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.

Abstract Much remains to be explored regarding the diversity of uncultured, host-associated microbes. Here, we describe rectangular bacterial structures (RBSs) in the mouths of bottlenose dolphins. DNA staining revealed multiple paired bands within RBSs, suggesting the presence of cells dividing along the longitudinal axis. Cryogenic transmission electron microscopy and tomography showed parallel membrane-bound segments that are likely cells, encapsulated by an S-layer-like periodic surface covering. RBSs displayed unusual pilus-like appendages with bundles of threads splayed at the tips. We present multiple lines of evidence, including genomic DNA sequencing of micromanipulated RBSs, 16S rRNA gene sequencing, and fluorescence in situ hybridization, suggesting that RBSs are bacterial and distinct from the genera Simonsiella and Conchiformibius (family Neisseriaceae), with which they share similar morphology and division patterning. Our findings highlight the diversity of novel microbial forms and lifestyles that await characterization using tools complementary to genomics such as microscopy.

Abstract The oral microbiota is emerging as an influential factor of host physiology and disease state. Factors influencing oral microbiota composition have not been well characterised. In particular, there is a lack of population-based studies. We undertook a large hypothesis-free study of the saliva microbiota, considering potential influential factors of host health (frailty; diet; periodontal disease), demographics (age; sex; BMI) and sample processing (storage time), in a sample (n = 679) of the TwinsUK cohort of adult twins. Alpha and beta diversity of the saliva microbiota was associated most strongly with frailty (alpha diversity: β = −0.16, Q = 0.003, Observed; β = −0.16, Q = 0.002, Shannon; β = −0.16, Q = 0.003, Simpson; Beta diversity: Q = 0.002, Bray Curtis dissimilarity) and age (alpha diversity: β = 0.15, Q = 0.006, Shannon; β = 0.12, Q = 0.003, Simpson; beta diversity: Q = 0.002, Bray Curtis dissimilarity; Q = 0.032, Weighted UniFrac) in multivariate models including age, frailty, sex, BMI, frailty and diet, and adjustment for multiple testing. Those with a more advanced age were more likely to be dissimilar in the saliva microbiota composition than younger participants (P = 5.125e−06, ANOVA). In subsample analyses, including consideration of periodontal disease (total n = 138, periodontal disease n = 66), the association with frailty remained for alpha diversity (Q = 0.002, Observed ASVs; Q = 0.04 Shannon Index), but not beta diversity, whilst age was not demonstrated to associate with alpha or beta diversity in this subsample, potentially due to insufficient statistical power. Length of time that samples were stored prior to sequencing was associated with beta diversity (Q = 0.002, Bray Curtis dissimilarity). Six bacterial taxa were associated with age after adjustment for frailty and diet. Of the factors studied, frailty and age emerged as the most influential with regards to saliva microbiota composition. Whilst age and frailty are correlates, the associations were independent of each other, giving precedence to both biological and chronological ageing as processes of potential importance when considering saliva microbiota composition.

Interactions between the immune system and adipose tissue contribute to the regulation of body weight, however, the underlying mechanisms remain incompletely understood. Here the authors dissect the role of two structurally and functionally similar immune mediators, BAFF and APRIL, in modifying diet-induced weight gain and adipocyte lipid handling.

ABSTRACT An outstanding question regarding the human gut microbiota is whether and how microbiota-directed interventions influence host phenotypic traits. Here, we employed a dietary intervention to probe this question in the context of lactose intolerance. To assess the effects of dietary dairy product elimination and (re)introduction on the microbiota and host phenotype, we studied 12 self-reported mildly lactose-intolerant adults with triweekly collection of fecal samples over a 12-week study period: 2 weeks of baseline diet, 4 weeks of dairy product elimination, and 6 weeks of gradual whole cow milk (re)introduction. Of the 12 subjects, 6 reported either no dairy or only lactose-free dairy product consumption. A clinical assay for lactose intolerance, the hydrogen breath test, was performed before and after each of these three study phases, and 16S rRNA gene amplicon sequencing was performed on all fecal samples. We found that none of the subjects showed change in a clinically defined measure of lactose tolerance. Similarly, fecal microbiota structure resisted modification. Although the mean fraction of the genus Bifidobacterium, a group known to metabolize lactose, increased slightly with milk (re)introduction (from 0.0125 to 0.0206; Wilcoxon P = 0.068), the overall structure of each subject’s gut microbiota remained highly individualized and largely stable in the face of diet manipulation. IMPORTANCE Lactose intolerance is a gastrointestinal disorder diagnosed with a lactose hydrogen breath test. Lifestyle changes such as diet interventions can impact the gut microbiome; however, the role of the microbiome in lactose intolerance is unclear. Our study assessed the effects of a 12-week dietary dairy product elimination and (re)introduction on the microbiome and clinical lactose intolerance status in 12 adult self-reported lactose-intolerant individuals. We found each subject’s gut microbiome remained highly individualized and largely stable in the face of this diet manipulation. We also report that none of the subjects showed change in a clinically defined measure of lactose tolerance.

Abstract For more than a century, the prenatal environment was considered sterile. Over the last few years, findings obtained with next-generation sequencing approaches from samples of the placenta, the amniotic fluid, meconium, and even fetal tissues have challenged the dogma of a sterile womb, and additional reports have emerged that used culture, microscopy, and quantitative PCR to support the presence of a low-biomass microbial community at prenatal sites. Given the substantial implications of prenatal exposure to microbes for the development and health of the host, the findings have gathered substantial interest from academics, high impact journals, the public press, and funding agencies. However, an increasing number of studies have challenged the prenatal microbiome identifying contamination as a major issue, and scientists that remained skeptical have pointed to inconsistencies with in utero colonization, the impact of c-sections on early microbiome assembly, and the ability to generate germ-free mammals. A lively academic controversy has emerged on the existence of the wider importance of prenatal microbial communities. Microbiome has asked experts to discuss these issues and provide their thoughts on the implications. To allow for a broader perspective of this discussion, we have specifically selected scientists, who have a long-standing expertise in microbiome sciences but who have not directly been involved in the debate so far.

A key security challenge with biosecurity threats is determining the responsible actor. In this Perspective, the authors review recent developments in using genetic sequence to assign a lab-of-origin and the potential protection it provides against misuse of synthetic biology.

Selective and neutral forces shape human microbiota assembly in early life. Here, Sprockett et al. study microbial community assembly in 47 infant-mother pairs from the Tsimane, an indigenous Bolivian population, highlighting the importance of neutral forces during microbiota assembly.

Abstract The southern sea otter (Enhydra lutris nereis) is a threatened sub‐species in coastal ecosystems. To understand better the role of diet, monitor health, and enhance management of this and other marine mammal species, we characterized the oral (gingival) and distal gut (rectal and fecal) microbiota of 158 wild southern sea otters living off the coast of central California, USA, and 12 captive sea otters, some of which were included in a diet shift experiment. We found that the sea otter fecal microbiota was distinct from that of three other otter species, and that captivity does not significantly alter the community structure of the sea otter gingival or distal gut microbiota. Metagenomic analysis unexpectedly revealed that the majority of sea otter fecal DNA is derived from prey, rather than from indigenous bacteria or host cells as with most other mammals. We speculate that a reduced bacterial biomass in the sea otter gut reflects rapid gut transit time and a particular strategy for foraging and energy harvest. This study establishes a reference for the healthy sea otter microbiota, highlights how a marine lifestyle may shape the mammalian microbiota, and may inform future health assessments and conservation management of sea otter populations.

Abstract Background The accuracy of microbial community surveys based on marker-gene and metagenomic sequencing (MGS) suffers from the presence of contaminants—DNA sequences not truly present in the sample. Contaminants come from various sources, including reagents. Appropriate laboratory practices can reduce contamination, but do not eliminate it. Here we introduce decontam (https://github.com/benjjneb/decontam), an open-source R package that implements a statistical classification procedure that identifies contaminants in MGS data based on two widely reproduced patterns: contaminants appear at higher frequencies in low-concentration samples and are often found in negative controls. Results Decontam classified amplicon sequence variants (ASVs) in a human oral dataset consistently with prior microscopic observations of the microbial taxa inhabiting that environment and previous reports of contaminant taxa. In metagenomics and marker-gene measurements of a dilution series, decontam substantially reduced technical variation arising from different sequencing protocols. The application of decontam to two recently published datasets corroborated and extended their conclusions that little evidence existed for an indigenous placenta microbiome and that some low-frequency taxa seemingly associated with preterm birth were contaminants. Conclusions Decontam improves the quality of metagenomic and marker-gene sequencing by identifying and removing contaminant DNA sequences. Decontam integrates easily with existing MGS workflows and allows researchers to generate more accurate profiles of microbial communities at little to no additional cost.

The decreasing effectiveness of antimicrobial agents is a global public health threat, yet risk factors for community-acquired antimicrobial resistance (CA-AMR) in low-income settings have not been clearly elucidated. Our aim was to identify risk factors for CA-AMR with extended-spectrum β-lactamase (ESBL)–producing organisms among urban-dwelling women in India. We collected microbiological and survey data in an observational study of primigravidae women in a public hospital in Hyderabad, India. We analyzed the data using multivariate logistic and linear regression and found that 7% of 1,836 women had bacteriuria; 48% of isolates were ESBL-producing organisms. Women in the bottom 50th percentile of income distribution were more likely to have bacteriuria (adjusted odds ratio 1.44, 95% CI 0.99–2.10) and significantly more likely to have bacteriuria with ESBL-producing organisms (adjusted odds ratio 2.04, 95% CI 1.17–3.54). Nonparametric analyses demonstrated a negative relationship between the prevalence of ESBL and income.

Oral microbial communities are excellent sites for examining the mechanisms of structural variation in the microbiota. Here, the authors use spatial and ecological analysis to examine over 9000 dental and mucosal microbiota samples. They find that saliva plays a key role in structuring bacterial communities in the mouth.

Abstract Genital mycoplasmas, which can be vertically transmitted, have been implicated in preterm birth, neonatal infections, and chronic lung disease of prematurity. Our prior work uncovered 16S rRNA genes belonging to a novel, as-yet-uncultivated mycoplasma (lineage ‘Mnola’) in the oral cavity of a premature neonate. Here, we characterize the organism’s associated community, growth status, metabolic potential, and population diversity. Sequencing of genomic DNA from the infant’s saliva yielded 1.44 Gbp of high-quality, non-human read data, from which we recovered three essentially complete (including ‘Mnola’) and three partial draft genomes (including Trichomonas vaginalis). The completed 629,409-bp ‘Mnola’ genome (Candidatus Mycoplasma girerdii str. UC-B3) was distinct at the strain level from its closest relative, vaginally-derived Ca. M. girerdii str. VCU-M1, which is also associated with T. vaginalis. Replication rate measurements indicated growth of str. UC-B3 within the infant. Genes encoding surface-associated proteins and restriction-modification systems were especially diverse within and between strains. In UC-B3, the population genetic underpinnings of phase variable expression were evident in vivo. Unique among mycoplasmas, Ca. M. girerdii encodes pyruvate-ferredoxin oxidoreductase and may be sensitive to metronidazole. This study reveals a metabolically unique mycoplasma colonizing a premature neonate, and establishes the value of genome-resolved metagenomics in tracking phase variation.

Iron overload disorder (IOD) affects many wildlife species cared for ex situ. Two of the four rhinoceros species in human care, Sumatran rhinoceros (Dicerorhinus sumatrensis) and black rhinoceros (Diceros bicornis), are susceptible, whereas the other two, white rhinoceros (Ceratotherium simum) and greater one-horned (GOH) rhinoceros (Rhinoceros unicornis), are relatively resistant to IOD. Complex interrelationships exist between mammalian hosts, their indigenous gut microbiota, metabolome, physical condition, and iron availability. The goal of this study was to gain insight into these relationships within the family Rhinocerotidae. Specific objectives were to (1) characterize the gut microbiome and metabolome of four rhinoceros species; (2) compare the microbiome and metabolome of IOD-susceptible and IOD-resistant rhinoceros species; and (3) identify variation in the microbiome and metabolome associated with compromised health or disease in IOD-susceptible rhinoceroses. Fecal samples were collected from 31 rhinoceroses (Sumatran rhinoceros, n = 3; black rhinoceros, n = 6; GOH rhinoceros, n = 9; white rhinoceros, n = 13) located at five facilities, and matched fecal aliquots were processed for microbiome and metabolome analyses using 16S rRNA gene sequencing and nuclear magnetic resonance spectroscopy, respectively. Despite the phylogenetic disparity and dissimilar zoo diets of the hosts, the structure of the fecal microbiota of the two IOD-susceptible rhinoceros species were more closely related to each other than to those of the two IOD-resistant species (Bray–Curtis dissimilarity; IOD-susceptible vs. IOD-resistant p-value < 0.001). In addition, IOD-susceptible rhinoceroses exhibited less microbial diversity than their IOD-resistant relatives (Shannon diversity; p-value < 0.001) which could have health implications. Of note, the black rhinoceros was distinct among the four rhinoceros species with the most divergent fecal metabolome; interestingly, it contained higher concentrations of short chain fatty acids. Neither age nor sex were associated with differences in microbial community composition (p = 0.253 and 0.488, respectively) or fecal metabolomic profile (p = 0.634 and 0.332, respectively). Differences in the distal gut microbiomes between IOD-resistant and IOD-susceptible rhinoceroses support hypotheses that gut microbes play a role in host iron acquisition, and further studies and experiments to test these hypotheses are warranted.

Little is known about the microbiota of marine mammals, despite the crucial ecological roles played by these animals. Here, Bik et al. describe the bacterial communities associated with various body sites in dolphins and sea lions, as well as the microbiota of their dietary fish and adjacent seawater.

Broad-range rDNA polymerase chain reaction (PCR) provides an alternative, cultivation-independent approach for identifying pathogens. In 1995, the Centers for Disease Control and Prevention initiated population-based surveillance for unexplained life-threatening infections (Unexplained Death and Critical Illness Project [UNEX]). To address the causes of UNEX cases, we examined 59 specimens from 46 cases by using broad-range bacterial 16S rDNA PCR and phylogenetic analysis of amplified sequences. Specimens from eight cases yielded sequences from Neisseria meningitidis (cerebrospinal fluid from two patients with meningitis), Streptococcus pneumoniae (cerebrospinal fluid from one patient with meningitis and pleural fluid from two patients with pneumonia), or Stenotrophomonas maltophilia (bone marrow aspirate from one patient with pneumonia). Streptococcus pneumoniae rDNA sequence microheterogeneity was found in one pleural fluid specimen, suggesting the presence of multiple strains. In conclusion, known bacterial pathogens cause some critical illnesses and deaths that fail to be explained with traditional diagnostic methods.

ABSTRACT The publication of scientific information that derives from dual use research of concern (DURC) poses major problems for journals because it brings into conflict the benefits of free access to data and the need to prevent misuse of that information by others. Recently, a group of authors and a major scientific journal addressed the issue of publishing information on a newly discovered, highly lethal toxin that can be delivered to large populations and for which there are no available countermeasures. The journal addressed this conflict by permitting the redaction of information that is normally considered essential for publication. This action establishes a precedent for redaction of sensitive data that also provides an example of responsible scientific publishing. However, this precedent leaves many questions unanswered and suggests a need for a discussion by all stakeholders of scientific information so as to derive normative standards for the publication of DURC.

Complete genomic sequences of microbial pathogens and hosts offer sophisticated new strategies for studying host-pathogen interactions. DNA microarrays exploit primary sequence data to measure transcript levels and detect sequence polymorphisms, for every gene, simultaneously. The design and construction of a DNA microarray for any given microbial genome are straightforward. By monitoring microbial gene expression, one can predict the functions of uncharacterized genes, probe the physiologic adaptations made under various environmental conditions, identify virulence-associated genes, and test the effects of drugs. Similarly, by using host gene microarrays, one can explore host response at the level of gene expression and provide a molecular description of the events that follow infection. Host profiling might also identify gene expression signatures unique for each pathogen, thus providing a novel tool for diagnosis, prognosis, and clinical management of infectious disease.

Rhinosporidium seeberi, a microorganism that can infect the mucosal surfaces of humans and animals, has been classified as a fungus on the basis of morphologic and histochemical characteristics. Using consensus polymerase chain reaction (PCR), we amplified a portion of the R. seeberi 18S rRNA gene directly from infected tissue. Analysis of the aligned sequence and inference of phylogenetic relationships showed that R. seeberi is a protist from a novel clade of parasites that infect fish and amphibians. Fluorescence in situ hybridization and R. seeberi-specific PCR showed that this unique 18S rRNA sequence is also present in other tissues infected with R. seeberi. Our data support the R. seeberi phylogeny recently suggested by another group. R. seeberi is not a classic fungus, but rather the first known human pathogen from the DRIPs clade, a novel clade of aquatic protistan parasites (Ichthyosporea).

Features of a number of important but poorly explained human clinical syndromes strongly indicate a microbial etiology. In these syndromes, the failure of cultivation-dependent microbial detection methods reveals our ignorance of microbial growth requirements. Sequence-based molecular methods, however, offer alternative approaches for microbial identification directly from host specimens found in the setting of unexplained acute illnesses, chronic inflammatory disease, and from anatomic sites that contain commensal microflora. The rapid expansion of genome sequence databases and advances in biotechnology present opportunities and challenges: identification of consensus sequences from which reliable, specific phylogenetic information can be inferred for all taxonomic groups of pathogens, broad-range pathogen identification on the basis of virulence-associated gene families, and use of host gene expression response profiles as specific signatures of microbial infection.

ABSTRACT The purpose of this study was to evaluate the composition and richness of bacterial communities associated with low-birthweight (LBW) infants in relation to host body site, individual, and age. Bacterial 16S rRNA genes from saliva samples, skin swabs, and stool samples collected on postnatal days 8, 10, 12, 15, 18, and 21 from six LBW (five premature) infants were amplified, pyrosequenced, and analyzed within a comparative framework that included analogous data from normal-birthweight (NBW) infants and healthy adults. We found that body site was the primary determinant of bacterial community composition in the LBW infants. However, site specificity depended on postnatal age: saliva and stool compositions diverged over time but were not significantly different until the babies were 15 days old. This divergence was primarily driven by progressive temporal turnover in the distal gut, which proceeded at a rate similar to that of age-matched NBW infants. Neonatal skin was the most adult-like in microbiota composition, while saliva and stool remained the least so. Compositional variation among infants was marked and depended on body site and age. Only the smallest, most premature infant received antibiotics during the study period; this heralded a coexpansion of Pseudomonas aeruginosa and a novel Mycoplasma sp. in the oral cavity of this vaginally delivered, intubated patient. We conclude that concurrent molecular surveillance of multiple body sites in LBW neonates reveals a delayed compositional differentiation of the oral cavity and distal gut microbiota and, in the case of one infant, an abundant, uncultivated oral Mycoplasma sp., recently detected in human vaginal samples. IMPORTANCE Complications of premature birth are the most common cause of neonatal mortality. Colonization by the indigenous microbiota, which begins at delivery, may predispose some high-risk newborns to invasive infection or necrotizing enterocolitis (NEC), and protect others, yet neonatal microbiome dynamics are poorly understood. Here, we present the first cultivation-independent time series tracking microbiota assembly across multiple body sites in a synchronous cohort of hospitalized low-birthweight (LBW) neonates. We take advantage of archived samples and publically available sequence data and compare our LBW infant findings to those from normal-birthweight (NBW) infants and healthy adults. Our results suggest potential windows of opportunity for the dispersal of microbes within and between hosts and support recent findings of substantial baseline spatiotemporal variation in microbiota composition among high-risk newborns.