Characterizing the human vaginal microbiome using high-throughput sequencing

The human vaginal microbiome undoubtedly has a significant role in reproductive health and for protection from infectious organisms. Recent efforts to characterize the bacterial species of the vagina using molecular techniques have uncovered an unexpected diversity. Using high-throughput sequencing I sought to describe the structure and function of the vaginal microbiome under different physiological states including healthy, bacterial vaginosis (BV), post-menopausal vaginal atrophy, and acute vulvovaginal candidiasis (VVC). Partial 16S rRNA gene sequencing revealed that healthy, asymptomatic women most often have vaginal biotas dominated by Lactobacillus iners or L. crispatus. In contrast, BV is a heterogeneous, highly diversified condition with reduced Lactobacillus abundance. Similar to BV, post-menopausal women experiencing vaginal dryness were depleted in lactobacilli and had a more diverse vaginal profile. In the case of VVC, the biotas were not significantly altered compared to healthy women despite the fungal overgrowth. One organism, Lactobacillus iners was ubiquitously present in all conditions, and became predominant following antibiotic and probiotic treatment of BV. To uncover the potential role of this bacterium, I used whole genome sequencing of vaginal isolate AB-1. The genome is predicted to be the smallest of any Lactobacillus at 1.3 Mbp, but having a higher proportion of horizontally acquired genes. These results, along with predicted adhesins and a cholesterol-dependent cytolysin, indicate L. iners is highly adapted for the vagina and could have an uncharacterized role in the etiology of BV. As BV is the most common vaginal ailment with severe implications on acquisition and transmission of sexually transmitted infections, and complications during pregnancy, I sought iito examine the functional contribution of the organisms during BV using meta-RNA sequencing. L. iners drastically modulates gene expression in response to BV, and notably increases expression of a cholesterol-dependent cytolysin, mucin and glycerol transport and metabolic enzymes, and genes belonging to a CRISPR system - suggestive of bacteriophage influence in the community. Although diverse in taxonomic membership, there is evidence of functional conservation in BV including preference for glycogen and glycerol as carbon sources, and predicted end products of metabolism including an abundance of succinate and short-chain fatty acids. These studies add significantly to our understanding of the role lactobacilli can play in vaginal and reproductive health

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

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

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

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

Spontaneous Preterm Birth Is Associated with Differential Expression of Vaginal Metabolites by Lactobacilli-Dominated Microflora

A major challenge in preventing preterm birth (PTB) is identifying women at greatest risk. This pilot study prospectively examined the differences in vaginal microbiota and metabolite profiles of women who delivered prematurely compared to their term counterparts in a cohort of asymptomatic (studied at 20–22, n = 80; and 26–28 weeks, n = 41) and symptomatic women (studied at 24–36 weeks, n = 37). Using 16S rRNA sequencing, the vaginal microbiota from cervicovaginal fluid samples was characterized into five Community State Types (CST) dominated by Lactobacillus spp.: CSTI (Lactobacillus crispatus), CSTII (Lactobacillus gasseri), CSTIII (Lactobacillus iners), CSTV (Lactobacillus jensenii); and mixed anaerobes—CSTIV. This was then related to the vaginal metabolite profile and pH determined by 1H-Nuclear Magnetic Resonance spectroscopy and pH indicator paper, respectively. At 20–22 weeks, the term-delivered women (TDW) indicated a proportion of CSTI-dominated microbiota >2-fold higher compared to the preterm-delivered women (PTDW) (40.3 vs. 16.7%, P = 0.0002), and a slightly higher proportion at 26–28 weeks (20.7 vs. 16.7%, P = 0.03). CSTV was >2-fold higher in the PTDW compared to TDW at 20–22 (22.2 vs. 9.7%, P = 0.0002) and 26–28 weeks (25.0 vs. 10.3%, P = 0.03). Furthermore, at 26–28 weeks no PTDW had a CSTII-dominated microbiome, in contrast to 28% of TDW (P < 0.0001). CSTI-dominated samples showed higher lactate levels than CSTV at 20–22 weeks (P < 0.01), and 26–28 weeks (P < 0.05), while CSTII-dominated samples indicated raised succinate levels over CSTV at 26–28 weeks (P < 0.05). These were supported by Principal coordinates analysis, which revealed strong clustering of metabolites according to CST. In addition, the CSTI-dominated samples had an average pH of 3.8, which was lower than those of CSTII—4.4, and CSTV—4.2 (P < 0.05). Elevated vaginal lactate and succinate were associated with predominance of CSTI and II over CSTV in women who delivered at term compared with their preterm counterparts. This suggests that L. jensenii-dominance and decreased lactate and/or succinate could increase the risk of PTB, while L. crispatus/gasseri may confer some protection against inflammation-associated PTB and highlight the need for further study in this area

Adhesion Forces and Coaggregation between Vaginal Staphylococci and Lactobacilli

Urogenital infections are the most common ailments afflicting women. They are treated with dated antimicrobials whose efficacy is diminishing. The process of infection involves pathogen adhesion and displacement of indigenous Lactobacillus crispatus and Lactobacillus jensenii. An alternative therapeutic approach to antimicrobial therapy is to reestablish lactobacilli in this microbiome through probiotic administration. We hypothesized that lactobacilli displaying strong adhesion forces with pathogens would facilitate coaggregation between the two strains, ultimately explaining the elimination of pathogens seen in vivo. Using atomic force microscopy, we found that adhesion forces between lactobacilli and three virulent toxic shock syndrome toxin 1-producing Staphylococcus aureus strains, were significantly stronger (2.2–6.4 nN) than between staphylococcal pairs (2.2–3.4 nN), especially for the probiotic Lactobacillus reuteri RC-14 (4.0–6.4 nN) after 120 s of bond-strengthening. Moreover, stronger adhesion forces resulted in significantly larger coaggregates. Adhesion between the bacteria occurred instantly upon contact and matured within one to two minutes, demonstrating the potential for rapid anti-pathogen effects using a probiotic. Coaggregation is one of the recognized mechanisms through which lactobacilli can exert their probiotic effects to create a hostile micro-environment around a pathogen. With antimicrobial options fading, it therewith becomes increasingly important to identify lactobacilli that bind strongly with pathogens