Microbiome datasets are compositional: and this is not optional

Datasets collected by high-throughput sequencing (HTS) of 16S rRNA gene amplimers, metagenomes or metatranscriptomes are commonplace and being used to study human disease states, ecological differences between sites, and the built environment. There is increasing awareness that microbiome datasets generated by HTS are compositional because they have an arbitrary total imposed by the instrument. However, many investigators are either unaware of this or assume specific properties of the compositional data. The purpose of this review is to alert investigators to the dangers inherent in ignoring the compositional nature of the data, and point out that HTS datasets derived from microbiome studies can and should be treated as compositions at all stages of analysis. We briefly introduce compositional data, illustrate the pathologies that occur when compositional data are analyzed inappropriately, and finally give guidance and point to resources and examples for the analysis of microbiome datasets using compositional data analysis.Peer ReviewedPostprint (published version

Unifying the analysis of high-throughput sequencing datasets: characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis.

BACKGROUND: Experimental designs that take advantage of high-throughput sequencing to generate datasets include RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), sequencing of 16S rRNA gene fragments, metagenomic analysis and selective growth experiments. In each case the underlying data are similar and are composed of counts of sequencing reads mapped to a large number of features in each sample. Despite this underlying similarity, the data analysis methods used for these experimental designs are all different, and do not translate across experiments. Alternative methods have been developed in the physical and geological sciences that treat similar data as compositions. Compositional data analysis methods transform the data to relative abundances with the result that the analyses are more robust and reproducible. RESULTS: Data from an in vitro selective growth experiment, an RNA-seq experiment and the Human Microbiome Project 16S rRNA gene abundance dataset were examined by ALDEx2, a compositional data analysis tool that uses Bayesian methods to infer technical and statistical error. The ALDEx2 approach is shown to be suitable for all three types of data: it correctly identifies both the direction and differential abundance of features in the differential growth experiment, it identifies a substantially similar set of differentially expressed genes in the RNA-seq dataset as the leading tools and it identifies as differential the taxa that distinguish the tongue dorsum and buccal mucosa in the Human Microbiome Project dataset. The design of ALDEx2 reduces the number of false positive identifications that result from datasets composed of many features in few samples. CONCLUSION: Statistical analysis of high-throughput sequencing datasets composed of per feature counts showed that the ALDEx2 R package is a simple and robust tool, which can be applied to RNA-seq, 16S rRNA gene sequencing and differential growth datasets, and by extension to other techniques that use a similar approach

A Multi-Platform Metabolomics Approach Identifies Highly Specific Biomarkers of Bacterial Diversity in the Vagina of Pregnant and Non-Pregnant Women

Bacterial vaginosis (BV) increases transmission of HIV, enhances the risk of preterm labour, and is associated with malodour. Clinical diagnosis often relies on microscopy, which may not reflect the microbiota composition accurately. We use an untargeted metabolomics approach, whereby we normalize the weight of samples prior to analysis, to obtained precise measurements of metabolites in vaginal fluid. We identify biomarkers for BV with high sensitivity and specificity (AUC = 0.99) in a cohort of 131 pregnant and non-pregnant Rwandan women, and demonstrate that the vaginal metabolome is strongly associated with bacterial diversity. Metabolites associated with high diversity and clinical BV include 2-hydroxyisovalerate and γ-hydroxybutyrate (GHB), but not succinate, which is produced by both Lactobacillus crispatus and BV-associated anaerobes in vitro. Biomarkers associated with high diversity and clinical BV are independent of pregnancy status, and were validated in a blinded replication cohort from Tanzania (n = 45), where we predicted clinical BV with 91% accuracy. Correlations between the metabolome and microbiota identified Gardnerella vaginalis as a putative producer of GHB, and we demonstrate production by this species in vitro. This work illustrates how changes in community structure alter the chemical composition of the vagina, and identifies highly specific biomarkers for a common condition

Vaginal Microbiome and Epithelial Gene Array in Post-Menopausal Women with Moderate to Severe Dryness

After menopause, many women experience vaginal dryness and atrophy of tissue, often attributed to the loss of estrogen. An understudied aspect of vaginal health in women who experience dryness due to atrophy is the role of the resident microbes. It is known that the microbiota has an important role in healthy vaginal homeostasis, including maintaining the pH balance and excluding pathogens. The objectives of this study were twofold: first to identify the microbiome of post-menopausal women with and without vaginal dryness and symptoms of atrophy; and secondly to examine any differences in epithelial gene expression associated with atrophy. The vaginal microbiome of 32 post-menopausal women was profiled using Illumina sequencing of the V6 region of the 16S rRNA gene. Sixteen subjects were selected for follow-up sampling every two weeks for 10 weeks. In addition, 10 epithelial RNA samples (6 healthy and 4 experiencing vaginal dryness) were acquired for gene expression analysis by Affymetrix Human Gene array. The microbiota abundance profiles were relatively stable over 10 weeks compared to previously published data on premenopausal women. There was an inverse correlation between Lactobacillus ratio and dryness and an increased bacterial diversity in women experiencing moderate to severe vaginal dryness. In healthy participants, Lactobacillus iners and L. crispatus were generally the most abundant, countering the long-held view that lactobacilli are absent or depleted in menopause. Vaginal dryness and atrophy were associated with down-regulation of human genes involved in maintenance of epithelial structure and barrier function, while those associated with inflammation were up-regulated consistent with the adverse clinical presentation

Microbiota of human breast tissue

In recent years, a greater appreciation for the microbes inhabiting human body sites has emerged. In the female mammary gland, milk has been shown to contain bacterial species, ostensibly reaching the ducts from the skin. We decided to investigate whether there is a microbiome within the mammary tissue. Using 16S rRNA sequencing and culture, we analyzed breast tissue from 81 women with and without cancer in Canada and Ireland. A diverse population of bacteria was detected within tissue collected from sites all around the breast in women aged 18 to 90, not all of whom had a history of lactation. The principal phylum was Proteobacteria. The most abundant taxa in the Canadian samples were Bacillus (11.4%), Acinetobacter (10.0%), Enterobacteriaceae (8.3%), Pseudomonas (6.5%), Staphylococcus (6.5%), Propionibacterium (5.8%), Comamonadaceae (5.7%), Gammaproteobacteria (5.0%), and Prevotella (5.0%). In the Irish samples the most abundant taxa were Enterobacteriaceae (30.8%), Staphylococcus (12.7%), Listeria welshimeri (12.1%), Propionibacterium (10.1%), and Pseudomonas (5.3%). None of the subjects had signs or symptoms of infection, but the presence of viable bacteria was confirmed in some samples by culture. The extent to which these organisms play a role in health or disease remains to be determined. © 2014, American Society for Microbiology

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

Deep Sequencing of the Vaginal Microbiota of Women with HIV

BACKGROUND: Women living with HIV and co-infected with bacterial vaginosis (BV) are at higher risk for transmitting HIV to a partner or newborn. It is poorly understood which bacterial communities constitute BV or the normal vaginal microbiota among this population and how the microbiota associated with BV responds to antibiotic treatment. METHODS AND FINDINGS: The vaginal microbiota of 132 HIV positive Tanzanian women, including 39 who received metronidazole treatment for BV, were profiled using Illumina to sequence the V6 region of the 16S rRNA gene. Of note, Gardnerella vaginalis and Lactobacillus iners were detected in each sample constituting core members of the vaginal microbiota. Eight major clusters were detected with relatively uniform microbiota compositions. Two clusters dominated by L. iners or L. crispatus were strongly associated with a normal microbiota. The L. crispatus dominated microbiota were associated with low pH, but when L. crispatus was not present, a large fraction of L. iners was required to predict a low pH. Four clusters were strongly associated with BV, and were dominated by Prevotella bivia, Lachnospiraceae, or a mixture of different species. Metronidazole treatment reduced the microbial diversity and perturbed the BV-associated microbiota, but rarely resulted in the establishment of a lactobacilli-dominated microbiota. CONCLUSIONS: Illumina based microbial profiling enabled high though-put analyses of microbial samples at a high phylogenetic resolution. The vaginal microbiota among women living with HIV in Sub-Saharan Africa constitutes several profiles associated with a normal microbiota or BV. Recurrence of BV frequently constitutes a different BV-associated profile than before antibiotic treatment

Adhesion of lactobacilli to polymer surfaces in vivo and in vitro

Lactobacillus iners is prominent in the human vagina and is able to persist despite development of bacterial vaginosis and treatment with antibiotics. A probable factor in its persistent survival is its ability to be retained in the vaginal epithelia. Genome sequencing of the strain showed an organism deplete of many metabolic pathways, yet equipped with fibronectin (Fn)-binding adhesins. The objective of the present study was to assess the ability of L iners AB-1 to bind immobilized Fn. Results showed that the organism superiorly bound the protein compared to other species of Lactobacillus and known binders such as Staphylococcus aureus. Treatment of L iners cells by protease rendered its binding abilities to Fn nonfunctional. The findings indicate a mechanism of vaginal persistence for a Lactobacillus species, with implications for reproductive health. © The Author(s) 2012