Nihar U. Sheth

The truth about metagenomics: quantifying and counteracting bias in 16S rRNA studies

BackgroundCharacterizing microbial communities via next-generation sequencing is subject to a number of pitfalls involving sample processing. The observed community composition can be a severe distortion of the quantities of bacteria actually present in the microbiome, hampering analysis and threatening the validity of conclusions from metagenomic studies. We introduce an experimental protocol using mock communities for quantifying and characterizing bias introduced in the sample processing pipeline. We used 80 bacterial mock communities comprised of prescribed proportions of cells from seven vaginally-relevant bacterial strains to assess the bias introduced in the sample processing pipeline. We created two additional sets of 80 mock communities by mixing prescribed quantities of DNA and PCR product to quantify the relative contribution to bias of (1) DNA extraction, (2) PCR amplification, and (3) sequencing and taxonomic classification for particular choices of protocols for each step. We developed models to predict the “true” composition of environmental samples based on the observed proportions, and applied them to a set of clinical vaginal samples from a single subject during four visits.ResultsWe observed that using different DNA extraction kits can produce dramatically different results but bias is introduced regardless of the choice of kit. We observed error rates from bias of over 85% in some samples, while technical variation was very low at less than 5% for most bacteria. The effects of DNA extraction and PCR amplification for our protocols were much larger than those due to sequencing and classification. The processing steps affected different bacteria in different ways, resulting in amplified and suppressed observed proportions of a community. When predictive models were applied to clinical samples from a subject, the predicted microbiome profiles were better reflections of the physiology and diagnosis of the subject at the visits than the observed community compositions.ConclusionsBias in 16S studies due to DNA extraction and PCR amplification will continue to require attention despite further advances in sequencing technology. Analysis of mock communities can help assess bias and facilitate the interpretation of results from environmental samples.

Differences in vaginal microbiome in African American women versus women of European ancestry.

Women of European ancestry are more likely to harbour a Lactobacillus-dominated microbiome, whereas African American women are more likely to exhibit a diverse microbial profile. African American women are also twice as likely to be diagnosed with bacterial vaginosis and are twice as likely to experience preterm birth. The objective of this study was to further characterize and contrast the vaginal microbial profiles in African American versus European ancestry women. Through the Vaginal Human Microbiome Project at Virginia Commonwealth University, 16S rRNA gene sequence analysis was used to compare the microbiomes of vaginal samples from 1268 African American women and 416 women of European ancestry. The results confirmed significant differences in the vaginal microbiomes of the two groups and identified several taxa relevant to these differences. Major community types were dominated by Gardnerella vaginalis and the uncultivated bacterial vaginosis-associated bacterium-1 (BVAB1) that were common among African Americans. Moreover, the prevalence of multiple bacterial taxa that are associated with microbial invasion of the amniotic cavity and preterm birth, including Mycoplasma, Gardnerella, Prevotella and Sneathia, differed between the two ethnic groups. We investigated the contributions of intrinsic and extrinsic factors, including pregnancy, body mass index, diet, smoking and alcohol use, number of sexual partners, and household income, to vaginal community composition. Ethnicity, pregnancy and alcohol use correlated significantly with the relative abundance of bacterial vaginosis-associated species. Trends between microbial profiles and smoking and number of sexual partners were observed; however, these associations were not statistically significant. These results support and extend previous findings that there are significant differences in the vaginal microbiome related to ethnicity and demonstrate that these differences are pronounced even in healthy women.

Species-level classification of the vaginal microbiome.

BackgroundThe application of next-generation sequencing to the study of the vaginal microbiome is revealing the spectrum of microbial communities that inhabit the human vagina. High-resolution identification of bacterial taxa, minimally to the species level, is necessary to fully understand the association of the vaginal microbiome with bacterial vaginosis, sexually transmitted infections, pregnancy complications, menopause, and other physiological and infectious conditions. However, most current taxonomic assignment strategies based on metagenomic 16S rDNA sequence analysis provide at best a genus-level resolution. While surveys of 16S rRNA gene sequences are common in microbiome studies, few well-curated, body-site-specific reference databases of 16S rRNA gene sequences are available, and no such resource is available for vaginal microbiome studies.ResultsWe constructed the Vaginal 16S rDNA Reference Database, a comprehensive and non-redundant database of 16S rDNA reference sequences for bacterial taxa likely to be associated with vaginal health, and we developed STIRRUPS, a new method that employs the USEARCH algorithm with a curated reference database for rapid species-level classification of 16S rDNA partial sequences. The method was applied to two datasets of V1-V3 16S rDNA reads: one generated from a mock community containing DNA from six bacterial strains associated with vaginal health, and a second generated from over 1,000 mid-vaginal samples collected as part of the Vaginal Human Microbiome Project at Virginia Commonwealth University. In both datasets, STIRRUPS, used in conjunction with the Vaginal 16S rDNA Reference Database, classified more than 95% of processed reads to a species-level taxon using a 97% global identity threshold for assignment.ConclusionsThis database and method provide accurate species-level classifications of metagenomic 16S rDNA sequence reads that will be useful for analysis and comparison of microbiome profiles from vaginal samples. STIRRUPS can be used to classify 16S rDNA sequence reads from other ecological niches if an appropriate reference database of 16S rDNA sequences is available.

Optimizing Read Mapping to Reference Genomes to Determine Composition and Species Prevalence in Microbial Communities

The Human Microbiome Project (HMP) aims to characterize the microbial communities of 18 body sites from healthy individuals. To accomplish this, the HMP generated two types of shotgun data: reference shotgun sequences isolated from different anatomical sites on the human body and shotgun metagenomic sequences from the microbial communities of each site. The alignment strategy for characterizing these metagenomic communities using available reference sequence is important to the success of HMP data analysis. Six next-generation aligners were used to align a community of known composition against a database comprising reference organisms known to be present in that community. All aligners report nearly complete genome coverage (>97%) for strains with over 6X depth of coverage, however they differ in speed, memory requirement and ease of use issues such as database size limitations and supported mapping strategies. The selected aligner was tested across a range of parameters to maximize sensitivity while maintaining a low false positive rate. We found that constraining alignment length had more impact on sensitivity than does constraining similarity in all cases tested. However, when reference species were replaced with phylogenetic neighbors, similarity begins to play a larger role in detection. We also show that choosing the top hit randomly when multiple, equally strong mappings are available increases overall sensitivity at the expense of taxonomic resolution. The results of this study identified a strategy that was used to map over 3 tera-bases of microbial sequence against a database of more than 5,000 reference genomes in just over a month.

An Emerging Mycoplasma Associated with Trichomoniasis, Vaginal Infection and Disease

Humans are colonized by thousands of bacterial species, but it is difficult to assess the metabolic and pathogenic potential of the majority of these because they have yet to be cultured. Here, we characterize an uncultivated vaginal mycoplasma tightly associated with trichomoniasis that was previously known by its 16S rRNA sequence as “Mnola.” In this study, the mycoplasma was found almost exclusively in women infected with the sexually transmitted pathogen Trichomonas vaginalis, but rarely observed in women with no diagnosed disease. The genomes of four strains of this species were reconstructed using metagenome sequencing and assembly of DNA from four discrete mid-vaginal samples, one of which was obtained from a pregnant woman with trichomoniasis who delivered prematurely. These bacteria harbor several putative virulence factors and display unique metabolic strategies. Genes encoding proteins with high similarity to potential virulence factors include two collagenases, a hemolysin, an O-sialoglycoprotein endopeptidase and a feoB-type ferrous iron transport system. We propose the name “Candidatus Mycoplasma girerdii” for this potential new pathogen.

Whole exome sequencing to estimate alloreactivity potential between donors and recipients in stem cell transplantation

Whole exome sequencing (WES) was performed on stem cell transplant donor‐recipient (D‐R) pairs to determine the extent of potential antigenic variation at a molecular level. In a small cohort of D‐R pairs, a high frequency of sequence variation was observed between the donor and recipient exomes independent of human leucocyte antigen (HLA) matching. Nonsynonymous, nonconservative single nucleotide polymorphisms were approximately twice as frequent in HLA‐matched unrelated, compared with related D‐R pairs. When mapped to individual chromosomes, these polymorphic nucleotides were uniformly distributed across the entire exome. In conclusion, WES reveals extensive nucleotide sequence variation in the exomes of HLA‐matched donors and recipients.

In silico Derivation of HLA-Specific Alloreactivity Potential from Whole Exome Sequencing of Stem-Cell Transplant Donors and Recipients: Understanding the Quantitative Immunobiology of Allogeneic Transplantation

Donor T-cell mediated graft versus host (GVH) effects may result from the aggregate alloreactivity to minor histocompatibility antigens (mHA) presented by the human leukocyte antigen (HLA) molecules in each donor–recipient pair undergoing stem-cell transplantation (SCT). Whole exome sequencing has previously demonstrated a large number of non-synonymous single nucleotide polymorphisms (SNP) present in HLA-matched recipients of SCT donors (GVH direction). The nucleotide sequence flanking each of these SNPs was obtained and the amino acid sequence determined. All the possible nonameric peptides incorporating the variant amino acid resulting from these SNPs were interrogated in silico for their likelihood to be presented by the HLA class I molecules using the Immune Epitope Database stabilized matrix method (SMM) and NetMHCpan algorithms. The SMM algorithm predicted that a median of 18,396 peptides weakly bound HLA class I molecules in individual SCT recipients, and 2,254 peptides displayed strong binding. A similar library of presented peptides was identified when the data were interrogated using the NetMHCpan algorithm. The bioinformatic algorithm presented here demonstrates that there may be a high level of mHA variation in HLA-matched individuals, constituting a HLA-specific alloreactivity potential.

Effects of combined oral contraceptives, depot medroxyprogesterone acetate and the levonorgestrel-releasing intrauterine system on the vaginal microbiome

OBJECTIVES Prior studies suggest that the composition of the vaginal microbiome may positively or negatively affect susceptibility to sexually transmitted infections (STIs) and bacterial vaginosis (BV). Some female hormonal contraceptive methods also appear to positively or negatively influence STI transmission and BV. Therefore, changes in the vaginal microbiome that are associated with different contraceptive methods may explain, in part, effects on STI transmission and BV. STUDY DESIGN We performed a retrospective study of 16S rRNA gene survey data of vaginal samples from a subset of participants from the Human Vaginal Microbiome Project at Virginia Commonwealth University. The subset included 682 women who reported using a single form of birth control that was condoms, combined oral contraceptives (COCs), depot medroxyprogesterone acetate (DMPA) or the levonorgestrel-releasing intrauterine system (LNG-IUS). RESULTS Women using COCs [adjusted odds ratio (aOR) 0.29, 95% confidence interval (CI) 0.13-0.64] and DMPA (aOR 0.34, 95% CI 0.13-0.89), but not LNG-IUS (aOR 1.55, 95% CI 0.72-3.35), were less likely to be colonized by BV-associated bacteria relative to women who used condoms. Women using COCs (aOR 1.94, 95% CI 1.25-3.02) were more likely to be colonized by beneficial H2O2-producing Lactobacillus species compared with women using condoms, while women using DMPA (aOR 1.09, 95% CI 0.63-1.86) and LNG-IUS (aOR 0.74, 95% CI 0.48-1.15) were not. CONCLUSIONS Use of COCs is significantly associated with increased vaginal colonization by healthy lactobacilli and reduced BV-associated taxa. IMPLICATIONS COC use may positively influence gynecologic health through an increase in healthy lactobacilli and a decrease in BV-associated bacterial taxa.

Skin-to-Skin Care and the Development of the Preterm Infant Oral Microbiome

OBJECTIVE The oral cavity represents an initial entry way for oral and gut indigenous colonization. Skin-to-skin (STS) care, in which the mother holds the diaper clad naked preterm (PT) infant between her breasts, is associated with improved digestive function, decreased stress, and improved survival. This study evaluated the development of oral microbial colonization repertoires and health characteristics in PT infants with or without STS exposure. METHODS Saliva from 42 PT infants (<32 weeks of gestation at birth) was collected prospectively at 1 month and/or at discharge. High-throughput 16S rRNA sequencing identified microbial diversity and prevalence of bacterial signatures correlated with clinical STS or non-STS care. RESULTS Corrected for gestational age (CGA) at sampling, bacterial taxa demonstrated increased Streptococcus as a signature of oral repertoire maturation. STS was associated with increased Streptococcus (p < 0.024), while non-STS was associated with greater Corynebacterium (p < 0.023) and Pseudomonas (p < 0.019) in infants ≤ 32 weeks CGA. In infants > 32 weeks CGA, Neisseria and Acinetobacter were more prevalent, 50 vs. 16.7% and 40 vs. 0%, respectively. STS care was associated with shorter hospitalization (p < 0.039). CONCLUSION STS care during earlier gestation was associated with a distinct microbial pattern and an accelerated pace of oral microbial repertoire maturity.

Stem cell transplantation as a dynamical system: are clinical outcomes deterministic?

Outcomes in stem cell transplantation (SCT) are modeled using probability theory. However, the clinical course following SCT appears to demonstrate many characteristics of dynamical systems, especially when outcomes are considered in the context of immune reconstitution. Dynamical systems tend to evolve over time according to mathematically determined rules. Characteristically, the future states of the system are predicated on the states preceding them, and there is sensitivity to initial conditions. In SCT, the interaction between donor T cells and the recipient may be considered as such a system in which, graft source, conditioning, and early immunosuppression profoundly influence immune reconstitution over time. This eventually determines clinical outcomes, either the emergence of tolerance or the development of graft versus host disease. In this paper, parallels between SCT and dynamical systems are explored and a conceptual framework for developing mathematical models to understand disparate transplant outcomes is proposed.