Extensive Unexplored Human Microbiome Diversity Revealed by Over 150,000 Genomes from Metagenomes Spanning Age, Geography, and Lifestyle.

The body-wide human microbiome plays a role in health, but its full diversity remains uncharacterized, particularly outside of the gut and in international populations. We leveraged 9,428 metagenomes to reconstruct 154,723 microbial genomes (45% of high quality) spanning body sites, ages, countries, and lifestyles. We recapitulated 4,930 species-level genome bins (SGBs), 77% without genomes in public repositories (unknown SGBs [uSGBs]). uSGBs are prevalent (in 93% of well-assembled samples), expand underrepresented phyla, and are enriched in non-Westernized populations (40% of the total SGBs). We annotated 2.85 M genes in SGBs, many associated with conditions including infant development (94,000) or Westernization (106,000). SGBs and uSGBs permit deeper microbiome analyses and increase the average mappability of metagenomic reads from 67.76% to 87.51% in the gut (median 94.26%) and 65.14% to 82.34% in the mouth. We thus identify thousands of microbial genomes from yet-to-be-named species, expand the pangenomes of human-associated microbes, and allow better exploitation of metagenomic technologies.

Macroevolutionary Dynamics in Micro-organisms: Generalists Give Rise to Specialists Across Biomes in the Ubiquitous Bacterial Phylum Myxococcota.

Prokaryotes dominate the Tree of Life, but our understanding of the macroevolutionary processes generating this diversity is still limited. Habitat transitions are thought to be a key driver of prokaryote diversity. However, relatively little is known about how prokaryotes successfully transition and persist across environments, and how these processes might vary between biomes and lineages. Here, we investigate biome transitions and specialization in natural populations of a focal bacterial phylum, the Myxococcota, sampled across a range of replicated soils and freshwater and marine sediments in Cornwall (UK). By targeted deep sequencing of the protein-coding gene rpoB, we found >2,000 unique Myxococcota lineages, with the majority (77%) classified as biome specialists and with only <5% of lineages distributed across the salt barrier. Discrete character evolution models revealed that specialists in one biome rarely transitioned into specialists in another biome. Instead, evolved generalism mediated transitions between biome specialists. State-dependent diversification models found variation in speciation rates across the tree, but this variation was independent of biome association or specialization. Our findings were robust to phylogenetic uncertainty, different levels of species delineation, and different assumed amounts of unsampled diversity resulting in an incomplete phylogeny. Overall, our results are consistent with a "jack-of-all-trades" tradeoff where generalists suffer a cost in any individual environment, resulting in rapid evolution of niche specialists and shed light on how bacteria could transition between biomes.

Swarm v3: towards tera-scale amplicon clustering.

MOTIVATION: Previously we presented swarm, an open-source amplicon clustering programme that produces fine-scale molecular operational taxonomic units (OTUs) that are free of arbitrary global clustering thresholds. Here, we present swarm v3 to address issues of contemporary datasets that are growing towards tera-byte sizes. RESULTS: When compared with previous swarm versions, swarm v3 has modernized C++ source code, reduced memory footprint by up to 50%, optimized CPU-usage and multithreading (more than 7 times faster with default parameters), and it has been extensively tested for its robustness and logic. AVAILABILITY AND IMPLEMENTATION: Source code and binaries are available at https://github.com/torognes/swarm. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Treatment of Active Crohn's Disease With an Ordinary Food-based Diet That Replicates Exclusive Enteral Nutrition.

BACKGROUND & AIMS: Exclusive enteral nutrition (EEN) is the only established dietary treatment for Crohn's disease (CD), but its acceptability is limited. There is a need for novel dietary treatments for CD. METHODS: We evaluated the effects of an individualized food-based diet (CD-TREAT), with similar composition to EEN, on the gut microbiome, inflammation, and clinical response in a rat model, healthy adults, and children with relapsing CD. Twenty-five healthy adults randomly received EEN or CD-TREAT for 7 days, followed by a 14-day washout period, followed by the alternate diet. Fecal microbiome and metabolome were assessed before and after each diet. HLA-B7 and HLA-B27 transgenic rats with gut inflammation received EEN, CD-TREAT, or standard chow for 4 weeks. Fecal, luminal, and tissue microbiome, fecal metabolites, and gut inflammation were assessed. Five children with active CD activity received CD-TREAT and their clinical activity and calprotectin were evaluated after 8 weeks of treatment. RESULTS: For healthy adults, CD-TREAT was easier to comply with and more acceptable than EEN. CD-TREAT induced similar effects to EEN (EEN vs CD-TREAT) on fecal microbiome composition, metabolome, mean total sulfide (increase 133.0 ± 80.5 vs 54.3 ± 47.0 nmol/g), pH (increase 1.3 ± 0.5 vs 0.9 ± 0.6), and the short-chain fatty acids (µmol/g) acetate (decrease 27.4 ± 22.6 vs 21.6 ± 20.4), propionate (decrease 5.7 ± 7.8 vs 5.2 ± 7.9), and butyrate (decrease 7.0 ± 7.4 vs 10.2 ± 8.5). In the rat model, CD-TREAT and EEN produced similar changes in bacterial load (decrease 0.3 ± 0.3 log10 16S rRNA gene copies per gram), short-chain fatty acids, microbiome, and ileitis severity (mean histopathology score decreases of 1.25 for EEN [P = .015] and 1.0 for CD-TREAT [P = .044] vs chow). In children receiving CD-TREAT, 4 (80%) had a clinical response and 3 (60%) entered remission, with significant concurrent decreases in fecal calprotectin (mean decrease 918 ± 555 mg/kg; P = .002). CONCLUSION: CD-TREAT replicates EEN changes in the microbiome, decreases gut inflammation, is well tolerated, and is potentially effective in patients with active CD. ClinicalTrials.gov, numbers NCT02426567 and NCT03171246.

Critical Assessment of Metagenome Interpretation-a benchmark of metagenomics software.

Methods for assembly, taxonomic profiling and binning are key to interpreting metagenome data, but a lack of consensus about benchmarking complicates performance assessment. The Critical Assessment of Metagenome Interpretation (CAMI) challenge has engaged the global developer community to benchmark their programs on highly complex and realistic data sets, generated from ∼700 newly sequenced microorganisms and ∼600 novel viruses and plasmids and representing common experimental setups. Assembly and genome binning programs performed well for species represented by individual genomes but were substantially affected by the presence of related strains. Taxonomic profiling and binning programs were proficient at high taxonomic ranks, with a notable performance decrease below family level. Parameter settings markedly affected performance, underscoring their importance for program reproducibility. The CAMI results highlight current challenges but also provide a roadmap for software selection to answer specific research questions.

Binning metagenomic contigs by coverage and composition.

Shotgun sequencing enables the reconstruction of genomes from complex microbial communities, but because assembly does not reconstruct entire genomes, it is necessary to bin genome fragments. Here we present CONCOCT, a new algorithm that combines sequence composition and coverage across multiple samples, to automatically cluster contigs into genomes. We demonstrate high recall and precision on artificial as well as real human gut metagenome data sets.

Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform.

With read lengths of currently up to 2 × 300 bp, high throughput and low sequencing costs Illumina's MiSeq is becoming one of the most utilized sequencing platforms worldwide. The platform is manageable and affordable even for smaller labs. This enables quick turnaround on a broad range of applications such as targeted gene sequencing, metagenomics, small genome sequencing and clinical molecular diagnostics. However, Illumina error profiles are still poorly understood and programs are therefore not designed for the idiosyncrasies of Illumina data. A better knowledge of the error patterns is essential for sequence analysis and vital if we are to draw valid conclusions. Studying true genetic variation in a population sample is fundamental for understanding diseases, evolution and origin. We conducted a large study on the error patterns for the MiSeq based on 16S rRNA amplicon sequencing data. We tested state-of-the-art library preparation methods for amplicon sequencing and showed that the library preparation method and the choice of primers are the most significant sources of bias and cause distinct error patterns. Furthermore we tested the efficiency of various error correction strategies and identified quality trimming (Sickle) combined with error correction (BayesHammer) followed by read overlapping (PANDAseq) as the most successful approach, reducing substitution error rates on average by 93%.

Fluvial network organization imprints on microbial co-occurrence networks.

Recent studies highlight linkages among the architecture of ecological networks, their persistence facing environmental disturbance, and the related patterns of biodiversity. A hitherto unresolved question is whether the structure of the landscape inhabited by organisms leaves an imprint on their ecological networks. We analyzed, based on pyrosequencing profiling of the biofilm communities in 114 streams, how features inherent to fluvial networks affect the co-occurrence networks that the microorganisms form in these biofilms. Our findings suggest that hydrology and metacommunity dynamics, both changing predictably across fluvial networks, affect the fragmentation of the microbial co-occurrence networks throughout the fluvial network. The loss of taxa from co-occurrence networks demonstrates that the removal of gatekeepers disproportionately contributed to network fragmentation, which has potential implications for the functions biofilms fulfill in stream ecosystems. Our findings are critical because of increased anthropogenic pressures deteriorating stream ecosystem integrity and biodiversity.

Illumina error profiles: resolving fine-scale variation in metagenomic sequencing data.

BACKGROUND: Illumina's sequencing platforms are currently the most utilised sequencing systems worldwide. The technology has rapidly evolved over recent years and provides high throughput at low costs with increasing read-lengths and true paired-end reads. However, data from any sequencing technology contains noise and our understanding of the peculiarities and sequencing errors encountered in Illumina data has lagged behind this rapid development. RESULTS: We conducted a systematic investigation of errors and biases in Illumina data based on the largest collection of in vitro metagenomic data sets to date. We evaluated the Genome Analyzer II, HiSeq and MiSeq and tested state-of-the-art low input library preparation methods. Analysing in vitro metagenomic sequencing data allowed us to determine biases directly associated with the actual sequencing process. The position- and nucleotide-specific analysis revealed a substantial bias related to motifs (3mers preceding errors) ending in "GG". On average the top three motifs were linked to 16 % of all substitution errors. Furthermore, a preferential incorporation of ddGTPs was recorded. We hypothesise that all of these biases are related to the engineered polymerase and ddNTPs which are intrinsic to any sequencing-by-synthesis method. We show that quality-score-based error removal strategies can on average remove 69 % of the substitution errors - however, the motif-bias remains. CONCLUSION: Single-nucleotide polymorphism changes in bacterial genomes can cause significant changes in phenotype, including antibiotic resistance and virulence, detecting them within metagenomes is therefore vital. Current error removal techniques are not designed to target the peculiarities encountered in Illumina sequencing data and other sequencing-by-synthesis methods, causing biases to persist and potentially affect any conclusions drawn from the data. In order to develop effective diagnostic and therapeutic approaches we need to be able to identify systematic sequencing errors and distinguish these errors from true genetic variation.

A comprehensive benchmarking study of protocols and sequencing platforms for 16S rRNA community profiling.

BACKGROUND: In the last 5 years, the rapid pace of innovations and improvements in sequencing technologies has completely changed the landscape of metagenomic and metagenetic experiments. Therefore, it is critical to benchmark the various methodologies for interrogating the composition of microbial communities, so that we can assess their strengths and limitations. The most common phylogenetic marker for microbial community diversity studies is the 16S ribosomal RNA gene and in the last 10 years the field has moved from sequencing a small number of amplicons and samples to more complex studies where thousands of samples and multiple different gene regions are interrogated. RESULTS: We assembled 2 synthetic communities with an even (EM) and uneven (UM) distribution of archaeal and bacterial strains and species, as metagenomic control material, to assess performance of different experimental strategies. The 2 synthetic communities were used in this study, to highlight the limitations and the advantages of the leading sequencing platforms: MiSeq (Illumina), The Pacific Biosciences RSII, 454 GS-FLX/+ (Roche), and IonTorrent (Life Technologies). We describe an extensive survey based on synthetic communities using 3 experimental designs (fusion primers, universal tailed tag, ligated adaptors) across the 9 hypervariable 16S rDNA regions. We demonstrate that library preparation methodology can affect data interpretation due to different error and chimera rates generated during the procedure. The observed community composition was always biased, to a degree that depended on the platform, sequenced region and primer choice. However, crucially, our analysis suggests that 16S rRNA sequencing is still quantitative, in that relative changes in abundance of taxa between samples can be recovered, despite these biases. CONCLUSION: We have assessed a range of experimental conditions across several next generation sequencing platforms using the most up-to-date configurations. We propose that the choice of sequencing platform and experimental design needs to be taken into consideration in the early stage of a project by running a small trial consisting of several hypervariable regions to quantify the discriminatory power of each region. We also suggest that the use of a synthetic community as a positive control would be beneficial to identify the potential biases and procedural drawbacks that may lead to data misinterpretation. The results of this study will serve as a guideline for making decisions on which experimental condition and sequencing platform to consider to achieve the best microbial profiling.

The skin microbiome in psoriatic arthritis: methodology development and pilot data.

BACKGROUND: Skin microbiota are likely to be important in the development of conditions such as psoriatic arthritis. Profiling the bacterial community in the psosriatic plaques will contribute to our understanding of the role of the skin microbiome in these conditions. The aim of this work was to determine the optimum study design for work on the skin microbiome with use of the MiSeq platform. The objectives were to compare data generated from two platforms for two primer pairs in a low density mock bacterial community. METHODS: DNA was obtained from two low density mock communities of 11 diverse bacterial strains (with and without human DNA supplementation) and from swabs taken from the skin of four healthy volunteers. The DNA was amplified with primer pairs covering hypervariable regions of the 16S rRNA gene: primers 63F and 519R (V1-V3), and 347F and 803R (V3-V4). The resultant libraries were indexed for the MiSeq and Roche454 platforms and sequenced. Both datasets were de-noised, cleaned of chimeras, and analysed by use of QIIME software (version 1.8.0). FINDINGS: No significant difference in the diversity indices at the phylum and the genus level between the platforms was seen. Comparison of the diversity indices for the mock community data for the two primer pairs demonstrated that the V3-V4 hypervariable region had significantly better capture of bacterial diversity than did the V1-V3 region. Amplification with the same primer pairs showed strong concordance within each platform (98·9-99·8%), with negligible effect of spiked human DNA contamination. Comparison at the family level classification between samples processed on the MiSeq and Roche454 platforms using the V3-V4 hypervariable region also showed a high level of concordance (87%), although less so for the V1-V3 primers (10%). The pilot data from healthy volunteers were similar. INTERPRETATION: Results obtained from the V3-V4 16S rRNA hypervariable region, sequencing on the MiSeq and Roche454 platforms, were concordant between replicates, and between each other. These findings suggest that the MiSeq platform, and these primers, is a comparable method for determining skin microbiota to the widely used Roche454 methodology. FUNDING: NIHR Manchester Musculoskeletal Biomedical Research Unit.

Benchmarking of viral haplotype reconstruction programmes: an overview of the capacities and limitations of currently available programmes.

Viral haplotype reconstruction from a set of observed reads is one of the most challenging problems in bioinformatics today. Next-generation sequencing technologies enable us to detect single-nucleotide polymorphisms (SNPs) of haplotypes-even if the haplotypes appear at low frequencies. However, there are two major problems. First, we need to distinguish real SNPs from sequencing errors. Second, we need to determine which SNPs occur on the same haplotype, which cannot be inferred from the reads if the distance between SNPs on a haplotype exceeds the read length. We conducted an independent benchmarking study that directly compares the currently available viral haplotype reconstruction programmes. We also present nine in silico data sets that we generated to reflect biologically plausible populations. For these data sets, we simulated 454 and Illumina reads and applied the programmes to test their capacity to reconstruct whole genomes and individual genes. We developed a novel statistical framework to demonstrate the strengths and limitations of the programmes. Our benchmarking demonstrated that all the programmes we tested performed poorly when sequence divergence was low and failed to recover haplotype populations with rare haplotypes.

Metabarcoding and metabolome analyses of copepod grazing reveal feeding preference and linkage to metabolite classes in dynamic microbial plankton communities.

In order to characterize copepod feeding in relation to microbial plankton community dynamics, we combined metabarcoding and metabolome analyses during a 22-day seawater mesocosm experiment. Nutrient amendment of mesocosms promoted the development of haptophyte (Phaeocystis pouchetii)- and diatom (Skeletonema marinoi)-dominated plankton communities in mesocosms, in which Calanus sp. copepods were incubated for 24 h in flow-through chambers to allow access to prey particles (<500 µm). Copepods and mesocosm water sampled six times spanning the experiment were analysed using metabarcoding, while intracellular metabolite profiles of mesocosm plankton communities were generated for all experimental days. Taxon-specific metabarcoding ratios (ratio of consumed prey to available prey in the surrounding seawater) revealed diverse and dynamic copepod feeding selection, with positive selection on large diatoms, heterotrophic nanoflagellates and fungi, while smaller phytoplankton, including P. pouchetii, were passively consumed or even negatively selected according to our indicator. Our analysis of the relationship between Calanus grazing ratios and intracellular metabolite profiles indicates the importance of carbohydrates and lipids in plankton succession and copepod-prey interactions. This molecular characterization of Calanus sp. grazing therefore provides new evidence for selective feeding in mixed plankton assemblages and corroborates previous findings that copepod grazing may be coupled to the developmental and metabolic stage of the entire prey community rather than to individual prey abundances.

Bioaugmentation Mitigates the Impact of Estrogen on Coliform-Grazing Protozoa in Slow Sand Filters.

Exposure to endocrine-disrupting chemicals (EDCs), such as estrogens, is a growing issue for human and animal health as they have been shown to cause reproductive and developmental abnormalities in wildlife and plants and have been linked to male infertility disorders in humans. Intensive farming and weather events, such as storms, flash flooding, and landslides, contribute estrogen to waterways used to supply drinking water. This paper explores the impact of estrogen exposure on the performance of slow sand filters (SSFs) used for water treatment. The feasibility and efficacy of SSF bioaugmentation with estrogen-degrading bacteria was also investigated, to determine whether removal of natural estrogens (estrone, estradiol, and estriol) and overall SSF performance for drinking water treatment could be improved. Strains for SSF augmentation were isolated from full-scale, municipal SSFs so as to optimize survival in the laboratory-scale SSFs used. Concentrations of the natural estrogens, determined by gas chromatography coupled with mass spectrometry (GC-MS), revealed augmented SSFs reduced the overall estrogenic potency of the supplied water by 25% on average and removed significantly more estrone and estradiol than nonaugmented filters. A negative correlation was found between coliform removal and estrogen concentration in nonaugmented filters. This was due to the toxic inhibition of protozoa, indicating that high estrogen concentrations can have functional implications for SSFs (such as impairing coliform removal). Consequently, we suggest that high estrogen concentrations could impact significantly on water quality production and, in particular, on pathogen removal in biological water filters.

Evolution of the plankton paleome in the Black Sea from the Deglacial to Anthropocene.

The complex interplay of climate shifts over Eurasia and global sea level changes modulates freshwater and saltwater inputs to the Black Sea. The dynamics of the hydrologic changes from the Late Glacial into the Holocene remain a matter of debate, and information on how these changes affected the ecology of the Black Sea is sparse. Here we used Roche 454 next-generation pyrosequencing of sedimentary 18S rRNA genes to reconstruct the plankton community structure in the Black Sea over the last ca. 11,400 y. We found that 150 of 2,710 species showed a statistically significant response to four environmental stages. Freshwater chlorophytes were the best indicator species for lacustrine conditions (>9.0 ka B.P.), although the copresence of previously unidentified marine taxa indicated that the Black Sea might have been influenced to some extent by the Marmara Sea since at least 9.6 ka calendar (cal) B.P. Dinoflagellates, cercozoa, eustigmatophytes, and haptophytes responded most dramatically to the gradual increase in salinity after the latest marine reconnection and during the warm and moist mid-Holocene climatic optimum. According to paired analysis of deuterium/hydrogen (D/H) isotope ratios in fossil alkenones, salinity increased rapidly with the onset of the dry Subboreal after ~5.2 ka B.P., leading to an increase in marine fungi and the first occurrence of marine copepods. A gradual succession of dinoflagellates, diatoms, and chrysophytes occurred during the refreshening after ~2.5 ka cal B.P. with the onset of the cool and wet Subatlantic climate and recent anthropogenic perturbations.

Extensive Modulation of the Fecal Metagenome in Children With Crohn's Disease During Exclusive Enteral Nutrition.

OBJECTIVES: Exploring associations between the gut microbiota and colonic inflammation and assessing sequential changes during exclusive enteral nutrition (EEN) may offer clues into the microbial origins of Crohn's disease (CD). METHODS: Fecal samples (n=117) were collected from 23 CD and 21 healthy children. From CD children fecal samples were collected before, during EEN, and when patients returned to their habitual diets. Microbiota composition and functional capacity were characterized using sequencing of the 16S rRNA gene and shotgun metagenomics. RESULTS: Microbial diversity was lower in CD than controls before EEN (P=0.006); differences were observed in 36 genera, 141 operational taxonomic units (OTUs), and 44 oligotypes. During EEN, the microbial diversity of CD children further decreased, and the community structure became even more dissimilar than that of controls. Every 10 days on EEN, 0.6 genus diversity equivalents were lost; 34 genera decreased and one increased during EEN. Fecal calprotectin correlated with 35 OTUs, 14 of which accounted for 78% of its variation. OTUs that correlated positively or negatively with calprotectin decreased during EEN. The microbiota of CD patients had a broader functional capacity than healthy controls, but diversity decreased with EEN. Genes involved in membrane transport, sulfur reduction, and nutrient biosynthesis differed between patients and controls. The abundance of genes involved in biotin (P=0.005) and thiamine biosynthesis decreased (P=0.017), whereas those involved in spermidine/putrescine biosynthesis (P=0.031), or the shikimate pathway (P=0.058), increased during EEN. CONCLUSIONS: Disease improvement following treatment with EEN is associated with extensive modulation of the gut microbiome.

Decreased gut microbiota diversity, delayed Bacteroidetes colonisation and reduced Th1 responses in infants delivered by caesarean section.

OBJECTIVE: The early intestinal microbiota exerts important stimuli for immune development, and a reduced microbial exposure as well as caesarean section (CS) has been associated with the development of allergic disease. Here we address how microbiota development in infants is affected by mode of delivery, and relate differences in colonisation patterns to the maturation of a balanced Th1/Th2 immune response. DESIGN: The postnatal intestinal colonisation pattern was investigated in 24 infants, born vaginally (15) or by CS (nine). The intestinal microbiota were characterised using pyrosequencing of 16S rRNA genes at 1 week and 1, 3, 6, 12 and 24 months after birth. Venous blood levels of Th1- and Th2-associated chemokines were measured at 6, 12 and 24 months. RESULTS: Infants born through CS had lower total microbiota diversity during the first 2 years of life. CS delivered infants also had a lower abundance and diversity of the Bacteroidetes phylum and were less often colonised with the Bacteroidetes phylum. Infants born through CS had significantly lower levels of the Th1-associated chemokines CXCL10 and CXCL11 in blood. CONCLUSIONS: CS was associated with a lower total microbial diversity, delayed colonisation of the Bacteroidetes phylum and reduced Th1 responses during the first 2 years of life.

Rare but active taxa contribute to community dynamics of benthic biofilms in glacier-fed streams.

Glaciers harbour diverse microorganisms, which upon ice melt can be released downstream. In glacier-fed streams microorganisms can attach to stones or sediments to form benthic biofilms. We used 454-pyrosequencing to explore the bulk (16S rDNA) and putatively active (16S rRNA) microbial communities of stone and sediment biofilms across 26 glacier-fed streams. We found differences in community composition between bulk and active communities among streams and a stronger congruence between biofilm types. Relative abundances of rRNA and rDNA were positively correlated across different taxa and taxonomic levels, but at lower taxonomic levels, the higher abundance in either the active or the bulk communities became more apparent. Here, environmental variables played a minor role in structuring active communities. However, we found a large number of rare taxa with higher relative abundances in rRNA compared with rDNA. This suggests that rare taxa contribute disproportionately to microbial community dynamics in glacier-fed streams. Our findings propose that high community turnover, where taxa repeatedly enter and leave the 'seed bank', contributes to the maintenance of microbial biodiversity in harsh ecosystems with continuous environmental perturbations, such as glacier-fed streams.

Loss of microbial diversity in soils is coincident with reductions in some specialized functions.

Loss of microbial diversity is considered a major threat because of its importance for ecosystem functions, but there is a lack of conclusive evidence that diversity itself is reduced under anthropogenic stress, and about the consequences of diversity loss. Heavy metals are one of the largest, widespread pollutant types globally, and these represent a significant environmental stressor for terrestrial microbial communities. Using combined metagenomics and functional assays, we show that the compositional and functional response of microbial communities to long-term heavy metal stress results in a significant loss of diversity. Our results indicate that even at a moderate loss of diversity, some key specialized functions (carried out by specific groups) may be compromised. Together with previous work, our data suggest disproportionate impact of contamination on microbes that carry out specialized, but essential, ecosystem functions. Based on these findings, we propose a conceptual framework to explicitly consider diversity of functions and microbial functional groups to test the relationship between biodiversity and soil functions.