Histone Deacetylase Inhibition by Gut Microbe-Generated Short-Chain Fatty Acids Entrains Intestinal Epithelial Circadian Rhythms.

BACKGROUND & AIMS: The circadian clock orchestrates ∼24-hour oscillations of gastrointestinal epithelial structure and function that drive diurnal rhythms in gut microbiota. Here, we use experimental and computational approaches in intestinal organoids to reveal reciprocal effects of gut microbial metabolites on epithelial timekeeping by an epigenetic mechanism. METHODS: We cultured enteroids in media supplemented with sterile supernatants from the altered Schaedler Flora (ASF), a defined murine microbiota. Circadian oscillations of bioluminescent PER2 and Bmal1 were measured in the presence or absence of individual ASF supernatants. Separately, we applied machine learning to ASF metabolomics to identify phase-shifting metabolites. RESULTS: Sterile filtrates from 3 of 7 ASF species (ASF360 Lactobacillus intestinalis, ASF361 Ligilactobacillus murinus, and ASF502 Clostridium species) induced minimal alterations in circadian rhythms, whereas filtrates from 4 ASF species (ASF356 Clostridium species, ASF492 Eubacterium plexicaudatum, ASF500 Pseudoflavonifactor species, and ASF519 Parabacteroides goldsteinii) induced profound, concentration-dependent phase shifts. Random forest classification identified short-chain fatty acid (SCFA) (butyrate, propionate, acetate, and isovalerate) production as a discriminating feature of ASF "shifters." Experiments with SCFAs confirmed machine learning predictions, with a median phase shift of 6.2 hours in murine enteroids. Pharmacologic or botanical histone deacetylase (HDAC) inhibitors yielded similar findings. Further, mithramycin A, an inhibitor of HDAC inhibition, reduced SCFA-induced phase shifts by 20% (P < .05) and conditional knockout of HDAC3 in enteroids abrogated butyrate effects on Per2 expression. Key findings were reproducible in human Bmal1-luciferase enteroids, colonoids, and Per2-luciferase Caco-2 cells. CONCLUSIONS: Gut microbe-generated SCFAs entrain intestinal epithelial circadian rhythms by an HDACi-dependent mechanism, with critical implications for understanding microbial and circadian network regulation of intestinal epithelial homeostasis.

Quantifying cumulative phenotypic and genomic evidence for procedural generation of metabolic network reconstructions.

Genome-scale metabolic network reconstructions (GENREs) are valuable tools for understanding microbial metabolism. The process of automatically generating GENREs includes identifying metabolic reactions supported by sufficient genomic evidence to generate a draft metabolic network. The draft GENRE is then gapfilled with additional reactions in order to recapitulate specific growth phenotypes as indicated with associated experimental data. Previous methods have implemented absolute mapping thresholds for the reactions automatically included in draft GENREs; however, there is growing evidence that integrating annotation evidence in a continuous form can improve model accuracy. There is a need for flexibility in the structure of GENREs to better account for uncertainty in biological data, unknown regulatory mechanisms, and context-specificity associated with data inputs. To address this issue, we present a novel method that provides a framework for quantifying combined genomic, biochemical, and phenotypic evidence for each biochemical reaction during automated GENRE construction. Our method, Constraint-based Analysis Yielding reaction Usage across metabolic Networks (CANYUNs), generates accurate GENREs with a quantitative metric for the cumulative evidence for each reaction included in the network. The structuring of CANYUNs allows for the simultaneous integration of three data inputs while maintaining all supporting evidence for biochemical reactions that may be active in an organism. CANYUNs is designed to maximize the utility of experimental and annotation datasets and to ultimately assist in the curation of the reference datasets used for the automatic construction of metabolic networks. We validated CANYUNs by generating an E. coli K-12 model and compared it to the manually curated reconstruction iML1515. Finally, we demonstrated the use of CANYUNs to build a model by generating an E. coli Nissle CANYUNs model using novel phenotypic data that we collected. This method may address key challenges for the procedural construction of metabolic networks by leveraging uncertainty and redundancy in biological data.

Medusa: Software to build and analyze ensembles of genome-scale metabolic network reconstructions.

Uncertainty in the structure and parameters of networks is ubiquitous across computational biology. In constraint-based reconstruction and analysis of metabolic networks, this uncertainty is present both during the reconstruction of networks and in simulations performed with them. Here, we present Medusa, a Python package for the generation and analysis of ensembles of genome-scale metabolic network reconstructions. Medusa builds on the COBRApy package for constraint-based reconstruction and analysis by compressing a set of models into a compact ensemble object, providing functions for the generation of ensembles using experimental data, and extending constraint-based analyses to ensemble scale. We demonstrate how Medusa can be used to generate ensembles and perform ensemble simulations, and how machine learning can be used in conjunction with Medusa to guide the curation of genome-scale metabolic network reconstructions. Medusa is available under the permissive MIT license from the Python Packaging Index (https://pypi.org) and from github (https://github.com/opencobra/Medusa), and comprehensive documentation is available at https://medusa.readthedocs.io/en/latest.

Transcriptome-guided parsimonious flux analysis improves predictions with metabolic networks in complex environments.

The metabolic responses of bacteria to dynamic extracellular conditions drives not only the behavior of single species, but also entire communities of microbes. Over the last decade, genome-scale metabolic network reconstructions have assisted in our appreciation of important metabolic determinants of bacterial physiology. These network models have been a powerful force in understanding the metabolic capacity that species may utilize in order to succeed in an environment. Increasingly, an understanding of context-specific metabolism is critical for elucidating metabolic drivers of larger phenotypes and disease. However, previous approaches to use network models in concert with omics data to better characterize experimental systems have met challenges due to assumptions necessary by the various integration platforms or due to large input data requirements. With these challenges in mind, we developed RIPTiDe (Reaction Inclusion by Parsimony and Transcript Distribution) which uses both transcriptomic abundances and parsimony of overall flux to identify the most cost-effective usage of metabolism that also best reflects the cell's investments into transcription. Additionally, in biological samples where it is difficult to quantify specific growth conditions, it becomes critical to develop methods that require lower amounts of user intervention in order to generate accurate metabolic predictions. Utilizing a metabolic network reconstruction for the model organism Escherichia coli str. K-12 substr. MG1655 (iJO1366), we found that RIPTiDe correctly identifies context-specific metabolic pathway activity without supervision or knowledge of specific media conditions. We also assessed the application of RIPTiDe to in vivo metatranscriptomic data where E. coli was present at high abundances, and found that our approach also effectively predicts metabolic behaviors of host-associated bacteria. In the setting of human health, understanding metabolic changes within bacteria in environments where growth substrate availability is difficult to quantify can have large downstream impacts on our ability to elucidate molecular drivers of disease-associated dysbiosis across the microbiota. Our results indicate that RIPTiDe may have potential to provide understanding of context-specific metabolism of bacteria within complex communities.

Gram-negative Microbiota Blooms in Premature Twins Discordant for Parenteral Nutrition-associated Cholestasis.

Parenteral nutrition-associated cholestasis (PNAC) causes serious morbidity in the neonatal intensive care unit. Infection with gut-associated bacteria is associated with cholestasis, but the role of intestinal microbiota in PNAC is poorly understood. We examined the composition of stool microbiota from premature twins discordant for PNAC as a strategy to reduce confounding from variables associated with both microbiota and cholestasis. Eighty-four serial stool samples were included from 4 twin sets discordant for PNAC. Random Forests was utilized to determine genera most discriminatory in classifying samples from infants with and without PNAC. In infants with PNAC, we detected a significant increase in the relative abundance of Klebsiella, Veillonella, Enterobacter, and Enterococcus (P < 0.05). Bray-Curtis dissimilarities in infants with PNAC were significantly different (P < 0.05) from infants without PNAC. Our findings warrant further exploration in larger cohorts and experimental models of PNAC to determine if a microbiota signature predicts PNAC, as a basis for future interventions to mitigate liver injury.

Fecal sphingolipids predict parenteral nutrition-associated cholestasis in the neonatal intensive care unit.

BACKGROUND: Parenteral nutrition-associated cholestasis (PNAC) in the neonatal intensive care unit (NICU) causes significant morbidity and associated healthcare costs. Laboratory detection of PNAC currently relies on elevated serum conjugated bilirubin levels in the aftermath of impaired bile flow. Here, we sought to identify fecal biomarkers, which when integrated with clinical data, would better predict risk for developing PNAC. METHODS: Using untargeted metabolomics in 200 serial stool samples from 60 infants, we applied statistical and machine learning approaches to identify clinical features and metabolic biomarkers with the greatest associative potential for risk of developing PNAC. Stools were collected prospectively from infants receiving PN with soybean oil-based lipid emulsion at a level IV NICU. RESULTS: Low birth weight, extreme prematurity, longer duration of PN, and greater number of antibiotic courses were all risk factors for PNAC (P < 0.05). We identified 78 stool biomarkers with early predictive potential (P < 0.05). From these 78 biomarkers, we further identified 12 sphingomyelin lipids with high association for the development of PNAC in precholestasis stool samples when combined with birth anthropometry. CONCLUSION: We demonstrate the potential for stool metabolomics to enhance early identification of PNAC risk. Earlier detection of high-risk infants would empower proactive mitigation with alterations to PN for at-risk infants and optimization of energy nutrition with PN for infants at lower risk.

Designing Consumer Health Information Technology to Support Biform and Articulation Work: A Qualitative Study of Diet and Nutrition Management as Patient Work.

BACKGROUND: Diet and nutrition management is an integral component of Crohn disease (CD) management. This type of management is highly variable and individualized and, thus, requires personalized approaches. Consumer health information technology (CHIT) designed to support CD management has typically supported this task as everyday life work and, not necessarily, as illness work. Moreover, CHIT has rarely supported the ways in which diet and nutrition management requires coordination between multiple forms of patient work. OBJECTIVE: The purpose of this study was to investigate diet and nutrition management as biform work, identify components of articulation work, and provide guidance on how to design CHIT to support this work. METHODS: We performed a qualitative study in which we recruited participants from CD-related Facebook pages and groups. RESULTS: Semistructured interviews with 21 individuals showed that diet and nutrition management strategies were highly individualized and variable. Four themes emerged from the data, emphasizing the interactions of diet and nutrition with physical, emotional, information, and technology-enabled management. CONCLUSIONS: This study shows that the extent to which diet and nutrition management is biform work fluctuates over time and that articulation work can be continuous and unplanned. The design guidance specifies the need for patient-facing technologies to support interactions among diet and nutrition and other management activities such as medication intake, stress reduction, and information seeking, as well as to respond to the ways in which diet and nutrition management needs change over time.

Novel co-culture plate enables growth dynamic-based assessment of contact-independent microbial interactions.

Interactions between microbes are central to the dynamics of microbial communities. Understanding these interactions is essential for the characterization of communities, yet challenging to accomplish in practice. There are limited available tools for characterizing diffusion-mediated, contact-independent microbial interactions. A practical and widely implemented technique in such characterization involves the simultaneous co-culture of distinct bacterial species and subsequent analysis of relative abundance in the total population. However, distinguishing between species can be logistically challenging. In this paper, we present a low-cost, vertical membrane, co-culture plate to quantify contact-independent interactions between distinct bacterial populations in co-culture via real-time optical density measurements. These measurements can be used to facilitate the analysis of the interaction between microbes that are physically separated by a semipermeable membrane yet able to exchange diffusible molecules. We show that diffusion across the membrane occurs at a sufficient rate to enable effective interaction between physically separate cultures. Two bacterial species commonly found in the cystic fibrotic lung, Pseudomonas aeruginosa and Burkholderia cenocepacia, were co-cultured to demonstrate how this plate may be implemented to study microbial interactions. We have demonstrated that this novel co-culture device is able to reliably generate real-time measurements of optical density data that can be used to characterize interactions between microbial species.

Systems-level metabolism of the altered Schaedler flora, a complete gut microbiota.

The altered Schaedler flora (ASF) is a model microbial community with both in vivo and in vitro relevance. Here we provide the first characterization of the ASF community in vitro, independent of a murine host. We compared the functional genetic content of the ASF to wild murine metagenomes and found that the ASF functionally represents wild microbiomes better than random consortia of similar taxonomic composition. We developed a chemically defined medium that supported growth of seven of the eight ASF members. To elucidate the metabolic capabilities of these ASF species-including potential for interactions such as cross-feeding-we performed a spent media screen and analyzed the results through dynamic growth measurements and non-targeted metabolic profiling. We found that cross-feeding is relatively rare (32 of 3570 possible cases), but is enriched between Clostridium ASF356 and Parabacteroides ASF519. We identified many cases of emergent metabolism (856 of 3570 possible cases). These data will inform efforts to understand ASF dynamics and spatial distribution in vivo, to design pre- and probiotics that modulate relative abundances of ASF members, and will be essential for validating computational models of ASF metabolism. Well-characterized, experimentally tractable microbial communities enable research that can translate into more effective microbiome-targeted therapies to improve human health.

Miniaturized plate readers for low-cost, high-throughput phenotypic screening.

We present a miniaturized plate reader for measuring optical density in 96-well plates. Our standalone reader fits in most incubators, environmental chambers, or biological containment suites, allowing users to leverage their existing laboratory infrastructure. The device contains no moving parts, allowing an entire 96-well plate to be read several times per second. We demonstrate how the fast sampling rate allows our reader to detect small changes in optical density, even when the device is placed in a shaking incubator. A wireless communication module allows remote monitoring of multiple devices in real time. These features allow easy assembly of multiple readers to create a scalable, accurate solution for high-throughput phenotypic screening.

Micropatterned dermal-epidermal regeneration matrices create functional niches that enhance epidermal morphogenesis.

Although tissue engineered skin substitutes have demonstrated some clinical success for the treatment of chronic wounds such as diabetic and venous ulcers, persistent graft take and stability remain concerns. Current bilayered skin substitutes lack the characteristic microtopography of the dermal-epidermal junction that gives skin enhanced mechanical stability and creates cellular microniches that differentially promote keratinocyte function to form skin appendages and enhance wound healing. We developed a novel micropatterned dermal-epidermal regeneration matrix (µDERM) which incorporates this complex topography and substantially enhances epidermal morphology. Here, we describe the use of this three-dimensional (3-D) in vitro culture model to systematically evaluate different topographical geometries and to determine their relationship to keratinocyte function. We identified three distinct keratinocyte functional niches: the proliferative niche (narrow geometries), the basement membrane protein synthesis niche (wide geometries) and the putative keratinocyte stem cell niche (narrow geometries and corners). Specifically, epidermal thickness and keratinocyte proliferation is significantly (p<0.05) increased in 50 and 100 µm channels while laminin-332 deposition is significantly (p<0.05) increased in 400 µm channels compared to flat controls. Additionally, ß1(bri)p63(+) keratinocytes, putative keratinocyte stem cells, preferentially cluster in channel geometries (similar to clustering observed in native skin) compared to a random distribution on flats. This study identifies specific target geometries to enhance skin regeneration and graft performance. Furthermore, these results suggest the importance of µDERM microtopography in designing the next generation of skin substitutes. Finally, we anticipate that 3-D organotypic cultures on µDERMS will provide a novel tissue engineered skin substitute for in vitro investigations of skin morphogenesis, wound healing and pathology.