Jean M. Macklaim

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 allowed an in-depth analysis of errors, and we found that PCR-induced errors composed the vast majority of non-organism derived species variants, an observation 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 bacterial 16S rRNA genes and the eukaryotic V9 region that is gaining popularity for sampling protist diversity.

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.

Comparative meta-RNA-seq of the vaginal microbiota and differential expression by Lactobacillus iners in health and dysbiosis

BackgroundBacterial vaginosis (BV), the most common vaginal condition of reproductive-aged women, is associated with a highly diverse and heterogeneous microbiota. Here we present a proof-of-principle analysis to uncover the function of the microbiota using meta-RNA-seq to uncover genes and pathways that potentially differentiate healthy vaginal microbial communities from those in the dysbiotic state of bacterial vaginosis (BV).ResultsThe predominant organism, Lactobacillus iners, was present in both conditions and showed a differing expression profile in BV compared to healthy. Despite its minimal genome, L. iners differentially expressed over 10% of its gene complement. Notably, in a BV environment L. iners increased expression of a cholesterol-dependent cytolysin, and of mucin and glycerol transport and related metabolic enzymes. Genes belonging to a CRISPR system were greatly upregulated suggesting that bacteriophage influence the community. Reflective of L. iners, the bacterial community as a whole demonstrated a preference for glycogen and glycerol as carbon sources under BV conditions. The predicted end-products of metabolism under BV conditions include an abundance of succinate and other short-chain fatty-acids, while healthy conditions are predicted to largely contain lactic acid.ConclusionsOur study underscores the importance of understanding the functional activity of the bacterial community in addition to characterizing the population structure when investigating the human microbiome.

At the crossroads of vaginal health and disease, the genome sequence of Lactobacillus iners AB-1

Lactobacilli have long been regarded as important constituents of the healthy human vagina. Lactobacillus iners is the most frequently detected bacterial species in the vagina, but little is known about its characteristics. We report a description of the whole-genome sequence of L. iners AB-1 along with comparative analysis of published genomes of closely related strains of lactobacilli. The genome is the smallest Lactobacillus reported to date, with a 1.3-Mbp single chromosome. The genome seems to have undergone one or more rapid evolution events that resulted in large-scale gene loss and horizontal acquisition of a number of genes for survival in the vagina. L. iners may exhibit specialized adaptation mechanisms to the vaginal environment, such as an iron–sulfur cluster assembly system, and several unique σ factors to regulate gene transcription in this fluctuating environment. A potentially highly expressed homolog of a cholesterol-binding lysin may also contribute to host cell adhesion or act as a defense mechanism against other microbes. Notably, there is a lack of apparent adhesion proteins, but several cell-anchor proteins were identified and may be important for interaction with the host mucosal tissues. L. iners is widely present in healthy females as well as those suffering from bacterial vaginosis or who have undergone antimicrobial therapy, suggesting that it is an important indigenous species of the vagina.

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

BackgroundExperimental 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.ResultsData 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.ConclusionStatistical 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.

Microbiota of Human Breast Tissue

ABSTRACT 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.

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.

ANOVA-like differential expression (ALDEx) analysis for mixed population RNA-Seq.

Experimental variance is a major challenge when dealing with high-throughput sequencing data. This variance has several sources: sampling replication, technical replication, variability within biological conditions, and variability between biological conditions. The high per-sample cost of RNA-Seq often precludes the large number of experiments needed to partition observed variance into these categories as per standard ANOVA models. We show that the partitioning of within-condition to between-condition variation cannot reasonably be ignored, whether in single-organism RNA-Seq or in Meta-RNA-Seq experiments, and further find that commonly-used RNA-Seq analysis tools, as described in the literature, do not enforce the constraint that the sum of relative expression levels must be one, and thus report expression levels that are systematically distorted. These two factors lead to misleading inferences if not properly accommodated. As it is usually only the biological between-condition and within-condition differences that are of interest, we developed ALDEx, an ANOVA-like differential expression procedure, to identify genes with greater between- to within-condition differences. We show that the presence of differential expression and the magnitude of these comparative differences can be reasonably estimated with even very small sample sizes.

Changes in vaginal microbiota following antimicrobial and probiotic therapy

Background The composition of the vaginal microbiota is known to be important for health. When infections occur, antimicrobial therapy is often poorly efficacious. Objective and design We used 16S rRNA gene sequencing to characterize changes in the bacterial microbiota following oral antimicrobial and probiotic interventions. Results While the bacterial vaginal profiles of women with vulvovaginal candidiasis were dominated by lactobacilli as in healthy women, and unchanged by therapy, Gardnerella vaginalis, Prevotella, Atopobium, Sneathia, and Megasphaera dominated the vagina of women with bacterial vaginosis (BV), and treatment with tinidazole plus Lactobacillus reuteri RC-14+L. rhamnosus GR-1 resulted in an increased relative abundance of indigenous L. iners or L. crispatus. Conclusions The ability to restore homeostasis provides a rationale for conjoint use of probiotics with antibiotic treatment of BV.

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.