Research article | Open | Published:
Geochip-based analysis of microbial communities in alpine meadow soils in the Qinghai-Tibetan plateau
BMC Microbiologyvolume 13, Article number: 72 (2013)
GeoChip 3.0, a microbial functional gene array, containing ~28,000 oligonucleotide probes and targeting ~57,000 sequences from 292 functional gene families, provided a powerful tool for researching microbial community structure in natural environments. The alpine meadow is a dominant plant community in the Qinghai-Tibetan plateau, hence it is important to profile the unique geographical flora and assess the response of the microbial communities to environmental variables. In this study, Geochip 3.0 was employed to understand the microbial functional gene diversity and structure, and metabolic potential and the major environmental factors in shaping microbial communities structure of alpine meadow soil in Qinghai-Tibetan Plateau.
A total of 6143 microbial functional genes involved in carbon degradation, carbon fixation, methane oxidation and production, nitrogen cycling, phosphorus utilization, sulphur cycling, organic remediation, metal resistance, energy process and other category were detected in six soil samples and high diversity was observed. Interestingly, most of the detected genes associated with carbon degradation were derived from cultivated organisms. To identify major environmental factors in shaping microbial communities, Mantel test and CCA Statistical analyses were performed. The results indicated that altitude, C/N, pH and soil organic carbon were significantly (P < 0.05) correlated with the microbial functional structure and a total of 80.97% of the variation was significantly explained by altitude, C/N and pH. The C/N contributed 38.2% to microbial functional gene variation, which is in accordance with the hierarchical clustering of overall microbial functional genes.
High overall functional genes and phylogenetic diversity of the alpine meadow soil microbial communities existed in the Qinghai-Tibetan Plateau. Most of the genes involved in carbon degradation were derived from characterized microbial groups. Microbial composition and structures variation were significantly impacted by local environmental conditions, and soil C/N is the most important factor to impact the microbial structure in alpine meadow in Qinghai-Tibetan plateau.
Microorganisms are the most abundant and diverse groups of organisms known on our planet, which play key roles in ecosystems and biogeochemical cycling of carbon, nitrogen, sulfur, phosphorus, and metals and biodegradation or stabilization of environmental contaminants [1–3]. Therefore, understanding microbial community structure, diversity, function and their relationships with environmental factors and ecosystem functioning is essential for the research of community formation and sustainability of life on our planet, which facilitates the management and protection of our natural environments [3, 4]. Numerous studies have been conducted to investigate the microbial community structure, diversity and their relationships with environments. Some studies showed that the microbial community is very sensitive to environmental changes, compared to plants and animals [5–8]. However, understanding is still limited on soil microbial communities in terms of structure, composition, and functional activity and their impact and response on environmental variations, especial for some special environments.
A large number of molecular approaches were developed and applied to analyze microbial diversity in the last two decades. Among them, high-throughput genomics technologies have shown great potential to study microbial diversity and the driving forces of different ecosystem processes as well as their response to different geological locations and environment changes [8–10]. GeoChip contains probes corresponding to genes encoding key enzymes involved in various biogeochemical cycling, thus it provided rapid, specific, sensitive and potentially quantitative analysis for microbial communities and was useful for studying the functional diversity and dynamics of microbial communities in different natural environments [8, 11–14]. Geochip 3.0, containing ~28,000 oligonucleotide probes and targeting ~57,000 sequences from 292 functional gene families, has been used to analyze microbial communities from different habitats of soils [14, 15], water , oil fields , marine sediments  and contaminated sites [10–19]. These studies showed that this Geochip served as a powerful tool for researching microbial community structure in natural environments .
The Qinghai-Tibetan Plateau, which extends over 2.5 million km2, is the youngest, highest and largest geo-morphological unit on the Eurasian continent , and was considered “The third pole of Earth”. However, this area also is a key region very sensitive to the impact of global warming. Therefore, the Qinghai-Tibetan Plateau has important significant values in scientific researches . The alpine meadow ecosystem, covering about 35% of the plateau area, is the dominant plant community type of the Qinghai-Tibetan Plateau . Kobresia, as one of the dominant genera of alpine meadows, is a typical vegetation on the Qinghai-Tibetan Plateau . At present, some studies found that the majority and diversity functional genes involved in nitrogen fixation and denitrifying existed in the alpine meadow in Qinghai-Tibetan plateau, and altitude and C/N ratio are the important environmental parameters affecting the activity of soil bacteria [20, 21]. However, little is known about the functional diversity and metabolic potential at the community level in the alpine meadow, especially for the Kobreasia, and the relationship between the functional gene structure of microbial communities and the surrounding environmental factors remains unclear .
In this study, Geochip 3.0 was employed to address two key questions. (i) what are the microbial functional gene diversity and structure, and metabolic potential of alpine meadow soil in Qinghai-Tibetan Plateau? (ii) what are the major environmental factors in shaping microbial communities structure in alpine meadow? To answer these questions, six soil samples were obtained and analyzed from the alpine meadow in the center part of the Qinghai-Tibetan Plateau, China.
Site description, sample collection, and geochemical analysis
The study sites were located in Sanjiangyuan Nature Reserve (89°24′-102°23′E, 31°39′-36°16′N), in the center of the Qinghai-Tibetan Plateau, China. Kobresia, as one of the dominant genera of alpine meadows, is a typical vegetation on the Qinghai-Tibetan Plateau. Six sites of typical Kobresia vegetation were selected in this study (Table 1). At each site, three 2 m × 2 m plots comprising typical vegetation were set up and the distance between nearly plots was about 20 m. Five to eight soil cores from the upper layer (0-15 cm) at a diameter of 1.5 cm were collected and mixed equally at each plot, and three plots were mixed and formed a soil sample at each site. Soil samples were stored at -20oC.
Soil samples were air-dried, sieved < 2 mm and analysed for pH (1:2 soil to H2O ratio), total organic carbon, total nitrogen, total phosphorus, total potassium, available potassium, available phosphorus as previously described .
Soil DNA extraction, purification and labeling
Microbial community genomic DNA was extracted directly from a 5 g soil sample by using a protocol that included liquid nitrogen grinding, freezing and thawing, and treatment with sodium dodecyl sulfate for cell lysis, which has been previously described . Then DNA was purified twice using 0.5% low melting point agarose gel followed by phenol-chloroform-butanol extraction. Purified DNA was quantified with an ND-1000 spectrophotometer (Nanodrop Inc.) and Quant-It PicoGreen (invitrogen, Carlsbd, CA). 3 μg of amplified DNA was labeled with a Cy5 fluorescent dye (GE Healthcare) by a random priming method .
DNA microarray hybridization, scanning and data processing
GeoChip 3.0 was used for DNA hybridization and this Geochip contains DNA probes targeting a total of 57,000 genes involved in key microbial processes . All hybridizations were carried out at 45°C for 10 h with 50% formamide using a TECAN HS4800. Arrays were scanned by using the ScanArray 5000 analysis system (Perkin-Elmer, Wellesley, MA). Signal intensities of each spot were measured with ImaGene 6.0 (Biodiscovery Inc., EI Segundo, CA, USA) and only the spots automatically scored as positive in the output of raw data were used for further data analysis . Spots with a signal-to-noise ratio [SNR = (signal intensity-background intensity)/standard deviation of the background] greater than 2.0 were used for further analysis.
Functional gene diversity was calculated by using Simpson’s reciprocal index (1/D) and Shannon-Weaver index (H’) using online software (http://www2.biology. ualberta.ca/jbrzusto/krebswin/html). Hierarchical clustering analysis of whole functional genes was performed using by the unweighted pairwise average-linkage clustering algorithm with CLUSTER (http://rana.lbl.gov/EisenSoftware.htm) and visualized by TREEVIEW software . The mantel tests were performed using R 2.9.1 (http://www.r-project.org/). Canonical correspondence analysis (CCA) is a multivariate ordination method and was performed to analysis the relationship between microbial community and environment variables by using the program package Canoco for Windows 4.5 (Biometris, The Nerherlands).
Geochemical properties in sampling sites
Soil characteristics of these six sampling sites are summarized in Table 1. pH in all those sites was neutral or close to alkali, and they were rich in organic carbon (C) and nitrogen (N), ranging from 91.99 g/kg to 209.19 g/kg and 1.50 g/kg to 15.50 g/kg, respectively. It was noted that C/N ratios displayed a decreasing trend as the elevation increased. For example, sample SJY-GH with the lowest elevation (3400 m) had the highest value of C/N ratio, whereas sample SJY-YS with the highest elevation (4813 m) had the lowest C/N ratio. In addition, sample SJY-GH had higher total C, N, P and K contents than the other samples.
Overview of functional gene diversity and structure of soil microbial communities
The examined microbial communities showed high diversity, as judged by the number of detected genes, overlapping genes between samples, unique genes and diversity indices (Table 2). The total number of detected genes ranged from 1,732 to 3,746 among the six study sites (Table 2). For instance, twice as many genes were detected in sample SJY-GH as in sample SJY-CD, SJY-ZD or SJY-YS. These samples had different community compositions, as shown by the unique and overlapped genes (Table 2). Sample SJY-GH and sample SJY-DR had the most overlapped genes (2029, 42.94%), while sample SJY-GH and sample SJY-YS had the fewest overlapped genes (1178, 27.22%). Simpson’s reciprocal diversity index (1/D) was the highest in sample SJY-GH and the lowest in sample SJY-CD (3716 and 1723, respectively). Similar results were obtained with Shannon-Weaver index (Table 2).
According to the phylogenetic analysis, the Proteobacteria group is the most dominant bacteria in all six samples, which account for over 56% (over 23% belong to α-proteobacteria, 13% belong to β-proteobacteria, 14% belong to γ-portecobacteria) among all the detected genes (Additional file 1: Table S1). The Actinobacteria (over 9.30%) and firmicutes (3.73%) are the second and third dominant bacteria groups. The genes derived from the Archaea and Eukaryota also were detected and accounted for 1.64% to 2.04% and 4.35% to 5.33% among all the detected genes in all samples, respectively. Although gene numbers belonging to different phylogenetic structure varied considerably in different samples, the proportions of genes number of different phylogenetic structure in all detected genes is similar. For example, the ratio of α-Proteobacteria ranged from 23.18% to 24.99% and the ratio of Actinobacteria ranged from 9.30% to 10.97% (Additional file 1: Table S1). Therefore, these results indicated the overall functional genes as well as the phylogenetic diversity of these alpine meadow soil microbial communities appeared to be quite high.
Analysis of detected functional genes
Among the 6143 genes detected in at least one sample, 567 were involved in carbon degradation, 202 in carbon fixation, 36 in methane oxidation, 18 in methane production, 754 in nitrogen cycling, 153 in phosphorus utilization, 279 in sulphur cycling, 2540 in organic remediation, 1275 in metal resistance, 126 in energy process, 193 in other category. Detected functional genes among these six alpine meadow soil samples were analyzed by hierarchical clustering (Additional file 1: Figure S1). A total of 39 different clusters of genes were observed. Genes in group 5, group 32 and group 35 are presented in all of the samples. The most obvious patterns were group 11 (1054 [17.16%]) and group 33 (373 [6.07%]); instead of, the genes in group 11 is only present at sample SJY-GH which is the lowest altitude sample and group 33 is only present at sample SJY-YS which is the highest altitude sample. The genes in group 11 were from functional categories involved in carbon degradation, carbon fixation, denitrification, nitrification, nitrogen fixation, phosphorus utilization, sulfite reductase, etc. Most of the genes in group 33 are involved in the carbon degradation, denitrification, nitrogen fixation, organic remediation, etc. These results showed that different microbial community structures existed in these samples and environment factors may influence them. To better understand microbial diversity involved in soil carbon cycling and nitrogen cycling, selected gene groups were further analyzed.
Functional genes involved in the carbon cycling
Microbe-mediated carbon cycling is one of the most important and complex process in the biogeochemical cycling. A total of 5196 gene probes belonging to carbon cycling were detected in the Geochip 3.0 . Among them, 823 gene probes were detected in all six soil samples (Table 3). Sample SJY-GH and SJY-CD have the most and least detected gene numbers, respectively. Carbon fixation and carbon degradation are the two most important gene categories in the carbon cycling in all samples (Table 3). The key enzymes (CODH, FTHFS, pcc and rubisco) involved in carbon fixation were detected in all samples and most of the detected genes (198/202) are derived from cultured bacteria. For example, fourteen genes were derived from Rhodopseudomonas palustris, four genes were derived from Xanthobacter autotrophicus, four genes were derived from Verminephrobacter eiseniae, three genes were derived from Roseiflexus Sp. and two genes were derived from Burkholderia xenovorans. However, only a few number of genes (10/202) involved in carbon fixation were shared by all six samples and Roseiflexus Sp. and Burkholderia xenovorans have high signal intensity in all of these soil samples.
Genes involved in the degradation of starch, cellulose, hemicellulose, chitin, lignin and pectin also were detected in Geochip and 161, 63, 105, 76, 90 and 12 gene probes were detected in all six samples (Table 3). All of the detected genes involved in the degradation of starch, cellulose and hemicellulose were derived from the cultured bacteria, and over 80% detected genes involved chitin, lignin and pectin (72/76, 85/90 and 10/12, respectively) were derived from cultured bacteria. However, only a few genes involved in the degradation of starch, cellulose, hemicellulose, chitin, lignin and pectin (14/161, 5/63, 6/105, 8/76, 8/90 and 0/12, respectively) were shared by all six samples. For methane cycle, a higher gene number and signal intensity of methane oxidation genes (mmoX and pmoA) were detected than that of methane production genes (mcrA) in all six samples. Most of the genes involved in methane oxidation and production (32/36 and 16/18) are derived from the uncultured microorganisms.
Most of shared genes involved in carbon cycling have high signal intensity in all the samples. For example, cellobiase gene involved in cellulose degradation derived from Roseiflexus castenholzii DSM 13941 was abundant in and shared by all six samples (Additional file 1: Figure S2), and gene derived from Rhodococcus sp. RHA1, Trichoderma harzianum and Arthrobacter sp. FB24 were also abundant. These results indicated that all of the processes involved in carbon cycling existed in the alpine meadow, and there were abundant genes diversity and most of the genes derived from the cultured bacteria in the Qinghai-Tibetan Plateau.
Functional genes involved in the nitrogen cycling
A total of 3763 gene probes belonging to different key gene categories involved in nitrogen fixation, denitrification, nitrification, dissimilatory N reduction, assimilatory N reduction and anaerobic ammonium oxidation are present in Geochip 3.0 . Among them, 754 gene probes were detected in all six soil samples (Table 3). 224, 372, 17, 51, 27 and 63 genes involved in nitrogen fixation, denitrification, nitrification, dissimilatory N reduction, assimilatory N reduction and anaerobic ammonium oxidation were detected in all samples, respectively (Table 3). Sample SJY-GH and SJY-CD have the most and least detected gene number, respectively.
Microbe-mediated nitrogen fixation and denitrification are the most important processes in nitrogen cycling. Microbe-mediated nitrogen fixation is the most important source of nitrogen in natural ecosystems, and occurs across a wide range of bacteria phyla, from Archaebacteria to Eubacteria . The majority of nifH genes (155/224) were derived from unidentified or uncultured organisms retrieved from different environments. Among nifH genes, 19 were shared by all samples. The shared gene 44829093 derived from an uncultured bacterium was dominant in samples SJY-GH and SJY-YS, and 780709 from an unidentified marine eubacterium was the most dominant gene in sample SJY-CD. These samples had a relatively high abundance of genes involved in nitrogen fixation.
Denitrification is a dissimilatory process of denitrifying bacteria where oxidized nitrogen compounds are used as alternative electron acceptors and nitrogen is transferred into the atmosphere in form of N2. Most of the detected genes involved in denitrification (320/372) were derived from the unidentified or uncultured organisms retrieved from different environments. These samples had a relatively high abundance of genes involved in denitrification (Table 3). 67 nosZ genes which encoding nitrous oxide reductase and it is considered a key enzyme in the denitrification process were detected. Few genes (13/67) were derived from the isolated bacteria. Four genes were shared and derived from the uncultured bacteria by all six soil samples (Additional file 1: Figure S3).
Together, these results indicated that all the processes involved in nitrogen cycling existed, and there were high gene diversity as well as high potential metabolic ability in nitrogen fixation and denitrification in all these samples.
Relationships between microbial community structure and environmental variables
To assess the relationships between microbial community structure and soil environmental variables, Mantel test and canonical correspondence analysis (CCA) were used. Mantel tests of all six soil samples were performed with 12 individual environmental variables. The environmental variables of altitude, C/N, pH and soil organic carbon were used to analyze the correlation with the microbial functional genes involved in carbon and nitrogen cycling (Additional file 1: Table S2), suggesting that these environmental variables play important roles in shaping the microbial community structure in these soil samples. Consistently, CCA results showed that the C/N and altitude were the most important factors when only significant environmental variables (altitude, C/N, pH and organic carbon) were included in the CCA biplot (Figure 1). Samples of SJY-DR, SJY-CD, SJY-ZD and SJY-QML clustered together which were separated from in SJY-GH and SJY-YS (Figure 1). On the basis of the relationship between environmental variables and microbial functional structure, altitude seemed to be the most important variable affecting the microbial functional structure. Notably, sample SJY-GH was collected at a low altitude (3400 m), while sample SJY-YS was collected at a high altitude (4813 m), while the altitude of Sample SJY-DR, SJY-CD, SJY-ZD and SJY-QML was 4000-4500 m.
Variance partitioning analysis was used to quantify the contributions of altitude (A), soil chemistry (S) and pH (p) to the microbial community variation. The total variation was partitioned into the independent effects of A, S and pH (when the effects of all other factors were removed), interactions between only two factors, common interactions of all three factors and the unexplained portion (Figure 2a). On the basis of Geochip data, a total of 80.97% of the variation was significantly explained by these three environmental variables (Figure 2b). Altitude, C/N and pH were able to independently explain 18.11%, 38.23% and 19.47% of the total variations observed, respectively. Interactions between any two factors or among the three factors seemed to have less effect than the individual factors. Only about 20% of the community variation could not be explained by these three environmental variables.
Analyzing microbial functional genes involved in major biogeochemical processes is important to link microbial community structure to their potential ecological functions . In this study, we used GeoChip 3.0 to analyze microbial functional gene diversity in alpine meadow soil samples from the Qinghai-Tibetan plateau. This report was one of the first ecological applications of an expanded functional gene microarray [13, 30], and it is the first application of this kind for studies in Qinghai-Tibetan plateau, China. These results indicated the overall functional genes as well as the phylogenetic diversity of these alpine meadow soil microbial communities is higher than in the Antarctic latitudinal transect or alpine soil in the Colorado Rocky Mountains [30, 31]. All the detected genes involved in the carbon degradation, carbon fixation, methane oxidation and production, nitrogen cycling, phosphorus utilization, sulphur cycling, organic remediation, metal resistance, energy process, and other category. According to the phylogenetic analysis, the proteobacteria group is the most dominant bacteria in all six samples, which account for over 56% among all the detected genes. Therefore, Proteobacteria maybe the most prevalent bacteria in Qinghai-Tibetan plateau.
Soil is the major reservoir of terrestrial organic carbon, and soil carbon degradation is largely controlled by the metabolic activities of the microorganisms present in the soil [32, 33]. The majority of microbial studies have monitored the relationship between organic carbon in soil, CO2 release, and microbial biomass in different soil types [34, 35]. In this study, metabolic genes involved in the degradation of starch, cellulose, hemicellulose, chitin, lignin and pectin were detected and the individual gene orthologs were abundant and diverse. For example, 76 genes related to lignin degradation were detected and the number of genes detected was 53, 37, 31, 23, 22 and 23 in SJY-GH, SJY-DR, SJY-QML, SJY-CD, SJY-ZD and SJY-YS, respectively. These detected genes related to lignin degradation belonged to 4 different gene families, including laccase, glyoxal oxidase, lignin peroxidase and manganese peroxidase, and most of the detected genes (94.59%) were derived from the isolated organisms (e.g., 17.57% from Phanerochaete sp.). Most of the shared genes were abundant in all the samples. For example, the cellobiase gene involved in cellulose degradation derived from Roseiflexus castenholzii DSM 13941 was shared by all of the six samples and had the highest signal intensity in all samples.
Understanding the environmental variables that affect microbial community structure is a key goal in microbial ecology . Different environmental variables affect the microbial structure and potential activity on ecosystem functions . He et al  found that the abundance of all detected genes was significantly (P < 0.05) and positively correlated with soil moisture and pH. Yergeau et al.  used the Geochip to examine soil microbial communities across an Antarctic latitudinal transect and revealed that cellulose degradation and denitrification genes were correlated with soil temperature. Our results showed that altitude, C/N, pH and available phosphorus had a significant impact on the microbial functional communities in alpine meadow soil, suggesting that these environmental variables play an important role in shaping microbial community structure. However, we know very little about how microbial distribution pattern varies along altitude gradients . This is a considerable gap in understanding microbial biodiversity and will likely be an important component of ecosystem response to global warming [37, 38].
Variation partitioning analysis in this study showed that a total of 80.97% of the variation was significantly explained by altitude, C/N and pH. The C/N contributed the most (38.2%) to microbial functional gene variation, which is in accordance with the hierarchical clustering of overall microbial functional genes, indicating a significant impact of local environmental conditions on the composition and structures of microbial communities. In this study, only 19.03% of the variation of microbial community structure could not be explained by of these three factors, which showed that considerable amounts of variations could be explained by environmental variables measured. However, some previous studies thought that most of the variation could be explained by environmental variables. For example, Zhou et al.  showed that more than 50% of variations in a forest soil community could not be explained by both environmental factors and geographic distance. Ramette and Tiedje  showed that 34-80% of microbial variations could not be explained by measured environmental variables in agricultural soils. Liang et al  indicated over 40% of the variations of microbial community could not be explained by geographic location, soil geochemical variables and oil contamination.
In summary, soil microbial functional gene diversity in alpine meadow in Qinghai-Tibetan plateau was examined by Geochip 3.0 and almost all genes involved in carbon, nitrogen and other element cycling were found, which showed that the microbial functional diversity in alpine meadow ecosystem was quietly high. Statistical analyses showed that the microbial communities may be shaped largely by the altitude, C/N, and pH. However, Geochip analyzed the distribution of metabolic genes may reflect the metabolic potential of the microbial community , but not necessarily the actual populations. For example, we detected many key enzyme genes involved in carbon degradation, which implied that the populations carrying those genes could exist in the alpine meadow ecosystem, but it does not mean that they express the enzymes of degradation organic carbon. Therefore, further analysis of the functional activity with different approaches such as mRNA-based microarray hybridization is needed to address it .
A highly overall functional genes and phylogenetic diversity of the alpine meadow soil microbial communities existed in the Qinghai-Tibetan Plateau. Most of the genes involved in carbon degradation were derived from characterized microbial groups. The considerable amounts of microbial the composition and structures variation was significant impacted by local environmental conditions, and the C/N is the most important factors to impact the microbial structure in alpine meadow in Qinghai-Tibetan plateau.
Availability of supporting data
The data set supporting the results of this article is available in the microarray data repository, unique persistent identifier and hyperlink to dataset(s) in http://ieg2.ou.edu/NimbleGen/analysis.cgi
Fierer N, Jackson RB: The diversity and biogeography of soil bacteria communities. PNAS. 2006, 103 (3): 626-631. 10.1073/pnas.0507535103.
Green JL, Bohannan BJM, Whitaker RJ: Microbial biogeography: from taxonomy to traits. Science. 2008, 320: 1039-1043. 10.1126/science.1153475.
He Z, Nostrnd JDV, Deng Y, Zhou J: Development and applications of functional gene microarrays in the analysis of the functional diversity, composition, and structure of microbial communities. Front Environ Sci Engin China. 2011, 5 (1): 1-20.
Meeteren MJ, Tietema A, Loon E, Verstraten J: Microbial dynamics and litter decomposition under a changed climate in a Dutch heathland. Appl Soil Ecol. 2008, 38: 119-127. 10.1016/j.apsoil.2007.09.006.
Willig MR, Kaufman DM, Stevens RD: Latitudinal gradients of biodiversity: pattern process, scale and synthesis. Annu Rev Eco Evol Sys. 2003, 34: 273-309. 10.1146/annurev.ecolsys.34.012103.144032.
Currie DJ, Francis AP: Regional versus climate effect on taxon richness in angiosperms: reply to Qian and Ricklefs. Am Nat. 2004, 163: 780-785. 10.1086/383596.
Qian H, Ricklefs RE: Taxon richness and climate in angiosperms: is there a globally consistent relationship that precludes region effect?. Am Nat. 2004, 163: 773-779. 10.1086/383097.
Zhou J, Kang S, Schadt CW, Garten CT: Spatial scaling of functional gene diversity across various microbial taxa. PNAS. 2008, 105 (22): 7768-7773. 10.1073/pnas.0709016105.
Waldron PJ, Wu L, Joy D, Schadt CW, He Z, Watson DB, Jardine PM, Palumbo AV, Hazen TC, Zhou J: Functional gene array-based analysis of microbial community structure in ground-waters with a gradient of contaminant levels. Environ Sci Technol. 2009, 43: 3529-3534. 10.1021/es803423p.
Wang F, Zhou H, Meng J, Peng X, Jiang L, Sun P, Zhang C, Joy D, Deng Y, He Z, Wu L, Zhou J, Xiao X: Geochip-based analysis of metabolic diversity of microbial communities at the Juan de Fuca ridge hybrothermal vent. PNAS. 2009, 106 (12): 4840-4845. 10.1073/pnas.0810418106.
Zhou J, Thompso DK: Challenges in applying microarrays to environmental studies. Curr Opin Biotechnol. 2002, 13: 204-207. 10.1016/S0958-1669(02)00319-1.
Rhee SK, Liu X, Wu L, Chong SC, Wan X, Zhou J: Detection of genes involved in biodegradation and biotransformation in microbial communities by using 50-mer oligonucleotide microarrays. Appl Environ Microbiol. 2004, 70: 4303-4317. 10.1128/AEM.70.7.4303-4317.2004.
He ZL, Gentry TJ, Schadt CW, Wu L, Liebich J, Chong SCZ: Geochip: a comprehensive microarray for investigating biogeochemical, ecological and environmental processes. ISME J. 2007, 1: 67-77. 10.1038/ismej.2007.2.
He ZL, Deng Y, Nostrand JD, Tu Q, Xu M, Hemme CL, Li X, Wu L, Gentry TJ, Yin Y, Liebich J, Hazen TC, Zhou J: GeoChip 3.0 as a high-throughput tool for analyzing microbial community composition, structure and functional activity. ISME J. 2010a, 4: 1-13. 10.1038/ismej.2009.83.
He Z, Xu M, Deng Y, Kang S, Kellogg L, Wu L, Nostrand JD, Hobbie SE, Reich PB, Zhou J: Metagenomic analysis reveals a marked divergence in the structure of belowground microbial communities at elevated CO2. Ecol Lett. 2010, 13: 564-575. 10.1111/j.1461-0248.2010.01453.x.
Lu Z, He Z, Parisi VA, Kang S, Deng Y, Van Nostrand JD, Masoner JR, Cozzarelli IM, Suflita JM, Zhou J: Geochip-based analysis ofmicrobial functional gene diversity in a landfill leachate-contaminated aquifer. Environ Sci Technol. 2012, 46: 5824-5833. 10.1021/es300478j.
Liang Y, Nostrand JD, Deng Y, He Z, Wu L, Zhang X, Li G, Zhou J: Functional gene diversity of soil microbial communities from five oil-contaminated fields in China. ISME J. 2011, 5: 403-413. 10.1038/ismej.2010.142.
Xiong JB, Wu LY, Tu SX, van Nostrand JD, He ZH, Zhou JZ, Wang GJ: Microbial communities and functional genes associated with soil arsenic contamination and the rhizosphere of the arsenic-hyperaccumulating plant Pteris vittata L. Appl Environ Microbiol. 2010, 76: 7277-7284. 10.1128/AEM.00500-10.
Xu M, Wu WM, Wu L, He Z, Van Nostrand JD, Deng Y, Luo J, Carley J, Ginder-Vogel M, Terry JG, Baouhua G, David W, Philip MJ, Terence LM, James MT, Terry H, Craig SC, Zhou J: Responses of microbial community functional structures to pilot-scale uranium in situ bioremediation. ISME J. 2010, 4: 1060-1070. 10.1038/ismej.2010.31.
Zhang Y, Li D, Wang H, Xiao Q, Liu X: Molecular diversity of nitrogen-fixing bacteria in Tibet plateau, China. FEMS Microbiol Lett. 2006, 260: 134-142. 10.1111/j.1574-6968.2006.00317.x.
Zhang Y, Li D, Wang H, Xiao Q, Liu X: The diversity of denitrifying bacteria in the alpine meadow soil of Sanjiangyuan natural reserve in Tibet Plateau. Chin Sci Bull. 2006, 51 (8): 1-8.
Cao G, Tang Y, Mo W, Wang Y, Li Y, Zhao X: Grazing intensity alters soil respiration in an alpine meadow on the Tibetan plateau. Soil Biol Biochem. 2004, 36 (2): 237-243. 10.1016/j.soilbio.2003.09.010.
China Vegetation Edits Commission: China vegetations. 1980, Beijing: Science Press
Nemergut DR, Costello EK, Meyer AF, Pescador MY, Weintraub MN, Schmidt SK: Structure and function of alpine and arctic soil microbial communities. Res Microbiol. 2005, 156: 775-784. 10.1016/j.resmic.2005.03.004.
Bao SD: Soil and agricultural chemistry analysis. 1999, Beijing: China Agriculture Press, 25-150.
Richard AH, Qiu XY, Wu LY, Roh Y, Palumbo AV, Tiedje JM, Zhou JZ: Simultaneous recovery of RNA and DNA from soils and sediments. Appl Environ Microbiol. 2001, 67: 4495-4503. 10.1128/AEM.67.10.4495-4503.2001.
Wu L, Kellogg L, Devol AH, Tiedje JM, Zhou J: Microarray-based characterization of microbial community functional structure and heterogeneity in marine sediments from the gulf of Mexico. Appl Environ Microbiol. 2008, 74 (14): 4516-4529. 10.1128/AEM.02751-07.
Young JPW: Phylogenetic classification of nitrogen fixing organisms. Biological nitrogen fixation. Edited by: Stacey G, Burries RH, Evans HJ. 1992, New York: Chapman and Hall, 43-86.
Torsvik V, Ovreas L: Microbial diversity and function in soil: from genes to ecosystems. Curr Opin Microbiol. 2002, 5: 240-245. 10.1016/S1369-5274(02)00324-7.
Yergeau E, Kang S, He Z, Zhou J, Kowalchuk GA: Functional microarray analysis of nitrogen and carbon cycling genes across an Antarctic latitudinal transect. ISME J. 2007, 1: 163-179. 10.1038/ismej.2007.24.
David AL, Steven KS: Seasonal changes in an alpine soil bacteria community in the Colorado rocky mountains. Appl Environ Microbiol. 2004, 70 (5): 2867-2879. 10.1128/AEM.70.5.2867-2879.2004.
Ross DJ, Tate KR, Scott NA, Feltham CW: Land-use change: effects on soil carbon, nitrogen and phosphorus pools and fluxes in three adjacent ecosystems. Soil Biol Biochem. 1999, 31: 803-813. 10.1016/S0038-0717(98)00180-1.
Zhang Y, Zhang X, Liu X, Xiao Y, Qu L, Wu L, Zhou J: Microarray-based analysis of changes in diversity of microbial genes involved in organic carbon decomposition following land use/cover changes. FEMS Microbiol Lett. 2007, 266: 144-151. 10.1111/j.1574-6968.2006.00511.x.
Wang WJ, Dalal RC, Moody PW, Smith CJ: Relationships of soil respiration to microbial biomass, substrate availability and clay content. Soil Biol Biochem. 2003, 35: 273-284. 10.1016/S0038-0717(02)00274-2.
Michelsen A, Andersson M, Jensen M, Kjoller A, Gashew M: Carbon stocks, soil respiration and microbial biomass if fire-phone tropical grassland, woodland and forest ecosystems. Soil Biol Biochem. 2004, 36: 1707-1717. 10.1016/j.soilbio.2004.04.028.
Bryant JA, Lamanna C, Morlon H, Kerkhoff AJ, Enquist BJ, Green JL: Microbes on mountainsides: Contrating elevational patterns of bacterial and plant diversity. PNAS. 2008, 105 (suppl.1): 11505-11511.
Carney KM, Hungate BA, Drake BG, Megonigal JP: Altered soil microbial community at elevated CO2 leads to loss of soil carbon. PNAS. 2007, 104 (12): 4990-4995. 10.1073/pnas.0610045104.
Monson RK: Winter forest soil respiration controlled by climate and microbial community composition. Nature. 2006, 439: 711-714. 10.1038/nature04555.
Ramette A, Tiedje J: Multiscale responses of microbial life to spatial distance and environmental heterogeneity in a patchy ecosystem. Proc Natl Acad Sci USA. 2007, 104: 2761-2766. 10.1073/pnas.0610671104.
This research was supported by the Public Welfare Project of the National Scientific Research Institution (CAFRIFEEP201101, CAFRIF200713) and National Natural Science Foundation of China (No. 30700018).
We declared that this manuscript have not any finical competing interests. We have not received reimbursements, fees, funding, or salary, or hold any stocks or shares from any organizations that may in any way gain or lose financially from the publication of this manuscript, either now or in the future. We also have not hold or apply any patents relating to content of the manuscript. No other financial competing interests are related to this manuscript. We declared that this manuscript have not any non-financial competing interests (political, personal, religious, ideological, academic, intellectual, commercial or any other).
Y Z carried out the lab design, sampling collecting, data analysis and the manuscript preparation. Z L carried out the soil microbial DNA extraction, microarray hybridization, scanning and data processing. S L participated the microarray data analysis. Y Y participated the microarray data analysis and manuscript preparation. Z R participated the sampling collecting and biogeochemical data analysis. J Z participated the lab design and data analysis. D L participated the lab design, data analysis and manuscript preparation. All authors read and approved the final manuscript.