Fungal communities in ancient peatlands developed from different periods in the Sanjiang Plain, China
Fungal communities in ancient peatlands developed from different periods in the Sanjiang Plain, China
Zhenqing Zhang 0 1 2
Xue Zhou 0 1 2
Lei Tian 0 1 2
Lina Ma 0 1 2
Shasha Luo 0 1 2
Jianfeng Zhang 0 1 2
Xiujun Li 0 1 2
Chunjie Tian 0 1 2
☯ These authors contributed equally to this work. 0 2
0 the Science Foundation of Chinese Academy of Sciences (XDB15030103), the National Basic Foundation (2016YFC0501202), Key Laboratory Foundation of Mollisols Agroecology (2016ZKHT- 05), the Key Technology Research and Development Program of CAS during the a13th Five-Year Plano (Y6H2043001) , the National
1 Key Laboratory of Mollisols Agroecology, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences , Changchun , China , 2 Key Laboratory of Wetland Ecology and Environment, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences , Changchun , China , 3 University of the Chinese Academy of Sciences , Beijing , China
2 Editor: Erika Kothe, Friedrich Schiller University , GERMANY
Peatlands in the Sanjiang Plain could be more vulnerable to global warming because they are located at the southernmost boundary of northern peatlands. Unlike bacteria, fungi are often overlooked, even though they play important roles in substance circulation in the peatland ecosystems. Accordingly, it is imperative that we deepen our understanding of fungal community structure and diversity in the peatlands. In this study, high-throughput Illumina sequencing was used to study the fungal communities in three fens in the Sanjiang Plain, located at the southern edge of northern peatlands. Peat soil was collected from the three fens which developed during different periods. A total of 463,198 fungal ITS sequences were obtained, and these sequences were classified into at least six phyla, 21 classes, more than 60 orders and over 200 genera. The fungal community structures were distinct in the three sites and were dominated by Ascomycota and Basidiomycota. However, there were no significant differences between these three fens in any α-diversity index (p > 0.05). Soil age and the carbon (C) accumulation rate, as well as total carbon (TC), total nitrogen (TN), C/N ratio, and bulk density were found to be closely related to the abundance of several dominant fungal taxa. We captured a rich fungal community and confirmed that the dominant taxa were those which were frequently detected in other northern peatlands. Soil age and the C accumulation rate were found to play important roles in shaping the fungal community structure.
Data Availability Statement: All relevant data are
within the paper and its Supporting Information
Peatlands, worldwide and particularly in the northern (boreal and subarctic) part, are shown
to be important participant in global C cycle in the recent past [
]. Despite covering only 6±8%
of the terrestrial ecosystems, northern peatlands store 550Pg C [
], accounting for 25±33% of
the world's soil C [
]. The rate of C accumulation is controlled more by changes in annual
Natural Science Foundation of China (41671094,
41571255, 31370144), and the Natural Science
Foundation of Jilin Province (20140101017JC).
plant productivity than by decomposition [
]. In particular, the reason for peatlands
functioning as a C sink is that more C is fixed by its vegetation than it is lost through the outflow of
dissolved organic C and emissions of CO2 and CH4 [
]. However, there is a growing concern
that the ongoing and future global warming will change the C cycling in these ecosystems,
making peatlands return the previously captured C to the atmosphere via releasing CO2 and/
or CH4, which would possibly accelerate the pace of global warming [7±10].
Many studies have shown that fungi play a crucial role in the soil ecosystem [11±13]. For
example, mycorrhizal fungi contribute to the growth of plants by providing root-available
nitrogen and phosphorus [
]. During the initial stage of plant litter decay, saprophagous
fungi are more important than bacteria in the degradation of fresh litter into less labile organic
]. The fungal communities, activities, and effects in the peatlands of America,
Europe, and Canada have been studied well. Lin et al. [
] found that dominant fungal taxa at
bog and fen sites in the peatlands of northwestern Minnesota were Ascomycota,
Basidiomycota, and Zygomycota, and the ammonia, dissolved organic carbon (DOC), and dissolved
organic nitrogen (DON) concentrations were more important than pH in structuring fungal
community structure. The same dominant fungal taxa were detected in the ombrotrophic
peatlands of northern England by Liu et al. [
]. They determined that fungal community
composition was significantly affected by organic matter (OM), moisture, phosphorus, and
ammonium. Sun et al. [
] reported that Ascomycota and Basidiomycota were the majority
groups in the boreal peatlands in Finland, and the dominant tree species and the availability of
certain nutrients (Ca, P, and Fe) had a strong effect on the composition of fungal communities.
Myers et al. [
] stated that fungal activity became increasingly important and was negatively
correlated with pH and total peat Mg concentrations in ombrotrophic peatlands located in
southern Canada. Soil age and the C accumulation rate were reportedly to be important factors
of fungal community structure in other ecosystems. Several studies found that fungal
communities shifted significantly across chronosequences in a forest ecosystem [20±22]. Zhang et al.
] found that C accumulation promoted macroaggregates formation and reduced the
effective diffusion coefficient of oxygen; thus, they identified changes in microhabitats and a shift
in microbial community when they investigated the linkage between organic C accumulation
and microbial community dynamics in a sandy loam soil.
However, there are few studies on the fungal community in the Sanjiang Plain, leading to
the lack of a high-quality curated database for taxonomic assignment. The Sanjiang Plain,
located in the temperate climate region, is the largest area of freshwater marshlands in China
]. These peatlands could be more vulnerable to global warming because they are primarily
located at the southern boundary of northern peatlands [
]. Thus, it is imperative that we
deepen our understanding of the microorganisms, particularly fungi, in this ecosystem. In this
study, high-throughput Illumina sequencing of ITS rRNA genes was used to study the fungal
communities in the Sanjiang Plain at the southern edge of northern peatlands. Peat soils were
collected from three fens which began development at different time periods. The greater
sequencing detail achieved by the high-throughput sequencing allowed the capture of the less
abundant and uncultured taxa, thus enabing a more thorough characterization of peatland
fungal diversity. The chronological sampling approach further facilitated the analysis of
potential linkages of fungal communities with soil age as well as the C accumulation rate.
Materials and methods
Study area and sampling description
Peat cores were sampled from three different minerotrophic fens Shenjiadian (S), Honghe
(H), and Qindelie (Q) in the Sanjiang Plain (129Ê110-135Ê050 E, 43Ê490-48Ê270 N), China
2 / 16
Fig 1. Map of the study region showing locations of the sampling sites in the Sanjiang Plain, which was
generated using ArcGIS 10.0 (http://www.esrichina.com.cn/).
(Fig 1). No specific permits were required for our field studies and our work did not involve
endangered or protected species. The three fens started to develop during different time
periods during the Holocene [
]. Over 70% of this region is dominated by fresh water wetlands
developing in ancient riverbeds and waterlogged depressions, and approximately 30% (or
nearly 3.3 × 104 ha) of this region is covered by peatlands, which developed under certain
topographic conditions during or before the Holocene [
]. The study area has a temperate humid
to sub-humid continental monsoon climate. The mean annual temperature ranges from 1.4ÊC
to 4.3ÊC, with an average maximum of 22ÊC in July and an average minimum of -18ÊC in
January. The mean annual precipitation is 500±650 mm, and 80% of rainfall occurs between May
and September [
]. The sampled fens receive water inputs from groundwater as well as
3 / 16
precipitation and are primarily covered with sedges (Carex lasiocarpa). In May 2012, triplicate
cores were sampled from each fen using a Russia peat core. The cores were subsampled for
chronological, microbiological, and physico-chemical analyses. The soil sample handling for
chronological analysis is described in more detail previously [
]. For microbiological and
physico-chemical analyses, surface peat soil (depth, 0±30 cm) was homogenized in sterile bags.
The peat soils were stored in the dark and kept chilled before they were transferred to a
laboratory. Chronological and physico-chemical characterization was conducted right after the
samples arrived. Samples for DNA extractions were immediately frozen at -80ÊC.
Chronological and physico-chemical characterization
Subsamples with a volume of 3 cm3 were used for loss-on-ignition (LOI) with sequential
combustion at 500ÊC to estimate organic matter content [
]. Bulk density was calculated at 1-cm
intervals based on the dry weight and volume of each subsample. Ash-free (organic matter)
bulk density was calculated from the measurements of the bulk density and organic matter
contents. All subsamples for accelerate mass spectrometry (AMS) dating were dated with an
accelerator mass spectrometry system at the Institute of Earth Environment, CAS. The AMS
14C dates were converted into calendar ages using the program Calib 7.02 based on the
INTCAL 13 calibration dataset [
]. The apparent C accumulation rates were calculated using
calibrated AMS 14C ages, ash-free bulk density measurements, and the C contents of peat
organic matter in peatlands (using 52% C in peat organic matter) [
]. Soil pH was measured
in a 1:5 soil/water suspension [
]. The total nitrogen (N) in the soil was determined by
dichromate oxidation using the continuous flow analytical system (SAN++, SKALAR,
DNA extraction and sequencing
Genomic DNA was extracted from 0.5 g of peat soil using a FASTDNA™ SPIN Kit for soil
(MPBio, Santa Ana, CA, USA), according to the manufacturer's instructions. DNA
concentration was measured using a NanoDrop 2000 spectrophotometer (NanoDrop Technologies,
Inc., Wilmington, DE, USA). The fungal ITS2 region was amplified using the primer pair
ITS3F (GCATCGATGAAGAACGCAGC) and ITS4R (TCCTCCGCTTATTGATATGC) combined
with Illumina adapter sequences and barcodes [
]. PCR reactions were performed in a 30-μL
mixture containing 3 μL of each primer (2 μM), 10 μL of template DNA (1 ng/μL), 15 uL of
Phusion1 High-Fidelity PCR Master Mix (BioLabs, Inc., New England, USA), and 2 μL of
water. The following thermal program was used for amplification: 95ÊC for 1 min, followed by
30 cycles of 98ÊC for 10 s, 50ÊC for 30 s, 72ÊCfor 30 s, and a final extension step at 72ÊC for 5
min. Each sample was amplified in triplicate, and the PCR products were pooled and purified
using the Qiagen Gel Extraction Kit (Qiagen, Hilden, Germany). Metagenomic sequencing
libraries were generated using the TruSeq1 DNA PCR Free Sample Preparation Kit (Illumina,
San Diego, CA, USA), according to the manufacturer's instructions and pooled at an
equimolar ratio. The 250-bp paired-end sequencing was performed on an Illumina HiSeq2000
platform at the Novogene Bioinformatics Technology Ltd.
Sequence data preprocessing and statistical analysis
Raw sequences were divided into sample libraries via sample-specific barcodes and were
truncated after cutting off the barcode and the primer sequence. Forward and reverse reads with at
least 10-bp overlaps and less than 5% mismatches were merged using FLASH [
filtering on the raw tags was performed according to the QIIME (V1.7.0, http://qiime.org/index.
html) quality control process [
], and all sequences shorter than 200 bp with an average
4 / 16
quality score lower than 25 in the raw reads were removed. The remaining sequences were
subjected to chimera removal using UCHIME Algorithm (http://www.drive5.com/usearch/
manual/uchime_algo.html). UPARSE (Version 7.0.1001, http://drive5.com/uparse/) was used
to classify the operational taxonomic units (OTUs) at the 97% similarity level [
taxonomic identity was annotated using a BLAST algorithm against sequences within the Unite
Database (https://unite.ut.ee/) using the QIIME software [
]. To study the phylogenetic
relationship of different OTUs and the differences between the dominant species in different
samples (groups), multiple sequence alignments were performed using the MUSCLE software
(Version 3.8.31, http://www.drive5.com/muscle/) [
]. Alpha and beta diversity analyses were
performed using harmonized data by random subsampling of 43,303 sequences. The alpha
diversity indices, including Chao1, Shannon [
], Simpson, ACE, and PieLou equitability [
were calculated. Beta diversity between microbial communities was evaluated using both
Unifrac (weighted and unweighted) and Jaccard distances. Principal coordinate analysis (PCoA)
and canonical correspondence analysis (CCA) were performed using the ape and vegan
packages in R, respectively. A linear discriminant analysis (LDA) effect size (LEfSe) algorithm was
used for high-dimensional biomarker discovery [
] using the non-parametric factorial
Kruskal-Wallis (KW) sum-rank test [
] to detect the features with having significant differential
abundance with respect to the class of interest; biological consistency was subsequently
investigated using with a set of pairwise tests among subclasses using the (unpaired) Wilcoxon
ranksum test [
]. As the final step, LEfSe used LDA to estimate the effect size of each differentially
abundant feature and perform dimension reduction, when necessary. Pearson correlation
analyses were used to correlate the relationships between the geochemical and microbial
parameters of the soil [
]. Differences in soil properties and alpha diversity indices across
samples were determined using ANOVA, which was followed by the least significant
difference (LSD) test performed in IBM SPSS (version 19.0, Chicago, IL, USA) [
beta-diversity community was compared by permutational MANOVA [
]. All the OTUs were assigned
to functional group by FUNGuild [
]. The Illumina sequencing data in the present study has
been deposited into NCBI SRA database with the accession number as SRP082472.
Physico-chemical and chronological characterization of peat
A total of nine peat cores were retrieved from S, H, and Q fens in the Sanjiang Plain,
northeastern China. The AMS dating results indicated that these fens had developed during different
periods and had different C accumulation rates. The peat was dated from 637 to 2,085 cal. yr
BP (Table 1). Peat cores from Q exhibited the most ancient chronological dates, being dated to
924±2085 cal. yr BP, and S was the youngest fen at 637±862 cal. yr BP. The C accumulation
rates of the peat cores ranged from 9.29 to 205.66 g C m−2 yr−1 (Table 1).
A summary of soil physico-chemical characteristics was presented in Table 1. Soils were
acidic; pH ranged from 4.80 to 5.44. Soil total C and N ranged from 286.88 to 430.12 g kg−1
and from 11.62 to 34.05 g kg−1, respectively.
Fungal community diversity
A total of 463,198 sequences targeting the ITS gene were obtained from nine surface (0±30
cm) soil samples using Illumina HiSeq sequencing, ranging from 43,303 to 57,269 reads per
sample. After OTU clustering at 97% sequence identity, a total of 989 OTUs were subsequently
generated after resampling with 43,303 sequences per sample. Shannon, Simpson, Chao1,
ACE index, and PieLou equitability were calculated to estimate microbial richness and
5 / 16
evenness (Table 2). However, there were no significant differences between these three fens in
all diversity indices (p > 0.05).
Microbial community composition
The obtained 463,198 fungal ITS sequences were classified into at least six phyla, 21 classes,
more than 60 orders, and over 200 genera, suggesting a rich fungal community in this
ecosystem. The dominant fungal phyla across all soil samples were Ascomycota and Basidiomycota,
with relative abundances ranging from 79.75% to 54.07% and 10.04% to 25.38%, respectively
(Fig 2). The relative abundances of the minor phyla Chytridiomycota, Zygomycota,
Glomeromycota, and Rozellomycota were all lower than 1%. In addition, numerous sequences could not be
classified into known fungi, with relative abundances varying from 1.65% to 19%. The 10 most
abundant fungal OTUs were found to belong to two different phyla (Ascomycota and
Basidiomycota) and three different orders (Agaricales, Helotiales, and Tremellales) (Fig 3).
Only 43.28% of the OTUs were assigned to a functional group by FUNGuild (S1 Table).
We detected 14 categories, broadly referred to as guilds; animal pathogens (1.65%), arbuscular
mycorrhizal fungi (1.42%), ectomycorrhizal fungi (3.30%), ericoid mycorrhizal fungi (0.35%),
lichenized fungi (0.94%), plant pathogens (4.25%), undefined saprotrophs (22.05%), wood
saprotrophs (1.18%), endophytes (2.36%), dung saprotrophs (0.83%), root associated
biotrophs (0.12%), soil saprotrophs (0.47%), fungal parasites (4.25%), and plant saprotrophs
(0.12%). Most of the OTUs were saprotrophs.
A Venn diagram was used to show a comparison of the similarities and differences between
the communities in the three fens (Fig 4a). The Shenjiadian, Honghe, and Qindelie fungal
Fig 2. Fungal community variation in the Shengjiadian (S), Honghe (H), and Qindelie (Q) fens. The mean relative abundance of fungi
from three cores taken from each site. Each bar height represents the relative abundance, and color represents a particular fungal phylum.
communities had 134 OTUs in common and 149, 228, 240 unique OTUs, respectively. The
unique OTUs accounted for 34%, 45% and 43% of the total detected OTUs in Shenjiadian,
Honghe, and Qindelie fens. The fungal diversity in these fens was rich, with OTU number,
Shannon index and Chao1 index ranging from 192 to 429, 1.516 to 5.027, 192.8 to 450.5,
respectively. The fungal community structures were distinct in the three sites (permutational
MANOVA, r = 0.375, p = 0.035) and dominated with sequences belonging to Ascomycota
To show the fungal community structures of S, H, and Q fens, PCoA based on both Unifrac
and Jaccard distance was performed (Fig 5). The highest relative abundance (80%) of
Ascomycota sequences was observed in Shenjiadian, the youngest fen. The three sites showed a similar
relative abundance of Zygomycota (0.35%-0.48%) sequences, while the Glomeromycota and
Chytridiomycota sequences were approximately five-fold and two-fold higher in Honghe than
in the other two sites.
The biomarkers explored using the LEfSe analysis were affiliated with Helotiaceae,
Tremellomycetes, and Eurotiomycetes, respectively (Fig 4b). According to the Pearson correlation
analysis, the relative abundance of Helotiaceae was positively correlated with the C accumulation
rate (r = 0.706, p = 0.033), and the relative abundance of Tremellomycetes was positively
correlated with soil age (r = 0.680, p = 0.044). However, the relative abundance of Eurotiomycetes
was not significantly correlated with any of measured soil parameter (p > 0.05).
Fungal distribution link to the soil properties
We used CCA to analyze the variation in fungal community structure and the associated
soil parameters. Soil parameters showing significant correlation (p < 0.05) with the fungal
7 / 16
Fig 3. Phylogenetic tree of the OTUs belonging to the top 10 genera. Each color represents a particular OTU. The second and the
outermost layer indicate relative abundance and confidence, respectively.
community were plotted as vectors (Fig 5). On the horizontal axis (CCA1, 19.30% of
constrained variability), the most influential constraining variable was C/N ratio (biplot
score = 0.82), followed by bulk density (biplot score = -0.36), TN (biplot score = -0.20), C
accumulation rate (biplot score = -0.11), soil age (biplot score = 0.07), and TC (biplot
score = 0.02). On the vertical axis (CCA2, 18.60% of constrained variability), the most
influential constraining variable was soil age (biplot score = 0.61), followed by TC (biplot
score = 0.52), C accumulation rate (biplot score = -0.45), bulk density (biplot score = -0.36),
C/N ratio (biplot score = 0.29), and TN (biplot score = -0.15).
8 / 16
Fig 4. Venn diagram (a) and LEfSe anaylsis (b and c) show the unique and shared OTUs between Shengjiadian (S), Honghe (H),
and Qindelie (Q) fens. In the cladogram (b), the circles radiating from the inside out represent fungal taxon from kingdom to family, and the
diameter of the circle is proportional to the relative abundance of each taxon. A taxon with significant difference is marked with the same
color as the sampling site where the taxon is ranked the highest, and the branch area is correspondingly shaded. A taxon without a
significant difference is marked in yellow. Graph (c) is a histogram of the LDA score, and fungal groups that are statistically significant among
sites are shown (p < 0.05).
Pearson correlation analyses were used to correlate the relationships between soil properties
and the relative abundance of the 10 most abundant fungal taxa at different levels (Fig 6). At
the phylum level, the relative abundance of only Basidiomycota was significantly (r = 0.708,
p = 0.033) correlated with the soil C/N ratio positively. Another dominant phylum exhibited
no correlation with measured soil properties. The correlation analysis of the dominant fungal
classes and soil properties revealed that the abundance of both Lecanoromycetes and
9 / 16
Fig 5. Similarity in fungal community structures indicated by PCoA plots based on weight UniFrac (a), unweight UniFrac (b) and
Jaccard (c) distances among sites.
Eurotiomycetes increased with increasing TN and C accumulation rate. At the order level, we
found that the relative abundance of Sebacinales was positively correlated with TN (r = 0.890,
p = 0.001) and the C accumulation rate (r = 0.847, p = 0.004). Pleosporales was negatively
correlated with soil bulk density (r = -0.727, p = 0.026) and positively correlated with soil age
(r = 0.873, p = 0.002). Filobasidiales also exhibited a positive correlation with soil age
(r = 0.681, p = 0.043). No significant correlation between soil properties and the relative
abundance of fungal taxa was observed at the genus level. Further, among all the measured soil
properties, soil pH was the only one which was found to be not correlated with fungal
abundance at any level (data not shown).
The fungal taxa observed in the present study primarily belonged to two phyla: Ascomycota
and Basidiomycota, which is consistent with previous studies of peat soils [46±48]. Previous
studies reported that Ascomycota and Basidiomycota were capable of aerobically degrading
dissolved organic matter (DOM), including cellulose and polyphenolic compounds [
predominance of Ascomycota and Basidiomycota at the surface-level peat soil was consistent with
this ability. The relative abundance of Ascomycota and Basidiomycota at the three sites were
varied, but the relative abundance of Zygomycota among these sites was similar (Fig 2),
suggesting that the physiology of Zygomycota may distinct from that of Ascomycota and
Basidiomycota. Zygomycota can survive over long periods of dormancy by producing thick-walled,
resistant spores [
]. Additionally, members of Zygomycota are not capable of using cellulose
and sucrose degradation products, but instead, they could use C substrates of animal and
fungal origin, such as fungal hyphae. For example, an important member of Zygomycete,
Mortierella sp., was observed in all samples. They can degrade chitin, the essential component of
fungal hyphae, as efficiently as chitinolytic actinomycetes [
]. We believe that Zygomycota
may play a role just as important as Ascomycota or Basidiomycota in the peatlands carbon
cycling. Saprotrophs were the dominant functional group. They play an important role in the
degradation of fresh litter to less labile organic matter during the initial stages of plant litter
decay . Biomarkers revealed by LEfSe were affiliated with Helotiaceae, Eurotiomycetes and
Tremellomycetes (Fig 4b). Eurotiomycetes are a class of Ascomycetes within the subphylum
Pezizomycotina. One subclass, Eurotiomycetidae, includes xerophiles and psychrophiles, which are
producers of toxic and useful secondary metabolites, such as fermentation agents used to
10 / 16
Fig 6. Canonical correspondence analysis (CCA) used to evaluate the effect of soil factors on the fungal community structure.
Soil factors indicated in blue text include TC (total carbon), TN (total nitrogen), C/N (C/N ratio), C.A.R. (carbon accumulation rate).
make food products and enzymes. It also includes the important genetic model Aspergillus
nidulans. Chaetothyriomycetidae, the other subclass, includes the common black yeasts, some
of which are pathogens of humans and animals, as well as some primarily lichen-forming
groups newly found to be phylogenetically associated with this group [
]. Helotiaceae is a
small group of fungicolous, lichenicolous, and discomycetes members [
]. Tremellomycetes, a
basal group in the Agaricomycotina, comprises mostly dimorphic taxa and species that form
hyphae and/or complex fruiting bodies [
]. We found that the relative abundance of
Tremellomycetes and Helotiaceae was positively correlated with soil age (r = 0.680, p = 0.044) and the
C accumulation rate (r = 0.706, p = 0.033), respectively. Peat samples collected from
Shenjiadian and Qindelie had relatively higher C accumulation rate and were older than those
11 / 16
collected from Honghe (Table 1). However, the relative abundance of Eurotiomycetes was not
significantly correlated with any of the measured soil parameters (p > 0.05). Vegetation can
affect the fungal community structures and may also be an important factor associated with
the biomarker. Although the S, H and Q fens were primarily covered with sedges (Carex
lasiocarpa), other vegetation and their proportion, such as the amount of shrubs and tree likely
differ. Faster C accumulation and older age likely mean that the surface of the Shenjiadan and
Qindelie fens were drier (with lower water tables) and that they had a greater percentage of
woody species, as a part of the minority vegetation. Further, Helotiaceae and Tremellomycetes
are dominate by terricolous or lichen-forming species, such as Gelatinopsis [
], two of the largest genera of Helotiaceae, and Tremella and Syzygospora
belonging to Tremellomycete [
]. The functional roles of these biomarkers in the three fens still
remain to be elucidated.
In this study, fungal communities is correlated with soil physico-chemical properties. Both
CCA and Pearson correlation analysis indicated that pH is less important as an environmental
force in shaping the fungal community structure. This could be attributed to the ability of
fungi to tolerate a wide range of pH and the fact that their optimal extracellular enzyme activity
is at low pH [
]. This study of only three acidic sites does not cover a broad enough pH range
to detect this relationship. Our results also suggested that TC, TN, and C/N ratio were
important in structuring fungal distribution, which was consistent with other studies, not only in
], but also in farmlands [
] and forests . All of these soil parameters were
related to nutrient availability, which may have obvious implications for fungi growth.
Furthermore, soil bulk density was identified as a driving force for fungal distribution in this
study. Soil bulk density is inversely related to porosity [
] as the pores serve as pathways for
water and oxygen, which are important growth factors of fungi [
]. We also found that the
fungal community was significantly correlated with soil age and C accumulation in peatlands.
The peat cores initiated from different ages, being affected by climatic conditions (e.g.,
temperature and precipitation), underwent different initiation and decomposition process, which
resulted in the respective soil properties [
] We speculate that the effect of soil age on fungal
community is due to its correlation with soil properties. Temperature, oxygen exposure time,
and vegetation have been proposed as important controls on the C accumulation rate of
peatlands and could yield a different C fraction, which may be used as potential substrate or cause
an inhibitory affect on microorganism growth [
]. Fungi, as an important participant of C
cycling, could change the C accumulation rate of peatlands. The correlation between C
accumulation and microbial community has been previously investigated. Mackay et al. found the
positive correlations between total organic C accumulation and microbial indicators (total
PLFA, fungal PLFA, bacterial PLFA and activities of decomposition enzymes) when they
surveyed a chronosequence (0±23 years) of reforesting riparian pastures [
]. Kaƚucka and
] reported that ectomycorrhizal fungi were either indirectly involved in the C
accumulation through their influence on tree biomass production and organic matter
decomposition in boreal and temperate forests or directly involved through C storage by the
mycelium. A mechanistic understanding of the role of important fungal taxa in peatland C cycling
requires further field experiments and ecophysiological studies in the laboratory. However,
several studies demonstrated that fungal community distribution pattern is mainly affected by
pH, peatland vegetation, DOC, or DON [
6, 12, 61
]. In addition, hydrology and the
hydroperiod dictates the above-ground vegetation, and this in turn is both influenced by and influences
fungal and bacterial communities [
]. The peatland fungal community distribution pattern
remains to be sorted out as the technology and methodology are developing and each wetland
is unique along a hydrological spectrum.
12 / 16
In conclusion, the composition and diversity of fungal communities in three different
minerotrophic fens distributed in the Sanjiang Plain, the southern edge of northern peatlands in
China, were investigated by high-throughput Illumina sequencing. We captured a rich fungal
community and confirmed that the dominating taxa have also been frequently detected in
other northern peatland ecosystems. The fungal community structures were distinct in the
three sites, however, there was no significant differences in α-diversity (p > 0.05). Unlike pH,
TC, TN, C/N ratio, and bulk density are determined to be important environmental
parameters shaping fungal community structure. Additionally, we found that the distribution patterns
of several abundant fungal taxa were closely related to the soil age and C accumulation rate.
However, because of limited sampling sites in the research, these results should be considered
preliminary. Data presented here may serve useful for future researchers in this field. Extensive
sampling of a larger range of peatland sites is required to further this research, and may reveal
a relationship between the fungal community and pH, which was not detected in this study.
S1 Fig. The relationship between the relative abundances of abundant fungal phyla (A),
fungal classes (B), fungal orders (C), fungal genera (D) and soil properties.
S1 Table. Annotation results assigned by FUNGuild.
This work was financially supported by the Science Foundation of Chinese Academy of
Sciences (XDB15030103), the National Basic Foundation (2016YFC0501202), Key Laboratory
Foundation of Mollisols Agroecology (2016ZKHT-05), the Key Technology Research and
Development Program of CAS during the ª13th Five-Year Planº (Y6H2043001), the National
Natural Science Foundation of China (41671094, 41571255, 31370144), and the Natural
Science Foundation of Jilin Province (20140101017JC).
Conceptualization: Chunjie Tian.
Data curation: Xue Zhou.
Formal analysis: Lei Tian, Jianfeng Zhang.
Funding acquisition: Chunjie Tian.
Investigation: Zhenqing Zhang, Lina Ma.
Methodology: Zhenqing Zhang, Shasha Luo.
Software: Lei Tian.
Supervision: Chunjie Tian.
Writing ± original draft: Xue Zhou.
Writing ± review & editing: Xue Zhou, Xiujun Li.
13 / 16
14 / 16
15 / 16
1. Yu Z , Beilman DW , Jones MC . Sensitivity of northern peatland carbon dynamics to Holocene climate change. Carbon cycling in northern peatlands . Geophysical Monograph . 2009 ; 184 : 55 ± 69 .
2. Philben M , Holmquist J , MacDonald G , Duan D , Kaiser K , Benner R. Temperature , oxygen, and vegetation controls on decomposition in a James Bay peatland . Global Biogeochemical Cycles . 2015 ; 29 ( 6 ): 729 ± 43 . https://doi.org/10.1002/2014gb004989
3. Utstol-Klein S , Halvorsen R , Ohlson M. Increase in carbon accumulation in a boreal peatland following a period of wetter climate and long-term decrease in nitrogen deposition . New Phytologist . 2015 ; 206 ( 4 ): 1238 ± 46 . https://doi.org/10.1111/nph.13311 PMID: 25678224 .
4. Xing W , Bao K , Gallego-Sala AV , Dan JC , Zhang Z , Gao C , et al. Climate controls on carbon accumulation in peatlands of Northeast China . Quaternary Science Reviews . 2015 ; 115 ( 9 ): 78 ± 88 .
5. Andersen R , Chapman SJ , Artz RRE . Microbial communities in natural and disturbed peatlands: A review . Soil Biology and Biochemistry . 2013 ; 57 : 979 ± 94 . https://doi.org/10.1016/j.soilbio. 2012 . 10 .003
6. Bragazza L , Parisod J , Buttler A , Bardgett RD . Biogeochemical plant-soil microbe feedback in response to climate warming in peatlands . Nature Climate Change . 2013 ; 3(3):273±7 . https://doi.org/10.1038/ nclimate1781
7. Christensen JH , Christensen OB . A summary of the PRUDENCE model projections of changes in European climate by the end of this century . Climatic Change . 2007 ; 81 :7± 30 . https://doi.org/10.1007/ s10584-006-9210-7
8. Davidson EA , Janssens IA . Temperature sensitivity of soil carbon decomposition and feedbacks to climate change . Nature . 2006 ; 440 ( 7081 ): 165 ± 73 . https://doi.org/10.1038/nature04514 PMID: 16525463
9. Jones MC , Yu Z. Rapid deglacial and early Holocene expansion of peatlands in Alaska . Proceedings of the National Academy of Sciences of the United States of America . 2010 ; 107 ( 16 ): 7347 ± 52 . https://doi. org/10.1073/pnas.0911387107 PMID: 20368451
10. Knorr W , Prentice IC , House JI , Holland EA . Long-term sensitivity of soil carbon turnover to warming . Nature . 2005 ; 433 ( 7023 ): 298 ± 301 . https://doi.org/10.1038/nature03226 PMID: 15662420
11. Kim S-Y , Lee S-H , Freeman C , Fenner N , Kang H . Comparative analysis of soil microbial communities and their responses to the short-term drought in bog, fen, and riparian wetlands . Soil Biology & Biochemistry . 2008 ; 40 ( 11 ): 2874 ± 80 . https://doi.org/10.1016/j.soilbio. 2008 . 08 .004
12. Lin X , Green S , Tfaily MM , Prakash O , Konstantinidis KT , Corbett JE , et al. Microbial community structure and activity linked to contrasting biogeochemical gradients in bog and fen environments of the Glacial Lake Agassiz Peatland . Applied and Environmental Microbiology . 2012 ; 78 ( 19 ): 7023 ± 31 . https:// doi.org/10.1128/AEM.01750-12 PMID: 22843538
13. Peltoniemi K , Strakova P , Fritze H , Alvira Iraizoz P , Pennanen T , Laiho R . How water-level drawdown modifies litter-decomposing fungal and actinobacterial communities in boreal peatlands . Soil Biology & Biochemistry . 2012 ; 51 : 20 ± 34 . https://doi.org/10.1016/j.soilbio. 2012 . 04 .013
14. Kiers ET , Duhamel M , Beesetty Y , Mensah JA , Franken O , Verbruggen E , et al. Reciprocal rewards stabilize cooperation in the mycorrhizal symbiosis . Science . 2011 ; 333 ( 6044 ): 880 ±2. https://doi.org/10. 1126/science.1208473 PMID: 21836016
15. Kuehn KA , Lemke MJ , Suberkropp K , Wetzel RG . Microbial biomass and production associated with decaying leaf litter of the emergent macrophyte Juncus effusus . Limnology and Oceanography . 2000 ; 45 ( 4 ): 862 ± 70 .
16. Lin X , Kennedy D , Fredrickson J , Bjornstad B , Konopka A . Vertical stratification of subsurface microbial community composition across geological formations at the Hanford Site . Environmental Microbiology . 2012 ; 14 ( 2 ): 414 ± 25 . https://doi.org/10.1111/j.1462- 2920 . 2011 . 02659 . x PMID : 22122741
17. Liu J , Sui Y , Yu Z , Shi Y , Chu H , Jin J , et al. Soil carbon content drives the biogeographical distribution of fungal communities in the black soil zone of northeast China . Soil Biology & Biochemistry . 2015 ; 83 : 29 ± 39 . https://doi.org/10.1016/j.soilbio. 2015 . 01 .009
18. Sun H , Terhonen E , Kovalchuk A , Tuovila H , Chen H , Oghenekaro AO , et al. Dominant tree species and soil type affect the fungal community structure in a boreal peatland forest . Applied and Environmental Microbiology . 2016 ; 82 ( 9 ): 2632 ± 43 . https://doi.org/10.1128/AEM.03858-15 PMID: 26896139
19. Myers B , Webster KL , McLaughlin JW , Basiliko N. Microbial activity across a boreal peatland nutrient gradient: the role of fungi and bacteria . Wetlands Ecology and Management . 2012 ; 20 ( 2 ): 77 ± 88 . https:// doi.org/10.1007/s11273-011-9242-2
20. Hui N , Jumpponen A , Francini G , Kotze DJ , Liu X , Romantschuk M , et al. Soil microbial communities are shaped by vegetation type and park age in cities under cold climate . Environmental Microbiology . 2017 ; 19 ( 3 ): 1281 ± 95 . https://doi.org/10.1111/ 1462 - 2920 .13660 PMID: 28063185
21. Kyaschenko J , Clemmensen KE , Hagenbo A , Karltun E , Lindahl BD . Shift in fungal communities and associated enzyme activities along an age gradient of managed Pinus sylvestris stands . Isme Journal . 2017 ; 11 ( 4 ): 863 ± 74 . https://doi.org/10.1038/ismej. 2016 .184 PMID: 28085155
22. Sun S , Li S , Avera BN , Strahm BD , Badgley BD . Soil bacterial and fungal communities show distinct recovery patterns during forest ecosystem restoration . Applied and Environmental Microbiology . 2017 , 83 ( 14 ): https://doi.org/10.1128/aem.00966-17 PMID: 28476769 .
23. Zhang H , Ding W , Yu H , He X . Linking organic carbon accumulation to microbial community dynamics in a sandy loam soil: result of 20 years compost and inorganic fertilizers repeated application experiment . Biology and Fertility of Soils . 2015 ; 51 ( 2 ): 137 ± 50 . https://doi.org/10.1007/s00374-014-0957-0
24. Xing W , Guo W , Liang H , Li X , Wang C , He J , et al. Holocene peatland initiation and carbon storage in temperate peatlands of the Sanjiang Plain, Northeast China . Holocene. 2016 ; 26 ( 1 ): 70 ±9. https://doi. org/10.1177/0959683615596824
25. Zhang Z , Xing W , Wang G , Tong S , Lv X , Sun J. The peatlands developing history in the Sanjiang Plain , NE China, and its response to East Asian monsoon variation . Scientific Reports . 2015 ; 5 . https://doi. org/10.1038/srep11316 PMID: 26076653
26. Ma X , Liu X , Wang R . China's wetlands and agro-ecological engineering . Ecological Engineering . 1993 ; 2 ( 3 ): 291 ± 301 .
27. Song C , Xu X , Tian H , Wang Y . Ecosystem-atmosphere exchange of CH4 and N2O and ecosystem respiration in wetlands in the Sanjiang Plain, Northeastern China . Global Change Biology . 2009 ; 15 ( 3 ): 692 ± 705 . https://doi.org/10.1111/j.1365- 2486 . 2008 . 01821 .x
28. D WE Jr . Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: comparison with other methods . Journal of Sedimentary Research . 1974 ; 44 ( 1 ): 242 ± 8 .
29. Stuiver M. Extended 14C data base and revised CALIB 3.0 14C age calibration program . Radiocarbon . 1993 ; 35 : 215 ± 30 .
30. Vitt DH , Halsey LA , Bauer IE , Campbell C . Spatial and temporal trends in carbon storage of peatlands of continental western Canada through the Holocene . Canadian Journal of Earth Sciences . 2000 ; 37 ( 5 ): 683 ± 93 . https://doi.org/10.1139/e99- 097
31. Zhang B , Chen S , He X , Liu W , Zhao Q , Zhao L , et al. Responses of soil microbial communities to experimental warming in alpine grasslands on the Qinghai-Tibet Plateau . Plos One . 2014 ; 9 ( 8 ). https:// doi.org/10.1371/journal.pone. 0103859 PMID: 25083904
32. Caporaso JG , Lauber CL , Walters WA , Berg-Lyons D , Huntley J , Fierer N , et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms . Isme Journal . 2012 ; 6 (8):1621±4 . https://doi.org/10.1038/ismej. 2012 .8 PMID: 22402401
33. Deng J , Gu Y , Zhang J , Xue K , Qin Y , Yuan M , et al. Shifts of tundra bacterial and archaeal communities along a permafrost thaw gradient in Alaska . Molecular Ecology . 2015 ; 24 ( 1 ): 222 ± 34 . https://doi.org/10. 1111/mec.13015 PMID: 25424441
34. Bokulich NA , Subramanian S , Faith JJ , Gevers D , Gordon JI , Knight R , et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing . Nature Methods . 2013 ; 10 ( 1 ):57± U11. https://doi.org/10.1038/nmeth.2276 PMID: 23202435
35. Edgar RC . MUSCLE: multiple sequence alignment with high accuracy and high throughput . Nucleic Acids Research . 2004 ; 32 ( 5 ): 1792 ±7. https://doi.org/10.1093/nar/gkh340 PMID: 15034147
36. Koljalg U , Nilsson RH , Abarenkov K , Tedersoo L , Taylor AFS , Bahram M , et al. Towards a unified paradigm for sequence-based identification of fungi . Molecular Ecology . 2013 ; 22 ( 21 ): 5271 ±7. https://doi. org/10.1111/mec.12481 PMID: 24112409
37. Chao A , Bunge J . Estimating the number of species in a stochastic abundance model . Biometrics . 2002 ; 58 ( 3 ): 531 ±9. https://doi.org/10.1111/j. 0006 - 341X . 2002 . 00531 . x PMID : 12229987
38. Simon M , Lopez-Garcia P , Deschamps P , Moreira D , Restoux G , Bertolino P , et al. Marked seasonality and high spatial variability of protist communities in shallow freshwater systems . Isme Journal . 2015 ; 9 ( 9 ): 1941 ±65. https://doi.org/10.1038/ismej. 2015 .6 PMID: 25853803
39. Segata N , Izard J , Waldron L , Gevers D , Miropolsky L , Garrett WS , et al. Metagenomic biomarker discovery and explanation . Genome Biology . 2011 ; 12 ( 6 ). https://doi.org/10.1186/gb-2011 -12-6-r60 PMID: 21702898
40. Kruskal WH , Wallis WA . Citation-classic-use of ranks in one-criterion varivance analysis . Current Contents/Social & Behavioral Sciences . 1987 ;(40): 20 -.
41. Wilcoxon F. Individual comparisons of grouped data by ranking methods . Journal of economic entomology . 1946 ; 39 : 269 -. PMID: 20983181.
42. Oksanen J , Kindt R , Legendre P , O'Hara B , Stevens M , Oksanen M. The vegan package . vegan: Community Ecology Package . 2010 .
43. Banerjee S , Kirkby CA , Schmutter D , Bissett A , Kirkegaard JA , Richardson AE . Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil . Soil Biology & Biochemistry . 2016 ; 97 : 188 ± 98 . https:// doi.org/10.1016/j.soilbio. 2016 . 03 .017
44. Jiao S , Liu Z , Lin Y , Yang J , Chen W , Wei G . Bacterial communities in oil contaminated soils: Biogeography and co-occurrence patterns . Soil Biology & Biochemistry . 2016 ; 98 : 64 ± 73 . https://doi.org/10. 1016/j.soilbio. 2016 . 04 .005
45. Nguyen NH , Song Z , Bates ST , Branco S , Tedersoo L , Menke J , et al. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild . Fungal Ecology . 2016 ; 20 : 241 ±8. https:// doi.org/10.1016/j.funeco. 2015 . 06 .006
46. Artz RRE , Anderson IC , Chapman SJ , Hagn A , Schloter M , Potts JM , et al. Changes in fungal community composition in response to vegetational succession during the natural regeneration of cutover peatlands . Microbial Ecology . 2007 ; 54 ( 3 ): 508 ± 22 . https://doi.org/10.1007/s00248-007 -9220-7 PMID: 17450394
47. Gilbert D , Mitchell EAD. Chapter 13 Microbial diversity in Sphagnum peatlands: Elsevier Science & Technology; 2006 . 287 ±318 p.
48. Thormann MN . Diversity and function of fungi in peatlands: A carbon cycling perspective . Canadian Journal of Soil Science . 2006 ; 86 ( 2 ): 281 ± 93 .
49. Lynd LR , Weimer PJ , van Zyl WH , Pretorius IS . Microbial cellulose utilization: Fundamentals and biotechnology (vol 66 , pg 506, 2002 ). Microbiology and Molecular Biology Reviews . 2002 ; 66 ( 4 ): 739 -. https://doi.org/10.1128/mmbr.66.4.739.2002
50. Young-Ju K , Zhao Y , Oh K-T , Nguyen V-N , Park R -D. Enzymatic deacetylation of chitin by extracellular chitin deacetylase from a newly screened Mortierella sp DY-52 . Journal of Microbiology and Biotechnology . 2008 ; 18 ( 4 ): 759 ± 66 . PMID: 18467873
51. De Boer W , Gerards S , Gunnewiek PJA , Modderman R . Response of the chitinolytic microbial community to chitin amendments of dune soils . Biology and Fertility of Soils . 1999 ; 29 ( 2 ): 170 ±7. https://doi.org/ 10.1007/s003740050541
52. Geiser DM , Gueidan C , Miadlikowska J , Lutzoni F , Kauff F , Hofstetter V , et al. Eurotiomycetes: Eurotiomycetidae and Chaetothyriomycetidae. Mycologia . 2006 ; 98 ( 6 ): 1053 ± 64 . https://doi.org/10.3852/ mycologia.98.6.1053 PMID: 17486980
53. Zhuang WY . Two new species of Unguiculariopsis (Helotiaceae, Encoelioideae) from China . Mycological Research . 2000 ; 104 : 507 ±9. https://doi.org/10.1017/s0953756299001902
54. Weiss M , Bauer R , Sampaio JP , Oberwinkler F. 12 Tremellomycetes and Related Groups2014 . 331 ±55 p.
55. KocourkovaÂ J , Knudsen K. Gelatinopsis acarosporicola (Helotiaceae), a new lichenicolous fungus on Acarospora socialis from western North America . Bryologist. 2009 ; 112 ( 2 ): 363 ± 7 .
56. Zheng HD , Zhuang WY . Five new species of Hymenoscyphus (Helotiaceae, Ascomycota) with notes on the phylogeny of the genus . Mycotaxon . 2016 ; 130 ( 4 ): 1017 ± 38 .
57. Millanes AM , Diederich P , Ekman S , Wedin M. Phylogeny and character evolution in the jelly fungi (Tremellomycetes, Basidiomycota , Fungi). Molecular Phylogenetics & Evolution . 2011 ; 61 ( 1 ): 12 .
58. Beales N. Adaptation of Microorganisms to Cold Temperatures, Weak acid preservatives, low pH, and osmotic stress: A review . Comprehensive Reviews in Food Science & Food Safety . 2004 ; 3 ( 1 ):1± 20 .
59. Lauber CL , Strickland MS , Bradford MA , Fierer N. The influence of soil properties on the structure of bacterial and fungal communities across land-use types . Soil Biology & Biochemistry . 2008 ; 40 ( 9 ): 2407 ± 15 .
60. Li CH , Ma BL , Zhang TQ . Soil bulk density effects on soil microbial populations and enzyme activities during the growth of maize (Zea mays L.) planted in large pots under field exposure . Canadian Journal of Soil Science . 2002 ; 82 ( 2 ): 147 ± 54 .
61. Elliott DR , Caporn SJM , Nwaishi F , Nilsson RH , Sen R . Bacterial and fungal communities in a degraded ombrotrophic peatland undergoing natural and managed re-vegetation . Plos One . 2015 ; 10 ( 5 ). https:// doi.org/10.1371/journal.pone. 0124726 PMID: 25969988
62. Mackay JE , Cunningham SC , Cavagnaro TR . Riparian reforestation: are there changes in soil carbon and soil microbial communities? Science of the Total Environment . 2016 ; 566 ± 567 : 960 .
63. Kaøucka I , Jagodziński AM . Ectomycorrhizal fungi and carbon dynamics in forest ecosystems . Sylwan . 2013 ; 157 ( 11 ): 817 ± 30 .