Comparison of microbial taxonomic and functional shift pattern along contamination gradient
© The Author(s). 2016
Received: 20 December 2015
Accepted: 6 June 2016
Published: 14 June 2016
The interaction mechanism between microbial communities and environment is a key issue in microbial ecology. Microbial communities usually change significantly under environmental stress, which has been studied both phylogenetically and functionally, however which method is more effective in assessing the relationship between microbial communities shift and environmental changes still remains controversial.
By comparing the microbial taxonomic and functional shift pattern along heavy metal contamination gradient, we found that both sedimentary composition and function shifted significantly along contamination gradient. For example, the relative abundance of Geobacter and Fusibacter decreased along contamination gradient (from high to low), while Janthinobacterium and Arthrobacter increased their abundances. Most genes involved in heavy metal resistance (e.g., metc, aoxb and mer) showed higher intensity in sites with higher concentration of heavy metals. Comparing the two shift patterns, there were correlations between them, because functional and phylogenetic β-diversities were significantly correlated, and many heavy metal resistance genes were derived from Geobacter, explaining their high abundance in heavily contaminated sites. However, there was a stronger link between functional composition and environmental drivers, while stochasticity played an important role in formation and succession of phylogenetic composition demonstrated by null model test.
Overall our research suggested that the responses of functional traits depended more on environmental changes, while stochasticity played an important role in formation and succession of phylogenetic composition for microbial communities. So profiling microbial functional composition seems more appropriate to study the relationship between microbial communities and environment, as well as explore the adaptation and remediation mechanism of microbial communities to heavy metal contamination.
Investigating microbial community composition, structure and function, as well as their response to environmental changes are key issues in microbial ecology. Nowadays, microbial communities have been studied both phylogenetically and functionally. Of them, 16S rRNA gene amplicons sequencing is the most effective and popular method to profile microbial community composition, and many uncultured species of nature environments have been identified [1, 2], but it provides little information about the metabolic potential of these species. GeoChip is a comprehensive microarray for investigating microbial community function, also provides direct linkages of microbial genes to ecosystem processes and functions . However, the relationship between community composition and function still remains controversial, which also makes it difficult to evaluate the interactions between microorganisms and environment [4, 5]. For example, a study of soil microbial communities in forest showed that microbial community structure and functional gene diversity were significantly linked to soil and plant factors, but some environmental factors correlated to community structure were not related to its function, such as temperature and tree diameter . Exploring the relationship between microbial community structure and function would help explaining such inconsistency.
Basically, microbial phylogeny and function are strongly correlated, and phylogenetic trees based on 16S usually closely resemble clusters obtained on the basis of shared gene content. Therefore, researchers often infer properties of uncultured organisms from cultured relatives . For example, using 16S information, PICRUSt recaptures key findings from the Human Microbiome Project and accurately predicts the abundance of gene families in host-associated and environmental communities . However, widespread cases of horizontal gene transfer impair the consistency between taxonomic and functional composition. And it was recently suggested that functional traits are more valuable to understand bacterial community assembly and explain shifts in microbial community composition across environmental gradients [9, 10], probably because the adaptation and ecological niche of microbial population was more directly determined by their functional composition.
Microbial communities usually change extensively under environmental stress, which provides us an opportunity to explore the correlation and difference between microbial phylogenetic and functional shift pattern. Recently, heavy metals such as mercury (Hg), chromium (Cr), lead (Pb), manganese (Mn), and arsenic (As) have induced serious diseases or even death of organisms through contaminated waters or soils, although heavy metals in trace amount are beneficial even significant to organisms [11–13]. The investigation of microbial populations distribution along heavy metal contamination gradients across spatial scales will help elucidate how natural communities respond to environmental changes . Many studies on the issue have been focused on functional and phylogenetical analyses. For example, in highly heavy metal contaminated sites, an overall lower gene diversity but higher abundance for specific functional genes, such as heavy metal homeostasis genes and sulfate-reducing genes were observed [11, 15, 16], and the dominant microbial groups included α-Proteobacteria, β-Proteobacteria and Firmicutes [17, 18].
In this study, we attempted to address two core issues (i) the taxonomic and functional shift pattern of sedimentary microbial communities to heavy metal contamination; and (ii) correlation and difference among the two shift patterns. To explore the taxonomic and functional response of microbial communities to heavy metal contamination, 12 sedimentary samples were taken from three sites in the Xiangjiang River with a gradient of contaminant levels (described before ), and analyzed by GeoChip 5.0 and 16S rRNA gene amplicons sequencing. The study provides us an insight into the shift pattern of microbial communities to heavy metal contamination, and demonstrates that functional profiling microbial communities is more effective in examining the interaction between microorganisms and environments.
Samples were collected from sediment of Xiangjiang River (Hunan, China), as previous described . In this study, we choose three groups of samples with different distance from drain outlet, 500 m, 1000 m, and 1500 m, separately. Geochemical properties of each sample were measured. The composition of heavy metals including Hg, As, Cr, Pb, Mn, cobalt (Co), cadmium (Cd), nickel (Ni), copper (Cu) and zinc (Zn) in the sediments was analyzed by ICP-AES . Total sedimentary organic nitrogen (N) was quantified by Kjeldahl distillation . The amount of total sedimentary organic carbon (C) was analyzed by potassium dichromate oxidation-ferrous sulphate titrimetry .
Illumina sequencing, GeoChip analysis and data processing
DNA was extracted using a TIANamp Bacterial DNA Kit (MO BIO Laboratories, Inc., Carlsbad, CA). The V4 region of the 16S rRNA genes was amplified with the primer pair 515 F (5’-GTGCCAGCMGCCGCGGTAA-3’) and 806R (5’- GGACTACHVGGGTWTCTAAT-3’). Sample libraries were generated from purified PCR products. The MiSeq 500 cycles kit was used for 2x250 bp paired-ends sequencing on MiSeq machine (Illumina, San Diego, CA). Sequences with perfect matches to barcodes were split to sample libraries, and trimmed. OTU clustering was performed through UCLUST at 97 % similarity level , and taxonomic assignment was through the RDP classifier  with a minimal 50 % confidence estimate. The above steps were conducted through the Galaxy pipeline (http://zhoulab5.rccc.ou.edu/) developed by Qin el al. Subsequent analyses were performed in R . Finally, samples were rarefied at 13,000 sequences per sample. All the 16S rRNA sequences were deposited in GenBank database and the accession number were KP784842 - KP788032.
For each sample, microbial community DNA was extracted and purified as described previously [15, 26]. Amplified DNA was labeled and hybridized with GeoChip 5.0, which is a powerful tool to study the functional diversity, composition, structure and metabolic potential of microbial communities . All GeoChip 5.0 hybridization data are available at the Institute for Environmental Genomics, University of Oklahoma (http://ieg.ou.edu/). The hybridized GeoChip 5.0 was analyzed as previously described . Software TMEV was used for hierarchical cluster analysis of sequencing and GeoChip data. Statistical differences between the functional microbial communities from the different sites were analyzed by analysis of variance (ANOVA).
Partial least squares path modeling (PLSPM) is a powerful structural equation modeling technique, which is used to elucidate the complex relationship among microbial community composition, structure and function of three groups of samples. Before performed in R v. 2.6.1 with the package plspm , principal component analysis (PCA) was conducted for 16S rRNA gene sequencing data, GeoChip hybridization data and environmental data respectively. Then PC1 and PC2 values were used for PLSPM. And α-diversity value could be used for PLSPM directly. Taxonomic composition and functional gene diversity was calculated using Shannon-Weiner’s H′ and evenness. Difference among three groups of microbial communities in composition and function was evaluated using dissimilarity test respectively. Mantel test was used to calculate correlations between microbial community diversity and environmental attributes . Null model analysis which assumes that a community is not structured by species interactions, was performed according to the method described by Zhou et al . In order to determine whether species/gene compositional differences among sites were caused by the forces causing communities to be different from the expectations by random chance or not, the permutational analysis of multivariate dispersions (PERMDISP) was used to test the significance of the differences of the sedimentary microbial communities of each group from null model expectations . All the analyses were performed online (http://ieg.ou.edu/).
Geochemical properties of samples were significantly (p < 0.05) different among three groups (Additional file 1: Table S1). The sample group nearest (500 m) to drain outlet had the highest concentration of heavy metals and organic C (1.51 ppm) and N (16.68 ppm), so we defined the group as Group H. Almost all the heavy metals detected in this study were most abundant in Group H, except for Pb and Zn. On the contrary, the sampling site farthest (1500 m) from drain outlet had the lowest concentration of heavy metals, thus was defined as Group L. The concentrations of Hg (0.18 ppm), Cd (3.0 ppm), Cu (34 ppm) and Zn (158 ppm) were the lowest in Group L. And the sample group with a sampling site between Group H and Group L was defined as Group M, which was about 1000 m away from the drain outlet. Group M had a moderate content of heavy metals, with similar Pb (103 and 124 ppm) and Zn (346 and 496 ppm) content to Group H, but similar As (70 and 73 ppm), Co (12.1 and 12.8 ppm), Cr (57 and 62 ppm), Ni (30.2 and 32.7 ppm) and Mn (788 and 1476 ppm) content to Group L.
Shift patterns of taxonomic and functional community composition
Dissimilarity test of three groups of microbial communities in phylogenetic composition
Dissimilarity test of three groups of microbial communities in functional composition
Comparison of taxonomic and functional shift pattern
Generally, a significant correlation (r = 0.711, p = 0.001) between phylogenetic β-diversity and functional β-diversity was observed. Although not all the microbial populations detected by 16S rRNA gene sequencing had their genes detected in GeoChip, there were still 10 microbial genera (relative abundance > 1 %) whose genes were detected by GeoChip. For example, 71 gene variants of Geobacter, 24 gene variants of Janthinobacterium and 211 gene variants of Arthrobacter were detected. More importantly, of all gene variants belonging to Geobacter 36.62 % of them were involved in heavy metal resistance, but only 4.17 % of Janthinobacterium genes were related to heavy metal resistance, indicating their different potential in adapting heavy metal contamination (Additional file 1: Table S3).
Significance test of the differences of centroids between the sedimentary microbial communities and null model simulations across different contamination gradient
Centroid of actual communities
Centroid of the null model
The interaction mechanism between microbial communities and environment is a key issue in microbial ecology. Microbial communities usually change significantly under environmental stress, which provides us an opportunity to study the microorganism-environment interaction mechanism. Here, by investigating how sedimentary microbial communities shifted along heavy metal contamination gradient both phylogenetically and functionally, we aimed to identify an effective method to profile the relationship between microbial communities and environment, and reveal the adaptation mechanism of sedimentary microbial communities to heavy metal contamination.
Taxonomic and functional shift pattern along contamination gradient
Generally, both phylogenetic and functional composition of microbial communities shifted significantly along contamination gradient. Previous studies showed that heavy metals would decrease the diversity of microbial community [15, 31], and we observed the same tendency as microbial communities had the lowest Shannon diversity index in H Group. All kinds of microbial populations were clustered into three major categories, according to their abundance pattern in three sites. They included microbial genera positively related to heavy metal contamination, negatively related to heavy metal contamination and showed no obvious change along contamination gradient. For example, Fusibacter and Geobacter were positively correlated to contamination, while Janthinobacterium showed a negative correlation. A larger percent of Fusibacter (35.21 %) were detected in Group H, probably because many Fusibacter species could reduce sulfur or thiosulfate [32, 33], which usually are of high metal resistance and play important roles in heavy metal bioremediation. Geobacter were more abundant in Group H and Group L, mainly resulting from their metal ion reduction ability and potential for use in bioremediation of radioactive metals [34, 35]. On the contrary, Janthinobacterium is an important genus of Betaproteobacteria , and it is found to be susceptible to heavy metals, such as Ag, Cu, Hg, Pb and Ni . In this study, the relative abundance of Janthinobacterium decreased significantly in highly contaminated samples, consistent with previous studies.
Meanwhile, functional composition of microbial communities also shifted along contamination gradient. It has been reported that prolonged exposure to high concentrations of heavy metals resulted in a significant loss of metabolic diversity . Although a decrease in metabolic diversity hasn’t been detected in heavily contaminated sites in this study, which might resulted from their long-term adaptation to contaminated environment , lots of genes involved in carbon and nitrogen metabolism showed a decrease in intensity, such as nifh and amyA, involved in nitrogen fixation and carbon degradation respectively. The results were supported by previous studies in which heavy metal pollution produced a dramatic reduce in nutrient metabolism [40, 41]. However, microorganisms have several mechanisms to survive heavy metal contamination, such as discharging toxic metals and enzymatic conversion [42, 43]. In this way, high intensity of heavy metal resistance genes was usually observed in heavy contaminated sites. Here, we detected higher abundance of heavy metal resistance genes (e.g., metc, aoxb, mer, merb and silaffin) in Group H while lower abundance in Group L, which was a vital part in the functional shift process of sedimentary microbial communities.
Taxonomic and functional compositions were correlated
Furthermore, we compared the phylogenetic and functional shift pattern in order to identify an effective method to study the interactions between microbial communities and environments as well as the adaptation mechanism of microorganisms to heavy metal contamination. Comparing the two shift patterns, correlation between them is one important point, because microbial populations respond to environment through their functional genes. For example, most of microorganisms in AMD are capable of oxidizing Fe or S, so we could approximately predict the Fe or S oxidation potential of a community based on 16S rRNA gene sequencing data . Here, genes related to heavy metal resistance were derived mainly from Actinobacteria, Gammaproteobacteria, Alphaproteobacteria and Betaproteobacteria, which were also the main classes in sedimentary microbial communities. However, not all the gene variants detected in GeoChip could be assigned to the specific microbial populations identified by sequencing, because probes in microarray are usually designed before knowing the microbial community composition . Nevertheless, the different potential in heavy metal resistance of Geobacter and Janthinobacterium detected by Geochip, was consistent with their abundance pattern in three contamination sites, demonstrating the correlation between taxonomic and functional shift pattern of microbial populations. Because of the limitation of GeoChip, a metagenomic insight into the microbial communities is needed to profile the relationship between microbial community composition and function in future study. Moreover, we found that functional and phylogenetic β-diversities were significantly (r = 0.711, p = 0.001) correlated, which was also reported in a previous study . We have acknowledged that microbial species closely related in phylogenetic trees based on 16S usually have similar functional composition, and researchers could almost accurately predict the abundance of gene families in environmental communities based on 16S information [7, 8]. In this way, evidences supported that phylogenetic and functional composition are correlated to a certain extent.
Functional composition was more environmentally dependent
However, since a broad range of functional variation may occur among closely related organisms, taxonomic distributions are assumed to be ambiguous in assessing the response of microbial communities to environmental changes . And horizontal gene transfer might be the other key factor resulting in the divergence between phylogenetic trees based on 16S information and functional composition . A recent study has documented that specific functions could be widely detected across a variety of taxa or phylogenetic groups . For example, sulphate-reducing bacterial species are widely distributed in various phyla, such as Deltaproteobacteria, Nitrospirae, Clostridia and Euryarchaeota . Such differences lead to their different relationship with environment. Researchers reported that functional traits were more reliable in assessing the relationship between microbial communities and ecological processes, as microbial assemblages were better predicted at the functional genes level rather than at taxonomic level . In this study, both mantel test and PLSPM analysis indicated that the responses of functional traits might depend more on environmental changes.
Historically contingent in taxonomic composition indicated that environmental conditions would determine the types of ecological niches available for specific functional groups, while species compositions with similar physiological fitness are stochastically influenced by the history . Particularly, a previous study showed that stochastic processes played important roles in controlling the assembly and succession of the groundwater microbial community . Therefore, we speculated that selection strength, mainly heavy metal contamination stress in this study, shaped and directed the functional shift pattern of sedimentary microbial communities, but their phylogenetic composition had various shift patterns to achieve the same function shift, because similar function genes are widely distributed. For example, heavy metal resistance related genes in this study were derived from 29 microbial phyla, including 57 classes, 107 orders and 502 genera. So each microbial population capable of heavy metal resistance had a chance to become more abundant in heavily contaminated sites in theory. Supporting the hypothesis, we demonstrated that stochasticity played a more significant role in phylogenetic composition than in functional composition of sedimentary microbial communities using null model test.
Collectively, by comparing the microbial taxonomic and functional shift pattern along heavy metal contamination gradient, the study demonstrated that: (i) the responses of functional traits depended more on environmental changes, and stochasticity played an important role in formation and succession of phylogenetic composition; (ii) taxonomic composition and functional composition were closely correlated, although taxonomically related populations neither sufficiently nor necessarily meant functional similarity. The study is of high significance in future metagenomic research, and also provides us an insight into the adaptation pattern of microbial communities to heavy metal contamination.
ANOVA, analysis of variance; As, arsenic; C, carbon; Cd, cadmium; Co, cobalt; Cr, chromium; Cu, copper; Gof, goodness of fit; Hg, mercury; Mn, manganesev; N, nitrogen; Ni, nickel; P, phosphorus; Pb, lead; PCA, principal component analysis; PERMDISP, permutational analysis of multivariate dispersions; PLSPM, partial least squares path modeling; S, sulfur; Zn, zinc
We are grateful to Ye Deng and the data analysis platform of Institute for Environmental Genomics.
This research was supported by Natural Science Foundation of China (41573072).
Availability of data and materials
All the 16S rRNA sequences were deposited in GenBank database and the accession number were KP784842 - KP788032.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
YR and JN: Data analysis and manuscript preparation. YR and HY: Design and conception of the experiments. YX and XZ: Execution of the experiments. WH, DP, YL, XL and HY: Participated in discussions and supervised overall experimental and theoretical works. All authors have read and approved the manuscript.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kim O-S, Cho Y-J, Lee K, Yoon S-H, Kim M, Na H, et al. Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species. Int J Syst Evol Microbiol. 2012;62(Pt 3):716–21.View ArticlePubMedGoogle Scholar
- Schmidt TM, DeLong E, Pace N. Analysis of a marine picoplankton community by 16S rRNA gene cloning and sequencing. J Bacteriol. 1991;173(14):4371–8.PubMedPubMed CentralGoogle Scholar
- He Z, Gentry TJ, Schadt CW, Wu L, Liebich J, Chong SC, et al. GeoChip: a comprehensive microarray for investigating biogeochemical, ecological and environmental processes. ISME J. 2007;1(1):67–77.View ArticlePubMedGoogle Scholar
- Fukami T, Martijn Bezemer T, Mortimer SR, Putten WH. Species divergence and trait convergence in experimental plant community assembly. Ecol Lett. 2005;8(12):1283–90.View ArticleGoogle Scholar
- Devictor V, Mouillot D, Meynard C, Jiguet F, Thuiller W, Mouquet N. Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world. Ecol Lett. 2010;13(8):1030–40.PubMedGoogle Scholar
- Zhang Y, Cong J, Lu H, Yang C, Yang Y, Zhou J, et al. An integrated study to analyze soil microbial community structure and metabolic potential in two forest types. PLoS One. 2014;9(4), e93773.View ArticlePubMedPubMed CentralGoogle Scholar
- Zaneveld JR, Lozupone C, Gordon JI, Knight R. Ribosomal RNA diversity predicts genome diversity in gut bacteria and their relatives. Nucleic Acids Res. 2010;38(12):3869–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Barberán A, Fernández‐Guerra A, Bohannan BJ, Casamayor EO. Exploration of community traits as ecological markers in microbial metagenomes. Mol Ecol. 2012;21(8):1909–17.View ArticlePubMedGoogle Scholar
- Edwards KF, Litchman E, Klausmeier CA. Functional traits explain phytoplankton responses to environmental gradients across lakes of the United States. Ecology. 2013;94(7):1626–35.View ArticlePubMedGoogle Scholar
- Waldron PJ, Wu L, Nostrand JDV, Schadt CW, He Z, Watson DB, et al. Functional gene array-based analysis of microbial community structure in groundwaters with a gradient of contaminant levels. Environ Sci Technol. 2009;43(10):3529–34.View ArticlePubMedGoogle Scholar
- Li J, Richter DD, Mendoza A, Heine P. Four-decade responses of soil trace elements to an aggrading old-field forest: B, Mn, Zn, Cu, and Fe. Ecology. 2008;89(10):2911–23.View ArticlePubMedGoogle Scholar
- Schlesinger WH, Cole JJ, Finzi AC, Holland EA. Introduction to coupled biogeochemical cycles. Front Ecol Environ. 2011;9(1):5–8.View ArticleGoogle Scholar
- Feris KP, Ramsey PW, Frazar C, Rillig M, Moore JN, Gannon JE, et al. Seasonal dynamics of shallow-hyporheic-zone microbial community structure along a heavy-metal contamination gradient. Appl Environ Microbiol. 2004;70(4):2323–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Kang S, Van Nostrand JD, Gough HL, He Z, Hazen TC, Stahl DA, et al. Functional gene array-based analysis of microbial communities in heavy metals-contaminated lake sediments. FEMS Microbiol Ecol. 2013;86(2):200–14.View ArticlePubMedGoogle Scholar
- Zhang X, Niu J, Liang Y, Liu X, Yin H. Metagenome-scale analysis yields insights into the structure and function of microbial communities in a copper bioleaching heap. BMC Genet. 2016;17(1):1.View ArticleGoogle Scholar
- Zhu J, Zhang J, Li Q, Han T, Xie J, Hu Y, et al. Phylogenetic analysis of bacterial community composition in sediment contaminated with multiple heavy metals from the Xiangjiang River in China. Mar Pollut Bull. 2013;70(1–2):134–9.View ArticlePubMedGoogle Scholar
- Xiao Y, Liu X, Ma L, Liang Y, Niu J, Gu Y, et al. Microbial communities from different subsystems in biological heap leaching system play different roles in iron and sulfur metabolisms. Appl Microbiol Biotechnol. 2016;99(23):1–10.Google Scholar
- Yin H, Niu J, Ren Y, Cong J, Zhang X, Fan F, et al. An integrated insight into the response of sedimentary microbial communities to heavy metal contamination. Scientific Reports. 2015;5:14266.View ArticlePubMedPubMed CentralGoogle Scholar
- Ramsey MH, Thompson M. High-accuracy analysis by inductively coupled plasma atomic emission spectrometry using the parameter-related internal standard method. J Anal At Spectrom. 1987;2(5):497–502.View ArticleGoogle Scholar
- McKenzie HA, Wallace HS. The Kjeldahl determination of nitrogen: a critical study of digestion conditions-temperature, catalyst, and oxidizing agent. Aust J Chem. 1954;7(1):55–70.View ArticleGoogle Scholar
- Liu C, Zhu Z, He X, Zhang B, Xia N. Rapid determination of organic carbon in marine sediment samples by potassium dichromate oxidation-ferrous sulphate titrimetry. Rock Miner Anal. 2007;6(3):205–8.Google Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.View ArticlePubMedGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Statistical Package R. R: A language of environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2009.Google Scholar
- Xiao Y, Xu Y, Dong W, Liang Y, Fan F, Zhang X, et al. The complicated substrates enhance the microbial diversity and zinc leaching efficiency in sphalerite bioleaching system. Appl Microbiol Biotechnol. 2015;99(23):10311–22.View ArticlePubMedGoogle Scholar
- Liang Y, He Z, Wu L, Deng Y, Li G, Zhou J. Development of a common oligonucleotide reference standard for microarray data normalization and comparison across different microbial communities. Appl Environ Microbiol. 2010;76(4):1088.View ArticlePubMedGoogle Scholar
- Smouse PE, Long JC, Sokal RR. Multiple regression and correlation extensions of the Mantel test of matrix correspondence. Syst Zool. 1986;35(4):627–32.View ArticleGoogle Scholar
- Zhou J, Deng Y, Zhang P, Xue K, Liang Y, Van Nostrand JD, et al. Stochasticity, succession, and environmental perturbations in a fluidic ecosystem. Proc Natl Acad Sci U S A. 2014;111(9):E836–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Chase JM. Stochastic community assembly causes higher biodiversity in more productive environments. Science. 2010;328(5984):1388–91.View ArticlePubMedGoogle Scholar
- Golebiewski M, Deja-Sikora E, Cichosz M, Tretyn A, Wrobel B. 16S rDNA pyrosequencing analysis of bacterial community in heavy metals polluted soils. Microb Ecol. 2014;67(3):635–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Fadhlaoui K, Hania WB, Postec A, Fauque G, Hamdi M, Ollivier B, et al. Fusibacter fontis sp. nov., a sulfur-reducing anaerobic bacterium, isolated from a mesothermic tunisian spring. Int J Syst Evol Microbiol. 2015;65:3501–6.View ArticlePubMedGoogle Scholar
- Ravot G, Magot M, Fardeau M-L, Patel BK, Thomas P, Garcia J-L, et al. Fusibacter paucivorans gen. nov., sp. nov., an anaerobic, thiosulfate-reducing bacterium from an oil-producing well. Int J Syst Bacteriol. 1999;49(3):1141–7.View ArticlePubMedGoogle Scholar
- Methe B, Nelson KE, Eisen J, Paulsen I, Nelson W, Heidelberg J, et al. Genome of Geobacter sulfurreducens: metal reduction in subsurface environments. Science. 2003;302(5652):1967–9.View ArticlePubMedGoogle Scholar
- Lovley DR, Giovannoni SJ, White DC, Champine JE, Phillips E, Gorby YA, et al. Geobacter metallireducens gen. nov. sp. nov., a microorganism capable of coupling the complete oxidation of organic compounds to the reduction of iron and other metals. Arch Microbiol. 1993;159(4):336–44.View ArticlePubMedGoogle Scholar
- Lincoln SP, Fermor TR, Tindall BJ. Janthinobacterium agaricidamnosum sp.nov.,a soft rot pathogen of Agaricus bisporus. Int Syst Bacteriol. 1999;49(4):1577–89.View ArticleGoogle Scholar
- Shi Y-L, Lu X-Z, Yu W-G. A new Beta-agarase from marine bacterium Janthinobacterium sp. SY12. Microbiol Biotechnol. 2008;24:2659–64.View ArticleGoogle Scholar
- Hemme CL, Deng Y, Gentry TJ, Fields MW, Wu L, Barua S, et al. Metagenomic insights into evolution of a heavy metal-contaminated groundwater microbial community. ISME J. 2010;4(5):660–72.View ArticlePubMedGoogle Scholar
- Chai L, Wang Z, Wang Y, Yang Z, Wang H, Wu X. Ingestion risks of metals in groundwater based on TIN model and dose-response assessment – a case study in the Xiangjiang watershed, central-south China. Total Environ. 2010;408(16):3118–24.View ArticleGoogle Scholar
- Llorens N, Arola L, Bladé C, Mas A. Effects of copper exposure upon nitrogen metabolism in tissue cultured Vitis vinifera. Plant Sci. 2000;160(1):159–63.View ArticlePubMedGoogle Scholar
- Boularbah A, Schwartz C, Bitton G, Morel JL. Heavy metal contamination from mining sites in South Morocco: 1. Use of a biotest to assess metal toxicity of tailings and soils. Chemosphere. 2006;63(5):802–10.View ArticlePubMedGoogle Scholar
- Valls M, Lorenzo VD. Exploiting the genetic and biochemical capacities of bacteria for the remediation of heavy metal pollution. FEMS Microbiol Rev. 2002;26:327–38.View ArticlePubMedGoogle Scholar
- Nies DH. Efflux-mediated heavy metal resistance in prokaryotes. FEMS Microbiol Rev. 2003;27(2–3):313–39.View ArticlePubMedGoogle Scholar
- Chen L-X, Hu M, Huang L-N, Hua Z-S, Kuang J-L, Li S-J, et al. Comparative metagenomic and metatranscriptomic analyses of microbial communities in acid mine drainage. The ISME journal. 2015;9(7):1579–92.View ArticlePubMedGoogle Scholar
- Kuang J, Huang L, He Z, Chen L, Hua Z, Jia P, et al. Predicting taxonomic and functional structure of microbial communities in acid mine drainage. The ISME journal. 2016.Google Scholar
- Green JL, Bohannan BJ, Whitaker RJ. Microbial biogeography: from taxonomy to traits. Science. 2008;320(5879):1039–43.View ArticlePubMedGoogle Scholar
- Gogarten JP, Townsend JP. Horizontal gene transfer, genome innovation and evolution. Nat Rev Microbiol. 2005;3(9):679–87.View ArticlePubMedGoogle Scholar
- Burke C, Steinberg P, Rusch D, Kjelleberg S, Thomas T. Bacterial community assembly based on functional genes rather than species. Proc Natl Acad Sci. 2011;108(34):14288–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Muyzer G, Stams AJ. The ecology and biotechnology of sulphate-reducing bacteria. Nat Rev Microbiol. 2008;6(6):441–54.PubMedGoogle Scholar