Interpreting 16S metagenomic data without clustering to achieve sub-OTU resolution

The standard approach to analyzing 16S tag sequence data, which relies on clustering reads by sequence similarity into Operational Taxonomic Units (OTUs), underexploits the accuracy of modern sequencing technology. We present a clustering-free approach to multi-sample Illumina datasets that can identify independent bacterial subpopulations regardless of the similarity of their 16S tag sequences. Using published data from a longitudinal time-series study of human tongue microbiota, we are able to resolve within standard 97% similarity OTUs up to 20 distinct subpopulations, all ecologically distinct but with 16S tags differing by as little as 1 nucleotide (99.2% similarity). A comparative analysis of oral communities of two cohabiting individuals reveals that most such subpopulations are shared between the two communities at 100% sequence identity, and that dynamical similarity between subpopulations in one host is strongly predictive of dynamical similarity between the same subpopulations in the other host. Our method can also be applied to samples collected in cross-sectional studies and can be used with the 454 sequencing platform. We discuss how the sub-OTU resolution of our approach can provide new insight into factors shaping community assembly.Comment: Updated to match the published version. 12 pages, 5 figures + supplement. Significantly revised for clarity, references added, results not change

Metagenomics reveals sediment microbial community response to Deepwater Horizon oil spill

The Deepwater Horizon (DWH) oil spill in the spring of 2010 resulted in an input of ∼4.1 million barrels of oil to the Gulf of Mexico; >22% of this oil is unaccounted for, with unknown environmental consequences. Here we investigated the impact of oil deposition on microbial communities in surface sediments collected at 64 sites by targeted sequencing of 16S rRNA genes, shotgun metagenomic sequencing of 14 of these samples and mineralization experiments using (14)C-labeled model substrates. The 16S rRNA gene data indicated that the most heavily oil-impacted sediments were enriched in an uncultured Gammaproteobacterium and a Colwellia species, both of which were highly similar to sequences in the DWH deep-sea hydrocarbon plume. The primary drivers in structuring the microbial community were nitrogen and hydrocarbons. Annotation of unassembled metagenomic data revealed the most abundant hydrocarbon degradation pathway encoded genes involved in degrading aliphatic and simple aromatics via butane monooxygenase. The activity of key hydrocarbon degradation pathways by sediment microbes was confirmed by determining the mineralization of (14)C-labeled model substrates in the following order: propylene glycol, dodecane, toluene and phenanthrene. Further, analysis of metagenomic sequence data revealed an increase in abundance of genes involved in denitrification pathways in samples that exceeded the Environmental Protection Agency (EPA)'s benchmarks for polycyclic aromatic hydrocarbons (PAHs) compared with those that did not. Importantly, these data demonstrate that the indigenous sediment microbiota contributed an important ecosystem service for remediation of oil in the Gulf. However, PAHs were more recalcitrant to degradation, and their persistence could have deleterious impacts on the sediment ecosystem

Development of the preterm gut microbiome in twins at risk of necrotising enterocolitis and sepsis

The preterm gut microbiome is a complex dynamic community influenced by genetic and environmental factors and is implicated in the pathogenesis of necrotising enterocolitis (NEC) and sepsis. We aimed to explore the longitudinal development of the gut microbiome in preterm twins to determine how shared environmental and genetic factors may influence temporal changes and compared this to the expressed breast milk (EBM) microbiome. Stool samples (n = 173) from 27 infants (12 twin pairs and 1 triplet set) and EBM (n = 18) from 4 mothers were collected longitudinally. All samples underwent PCR-DGGE (denaturing gradient gel electrophoresis) analysis and a selected subset underwent 454 pyrosequencing. Stool and EBM shared a core microbiome dominated by Enterobacteriaceae, Enterococcaceae, and Staphylococcaceae. The gut microbiome showed greater similarity between siblings compared to unrelated individuals. Pyrosequencing revealed a reduction in diversity and increasing dominance of Escherichia sp. preceding NEC that was not observed in the healthy twin. Antibiotic treatment had a substantial effect on the gut microbiome, reducing Escherichia sp. and increasing other Enterobacteriaceae. This study demonstrates related preterm twins share similar gut microbiome development, even within the complex environment of neonatal intensive care. This is likely a result of shared genetic and immunomodulatory factors as well as exposure to the same maternal microbiome during birth, skin contact and exposure to EBM. Environmental factors including antibiotic exposure and feeding are additional significant determinants of community structure, regardless of host genetics

Lactobacillus rhamnosus GG-supplemented formula expands butyrate-producing bacterial strains in food allergic infants.

Dietary intervention with extensively hydrolyzed casein formula supplemented with Lactobacillus rhamnosus GG (EHCF+LGG) accelerates tolerance acquisition in infants with cow's milk allergy (CMA). We examined whether this effect is attributable, at least in part, to an influence on the gut microbiota. Fecal samples from healthy controls (n=20) and from CMA infants (n=19) before and after treatment with EHCF with (n=12) and without (n=7) supplementation with LGG were compared by 16S rRNA-based operational taxonomic unit clustering and oligotyping. Differential feature selection and generalized linear model fitting revealed that the CMA infants have a diverse gut microbial community structure dominated by Lachnospiraceae (20.5±9.7%) and Ruminococcaceae (16.2±9.1%). Blautia, Roseburia and Coprococcus were significantly enriched following treatment with EHCF and LGG, but only one genus, Oscillospira, was significantly different between infants that became tolerant and those that remained allergic. However, most tolerant infants showed a significant increase in fecal butyrate levels, and those taxa that were significantly enriched in these samples, Blautia and Roseburia, exhibited specific strain-level demarcations between tolerant and allergic infants. Our data suggest that EHCF+LGG promotes tolerance in infants with CMA, in part, by influencing the strain-level bacterial community structure of the infant gut

Microbiome profiling by Illumina sequencing of combinatorial sequence-tagged PCR products

We developed a low-cost, high-throughput microbiome profiling method that uses combinatorial sequence tags attached to PCR primers that amplify the rRNA V6 region. Amplified PCR products are sequenced using an Illumina paired-end protocol to generate millions of overlapping reads. Combinatorial sequence tagging can be used to examine hundreds of samples with far fewer primers than is required when sequence tags are incorporated at only a single end. The number of reads generated permitted saturating or near-saturating analysis of samples of the vaginal microbiome. The large number of reads al- lowed an in-depth analysis of errors, and we found that PCR-induced errors composed the vast majority of non-organism derived species variants, an ob- servation that has significant implications for sequence clustering of similar high-throughput data. We show that the short reads are sufficient to assign organisms to the genus or species level in most cases. We suggest that this method will be useful for the deep sequencing of any short nucleotide region that is taxonomically informative; these include the V3, V5 regions of the bac- terial 16S rRNA genes and the eukaryotic V9 region that is gaining popularity for sampling protist diversity.Comment: 28 pages, 13 figure

Transkingdom Networks: A Systems Biology Approach to Identify Causal Members of Host-Microbiota Interactions

Improvements in sequencing technologies and reduced experimental costs have resulted in a vast number of studies generating high-throughput data. Although the number of methods to analyze these "omics" data has also increased, computational complexity and lack of documentation hinder researchers from analyzing their high-throughput data to its true potential. In this chapter we detail our data-driven, transkingdom network (TransNet) analysis protocol to integrate and interrogate multi-omics data. This systems biology approach has allowed us to successfully identify important causal relationships between different taxonomic kingdoms (e.g. mammals and microbes) using diverse types of data

Prospecting environmental mycobacteria: combined molecular approaches reveal unprecedented diversity

Background: Environmental mycobacteria (EM) include species commonly found in various terrestrial and aquatic environments, encompassing animal and human pathogens in addition to saprophytes. Approximately 150 EM species can be separated into fast and slow growers based on sequence and copy number differences of their 16S rRNA genes. Cultivation methods are not appropriate for diversity studies; few studies have investigated EM diversity in soil despite their importance as potential reservoirs of pathogens and their hypothesized role in masking or blocking M. bovis BCG vaccine. Methods: We report here the development, optimization and validation of molecular assays targeting the 16S rRNA gene to assess diversity and prevalence of fast and slow growing EM in representative soils from semi tropical and temperate areas. New primer sets were designed also to target uniquely slow growing mycobacteria and used with PCR-DGGE, tag-encoded Titanium amplicon pyrosequencing and quantitative PCR. Results: PCR-DGGE and pyrosequencing provided a consensus of EM diversity; for example, a high abundance of pyrosequencing reads and DGGE bands corresponded to M. moriokaense, M. colombiense and M. riyadhense. As expected pyrosequencing provided more comprehensive information; additional prevalent species included M. chlorophenolicum, M. neglectum, M. gordonae, M. aemonae. Prevalence of the total Mycobacterium genus in the soil samples ranged from 2.3×107 to 2.7×108 gene targets g−1; slow growers prevalence from 2.9×105 to 1.2×107 cells g−1. Conclusions: This combined molecular approach enabled an unprecedented qualitative and quantitative assessment of EM across soil samples. Good concordance was found between methods and the bioinformatics analysis was validated by random resampling. Sequences from most pathogenic groups associated with slow growth were identified in extenso in all soils tested with a specific assay, allowing to unmask them from the Mycobacterium whole genus, in which, as minority members, they would have remained undetected

Forensic analysis of the microbiome of phones and shoes

BACKGROUND: Microbial interaction between human-associated objects and the environments we inhabit may have forensic implications, and the extent to which microbes are shared between individuals inhabiting the same space may be relevant to human health and disease transmission. In this study, two participants sampled the front and back of their cell phones, four different locations on the soles of their shoes, and the floor beneath them every waking hour over a 2-day period. A further 89 participants took individual samples of their shoes and phones at three different scientific conferences. RESULTS: Samples taken from different surface types maintained significantly different microbial community structures. The impact of the floor microbial community on that of the shoe environments was strong and immediate, as evidenced by Procrustes analysis of shoe replicates and significant correlation between shoe and floor samples taken at the same time point. Supervised learning was highly effective at determining which participant had taken a given shoe or phone sample, and a Bayesian method was able to determine which participant had taken each shoe sample based entirely on its similarity to the floor samples. Both shoe and phone samples taken by conference participants clustered into distinct groups based on location, though much more so when an unweighted distance metric was used, suggesting sharing of low-abundance microbial taxa between individuals inhabiting the same space. CONCLUSIONS: Correlations between microbial community sources and sinks allow for inference of the interactions between humans and their environment. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s40168-015-0082-9) contains supplementary material, which is available to authorized users

Diet rapidly and reproducibly alters the human gut microbiome

Long-term diet influences the structure and activity of the trillions of microorganisms residing in the human gut1–5, but it remains unclear how rapidly and reproducibly the human gut microbiome responds to short-term macronutrient change. Here, we show that the short-term consumption of diets composed entirely of animal or plant products alters microbial community structure and overwhelms inter-individual differences in microbial gene expression. The animal-based diet increased the abundance of bile-tolerant microorganisms (Alistipes, Bilophila, and Bacteroides) and decreased the levels of Firmicutes that metabolize dietary plant polysaccharides (Roseburia, Eubacterium rectale, and Ruminococcus bromii). Microbial activity mirrored differences between herbivorous and carnivorous mammals2, reflecting trade-offs between carbohydrate and protein fermentation. Foodborne microbes from both diets transiently colonized the gut, including bacteria, fungi, and even viruses. Finally, increases in the abundance and activity of Bilophila wadsworthia on the animal-based diet support a link between dietary fat, bile acids, and the outgrowth of microorganisms capable of triggering inflammatory bowel disease6. In concert, these results demonstrate that the gut microbiome can rapidly respond to altered diet, potentially facilitating the diversity of human dietary lifestyles

Human oral viruses are personal, persistent and gender-consistent.

Viruses are the most abundant members of the human oral microbiome, yet relatively little is known about their biodiversity in humans. To improve our understanding of the DNA viruses that inhabit the human oral cavity, we examined saliva from a cohort of eight unrelated subjects over a 60-day period. Each subject was examined at 11 time points to characterize longitudinal differences in human oral viruses. Our primary goals were to determine whether oral viruses were specific to individuals and whether viral genotypes persisted over time. We found a subset of homologous viral genotypes across all subjects and time points studied, suggesting that certain genotypes may be ubiquitous among healthy human subjects. We also found significant associations between viral genotypes and individual subjects, indicating that viruses are a highly personalized feature of the healthy human oral microbiome. Many of these oral viruses were not transient members of the oral ecosystem, as demonstrated by the persistence of certain viruses throughout the entire 60-day study period. As has previously been demonstrated for bacteria and fungi, membership in the oral viral community was significantly associated with the sex of each subject. Similar characteristics of personalized, sex-specific microflora could not be identified for oral bacterial communities based on 16S rRNA. Our findings that many viruses are stable and individual-specific members of the oral ecosystem suggest that viruses have an important role in the human oral ecosystem