Host-Microbe-Drug-Nutrient Screen Identifies Bacterial Effectors of Metformin Therapy.

Metformin is the first-line therapy for treating type 2 diabetes and a promising anti-aging drug. We set out to address the fundamental question of how gut microbes and nutrition, key regulators of host physiology, affect the effects of metformin. Combining two tractable genetic models, the bacterium E. coli and the nematode C. elegans, we developed a high-throughput four-way screen to define the underlying host-microbe-drug-nutrient interactions. We show that microbes integrate cues from metformin and the diet through the phosphotransferase signaling pathway that converges on the transcriptional regulator Crp. A detailed experimental characterization of metformin effects downstream of Crp in combination with metabolic modeling of the microbiota in metformin-treated type 2 diabetic patients predicts the production of microbial agmatine, a regulator of metformin effects on host lipid metabolism and lifespan. Our high-throughput screening platform paves the way for identifying exploitable drug-nutrient-microbiome interactions to improve host health and longevity through targeted microbiome therapies. VIDEO ABSTRACT.

Dynamic optimization reveals alveolar epithelial cells as key mediators of host defense in invasive aspergillosis.

Aspergillus fumigatus is an important human fungal pathogen and its conidia are constantly inhaled by humans. In immunocompromised individuals, conidia can grow out as hyphae that damage lung epithelium. The resulting invasive aspergillosis is associated with devastating mortality rates. Since infection is a race between the innate immune system and the outgrowth of A. fumigatus conidia, we use dynamic optimization to obtain insight into the recruitment and depletion of alveolar macrophages and neutrophils. Using this model, we obtain key insights into major determinants of infection outcome on host and pathogen side. On the pathogen side, we predict in silico and confirm in vitro that germination speed is an important virulence trait of fungal pathogens due to the vulnerability of conidia against host defense. On the host side, we found that epithelial cells, which have been underappreciated, play a role in fungal clearance and are potent mediators of cytokine release. Both predictions were confirmed by in vitro experiments on established cell lines as well as primary lung cells. Further, our model affirms the importance of neutrophils in invasive aspergillosis and underlines that the role of macrophages remains elusive. We expect that our model will contribute to improvement of treatment protocols by focusing on the critical components of immune response to fungi but also fungal virulence traits.

Metabolic Functions of Gut Microbes Associate With Efficacy of Tumor Necrosis Factor Antagonists in Patients With Inflammatory Bowel Diseases.

BACKGROUND & AIMS: Altered interactions between the mucosal immune system and intestinal microbiota contribute to pathogenesis of inflammatory bowel diseases (IBD). It is not clear how inhibitors of cytokines, such as antagonists of tumor necrosis factor (anti-TNF), affect the intestinal microbiome. We investigated the effects of anti-TNF agents on gut microbe community structure and function in a longitudinal 2-step study of patients with IBD. We correlated our findings with outcomes of treatment and investigated patterns of metabolites in fecal samples before and after anti-TNF therapy. METHODS: We performed a prospective study of 2 cohorts of patients in Germany; the discovery cohort comprised 12 patients with IBD, 17 patients with rheumatic disease, and 19 healthy individuals (controls); fecal samples were collected at baseline and 2, 6, and 30 weeks after induction of anti-TNF therapy. The validation cohort comprised 23 patients with IBD treated with anti-TNF or vedolizumab (anti-α4ß7 integrin) and 99 healthy controls; fecal samples were collected at baseline and at weeks 2, 6, and 14. Fecal microbiota were analyzed by V3-V4 16S ribosomal RNA gene amplicon sequencing. Clinical response and remission were determined by clinical disease activity scores. Metabolic network reconstruction and associated fecal metabolite level inference was performed in silico using the AGORA (Assembly of Gut Organisms through Reconstruction and Analysis) resource. Metabolomic analyses of fecal samples from a subset of patients were performed to validate metabolites associated with treatment outcomes. RESULTS: Anti-TNF therapy shifted the diversity of fecal microbiota in patients with IBD, but not with rheumatic disease, toward that of controls. Across timepoints, diversity indices did not vary significantly between patients with IBD who did or did not achieve clinical remission after therapy. In contrast, in silico modeling of metabolic interactions between gut microbes found metabolite exchange to be significantly reduced at baseline in fecal samples from patients with IBD and to be associated with later clinical remission. Predicted levels of butyrate and substrates involved in butyrate synthesis (ethanol or acetaldehyde) were significantly associated with clinical remission following anti-TNF therapy, verified by fecal metabolomic analyses. CONCLUSIONS: Metabolic network reconstruction and assessment of metabolic profiles of fecal samples might be used to identify patients with IBD likely to achieve clinical remission following anti-TNF therapy and increase our understanding of the heterogeneity of IBD.

Epithelial RNase H2 Maintains Genome Integrity and Prevents Intestinal Tumorigenesis in Mice.

BACKGROUND & AIMS: RNase H2 is a holoenzyme, composed of 3 subunits (ribonuclease H2 subunits A, B, and C), that cleaves RNA:DNA hybrids and removes mis-incorporated ribonucleotides from genomic DNA through ribonucleotide excision repair. Ribonucleotide incorporation by eukaryotic DNA polymerases occurs during every round of genome duplication and produces the most frequent type of naturally occurring DNA lesion. We investigated whether intestinal epithelial proliferation requires RNase H2 function and whether RNase H2 activity is disrupted during intestinal carcinogenesis. METHODS: We generated mice with epithelial-specific deletion of ribonuclease H2 subunit B (H2bΔIEC) and mice that also had deletion of tumor-suppressor protein p53 (H2b/p53ΔIEC); we compared phenotypes with those of littermate H2bfl/fl or H2b/p53fl/fl (control) mice at young and old ages. Intestinal tissues were collected and analyzed by histology. We isolated epithelial cells, generated intestinal organoids, and performed RNA sequence analyses. Mutation signatures of spontaneous tumors from H2b/p53ΔIEC mice were characterized by exome sequencing. We collected colorectal tumor specimens from 467 patients, measured levels of ribonuclease H2 subunit B, and associated these with patient survival times and transcriptome data. RESULTS: The H2bΔIEC mice had DNA damage to intestinal epithelial cells and proliferative exhaustion of the intestinal stem cell compartment compared with controls and H2b/p53ΔIEC mice. However, H2b/p53ΔIEC mice spontaneously developed small intestine and colon carcinomas. DNA from these tumors contained T>G base substitutions at GTG trinucleotides. Analyses of transcriptomes of human colorectal tumors associated lower levels of RNase H2 with shorter survival times. CONCLUSIONS: In analyses of mice with disruption of the ribonuclease H2 subunit B gene and colorectal tumors from patients, we provide evidence that RNase H2 functions as a colorectal tumor suppressor. H2b/p53ΔIEC mice can be used to study the roles of RNase H2 in tissue-specific carcinogenesis.

Optimality principles reveal a complex interplay of intermediate toxicity and kinetic efficiency in the regulation of prokaryotic metabolism.

A precise and rapid adjustment of fluxes through metabolic pathways is crucial for organisms to prevail in changing environmental conditions. Based on this reasoning, many guiding principles that govern the evolution of metabolic networks and their regulation have been uncovered. To this end, methods from dynamic optimization are ideally suited since they allow to uncover optimality principles behind the regulation of metabolic networks. We used dynamic optimization to investigate the influence of toxic intermediates in connection with the efficiency of enzymes on the regulation of a linear metabolic pathway. Our results predict that transcriptional regulation favors the control of highly efficient enzymes with less toxic upstream intermediates to reduce accumulation of toxic downstream intermediates. We show that the derived optimality principles hold by the analysis of the interplay between intermediate toxicity and pathway regulation in the metabolic pathways of over 5000 sequenced prokaryotes. Moreover, using the lipopolysaccharide biosynthesis in Escherichia coli as an example, we show how knowledge about the relation of regulation, kinetic efficiency and intermediate toxicity can be used to identify drug targets, which control endogenous toxic metabolites and prevent microbial growth. Beyond prokaryotes, we discuss the potential of our findings for the development of antifungal drugs.

BacArena: Individual-based metabolic modeling of heterogeneous microbes in complex communities.

Recent advances focusing on the metabolic interactions within and between cellular populations have emphasized the importance of microbial communities for human health. Constraint-based modeling, with flux balance analysis in particular, has been established as a key approach for studying microbial metabolism, whereas individual-based modeling has been commonly used to study complex dynamics between interacting organisms. In this study, we combine both techniques into the R package BacArena (https://cran.r-project.org/package=BacArena) to generate novel biological insights into Pseudomonas aeruginosa biofilm formation as well as a seven species model community of the human gut. For our P. aeruginosa model, we found that cross-feeding of fermentation products cause a spatial differentiation of emerging metabolic phenotypes in the biofilm over time. In the human gut model community, we found that spatial gradients of mucus glycans are important for niche formations which shape the overall community structure. Additionally, we could provide novel hypothesis concerning the metabolic interactions between the microbes. These results demonstrate the importance of spatial and temporal multi-scale modeling approaches such as BacArena.

Uncoupling of mucosal gene regulation, mRNA splicing and adherent microbiota signatures in inflammatory bowel disease.

OBJECTIVE: An inadequate host response to the intestinal microbiota likely contributes to the manifestation and progression of human inflammatory bowel disease (IBD). However, molecular approaches to unravelling the nature of the defective crosstalk and its consequences for intestinal metabolic and immunological networks are lacking. We assessed the mucosal transcript levels, splicing architecture and mucosa-attached microbial communities of patients with IBD to obtain a comprehensive view of the underlying, hitherto poorly characterised interactions, and how these are altered in IBD. DESIGN: Mucosal biopsies from Crohn's disease and patients with UC, disease controls and healthy individuals (n=63) were subjected to microbiome, transcriptome and splicing analysis, employing next-generation sequencing. The three data levels were integrated by different bioinformatic approaches, including systems biology-inspired network and pathway analysis. RESULTS: Microbiota, host transcript levels and host splicing patterns were influenced most strongly by tissue differences, followed by the effect of inflammation. Both factors point towards a substantial disease-related alteration of metabolic processes. We also observed a strong enrichment of splicing events in inflamed tissues, accompanied by an alteration of the mucosa-attached bacterial taxa. Finally, we noted a striking uncoupling of the three molecular entities when moving from healthy individuals via disease controls to patients with IBD. CONCLUSIONS: Our results provide strong evidence that the interplay between microbiome and host transcriptome, which normally characterises a state of intestinal homeostasis, is drastically perturbed in Crohn's disease and UC. Consequently, integrating multiple OMICs levels appears to be a promising approach to further disentangle the complexity of IBD.

Deciphering the regulation of metabolism with dynamic optimization: an overview of recent advances.

Understanding optimality principles shaping the evolution of regulatory networks controlling metabolism is crucial for deriving a holistic picture of how metabolism is integrated into key cellular processes such as growth, adaptation and pathogenicity. While in the past the focus of research in pathway regulation was mainly based on stationary states, more recently dynamic optimization has proved to be an ideal tool to decipher regulatory strategies for metabolic pathways in response to environmental cues. In this short review, we summarize recent advances in the elucidation of optimal regulatory strategies and identification of optimal control points in metabolic pathways. We discuss biological implications of the discovered optimality principles on genome organization and provide examples how the derived knowledge can be used to identify new treatment strategies against pathogens. Furthermore, we briefly discuss the variety of approaches for solving dynamic optimization problems and emphasize whole-cell resource allocation models as an important emerging area of research that will allow us to study the regulation of metabolism on the whole-cell level.

Pervasive Selection for Cooperative Cross-Feeding in Bacterial Communities.

Bacterial communities are taxonomically highly diverse, yet the mechanisms that maintain this diversity remain poorly understood. We hypothesized that an obligate and mutual exchange of metabolites, as is very common among bacterial cells, could stabilize different genotypes within microbial communities. To test this, we developed a cellular automaton to model interactions among six empirically characterized genotypes that differ in their ability and propensity to produce amino acids. By systematically varying intrinsic (i.e. benefit-to-cost ratio) and extrinsic parameters (i.e. metabolite diffusion level, environmental amino acid availability), we show that obligate cross-feeding of essential metabolites is selected for under a broad range of conditions. In spatially structured environments, positive assortment among cross-feeders resulted in the formation of cooperative clusters, which limited exploitation by non-producing auxotrophs, yet allowed them to persist at the clusters' periphery. Strikingly, cross-feeding helped to maintain genotypic diversity within populations, while amino acid supplementation to the environment decoupled obligate interactions and favored auxotrophic cells that saved amino acid production costs over metabolically autonomous prototrophs. Together, our results suggest that spatially structured environments and limited nutrient availabilities should facilitate the evolution of metabolic interactions, which can help to maintain genotypic diversity within natural microbial populations.

Computing autocatalytic sets to unravel inconsistencies in metabolic network reconstructions.

MOTIVATION: Genome-scale metabolic network reconstructions have been established as a powerful tool for the prediction of cellular phenotypes and metabolic capabilities of organisms. In recent years, the number of network reconstructions has been constantly increasing, mostly because of the availability of novel (semi-)automated procedures, which enabled the reconstruction of metabolic models based on individual genomes and their annotation. The resulting models are widely used in numerous applications. However, the accuracy and predictive power of network reconstructions are commonly limited by inherent inconsistencies and gaps. RESULTS: Here we present a novel method to validate metabolic network reconstructions based on the concept of autocatalytic sets. Autocatalytic sets correspond to collections of metabolites that, besides enzymes and a growth medium, are required to produce all biomass components in a metabolic model. These autocatalytic sets are well-conserved across all domains of life, and their identification in specific genome-scale reconstructions allows us to draw conclusions about potential inconsistencies in these models. The method is capable of detecting inconsistencies, which are neglected by other gap-finding methods. We tested our method on the Model SEED, which is the largest repository for automatically generated genome-scale network reconstructions. In this way, we were able to identify a significant number of missing pathways in several of these reconstructions. Hence, the method we report represents a powerful tool to identify inconsistencies in large-scale metabolic networks. AVAILABILITY AND IMPLEMENTATION: The method is available as source code on http://users.minet.uni-jena.de/∼m3kach/ASBIG/ASBIG.zip. CONTACT: christoph.kaleta@uni-jena.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Optimal programs of pathway control: dissecting the influence of pathway topology and feedback inhibition on pathway regulation.

BACKGROUND: Adjusting the capacity of metabolic pathways in response to rapidly changing environmental conditions is an important component of microbial adaptation strategies to stochastic environments. In this work, we use advanced dynamic optimization techniques combined with theoretical models to study which reactions in pathways are optimally targeted by regulatory interactions in order to minimize the regulatory effort that is required to adjust the flux through a complex metabolic network. Moreover, we analyze how constraints in the speed at which an organism can respond on a proteomic level influences these optimal targets of pathway control. RESULTS: We find that limitations in protein biosynthetic rates have a strong influence. With increasing protein biosynthetic rates the regulatory effort targeting the initial enzyme in a pathway is reduced while the regulatory effort in the terminal enzyme is increased. Studying the impact of allosteric regulation for different pathway topologies, we find that the presence of feedback inhibition by products of metabolic pathways allows organisms to reduce the regulatory effort that is required to control a metabolic pathway in all cases. In a linear pathway this even leads to the case where the sole transcriptional regulatory control of the terminal enzyme is sufficient to control flux through the entire pathway. We confirm the utilization of these pathway regulation strategies through the large-scale analysis of transcriptional regulation in several hundred prokaryotes. CONCLUSIONS: This work expands our knowledge about optimal programs of pathway control. Optimal targets of pathway control strongly depend on the speed at which proteins can be synthesized. Moreover, post-translational regulation such as allosteric regulation allows to strongly reduce the number of transcriptional regulatory interactions required to control a metabolic pathway across different pathway topologies.

Mathematical models for explaining the Warburg effect: a review focussed on ATP and biomass production.

For producing ATP, tumour cells rely on glycolysis leading to lactate to about the same extent as on respiration. Thus, the ATP synthesis flux from glycolysis is considerably higher than in the corresponding healthy cells. This is known as the Warburg effect (named after German biochemist Otto H. Warburg) and also applies to striated muscle cells, activated lymphocytes, microglia, endothelial cells and several other cell types. For similar phenomena in several yeasts and many bacteria, the terms Crabtree effect and overflow metabolism respectively, are used. The Warburg effect is paradoxical at first sight because the molar ATP yield of glycolysis is much lower than that of respiration. Although a straightforward explanation is that glycolysis allows a higher ATP production rate, the question arises why cells do not re-allocate protein to the high-yield pathway of respiration. Mathematical modelling can help explain this phenomenon. Here, we review several models at various scales proposed in the literature for explaining the Warburg effect. These models support the hypothesis that glycolysis allows for a higher proliferation rate due to increased ATP production and precursor supply rates.

Optimality in the zonation of ammonia detoxification in rodent liver.

The rodent liver eliminates toxic ammonia. In mammals, three enzymes (or enzyme systems) are involved in this process: glutaminase, glutamine synthetase and the urea cycle enzymes, represented by carbamoyl phosphate synthetase. The distribution of these enzymes for optimal ammonia detoxification was determined by numerical optimization. This in silico approach predicted that the enzymes have to be zonated in order to achieve maximal removal of toxic ammonia and minimal changes in glutamine concentration. Using 13 compartments, representing hepatocytes, the following predictions were generated: glutamine synthetase is active only within a narrow pericentral zone. Glutaminase and carbamoyl phosphate synthetase are located in the periportal zone in a non-homogeneous distribution. This correlates well with the paradoxical observation that in a first step glutamine-bound ammonia is released (by glutaminase) although one of the functions of the liver is detoxification by ammonia fixation. The in silico approach correctly predicted the in vivo enzyme distributions also for non-physiological conditions (e.g. starvation) and during regeneration after tetrachloromethane (CCl4) intoxication. Metabolite concentrations of glutamine, ammonia and urea in each compartment, representing individual hepatocytes, were predicted. Finally, a sensitivity analysis showed a striking robustness of the results. These bioinformatics predictions were validated experimentally by immunohistochemistry and are supported by the literature. In summary, optimization approaches like the one applied can provide valuable explanations and high-quality predictions for in vivo enzyme and metabolite distributions in tissues and can reveal unknown metabolic functions.

Recurrent Phases of Strict Protein Limitation Inhibit Tumor Growth and Restore Lifespan in A Drosophila Intestinal Cancer Model.

Diets that restrict caloric or protein intake offer a variety of benefits, including decreasing the incidence of cancer. However, whether such diets pose a substantial therapeutic benefit as auxiliary cancer treatments remains unclear. We determined the effects of severe protein depletion on tumorigenesis in a Drosophila melanogaster intestinal tumor model, using a human RAF gain-of-function allele. Severe and continuous protein restriction significantly reduced tumor growth but resulted in premature death. Therefore, we developed a diet in which short periods of severe protein restriction alternated cyclically with periods of complete feeding. This nutritional regime reduced tumor mass, restored gut functionality, and rescued the lifespan of oncogene-expressing flies to the levels observed in healthy flies on a continuous, fully nutritious diet. Furthermore, this diet reduced the chemotherapy-induced stem cell activity associated with tumor recurrence. Transcriptome analysis revealed long-lasting changes in the expression of key genes involved in multiple major developmental signaling pathways. Overall, the data suggest that recurrent severe protein depletion effectively mimics the health benefits of continuous protein restriction, without undesired nutritional shortcomings. This provides seminal insights into the mechanisms of the memory effect required to maintain the positive effects of protein restriction throughout the phases of a full diet. Finally, the repetitive form of strict protein restriction is an ideal strategy for adjuvant cancer therapy that is useful in many tumor contexts.

Gut-associated functions are favored during microbiome assembly across a major part of C. elegans life.

The microbiome expresses a variety of functions that influence host biology. The range of functions depends on the microbiome's composition, which can change during the host's lifetime due to neutral assembly processes, host-mediated selection, and environmental conditions. To date, the exact dynamics of microbiome assembly, the underlying determinants, and the effects on host-associated functions remain poorly understood. Here, we used the nematode Caenorhabditis elegans and a defined community of fully sequenced, naturally associated bacteria to study microbiome dynamics and functions across a major part of the worm's lifetime of hosts under controlled experimental conditions. Bacterial community composition initially shows strongly declining levels of stochasticity, which increases during later time points, suggesting selective effects in younger animals as opposed to more random processes in older animals. The adult microbiome is enriched in genera Ochrobactrum and Enterobacter compared to the direct substrate and a host-free control environment. Using pathway analysis, metabolic, and ecological modeling, we further find that the lifetime assembly dynamics increase competitive strategies and gut-associated functions in the host-associated microbiome, indicating that the colonizing bacteria benefit the worm. Overall, our study introduces a framework for studying microbiome assembly dynamics based on stochastic, ecological, and metabolic models, yielding new insights into the processes that determine host-associated microbiome composition and function. IMPORTANCE: The microbiome plays a crucial role in host biology. Its functions depend on the microbiome composition that can change during a host's lifetime. To date, the dynamics of microbiome assembly and the resulting functions still need to be better understood. This study introduces a new approach to characterize the functional consequences of microbiome assembly by modeling both the relevance of stochastic processes and metabolic characteristics of microbial community changes. The approach was applied to experimental time-series data obtained for the microbiome of the nematode Caenorhabditis elegans across the major part of its lifetime. Stochastic processes played a minor role, whereas beneficial bacteria as well as gut-associated functions enriched in hosts. This indicates that the host might actively shape the composition of its microbiome. Overall, this study provides a framework for studying microbiome assembly dynamics and yields new insights into C. elegans microbiome functions.

Changes in Bacterial Gut Composition in Parkinson's Disease and Their Metabolic Contribution to Disease Development: A Gut Community Reconstruction Approach.

Parkinson's disease (PD) is a chronic and progressive neurodegenerative disease with the major symptoms comprising loss of movement coordination (motor dysfunction) and non-motor dysfunction, including gastrointestinal symptoms. Alterations in the gut microbiota composition have been reported in PD patients vs. controls. However, it is still unclear how these compositional changes contribute to disease etiology and progression. Furthermore, most of the available studies have focused on European, Asian, and North American cohorts, but the microbiomes of PD patients in Latin America have not been characterized. To address this problem, we obtained fecal samples from Colombian participants (n = 25 controls, n = 25 PD idiopathic cases) to characterize the taxonomical community changes during disease via 16S rRNA gene sequencing. An analysis of differential composition, diversity, and personalized computational modeling was carried out, given the fecal bacterial composition and diet of each participant. We found three metabolites that differed in dietary habits between PD patients and controls: carbohydrates, trans fatty acids, and potassium. We identified six genera that changed significantly in their relative abundance between PD patients and controls, belonging to the families Lachnospiraceae, Lactobacillaceae, Verrucomicrobioaceae, Peptostreptococcaceae, and Streptococcaceae. Furthermore, personalized metabolic modeling of the gut microbiome revealed changes in the predicted production of seven metabolites (Indole, tryptophan, fructose, phenylacetic acid, myristic acid, 3-Methyl-2-oxovaleric acid, and N-Acetylneuraminic acid). These metabolites are associated with the metabolism of aromatic amino acids and their consumption in the diet. Therefore, this research suggests that each individual's diet and intestinal composition could affect host metabolism. Furthermore, these findings open the door to the study of microbiome-host interactions and allow us to contribute to personalized medicine.

Elevated fecal calprotectin is associated with gut microbial dysbiosis, altered serum markers and clinical outcomes in older individuals.

Fecal calprotectin is an established marker of gut inflammation in inflammatory bowel disease (IBD). Elevated levels of fecal calprotectin as well as gut microbial dysbiosis have also been observed in other clinical conditions. However, systemic and multi-omics alterations linked to elevated fecal calprotectin in older individuals remain unclear. This study comprehensively investigated the relationship between fecal calprotectin levels, gut microbiome composition, serum inflammation and targeted metabolomics markers, and relevant lifestyle and medical data in a large sample of older individuals (n = 735; mean age ± SD: 68.7 ± 6.3) from the TREND cohort study. Low (0-50 µg/g; n = 602), moderate (> 50-100 µg/g; n = 64) and high (> 100 µg/g; n = 62) fecal calprotectin groups were stratified. Several pro-inflammatory gut microbial genera were significantly increased and short-chain fatty acid producing genera were decreased in high vs. low calprotectin groups. In serum, IL-17C, CCL19 and the toxic metabolite indoxyl sulfate were increased in high vs. low fecal calprotectin groups. These changes were partially mediated by the gut microbiota. Moreover, the high fecal calprotectin group showed increased BMI and a higher disease prevalence of heart attack and obesity. Our findings contribute to the understanding of fecal calprotectin as a marker of gut dysbiosis and its broader systemic and clinical implications in older individuals.

Nonlinear DNA methylation trajectories in aging male mice.

Although DNA methylation data yields highly accurate age predictors, little is known about the dynamics of this quintessential epigenomic biomarker during lifespan. To narrow the gap, we investigate the methylation trajectories of male mouse colon at five different time points of aging. Our study indicates the existence of sudden hypermethylation events at specific stages of life. Precisely, we identify two epigenomic switches during early-to-midlife (3-9 months) and mid-to-late-life (15-24 months) transitions, separating the rodents' life into three stages. These nonlinear methylation dynamics predominantly affect genes associated with the nervous system and enrich in bivalently marked chromatin regions. Based on groups of nonlinearly modified loci, we construct a clock-like classifier STageR (STage of aging estimatoR) that accurately predicts murine epigenetic stage. We demonstrate the universality of our clock in an independent mouse cohort and with publicly available datasets.

Metabolic model predictions enable targeted microbiome manipulation through precision prebiotics.

While numerous health-beneficial interactions between host and microbiota have been identified, there is still a lack of targeted approaches for modulating these interactions. Thus, we here identify precision prebiotics that specifically modulate the abundance of a microbiome member species of interest. In the first step, we show that defining precision prebiotics by compounds that are only taken up by the target species but no other species in a community is usually not possible due to overlapping metabolic niches. Subsequently, we use metabolic modeling to identify precision prebiotics for a two-member Caenorhabditis elegans microbiome community comprising the immune-protective target species Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71. We experimentally confirm four of the predicted precision prebiotics, L-serine, L-threonine, D-mannitol, and γ-aminobutyric acid, to specifically increase the abundance of MYb11. L-serine was further assessed in vivo, leading to an increase in MYb11 abundance also in the worm host. Overall, our findings demonstrate that metabolic modeling is an effective tool for the design of precision prebiotics as an important cornerstone for future microbiome-targeted therapies.IMPORTANCEWhile various mechanisms through which the microbiome influences disease processes in the host have been identified, there are still only few approaches that allow for targeted manipulation of microbiome composition as a first step toward microbiome-based therapies. Here, we propose the concept of precision prebiotics that allow to boost the abundance of already resident health-beneficial microbial species in a microbiome. We present a constraint-based modeling pipeline to predict precision prebiotics for a minimal microbial community in the worm Caenorhabditis elegans comprising the host-beneficial Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71 with the aim to boost the growth of MYb11. Experimentally testing four of the predicted precision prebiotics, we confirm that they are specifically able to increase the abundance of MYb11 in vitro and in vivo. These results demonstrate that constraint-based modeling could be an important tool for the development of targeted microbiome-based therapies against human diseases.