ORIGINAL RESEARCH article
Dramatic Increases of Soil Microbial Functional Gene Diversity at the Treeline Ecotone of Changbai Mountain
- 1State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing, China
- 2State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China
- 3University of the Chinese Academy of Sciences, Beijing, China
- 4CAS Key Laboratory of Environmental Biotechnology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing, China
- 5Institute for Environmental Genomics and Department of Microbiology and Plant Biology, University of Oklahoma, Norman, OK, USA
- 6State Key Laboratory of Environment Simulation and Pollution Control, School of Environment, Tsinghua University, Beijing, China
- 7Earth and Environmental Sciences, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
The elevational and latitudinal diversity patterns of microbial taxa have attracted great attention in the past decade. Recently, the distribution of functional attributes has been in the spotlight. Here, we report a study profiling soil microbial communities along an elevation gradient (500–2200 m) on Changbai Mountain. Using a comprehensive functional gene microarray (GeoChip 5.0), we found that microbial functional gene richness exhibited a dramatic increase at the treeline ecotone, but the bacterial taxonomic and phylogenetic diversity based on 16S rRNA gene sequencing did not exhibit such a similar trend. However, the β-diversity (compositional dissimilarity among sites) pattern for both bacterial taxa and functional genes was similar, showing significant elevational distance-decay patterns which presented increased dissimilarity with elevation. The bacterial taxonomic diversity/structure was strongly influenced by soil pH, while the functional gene diversity/structure was significantly correlated with soil dissolved organic carbon (DOC). This finding highlights that soil DOC may be a good predictor in determining the elevational distribution of microbial functional genes. The finding of significant shifts in functional gene diversity at the treeline ecotone could also provide valuable information for predicting the responses of microbial functions to climate change.
With the rapid development of molecular tools, great progress on microbial biogeography has been made over the past decade (Martiny et al., 2006; Green et al., 2008; Hanson et al., 2012; Zhou et al., 2016). Researchers examining biogeographic patterns are now focusing on functional traits rather than just taxa (Fierer et al., 2012; Hanson et al., 2012; Fierer et al., 2013). Patterns of functional traits provide a better understanding of ecosystem functions, biogeochemical cycles and microbial responses to environmental change, which can shed light on fundamental questions in biology (Petchey and Gaston, 2002; Green et al., 2008; Fuhrman, 2009; Reiss et al., 2009). Recent studies have focused on functional gene distribution in soil microbial communities across large-scale natural ecosystems or latitudinal gradients (Fierer et al., 2012, 2013; Shi et al., 2015). There were also studies about the functional gene distribution along elevation, but the range and scale of elevational gradients were relatively small (Yang et al., 2014; Ding et al., 2015).
The high altitude limit of forests, commonly referred to as the treeline, timberline or forest line, has obtained extensive attention due to its sensitivity to climate change (Körner and Paulsen, 2004; Holtmeier and Broll, 2007). In reality, the transition from the uppermost closed montane forests to the treeless alpine vegetation is not really a line, but rather a steep gradient of increasing tree stands fragmentation and stunted growth, often called the treeline ecotone (Körner and Paulsen, 2004). However, little research has focused on the soil microbial communities at the treeline in elevational studies, until recently when Thébault et al. (2014) examined these communities with phospholipid fatty acid analysis (PLFA) technique and observed close correlations between microbial taxonomic communities and nutrient availability. Ding et al. (2015) examined the functional gene structure at the forest timberline using GeoChip technology, and they found that functional gene compositions were significantly different between shrubland and coniferous forest. Nonetheless, we do not know if microbial elevational diversity patterns will change in treeline ecotone and if microbial taxa and functional genes show the same pattern along an elevation gradient.
The Changbai Mountain is the highest mountain in Northeast China and is one of the most well-protected and conserved natural ecosystems on Earth (He et al., 2005). The elevational distribution of vegetation along the mountainside mirrors the latitudinal vegetation gradient from temperate to frigid zones on the Eurasian continent (Xu et al., 2004; Zhang et al., 2011). The treeline ecotone is at 1900–2000 m at the transition from birch forest to tundra. In this study, we hypothesized that soil microbial community would show different elevational diversity patterns due to different environmental controls. We used GeoChip 5.0, a functional gene microarray, and 16S rRNA gene sequencing to examine microbial functional gene and taxonomic diversity along an elevation gradient on Changbai Mountain. We found that microbial functional gene richness exhibited a dramatic increase at the treeline ecotone, but the bacterial taxonomic and phylogenetic diversity did not exhibit such a similar trend. Our results indicated that the elevational α-diversity patterns for bacterial taxa and functional genes were quite different, but that β-diversity patterns were similar. Our results also indicated that soil DOC might be a good predictor of microbial functional gene elevational distribution. This work provides valuable information for predicting the response of microbial functions to climate change.
Materials and Methods
Site Selection and Soil Sampling
As the highest mountain in north-eastern China and the head of three large rivers (the Songhua, Yalu and Tumen), Changbai Mountain (126°55′–129°00′E; 41°23′–42°36′N) is located in Jilin Province, which extends along the border of China and North Korea. Topographic and climatic variations result in a distinct vertical zonation of major forest types, especially along the northern slope. It has a typical continental temperate monsoon climate. Along the elevational gradient from 530 to 2200 m, the mean annual temperature decreases from 2.9 to -4.8°C, and the mean annual precipitation increases from 632 to 1154 mm (Tong et al., 2010). Below 1100 m there is a typical temperate forest, mainly composed of Korean pine and hardwood species; from 1100 to 1700 m, evergreen coniferous forest, dominated by spruce and fir, is typical; from 1700 to 2000 m is mountain birch forest; and above 2000 m is unique alpine tundra, which marks the southernmost occurrence of this ecosystem type on the eastern Eurasian continent. The forest-tundra treeline ecotone is present from 1900 to 2000 m (Supplementary Figure S1). The main climatic and ecological characteristics along the elevational gradient are summarized in Supplementary Table S1.
We collected soil samples from the northern slope of Changbai Mountain on 30 August 2013 at seven elevations (500, 700, 1000, 1300, 1600, 1900, and 2200 m). At each elevation, soil samples were collected from 5 plots (10 m × 10 m) as five independent replicates. In each plot, samples of the soil organic layer (∼10 cm × 10 cm in area, and 0–5 cm depth directly below the litter layer) were collected at ten random points using a sterile blade and composited together as a single sample. Visible roots and residue were removed from each sample prior to homogenizing the soil fraction. The fresh soil samples were sieved through 2-mm meshes and subdivided into two subsamples. One was stored at 4°C to determine the physical and chemical properties, and the other was stored at -20°C for subsequent DNA extraction.
Soil Physicochemical Analyses
Soil pH was measured using a pH meter (FE20-FiveEasyTM pH, Mettler Toledo, Germany) after shaking a soil water (1:5 w/v) suspension in the shaker for 30 min. Soil moisture was measured gravimetrically. Total carbon (TC) and total nitrogen (TN) contents were measured by elemental analyzer (Vario MAX, Elementar, Germany). Ammonium (NH4+-N), nitrate (-N), dissolved organic carbon (DOC) and dissolved total nitrogen (DTN) were extracted at a ratio of 10 g fresh soil to 100 mL 2 M KCl. After shaking for 1 h, -N, -N and DTN in the extracts (filtering through 10 cm diameter quantitative filter paper) were analyzed using a continuous flow analytical system (San++ System, Skalar, Holland), and DOC was determined using a TOC analyzer (Multi N/C 3000, Analytik Jena, Germany). Dissolved organic nitrogen (DON) was calculated as follows: DON = DTN --N --N. Available P (AP) in soil was extracted by sodium bicarbonate and determined using the molybdenum blue method (Olsen et al., 1954). Available K (AK) in soil was extracted by ammonium acetate and determined by flame photometry (Carson, 1980).
Soil DNA was extracted using a PowerSoil DNA kit (MO BIO, Carlsbad, CA, USA), and then purified with an UltraClean 15 DNA purification kit (MO BIO, Carlsbad, CA, USA). DNA concentrations were estimated by electrophoresis on 1% agarose gels, and DNA was diluted to 20 ng/μl before use in PCR.
Illumina MiSeq Sequencing and Data Processing
16S rRNA genes were amplified using a common primer set (515F, 5′-GTGCCAGCMGCCGCGGTAA-3′; 806R, 5′-GGACTACHVGGGTWTCTAAT-3′) combined with adapter sequences and barcode sequences (Caporaso et al., 2011). PCR amplifications were conducted with 25 μl 2× premix (TaKaRa, Otsu, Japan), 0.5 μl 20 mM each forward and reverse primer, and 50 ng of DNA, and the volume was completed to 50 μl with double-distilled water. Each sample was amplified in triplicate using 30 cycles (94°C for 30 s, 55°C for 30 s, and 72°C for 30 s) with a final extension at 72°C for 10 min. The three reaction products were pooled and purified using a QIAquick PCR purification kit (QIAGEN, Shenzhen, China) and then quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, NC, USA). PCR products were sequenced on two lanes of a 2 × 151 bp sequencing run on an Illumina MiSeq (Caporaso et al., 2012).
Raw data were processed and analyzed as previously described (Caporaso et al., 2012), using the QIIME software and following the workflow1. Briefly, sequences were quality filtered (max value of 0.5) and clustered into 97% similar phylotypes after removing singleton sequences. Taxonomy was identified using the Ribosomal Database Project classifier (Wang et al., 2007) which was trained on the Greengenes 13_8 16S rRNA database (McDonald et al., 2012). To rarify all the data sets to the same level of sampling effort, 4000 sequences were randomly selected. Sequences were deposited to the MG-RAST metagenomics analysis server2 and are available to the public (accession numbers from 4706639.3 to 4706673.3).
GeoChip Hybridization and Data Processing
DNA was extracted using a PowerSoil DNA kit as described above. GeoChip 5.0 was used to analyze DNA samples as described previously (Wang et al., 2014). Briefly, DNA (1 μg) was labeled with the fluorescent dye Cy-3 using a random priming method and then purified with the QIA quick purification kit (QIAGEN, Shenzhen, China) according to the manufacturer’s instructions. After purification, the labeled DNAs were hybridized with the Agilent platform-based GeoChip 5.0 arrays at 67°C plus 10% formamide for 24 h. GeoChip microarrays were scanned with a 100% laser power and 100% photomultiplier tube with a NimbleGen MS 200 Microarray Scanner (Roche, Basel, Switzerland). Spots with a signal-to-noise ratio <2.0, signals <150, or <1.3 times the background were discarded prior to statistical analyses. Raw data were quantified and processed using the analysis pipeline as previously described (He et al., 2010; Tu et al., 2014). Then processed GeoChip data were analyzed using the following steps: (i) removing genes detected in less than 3 of 5 samples from the same elevation; (ii) normalizing the signal intensity of each spot by dividing the signal intensity by the total intensity of the microarray followed by multiplying by a constant; and (iii) transforming data to the natural logarithm.
Bray–Curtis similarity distance based on bacterial OTU table and functional gene table was used to calculate bacterial dissimilarities and microbial functional gene dissimilarities. Principal co-ordinates analysis (PCoA) were performed on the basis of Bray–Curtis dissimilarities. The variations in beta diversity were plotted against changes in elevational distance to discern any patterns. The distance–decay relationship (which measures how similarity decays with increasing distance between pairwise sites) was analyzed using a Gaussian generalized linear model, and the significance was determined using Pearson’s correlation. In order to identify environmental and biogeochemical factors associated with functional gene diversity, correlations between diversity metrics and GeoChip data were conducted by SPSS software (version 20.0). Mantel tests of Bray-Curtis similarity distance values were calculated and carried out within vegan package in R (R Development Core Team, 2010; Dixon, 2003). The multivariate regression tree was created within the mvpart package in R (De’Ath, 2002). Distance-based linear model multivariate analysis (DistLM) (McArdle and Anderson, 2001) was implemented to determine the influence of environmental variables on microbial functional gene composition. Marginal tests were applied to assess the statistical significance and percentage contribution of each variable in isolation, and sequential tests were performed to evaluate the cumulative effect of environmental variables.
General Characteristics of Soil Microbial Communities
We obtained 329,598 high quality sequences from 16S rDNA amplicon sequencing for all soil samples (4,321 to 20,159 sequences per sample). A total of 30,302 unique OTUs were identified and were assigned to 43 bacterial phyla. Among the identified phyla, Proteobacteria (23.1%), Acidobacteria (20.8%), and Verrucomicrobia (17.29%) were the most abundant phyla across the seven elevations (Supplementary Table S4). Other dominant phyla across all the Changbai Mountain soils were Actinobacteria (5.92%), Bacteroidetes (5.76%), Chloroflexi (4.27), Planctomycetes (4.47%), Firmicutes (7.21%), accounting for more than 27% of the bacterial sequences from each of the soils (Supplementary Figure S3). In addition, Nitrospirae, Gemmatimonadetes, Armatimonadetes and other 29 phyla were present in most soils but at relatively low abundances (Supplementary Table S4).
A total of 55,874 functional genes in 35 samples were detected by GeoChip 5.0. These genes were mainly involved in six functional categories: carbon (C)-cycling (15.9%), nitrogen (N)-cycling (4.2%), environmental stress response (15.3%), metal homeostasis (26.2%), organic remediation (8.9%) and virulence (16.2%). Other functional categories involved in sulfur and phosphorus cycling, and electron transfer accounted for a small proportion (Supplementary Table S2). Functional genes shared among all elevations accounted for 15.57% of the detected genes, while those shared among 2–6 elevations accounted for 83.47%. The greatest number of unique functional genes (functional genes found exclusively in one elevation) were detected at 2200 m (456 functional genes), suggesting that the tundra ecosystem has more unique functional genes than forest ecosystems (Figure 1). Specific genes involved in carbon and nitrogen cycling like nifH, nosZ, ureC, CODH genes, mcrA, and pmoA, all show a higher gene richness at 1900 and 2200 m (Supplementary Table S3). Response ratio analysis showed that specific functional gene categories that exhibited significant (95% CI) higher gene abundances from coniferous forest to birch forest at the treeline ecotone are mainly involved in carbon degradation and fixation (Figure 2).
FIGURE 1. A network map showing the interactions of microbial functional genes among all the samples from different elevations. Each point represents one independent functional gene. Functional genes in the left column were unique to one site, while those in the right belonged to multiple elevation sites.
FIGURE 2. Specific functional gene categories that exhibited significant changes in abundance from coniferous forest (1600 m) to bich forest at the treeline ecotone (1900 m). Significance was determined using the response ratio analysis at a 95% confidence interval (CI).
Alpha and Beta Diversity Patterns
At the taxonomic and phylogenetic level, the alpha diversity of bacterial communities (randomly selected 4,000 sequences per sample) did not show an apparent trend along the elevation gradient (Figure 3), which is consistent with our previous results (454 data in Shen et al., 2013). The number of functional genes detected was significantly (p < 0.05) higher at 1900 and 2200 m than at lower elevations, resulting in a sharp increase in functional gene richness at the treeline ecotone (Figure 3). These results suggest that the treeline ecotone has a strong influence on microbial functional gene richness, rather than bacterial diversity.
FIGURE 3. The elevational alpha diversity patterns for bacterial communities and microbial functional genes. Bacterial alpha diversity (OTU richness) exhibited no apparent elevational gradient, whereas microbial functional gene alpha diversity (functional gene richness) showed a sharp increase at 1900 m (treeline ecotone). Different letters indicate significant differences (p < 0.05).
Principal co-ordinates analysis analysis based on 16S rRNA gene sequencing and GeoChip datasets showed that both bacterial and functional gene composition tended to be relatively similar among samples within the same elevation and distinctly different among the different elevations (Figure 4). Dissimilarity tests revealed that bacterial composition and functional gene composition differed significantly (p < 0.05) among elevations (Supplementary Table S6; Table 1). Furthermore, both the bacterial (r = 0.409, P < 0.001) and microbial functional gene (r = 0.708, P < 0.001) compositional dissimilarities significantly increased with increasing elevations, which showed significant elevational distance-decay patterns (Figure 5). It is noteworthy that functional gene data separated into two groups based on which side of the treeline ecotone the samples were collected (Figure 5). This suggests that effect of treeline ecotone on functional gene composition is bigger than that on bacterial composition. Nevertheless, these results suggest that the elevational alpha diversity patterns for bacterial taxa and functional genes are different, but they could exhibit similar elevational beta diversity patterns.
FIGURE 4. Principal co-ordinates analysis (PCoA) of the bacterial community composition (A) and microbial functional gene composition (B).
TABLE 1. Results from analysis of microbial functional gene composition dissimilarities (MRPP, ADONIS, and ANOSIM) among different elevations.
FIGURE 5. The relationships between the community dissimilarity (bacterial communities and microbial functional genes) and elevational change. Both the bacterial communities and microbial functional genes showed significant elevational distance-decay patterns (r = 0.409, p < 0.001; r = 0.708, p < 0.001).
The Key Environmental Factors Influencing the Diversity and Composition of Soil Bacterial Community and Microbial Functional Genes
In terms of OTU richness and Faith’s PD, bacterial OTU richness was significantly (p < 0.05) correlated with soil pH, moisture, -N, TC and TN, while bacterial phylogenetic diversity was significantly correlated with soil pH and AK (Supplementary Table S5). Mantel tests showed that pH had the highest correlation with bacterial community composition (r = 0.36, P < 0.001, Table 2). Specifically, soil pH showed a significant correlation with the relative abundance of the four dominant phyla (Supplementary Figure S4).
TABLE 2. Mantel test results for the correlation between bacterial community composition and microbial functional gene composition and environmental variables along the elevational gradient.
In terms of functional gene richness and Shannon index, soil DOC, DON, TC, moisture, AP and -N showed significant (p < 0.05) correlations with functional gene diversity (Supplementary Table S5). Mantel tests showed that the microbial functional gene composition was significantly (p < 0.05) correlated with soil DOC, DON and TC (Table 2). Of all of the environmental variables examined, DOC showed the highest correlation with the functional structure based on Multivariate Regression Trees (MRT) analyses (Figure 6). Sequential tests of DistLM analysis indicated four significant variables (DOC, AP, AK, -N) that explained 62.3% of the total variation of microbial functional gene composition, with soil DOC providing the greatest explanatory power (43.5% of the total variation; Table 3). These results suggest that the bacterial community composition and microbial functional gene diversity were influenced by different drivers with the bacterial diversity and community composition being strongly influenced by soil pH, while the microbial functional gene diversity and composition were significantly correlated with soil DOC.
FIGURE 6. Multivariate regression tree (MRT) showing the most important factor with microbial functional gene composition. DOC, dissolved organic carbon; NH4, ammonium; P, available P.
TABLE 3. Results of distance-based multivariate linear model (DistLM) for microbial functional gene composition showing the % variation explained by environmental variables.
The elevational diversity gradient is one of the most fundamental patterns in biogeography (Lomolino, 2001; Bryant et al., 2008). There are many lines of evidence showing that microbes do not follow the classic decreasing or unimodal elevational diversity patterns that plants and animals follow. As expected, we observed that soil bacteria exhibited no apparent elevational trend based on taxonomic diversity, which is consistent with previous reports (Fierer et al., 2011; Shen et al., 2013; Xu et al., 2014; Yuan et al., 2014). In contrast, microbial functional gene richness showed a dramatic shift at the treeline ecotone along this elevation gradient. To the best of our knowledge, this finding has not been reported in other elevational studies, suggesting that the primary controls on microbial functional gene diversity may be different from those for taxa. Moreover, microbial functional gene composition significantly differed across elevations, which is in line with a previous study of the Tibetan Plateau (Yang et al., 2014). Given that many researchers have found a clear bacterial or micro-eukaryotic compositional difference among different elevations (Singh et al., 2012; Shen et al., 2014, 2015; Yuan et al., 2014; Jarvis et al., 2015), these results suggest that elevation could be a major driver of variation in both microbial taxonomic and functional gene composition. Beta diversity (compositional dissimilarity among sites) is central to the ecological and evolutionary processes shaping the distribution of species (Condit et al., 2002). Yet elevational patterns of beta diversity have not received as much attention as that of alpha diversity (Wang et al., 2012). Here we found that bacterial community dissimilarities significantly increased with increasing elevation, and that microbial functional genes also exhibited a significant elevational distance-decay pattern, though the treeline ecotone had a bigger influence on functional gene composition than bacterial taxa composition. However, this elevational distance-decay pattern may not be observed in all environments. Bryant et al. (2008) observed that the compositional similarity of soil Acidobacteria significantly decreased with elevational distance, whereas Wang et al. (2012) found that benthic bacteria did not show a significant elevational distance-decay relationship. The elevational distance-decay pattern may be more easily observed in soils than aquatic environments or it could be caused by the greater spatial heterogeneity in soils than that in water. Together, these results refute our hypothesis, and suggest that the elevational alpha diversity patterns for bacterial communities and microbial functional genes are different, although they could have similar elevational beta diversity patterns.
The relationship between taxonomic diversity and functional diversity in soils is largely unknown (Petchey and Gaston, 2002; Rosenfeld, 2002; Allison and Martiny, 2008). Scientists are trying to test this relationship at different scales by combining 16S rRNA gene sequencing with whole genome shotgun sequencing or GeoChip analysis (Fierer et al., 2012, 2013; Ding et al., 2015; Leff et al., 2015; Shi et al., 2015). For example, Fierer et al. (2012) found significant correlations between bacterial diversity and functional gene diversity (include alpha and beta) in different ecosystems across the globe. Furthermore, this correlation seems to become stronger within the single tallgrass prairie ecosystem in the midwestern United States (Fierer et al., 2013). The above-mentioned findings highlight that the overall diversity of functional gene categories found in a given sample is, to some degree, predictable from the taxonomic diversity of the microbial communities. Actually, while the taxonomic and phylogenetic beta diversity is significantly correlated with functional gene beta diversity, the alpha diversity may not show identical manner. For instance, Shi et al. (2015) found no correlation between bacterial richness and microbial functional gene richness across an Arctic tundra ecosystem. Similarly, our study here together with a study in a forest timberline, did not detect this relationship (Ding et al., 2015). The contradictory results could be caused by differences in terms of technique (shotgun metagenome sequencing vs. GeoChip, see Zhou et al., 2015), scale or ecosystem examined. Therefore, more research using multiple comparable techniques are needed to test this relationship (especially for alpha diversity) across different ecosystems. Functional redundancy of species is assumed to be a common feature in soils (Heemsbergen et al., 2004). This concept is based on the observation that some species perform similar roles in communities and ecosystems, and may therefore be substitutable with little impact on ecosystem processes (Lawton and Brown, 1993). Given the non-significant bacterial richness pattern observed along the elevation gradient and the sharp increase of functional gene richness at the treeline ecotone observed in this study, soil bacterial taxa appear to exhibit a high degree of functional redundancy. Of course, functional redundancy cannot be determined simply based on the relationship of taxonomic and functional gene alpha diversity, as redundancy is difficult to establish because it requires detailed knowledge of the microbial populations that perform a specific process (Allison and Martiny, 2008).
Identifying the environmental factors that best explain microbial community variation is a fundamental goal in microbial ecology (Martiny et al., 2006; Hanson et al., 2012). Microbial taxa display spatial patterns linked to geographic distance (Cho and Tiedje, 2000; Green and Bohannan, 2006; Martiny et al., 2006), vegetation type (Knelman et al., 2012; Shen et al., 2013), and soil characteristics (Fierer and Jackson, 2006; Lauber et al., 2009; Chu et al., 2010; Liu et al., 2014). Previous studies have demonstrated that soil pH drives both bacterial horizontal and elevational distribution (Fierer and Jackson, 2006; Chu et al., 2010; Griffiths et al., 2011; Shen et al., 2013). In this study, we found a significant relationship between soil pH and bacterial diversity, which is in line with our previous results based on pyrosequencing (Shen et al., 2013). However, no significant correlations were found here between pH and functional gene diversity and gene composition, suggesting a weak role of soil pH in influencing microbial functional genes.
The sensitivity and response of the treeline to climate change has been increasingly discussed from global to regional and smaller scales (Körner and Paulsen, 2004; Holtmeier and Broll, 2007). In this study, our results showed a dramatic shift of functional gene richness at the treeline ecotone. To our knowledge, this is the first reported observation of a significant shift of functional gene richness from a forest ecosystem to a tundra ecosystem along an elevation gradient. Climate warming leads to movement of the treeline toward a higher elevation, but little is known about the dynamics of microbial functions under this change. Thus, this finding provides a great value in predicting the response of microbial functions to climate change. We concluded this pattern was largely driven by soil DOC, since it was the most significant factor correlated with microbial functional gene diversity. Our study showed that soil DOC content significantly increased with elevation, particularly was higher at the treeline ecotone and tundra (Supplementary Figure S2). Variations in DOC content might caused by differences in vegetation types (see Supplementary Table S1). On one hand, litter input produced by different vegetations can lead to difference in soil organic matter composition (Quideau et al., 2001). On the other hand, vegetation type greatly influences soil physicochemical properties, which forms different soil microclimate (Knelman et al., 2012). DOC is a complex and heterogeneous mixture of C compounds, including humic acids, fulvic acids, hydrophobic neutrals, and hydrophilic compounds, which may be both a substrate for microbial activity and a byproduct of the subsequent microbial metabolic processes, thus providing strong evidence of links between microbial activity and DOC concentrations (Marschner and Kalbitz, 2003; van Hees et al., 2005; Bolan et al., 2011). For instance, previous studies on alpine ecosystems showed that tundra soils at higher elevations had more labile carbon and greater microbial activity than forest soils (Neff and Hooper, 2002; Sjögersten et al., 2003). Using the Biolog method, Tian et al. (2015) found that labile DOC accounted for the largest amount of variation in microbial metabolic functional diversity. Importantly, Straathof (2015) studied the relationship between composition of DOC and microbial respiration, and concluded that the quality of DOC, via regulation of microbial processes, might be an important indicator of soil functions. However, while above-mentioned studies mainly focused on microbial respiration and substrate utilization, we first found significant correlations between functional gene diversity and DOC. Meanwhile, specific functional gene categories that exhibited significant higher gene abundances at the treeline ecotone are mainly involved in carbon degradation and fixation. This suggests that these functional genes associated with C cycling are closely linked with soil C dynamics, which may result in significant differences in microbial metabolic potential between the two sides of treeline. In addition, it should be noted that other factors which were not detected in our study may also have fundamental contributions for microbial functional gene variations. For example, Ding et al. (2015) found soil temperature was the best predictor for microbial functional gene variations. At the same time, the authors noted that the contribution of soil DOC and other nutrient contents to the functional gene variations cannot be overlooked, since temperature is likely an indirect factor affecting the microbial community by controlling soil nutrient availability (Mitchell et al., 2012). Taken together, our results suggest that soil DOC could be a good predictor for microbial functional gene diversity and composition, and this finding provides strong evidence for DOC as an important indicator of soil functions at the functional gene level.
In summary, we observed the distribution of soil microbial functional gene along a complete elevational gradient on Changbai Mountain. Soil microbial functional gene richness exhibited a dramatic increase at the treeline ecotone, but bacterial diversity did not exhibit a similar elevational trend. Both bacterial taxa and functional genes showed significant elevational distance-decay patterns as their dissimilarities significantly increased with increased elevation. This suggests that the elevational α-diversity patterns for microbial taxonomic and functional genes are quite different, but that they may have similar elevational β-diversity patterns. While bacterial diversity/composition was strongly influenced by soil pH, microbial functional gene diversity/composition was significantly correlated with soil DOC. This finding highlights that soil DOC may be a good predictor of microbial functional gene elevational distribution. These results provide valuable information for our understanding of elevational diversity patterns and predicting the response of microbial functions to climate change. Coupling community and function is fundamental to advancing our understanding of ecology. We anticipate that our results will lead to more investigations of soil microbial communities from multiple viewpoints (taxonomic/phylogenetic/functional) to test existing or novel relationships and hypotheses.
HC conceived the idea. CS and YN conducted all the experiments for soil physicochemical analysis. CS collected soil samples. CS and YS conducted the data analysis. CS wrote the first draft and HC finalized the manuscript with assistance from all co-authors.
This work was supported by the National Natural Science Foundation of China (41371254), the National Program on Key Basic Research Project (2014CB954002), the Strategic Priority Research Program (XDB15010101, XDB15010302) of Chinese Academy of Sciences, Plan and Frontiers Projects of Institute of Soil Science (ISSASIP1641), the National Basic Research Program of China (2015FY110100), and the State Key Laboratory of Forest and Soil Ecology (LFSE2014-02). CS was also supported by the China Scholarship Council (CSC).
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.
We thank Wenju Liang, Huaibo Sun, Teng Yang, and Fang Liu for assistance with soil sampling. We also thank Jonathan Adams for language editing.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.01184
- ^ http://nbviewer.ipython.org/github/biocore/qiime/blob/1.9.1/examples/ipynb/illumina_overview_tutorial.ipynb.
- ^ http://metagenomics.anl.gov/
Bolan, N. S., Adriano, D. C., Kunhikrishnan, A., James, T., McDowell, R., and Senesi, N. (2011). Dissolved organic matter: biogeochemistry, dynamics, and environmental significance in soils. Adv. Agron. 110, 1–75. doi: 10.1016/B978-0-12-385531-2.00001-3
Bryant, J. A., Lamanna, C., Morlon, H., Kerkhoff, A. J., Enquist, B. J., and Green, J. L. (2008). Microbes on mountainsides, Contrasting elevational patterns of bacterial and plant diversity. Proc. Natl. Acad. Sci. U.S.A. 105, 11505–11511. doi: 10.1073/pnas.0801920105
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Huntley, J., Fierer, N., et al. (2012). Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. doi: 10.1038/ismej.2012.8
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U.S.A. 108, 4516–4522. doi: 10.1073/pnas.1000080107
Carson, P. L. (1980). “Recommended potassium test,” in Recommended Chemical Soil Test Procedures for the North Central Region. Bulletin 499, ed. W. C. Dahnke (Fargo: North Dakota Agricultural Experiment Station), 17–18.
Chu, H. Y., Fierer, N., Lauber, C. L., Caporaso, J. G., Knight, R., and Grogan, P. (2010). Soil bacterial diversity in the Arctic is not fundamentally different from that found in other biomes. Environ. Microbiol. 12, 2998–3006. doi: 10.1111/j.1462-2920.2010.02277.x
Ding, J. J., Zhang, Y. G., Deng, Y., Cong, J., Lu, H., Sun, X., et al. (2015). Integrated metagenomics and network analysis of soil microbial community of the forest timberline. Sci. Rep. 5, 7994. doi: 10.1038/srep07994
Fierer, N., Ladau, J., Clemente, J. C., Leff, J. W., Owens, S. M., Pollard, K. S., et al. (2013). Reconstructing the Microbial diversity and function of pre-agricultural tallgrass prairie soils in the United States. Science 342, 621–624. doi: 10.1126/science.1243768
Fierer, N., Leff, J. W., Adams, B. J., Nielsen, U. N., Bates, S. T., Lauber, C. L., et al. (2012). Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. Proc. Natl. Acad. Sci. U.S.A. 109, 21390–21395. doi: 10.1073/pnas.1215210110
Fierer, N., McCain, C. M., Meir, P., Zimmermann, M., Rapp, J. M., Silman, M. R., et al. (2011). Microbes do not follow the elevational diversity patterns of plants and animals. Ecology 92, 797–804. doi: 10.1890/10-1170.1
Griffiths, R. I., Thomson, B. C., James, P., Bell, T., Bailey, M., and Whiteley, A. S. (2011). The bacterial biogeography of British soils. Environ. Microbiol. 13, 1642–1654. doi: 10.1111/j.1462-2920.2011.02480.x
Hanson, C. A., Fuhrman, J. A., Horner-Devine, M. C., and Martiny, J. B. H. (2012). Beyond biogeographic patterns: processes shaping the microbial landscape. Nat. Rev. Microbiol. 10, 497–506. doi: 10.1038/nrmicro2795
He, H. S., Hao, Z. Q., Mladenoff, D. J., Shao, G. F., Hu, Y. M., and Chang, Y. (2005). Simulating forest ecosystem response to climate warming incorporating spatial effects in north-eastern China. J. Biogeogr. 32, 2043–2056. doi: 10.1111/j.1365-2699.2005.01353.x
He, Z. L., Deng, Y., Van Nostrand, J. D., Tu, Q. C., Xu, M. Y., Hemme, C. L., et al. (2010). GeoChip 3.0 as a high-throughput tool for analyzing microbial community composition, structure and functional activity. ISME J. 4, 1167–1179. doi: 10.1038/ismej.2010.46
Heemsbergen, D. A., Berg, M. P., Loreau, M., van Haj, J. R., Faber, J. H., and Verhoef, H. A. (2004). Biodiversity effects on soil processes explained by interspecific functional dissimilarity. Science 306, 1019–1020. doi: 10.1126/science.1101865
Jarvis, S. G., Woodward, S., and Taylor, A. F. S. (2015). Strong altitudinal partitioning in the distributions of ectomycorrhizal fungi along a short (300 m) elevation gradient. New Phytol. 206, 1145–1155. doi: 10.1111/nph.13315
Knelman, J. E., Legg, T. M., O’Neill, S. P., Washenberger, C. L., González, A., Cleveland, C. C., et al. (2012). Bacterial community structure and function change in association with colonizer plants during early primary succession in a glacier forefield. Soil Biol. Biochem. 46, 172–180. doi: 10.1016/j.soilbio.2011.12.001
Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75, 5111–5120. doi: 10.1128/AEM.00335-09
Leff, J. W., Jones, S. E., Prober, S. M., Barberan, A., Borer, E. T., Firn, J. L., et al. (2015). Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe. Proc. Natl. Acad. Sci. U.S.A. 112, 10967–10972. doi: 10.1073/pnas.1508382112
Liu, J. J., Sui, Y. Y., Yu, Z. H., Shi, Y., Chu, H. Y., Jin, J., et al. (2014). High throughput sequencing analysis of biogeographical distribution of bacterial communities in the black soils of northeast China. Soil Biol. Biochem. 70, 113–122. doi: 10.1016/j.soilbio.2013.12.014
Martiny, J. B. H., Bohannan, B. J. M., Brown, J. H., Colwell, R. K., Fuhrman, J. A., Green, J. L., et al. (2006). Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112. doi: 10.1038/nrmicro1341
McArdle, B. H., and Anderson, M. J. (2001). Fitting multivariate models to community data: a comment on distance-based redundancy analysis. Ecology 82, 290–297. doi: 10.1890/0012-9658(2001)082[0290:FMMTCD]2.0.CO;2
McDonald, D., Price, M. N., Goodrich, J., Nawrocki, E. P., DeSantis, T. Z., Probst, A., et al. (2012). An improved greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of Bacteria and Archaea. ISME J. 6, 610–618. doi: 10.1038/ismej.2011.139
Mitchell, R. J., Hester, A. J., Campbell, C. D., Chapman, S. J., Cameron, C. M., Hewison, R. L., et al. (2012). Explaining the variation in the soil microbial community: do vegetation composition and soil chemistry explain the same or different parts of the microbial variation? Plant Soil 351, 355–362. doi: 10.1007/s11104-011-0968-7
Olsen, S. R., Cole, C. V., Watanabe, F. S., and Dean, L. A. (1954). Estimation of Available Phosphorus in Soils by Extraction with Sodium Bicarbonate. Washington, DC: U.S. Department of Agriculture, 19.
Quideau, S. A., Chadwick, O. A., Benesi, A., Graham, R. C., and Anderson, M. A. (2001). A direct link between forest vegetation type and soil organic matter composition. Geoderma 104, 41–60. doi: 10.1016/S0016-7061(01)00055-6
Reiss, J., Bridle, J. R., Montoya, J. M., and Woodward, G. (2009). Emerging horizons in biodiversity and ecosystem functioning research. Trends Ecol. Evol. 24, 505–514. doi: 10.1016/j.tree.2009.03.018
Shen, C. C., Liang, W. J., Shi, Y., Lin, X. G., Zhang, H. Y., Wu, X., et al. (2014). Contrasting elevational diversity patterns between eukaryotic soil microbes and plants. Ecology 95, 3190–3202. doi: 10.1890/14-0310.1
Shen, C. C., Ni, Y. Y., Liang, W. J., Wang, J. J., and Chu, H. Y. (2015). Distinct soil bacterial communities along a small-scale elevational gradient in alpine tundra. Front. Microbiol. 6:582. doi: 10.3389/fmicb.2015.00582
Shen, C. C., Xiong, J. B., Zhang, H. Y., Feng, Y. Z., Lin, X. G., Li, X. Y., et al. (2013). Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem. 57, 204–211. doi: 10.1016/j.soilbio.2012.07.013
Shi, Y., Grogan, P., Sun, H. B., Xiong, J. B., Yang, Y. F., Zhou, J. Z., et al. (2015). Multi-scale variability analysis reveals the importance of spatial distance in shaping Arctic soil microbial functional communities. Soil Biol Biochem. 86, 126–134. doi: 10.1016/j.soilbio.2015.03.028
Singh, D., Takahashi, K., Kim, M., Chun, J., and Adams, J. M. (2012). A Hump-Backed trend in bacterial diversity with elevation on Mount Fuji. Jpn. Microb. Ecol. 63, 429–437. doi: 10.1007/s00248-011-9900-1
Sjögersten, S., Turner, B. L., Mahieu, N., Condron, L. M., and Wookey, P. A. (2003). Soil organic matter biochemistry and potential susceptibility to climatic change across the forest-tundra ecotone in the Fennoscandian mountains. Glob. Change Biol. 9, 759–772. doi: 10.1046/j.1365-2486.2003.00598.x
Thébault, A., Clement, J. C., Ibanez, S., Roy, J., Geremia, R. A., Perez, C. A., et al. (2014). Nitrogen limitation and microbial diversity at the treeline. Oikos 123, 729–740. doi: 10.1111/j.1600-0706.2013.00860.x
Tian, J., McCormack, L., Wang, J. Y., Guo, D., Wang, Q. F., Zhang, X. Y., et al. (2015). Linkages between the soil organic matter fractions and the microbial metabolic functional diversity within a broad-leaved Korean pine forest. Eur. J. Soil Biol. 66, 57–64. doi: 10.1016/j.ejsobi.2014.12.001
Tu, Q. C., Yu, H., He, Z. L., Deng, Y., Wu, L. Y., Van Nostrand, J. D., et al. (2014). GeoChip 4: a functional gene-array-based high-throughput environmental technology for microbial community analysis. Mol. Ecol. Resour. 14, 914–928. doi: 10.1111/1755-0998.12239
van Hees, P. A. W., Jones, D. L., Finlay, R., Godbold, D. L., and Lundström, U. S. (2005). The carbon we do not see—the impact of low molecular weight compounds on carbon dynamics and respiration in forest soils: a review. Soil Biol. Biochem. 37, 1–13. doi: 10.1016/j.soilbio.2004.06.010
Wang, C., Wang, X. B., Liu, D. W., Wu, H. H., Lu, X. T., Fang, Y. T., et al. (2014). Aridity threshold in controlling ecosystem nitrogen cycling in arid and semi-arid grasslands. Nat. Commun. 5:4799. doi: 10.1038/ncomms5799
Wang, J. J., Soininen, J., Zhang, Y., Wang, B. X., Yang, X. D., and Shen, J. (2012). Patterns of elevational beta diversity in micro- and macroorganisms. Glob. Ecol. Biogeogr. 21, 743–750. doi: 10.1111/j.1466-8238.2011.00718.x
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Xu, M., Li, X. L., Cai, X. B., Gai, J. P., Li, X. L., Christie, P., et al. (2014). Soil microbial community structure and activity along a montane elevational gradient on the Tibetan Plateau. Eur. J. Soil Biol. 64, 6–14. doi: 10.1016/j.ejsobi.2014.06.002
Yang, Y. F., Gao, Y., Wang, S. P., Xu, D. P., Yu, H., Wu, L. W., et al. (2014). The microbial gene diversity along an elevation gradient of the Tibetan grassland. ISME J. 8, 430–440. doi: 10.1038/ismej.2013.146
Yuan, Y. L., Si, G. C., Wang, J., Luo, T. X., and Zhang, G. X. (2014). Bacterial community in alpine grasslands along an altitudinal gradient on the Tibetan Plateau. FEMS Microbiol. Ecol. 87, 121–132. doi: 10.1111/1574-6941.12197
Zhang, M., Zhang, X. K., Liang, W. J., Jiang, Y., Dai, G. H., Wang, X. G., et al. (2011). Distribution of soil organic carbon fractions along the altitudinal gradient in Changbai Mountain, China. Pedosphere 21, 615–620. doi: 10.1016/S1002-0160(11)60163-X
Zhou, J. Z., Deng, Y., Shen, L. N., Wen, C. Q., Yan, Q. Y., Ning, D. L., et al. (2016). Temperature mediates continental-scale diversity of microbes in forest soils. Nat. Commun. 7:12083. doi: 10.1038/ncomms12083
Keywords: metagenomics, GeoChip, microbial functional genes, bacterial taxonomic and phylogenetic diversity, alpha and beta diversity patterns, treeline ecotone, elevation gradient, soil dissolved organic carbon
Citation: Shen C, Shi Y, Ni Y, Deng Y, Van Nostrand JD, He Z, Zhou J and Chu H (2016) Dramatic Increases of Soil Microbial Functional Gene Diversity at the Treeline Ecotone of Changbai Mountain. Front. Microbiol. 7:1184. doi: 10.3389/fmicb.2016.01184
Received: 16 June 2016; Accepted: 18 July 2016;
Published: 29 July 2016.
Edited by:Hongchen Jiang, Miami University, USA
Copyright © 2016 Shen, Shi, Ni, Deng, Van Nostrand, He, Zhou and Chu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Haiyan Chu, email@example.com
†The authors have contributed equally to this work.