MeFiT: merging and filtering tool for illumina paired-end reads for 16S rRNA amplicon sequencing.

BACKGROUND: Recent advances in next-generation sequencing have revolutionized genomic research. 16S rRNA amplicon sequencing using paired-end sequencing on the MiSeq platform from Illumina, Inc., is being used to characterize the composition and dynamics of extremely complex/diverse microbial communities. For this analysis on the Illumina platform, merging and quality filtering of paired-end reads are essential first steps in data analysis to ensure the accuracy and reliability of downstream analysis. RESULTS: We have developed the Merging and Filtering Tool (MeFiT) to combine these pre-processing steps into one simple, intuitive pipeline. MeFiT invokes CASPER (context-aware scheme for paired-end reads) for merging paired-end reads and provides users the option to quality filter the reads using the traditional average Q-score metric or using a maximum expected error cut-off threshold. CONCLUSIONS: MeFiT provides an open-source solution that permits users to merge and filter paired end illumina reads. The tool has been implemented in python and the source-code is freely available at https://github.com/nisheth/MeFiT .

Comparison of Lactobacillus crispatus isolates from Lactobacillus-dominated vaginal microbiomes with isolates from microbiomes containing bacterial vaginosis-associated bacteria.

Vaginal lactobacilli can inhibit colonization by and growth of other bacteria, thereby preventing development of bacterial vaginosis (BV). Amongst the lactobacilli, Lactobacillus crispatus appears to be particularly effective at inhibiting growth of BV-associated bacteria. Nonetheless, some women who are colonized with this species can still develop clinical BV. Therefore, we sought to determine whether strains of L. crispatus that colonize women with lactobacilli-dominated vaginal microbiomes are distinct from strains that colonize women who develop BV. The genomes of L. crispatus isolates from four women with lactobacilli-dominated vaginal microbiomes ( <1% 16S rRNA reads above threshold from genera other than Lactobacillus) and four women with microbiomes containing BV-associated bacteria (>12% 16S rRNA reads from bacterial taxa associated with BV) were sequenced and compared. Lactic acid production by the different strains was quantified. Phage induction in the strains was also analysed. There was considerable genetic diversity between strains, and several genes were exclusive to either the strains from Lactobacillus-dominated microbiomes or those containing BV-associated bacteria. Overall, strains from microbiomes dominated by lactobacilli did not differ from strains from microbiomes containing BV-associated bacteria with respect to lactic acid production. All of the strains contained multiple phage, but there was no clear distinction between the presence or absence of BV-associated bacteria with respect to phage-induced lysis. Genes found to be exclusive to the Lactobacillus-dominated versus BV-associated bacteria-containing microbiomes could play a role in the maintenance of vaginal health and the development of BV, respectively.

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

BACKGROUND: Characterizing 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. RESULTS: We 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. CONCLUSIONS: Bias 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.

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.

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.

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.

BACKGROUND: The 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. RESULTS: We 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. CONCLUSIONS: This 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.

Racioethnic diversity in the dynamics of the vaginal microbiome during pregnancy.

The microbiome of the female reproductive tract has implications for women's reproductive health. We examined the vaginal microbiome in two cohorts of women who experienced normal term births: a cross-sectionally sampled cohort of 613 pregnant and 1,969 non-pregnant women, focusing on 300 pregnant and 300 non-pregnant women of African, Hispanic or European ancestry case-matched for race, gestational age and household income; and a longitudinally sampled cohort of 90 pregnant women of African or non-African ancestry. In these women, the vaginal microbiome shifted during pregnancy toward Lactobacillus-dominated profiles at the expense of taxa often associated with vaginal dysbiosis. The shifts occurred early in pregnancy, followed predictable patterns, were associated with simplification of the metabolic capacity of the microbiome and were significant only in women of African or Hispanic ancestry. Both genomic and environmental factors are likely contributors to these trends, with socioeconomic status as a likely environmental influence.

The vaginal microbiome and preterm birth.

The incidence of preterm birth exceeds 10% worldwide. There are significant disparities in the frequency of preterm birth among populations within countries, and women of African ancestry disproportionately bear the burden of risk in the United States. In the present study, we report a community resource that includes 'omics' data from approximately 12,000 samples as part of the integrative Human Microbiome Project. Longitudinal analyses of 16S ribosomal RNA, metagenomic, metatranscriptomic and cytokine profiles from 45 preterm and 90 term birth controls identified harbingers of preterm birth in this cohort of women predominantly of African ancestry. Women who delivered preterm exhibited significantly lower vaginal levels of Lactobacillus crispatus and higher levels of BVAB1, Sneathia amnii, TM7-H1, a group of Prevotella species and nine additional taxa. The first representative genomes of BVAB1 and TM7-H1 are described. Preterm-birth-associated taxa were correlated with proinflammatory cytokines in vaginal fluid. These findings highlight new opportunities for assessment of the risk of preterm birth.

Ethanol-Induced Behavioral Sensitization Alters the Synaptic Transcriptome and Exon Utilization in DBA/2J Mice.

Alcoholism is a complex behavioral disorder characterized by loss of control in limiting intake, and progressive compulsion to seek and consume ethanol. Prior studies have suggested that the characteristic behaviors associated with escalation of drug use are caused, at least in part, by ethanol-evoked changes in gene expression affecting synaptic plasticity. Implicit in this hypothesis is a dependence on new protein synthesis and remodeling at the synapse. It is well established that mRNA can be transported to distal dendritic processes, where it can undergo localized translation. It is unknown whether such modulation of the synaptic transcriptome might contribute to ethanol-induced synaptic plasticity. Using ethanol-induced behavioral sensitization as a model of neuroplasticity, we investigated whether repeated exposure to ethanol altered the synaptic transcriptome, contributing to mechanisms underlying subsequent increases in ethanol-evoked locomotor activity. RNAseq profiling of DBA/2J mice subjected to acute ethanol or ethanol-induced behavioral sensitization was performed on frontal pole synaptoneurosomes to enrich for synaptic mRNA. Genomic profiling showed distinct functional classes of mRNA enriched in the synaptic vs. cytosolic fractions, consistent with their role in synaptic function. Ethanol sensitization regulated more than twice the number of synaptic localized genes compared to acute ethanol exposure. Synaptic biological processes selectively perturbed by ethanol sensitization included protein folding and modification as well as and mitochondrial respiratory function, suggesting repeated ethanol exposure alters synaptic energy production and the processing of newly translated proteins. Additionally, marked differential exon usage followed ethanol sensitization in both synaptic and non-synaptic cellular fractions, with little to no perturbation following acute ethanol exposure. Altered synaptic exon usage following ethanol sensitization strongly affected genes related to RNA processing and stability, translational regulation, and synaptic function. These genes were also enriched for targets of the FMRP RNA-binding protein and contained consensus sequence motifs related to other known RNA binding proteins, suggesting that ethanol sensitization altered selective mRNA trafficking mechanisms. This study provides a foundation for investigating the role of ethanol in modifying the synaptic transcriptome and inducing changes in synaptic plasticity.

Rare mutations and potentially damaging missense variants in genes encoding fibrillar collagens and proteins involved in their production are candidates for risk for preterm premature rupture of membranes.

Preterm premature rupture of membranes (PPROM) is the leading identifiable cause of preterm birth with ~ 40% of preterm births being associated with PPROM and occurs in 1% - 2% of all pregnancies. We hypothesized that multiple rare variants in fetal genes involved in extracellular matrix synthesis would associate with PPROM, based on the assumption that impaired elaboration of matrix proteins would reduce fetal membrane tensile strength, predisposing to unscheduled rupture. We performed whole exome sequencing (WES) on neonatal DNA derived from pregnancies complicated by PPROM (49 cases) and healthy term deliveries (20 controls) to identify candidate mutations/variants. Genotyping for selected variants from the WES study was carried out on an additional 188 PPROM cases and 175 controls. All mothers were self-reported African Americans, and a panel of ancestry informative markers was used to control for genetic ancestry in all genetic association tests. In support of the primary hypothesis, a statistically significant genetic burden (all samples combined, SKAT-O p-value = 0.0225) of damaging/potentially damaging rare variants was identified in the genes of interest-fibrillar collagen genes, which contribute to fetal membrane strength and integrity. These findings suggest that the fetal contribution to PPROM is polygenic, and driven by an increased burden of rare variants that may also contribute to the disparities in rates of preterm birth among African Americans.

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.

Association between statin use, the vaginal microbiome, and Gardnerella vaginalis vaginolysin-mediated cytotoxicity.

BACKGROUND: Bacterial vaginosis (BV) is the leading dysbiosis of the vaginal microbiome. The pathways leading towards the development of BV are not well understood. Gardnerella vaginalis is frequently associated with BV. G. vaginalis produces the cholesterol-dependent cytolysin (CDC), vaginolysin, which can lyse a variety of human cells and is thought to play a role in pathogenesis. Because membrane cholesterol is required for vaginolysin to function, and because HMG-CoA reductase inhibitors (statins) affect not only serum levels of cholesterol but membrane levels as well, we hypothesized that statins might affect the vaginal microbiome. METHODS: To investigate the relationship between use of the statins and the vaginal microbiome, we analyzed 16S rRNA gene taxonomic surveys performed on vaginal samples from 133 women who participated in the Vaginal Human Microbiome Project and who were taking statins at the time of sampling, 152 women who reported high cholesterol levels but were not taking statins, and 316 women who did not report high cholesterol. To examine the effect of statins on the cytolytic effect of vaginolysin, the cholesterol-dependent cytolysin (CDC) produced by Gardnerella vaginalis, we assessed the effect of simvastatin pretreatment of VK2E6/E7 vaginal epithelial cells on vaginolysin-mediated cytotoxicity. RESULTS: The mean proportion of G. vaginalis among women taking statins was significantly lower relative to women not using statins. Women using statins had higher mean proportions of Lactobacillus crispatus relative to women with normal cholesterol levels, and higher levels of Lactobacillus jensenii relative to women with high cholesterol but not taking statins. In vitro, vaginal epithelial cells pretreated with simvastatin were relatively resistant to vaginolysin and this effect was inhibited by cholesterol. CONCLUSIONS: In this cross-sectional study, statin use was associated with reduced proportions of G. vaginalis and greater proportions of beneficial lactobacilli within the vaginal microbiome. The negative association between statin use and G. vaginalis may be related to inhibition of vaginolysin function.

BOTUX: bayesian-like operational taxonomic unit examiner.

Bayesian-like operational taxonomic unit examiner (BOTUX) is a new tool for the classification of 16S rRNA gene sequences into operational taxonomic units (OTUs) that addresses the problem of overestimation caused by errors introduced during PCR amplification and DNA sequencing steps. BOTUX utilises a grammar-based assignment strategy, where Bayesian models are built from each word of a given length (e.g., 8-mers). de novo analysis is possible with BOTUX as it does not require a training set, and updates probabilistic models as new sequences are recruited to an OTU. In benchmarking tests performed with real and simulated datasets of 16S rDNA sequences, BOTUX accurately identifies OTUs with comparable or better clustering efficiency and lower execution times than other OTU algorithms tested. BOTUX is the only OTU classifier, which allows incremental analysis of large datasets, and is also adept in clustering both 454 and Illumina datasets in a reasonable timeframe.

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.

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.