Search

Abstract Background Lifestyle plays an important role in shaping the gut microbiome. However, its contributions to the oral microbiome remain less clear, due to the confounding effects of geography and methodology in investigations of populations studied to date. Furthermore, while the oral microbiome seems to differ between foraging and industrialized populations, we lack insight into whether transitions to and away from agrarian lifestyles shape the oral microbiota. Given the growing interest in so-called “vanishing microbiomes” potentially being a risk factor for increased disease prevalence in industrialized populations, it is important that we distinguish lifestyle from geography in the study of microbiomes across populations. Results Here, we investigate salivary microbiomes of 63 Nepali individuals representing a spectrum of lifestyles: foraging, subsistence farming (individuals that transitioned from foraging to farming within the last 50 years), agriculturalists (individuals that have transitioned to farming for at least 300 years), and industrialists (expatriates that immigrated to the USA within the last 20 years). We characterize the role of lifestyle in microbial diversity, identify microbes that differ between lifestyles, and pinpoint specific lifestyle factors that may be contributing to differences in the microbiomes across populations. Contrary to prevailing views, when geography is controlled for, oral microbiome alpha diversity does not differ significantly across lifestyles. Microbiome composition, however, follows the gradient of lifestyles from foraging through agrarianism to industrialism, supporting the notion that lifestyle indeed plays a role in the oral microbiome. Relative abundances of several individual taxa, including Streptobacillus and an unclassified Porphyromonadaceae genus, also mirror lifestyle. Finally, we identify specific lifestyle factors associated with microbiome composition across the gradient of lifestyles, including smoking and grain sources. Conclusion Our findings demonstrate that by studying populations within Nepal, we can isolate an important role of lifestyle in determining oral microbiome composition. In doing so, we highlight the potential contributions of several lifestyle factors, underlining the importance of carefully examining the oral microbiome across lifestyles to improve our understanding of global microbiomes. Video Abstract

[This corrects the article DOI: 10.1371/journal.ppat.1005943.].

Our work focuses on the stability, resilience, and response to perturbation of the bacterial communities in the human gut. Informative flash flood-like disturbances that eliminate most gastrointestinal biomass can be induced using a clinically-relevant iso-osmotic agent. We designed and executed such a disturbance in human volunteers using a dense longitudinal sampling scheme extending before and after induced diarrhea. This experiment has enabled a careful multidomain analysis of a controlled perturbation of the human gut microbiota with a new level of resolution. These new longitudinal multidomain data were analyzed using recently developed statistical methods that demonstrate improvements over current practices. By imposing sparsity constraints we have enhanced the interpretability of the analyses and by employing a new adaptive generalized principal components analysis, incorporated modulated phylogenetic information and enhanced interpretation through scoring of the portions of the tree most influenced by the perturbation. Our analyses leverage the taxa-sample duality in the data to show how the gut microbiota recovers following this perturbation. Through a holistic approach that integrates phylogenetic, metagenomic and abundance information, we elucidate patterns of taxonomic and functional change that characterize the community recovery process across individuals. We provide complete code and illustrations of new sparse statistical methods for high-dimensional, longitudinal multidomain data that provide greater interpretability than existing methods.

The WHO estimates around a million children contract tuberculosis (TB) annually with over 80 000 deaths from dissemination of infection outside of the lungs. The insidious onset and association with skin test anergy suggests failure of the immune system to both recognise and respond to infection. To understand the immune mechanisms, we studied genome-wide whole blood RNA expression in children with TB meningitis (TBM). Findings were validated in a second cohort of children with TBM and pulmonary TB (PTB), and functional T-cell responses studied in a third cohort of children with TBM, other extrapulmonary TB (EPTB) and PTB. The predominant RNA transcriptional response in children with TBM was decreased abundance of multiple genes, with 140/204 (68%) of all differentially regulated genes showing reduced abundance compared to healthy controls. Findings were validated in a second cohort with concordance of the direction of differential expression in both TBM (r2 = 0.78 p = 2x10-16) and PTB patients (r2 = 0.71 p = 2x10-16) when compared to a second group of healthy controls. Although the direction of expression of these significant genes was similar in the PTB patients, the magnitude of differential transcript abundance was less in PTB than in TBM. The majority of genes were involved in activation of leucocytes (p = 2.67E-11) and T-cell receptor signalling (p = 6.56E-07). Less abundant gene expression in immune cells was associated with a functional defect in T-cell proliferation that recovered after full TB treatment (p

Immunocompromised individuals are at high risk of developing Clostridium difficile-associated disease (CDAD). Fecal microbiota transplantation (FMT) is a highly effective therapy for refractory or recurrent CDAD and, despite safety concerns, has recently been offered to immunocompromised patients. We investigated the genomics of bacterial composition following FMT in immunocompromised patients over a 1-year period. Metagenomic, strain and gene-level bacterial dynamics were characterized in two CDAD-affected hematopoietic stem cell (HCT) recipients following FMT. We found alterations in gene content, including loss of virulence and antibiotic resistance genes. These alterations were accompanied by long-term bacterial divergence at the species and strain levels. Our findings suggest limited durability of the specific bacterial consortium introduced with FMT and indicate that alterations of the functional potential of the microbiome are more complex than can be inferred by taxonomic information alone. Our observation that FMT alone cannot induce long-term donor-like alterations of the microbiota of HCT recipients suggests that FMT cannot indefinitely supersede environmental and/or host factors in shaping bacterial composition.

Leptospirosis causes significant morbidity and mortality worldwide; however, the role of the host immune response in disease progression and high case fatality (>10-50%) is poorly understood. We conducted a multi-parameter investigation of patients with acute leptospirosis to identify mechanisms associated with case fatality. Whole blood transcriptional profiling of 16 hospitalized Brazilian patients with acute leptospirosis (13 survivors, 3 deceased) revealed fatal cases had lower expression of the antimicrobial peptide, cathelicidin, and chemokines, but more abundant pro-inflammatory cytokine receptors. In contrast, survivors generated strong adaptive immune signatures, including transcripts relevant to antigen presentation and immunoglobulin production. In an independent cohort (23 survivors, 22 deceased), fatal cases had higher bacterial loads (P = 0.0004) and lower anti-Leptospira antibody titers (P = 0.02) at the time of hospitalization, independent of the duration of illness. Low serum cathelicidin and RANTES levels during acute illness were independent risk factors for higher bacterial loads (P = 0.005) and death (P = 0.04), respectively. To investigate the mechanism of cathelicidin in patients surviving acute disease, we administered LL-37, the active peptide of cathelicidin, in a hamster model of lethal leptospirosis and found it significantly decreased bacterial loads and increased survival. Our findings indicate that the host immune response plays a central role in severe leptospirosis disease progression. While drawn from a limited study size, significant conclusions include that poor clinical outcomes are associated with high systemic bacterial loads, and a decreased antibody response. Furthermore, our data identified a key role for the antimicrobial peptide, cathelicidin, in mounting an effective bactericidal response against the pathogen, which represents a valuable new therapeutic approach for leptospirosis.

The development of a vaccine against dengue faces unique challenges, including the complexity of the immune responses to the four antigenically distinct serotypes. Genome-wide transcriptional profiling provides insight into the pathways and molecular features that underlie responses to immune system stimulation, and may facilitate predictions of immune protection.In this study, we measured early transcriptional responses in the peripheral blood of cynomolgus macaques following vaccination with a live, attenuated tetravalent dengue vaccine candidate, TDV, which is based on a DENV-2 backbone. Different doses and routes of vaccine administration were used, and viral load and neutralizing antibody titers were measured at different time-points following vaccination. All 30 vaccinated animals developed a neutralizing antibody response to each of the four dengue serotypes, and only 3 of these animals had detectable serum viral RNA after challenge with wild-type dengue virus (DENV), suggesting protection of vaccinated animals to DENV infection. The vaccine induced statistically significant changes in 595 gene transcripts on days 1, 3, 5 and 7 as compared with baseline and placebo-treated animals. Genes involved in the type I interferon (IFN) response, including IFI44, DDX58, MX1 and OASL, exhibited the highest fold-change in transcript abundance, and this response was strongest following double dose and subcutaneous (versus intradermal) vaccine administration. In addition, modules of genes involved in antigen presentation, dendritic cell activation, and T cell activation and signaling were enriched following vaccination. Increased abundance of gene transcripts related to T cell activation on day 5, and the type I IFN response on day 7, were significantly correlated with the development of high neutralizing antibody titers on day 30.These results suggest that early transcriptional responses may be predictive of development of adaptive immunity to TDV vaccination in cynomolgus macaques, and will inform studies of human responses to dengue vaccines.

Genetics, diet, and other environmental exposures are thought to be major factors in the development and composition of the intestinal microbiota of animals. However, the relative contributions of these factors in adult animals, as well as variation with time in a variety of important settings, are still not fully understood. We studied a population of inbred, female mice fed the same diet and housed under the same conditions. We collected fecal samples from 46 individual mice over two weeks, sampling four of these mice for periods as long as 236 days for a total of 190 samples, and determined the phylogenetic composition of their microbial communities after analyzing 1,849,990 high-quality pyrosequencing reads of the 16S rRNA gene V3 region. Even under these controlled conditions, we found significant inter-individual variation in community composition, as well as variation within an individual over time, including increases in alpha diversity during the first 2 months of co-habitation. Some variation was explained by mouse membership in different cage and vendor shipment groups. The differences among individual mice from the same shipment group and cage were still significant. Overall, we found that 23% of the variation in intestinal microbiota composition was explained by changes within the fecal microbiota of a mouse over time, 12% was explained by persistent differences among individual mice, 14% by cage, and 18% by shipment group. Our findings suggest that the microbiota of controlled populations of inbred laboratory animals may not be as uniform as previously thought, that animal rearing and handling may account for some variation, and that as yet unidentified factors may explain additional components of variation in the composition of the microbiota within populations and individuals over time. These findings have implications for the design and interpretation of experiments involving laboratory animals.

Our current understanding of the composition and stability of the human distal gut microbiota is based largely on studies of infants and adults living in developed countries. In contrast, little is known about the gut microbiota and its variation over time in older children and adolescents, especially in developing countries.We compared the diversity, composition, and temporal stability of the fecal microbiota of healthy children, ages 9 to 14 years, living in an urban slum in Bangladesh with that of children of the same age range in an upper-middle class suburban community in the United States. We analyzed >8,000 near full-length 16S rRNA gene sequences and over 845,000 pyrosequencing reads of the 16S rRNA V1-V3 region. The distal gut of Bangladeshi children harbored significantly greater bacterial diversity than that of U.S. children, including novel lineages from several bacterial phyla. Bangladeshi and U.S. children had distinct fecal bacterial community membership and structure; the microbiota of Bangladeshi children was enriched in Prevotella, Butyrivibrio, and Oscillospira and depleted in Bacteroides relative to U.S. children (although similar to Bangladeshi adults). Furthermore, community membership and structure in Bangladeshi children was significantly less stable month-to-month than U.S. children.Together, these results suggest that differing environmental or genetic factors may shape the microbiota of healthy children in the two countries. Further investigation is necessary to understand the mechanisms and factors that underlie these differences, and to incorporate these findings into new strategies for the prevention and treatment of childhood and adolescent diseases.

Dengue is the most prevalent mosquito-borne human illness worldwide. The ability to predict disease severity during the earliest days of the illness is a long-sought, but unachieved goal.We examined human genome-wide transcript abundance patterns in daily peripheral blood mononuclear cell (PBMC) samples from 41 children hospitalized with dengue virus (DENV) infection in Nicaragua, as well as 8 healthy control subjects. Nine patients had primary dengue fever (DF1), 11 had dengue fever with serologic evidence of prior DENV infection, i.e., secondary dengue fever (DF2), 12 had dengue hemorrhagic fever (DHF), and 9 had dengue shock syndrome (DSS). We identified 2,092 genes for which transcript abundance differed significantly between patients on days 3-6 of fever and healthy subjects (FDR

Bordetella pertussis filamentous hemagglutinin (FHA) is a surface-associated and secreted protein that serves as a crucial adherence factor, and displays immunomodulatory activity in human peripheral blood mononuclear cells (PBMCs). In order to appreciate more fully the role of secreted FHA in pathogenesis, we analyzed FHA-induced changes in genome-wide transcript abundance in human PBMCs. Among the 683 known unique genes with greater than 3-fold change in transcript abundance following FHA treatment, 125 (18.3%) were identified as interferon (IFN)-regulated. Among the latter group were genes encoding several members of the IFN type I response, as well as 3 key components of the ISGylation pathway. Using real-time RT-PCR, we confirmed FHA-associated increases in transcript abundance for the genes encoding ubiquitin-like protein, ISG15, and its specific protease USP18. Western-blot analysis demonstrated the presence of both, free ISG15 and several ISGylated conjugates in FHA-stimulated PBMC lysates, but not in unstimulated cells. Intracellular FACS analysis provided evidence that monocytes and a natural killer-enriched cell population were the primary producers of ISG15 in PBMCs after FHA stimulation. Our data reveal previously-unrecognized effects of B. pertussis FHA on host IFN and ISGylation responses, and suggest previously-unsuspected mechanisms by which FHA may alter the outcome of the host-pathogen interaction.

Poxviruses use an arsenal of molecular weapons to evade detection and disarm host immune responses. We used DNA microarrays to investigate the gene expression responses to infection by monkeypox virus (MPV), an emerging human pathogen, and Vaccinia virus (VAC), a widely used model and vaccine organism, in primary human macrophages, primary human fibroblasts and HeLa cells. Even as the overwhelmingly infected cells approached their demise, with extensive cytopathic changes, their gene expression programs appeared almost oblivious to poxvirus infection. Although killed (gamma-irradiated) MPV potently induced a transcriptional program characteristic of the interferon response, no such response was observed during infection with either live MPV or VAC. Moreover, while the gene expression response of infected cells to stimulation with ionomycin plus phorbol 12-myristate 13-acetate (PMA), or poly (I-C) was largely unimpaired by infection with MPV, a cluster of pro-inflammatory genes were a notable exception. Poly(I-C) induction of genes involved in alerting the innate immune system to the infectious threat, including TNF-alpha, IL-1 alpha and beta, CCL5 and IL-6, were suppressed by infection with live MPV. Thus, MPV selectively inhibits expression of genes with critical roles in cell-signaling pathways that activate innate immune responses, as part of its strategy for stealthy infection.

Interferons are key modulators of the immune system, and are central to the control of many diseases. The response of immune cells to stimuli in complex populations is the product of direct and indirect effects, and of homotypic and heterotypic cell interactions. Dissecting the global transcriptional profiles of immune cell populations may provide insights into this regulatory interplay. The host transcriptional response may also be useful in discriminating between disease states, and in understanding pathophysiology. The transcriptional programs of cell populations in health therefore provide a paradigm for deconvoluting disease-associated gene expression profiles.We used human cDNA microarrays to (1) compare the gene expression programs in human peripheral blood mononuclear cells (PBMCs) elicited by 6 major mediators of the immune response: interferons alpha, beta, omega and gamma, IL12 and TNFalpha; and (2) characterize the transcriptional responses of purified immune cell populations (CD4+ and CD8+ T cells, B cells, NK cells and monocytes) to IFNgamma stimulation. We defined a highly stereotyped response to type I interferons, while responses to IFNgamma and IL12 were largely restricted to a subset of type I interferon-inducible genes. TNFalpha stimulation resulted in a distinct pattern of gene expression. Cell type-specific transcriptional programs were identified, highlighting the pronounced response of monocytes to IFNgamma, and emergent properties associated with IFN-mediated activation of mixed cell populations. This information provides a detailed view of cellular activation by immune mediators, and contributes an interpretive framework for the definition of host immune responses in a variety of disease settings.

The human intestinal microbiota is essential to the health of the host and plays a role in nutrition, development, metabolism, pathogen resistance, and regulation of immune responses. Antibiotics may disrupt these coevolved interactions, leading to acute or chronic disease in some individuals. Our understanding of antibiotic-associated disturbance of the microbiota has been limited by the poor sensitivity, inadequate resolution, and significant cost of current research methods. The use of pyrosequencing technology to generate large numbers of 16S rDNA sequence tags circumvents these limitations and has been shown to reveal previously unexplored aspects of the "rare biosphere." We investigated the distal gut bacterial communities of three healthy humans before and after treatment with ciprofloxacin, obtaining more than 7,000 full-length rRNA sequences and over 900,000 pyrosequencing reads from two hypervariable regions of the rRNA gene. A companion paper in PLoS Genetics (see Huse et al., doi: 10.1371/journal.pgen.1000255) shows that the taxonomic information obtained with these methods is concordant. Pyrosequencing of the V6 and V3 variable regions identified 3,300-5,700 taxa that collectively accounted for over 99% of the variable region sequence tags that could be obtained from these samples. Ciprofloxacin treatment influenced the abundance of about a third of the bacterial taxa in the gut, decreasing the taxonomic richness, diversity, and evenness of the community. However, the magnitude of this effect varied among individuals, and some taxa showed interindividual variation in the response to ciprofloxacin. While differences of community composition between individuals were the largest source of variability between samples, we found that two unrelated individuals shared a surprising degree of community similarity. In all three individuals, the taxonomic composition of the community closely resembled its pretreatment state by 4 weeks after the end of treatment, but several taxa failed to recover within 6 months. These pervasive effects of ciprofloxacin on community composition contrast with the reports by participants of normal intestinal function and with prior assumptions of only modest effects of ciprofloxacin on the intestinal microbiota. These observations support the hypothesis of functional redundancy in the human gut microbiota. The rapid return to the pretreatment community composition is indicative of factors promoting community resilience, the nature of which deserves future investigation.

Massively parallel pyrosequencing of hypervariable regions from small subunit ribosomal RNA (SSU rRNA) genes can sample a microbial community two or three orders of magnitude more deeply per dollar and per hour than capillary sequencing of full-length SSU rRNA. As with full-length rRNA surveys, each sequence read is a tag surrogate for a single microbe. However, rather than assigning taxonomy by creating gene trees de novo that include all experimental sequences and certain reference taxa, we compare the hypervariable region tags to an extensive database of rRNA sequences and assign taxonomy based on the best match in a Global Alignment for Sequence Taxonomy (GAST) process. The resulting taxonomic census provides information on both composition and diversity of the microbial community. To determine the effectiveness of using only hypervariable region tags for assessing microbial community membership, we compared the taxonomy assigned to the V3 and V6 hypervariable regions with the taxonomy assigned to full-length SSU rRNA sequences isolated from both the human gut and a deep-sea hydrothermal vent. The hypervariable region tags and full-length rRNA sequences provided equivalent taxonomy and measures of relative abundance of microbial communities, even for tags up to 15% divergent from their nearest reference match. The greater sampling depth per dollar afforded by massively parallel pyrosequencing reveals many more members of the "rare biosphere" than does capillary sequencing of the full-length gene. In addition, tag sequencing eliminates cloning bias and the sequences are short enough to be completely sequenced in a single read, maximizing the number of organisms sampled in a run while minimizing chimera formation. This technique allows the cost-effective exploration of changes in microbial community structure, including the rare biosphere, over space and time and can be applied immediately to initiatives, such as the Human Microbiome Project.

Preterm delivery causes substantial neonatal mortality and morbidity. Unrecognized intra-amniotic infections caused by cultivation-resistant microbes may play a role. Molecular methods can detect, characterize and quantify microbes independently of traditional culture techniques. However, molecular studies that define the diversity and abundance of microbes invading the amniotic cavity, and evaluate their clinical significance within a causal framework, are lacking.In parallel with culture, we used broad-range end-point and real-time PCR assays to amplify, identify and quantify ribosomal DNA (rDNA) of bacteria, fungi and archaea from amniotic fluid of 166 women in preterm labor with intact membranes. We sequenced up to 24 rRNA clones per positive specimen and assigned taxonomic designations to approximately the species level. Microbial prevalence, diversity and abundance were correlated with host inflammation and with gestational and neonatal outcomes. Study subjects who delivered at term served as controls. The combined use of molecular and culture methods revealed a greater prevalence (15% of subjects) and diversity (18 taxa) of microbes in amniotic fluid than did culture alone (9.6% of subjects; 11 taxa). The taxa detected only by PCR included a related group of fastidious bacteria, comprised of Sneathia sanguinegens, Leptotrichia amnionii and an unassigned, uncultivated, and previously-uncharacterized bacterium; one or more members of this group were detected in 25% of positive specimens. A positive PCR was associated with histologic chorioamnionitis (adjusted odds ratio [OR] 20; 95% CI, 2.4 to 172), and funisitis (adjusted OR 18; 95% CI, 3.1 to 99). The positive predictive value of PCR for preterm delivery was 100 percent. A temporal association between a positive PCR and delivery was supported by a shortened amniocentesis-to-delivery interval (adjusted hazard ratio 4.6; 95% CI, 2.2 to 9.5). A dose-response association was demonstrated between bacterial rDNA abundance and gestational age at delivery (r(2) = 0.42; P

Poxviruses engage in a complex and intricate dialogue with host cells as part of their strategy for replication. However, relatively little molecular detail is available with which to understand the mechanisms behind this dialogue.We designed a specialized microarray that contains probes specific to all predicted ORFs in the Monkeypox Zaire (MPXV) and Vaccinia Western Reserve (VACV) genomes, as well as >18,000 human genes, and used this tool to characterize MPXV and VACV gene expression responses in vitro during the course of primary infection of human monocytes, primary human fibroblasts and HeLa cells. The two viral transcriptomes show distinct features of temporal regulation and species-specific gene expression, and provide an early foundation for understanding global gene expression responses during poxvirus infection.The results provide a temporal map of the transcriptome of each virus during infection, enabling us to compare viral gene expression across species, and classify expression patterns of previously uncharacterized ORFs.

Filamentous hemagglutinin (FHA) is a cell-associated and secreted adhesin produced by Bordetella pertussis with pro-apoptotic and pro-inflammatory activity in host cells. Given the importance of the NF-kappaB transcription factor family in these host cell responses, we examined the effect of FHA on NF-kappaB activation in macrophages and bronchial epithelial cells, both of which are relevant cell types during natural infection.Exposure to FHA of primary human monocytes and transformed U-937 macrophages, but not BEAS-2B epithelial cells, resulted in early activation of the NF-kappaB pathway, as manifested by the degradation of cytosolic IkappaB alpha, by NF-kappaB DNA binding, and by the subsequent secretion of NF-kappaB-regulated inflammatory cytokines. However, exposure of macrophages and human monocytes to FHA for two hours or more resulted in the accumulation of cytosolic IkappaB alpha, and the failure of TNF-alpha to activate NF-kappaB. Proteasome activity was attenuated following exposure of cells to FHA for 2 hours, as was the nuclear translocation of RelA in BEAS-2B cells.These results reveal a complex temporal dynamic, and suggest that despite short term effects to the contrary, longer exposures of host cells to this secreted adhesin may block NF-kappaB activation, and perhaps lead to a compromised immune response to this bacterial pathogen.

Almost immediately after a human being is born, so too is a new microbial ecosystem, one that resides in that person's gastrointestinal tract. Although it is a universal and integral part of human biology, the temporal progression of this process, the sources of the microbes that make up the ecosystem, how and why it varies from one infant to another, and how the composition of this ecosystem influences human physiology, development, and disease are still poorly understood. As a step toward systematically investigating these questions, we designed a microarray to detect and quantitate the small subunit ribosomal RNA (SSU rRNA) gene sequences of most currently recognized species and taxonomic groups of bacteria. We used this microarray, along with sequencing of cloned libraries of PCR-amplified SSU rDNA, to profile the microbial communities in an average of 26 stool samples each from 14 healthy, full-term human infants, including a pair of dizygotic twins, beginning with the first stool after birth and continuing at defined intervals throughout the first year of life. To investigate possible origins of the infant microbiota, we also profiled vaginal and milk samples from most of the mothers, and stool samples from all of the mothers, most of the fathers, and two siblings. The composition and temporal patterns of the microbial communities varied widely from baby to baby. Despite considerable temporal variation, the distinct features of each baby's microbial community were recognizable for intervals of weeks to months. The strikingly parallel temporal patterns of the twins suggested that incidental environmental exposures play a major role in determining the distinctive characteristics of the microbial community in each baby. By the end of the first year of life, the idiosyncratic microbial ecosystems in each baby, although still distinct, had converged toward a profile characteristic of the adult gastrointestinal tract.

ABSTRACT The vaginal microbiome has been linked to negative health outcomes including preterm birth. Specific taxa, including Gardnerella spp., have been identified as risk factors for these conditions. Historically, microbiome analysis methods have treated all Gardnerella spp. as one species, but the broad diversity of Gardnerella has become more apparent. We explore the diversity of Gardnerella clades and genomic species in the vaginal microbiome of pregnant women and their associations with microbiome composition and preterm birth. Relative abundance of Gardnerella clades and genomic species and other taxa was quantified in shotgun metagenomic sequencing data from three distinct cohorts of pregnant women. We also assessed the diversity and abundance of Gardnerella variants in 16S rRNA gene amplicon sequencing data from seven previously conducted studies in differing populations. Individual microbiomes often contained multiple Gardnerella variants, and the number of clades was associated with increased microbial load, or the ratio of non-human reads to human reads. Taxon co-occurrence patterns were largely consistent across Gardnerella clades and among cohorts. Some variants previously described as rare were prevalent in other cohorts, highlighting the importance of surveying a diverse set of populations to fully capture the diversity of Gardnerella. The diversity of Gardnerella both across populations and within individual vaginal microbiomes has long been unappreciated, as has been the intra-species diversity of many other members of the vaginal microbiome. The broad genomic diversity of Gardnerella has led to its reclassification as multiple species; here we demonstrate the diversity of Gardnerella found within and between vaginal microbiomes.IMPORTANCEThe present study shows that single microbiomes can contain all currently known species of Gardnerella and that multiple similar species can exist within the same environment. Furthermore, surveys of demographically distinct populations suggest that some species appear more commonly in certain populations. Further studies in broad and diverse populations will be necessary to fully understand the ecological roles of each Gardnerella sp., how they can co-exist, and their distinct impacts on microbial communities, preterm birth, and other health outcomes.

Bordetella pertussis, B. bronchiseptica, B. parapertussis(hu), and B. parapertussis(ov) are closely related respiratory pathogens that infect mammalian species. B. pertussis and B. parapertussis(hu) are exclusively human pathogens and cause whooping cough, or pertussis, a disease that has resurged despite vaccination. Although it most often infects animals, infrequently B. bronchiseptica is isolated from humans, and these infections are thought to be zoonotic. B. pertussis and B. parapertussis(hu) are assumed to have evolved from a B. bronchiseptica-like ancestor independently. To determine the phylogenetic relationships among these species, housekeeping and virulence genes were sequenced, comparative genomic hybridizations were performed using DNA microarrays, and the distribution of insertion sequence elements was determined, using a collection of 132 strains. This multifaceted approach distinguished four complexes, representing B. pertussis, B. parapertussis(hu), and two distinct B. bronchiseptica subpopulations, designated complexes I and IV. Of the two B. bronchiseptica complexes, complex IV was more closely related to B. pertussis. Of interest, while only 32% of the complex I strains were isolated from humans, 80% of the complex IV strains were human isolates. Comparative genomic hybridization analysis identified the absence of the pertussis toxin locus and dermonecrotic toxin gene, as well as a polymorphic lipopolysaccharide biosynthesis locus, as associated with adaptation of complex IV strains to the human host. Lipopolysaccharide structural diversity among these strains was confirmed by gel electrophoresis. Thus, complex IV strains may comprise a human-associated lineage of B. bronchiseptica from which B. pertussis evolved. These findings will facilitate the study of pathogen host-adaptation. Our results shed light on the origins of the disease pertussis and suggest that the association of B. pertussis with humans may be more ancient than previously assumed.

Abstract The vaginal ecosystem is closely tied to human health and reproductive outcomes, yet its dynamics in the wake of childbirth remain poorly characterized. Here, we profile the vaginal microbiota and cytokine milieu of participants sampled longitudinally throughout pregnancy and for at least one year postpartum. We show that delivery, regardless of mode, is associated with a vaginal pro-inflammatory cytokine response and the loss of Lactobacillus dominance. By contrast, neither the progression of gestation nor the approach of labor strongly altered the vaginal ecosystem. At 9.5-months postpartum—the latest timepoint at which cytokines were assessed—elevated inflammation coincided with vaginal bacterial communities that had remained perturbed (highly diverse) from the time of delivery. Time-to-event analysis indicated a one-year postpartum probability of transitioning to Lactobacillus dominance of 49.4%. As diversity and inflammation declined during the postpartum period, dominance by L. crispatus, the quintessential health-associated commensal, failed to return: its prevalence before, immediately after, and one year after delivery was 41%, 4%, and 9%, respectively. Revisiting our pre-delivery data, we found that a prior live birth was associated with a lower odds of L. crispatus dominance in pregnant participants—an outcome modestly tempered by a longer ( > 18-month) interpregnancy interval. Our results suggest that reproductive history and childbirth in particular remodel the vaginal ecosystem and that the timing and degree of recovery from delivery may help determine the subsequent health of the woman and of future pregnancies.