Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments.

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

Group B Streptococcus and the Vaginal Microbiota.

Background: Streptococcus agalactiae (group B Streptococcus [GBS]) is an important neonatal pathogen and emerging cause of disease in adults. The major risk factor for neonatal disease is maternal vaginal colonization. However, little is known about the relationship between GBS and vaginal microbiota. Methods: Vaginal lavage samples from nonpregnant women were tested for GBS, and amplicon-based sequencing targeting the 16S ribosomal RNA V3-V4 region was performed. Results: Four hundred twenty-eight of 432 samples met the high-quality read threshold. There was no relationship between GBS carriage and demographic characteristics, α-diversity, or overall vaginal microbiota community state type (CST). Within the non-Lactobacillus-dominant CST IV, GBS positive status was significantly more prevalent in CST IV-A than CST IV-B. Significant clustering by GBS status was noted on principal coordinates analysis, and 18 individual taxa were found to be significantly associated with GBS carriage by linear discriminant analysis. After adjusting for race/ethnicity, 4 taxa were positively associated with GBS, and 6 were negatively associated. Conclusions: Vaginal microbiota CST and α-diversity are not related to GBS status. However, specific microbial taxa are associated with colonization of this important human pathogen, highlighting a potential role for the microbiota in promotion or inhibition of GBS colonization.

Insight into the ecology of vaginal bacteria through integrative analyses of metagenomic and metatranscriptomic data.

BACKGROUND: Vaginal bacterial communities dominated by Lactobacillus species are associated with a reduced risk of various adverse health outcomes. However, somewhat unexpectedly, many healthy women have microbiota that are not dominated by lactobacilli. To determine the factors that drive vaginal community composition we characterized the genetic composition and transcriptional activities of vaginal microbiota in healthy women. RESULTS: We demonstrate that the abundance of a species is not always indicative of its transcriptional activity and that impending changes in community composition can be predicted from metatranscriptomic data. Functional comparisons highlight differences in the metabolic activities of these communities, notably in their degradation of host produced mucin but not glycogen. Degradation of mucin by communities not dominated by Lactobacillus may play a role in their association with adverse health outcomes. Finally, we show that the transcriptional activities of L. crispatus, L. iners, and Gardnerella vaginalis vary with the taxonomic composition of the communities in which they reside. Notably, L. iners and G. vaginalis both demonstrate lower expression of their cholesterol-dependent cytolysins when co-resident with Lactobacillus spp. and higher expression when co-resident with other facultative and obligate anaerobes. The pathogenic potential of these species may depend on the communities in which they reside and thus could be modulated by interventional strategies. CONCLUSIONS: Our results provide insight to the functional ecology of the vaginal microbiota, demonstrate the diagnostic potential of metatranscriptomic data, and reveal strategies for the management of these ecosystems.

Dual RNA-seq analysis of in vitro infection multiplicity and RNA depletion methods in Chlamydia-infected epithelial cells.

Dual RNA-seq experiments examining viral and bacterial pathogens are increasing, but vary considerably in their experimental designs, such as infection rates and RNA depletion methods. Here, we have applied dual RNA-seq to Chlamydia trachomatis infected epithelial cells to examine transcriptomic responses from both organisms. We compared two time points post infection (1 and 24 h), three multiplicity of infection (MOI) ratios (0.1, 1 and 10) and two RNA depletion methods (rRNA and polyA). Capture of bacterial-specific RNA were greatest when combining rRNA and polyA depletion, and when using a higher MOI. However, under these conditions, host RNA capture was negatively impacted. Although it is tempting to use high infection rates, the implications on host cell survival, the potential reduced length of infection cycles and real world applicability should be considered. This data highlights the delicate nature of balancing host-pathogen RNA capture and will assist future transcriptomic-based studies to achieve more specific and relevant infection-related biological insights.

NsrR targets in the Escherichia coli genome: new insights into DNA sequence requirements for binding and a role for NsrR in the regulation of motility.

The Escherichia coli NsrR protein is a nitric oxide-sensitive repressor of transcription. The NsrR-binding site is predicted to comprise two copies of an 11 bp motif arranged as an inverted repeat with 1 bp spacing. By mutagenesis we confirmed that both 11 bp motifs are required for maximal NsrR repression of the ytfE promoter. We used chromatin immunoprecipitation and microarray analysis (ChIP-chip) to show that NsrR binds to 62 sites close to the 5' ends of genes. Analysis of the ChIP-chip data suggested that a single 11 bp motif (with the consensus sequence AANATGCATTT) can function as an NsrR-binding site in vivo. NsrR binds to sites in the promoter regions of the fliAZY, fliLMNOPQR and mqsR-ygiT transcription units, which encode proteins involved in motility and biofilm development. Reporter fusion assays confirmed that NsrR negatively regulates the fliA and fliL promoters. A mutation in the predicted 11 bp NsrR-binding site in the fliA promoter impaired repression by NsrR and prevented detectable binding in vivo. Assays on soft-agar confirmed that NsrR is a negative regulator of motility in E. coli K12 and in a uropathogenic strain; surface attachment assays revealed decreased levels of attached growth in the absence of NsrR.

Chromatin accessibility dynamics of Chlamydia-infected epithelial cells.

Chlamydia are Gram-negative, obligate intracellular bacterial pathogens responsible for a broad spectrum of human and animal diseases. In humans, Chlamydia trachomatis is the most prevalent bacterial sexually transmitted infection worldwide and is the causative agent of trachoma (infectious blindness) in disadvantaged populations. Over the course of its developmental cycle, Chlamydia extensively remodels its intracellular niche and parasitises the host cell for nutrients, with substantial resulting changes to the host cell transcriptome and proteome. However, little information is available on the impact of chlamydial infection on the host cell epigenome and global gene regulation. Regions of open eukaryotic chromatin correspond to nucleosome-depleted regions, which in turn are associated with regulatory functions and transcription factor binding. We applied formaldehyde-assisted isolation of regulatory elements enrichment followed by sequencing (FAIRE-Seq) to generate temporal chromatin maps of C. trachomatis-infected human epithelial cells in vitro over the chlamydial developmental cycle. We detected both conserved and distinct temporal changes to genome-wide chromatin accessibility associated with C. trachomatis infection. The observed differentially accessible chromatin regions include temporally-enriched sets of transcription factors, which may help shape the host cell response to infection. These regions and motifs were linked to genomic features and genes associated with immune responses, re-direction of host cell nutrients, intracellular signalling, cell-cell adhesion, extracellular matrix, metabolism and apoptosis. This work provides another perspective to the complex response to chlamydial infection, and will inform further studies of transcriptional regulation and the epigenome in Chlamydia-infected human cells and tissues.

A comprehensive non-redundant gene catalog reveals extensive within-community intraspecies diversity in the human vagina.

Analysis of metagenomic and metatranscriptomic data is complicated and typically requires extensive computational resources. Leveraging a curated reference database of genes encoded by members of the target microbiome can make these analyses more tractable. In this study, we assemble a comprehensive human vaginal non-redundant gene catalog (VIRGO) that includes 0.95 million non-redundant genes. The gene catalog is functionally and taxonomically annotated. We also construct a vaginal orthologous groups (VOG) from VIRGO. The gene-centric design of VIRGO and VOG provides an easily accessible tool to comprehensively characterize the structure and function of vaginal metagenome and metatranscriptome datasets. To highlight the utility of VIRGO, we analyze 1,507 additional vaginal metagenomes, and identify a high degree of intraspecies diversity within and across vaginal microbiota. VIRGO offers a convenient reference database and toolkit that will facilitate a more in-depth understanding of the role of vaginal microorganisms in women's health and reproductive outcomes.

VALENCIA: a nearest centroid classification method for vaginal microbial communities based on composition.

BACKGROUND: Taxonomic profiles of vaginal microbial communities can be sorted into a discrete number of categories termed community state types (CSTs). This approach is advantageous because collapsing a hyper-dimensional taxonomic profile into a single categorical variable enables efforts such as data exploration, epidemiological studies, and statistical modeling. Vaginal communities are typically assigned to CSTs based on the results of hierarchical clustering of the pairwise distances between samples. However, this approach is problematic because it complicates between-study comparisons and because the results are entirely dependent on the particular set of samples that were analyzed. We sought to standardize and advance the assignment of samples to CSTs. RESULTS: We developed VALENCIA (VAginaL community state typE Nearest CentroId clAssifier), a nearest centroid-based tool which classifies samples based on their similarity to a set of reference centroids. The references were defined using a comprehensive set of 13,160 taxonomic profiles from 1975 women in the USA. This large dataset allowed us to comprehensively identify, define, and characterize vaginal CSTs common to reproductive age women and expand upon the CSTs that had been defined in previous studies. We validated the broad applicability of VALENCIA for the classification of vaginal microbial communities by using it to classify three test datasets which included reproductive age eastern and southern African women, adolescent girls, and a racially/ethnically and geographically diverse sample of postmenopausal women. VALENCIA performed well on all three datasets despite the substantial variations in sequencing strategies and bioinformatics pipelines, indicating its broad application to vaginal microbiota. We further describe the relationships between community characteristics (vaginal pH, Nugent score) and participant demographics (race, age) and the CSTs defined by VALENCIA. CONCLUSION: VALENCIA provides a much-needed solution for the robust and reproducible assignment of vaginal community state types. This will allow unbiased analysis of both small and large vaginal microbiota datasets, comparisons between datasets and meta-analyses that combine multiple datasets. Video abstract.

Cervicovaginal microbiota and local immune response modulate the risk of spontaneous preterm delivery.

Failure to predict and understand the causes of preterm birth, the leading cause of neonatal morbidity and mortality, have limited effective interventions and therapeutics. From a cohort of 2000 pregnant women, we performed a nested case control study on 107 well-phenotyped cases of spontaneous preterm birth (sPTB) and 432 women delivering at term. Using innovative Bayesian modeling of cervicovaginal microbiota, seven bacterial taxa were significantly associated with increased risk of sPTB, with a stronger effect in African American women. However, higher vaginal levels of ß-defensin-2 lowered the risk of sPTB associated with cervicovaginal microbiota in an ethnicity-dependent manner. Surprisingly, even in Lactobacillus spp. dominated cervicovaginal microbiota, low ß-defensin-2 was associated with increased risk of sPTB. These findings hold promise for diagnostics to accurately identify women at risk for sPTB early in pregnancy. Therapeutic strategies could include immune modulators and microbiome-based therapeutics to reduce this significant health burden.

Dual RNA-Seq of Chlamydia and Host Cells.

During the infection of a host cell by a bacterial pathogen, a cascading series of gene expression changes occurs as each organism manipulates or responds to the other via defense or survival strategies. Unraveling this complex interplay is key for our understanding of bacterial virulence and host response pathways for the development of novel therapeutics. Dual RNA sequencing (dual RNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. Leveraging the sensitivity and resolution allowed by RNA-seq, dual RNA-Seq can be applied to any bacteria-eukaryotic host interaction. We pioneered dual RNA-Seq to simultaneously capture Chlamydia and host expression profiles during an in vitro infection as proof of principle. Here we provide a detailed laboratory protocol and bioinformatics analysis guidelines for dual RNA-seq experiments focusing on Chlamydia as the organism of interest.

Early Transcriptional Landscapes of Chlamydia trachomatis-Infected Epithelial Cells at Single Cell Resolution.

Chlamydia are Gram-negative obligate intracellular bacterial pathogens responsible for a variety of disease in humans and animals worldwide. Chlamydia trachomatis causes trachoma in disadvantaged populations, and is the most common bacterial sexually transmitted infection in humans, causing reproductive tract disease. Antibiotic therapy successfully treats diagnosed chlamydial infections, however asymptomatic infections are common. High-throughput transcriptomic approaches have explored chlamydial gene expression and infected host cell gene expression. However, these were performed on large cell populations, averaging gene expression profiles across all cells sampled and potentially obscuring biologically relevant subsets of cells. We generated a pilot dataset, applying single cell RNA-Seq (scRNA-Seq) to C. trachomatis infected and mock-infected epithelial cells to assess the utility, pitfalls and challenges of single cell approaches applied to chlamydial biology, and to potentially identify early host cell biomarkers of chlamydial infection. Two hundred sixty-four time-matched C. trachomatis-infected and mock-infected HEp-2 cells were collected and subjected to scRNA-Seq. After quality control, 200 cells were retained for analysis. Two distinct clusters distinguished 3-h cells from 6- and 12-h. Pseudotime analysis identified a possible infection-specific cellular trajectory for Chlamydia-infected cells, while differential expression analyses found temporal expression of metallothioneins and genes involved with cell cycle regulation, innate immune responses, cytoskeletal components, lipid biosynthesis and cellular stress. We find that changes to the host cell transcriptome at early times of C. trachomatis infection are readily discernible by scRNA-Seq, supporting the utility of single cell approaches to identify host cell biomarkers of chlamydial infection, and to further deconvolute the complex host response to infection.

Ultrahigh-Throughput Multiplexing and Sequencing of >500-Base-Pair Amplicon Regions on the Illumina HiSeq 2500 Platform.

Amplification, sequencing, and analysis of the 16S rRNA gene affords characterization of microbial community composition. As this tool has become more popular and amplicon-sequencing applications have grown in the total number of samples, growth in sample multiplexing is becoming necessary while maintaining high sequence quality and sequencing depth. Here, modifications to the Illumina HiSeq 2500 platform are described which produce greater multiplexing capabilities and 300-bp paired-end reads of higher quality than those produced by the current Illumina MiSeq platform. To improve the feasibility and flexibility of this method, a 2-step PCR amplification protocol is also described that allows for targeting of different amplicon regions, and enhances amplification success from samples with low bacterial bioburden. IMPORTANCE Amplicon sequencing has become a popular and widespread tool for surveying microbial communities. Lower overall costs associated with high-throughput sequencing have made it a widely adopted approach, especially for projects that necessitate sample multiplexing to eliminate batch effect and reduced time to acquire data. The method for amplicon sequencing on the Illumina HiSeq 2500 platform described here provides improved multiplexing capabilities while simultaneously producing greater quality sequence data and lower per-sample cost relative to those of the Illumina MiSeq platform without sacrificing amplicon length. To make this method more flexible for various amplicon-targeted regions as well as improve amplification from low-biomass samples, we also present and validate a 2-step PCR library preparation method.

The Cervicovaginal Microbiota-Host Interaction Modulates Chlamydia trachomatis Infection.

The mechanism(s) by which Lactobacillus-dominated cervicovaginal microbiota provide a barrier to Chlamydia trachomatis infection remain(s) unknown. Here we evaluate the impact of different Lactobacillus spp. identified via culture-independent metataxonomic analysis of C. trachomatis-infected women on C. trachomatis infection in a three-dimensional (3D) cervical epithelium model. Lactobacillus spp. that specifically produce d(-) lactic acid were associated with long-term protection against C. trachomatis infection, consistent with reduced protection associated with Lactobacillus iners, which does not produce this isoform, and with decreased epithelial cell proliferation, consistent with the observed prolonged protective effect. Transcriptomic analysis revealed that epigenetic modifications involving histone deacetylase-controlled pathways are integral to the cross talk between host and microbiota. These results highlight a fundamental mechanism whereby the cervicovaginal microbiota modulates host functions to protect against C. trachomatis infection.IMPORTANCE The vaginal microbiota is believed to protect women against Chlamydia trachomatis, the etiologic agent of the most prevalent sexually transmitted infection (STI) in developed countries. The mechanism underlying this protection has remained elusive. Here, we reveal the comprehensive strategy by which the cervicovaginal microbiota modulates host functions to protect against chlamydial infection, thereby providing a novel conceptual mechanistic understanding. Major implications of this work are that (i) the impact of the vaginal microbiota on the epithelium should be considered in future studies of chlamydial infection and other STIs and (ii) a fundamental understanding of the cervicovaginal microbiota's role in protection against STIs may enable the development of novel microbiome-based therapeutic strategies to protect women from infection and improve vaginal and cervical health.

A Laboratory Methodology for Dual RNA-Sequencing of Bacteria and their Host Cells In Vitro.

Dual RNA-Sequencing leverages established next-generation sequencing (NGS)-enabled RNA-Seq approaches to measure genome-wide transcriptional changes of both an infecting bacteria and host cells. By simultaneously investigating both organisms from the same biological sample, dual RNA-Seq can provide unique insight into bacterial infection processes and reciprocal host responses at once. However, the difficulties involved in handling both prokaryotic and eukaryotic material require distinct, optimized procedures. We previously developed and applied dual RNA-Seq to measure prokaryotic and eukaryotic expression profiles of human cells infected with bacteria, using in vitro Chlamydia-infected epithelial cells as proof of principle. Here we provide a detailed laboratory protocol for in vitro dual RNA-Seq that is readily adaptable to any host-bacteria system of interest.

The effect of cigarette smoking on the oral and nasal microbiota.

BACKGROUND: The goal of the study was to investigate whether cigarette smoking alters oral and nasal microbial diversity, composition, and structure. Twenty-three current smokers and 20 never smokers were recruited. From each subject, nine samples including supra and subgingiva plaque scrapes, saliva, swabs from five soft oral tissue sites, and one nasal swab from both the anterior nares were collected. 16S rRNA V3-V4 region was sequenced for microbial profiles. RESULTS: We found that alpha diversity was lower in smokers than in nonsmokers in the buccal mucosa, but in other sample sites, microbial diversity and composition were not significantly different by smoking status. Microbial profiles differed significantly among eight oral sites. CONCLUSIONS: This study investigates the effect of cigarette smoking on different sites of the oral cavity and shows a potential effect of cigarette smoking on the buccal mucosa microbiota. The marked heterogeneity of the oral microbial ecosystem that we found may contribute to the stability of the oral microbiota in most sites when facing environmental perturbations such as that caused by cigarette smoking.

Molecular Characterization of the Human Stomach Microbiota in Gastric Cancer Patients.

Helicobacter pylori (Hp) is the primary cause of gastric cancer but we know little of its relative abundance and other microbes in the stomach, especially at the time of gastric cancer diagnosis. Here we characterized the taxonomic and derived functional profiles of gastric microbiota in two different sets of gastric cancer patients, and compared them with microbial profiles in other body sites. Paired non-malignant and tumor tissues were sampled from 160 gastric cancer patients with 80 from China and 80 from Mexico. The 16S rRNA gene V3-V4 region was sequenced using MiSeq platform for taxonomic profiles. PICRUSt was used to predict functional profiles. Human Microbiome Project was used for comparison. We showed that Hp is the most abundant member of gastric microbiota in both Chinese and Mexican samples (51 and 24%, respectively), followed by oral-associated bacteria. Taxonomic (phylum-level) profiles of stomach microbiota resembled oral microbiota, especially when the Helicobacter reads were removed. The functional profiles of stomach microbiota, however, were distinct from those found in other body sites and had higher inter-subject dissimilarity. Gastric microbiota composition did not differ by Hp colonization status or stomach anatomic sites, but did differ between paired non-malignant and tumor tissues in either Chinese or Mexican samples. Our study showed that Hp is the dominant member of the non-malignant gastric tissue microbiota in many gastric cancer patients. Our results provide insights on the gastric microbiota composition and function in gastric cancer patients, which may have important clinical implications.

Enhanced Trapping of HIV-1 by Human Cervicovaginal Mucus Is Associated with Lactobacillus crispatus-Dominant Microbiota.

UNLABELLED: Cervicovaginal mucus (CVM) can provide a barrier that precludes HIV and other sexually transmitted virions from reaching target cells in the vaginal epithelium, thereby preventing or reducing infections. However, the barrier properties of CVM differ from woman to woman, and the causes of these variations are not yet well understood. Using high-resolution particle tracking of fluorescent HIV-1 pseudoviruses, we found that neither pH nor Nugent scores nor total lactic acid levels correlated significantly with virus trapping in unmodified CVM from diverse donors. Surprisingly, HIV-1 was generally trapped in CVM with relatively high concentrations of d-lactic acid and a Lactobacillus crispatus-dominant microbiota. In contrast, a substantial fraction of HIV-1 virions diffused rapidly through CVM with low concentrations of d-lactic acid that had a Lactobacillus iners-dominant microbiota or significant amounts of Gardnerella vaginalis, a bacterium associated with bacterial vaginosis. Our results demonstrate that the vaginal microbiota, including specific species of Lactobacillus, can alter the diffusional barrier properties of CVM against HIV and likely other sexually transmitted viruses and that these microbiota-associated changes may account in part for the elevated risks of HIV acquisition linked to bacterial vaginosis or intermediate vaginal microbiota. IMPORTANCE: Variations in the vaginal microbiota, especially shifts away from Lactobacillus-dominant microbiota, are associated with differential risks of acquiring HIV or other sexually transmitted infections. However, emerging evidence suggests that Lactobacillus iners frequently colonizes women with recurring bacterial vaginosis, raising the possibility that L. iners may not be as protective as other Lactobacillus species. Our study was designed to improve understanding of how the cervicovaginal mucus barrier against HIV may vary between women along with the vaginal microbiota and led to the finding that the vaginal microbiota, including specific species of Lactobacillus, can directly alter the diffusional barrier properties of cervicovaginal mucus. This work advances our understanding of the complex barrier properties of mucus and highlights the differential protective ability of different species of Lactobacillus, with Lactobacillus crispatus and possibly other species playing a key role in protection against HIV and other sexually transmitted infections. These findings could lead to the development of novel strategies to protect women against HIV.

Chlamydia caviae infection alters abundance but not composition of the guinea pig vaginal microbiota.

In humans, the vaginal microbiota is thought to be the first line of defense again pathogens including Chlamydia trachomatis. The guinea pig has been extensively used as a model to study chlamydial infection because it shares anatomical and physiological similarities with humans, such as a squamous vaginal epithelium as well as some of the long-term outcomes caused by chlamydial infection. In this study, we aimed to evaluate the guinea pig-C. caviae model of genital infection as a surrogate for studying the role of the vaginal microbiota in the early steps of C. trachomatis infection in humans. We used culture-independent molecular methods to characterize the relative and absolute abundance of bacterial phylotypes in the guinea pig vaginal microbiota in animals non-infected, mock-infected or infected by C. caviae. We showed that the guinea pig and human vaginal microbiotas are of different bacterial composition and abundance. Chlamydia caviae infection had a profound effect on the absolute abundance of bacterial phylotypes but not on the composition of the guinea pig vaginal microbiota. Our findings compromise the validity of the guinea pig-C. caviae model to study the role of the vaginal microbiota during the early steps of sexually transmitted infection.

Environmental surveillance for toxigenic Vibrio cholerae in surface waters of Haiti.

Epidemic cholera was reported in Haiti in 2010, with no information available on the occurrence or geographic distribution of toxigenic Vibrio cholerae in Haitian waters. In a series of field visits conducted in Haiti between 2011 and 2013, water and plankton samples were collected at 19 sites. Vibrio cholerae was detected using culture, polymerase chain reaction, and direct viable count methods (DFA-DVC). Cholera toxin genes were detected by polymerase chain reaction in broth enrichments of samples collected in all visits except March 2012. Toxigenic V. cholerae was isolated from river water in 2011 and 2013. Whole genome sequencing revealed that these isolates were a match to the outbreak strain. The DFA-DVC tests were positive for V. cholerae O1 in plankton samples collected from multiple sites. Results of this survey show that toxigenic V. cholerae could be recovered from surface waters in Haiti more than 2 years after the onset of the epidemic.

Draft Genome Sequences of Nine Enteropathogenic Escherichia coli Strains from Kenya.

We report here the draft genome sequences of nine enteropathogenic Escherichia coli (EPEC) strains isolated from children in Kenya who died during hospitalization with diarrhea. Each of the isolates possess the EPEC adherence factor (EAF) plasmid encoding the bundle-forming pilus, which is characteristic of EPEC. These isolates represent diverse serogroups and EPEC phylogenomic lineages.