Longitudinal profiling of the microbiome at four body sites reveals core stability and individualized dynamics during health and disease.

To understand the dynamic interplay between the human microbiome and host during health and disease, we analyzed the microbial composition, temporal dynamics, and associations with host multi-omics, immune, and clinical markers of microbiomes from four body sites in 86 participants over 6 years. We found that microbiome stability and individuality are body-site specific and heavily influenced by the host. The stool and oral microbiome are more stable than the skin and nasal microbiomes, possibly due to their interaction with the host and environment. We identify individual-specific and commonly shared bacterial taxa, with individualized taxa showing greater stability. Interestingly, microbiome dynamics correlate across body sites, suggesting systemic dynamics influenced by host-microbial-environment interactions. Notably, insulin-resistant individuals show altered microbial stability and associations among microbiome, molecular markers, and clinical features, suggesting their disrupted interaction in metabolic disease. Our study offers comprehensive views of multi-site microbial dynamics and their relationship with host health and disease.

Longitudinal profiling of the microbiome at four body sites reveals core stability and individualized dynamics during health and disease.

To understand dynamic interplay between the human microbiome and host during health and disease, we analyzed the microbial composition, temporal dynamics, and associations with host multi-omics, immune and clinical markers of microbiomes from four body sites in 86 participants over six years. We found that microbiome stability and individuality are body-site-specific and heavily influenced by the host. The stool and oral microbiome were more stable than the skin and nasal microbiomes, possibly due to their interaction with the host and environment. Also, we identified individual-specific and commonly shared bacterial taxa, with individualized taxa showing greater stability. Interestingly, microbiome dynamics correlated across body sites, suggesting systemic coordination influenced by host-microbial-environment interactions. Notably, insulin-resistant individuals showed altered microbial stability and associations between microbiome, molecular markers, and clinical features, suggesting their disrupted interaction in metabolic disease. Our study offers comprehensive views of multi-site microbial dynamics and their relationship with host health and disease. Study Highlights: The stability of the human microbiome varies among individuals and body sites.Highly individualized microbial genera are more stable over time.At each of the four body sites, systematic interactions between the environment, the host and bacteria can be detected.Individuals with insulin resistance have lower microbiome stability, a more diversified skin microbiome, and significantly altered host-microbiome interactions.

Microbial glutamate metabolism predicts intravenous cocaine self-administration in diversity outbred mice.

The gut microbiome is thought to play a critical role in the onset and development of psychiatric disorders, including depression and substance use disorder (SUD). To test the hypothesis that the microbiome affects addiction predisposing behaviors and cocaine intravenous self-administration (IVSA) and to identify specific microbes involved in the relationship, we performed 16S rRNA gene sequencing on feces from 228 diversity outbred mice. Twelve open field measures, two light-dark assay measures, one hole board and novelty place preference measure significantly differed between mice that acquired cocaine IVSA (ACQ) and those that failed to acquire IVSA (FACQ). We found that ACQ mice are more active and exploratory and display decreased fear than FACQ mice. The microbial abundances that differentiated ACQ from FACQ mice were an increased abundance of Barnesiella, Ruminococcus, and Robinsoniella and decreased Clostridium IV in ACQ mice. There was a sex-specific correlation between ACQ and microbial abundance, a reduced Lactobacillus abundance in ACQ male mice, and a decreased Blautia abundance in female ACQ mice. The abundance of Robinsoniella was correlated, and Clostridium IV inversely correlated with the number of doses of cocaine self-administered during acquisition. Functional analysis of the microbiome composition of a subset of mice suggested that gut-brain modules encoding glutamate metabolism genes are associated with the propensity to self-administer cocaine. These findings establish associations between the microbiome composition and glutamate metabolic potential and the ability to acquire cocaine IVSA thus indicating the potential translational impact of targeting the gut microbiome or microbial metabolites for treatment of SUD. This article is part of the Special Issue on "Microbiome & the Brain: Mechanisms & Maladies".

Exploratory studies of oral and fecal microbiome in healthy human aging.

Growing evidence has linked an altered host fecal microbiome composition with health status, common chronic diseases, and institutionalization in vulnerable older adults. However, fewer studies have described microbiome changes in healthy older adults without major confounding diseases or conditions, and the impact of aging on the microbiome across different body sites remains unknown. Using 16S ribosomal RNA gene sequencing, we reconstructed the composition of oral and fecal microbiomes in young (23-32; mean = 25 years old) and older (69-94; mean = 77 years old) healthy community-dwelling research subjects. In both body sites, we identified changes in minor bacterial operational taxonomic units (OTUs) between young and older subjects. However, the composition of the predominant bacterial species of the healthy older group in both microbiomes was not significantly different from that of the young cohort, which suggests that dominant bacterial species are relatively stable with healthy aging. In addition, the relative abundance of potentially pathogenic genera, such as Rothia and Mycoplasma, was enriched in the oral microbiome of the healthy older group relative to the young cohort. We also identified several OTUs with a prevalence above 40% and some were more common in young and others in healthy older adults. Differences with aging varied for oral and fecal samples, which suggests that members of the microbiome may be differentially affected by aging in a tissue-specific fashion. This is the first study to investigate both oral and fecal microbiomes in the context of human aging, and provides new insights into interactions between aging and the microbiome within two different clinically relevant sites.

Murine Model for Measuring Effects of Humanized-Dosing of Antibiotics on the Gut Microbiome.

There is a current need for enhancing our insight in the effects of antimicrobial treatment on the composition of human microbiota. Also, the spontaneous restoration of the microbiota after antimicrobial treatment requires better understanding. This is best addressed in well-defined animal models. We here present a model in which immune-competent or neutropenic mice were administered piperacillin-tazobactam (TZP) according to human treatment schedules. Before, during and after the TZP treatment, fecal specimens were longitudinally collected at established intervals over several weeks. Gut microbial taxonomic distribution and abundance were assessed through culture and molecular means during all periods. Non-targeted metabolomics analyses of stool samples using Quadrupole Time of Flight mass spectrometry (QTOF MS) were also applied to determine if a metabolic fingerprint correlated with antibiotic use, immune status, and microbial abundance. TZP treatment led to a 5-10-fold decrease in bacterial fecal viability counts which were not fully restored during post-antibiotic follow up. Two distinct, relatively uniform and reproducible restoration scenarios of microbiota changes were seen in post TZP-treatment mice. Post-antibiotic flora could consist of predominantly Firmicutes or, alternatively, a more diverse mix of taxa. In general, the pre-treatment microbial communities were not fully restored within the screening periods applied. A new species, closely related to Eubacterium siraeum, Mageeibacillus indolicus, and Saccharofermentans acetigenes, became predominant post-treatment in a significant proportion of mice, identified by 16S rRNA gene sequencing. Principal component analysis of QTOF MS of mouse feces successfully distinguished treated from non-treated mice as well as immunocompetent from neutropenic mice. We observe dynamic but distinct and reproducible responses in the mouse gut microbiota during and after TZP treatment and propose the current murine model as a useful tool for defining the more general post-antibiotic effects in the gastro-intestinal ecosystem where humanized antibiotic dosing may ultimately facilitate extrapolation to humans.

Differences in Gut Microbiome in Hospitalized Immunocompetent vs. Immunocompromised Children, Including Those With Sickle Cell Disease.

Background: Gut microbial diversity and composition play important roles in health. This cross-sectional study was designed to test the hypothesis that hospitalized children who may be relatively immunocompromised (IC), defined as those with cancer, sickle cell disease (SCD), transplantation, or receiving immunosuppressive therapy) would have decreased microbial diversity, increased Clostridioides difficile colonization and different species composition compared to non-immunocompromised (Non-IC) children admitted to the same pediatric unit. Methods: A stool sample was obtained within 72 h of admission to a single unit at The Children's Hospital at Montefiore, Bronx, NY from March 2016 to February 2017 and the microbiome assessed by 16S rRNA sequencing. C. difficile colonization was assessed by glutamate dehydrogenase antigen and toxin polymerase chain reaction assays. Results: Stool samples were obtained from 69 IC (32 SCD, 19 cancer, 9 transplantation and 9 other) and 37 Non-IC patients. There were no significant differences in microbial alpha diversity and C. difficile colonization comparing IC vs. non-IC patients. Lower alpha diversity, however, was independently associated with the use of proton pump inhibitors or antibiotics, including prophylactic penicillin in patients with SCD. Differences in specific species abundances were observed when comparing IC vs. non-IC patients, particularly children with SCD. Non-IC patients had increased abundance of commensals associated with health including Alistipes putredinis, Alistipes ihumii, Roseburia inulinivorans, Roseburia intestinalis, and Ruminococcus albus (p < 0.005). Conclusions: Antibiotics and proton pump inhibitors, which were more commonly used in IC children, were identified as risk factors for lower microbial diversity. Non-IC patients had higher abundance of several bacterial species associated with health. Longitudinal studies are needed to determine the clinical significance of these differences in gut microbiome.

Longitudinal Analysis of Serum Cytokine Levels and Gut Microbial Abundance Links IL-17/IL-22 With Clostridia and Insulin Sensitivity in Humans.

Recent studies using mouse models suggest that interaction between the gut microbiome and IL-17/IL-22-producing cells plays a role in the development of metabolic diseases. We investigated this relationship in humans using data from the prediabetes study of the Integrated Human Microbiome Project (iHMP). Specifically, we addressed the hypothesis that early in the onset of metabolic diseases there is a decline in serum levels of IL-17/IL-22, with concomitant changes in the gut microbiome. Clustering iHMP study participants on the basis of longitudinal IL-17/IL-22 profiles identified discrete groups. Individuals distinguished by low levels of IL-17/IL-22 were linked to established markers of metabolic disease, including insulin sensitivity. These individuals also displayed gut microbiome dysbiosis, characterized by decreased diversity, and IL-17/IL-22-related declines in the phylum Firmicutes, class Clostridia, and order Clostridiales This ancillary analysis of the iHMP data therefore supports a link between the gut microbiome, IL-17/IL-22, and the onset of metabolic diseases. This raises the possibility for novel, microbiome-related therapeutic targets that may effectively alleviate metabolic diseases in humans as they do in animal models.

Colon Cancer Prevention with Walnuts: A Longitudinal Study in Mice from the Perspective of a Gut Enterotype-like Cluster.

There is limited understanding of how walnut consumption inhibits the development of colorectal cancer. A possible mechanism may involve alterations to the gut microbiota. In this study, the effects of walnut on gut microbiota were tested in a mouse tumor bioassay using the colonotropic carcinogen, azoxymethane (AOM) added to the total Western diet (TWD). 16S rRNA pyrosequencing identified three enterotype-like clusters (E1, E2, and E3) in this murine model. E1, E2, and E3 are associated with AOM exposure, walnut consumption, and TWD diet, respectively. E2 and E3 showed distinct taxonomic and functional characteristics, while E1 represented an intermediate state. At the family level, E1 and E3 were both enriched with Bacteroidaceae, but driven by two different operational taxonomic units (OTU; OTU-2 for E1, OTU-4 for E3). E2 was overrepresented with Porphyromonadaceae and Lachnospiraceae, with OTU-3 (family Porphyromonadaceae) as the "driver" OTU for this cluster. Functionally, E3 is overrepresented with genes of glycan biosynthesis and metabolism, xenobiotic metabolism, and lipid metabolism. E2 is enriched with genes associated with cell motility, replication and repair, and amino acid metabolism. Longitudinally, E2 represents the gut microbial status of early life in these mice. In comparison with E1 and E3, E2 is associated with a moderate lower tumor burden (P = 0.12). Our results suggest that walnuts may reduce the risk of colorectal cancer within a Western diet by altering the gut microbiota. Our findings provide further evidence that colorectal cancer risk is potentially modifiable by diet via alterations to the microbiota.

New statistical method identifies cytokines that distinguish stool microbiomes.

Regressing an outcome or dependent variable onto a set of input or independent variables allows the analyst to measure associations between the two so that changes in the outcome can be described by and predicted by changes in the inputs. While there are many ways of doing this in classical statistics, where the dependent variable has certain properties (e.g., a scalar, survival time, count), little progress on regression where the dependent variable are microbiome taxa counts has been made that do not impose extremely strict conditions on the data. In this paper, we propose and apply a new regression model combining the Dirichlet-multinomial distribution with recursive partitioning providing a fully non-parametric regression model. This model, called DM-RPart, is applied to cytokine data and microbiome taxa count data and is applicable to any microbiome taxa count/metadata, is automatically fit, and intuitively interpretable. This is a model which can be applied to any microbiome or other compositional data and software (R package HMP) available through the R CRAN website.

Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis.

The 16S rRNA gene has been a mainstay of sequence-based bacterial analysis for decades. However, high-throughput sequencing of the full gene has only recently become a realistic prospect. Here, we use in silico and sequence-based experiments to critically re-evaluate the potential of the 16S gene to provide taxonomic resolution at species and strain level. We demonstrate that targeting of 16S variable regions with short-read sequencing platforms cannot achieve the taxonomic resolution afforded by sequencing the entire (~1500 bp) gene. We further demonstrate that full-length sequencing platforms are sufficiently accurate to resolve subtle nucleotide substitutions (but not insertions/deletions) that exist between intragenomic copies of the 16S gene. In consequence, we argue that modern analysis approaches must necessarily account for intragenomic variation between 16S gene copies. In particular, we demonstrate that appropriate treatment of full-length 16S intragenomic copy variants has the potential to provide taxonomic resolution of bacterial communities at species and strain level.

Microbiome Signatures Associated With Steatohepatitis and Moderate to Severe Fibrosis in Children With Nonalcoholic Fatty Liver Disease.

BACKGROUND & AIMS: The intestinal microbiome might affect the development and severity of nonalcoholic fatty liver disease (NAFLD). We analyzed microbiomes of children with and without NAFLD. METHODS: We performed a prospective, observational, cross-sectional study of 87 children (age range, 8-17 years) with biopsy-proven NAFLD and 37 children with obesity without NAFLD (controls). Fecal samples were collected and microbiome composition and functions were assessed using 16S ribosomal RNA amplicon sequencing and metagenomic shotgun sequencing. Microbial taxa were identified using zero-inflated negative binomial modeling. Genes contributing to bacterial pathways were identified using gene set enrichment analysis. RESULTS: Fecal microbiomes of children with NAFLD had lower α-diversity than those of control children (3.32 vs 3.52, P = .016). Fecal microbiomes from children with nonalcoholic steatohepatitis (NASH) had the lowest α-diversity (control, 3.52; NAFLD, 3.36; borderline NASH, 3.37; NASH, 2.97; P = .001). High abundance of Prevotella copri was associated with more severe fibrosis (P = .036). Genes for lipopolysaccharide biosynthesis were enriched in microbiomes from children with NASH (P < .001). Classification and regression tree model with level of alanine aminotransferase and relative abundance of the lipopolysaccharide pathway gene encoding 3-deoxy-d-manno-octulosonate 8-phosphate-phosphatase identified patients with NASH with an area under the receiver operating characteristic curve value of 0.92. Genes involved in flagellar assembly were enriched in the fecal microbiomes of patients with moderate to severe fibrosis (P < .001). Classification and regression tree models based on level of alanine aminotransferase and abundance of genes encoding flagellar biosynthesis protein had good accuracy for identifying case children with moderate to severe fibrosis (area under the receiver operating characteristic curve, 0.87). CONCLUSIONS: In an analysis of fecal microbiomes of children with NAFLD, we associated NAFLD and NASH with intestinal dysbiosis. NAFLD and its severity were associated with greater abundance of genes encoding inflammatory bacterial products. Alterations to the intestinal microbiome might contribute to the pathogenesis of NAFLD and be used as markers of disease or severity.

Longitudinal multi-omics of host-microbe dynamics in prediabetes.

Type 2 diabetes mellitus (T2D) is a growing health problem, but little is known about its early disease stages, its effects on biological processes or the transition to clinical T2D. To understand the earliest stages of T2D better, we obtained samples from 106 healthy individuals and individuals with prediabetes over approximately four years and performed deep profiling of transcriptomes, metabolomes, cytokines, and proteomes, as well as changes in the microbiome. This rich longitudinal data set revealed many insights: first, healthy profiles are distinct among individuals while displaying diverse patterns of intra- and/or inter-personal variability. Second, extensive host and microbial changes occur during respiratory viral infections and immunization, and immunization triggers potentially protective responses that are distinct from responses to respiratory viral infections. Moreover, during respiratory viral infections, insulin-resistant participants respond differently than insulin-sensitive participants. Third, global co-association analyses among the thousands of profiled molecules reveal specific host-microbe interactions that differ between insulin-resistant and insulin-sensitive individuals. Last, we identified early personal molecular signatures in one individual that preceded the onset of T2D, including the inflammation markers interleukin-1 receptor agonist (IL-1RA) and high-sensitivity C-reactive protein (CRP) paired with xenobiotic-induced immune signalling. Our study reveals insights into pathways and responses that differ between glucose-dysregulated and healthy individuals during health and disease and provides an open-access data resource to enable further research into healthy, prediabetic and T2D states.

A longitudinal big data approach for precision health.

Precision health relies on the ability to assess disease risk at an individual level, detect early preclinical conditions and initiate preventive strategies. Recent technological advances in omics and wearable monitoring enable deep molecular and physiological profiling and may provide important tools for precision health. We explored the ability of deep longitudinal profiling to make health-related discoveries, identify clinically relevant molecular pathways and affect behavior in a prospective longitudinal cohort (n = 109) enriched for risk of type 2 diabetes mellitus. The cohort underwent integrative personalized omics profiling from samples collected quarterly for up to 8 years (median, 2.8 years) using clinical measures and emerging technologies including genome, immunome, transcriptome, proteome, metabolome, microbiome and wearable monitoring. We discovered more than 67 clinically actionable health discoveries and identified multiple molecular pathways associated with metabolic, cardiovascular and oncologic pathophysiology. We developed prediction models for insulin resistance by using omics measurements, illustrating their potential to replace burdensome tests. Finally, study participation led the majority of participants to implement diet and exercise changes. Altogether, we conclude that deep longitudinal profiling can lead to actionable health discoveries and provide relevant information for precision health.

Symbiotic organs shaped by distinct modes of genome evolution in cephalopods.

Microbes have been critical drivers of evolutionary innovation in animals. To understand the processes that influence the origin of specialized symbiotic organs, we report the sequencing and analysis of the genome of Euprymna scolopes, a model cephalopod with richly characterized host-microbe interactions. We identified large-scale genomic reorganization shared between E. scolopes and Octopus bimaculoides and posit that this reorganization has contributed to the evolution of cephalopod complexity. To reveal genomic signatures of host-symbiont interactions, we focused on two specialized organs of E. scolopes: the light organ, which harbors a monoculture of Vibrio fischeri, and the accessory nidamental gland (ANG), a reproductive organ containing a bacterial consortium. Our findings suggest that the two symbiotic organs within E. scolopes originated by different evolutionary mechanisms. Transcripts expressed in these microbe-associated tissues displayed their own unique signatures in both coding sequences and the surrounding regulatory regions. Compared with other tissues, the light organ showed an abundance of genes associated with immunity and mediating light, whereas the ANG was enriched in orphan genes known only from E. scolopes Together, these analyses provide evidence for different patterns of genomic evolution of symbiotic organs within a single host.

Diet during Pregnancy and Infancy and the Infant Intestinal Microbiome.

OBJECTIVES: To determine the association between diet during pregnancy and infancy, including breastfeeding vs formula feeding, solid food introduction, and the infant intestinal microbiome. STUDY DESIGN: Infants participating in the Vitamin D Antenatal Asthma Reduction Trial were included in this study (n = 323). Maternal and infant diets were assessed by questionnaire. Infant stool samples were collected at age 3-6 months. Stool sequencing was performed using the Roche 454 platform. Analyses were stratified by race/ethnicity. RESULTS: Breastfeeding, compared with formula feeding, was independently associated with infant intestinal microbial diversity. Breastfeeding also had the most consistent associations with individual taxa that have been previously linked to early-life diet and health outcomes (eg, Bifidobacterium). Maternal diet during pregnancy and solid food introduction were less associated with the infant gut microbiome than breastfeeding status. We found evidence of a possible interaction between breastfeeding and child race/ethnicity on microbial composition. CONCLUSIONS: Breastfeeding vs formula feeding is the dietary factor that is most consistently independently associated with the infant intestinal microbiome. The relationship between breastfeeding status and intestinal microbiome composition varies by child race/ethnicity. Future studies will need to investigate factors, including genomic factors, which may influence the response of the microbiome to diet. TRIAL REGISTRATION: ClinicalTrials.gov: NCT00920621.

Extensive Gene Amplification as a Mechanism for Piperacillin-Tazobactam Resistance in Escherichia coli.

Although the TEM-1 ß-lactamase (BlaTEM-1) hydrolyzes penicillins and narrow-spectrum cephalosporins, organisms expressing this enzyme are typically susceptible to ß-lactam/ß-lactamase inhibitor combinations such as piperacillin-tazobactam (TZP). However, our previous work led to the discovery of 28 clinical isolates of Escherichia coli resistant to TZP that contained only blaTEM-1 One of these isolates, E. coli 907355, was investigated further in this study. E. coli 907355 exhibited significantly higher ß-lactamase activity and BlaTEM-1 protein levels when grown in the presence of subinhibitory concentrations of TZP. A corresponding TZP-dependent increase in blaTEM-1 copy number was also observed, with as many as 113 copies of the gene detected per cell. These results suggest that TZP treatment promotes an increase in blaTEM-1 gene dosage, allowing BlaTEM-1 to reach high enough levels to overcome inactivation by the available tazobactam in the culture. To better understand the nature of the blaTEM-1 copy number proliferation, whole-genome sequence (WGS) analysis was performed on E. coli 907355 in the absence and presence of TZP. The WGS data revealed that the blaTEM-1 gene is located in a 10-kb genomic resistance module (GRM) that contains multiple resistance genes and mobile genetic elements. The GRM was found to be tandemly repeated at least 5 times within a p1ESCUM/p1ECUMN-like plasmid when bacteria were grown in the presence of TZP.IMPORTANCE Understanding how bacteria acquire resistance to antibiotics is essential for treating infected patients effectively, as well as preventing the spread of resistant organisms. In this study, a clinical isolate of E. coli was identified that dedicated more than 15% of its genome toward tandem amplification of a ~10-kb resistance module, allowing it to escape antibiotic-mediated killing. Our research is significant in that it provides one possible explanation for clinical isolates that exhibit discordant behavior when tested for antibiotic resistance by different phenotypic methods. Our research also shows that GRM amplification is difficult to detect by short-read WGS technologies. Analysis of raw long-read sequence data was required to confirm GRM amplification as a mechanism of antibiotic resistance.

Integrative Personal Omics Profiles during Periods of Weight Gain and Loss.

Advances in omics technologies now allow an unprecedented level of phenotyping for human diseases, including obesity, in which individual responses to excess weight are heterogeneous and unpredictable. To aid the development of better understanding of these phenotypes, we performed a controlled longitudinal weight perturbation study combining multiple omics strategies (genomics, transcriptomics, multiple proteomics assays, metabolomics, and microbiomics) during periods of weight gain and loss in humans. Results demonstrated that: (1) weight gain is associated with the activation of strong inflammatory and hypertrophic cardiomyopathy signatures in blood; (2) although weight loss reverses some changes, a number of signatures persist, indicative of long-term physiologic changes; (3) we observed omics signatures associated with insulin resistance that may serve as novel diagnostics; (4) specific biomolecules were highly individualized and stable in response to perturbations, potentially representing stable personalized markers. Most data are available open access and serve as a valuable resource for the community.

A prospective microbiome-wide association study of food sensitization and food allergy in early childhood.

BACKGROUND: Alterations in the intestinal microbiome are prospectively associated with the development of asthma; less is known regarding the role of microbiome alterations in food allergy development. METHODS: Intestinal microbiome samples were collected at age 3-6 months in children participating in the follow-up phase of an interventional trial of high-dose vitamin D given during pregnancy. At age 3, sensitization to foods (milk, egg, peanut, soy, wheat, walnut) was assessed. Food allergy was defined as caretaker report of healthcare provider-diagnosed allergy to the above foods prior to age 3 with evidence of IgE sensitization. Analysis was performed using Phyloseq and DESeq2; P-values were adjusted for multiple comparisons. RESULTS: Complete data were available for 225 children; there were 87 cases of food sensitization and 14 cases of food allergy. Microbial diversity measures did not differ between food sensitization and food allergy cases and controls. The genera Haemophilus (log2 fold change -2.15, P=.003), Dialister (log2 fold change -2.22, P=.009), Dorea (log2 fold change -1.65, P=.02), and Clostridium (log2 fold change -1.47, P=.002) were underrepresented among subjects with food sensitization. The genera Citrobacter (log2 fold change -3.41, P=.03), Oscillospira (log2 fold change -2.80, P=.03), Lactococcus (log2 fold change -3.19, P=.05), and Dorea (log2 fold change -3.00, P=.05) were underrepresented among subjects with food allergy. CONCLUSIONS: The temporal association between bacterial colonization and food sensitization and allergy suggests that the microbiome may have a causal role in the development of food allergy. Our findings have therapeutic implications for the prevention and treatment of food allergy.