Fungal Communities Along a Small-Scale Elevational Gradient in an Alpine Tundra Are Determined by Soil Carbon Nitrogen Ratios.

Elevational gradients are associated not only with variations in temperature and precipitation, but also with shifts in vegetation types and changes in soil physicochemical properties. While large-scale elevational patterns of soil microbial diversity, such as monotonic declines and hump-shaped models, have been reported, it is unclear whether within-ecosystem elevational distribution patterns exist for soil fungal communities at the small scale. Using Illumina Miseq DNA sequencing, we present a comprehensive analysis of soil fungal diversity and community compositions in an alpine tundra ecosystem at elevations ranging from 2000 to 2500 m on the Changbai Mountain, China. Soil fungal community composition differed among elevations, and the fungal diversity (i.e., species richness and Chao1) increased along elevations. Soil fungal richness was negatively correlated with soil carbon/nitrogen (C/N) ratio, and community composition varied according to the C/N ratio. In addition, the relative abundances of Basidiomycota and Leotiomycetes were similarly negatively correlated with C/N ratio. For functional guilds, our data showed that mycoparasite and foliar epiphyte abundances were also influenced by C/N ratio. These results indicated that soil C/N ratio might be a key factor in determining soil fungal distribution at small-scale elevational gradients.

Introduction

The study of elevational diversity patterns contributes to important insights for developing the general theory of species diversity (Lomolino, 2001). It also provides novel perspectives on the potential responses of biological communities and ecosystems to the changing environments (Sundqvist et al., 2013). While there is a long history of mapping the elevational distribution and abundance of macroorganisms across terrestrial ecosystems (Kessler et al., 2011; Descombes et al., 2017; Dvorsky et al., 2017; Forister et al., 2010), less is known about patterns of microbial elevational diversity. Some studies on bacterial diversity along elevational gradients have showed a range of microbial responses to increasing elevation, including decreasing (Bryant et al., 2008), unimodal (Singh et al., 2012), and inconsistent trends (Fierer et al., 2011; Shen et al., 2013; Xu et al., 2014; Wang et al., 2017). It has been demonstrated that bacterial diversity and community composition are affected by various factors, such as soil pH (Shen et al., 2013; Zhang et al., 2013), soil carbon and n class="Chemical">nitrogen content (Lin et al., 2015; Shen et al., 2015), and spatial distribution (Hayden and Beman, 2016), however, similar studies of soil fungi remain limited.

Soil is a vital habitat for morphologically and functionally diverse fungal groups, including decomposers, mutualists, and pathogens of plants and animals (Tedersoo et al., 2014, 2017). Fungal activity is known to influence the structures of plant and animal communities, as well as rates of ecosystem processes (Peay et al., 2016), but despite these critical roles in the ecosystem, soil fungal diversity, and biogeographic patterns remain generally under-studied and undefined (Pautasso, 2013). It is thought that drivers of soil fungal spatial distribution patterns vary according to ecosystem. For example, soil fungal communities were predominantly determined by the aboveground plant community in natural grasslands of Tibetan Plateau (Yang et al., 2017a), while in a high elevation desert, the fungal communities were determined by highly stochastic dispersal of fungal species (Yang et al., 2017b). In addition, a temperature was found to be the strongest predictor of soil fungal diversity in the maritime Antarctic (Newsham et al., 2016).

Altitude (elevation) has been shown to affect fungal communities, where, for example, lower phylogenetic diversity of alpine arbuscular mycorrhizal fungal communities were recorded at higher elevations in the North American Rocky Mountains (Egan et al., 2017). Drivers of soil fungal communities at varying elevations include climate variables: Bahram et al. (2012) found that monotonically declining species richness of ectomycorrhizal fungi with increasing elevation in temperate old-growth forests was constrained by the variation in mean annual precipitation and temperature. The effects of soil properties on elevational differences in microbial communities, however, are unclear. At a large-scale elevational mountain gradient, Shen et al. (2014) found that although soil eukaryotic microbial richness was not correlated with elevation, it was strongly correlated with soil pH on the Changbai Mountain, China. However, Ping et al. (2017) observed on the same mountain that soil fungal richness and evenness decreased with the increasing elevation gradient, and the community composition was significantly influenced by soil C/N ratio, moisture, pH, and total n class="Chemical">nitrogen (TN), suggesting that this elevational diversity pattern was scale dependent.

To date, most studies examining vertical fungal diversity patterns have focused on contrasting habitats and vegetation types over large-scale elevational gradients, with little attention to dynamics over smaller scales or within biomes (Guo et al., 2013; Sundqvist et al., 2013). Tundra biomes are under threat from climate change, since they are highly sensitive to the global warming, but the fungal communities in Changbai Mountain Tundra are little studied. Shi et al. (2015) found that the soil fungal community composition in an n class="Disease">arctic tundra ecosystem was significantly affected by pH, ammonium concentration, C/N ratio, dissolved organic nitrogen (DON) content, and soil moisture content, however, it is not known whether these drivers of fungal community composition also apply in alpine tundra ecosystems.

The Changbai Natural Reserve, which encompasses within Changbai Mountain, is one of the largest biosphere reserves in China and has been spared from logging and other damaging human impacts (He et al., 2005). The vertical zonation of vegetation on Changbai Mountain transitions from broadleaf forest through to mixed, coniferous, and birch forests to n class="Species">alpine tundra (Dai et al., 2002). To increase our understanding of the structures of alpine tundra soil fungal communities, we investigated soil fungal diversity and community composition and soil physicochemical properties in the Changbai Mountain alpine tundra, along a small elevation gradient from 2000 to 2500 m. Specifically, we want to answer whether fungal diversity and community composition would change along the small-scale elevational gradient, and if so, what are the key drivers of the trend.

Materials and Methods

Study Site and Soil Sampling

The study site was located on Changbai Mountain (41°23′N–42°36′N, 126°55′E–129°00′E,) in northeast China, where alpine tundra occurs above 2000 m. The site and soil sampling methods have previously been fully described by Shen et al. (2013, 2015), so we provide a brief description here. On July 29, 2011, we selected n class="Species">alpine tundra study sites at six elevations between 2000 and 2500 m, separated by 100 m from the northern slope of Changbai Mountain. At each site, we created four 10 m × 10 m plots and from each of these, we randomly collected and combined six organic layer soil samples from an area of c. 10 cm × 10 cm and depth of 0–5 cm below the litter layer (N = 24). The fresh soil samples were sieved through a 2 mm sieve and visible roots and other residues were removed. Each of the 24 samples was divided into two subsamples, one was stored at -40°C for DNA extraction and the other was stored at 4°C for measurement of soil physicochemical properties. The main soil characteristics at the different elevations have been shown in the study of Shen et al. (2015).

Soil Physicochemical Properties Analyses and Soil DNA Extraction

Soil pH, moisture, dissolved organic carbon (n class="Chemical">DOC), Nitrate (NO3--N), ammonium (NH4+-N), DON, total carbon (TC), and TN were measured from each of the soil samples using methods described in detail by Shen et al. (2015). In brief, soil pH was measured after vigorous shaking water (1:5 w/v) supernatants for 30 min and soil moisture was measured gravimetrically. DOC was analyzed using a TOC analyzer (Analytik Jena, Multi N/C 3000, Germany) and NO3--N, NH4+-N, and DTN were determined using continuous flow analytical system (Skalar, Holland). DON was calculated according the content of NO3--N, NH4+-N, and DTN. Elemental analyzer (Vario MAX, Germany) was used to measure total TC and TN contents.

Soil DNA was extracted from 0.5 g fresh soil using a Fast DNA® SPIN Kit for soil (MP Biomedicals, Santa Ana, CA, United States) according to the manufacturer’s instructions and crude DNA was purified by Power Clean® Pro DNA Kit (MO BIO Laboratories, Inc.). The purified DNA was quantified by a Nano Drop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, NC, United States) and then stored at -20°C for use.

PCR Amplification and Illumina Miseq Sequencing

Primers ITS3F (5′-GCATCGATGAAGAACGCAGC-3′) and ITS4F (5′-TCCTCCGCTTATTGATATGC-3′) (Balajee et al., 2007) were used to amplify the ITS2 region. PCR was performed in 50 μl of reaction mixture containing 25 μl of Premix Taq DNA polymerase (TaKaRa, Japan), 0.5 μl of forward and reverse primers (20 μM), 23 μl of double distilled water (ddH2O), and 1 μl of DNA template (20 ng total soil DNA). The PCR cycling conditions were 94°C for 5 min, 32 cycles of 94°C for 30 s, 54°C for 30 s, and 72°C for 1 min. The PCR reaction was terminated at 72°C for 10 min and then cooled to 4°C. The PCR products were purified using EasyPure Quick Gel Extraction Kits (TransGen Biotech, Beijing, China) and quantified by the NanoDrop ND-1000 and then sequenced using the Illumina MiSeq PE 250 platform (Caporaso et al., 2012).

Processing of Fungal ITS Sequencing Date

Sequences obtained by Illumina Miseq sequencing were processed and analyzed using the QIIME software package (version 1.9.0) following the default settings (Caporaso et al., 2010). Low quality sequences shorter than 200 bp were removed during filtering and chimera checking was performed using UCHIME algorithm (Edgar et al., 2011) by the USEARCH tool. Sequences were assigned to operational taxonomic units (OTUs) at 97% similarity, using the UCLUST algorithm (Edgar, 2010), which were used as measures of fungal species richness. Taxonomy was assigned to all representative sequences using the ribosomal database project (RDP) classifier based on the UNITE fungal ITS database (QIIME release, version 7.0) (Bengtsson-Palme et al., 2013; Nilsson et al., 2015). Singletons and non-fungal sequences were discarded and 39,900 sequences per sample were randomly subsampled to rarify the data sets to the same level for further analysis.

Statistical Analyses

Alpha diversity was estimated using both taxonomic metrics (observed numbers of OTUs and Chao1 index). Chao1 was used to estimate the OTU richness based on frequencies of doubletons and singletons (Chao et al., 2009). Pearson’s correlation analysis in SPSS Statistin class="Chemical">cs Version21 for Windows (IBM SPSS, United States) was used to test for correlations (P < 0.05) between fungal diversity and relative fungal abundance and elevation and soil physicochemical variable.

Principal Co-ordinates Analysis (PCoA) was used to compare beta diversity between samples along the elevations based on the Bray-Curtis distance matrix and visualized using the ggplot2 package in R statistical software (version 3.1.2). In addition, we tested for differences in community composition among the elevations using analysis of similarities (ANOSIM) (Clarke, 1993; Warton et al., 2012), based on Bray-Curtis distances, in the vegan package of R. Meanwhile, a permutational multivariate analysis of variance (PERMANOVA) was also used to determine the effects of elevation. We then analyzed beta diversity of fungi along the elevation gradient using non-metric multidimensional scaling (NMDS), based on Bray-Curtis distance, in the “vegan” package of R; NMDS analysis was performed with the metaMDS function (Oksanen et al., 2013).

We used distance-based multivariate linear model analysis (DistLM) (McArdle and Anderson, 2001) to assess the influence of the soil physicochemical variables on fungal community composition and marginal tests were used to assess the statistical significance and contribution of each variable separately. These tests were performed using the computer program DISTLM_forward3 (Anderson, 2003). Fungal functional guilds were assigned according to Tedersoo et al. (2014) and Nguyen et al. (2016a) using an open annotation tool (FUNGuild)[1]. Here, we only accepted the guild assignment that confidence ranking were “probable” and “highly probable,” which was recommended by Nguyen et al. (2016a). Fourteen fungal functional groups were detected including Ectomycorrhizal fungi, Plant pathogens, Mycoparasites, Wood saprotrophs, Dung saprotrophs, Animal pathogens, Ericoid mycorrhizal fungi, Lichenized, Lichenicolous, Foliar epiphytes, Arbuscular mycorrhizal fungi, Endophytes, Orchid mycorrhizal fungi, and Soil saprotrophs.

Nucleotide Sequence Accession Numbers

Sequences were deposited in the NCBI Sequence Read Archive (SRA) under the accession number SRP110350[2].

Results

Soil Physicochemical Properties

Soil physicochemical properties including soil pH, soil moisture, DOC, n class="Chemical">DON, NO3--N, NH4+-N, TC, TN, and C/N ratio have been described in our previous study (Shen et al., 2015). Duncan’s multiple range tests showed that all the measured soil properties were significantly different between some elevations, except for pH (Shen et al., 2015). In addition, soil C/N ratio, the contents of TC and DOC were significantly negatively correlated with elevation (Supplementary Table ).

Fungal Diversity and Community Composition Along the Elevation Gradient

A total of 1,173,608 quality sequences were obtained from the 24 soil samples, ranging from 39,905 to 60,636 per soil sample (average = 48,900). We found some differences in fungal richness between the different elevations (Supplementary Table ), where species richness and Chao1 index of soil fungi increased with increasing elevation (r = 0.522, P = 0.009 and r = 0.557, P = 0.005, respectively; Figure ).

Relationship between fungal Chao1 (A) value and richness (B) and elevation.

A total of 4493 unique OTUs were identified and assigned to six phyla across all soil samples at the six elevations. The Ascomycota, Zygomycota, and Basidiomycota were the dominant phyla, representing 52.34, 32.18, and 15.05% of the sequences, respectively, while relative abundance of the Chytridiomycota and Glomeromycota ranged from 0.05 to 0.17% and from 0.003 to 0.02%, respectively (Figure and Supplementary Table ). Nine classes of fungi accounted for 91.97% of the total sequences, including Leotiomycetes, Dothideomycetes, Eurotiomycetes, Sordariomycetes, n class="Species">Ascomycota Incertae sedis, Lecanoromycetes, Pezizomycetes, Agaricomycetes, and Zygomycota Incertae_sedis (Supplementary Table ). We allocated OTUs to 14 functional groups, among which, the most dominant were the ectomycorrhiza, plant pathogens, and mycoparasites with mean relative abundances >1% (Supplementary Table ).

Relative abundances of the dominant fungal phyla and classes in soils at different elevations. The phylum and classes with relative abundances of <1% were assigned to “other.”

We found that while there was a difference in fungal community composition between 2500 m and the other elevations, there were no differences between 2000 and 2100 m, 2100 and 2200 m, 2200 and 2300 m, and 2300 and 2400 m using the ANOSIM test (Table ). The first two axes of the PCoA explained 37.38% of the observed variation in community composition, and the biplot indicated that fungal communities tended to be more similar within than between elevations (Supplementary Figure ). The significant difference in community composition between different elevations was also shown by PERMANOVA (R2 = 0.125, P = 0.001).

Correlations Between Soil Variables and Fungal Diversity and Community Composition

Fungal species richness was negatively correlated with C/N ratio (r = -0.560, P = 0.004) (Figure ) and TC (r = -0.410, P = 0.047), but was not associated with the other soil parameters (Supplementary Table ). The distance-based multivariate showed that elevation, C/N ratio, n class="Chemical">TC, TN, and pH contributed to the variation in fungal community composition individually (Table ). Among the soil physicochemical variables, C/N ratio provided the greatest explanatory power of the variation in community composition (11.2%) and it was selected as one of best environmental drivers in the multivariate model (Table ). RDA analysis explained 36.35% of the total variation in soil fungal community composition, and the first two components explained 28.65% of the variation (Supplementary Figure ). Soil C/N ratio was the strongest environmental driver of fungal community composition except elevation (Table and Supplementary Figure ). The NMDS plot also exhibited the best separation of soil samples along the C/N ratio gradient (Figure ).

Effects of C/N ratio on the alpha and beta diversity of soil fungi. Alpha diversity index (observed species) decreased with increasing C/n class="Chemical">N ratio (A). NMDS plot showing the best separation of samples along the C/N ratio gradients (B).

We found that the relative abundance of some fungal taxonomic groups, such as Basidiomycota, Ascomycota Incertae sedis, and Leotiomycetes, was significantly negatively correlated with soil C/N ratio (Supplementary Tables ). We also found a negative correlation between the relative abundance of Leotiomycetes and content of soil n class="Chemical">TN (r = -0.524, P = 0.009) and TC (r = -0.621, P = 0.001) (Supplementary Table ). Among the functional guilds of fungi, the relative abundance of mycoparasites and foliar epiphytes were negatively correlated with soil C/N ratio (r = -0.405, P = 0.049 and r = -0.427, P = 0.037, respectively), and the relative abundance of foliar epiphytes was correlated with other soil properties, such as DOC, DON, TN, and TC (Supplementary Table ).

Discussion

Our results show that elevation strongly influenced the diversity of fungal communities in the tundra soils of Changbai Mountain, where diversity increased linearly with increasing elevation (Figure ). Previous studies have observed different elevational patterns in soil fungal diversity. For example, Ping et al. (2017) found that soil fungal diversity decreased between the elevations from 699 to 937 m but increased between 937 and 1044 m in Korean pine forests on Changbai Mountain, showing a hollow curve’s pattern. Yang et al. (2017c) observed a monotonic decrease in soil fungal diversity along elevation from 700 to 2600 m by the culture-dependent method on Changbai Mountain. Our previous study showed that there was a lack of elevational pattern in soil fungal diversity at an altitudinal gradient from 530 to 2200 m on Changbai Mountain (Shen et al., 2014). The larger scale in elevational gradient contains different vegetation types, which is assumed to impact the elevational pattern in soil fungal diversity (He et al., 2017; Yang et al., 2017a). In addition, different microbial community analysis methods may also influence the observations on fungal distribution patterns along elevation.

To the best of our knowledge, our study is the first observation of a significant linear increase in soil fungal diversity along elevation in alpine tundra ecosystem. Similarly, a linear increase in diversity of foliar n class="Species">fungal endophytes along the increasing elevational gradient was observed within an Ermans birch forest on Changbai Mountain, which was strongly driven by leaf carbon content (Yang et al., 2016). Some environmental factors that differ among elevations may account for the diversity pattern along elevation. For example, soil pH has been found to drive the diversity pattern of soil fungi at the large-scale elevational gradient on Changbai Mountain (Shen et al., 2014), and climatic variables (temperature and precipitation) were assumed to be the potential drivers on the diversity pattern of ectomycorrhizal fungi along the elevational gradient in Hyrcanian forests of northern Iran (Bahram et al., 2012).

In this study, soil C/N ratio was significantly correlated with elevation and decreased with increasing elevations (r = -0.558, P = 0.005). Within the alpine tundra of Changbai Mountain, the plant communities at the lower elevations were dominated by dwarf shrubs (e.g., n class="Species">Rhododendron aureum and Vaccinium uliginosum), whereas the communities at the higher elevations were dominated by herbs/sedges, such as Sanguisorba stipulata, Sanguisorba parviflora, and Carex atrata (Wei et al., 2004; Zhang et al., 2017). Different dominant plant species at different elevations within alpine tundra may determine the variations of soil C/N ratio at the small-scale elevational gradient, which might be regulated by the components of root exudates and plant residues (Quideau et al., 2001; Knelman et al., 2012; Waring et al., 2015). In our study, soil fungal diversity significantly decreased with the increasing C/N ratio, and soil fungal community composition was mostly correlated with the C/N ratio among all soil variables. The above results indicated that C/N ratio was a key factor constraining both diversity pattern and community composition along elevation within alpine tundra on Changbai Mountain.

Similarly, the effect of C/N ratio on fungal diversity was also reported in other high and cold environments, such as the maritime Antarctic (n class="Chemical">Newsham et al., 2016) and the alpine grasslands of the Tibetan Plateau (Yang et al., 2017a). In addition, the effect of C/N ratio on soil fungal community composition was found in arctic tundra (Shi et al., 2015), permafrost-affected soils on the Tibetan Plateau (Zhang et al., 2014) and Korean pine forests on Changbai Mountain (Ping et al., 2017). Soil C/N ratio is generally thought as a proxy for nutrient availability (Cleveland and Liptzin, 2007) and as such, the change in ambient C/N ratio can substantially influence the fungal anabolism and foraging strategies (Prevost-Boure et al., 2011; Drake et al., 2013; Grosso et al., 2016). In terms of stoichiometry, the C/N ratio of a typical fungal biomass in soils was at 10–15 (Said et al., 2014), while the highest soil C/N ratio reached 18.44 in this study. High C/N ratio may broke the stoichiometric balance between soil and mycelia, restraining the activity of exoenzymes and accumulation of fungal biomass (Sinsabaugh et al., 2013). Of note, soil C/N ratio co-varied with the community composition of aboveground plants. On the one hand, rooting architectures and root exudates were quite different between dwarf shrubs and sedges/herbs, which dominated different elevations in this study. On the other hand, different dominant species at different elevations may select different fungal partners, given that mycorrhizal fungi is generally host specificity (Ishida et al., 2007; Tedersoo et al., 2008). Consequently, the observed pattern driven by C/N ratio in this study may be partly attributed to the variation of aboveground plant communities. Previously, many studies reported that plant species identity strongly affected soil fungal distributions (Tedersoo et al., 2016; Dassen et al., 2017; Gao et al., 2017). Nonetheless, the relative effects of plant dominant species and soil C/N ratio on soil fungal communities within alpine tundra still warrant further investigations.

In this study, the phylum of Zygomycota was more dominant than Basidiomycota, which was different from the observations in other tundra studies (Nemergut et al., 2005; Timling et al., 2014). The fungi in Zygomycota are generally regarded as the degraders of sucrose and n class="Chemical">cellulose, while most fungi in Basidiomycota are thought as lignin-degrading fungi (Waksman et al., 1928). We also found that ectomycorrhizal fungi, plant pathogens, and mycoparasites were the dominant functional guilds. The dominant plant species (e.g., dwarf shrubs and sedges) in alpine tundra were the host of ectomycorrhizal fungi (Wang and Qiu, 2006), which provided lots of living space for ectomycorrhizal fungi. The ectomycorrhizal fungi and plant pathogens have also been found as the dominant functional guilds in the arctic (Timling et al., 2014). Here, we also found that different taxonomic groups and functional guilds responded to different environmental variables (Supplementary Tables ). For example, the relative abundance of Ascomycota was significantly correlated with DON, while the relative abundance of Basidiomycota was significantly correlated with soil C/N ratio. The relative abundance of ectomycorrhizal fungi was mostly related with NH4+-N, while the mycoparasite and foliar epiphytes were mostly related with soil C/N ratio. In a field experiment, Nguyen et al. (2016b) found that the compositions of ectomycorrhizal fungi and saprotrophs were affected by host plant phylogenetic diversity and plant community structure, respectively. Vasutova et al. (2017) observed tree height was an important driver for the mycorrhizal fungi, whereas, the community composition of saprotrophs were mainly shaped by vegetation and edaphic variables. Solly et al. (2017) revealed that the relative abundance of ectomycorrhizal fungi was positively correlated with NH4+ concentrations but negatively correlated with soil C/N ratio in an alpine treeline. The responses of ectomycorrhizal fungal communities to some environmental factors under warming conditions may lead to major shifts in dominant taxa and mycelia production (Peter et al., 2001; Lilleskov et al., 2002; Solly et al., 2017). Therefore, in the context of climate changing, the studies on the shifts of functional guilds in alpine tundra should be more cautious.

Conclusion

In summary, soil fungal community composition and diversity changed along a small-scale elevational gradient in an alpine tundra on Changbai Mountain. Soil C/N ratio was not only negatively correlated with fungal diversity, but also significantly correlated with fungal community composition, and the relative abundances of specific phyla, class, and functional guilds. The results suggest that soil C/N ratio is a key factor for determining fungal elevational distribution pattern in n class="Species">alpine tundra ecosystems.

Author Contributions

HC designed the study. CS collected soil samples. YN and n class="Chemical">CS conducted the experiment. YN, TY, and KZ analyzed the data. YN, TY, and HC wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Table 1

Dissimilarities in fungal OTU community composition between elevations on Changbai Mountain as determined by ANOSIM.

Elevation (m)	2100	2200	2300	2400	2500
2000	0.29	0.11	0.72	0.94	1.00
2100		0.23	0.91	0.69	1.00
2200			0.2	0.3	0.91
2300				0.2	1.00
2400					0.50

Table 2

Distance-based multivariate linear model analysis (DistLM) for microbial fungal community composition.

Variable	% Var	Pseudo-F	P	Cum (%)
Variables individually
Elevation	20.8	5.7784	0.001
C/N	11.2	2.7701	0.003
TC	9.3	2.2579	0.010
ph	7.7	1.8408	0.023
TN	8.0	1.9103	0.033
DON	7.2	1.7024	0.053
NO3--N	6.2	1.4515	0.099
DOC	6.4	1.4958	0.108
NH4+-N	4.8	1.1016	0.315
Moisture	2.6	0.5840	0.932
Variables fitted sequentially
Elevation	20.8	5.7784	0.001	20.8
C/N	5.7	1.6325	0.031	26.5
DON	5.7	1.6898	0.036	32.2
NH4+-N	5.2	1.5937	0.051	37.4
NO3--N	4.5	1.4007	0.088	41.9
TN	4.0	1.2694	0.175	45.9
DOC	3.7	1.1803	0.246	49.6
Moisture	3.4	1.0926	0.334	53.0
pH	2.8	0.8813	0.630	55.8
TC	1.9	0.5787	0.901	57.7