Soil bacterial communities shaped by geochemical factors and land use in a less-explored area, Tibetan Plateau
- Xiangyu Guan†1, 3,
- Jinfeng Wang†2,
- Hui Zhao2,
- Jianjun Wang4,
- Ximing Luo3,
- Fei Liu1Email author and
- Fangqing Zhao2Email author
© Guan et al.; licensee BioMed Central Ltd. 2013
Received: 9 June 2013
Accepted: 20 November 2013
Published: 22 November 2013
As the largest low-latitude permafrost region, the Tibetan Plateau (TP) is an important part of the earth’s terrestrial ecosystem and one of the most vulnerable areas to climate change and human activities. However, to the best of our knowledge, the bacterial communities in TP soils and their roles in biogeochemical cycles remain limited.
In this study, we report the bacterial community structure and function as well as their correlation with environmental factors in TP major ecosystems (farmland, alpine meadow and oligosaline lake) by using metagenomic approaches. Compared with other soil samples in various environments, TP soils share a core set of microorganisms with a distinct abundance and composition. Among TP soil samples, the taxonomic and functional composition of bacterial communities among the upper (3-5 cm) and lower (18-20 cm) soils of farmland sites were highly similar, whereas the dissimilarities within alpine meadow samples were significantly greater than among farmland samples. A similar pattern was observed in elements cycles and pathways associated with adaption to environment and land use types. Canonical correlation analysis revealed that the bacterial communities in most of farmland and alpine meadow soil samples were also significantly correlated with geogenic variables. Specifically, the root-nodule bacteria are negatively correlated with the soil moisture and pH, while Thiobacillus associated with sulfur cycles show potential responses to low temperature and intense UV radiation.
These findings indicate that the bacterial community structure and functions in TP soils were influenced by both human activities and soil environmental properties, and that the bacterial communities appeared to be more homogenized in the farmland soils compared with pristine alpine meadows.
Soil is an essential ecosystem for organisms living on land. Soil microorganisms contribute to the breakdown of organic residues, element cycling and soil mineralization. Recent advances in high-throughput sequencing technologies have enabled comprehensive analysis of the microbial composition and function in a variety of soil environments, including farmland, forest, grassland, tundra, permafrost and desert [1–8]. These studies have revealed the high diversity of microbes in different niches and improved our understanding of new taxa, biogeographic distribution, and the association of specific microbial groups with geochemical factors.
The distribution of bacterial communities in different soil environments is controlled by a number of variables, including physical and chemical characteristics (e.g. depth, geogenic factors, pH, salinity and temperature), microenvironment structure (e.g. pore spaces, water films and dead organic matter), and other living organisms (e.g. arthropods, fungi). Geogenic factors and landforms have long been recognized as the drivers of soil geochemistry, but geogenic factors are fundamental determinants of microbial communities in naturally metal-rich soils in Australia [9, 10]. Nitrogen fertilization was shown to significantly influence bacterial community composition in a grassland soil and in agricultural soil in the United States while stimulating the distribution of Proteobacteria and Bacteroidetes . The relative abundance of Acidobacteria, Actinobacteria, and Bacteroidetes changed with soil pH, which is believed to represent a prominent factor in bacterial community composition and proportion [6, 12]. The Tibetan Plateau (TP) has unique climatic and geochemical characteristics including high elevation (exceeding 4,500 meters on average), low temperature, limited nutrients and intense UV radiation [13–15]. Only a small number of studies have examined the microbial communities in glacier, snow, and lake sediments by sequencing of 16S rRNA in TP [16–19]. Consequently, the composition of microbial communities and the associated environmental factors in TP soils remain poorly explored.
Deep metagenomic analyses provide not only an understanding of microbial composition, but also comprehensive functional diversity of soil samples [1, 20]. Some integrated studies of microbial composition and function have provided insight into the relationship of taxonomic and functional structure and their different potential responses to environments across major global biomes such as permafrost thaw . Composition and function of soil microbial communities from cold deserts, hot deserts, forests, grasslands and tundra have been compared, and some functional genes associated with osmoregulation, dormancy and antibiotic resistance have been shown to be related to environments . Additionally, soil samples from different sampling times and depths in the same grassland and their bacterial structure and functional fluctuations have been sequenced and analyzed . Metagenomic approaches can be used to explain how particular pathways reflect the adaptation of microbial communities to environments. For example, the abundance of genes for amino acid synthesis, vitamin and cofactor metabolism, and lipid and glycan metabolism can indicate bacterial community adaptation to specific environmental challenges, nutrient-limited conditions and physicochemical conditions respectively .
As the third pole of the world, although TP is still in a relatively pristine state, pristine alpine meadows have been converted to farmlands in recent years [13–15]. We collected 20 soil samples from ten different locations to cover the major ecosystems of alpine meadow, farmland and lake in the TP. Among these ecosystems, alpine meadow soil was found to have significantly higher organic carbon content than soils from other types of ecosystems [14, 23, 24], which account for more than 60% of the TP [25–27]. Farmland is the major ecosystem affected by human activity in this region . Namco Lake (oligosaline water), which is one of the three largest lakes in TP, has shown increasing water levels owing to glacial melt [29, 30]. Previous studies have shown that bacterial communities and functions were not only affected by environmental heterogeneity, but can also be influenced by the strength of human activities and land use [11, 31, 32]. Moreover, the top 10 cm of soil layer was found to be more easily influenced and variable than the lower soil layer . Therefore, in this study, samples were collected from the upper (3-5 cm) and lower (18-20 cm) soil layers for assessment of bacterial community diversity and anthropogenic impact on them. Additionally, we compared whether the degree of similarity in bacterial communities or functions between surface and lower layers in farmlands differed from those in pristine alpine meadows. To accomplish this, we measured a number of physical and chemical properties of the TP soil samples and then performed metagenomic sequencing of the total DNA extracted from them. We then integrated this data to identify taxonomic and functional diversities of microbial communities in TP soils and to assess potential environmental factors and land use that may shape the bacterial communities.
Sample collection and soil analyses
Soil pH value was determined in a soil/water (1:5 w/v) suspension by using pH meter (Sartorius PB-10, Goettingen, Germany). Total organic carbon (TOC) was determined using the K2Cr2O7 oxidation method , and total nitrogen (TN) was analyzed by Eurovector elemental analyser (Isoprime-EuroEA 3000, Milan, Italy). Total sulfur (TS) was measured by an infrared absorption method after high frequency combustion (High-speed analyzer HWF-900A, Jinan, China). Other trace metals were analyzed by ICP-MS (Thermo Fisher X-series, Franklin, MA, USA) and ICP-AES (TJA IRIS-Advantage, Franklin, MA, USA). Analysis on volatile organic compounds (VOCs) of soil samples was performed with an automatic static headspace sampler. The samples were analyzed by a GC/MSD System (Agilent 6890/5973 N, Palo Alto, CA, USA). The reporting limits were 1.0 μg/L for all the analytes.
DNA extraction and metagenomic sequencing
DNA of homogenized samples (0.5 g) was extracted by the PowerSoil DNA extraction kit (MoBio Laboratories, Carlsbad, CA, USA), and stored at -20°C for further use and -80°C for permanent preservation. The quantity and quality of isolated DNA were evaluated using a Nano Drop spectrophotometer (Thermo Fisher Scientific ND-1000, Waltham, MA, USA) and agarose gel electrophoresis (Bio-Rad, Hercules, CA, USA), respectively. For each sample, 0.5 μg of purified metagenomic DNA was sheared into fragments of ~180 bp in length, and a library was constructed according to a standard protocol provided by Illumina, Inc. (San Diego, CA, USA). Quantification was performed using a Qubit Fluorometer (Invitrogen, Life Technologies, Grand Island, NY, USA) and a Stratagene Real-time PCR Cycler (Agilent Mx3000P, Santa Clara, CA, USA) prior to cluster generation in a c-Bot automated sequencing system (Illumina, Inc., San Diego, CA, USA). Eight libraries with different indexes were pooled together and sequenced in one lane using an Illumina HiSeq 2000 high-throughput sequencing instrument with 2 × 100 bp paired-end (PE) sequencing. A total of two and a half lanes were sequenced for the 20 libraries.
Quality filtering and PE reads merging
Firstly, low quality PE reads were removed before further analysis by custom scripts. Secondly, quality-filtered reads were merged based on the overlap of PE reads. Briefly, we iteratively aligned a pair of reads (read 1 and read 2) with an overlap length ranging from 6 to 40 bp. In each iteration, the overlap score was calculated as the number of mismatches divided by the overlap length. If the score of the best overlap was smaller than the mismatch threshold (0.15), read 1 and read 2 were merged into a long read. The WGS sequencing data were deposited in the Sequence Read Archive (SRA) database with an accession number as SRP032429.
Sequences were aligned and annotated using the Basic Local Alignment Search Tool (BLASTX) algorithm and the National Center for Biotechnology Information (NCBI) non-redundant (NR) sequence database with an E-score cutoff value of <10-4, and the BLASTX hits were further processed by the Metagenome Analyzer program (MEGAN, V4) to statistically analyze the abundance of taxa in each sample . After normalizing sequence counts of each taxonomic group by the total number of reads, statistical analysis was performed on the bacterial composition and abundance at the phylum and genus levels. To compare the taxonomic communities between TP and other soil environments, metagenomic datasets were obtained from MG-RAST , and similar annotation analyses were carried out as below.
Metagenomic reads were assembled using SOAPdenovo2 v2.04 , and certain assembled contigs were further scaffolded using PGA . The MetaGeneMark v2.8 gene prediction tool  was used to predict genes from assembled contigs. The predicted proteins were compared against the NR database using BLASTP with an E-score cutoff value of <10-2 and were taxonomically and functionally annotated using a combination of lowest common ancestor and consensus approaches . Metagenomic sequences were functionally annotated using the Kyoto Encyclopedia of Genes and Genomes (KEGG) bioinformatics database  and MetaCV . Functional categories and genes in the KEGG pathway were counted in each sample for further functional analysis. Certain reads were extracted from KEGG sulfur metabolism pathway and BLASTXed against the NCBI NR protein database with an E-score cutoff value of <10-2. The alignment results were parsed by MEGAN to determine the abundance of each bacterium in contributing to sulfur metabolism. Complete genome of Thiobacillus denitrificans ATCC 25259 (NC_007404.1, 2.91 Mb) was chosen as a reference to calculate the sequence frequency of each genotype of Thiobacillus reads from samples of G5U, and Burrows-Wheeler Aligner (BWA) , SAMtools  and the Integrated Next-gen Genome Analysis Platform (inGAP) [45, 46] were used as mapping, calculating and visualization tools, respectively.
Statistical analyses and visualization
Alpha diversity was calculated by using the taxonomic and functional metrics (Simpson Diversity Index, 1/D). The dissimilarity in taxonomic and functional composition between samples was measured with Bray-Curtis distance . The variation of dissimilarity was compared between farm and meadow ecosystems by focusing on different sample groups, that is the upper, the lower and the lower vs. upper samples. The significance between the sample groups was tested using t-test with permutations (999). All the analyses were carried out in R.
Principal component analysis (PCA) implemented in imDEV  was performed to evaluate the similarity among the twenty metagenomic communities from TP based on the metadata of taxonomic composition, functional annotation and environmental factors. We further use Canonical Correlation Analysis (CCA) to measure the relationship between environmental factors and sampling sites in CANOCO 4.5 (Biometrics Wageningen, The Netherlands). The number of genera and functional reads were counted using the annotation results from MEGAN.
Geochemical measurements of 20 soil samples
Geochemical and field data for the sampling sites are shown in Additional file 1: Table S1 and S2. TN was the highest in alpine meadow G2, TS was significantly higher in alpine meadow samples G4U/L and G5U/L than all other samples, and G2 and G/L4 had higher TOC than other studied samples. The C: N and C: S ratios are important measures of geochemical cycling, and the C: N ratio can strongly influence the bacterial community structure [4, 49]. The C: N ratio in samples L1U/L, G5U/L, and G3U was more than 30:1. In samples G3L, G6U, and G4U/L, the C: N ratio was between 20 and 30, and nitrogen was neither fixed nor mineralized. The C: N ratio in the remaining samples was lower than 20, indicating that nitrogen was likely mineralized. The C: S ratio was below 200 in all samples, indicating sulfur mineralization. The phosphorous content of the samples was ranked as follows: farmland (~0.1%) > alpine meadow (<0.1%) > saline lake (<0.001%). The K: N and K: P ratios varied among soil samples, but all were more than 14 and 134, respectively, indicating they were rich in potassium. The VOCs were below the limit of detection in all samples.
Bacterial community composition and diversity in TP major ecosystems
Sequencing of both the upper (3-5 cm) and lower (18-20 cm) soil samples from 10 sampling sites in TP resulted in a total of 336 million read pairs. The paired reads were merged into 61 million long sequences (160-197 bp each), and these PE-merged sequences were compared against the NCBI NR protein database by BLASTX. In different samples, 48.2-58.3% of these PE-merged sequences could be assigned to eubacteria, 0.8-3.3% to eukarya, ~0.1% to both archaea and viruses, while the remaining reads could not be assigned.
New bacterial resources and similar bacterial core
To investigate the taxonomic diversity among samples, the similarity between assigned metagenomic reads and their reference sequences were compared and analyzed (Additional file 2: Figure S1). We found that the vast majority of reads shared 65-85% sequence identity with known species. One genus, Thiobacillus, had higher sequence identity with soil metagenomic sequences. The majority of metagenomic sequences classified as Thiobacillus had ≥ 90% identity at the amino acid level with known sequences from T. denitrificans ATCC 25259, which is a facultatively anaerobic, chemolithotrophic, sulfur-oxidizing bacterium with an ability to couple the oxidation of inorganic sulfur to the reduction of oxidized nitrogen compounds (such as nitrate, nitrite) to dinitrogen. Interestingly, we found that a small fraction of reads in F1 could be assigned to an insect, Acyrthosiphon pisum, with 90-95% sequence identity.
The distribution of the top 50 most abundant genera revealed a core set of taxa shared by the 20 samples and was different from other environmental samples such as tundra, grasslands, farms, forests, lakes, creeks and marine systems (Additional file 3: Figure S2). The core microbiota of TP samples included Candidatus Solibacter, Gemmatimonas, Sorangium, Nitrospira, Chthoniobacter and others. Samples from farmlands and alpine meadows shared a more similar core of genera than other environmental samples. It should be noted that samples from alpine (swamp) meadow G4U/L and its adjacent site G5U/L had the largest fraction of Thiobacillus. The bacterial composition and abundance of saline lake samples L1U/L were obviously different from other TP samples. The proportion of core genera in L1U/L differed from other TP samples in that it contained Thiobacillus, Conexibacter, Nitrospira, Burkholderia, Chitinophaga, and a few others.
Comparison of the bacterial communities between upper and lower layers of TP soils is shown (Figure 3D). Most of genera showed an unobvious diversity (grey color), and the difference between upper and lower layers samples were mostly low abundance genera (less 0.1%). The bacterial communities in both layers of F1, F2, and F3 were very similar, and a slight difference was only three microorganisms including Chitinophaga, Geobacter and Acyrthosiphon. Proglacial alpine meadow G2 was lowest diversity among alpine meadow samples, in which Nitrospira as nitrifying bacteria was the common genus at more than twice the amount of bacteria in the lower layer than in the upper layer except G2. Upper and lower layer samples from G5 were the more divergent than other sites samples, and the differences were mainly from genera including Thiobacillus, Geobacter, Gemmatimonas, Opitutus and so on. These observations showed that the bacterial communities at depth of 3-5 and 18-20 cm in soil from various farmland and alpine meadow soil environments were really similar with minor differences in root-associated bacteria and nitrifying bacteria (more abundant in upper soil than lower soil) and denitrifying bacteria (less abundant in upper than lower soil).
PCA and CCA based on metagenomic data
In addition to bacterial community composition, we also performed PCA of both environmental factors and functional annotations. Functional annotations were assigned to metagenomic sequences by a TBLASTX comparison to the non-redundant KEGG database. In contrast to the PCA of bacterial community composition, samples from sites G1 and G6 were clustered within the farmlands, and G2 was separated from the cluster of farmland samples, mainly correlated with moisture, soil carbon, nitrogen, and sulfur concentrations (Figure 5B). The sites of bacterial community functional profiles showed the greatest dispersion among the three PCA (Figure 5D). Although functional profiles of samples from F1-3 showed a closer relationship than samples from farmlands and lake, they could not be distinguished completely along PC1 (84%) (Figure 5D).
CCA was performed to explore the relationship of environmental factors and bacterial community composition. The results revealed a strong association between moisture, magnesium, calcium and saline water-saturated samples (L1U/L) (Figure 5C). G4 and G5 were clearly controlled by the TS and TOC concentration, while the other samples from farmland and meadow sites were clearly correlated with all minor elements and lithophile major elements except magnesium and calcium.
Nitrogen and sulfur-associated bacteria
Reads assigned to Thiobacillus in G5U were aligned to the complete genome sequence of T. denitrificans ATCC 25259, which revealed an uneven distribution of the depth of coverage (Additional file 5: Figure S3). High frequency genes included sulfur metabolism-related genes, the NADH dehydrogenase subunit, a heavy metal-related gene and a cold-shock DNA-binding protein family, whereas nitric oxide reductase occurred at a low frequency (Additional file 4: Table S6).
TP is the highest and largest plateau on earth and is considered one of the most fragile ecosystems . However, only a very limited number of studies have explored the microbial communities in TP soils by using 16 s rRNA gene clone library approaches to date [16, 50]. Therefore, in this study, we generated and analyzed 67.2 Gb metagenomic sequences from 20 soil samples to assess and quantify the variation in bacterial communities between major ecosystems in TP soils.
Comparison of the TP bacterial communities with other publicly available metagenomes (MG-RAST) showed that the Tibetan samples formed a distinct cluster (Figure 2B). The TP soil communities were distinct despite the fact that the soils share some physicochemical characteristics and land-use with soils from other environments such as arctic areas or farmland. Although phyla are detected in all compiled soils (including TP), the relative frequencies of bacterial phyla vary in different niches, which are generally thought to be controlled by local environmental variables including thawing, nitrogen gradients, pH gradients and other geochemical factors [11, 51]. TP has created a separate group from other environmental samples with similar phyla but different abundance (Figure 2A). Proteobacteria was the dominant phylum in high pH soil from shrub-covered permafrost and farmland . In metal-rich soils with low carbon and nitrogen, some genera of Alpha-proteobacteria, which are capable of nitrogen fixation, are present in high levels . In tundra soil, Actinobacteria was the dominant phylum due to its capacity to maintain metabolism and DNA repair at a low temperature and to degrade complex organic matter in soil [5, 52]. These observations are consistent with our findings in TP soils owing to the low-nitrogen and low-temperature of soils in these regions .
Soil samples from the TP shared a similar core microbiota compared with other environmental samples (Additional file 3: Figure S2). The core microbiota in samples from L1 differed from other TP samples, but the most abundant were relatively similar across all TP soil samples when compared with other environmental water, soil and sediment samples. These predominant genera are involved in nutrient element supplies, and phosphorus, sulfur, nitrogen cycles and these soil bacteria depend on and interact with each other while metabolizing the many raw materials present in soil. Gemmatimonas is a polyphosphate-accumulating genus for plant absorbing essential phosphorus in nutrient deficient TP areas . Thiobacillus has the ability to couple denitrification to sulfur-compound oxidation . Genera participating in the nitrogen cycle were also abundant in TP soil. Nitrospira, a nitrite-oxidizing bacteria, converts nitrite to nitrate as central step in the nitrogen cycle . Chthoniobacte r is involved in the transformation of organic carbon compounds in the soil, but does not have the capacity to fix nitrogen . Anaeromyxobacter can reduce nitrate and nitrite to ammonia, and nitrous oxide to nitrogen gas under both oxic and anoxic redox conditions . Opitutus species can produce acetate and propionate and hydrogen, which are important in methanogenesis, and can reduce nitrate to nitrite . These nitrogen cycle associated genera also obviously differed in abundance in upper and lower layer soil samples.
Alpha diversity of the bacterial communities from Namco Lake was lower than all farmland and meadow sites in TP due to the oligotrophic environment of lake. Farmlands investigated in this study had been converted from alpine meadows within a decade. Although anthropogenic disturbance of farmlands did not increase taxonomic alpha diversity, taxonomic beta diversity decreased significantly relative to alpine meadows (Figure 3A and 3B). A similar pattern has been reported in other studies of microbial communities in Netherlands pine forest soils and Amazon agriculture soils [4, 59]. It is worth noting that functional dissimilarities of alpine meadow sites were also higher than those of farmlands in both upper and lower layers of soil samples (Figure 3C). The metabolic pathways involved in amino acid synthesis, energy metabolism, vitamin and cofactor metabolism and lipid and glycan metabolism were considered to be responsible for bacterial adaptation to special environmental conditions . We found that the dissimilarities of these pathways in farmlands were lower than those of alpine meadows, indicating that human activities might drive the homogenization of bacterial communities and their functional gene pools in TP farmland soils (Figure 4). Similarly, dissimilarities in elemental metabolism cycles such as carbon, nitrogen, sulfur and methane among alpine meadows were also higher than those among farmlands, which further confirmed that the bacterial communities in alpine meadows were more diverse with respect to elements utilization and energy conversion without human disturbance. Loss of diversity with persistent human disturbance will increase the risk of abrupt and potentially irreversible ecosystem collapse . Homogenization of soil bacterial communities of farmlands converted from alpine meadow is also likely to make TP ecosystems vulnerable to environmental changes.
Differences between bacterial communities at depth of 3-5 and 18-20 cm in TP soil were mostly observed among root-associated bacteria and nitrifying bacteria (more abundant in upper soil than lower soil) and denitrifying bacteria (less abundant in upper than lower soil) (Figure 3D). Nitrifying bacteria prefer aerobic conditions, whereas denitrification is likely to occur in anaerobic environments . Some genera are indicators of different environmental factors. For example, Haliangium species are obligate halophiles  that differed in abundance between upper and lower layers of the alpine meadow G4. These genera varied between layers owing to variations in soil oxygen and pH. Taken together, these results highlight that human disturbance homogenized the bacterial communities in both composition and function at depths above 20 cm.
Bacterial community composition in TP samples was influenced by a number of geochemical factors including carbon, nitrogen, sulfur, temperature and moisture. In all PCA based on environmental factors, taxonomic composition, and functional composition, samples from three farmland sites (far or near residential area) were clustered together, and all other samples were dispersed, that they varied greatly (Figure 5A-C). CCA plots of sampling sites and environmental factors suggested that TS may control the taxonomic composition of G4 and G5 samples, while saline water saturation may affect L1 samples (Figure 5D). All remaining samples were positively correlated with most of the major and minor elements. Changes in pH value in the studied soil samples were not relevant to taxonomic composition. Viles  reported that vegetation changes, land use, and anthropogenic activities influenced soil microbial communities. Soil samples collected from four regions of Australia showed that their microbial community composition was primarily affected by geochemistry and secondarily by other factors . The bacterial community composition of F1 (far from residential area) and F2/3 (near residential area with fertilization) were nearly the same, indicating that these geochemical factors and land use both shaped bacterial communities.
Evaluation of the main element cycle and functions of bacterial communities revealed that the availability of nitrogen and sulfur and moisture of samples influenced the composition of different TP bacterial communities. Nitrogen fixation is very important for plant growth in the TP. Abundant bacteria in the Tibetan samples included Burkholderia, a β-proteobacteria genus associated with root-nodules. Root-nodules bacteria are capable of establishing a nitrogen-fixing symbiosis on legumes or other host plant roots . The moisture differences among the twenty Tibetan samples were negatively correlated with five Rhizobia genera (Figure 6). The abundance of Bradyrhizobium, Mesorhizobium, Rhizobium and Sinorhizobium influenced by moisture was also found in the north China Plain, but the relationship between Rhizobia and pH values differed between them [6, 64]. Due to climate warming and glacial melting, deposition of reactive nitrogen in soils is increasing . Accordingly, the impact of such environmental changes on soil microbial communities should be further investigated in alpine meadows associated with glaciers.
The sulfur cycle is another important biogeochemical cycle, because sulfur is a constituent of many proteins and cofactors . Thiobacillus was the dominant genus in alpine meadow G4 and G5, which had the highest sulfur compared to other samples. Thiobacillus is a chemolithoautotrophic sulfur-oxidizing β-Proteobacteria capable of nitrate-dependent Fe (II) and U (IV) oxidation [54, 67, 68]. The different number of reads assigned to various T. denitrificans protein-coding genes further reflected its genomic potential adaption to local environmental factors. The elevated sequencing depth coverage in sulfite reductase genes and cold-shock genes reflect the adaption to a high sulfur concentration and low temperature (Figure 6 and Additional file 4: Table S7). A higher frequency of NADH dehydrogenase genes, which catalyze the transfer of electrons from NADH to coenzyme Q, contributes to UV radiation resistance [69, 70] (Additional file 5: Figure 3).
This study presented the first comprehensive assessment of bacterial community structure and function in TP soils. However, considering the complexity of soil microbial community and the less exploration of TP soils, there were still a lot of difficulties in decoding these metagenomic data. The microbial community in soils especially in pristine environments is dominated by novel microorganisms, which are lack of closely related reference sequences in known public genomic databases. Homology based approaches, such as NR-BLAST method plus the lowest common ancestor (LCA) algorithm, tend to classify metagenomic sequences into higher taxonomic levels. A similar problem also occurs in functional annotation, which may lead to inaccurate annotation or the absence of some important and key functional genes. This is a universal problem that puzzles all metagenomic studies in microbial ecology. As the number of sequenced species grows, the situation will be improved.
TP soil bacterial community structure and functional profiles exhibit several distinct features from other environments, which are associated with the adaptation to environmental factors and land use. The beta diversity of bacterial communities in farmland decreased due to anthropogenic activities by homogenization even at different soil depths. As a much less explored environment, the metagenomic sequences and associated geochemical data generated in this study provide a foundation for further exploration of microbial resources and environmental adaptation in TP.
This research was supported by The Fundamental Research Funds for the Central Universities (2652012071 to XG), National Natural Science Foundation (40972162 to FL, 91131013 to FZ), and Funds for China Geological Survey (1212011121171). The authors would like to thank Cidan ZX, Yuxuan X, Minghan W and Yanming Z for their assistance with this study.
- Delmont TO, Prestat E, Keegan KP, Faubladier M, Robe P, Clark IM, Pelletier E, Hirsch PR, Meyer F, Gilbert JA: Structure, fluctuation and magnitude of a natural grassland soil metagenome. Isme J. 2012, 6 (9): 1677-1687. 10.1038/ismej.2011.197.PubMed CentralView ArticlePubMedGoogle Scholar
- Wallenstein MD, McMahon S, Schimel J: Bacterial and fungal community structure in Arctic tundra tussock and shrub soils. Fems Microbiol Ecol. 2007, 59 (2): 428-435. 10.1111/j.1574-6941.2006.00260.x.View ArticlePubMedGoogle Scholar
- Pointing SB, Chan Y, Lacap DC, Lau MC, Jurgens JA, Farrell RL: Highly specialized microbial diversity in hyper-arid polar desert. Proc Natl Acad Sci U S A. 2009, 106 (47): 19964-19969. 10.1073/pnas.0908274106.PubMed CentralView ArticlePubMedGoogle Scholar
- Kuramae EE, Yergeau E, Wong LC, Pijl AS, van Veen JA, Kowalchuk GA: Soil characteristics more strongly influence soil bacterial communities than land-use type. Fems Microbiol Ecol. 2012, 79 (1): 12-24. 10.1111/j.1574-6941.2011.01192.x.View ArticlePubMedGoogle Scholar
- Yergeau E, Hogues H, Whyte LG, Greer CW: The functional potential of high Arctic permafrost revealed by metagenomic sequencing, qPCR and microarray analyses. Isme J. 2010, 4 (9): 1206-1214. 10.1038/ismej.2010.41.View ArticlePubMedGoogle Scholar
- Fierer N, Jackson RB: The diversity and biogeography of soil bacterial communities. Proc Natl Acad Sci U S A. 2006, 103 (3): 626-631. 10.1073/pnas.0507535103.PubMed CentralView ArticlePubMedGoogle Scholar
- Nemergut DR, Costello EK, Hamady M, Lozupone C, Jiang L, Schmidt SK, Fierer N, Townsend AR, Cleveland CC, Stanish L: Global patterns in the biogeography of bacterial taxa. Environ Microbiol. 2011, 13 (1): 135-144. 10.1111/j.1462-2920.2010.02315.x.View ArticlePubMedGoogle Scholar
- Youssef NH, Elshahed MS: Diversity rankings among bacterial lineages in soil. Isme J. 2009, 3: 305-313. 10.1038/ismej.2008.106.View ArticlePubMedGoogle Scholar
- Reith F, Brugger J, Zammit CM, Gregg AL, Goldfarb KC, Andersen GL, Desantis TZ, Piceno YM, Brodie EL, Lu Z: Influence of geogenic factors on microbial communities in metallogenic Australian soils. Isme J. 2012, 6: 2107-2118. 10.1038/ismej.2012.48.PubMed CentralView ArticlePubMedGoogle Scholar
- Viles HA, Naylor LA, Carter NEA, Chaput D: Biogeomorphological disturbance regimes: progress in linking ecological and geornorphological systems. Earth Surf Proc Land. 2008, 33 (9): 1419-1435. 10.1002/esp.1717.View ArticleGoogle Scholar
- Fierer N, Lauber CL, Ramirez KS, Zaneveld J, Bradford MA, Knight R: Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. Isme J. 2012, 6 (5): 1007-1017. 10.1038/ismej.2011.159.PubMed CentralView ArticlePubMedGoogle Scholar
- Lauber CL, Hamady M, Knight R, Fierer N: Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microbiol. 2009, 75 (15): 5111-5120. 10.1128/AEM.00335-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Ni J: A simulation of biomes on the Tibetan Plateau and their responses to global climate change. Mt Res Dev. 2000, 20 (1): 80-89. 10.1659/0276-4741(2000)020[0080:ASOBOT]2.0.CO;2.View ArticleGoogle Scholar
- Kato T, Tang YH, Gu S, Hirota M, Du MY, Li YN, Zhao XQ: Temperature and biomass influences on interannual changes in CO2 exchange in an alpine meadow on the Qinghai-Tibetan Plateau. Global Change Biol. 2006, 12 (7): 1285-1298. 10.1111/j.1365-2486.2006.01153.x.View ArticleGoogle Scholar
- Qiu J: China: the third pole. Nature. 2008, 454 (7203): 393-396. 10.1038/454393a.View ArticlePubMedGoogle Scholar
- Liu Y, Yao T, Jiao N, Kang S, Xu B, Zeng Y, Huang S, Liu X: Bacterial diversity in the snow over Tibetan Plateau Glaciers. Extremophiles. 2009, 13 (3): 411-423. 10.1007/s00792-009-0227-5.View ArticlePubMedGoogle Scholar
- Jiang H, Dong H, Yu B, Liu X, Li Y, Ji S, Zhang CL: Microbial response to salinity change in Lake Chaka, a hypersaline lake on Tibetan plateau. Environ Microbiol. 2007, 9 (10): 2603-2621. 10.1111/j.1462-2920.2007.01377.x.View ArticlePubMedGoogle Scholar
- Wong FK, Lacap DC, Lau MC, Aitchison JC, Cowan DA, Pointing SB: Hypolithic microbial community of quartz pavement in the high-altitude tundra of central Tibet. Microb Ecol. 2010, 60 (4): 730-739. 10.1007/s00248-010-9653-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong FK, Lau MC, Lacap DC, Aitchison JC, Cowan DA, Pointing SB: Endolithic microbial colonization of limestone in a high-altitude arid environment. Microb Ecol. 2010, 59 (4): 689-699. 10.1007/s00248-009-9607-8.View ArticlePubMedGoogle Scholar
- Fierer N, Leff JW, Adams BJ, Nielsen UN, Bates ST, Lauber CL, Owens S, Gilbert JA, Wall DH, Caporaso JG: Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. Proc Natl Acad Sci U S A. 2012, 109 (52): 21390-21395. 10.1073/pnas.1215210110.PubMed CentralView ArticlePubMedGoogle Scholar
- Mackelprang R, Waldrop MP, DeAngelis KM, David MM, Chavarria KL, Blazewicz SJ, Rubin EM, Jansson JK: Metagenomic analysis of a permafrost microbial community reveals a rapid response to thaw. Nature. 2011, 480 (7377): 368-371. 10.1038/nature10576.View ArticlePubMedGoogle Scholar
- Gianoulis TA, Raes J, Patel PV, Bjornson R, Korbel JO, Letunic I, Yamada T, Paccanaro A, Jensen LJ, Snyder M: Quantifying environmental adaptation of metabolic pathways in metagenomics. Proc Natl Acad Sci U S A. 2009, 106 (5): 1374-1379. 10.1073/pnas.0808022106.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao GM, Tang YH, Mo WH, Wang YA, Li YN, Zhao XQ: 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.View ArticleGoogle Scholar
- Wang GX, Qian J, Cheng GD, Lai YM: Soil organic carbon pool of grassland soils on the Qinghai-Tibetan Plateau and its global implication. Sci Total Environ. 2002, 291 (1–3): 207-217.Google Scholar
- Zheng D, Zhang Q: Mountain Geoecology and Sustainable Development of the TibetanPlateau. 2000, Dordrecht, The Netherlands: Kluwer Academic PublishersView ArticleGoogle Scholar
- Fan JW, Zhong HP, Harris W, Yu GR, Wang SQ, Hu ZM, Yue YZ: Carbon storage in the grasslands of China based on field measurements of above- and below-ground biomass. Clim Change. 2008, 86 (3–4): 375-396.View ArticleGoogle Scholar
- Genxu W, Ju Q, Guodong C, Yuanmin L: Soil organic carbon pool of grassland soils on the Qinghai-Tibetan Plateau and its global implication. Sci Total Environ. 2002, 291 (1–3): 207-217.View ArticlePubMedGoogle Scholar
- Pan ZF, Deng GB, Zhai XG, Wu F, Yu MQ: Genetic diversity of Acid-PAGE monomeric prolamins in cultivated hulless barley (Hordeum vulgare L.) from Qinghai-Tibet plateau in China. Genet Resour Crop Ev. 2007, 54 (8): 1691-1699. 10.1007/s10722-006-9177-2.View ArticleGoogle Scholar
- Xu JQ, Yu SM, Liu JS, Haginoya S, Ishigooka Y, Kuwagata T, Hara M, Yasunari T: The implication of heat and water balance changes in a lake basin on the tibetan plateau. Hydrol Res Lett. 2009, 3: 1-5.View ArticleGoogle Scholar
- Wu YH, Zhu LP, Q H Y, L B W: The response of lake-glacier area change to climate variations in Namco Basin, Central Tibetan Plateau, during the last three decades. Acta Geograph Sin. 2007, 3: 301-311.Google Scholar
- Hallin S, Jones CM, Schloter M, Philippot L: Relationship between N-cycling communities and ecosystem functioning in a 50-year-old fertilization experiment. Isme J. 2009, 3 (5): 597-605. 10.1038/ismej.2008.128.View ArticlePubMedGoogle Scholar
- Bowers RM, McLetchie S, Knight R, Fierer N: Spatial variability in airborne bacterial communities across land-use types and their relationship to the bacterial communities of potential source environments. Isme J. 2011, 5 (4): 601-612. 10.1038/ismej.2010.167.PubMed CentralView ArticlePubMedGoogle Scholar
- Eilers KG, Debenport S, Anderson S, Fierer N: Digging deeper to find unique microbial communities: the strong effect of depth on the structure of bacterial and archaeal communities in soil. Soil Biol Biochem. 2012, 50: 58-65.View ArticleGoogle Scholar
- Islam KR, Weil RR: A rapid microwave digestion method for colorimetric measurement of soil organic carbon. Commun Soil Sci Plant Anal. 1998, 29: 2269-2284. 10.1080/00103629809370110.View ArticleGoogle Scholar
- Huson DH, Mitra S, Ruscheweyh HJ, Weber N, Schuster SC: Integrative analysis of environmental sequences using MEGAN4. Genome Res. 2011, 21 (9): 1552-1560. 10.1101/gr.120618.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A: The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinforma. 2008, 9: 386-10.1186/1471-2105-9-386.View ArticleGoogle Scholar
- Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, He G, Chen Y, Pan Q, Liu Y: SOAPdenovo2: An Empirically Improved Memory-Efficient Short-Read De Novo Assembler. 2012Google Scholar
- Zhao F, Zhao F, Li T, Bryant DA: A new pheromone trail-based genetic algorithm for comparative genome assembly. Nucleic Acids Res. 2008, 36 (10): 3455-3462. 10.1093/nar/gkn168.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu W, Lomsadze A, Borodovsky M: Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010, 38 (12): e132-10.1093/nar/gkq275.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang J, Qi J, Zhao H, He S, Zhang Y, Wei S, Zhao F: Metagenomic sequencing reveals microbiota and its functional potential associated with periodontal disease. Sci Rep. 2013, 3: 1843-PubMed CentralPubMedGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36: D480-D484.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu J, Wang H, Yang H, Zhang Y, Wang J, Zhao F, Qi J: Composition-based classification of short metagenomic sequences elucidates the landscapes of taxonomic and functional enrichment of microorganisms. Nucleic Acids Res. 2013, 41 (1): e3-10.1093/nar/gks828.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: Genome project data P: the sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Qi J, Zhao F, Buboltz A, Schuster SC: inGAP: an integrated next-generation genome analysis pipeline. Bioinformatics. 2010, 26 (1): 127-129. 10.1093/bioinformatics/btp615.PubMed CentralView ArticlePubMedGoogle Scholar
- Qi J, Zhao F: inGAP-sv: a novel scheme to identify and visualize structural variation from paired end mapping data. Nucleic Acids Res. 2011, 39: W567-W575. 10.1093/nar/gkr506.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang J, Shen J, Wu Y, Tu C, Soininen J, Stegen JC, He J, Liu X, Zhang L, Zhang E: Phylogenetic beta diversity in bacterial assemblages across ecosystems: deterministic versus stochastic processes. Isme J. 2013, 7: 1310-1321. 10.1038/ismej.2013.30.PubMed CentralView ArticlePubMedGoogle Scholar
- Grapov D, Newman JW: imDEV: a graphical user interface to R multivariate analysis tools in microsoft excel. Bioinformatics. 2012, 28 (17): 2288-2290. 10.1093/bioinformatics/bts439.PubMed CentralView ArticlePubMedGoogle Scholar
- Parton WJ, Stewart JWB, Cole CV: Dynamics of C, N, P and S in grassland soils: a model. Biogeochemistry. 1988, 5: 109-131. 10.1007/BF02180320.View ArticleGoogle Scholar
- Liu Y, Yao T, Jiao N, Kang S, Zeng Y, Huang S: Microbial community structure in moraine lakes and glacial meltwaters, Mount Everest. FEMS Microbiol Lett. 2006, 265 (1): 98-105. 10.1111/j.1574-6968.2006.00477.x.View ArticlePubMedGoogle Scholar
- Rousk J, Baath E, Brookes PC, Lauber CL, Lozupone C, Caporaso JG, Knight R, Fierer N: Soil bacterial and fungal communities across a pH gradient in an arable soil. Isme J. 2010, 4 (10): 1340-1351. 10.1038/ismej.2010.58.View ArticlePubMedGoogle Scholar
- Campbell BJ, Polson SW, Hanson TE, Mack MC, Schuur EA: The effect of nutrient deposition on bacterial communities in Arctic tundra soil. Environ Microbiol. 2010, 12 (7): 1842-1854. 10.1111/j.1462-2920.2010.02189.x.View ArticlePubMedGoogle Scholar
- Zhang H, Sekiguchi Y, Hanada S, Hugenholtz P, Kim H, Kamagata Y, Nakamura K: Gemmatimonas aurantiaca gen. nov., sp. nov., a gram-negative, aerobic, polyphosphate-accumulating micro-organism, the first cultured representative of the new bacterial phylum Gemmatimonadetes phyl. nov. Int J Syst Evol Microbiol. 2003, 53 (4): 1155-1163. 10.1099/ijs.0.02520-0.View ArticlePubMedGoogle Scholar
- Beller HR, Chain PS, Letain TE, Chakicherla A, Larimer FW, Richardson PM, Coleman MA, Wood AP, Kelly DP: The genome sequence of the obligately chemolithoautotrophic, facultatively anaerobic bacterium Thiobacillus denitrificans. J Bacteriol. 2006, 188 (4): 1473-1488. 10.1128/JB.188.4.1473-1488.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Off S, Alawi M, Spieck E: Enrichment and physiological characterization of a novel Nitrospira-like bacterium obtained from a marine sponge. Appl Environ Microbiol. 2010, 76 (14): 4640-4646. 10.1128/AEM.00320-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Kant R, van Passel MW, Palva A, Lucas S, Lapidus A, Glavina del Rio T, Dalin E, Tice H, Bruce D, Goodwin L: Genome sequence of Chthoniobacter flavus Ellin428, an aerobic heterotrophic soil bacterium. J Bacteriol. 2011, 193 (11): 2902-2903. 10.1128/JB.00295-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanford RA, Cole JR, Tiedje JM: Characterization and description of Anaeromyxobacter dehalogenans gen. nov., sp. nov., an aryl-halorespiring facultative anaerobic myxobacterium. Appl Environ Microbiol. 2002, 68 (2): 893-900. 10.1128/AEM.68.2.893-900.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Chin KJ, Janssen PH: Propionate formation by Opitutus terrae in pure culture and in mixed culture with a hydrogenotrophic methanogen and implications for carbon fluxes in anoxic rice paddy soil. Appl Environ Microbiol. 2002, 68 (4): 2089-2092. 10.1128/AEM.68.4.2089-2092.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodrigues JL, Pellizari VH, Mueller R, Baek K, Jesus Eda C, Paula FS, Mirza B, Hamaoui GS, Tsai SM, Feigl B: Conversion of the Amazon rainforest to agriculture results in biotic homogenization of soil bacterial communities. Proc Natl Acad Sci U S A. 2013, 110 (3): 988-993. 10.1073/pnas.1220608110.PubMed CentralView ArticlePubMedGoogle Scholar
- MacDougall AS, McCann KS, Gellner G, Turkington R: Diversity loss with persistent human disturbance increases vulnerability to ecosystem collapse. Nature. 2013, 494 (7435): 86-89. 10.1038/nature11869.View ArticlePubMedGoogle Scholar
- Galloway JN, Townsend AR, Erisman JW, Bekunda M, Cai Z, Freney JR, Martinelli LA, Seitzinger SP, Sutton MA: Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science. 2008, 320 (5878): 889-892. 10.1126/science.1136674.View ArticlePubMedGoogle Scholar
- Zhang YQ, Li YZ, Wang B, Wu ZH, Zhang CY, Gong X, Qiu ZJ, Zhang Y: Characteristics and living patterns of marine myxobacterial isolates. Appl Environ Microbiol. 2005, 71 (6): 3331-3336. 10.1128/AEM.71.6.3331-3336.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez-Romero E: Coevolution in Rhizobium-legume symbiosis?. DNA Cell Biol. 2009, 28 (8): 361-370. 10.1089/dna.2009.0863.View ArticlePubMedGoogle Scholar
- Zhang YM, Li Y, Chen WF, Wang ET, Tian CF, Li QQ, Zhang YZ, Sui XH, Chen WX: Biodiversity and biogeography of rhizobia associated with soybean plants grown in the North China Plain. Appl Environ Microbiol. 2011, 77 (18): 6331-6342. 10.1128/AEM.00542-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Saros JE, Rose KC, Clow DW, Stephens VC, Nurse AB, Arnett HA, Stone JR, Williamson CE, Wolfe AP: Melting Alpine glaciers enrich high-elevation lakes with reactive nitrogen. Environ Sci Technol. 2010, 44 (13): 4891-4896. 10.1021/es100147j.View ArticlePubMedGoogle Scholar
- Klimmek O: The biological cycle of sulfur. Met Ions Biol Syst. 2005, 43: 105-130.View ArticlePubMedGoogle Scholar
- Beller HR: Anaerobic, nitrate-dependent oxidation of U(IV) oxide minerals by the chemolithoautotrophic bacterium Thiobacillus denitrificans. Appl Environ Microbiol. 2005, 71 (4): 2170-2174. 10.1128/AEM.71.4.2170-2174.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Straub KL, Benz M, Schink B, Widdel F: Anaerobic, nitrate-dependent microbial oxidation of ferrous iron. Appl Environ Microbiol. 1996, 62 (4): 1458-1460.PubMed CentralPubMedGoogle Scholar
- Hoppe U, Bergemann J, Diembeck W, Ennen J, Gohla S, Harris I, Jacob J, Kielholz J, Mei W, Pollet D: Coenzyme Q10, a cutaneous antioxidant and energizer. Biofactors. 1999, 9 (2–4): 371-378.View ArticlePubMedGoogle Scholar
- Dunlap WC, Fujisawa A, Yamamoto Y: UV radiation increases the reduced coenzyme Q ratio in marine bacteria. Redox Rep. 2002, 7 (5): 320-323. 10.1179/135100002125000893.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.