RESUMEN
Negative psychological states impact immunity by altering the gut microbiome. However, the relationship between brain states and microbiome composition remains unclear. We show that Brunner's glands in the duodenum couple stress-sensitive brain circuits to bacterial homeostasis. Brunner's glands mediated the enrichment of gut Lactobacillus species in response to vagus nerve stimulation. Cell-specific ablation of the glands markedly suppressed Lactobacilli counts and heightened vulnerability to infection. In the forebrain, we mapped a vagally mediated, polysynaptic circuit connecting the central nucleus of the amygdala to Brunner's glands. Chronic stress suppressed central amygdala activity and phenocopied the effects of gland lesions. Conversely, excitation of either the central amygdala or parasympathetic vagal neurons activated Brunner's glands and reversed the effects of stress on the gut microbiome and immunity. The findings revealed a tractable brain-body mechanism linking psychological states to host defense.
RESUMEN
Clostridia are abundant in the human gut and comprise families associated with host health such as Oscillospiraceae, which has been correlated with leanness. However, culturing bacteria within this family is challenging, leading to their detection primarily through 16S rRNA amplicon sequencing, which has a limited ability to unravel diversity at low taxonomic levels, or by shotgun metagenomics, which is hindered by its high costs and complexity. In this cross-sectional study involving 114 Colombian adults, we used an amplicon-based sequencing strategy with alternative markers-gyrase subunit B (gyrB) and DNA K chaperone heat protein 70 (dnaK)-that evolve faster than the 16S rRNA gene. Comparing the diversity and abundance observed with the three markers in our cohort, we found a reduction in the diversity of Clostridia, particularly within Lachnospiraceae and Oscillospiraceae among obese individuals [as measured by the body mass index (BMI)]. Within Lachnospiraceae, the diversity of Ruminococcus_A negatively correlated with BMI. Within Oscillospiraceae, the genera CAG-170 and Vescimonas also exhibited this negative correlation. In addition, the abundance of Vescimonas was negatively correlated with BMI. Leveraging shotgun metagenomic data, we conducted a phylogenetic and genomic characterization of 120 metagenome-assembled genomes from Vescimonas obtained from a larger sample of the same cohort. We identified 17 of the 72 reported species. The functional annotation of these genomes showed the presence of multiple carbohydrate-active enzymes, particularly glycosyl transferases and glycoside hydrolases, suggesting potential beneficial roles in fiber degradation, carbohydrate metabolism, and butyrate production. IMPORTANCE: The gut microbiota is diverse across various taxonomic levels. At the intra-species level, it comprises multiple strains, some of which may be host-specific. However, our understanding of fine-grained diversity has been hindered by the use of the conserved 16S rRNA gene. While shotgun metagenomics offers higher resolution, it remains costly, may fail to identify specific microbes in complex samples, and requires extensive computational resources and expertise. To address this, we employed a simple and cost-effective analysis of alternative genetic markers to explore diversity within Clostridia, a crucial group within the human gut microbiota whose diversity may be underestimated. We found high intra-species diversity for certain groups and associations with obesity. Notably, we identified Vescimonas, an understudied group. Making use of metagenomic data, we inferred functionality, uncovering potential beneficial roles in dietary fiber and carbohydrate degradation, as well as in short-chain fatty acid production.
RESUMEN
Microbial species diversify into strains through single-nucleotide mutations and structural changes, such as recombination, insertions and deletions. Most strain-comparison methods quantify differences in single-nucleotide polymorphisms (SNPs) and are insensitive to structural changes. However, recombination is an important driver of phenotypic diversification in many species, including human pathogens. We introduce SynTracker, a tool that compares microbial strains using genome synteny-the order of sequence blocks in homologous genomic regions-in pairs of metagenomic assemblies or genomes. Genome synteny is a rich source of genomic information untapped by current strain-comparison tools. SynTracker has low sensitivity to SNPs, has no database requirement and is robust to sequencing errors. It outperforms existing tools when tracking strains in metagenomic data and is particularly suited for phages, plasmids and other low-data contexts. Applied to single-species datasets and human gut metagenomes, SynTracker, combined with an SNP-based tool, detects strains enriched in either point mutations or structural changes, providing insights into microbial evolution in situ.
RESUMEN
The extent to which bacterial lipids produced by the gut microbiota penetrate host tissues is unclear. Here, we combined mass spectrometry approaches to identify lipids produced by the human gut symbiont Bacteroides thetaiotaomicron (B. theta) and spatially track these bacterial lipids in the mouse colon. We characterize 130 B. theta lipids by liquid chromatography-tandem mass spectrometry (LC-MS/MS), using wild-type and mutant B. theta strains to confidently identify lipid structures and their interconnected pathways in vitro. Of these, 103 B. theta lipids can be detected and spatially mapped in a single MALDI mass spectrometry imaging run. We map unlabeled bacterial lipids across colon sections of germ-free and specific-pathogen-free (SPF) mice and mice mono-colonized with wild-type or sphingolipid-deficient (BTMUT) B. theta. We observe co-localization of bacterially derived phosphatidic acid with host tissues in BTMUT mice, consistent with lipid penetration into host tissues. These results indicate limited and selective transfer of bacterial lipids to the host.
Asunto(s)
Bacteroides thetaiotaomicron , Colon , Microbioma Gastrointestinal , Lipidómica , Animales , Ratones , Bacteroides thetaiotaomicron/metabolismo , Microbioma Gastrointestinal/fisiología , Colon/microbiología , Colon/metabolismo , Lípidos/análisis , Espectrometría de Masas en Tándem , Cromatografía Liquida , Metabolismo de los Lípidos , Vida Libre de Gérmenes , Organismos Libres de Patógenos Específicos , Ácidos Fosfatidicos/metabolismo , Humanos , Espectrometría de Masa por Láser de Matriz Asistida de Ionización Desorción , Esfingolípidos/metabolismo , Ratones Endogámicos C57BL , FemeninoRESUMEN
The incidence of antibiotic-resistant bacterial infections is increasing, and development of new antibiotics has been deprioritised by the pharmaceutical industry. Interdisciplinary research approaches, based on the ecological principles of bacterial fitness, competition, and transmission, could open new avenues to combat antibiotic-resistant infections. Many facultative bacterial pathogens use human mucosal surfaces as their major reservoirs and induce infectious diseases to aid their lateral transmission to new host organisms under some pathological states of the microbiome and host. Beneficial bacterial commensals can outcompete specific pathogens, thereby lowering the capacity of the pathogens to spread and cause serious infections. Despite the clinical relevance, however, the understanding of commensal-pathogen interactions in their natural habitats remains poor. In this Personal View, we highlight directions to intensify research on the interactions between bacterial pathogens and commensals in the context of human microbiomes and host biology that can lead to the development of innovative and sustainable ways of preventing and treating infectious diseases.
Asunto(s)
Bacterias , Infecciones Bacterianas , Microbiota , Simbiosis , Humanos , Infecciones Bacterianas/microbiología , Infecciones Bacterianas/transmisión , Bacterias/patogenicidad , Bacterias/efectos de los fármacos , Animales , Interacciones Huésped-Patógeno , Antibacterianos/farmacología , Antibacterianos/uso terapéuticoRESUMEN
This study investigated the relationship between gut microbiota and neuropsychiatric disorders (NPDs), specifically anxiety disorder (ANXD) and/or major depressive disorder (MDD), as defined by DSM-IV or V criteria. The study also examined the influence of medication use, particularly antidepressants and/or anxiolytics, classified through the Anatomical Therapeutic Chemical (ATC) Classification System, on the gut microbiota. Both 16S rRNA gene amplicon sequencing and shallow shotgun sequencing were performed on DNA extracted from 666 fecal samples from the Tulsa-1000 and NeuroMAP CoBRE cohorts. The results highlight the significant influence of medication use; antidepressant use is associated with significant differences in gut microbiota beta diversity and has a larger effect size than NPD diagnosis. Next, specific microbes were associated with ANXD and MDD, highlighting their potential for non-pharmacological intervention. Finally, the study demonstrated the capability of Random Forest classifiers to predict diagnoses of NPD and medication use from microbial profiles, suggesting a promising direction for the use of gut microbiota as biomarkers for NPD. The findings suggest that future research on the gut microbiota's role in NPD and its interactions with pharmacological treatments are needed.
RESUMEN
BACKGROUND: Standard pediatric growth curves cannot be used to impute missing height or weight measurements in individual children. The Michaelis-Menten equation, used for characterizing substrate-enzyme saturation curves, has been shown to model growth in many organisms including nonhuman vertebrates. We investigated whether this equation could be used to interpolate missing growth data in children in the first three years of life and compared this interpolation to several common interpolation methods and pediatric growth models. METHODS: We developed a modified Michaelis-Menten equation and compared expected to actual growth, first in a local birth cohort (N = 97) then in a large, outpatient, pediatric sample (N = 14,695). RESULTS: The modified Michaelis-Menten equation showed excellent fit for both infant weight (median RMSE: boys: 0.22 kg [IQR:0.19; 90% < 0.43]; girls: 0.20 kg [IQR:0.17; 90% < 0.39]) and height (median RMSE: boys: 0.93 cm [IQR:0.53; 90% < 1.0]; girls: 0.91 cm [IQR:0.50;90% < 1.0]). Growth data were modeled accurately with as few as four values from routine well-baby visits in year 1 and seven values in years 1-3; birth weight or length was essential for best fit. Interpolation with this equation had comparable (for weight) or lower (for height) mean RMSE compared to the best performing alternative models. CONCLUSIONS: A modified Michaelis-Menten equation accurately describes growth in healthy babies aged 0-36 months, allowing interpolation of missing weight and height values in individual longitudinal measurement series. The growth pattern in healthy babies in resource-rich environments mirrors an enzymatic saturation curve.
Asunto(s)
Cinética , Masculino , Lactante , Femenino , Humanos , Niño , Peso al NacerAsunto(s)
Microbioma Gastrointestinal , Microbiota , Humanos , Carboximetilcelulosa de Sodio , Intestinos , DietaRESUMEN
The gut bacteria of the family Christensenellaceae are consistently associated with metabolic health, but their role in promoting host health is not fully understood. Here, we explored the effect of Christensenella minuta amendment on voluntary physical activity and the gut microbiome. We inoculated male and female germ-free mice with an obese human donor microbiota together with live or heat-killed C. minuta for 28 days and measured physical activity in respirometry cages. Compared to heat-killed, the live-C. minuta treatment resulted in reduced feed efficiency and higher levels of physical activity, with significantly greater distance traveled for males and higher levels of small movements and resting metabolic rate in females. Sex-specific effects of C. minuta treatment may be in part attributable to different housing conditions for males and females. Amendment with live C. minuta boosted gut microbial biomass in both sexes, immobilizing dietary carbon in the microbiome, and mice with high levels of C. minuta lose more energy in stool. Live C. minuta also reduced within and between-host gut microbial diversity. Overall, our results showed that C. minuta acts as a keystone species: despite low relative abundance, it has a large impact on its ecosystem, from the microbiome to host energy homeostasis.IMPORTANCEThe composition of the human gut microbiome is associated with human health. Within the human gut microbiome, the relative abundance of the bacterial family Christensenellaceae has been shown to correlate with metabolic health and a lean body type. The mechanisms underpinning this effect remain unclear. Here, we show that live C. minuta influences host physical activity and metabolic energy expenditure, accompanied by changes in murine metabolism and the gut microbial community in a sex-dependent manner in comparison to heat-killed C. minuta. Importantly, live C. minuta boosts the biomass of the microbiome in the gut, and a higher level of C. minuta is associated with greater loss of energy in stool. These observations indicate that modulation of activity levels and changes to the microbiome are ways in which the Christensenellaceae can influence host energy homeostasis and health.
Asunto(s)
Clostridiales , Microbioma Gastrointestinal , Microbiota , Humanos , Masculino , Femenino , Animales , Ratones , Biomasa , Heces/microbiología , Bacterias/metabolismoRESUMEN
Methanogenesis allows methanogenic archaea to generate cellular energy for their growth while producing methane. Thermophilic hydrogenotrophic species of the genus Methanothermobacter have been recognized as robust biocatalysts for a circular carbon economy and are already applied in power-to-gas technology with biomethanation, which is a platform to store renewable energy and utilize captured carbon dioxide. Here, we generated curated genome-scale metabolic reconstructions for three Methanothermobacter strains and investigated differences in the growth performance of these same strains in chemostat bioreactor experiments with hydrogen and carbon dioxide or formate as substrates. Using an integrated systems biology approach, we identified differences in formate anabolism between the strains and revealed that formate anabolism influences the diversion of carbon between biomass and methane. This finding, together with the omics datasets and the metabolic models we generated, can be implemented for biotechnological applications of Methanothermobacter in power-to-gas technology, and as a perspective, for value-added chemical production.
RESUMEN
Obesity (OB) and cardiometabolic disease are major public health issues linked to changes in the gut microbiome. OB and poor cardiometabolic health status (CHS) are often comorbid, which hinders efforts to identify components of the microbiome uniquely linked to either one. Here, we used a deeply phenotyped cohort of 408 adults from Colombia, including subjects with OB, unhealthy CHS, or both, to validate previously reported features of gut microbiome function and diversity independently correlated with OB or CHS using fecal metagenomes. OB was defined by body mass index, waist circumference, and body fat; CHS as healthy or unhealthy according to blood biochemistry and anthropometric data. We found that OB, more so than metabolic status, drove associations with gut microbiome structure and functions. The microbiome of obese individuals with and without co-existing unhealthy CHS was characterized by reduced metagenomic diversity, reduced fermentative potential and elevated capacity to respond to oxidative stress and produce bacterial antigens. Disease-linked features were correlated with increased host blood pressure and inflammatory markers, and were mainly contributed by members of the family Enterobacteriaceae. Our results link OB with a microbiome able to tolerate an inflammatory and oxygenated gut state, and suggest that OB is the main driver of microbiome functional differences when poor CHS is a comorbidity.
Asunto(s)
Microbioma Gastrointestinal , Microbiota , Adulto , Humanos , Obesidad , Tejido Adiposo , AntropometríaRESUMEN
Genetic manipulation is necessary to interrogate the functions of microbes in their environments, such as the human gut microbiome. Yet, the vast majority of human gut microbiome species are not genetically tractable. Here, we review the hurdles to seizing genetic control of more species. We address the barriers preventing the application of genetic techniques to gut microbes and report on genetic systems currently under development. While methods aimed at genetically transforming many species simultaneously in situ show promise, they are unable to overcome many of the same challenges that exist for individual microbes. Unless a major conceptual breakthrough emerges, the genetic tractability of the microbiome will remain an arduous task. Increasing the list of genetically tractable organisms from the human gut remains one of the highest priorities for microbiome research and will provide the foundation for microbiome engineering.
Asunto(s)
Microbioma Gastrointestinal , Microbiota , HumanosRESUMEN
The number of published metagenome assemblies is rapidly growing due to advances in sequencing technologies. However, sequencing errors, variable coverage, repetitive genomic regions, and other factors can produce misassemblies, which are challenging to detect for taxonomically novel genomic data. Assembly errors can affect all downstream analyses of the assemblies. Accuracy for the state of the art in reference-free misassembly prediction does not exceed an AUPRC of 0.57, and it is not clear how well these models generalize to real-world data. Here, we present the Residual neural network for Misassembled Contig identification (ResMiCo), a deep learning approach for reference-free identification of misassembled contigs. To develop ResMiCo, we first generated a training dataset of unprecedented size and complexity that can be used for further benchmarking and developments in the field. Through rigorous validation, we show that ResMiCo is substantially more accurate than the state of the art, and the model is robust to novel taxonomic diversity and varying assembly methods. ResMiCo estimated 7% misassembled contigs per metagenome across multiple real-world datasets. We demonstrate how ResMiCo can be used to optimize metagenome assembly hyperparameters to improve accuracy, instead of optimizing solely for contiguity. The accuracy, robustness, and ease-of-use of ResMiCo make the tool suitable for general quality control of metagenome assemblies and assembly methodology optimization.
Asunto(s)
Aprendizaje Profundo , Metagenoma , Metagenoma/genética , Genómica/métodos , Análisis de Secuencia de ADN/métodos , Metagenómica , Programas InformáticosRESUMEN
BACKGROUND & AIMS: Fiber-rich foods promote health, but mechanisms by which they do so remain poorly defined. Screening fiber types, in mice, revealed psyllium had unique ability to ameliorate 2 chronic inflammatory states, namely, metabolic syndrome and colitis. We sought to determine the mechanism of action of the latter. METHODS: Mice were fed grain-based chow, which is naturally rich in fiber or compositionally defined diets enriched with semi-purified fibers. Mice were studied basally and in models of chemical-induced and T-cell transfer colitis. RESULTS: Relative to all diets tested, mice consuming psyllium-enriched compositionally defined diets were markedly protected against both dextran sulfate sodium- and T-cell transfer-induced colitis, as revealed by clinical-type, histopathologic, morphologic, and immunologic parameters. Such protection associated with stark basal changes in the gut microbiome but was independent of fermentation and, moreover, maintained in mice harboring a minimal microbiota (ie, Altered Schaedler Flora). Transcriptomic analysis revealed psyllium induced expression of genes mediating bile acids (BA) secretion, suggesting that psyllium's known ability to bind BA might contribute to its ability to prevent colitis. As expected, psyllium resulted in elevated level of fecal BA, reflecting their removal from enterohepatic circulation but, in stark contrast to the BA sequestrant cholestyramine, increased serum BA levels. Moreover, the use of BA mimetics that activate the farnesoid X receptor (FXR), as well as the use of FXR-knockout mice, suggested that activation of FXR plays a central role in psyllium's protection against colitis. CONCLUSIONS: Psyllium protects against colitis via altering BA metabolism resulting in activation of FXR, which suppresses pro-inflammatory signaling.
Asunto(s)
Colitis , Psyllium , Ratones , Animales , Psyllium/efectos adversos , Ácidos y Sales Biliares , Promoción de la Salud , Colitis/inducido químicamente , Colitis/prevención & control , Colitis/metabolismo , Inflamación , Ratones NoqueadosRESUMEN
The human gut virome and its early life development are poorly understood. Prior studies have captured single-point assessments with the evolution of the infant virome remaining largely unexplored. We performed viral metagenomic sequencing on stool samples collected longitudinally from a cohort of 53 infants from age 2 weeks to 3 years (80.7 billion reads), and from their mothers (9.8 billion reads) to examine and compare viromes. The asymptomatic infant virome consisted of bacteriophages, nonhuman dietary/environmental viruses, and human-host viruses, predominantly picornaviruses. In contrast, human-host viruses were largely absent from the maternal virome. Previously undescribed, sequence-divergent vertebrate viruses were detected in the maternal but not infant virome. As infants aged, the phage component evolved to resemble the maternal virome, but by age 3, the human-host component remained dissimilar from the maternal virome. Thus, early life virome development is determined predominantly by dietary, infectious, and environmental factors rather than direct maternal acquisition.
Asunto(s)
Bacteriófagos , Virus , Femenino , Humanos , Viroma/genética , Virus/genética , Bacteriófagos/genética , Madres , Metagenoma , MetagenómicaRESUMEN
Flagellin, the protein subunit of the bacterial flagellum, stimulates the innate immune receptor Toll-like receptor 5 (TLR5) after pattern recognition or evades TLR5 through lack of recognition. This binary response fails to explain the weak agonism of flagellins from commensal bacteria, raising the question of how TLR5 response is tuned. Here, we screened abundant flagellins present in metagenomes from human gut for both TLR5 recognition and activation and uncovered a class of flagellin-TLR5 interaction termed silent recognition. Silent flagellins were weak TLR5 agonists despite pattern recognition. Receptor activity was tuned by a TLR5-flagellin interaction distal to the site of pattern recognition that was present in Salmonella flagellin but absent in silent flagellins. This interaction enabled flagellin binding to preformed TLR5 dimers and increased TLR5 signaling by several orders of magnitude. Silent recognition by TLR5 occurred in human organoids and mice, and silent flagellin proteins were present in human stool. These flagellins were produced primarily by the abundant gut bacteria Lachnospiraceae and were enriched in nonindustrialized populations. Our findings provide a mechanism for the innate immune system to tolerate commensal-derived flagellins while remaining vigilant to the presence of flagellins produced by pathogens.
Asunto(s)
Flagelina , Receptor Toll-Like 5 , Animales , Humanos , Ratones , Bacterias , Flagelina/metabolismo , Transducción de Señal , IntestinosRESUMEN
Whether the human fetus and the prenatal intrauterine environment (amniotic fluid and placenta) are stably colonized by microbial communities in a healthy pregnancy remains a subject of debate. Here we evaluate recent studies that characterized microbial populations in human fetuses from the perspectives of reproductive biology, microbial ecology, bioinformatics, immunology, clinical microbiology and gnotobiology, and assess possible mechanisms by which the fetus might interact with microorganisms. Our analysis indicates that the detected microbial signals are likely the result of contamination during the clinical procedures to obtain fetal samples or during DNA extraction and DNA sequencing. Furthermore, the existence of live and replicating microbial populations in healthy fetal tissues is not compatible with fundamental concepts of immunology, clinical microbiology and the derivation of germ-free mammals. These conclusions are important to our understanding of human immune development and illustrate common pitfalls in the microbial analyses of many other low-biomass environments. The pursuit of a fetal microbiome serves as a cautionary example of the challenges of sequence-based microbiome studies when biomass is low or absent, and emphasizes the need for a trans-disciplinary approach that goes beyond contamination controls by also incorporating biological, ecological and mechanistic concepts.
Asunto(s)
Biomasa , Contaminación de ADN , Feto , Microbiota , Animales , Femenino , Humanos , Embarazo , Líquido Amniótico/inmunología , Líquido Amniótico/microbiología , Mamíferos , Microbiota/genética , Placenta/inmunología , Placenta/microbiología , Feto/inmunología , Feto/microbiología , Reproducibilidad de los ResultadosRESUMEN
OBJECTIVE: Gestational diabetes mellitus (GDM) is a condition in which women without diabetes are diagnosed with glucose intolerance during pregnancy, typically in the second or third trimester. Early diagnosis, along with a better understanding of its pathophysiology during the first trimester of pregnancy, may be effective in reducing incidence and associated short-term and long-term morbidities. DESIGN: We comprehensively profiled the gut microbiome, metabolome, inflammatory cytokines, nutrition and clinical records of 394 women during the first trimester of pregnancy, before GDM diagnosis. We then built a model that can predict GDM onset weeks before it is typically diagnosed. Further, we demonstrated the role of the microbiome in disease using faecal microbiota transplant (FMT) of first trimester samples from pregnant women across three unique cohorts. RESULTS: We found elevated levels of proinflammatory cytokines in women who later developed GDM, decreased faecal short-chain fatty acids and altered microbiome. We next confirmed that differences in GDM-associated microbial composition during the first trimester drove inflammation and insulin resistance more than 10 weeks prior to GDM diagnosis using FMT experiments. Following these observations, we used a machine learning approach to predict GDM based on first trimester clinical, microbial and inflammatory markers with high accuracy. CONCLUSION: GDM onset can be identified in the first trimester of pregnancy, earlier than currently accepted. Furthermore, the gut microbiome appears to play a role in inflammation-induced GDM pathogenesis, with interleukin-6 as a potential contributor to pathogenesis. Potential GDM markers, including microbiota, can serve as targets for early diagnostics and therapeutic intervention leading to prevention.
Asunto(s)
Diabetes Gestacional , Microbiota , Embarazo , Femenino , Humanos , Diabetes Gestacional/diagnóstico , Tercer Trimestre del Embarazo , Inflamación , CitocinasRESUMEN
Background and Objectives: Standard pediatric growth curves cannot be used to impute missing height or weight measurements in individual children. The Michaelis-Menten equation, used for characterizing substrate-enzyme saturation curves, has been shown to model growth in many organisms including nonhuman vertebrates. We investigated this equation could be used to interpolate missing growth data in children in the first three years of life. Methods: We developed a modified Michaelis-Menten equation and compared expected to actual growth, first in a local birth cohort (N=97) then in a large, outpatient, pediatric sample (N=14,695). Results: The modified Michaelis-Menten equation showed excellent fit for both infant weight (median RMSE: boys: 0.22kg [IQR:0.19; 90%<0.43]; girls: 0.20kg [IQR:0.17; 90%<0.39]) and height (median RMSE: boys: 0.93cm [IQR:0.53; 90%<1.0]; girls: 0.91cm [IQR:0.50;90%<1.0]). Growth data were modeled accurately with as few as four values from routine well-baby visits in year 1 and seven values in years 1-3; birth weight or length was essential for best fit. Conclusions: A modified Michaelis-Menten equation accurately describes growth in healthy babies aged 0-36 months, allowing interpolation of missing weight and height values in individual longitudinal measurement series. The growth pattern in healthy babies in resource-rich environments mirrors an enzymatic saturation curve.
RESUMEN
Tree ensemble machine learning models are increasingly used in microbiome science as they are compatible with the compositional, high-dimensional, and sparse structure of sequence-based microbiome data. While such models are often good at predicting phenotypes based on microbiome data, they only yield limited insights into how microbial taxa may be associated. We developed endoR, a method to interpret tree ensemble models. First, endoR simplifies the fitted model into a decision ensemble. Then, it extracts information on the importance of individual features and their pairwise interactions, displaying them as an interpretable network. Both the endoR network and importance scores provide insights into how features, and interactions between them, contribute to the predictive performance of the fitted model. Adjustable regularization and bootstrapping help reduce the complexity and ensure that only essential parts of the model are retained. We assessed endoR on both simulated and real metagenomic data. We found endoR to have comparable accuracy to other common approaches while easing and enhancing model interpretation. Using endoR, we also confirmed published results on gut microbiome differences between cirrhotic and healthy individuals. Finally, we utilized endoR to explore associations between human gut methanogens and microbiome components. Indeed, these hydrogen consumers are expected to interact with fermenting bacteria in a complex syntrophic network. Specifically, we analyzed a global metagenome dataset of 2203 individuals and confirmed the previously reported association between Methanobacteriaceae and Christensenellales. Additionally, we observed that Methanobacteriaceae are associated with a network of hydrogen-producing bacteria. Our method accurately captures how tree ensembles use features and interactions between them to predict a response. As demonstrated by our applications, the resultant visualizations and summary outputs facilitate model interpretation and enable the generation of novel hypotheses about complex systems.