Background: Humans live in constant and vital symbiosis with a closely linked bacterial ecosystem called the microbiome, which influences many aspects of human health. When this microbial ecosystem becomes disrupted, the health of the human host can suffer; a condition called dysbiosis. However, the community compositions of human microbiomes also vary dramatically from individual to individual, and over time, making it difficult to uncover the underlying mechanisms linking the microbiome to human health. We propose that a microbiome's interaction with its human host is not necessarily dependent upon the presence or absence of particular bacterial species, but instead is dependent on its community metabolome; an emergent property of the microbiome. Results: Using data from a previously published, longitudinal study of microbiome populations of the human gut, we extrapolated information about microbiome community enzyme profiles and metabolome models. Using machine learning techniques, we demonstrated that the aggregate predicted community enzyme function profiles and modeled metabolomes of a microbiome are more predictive of dysbiosis than either observed microbiome community composition or predicted enzyme function profiles. Conclusions: Specific enzyme functions and metabolites predictive of dysbiosis provide insights into the molecular mechanisms of microbiome–host interactions. The ability to use machine learning to predict dysbiosis from microbiome community interaction data provides a potentially powerful tool for understanding the links between the human microbiome and human health, pointing to potential microbiome-based diagnostics and therapeutic interventions.

Accurate prediction of the host phenotype from a metagenomic sample and identification of the associated microbial markers are important in understanding potential host-microbiome interactions related to disease initiation and progression. We introduce PopPhy-CNN, a novel convolutional neural network (CNN) learning framework that effectively exploits phylogenetic structure in microbial taxa for host phenotype prediction. Our approach takes an input format of a 2D matrix representing the phylogenetic tree populated with the relative abundance of microbial taxa in a metagenomic sample. This conversion empowers CNNs to explore the spatial relationship of the taxonomic annotations on the tree and their quantitative characteristics in metagenomic data. We show the competitiveness of our model compared to other available methods using nine metagenomic datasets of moderate size for binary classification. With synthetic and biological datasets, we show the superior and robust performance of our model for multi- class classification. Furthermore, we design a novel scheme for feature extraction from the learned CNN models and demonstrate improved performance when the extracted features. PopPhy-CNN is a practical deep learning framework for the prediction of host phenotype with the ability of facilitating the retrieval of predictive microbial taxa.

Senescent cells accumulate in most tissues with organismal aging, exposure to stressors, or disease progression. It is challenging to identify senescent cells because cellular senescence signatures and phenotypes vary widely across distinct cell types and tissues. Here we developed an analytical algorithm that defines cell-type-specific and universal signatures of cellular senescence across a wide range of cell types and tissues. We utilize 72 mouse and 64 human weighted single-cell transcriptomic signatures of cellular senescence to create the SenePy scoring platform. SenePy signatures better recapitulate in vivo cellular senescence than signatures derived from in vitro senescence studies. We use SenePy to map the kinetics of senescent cell accumulation in healthy aging as well as multiple disease contexts, including tumorigenesis, inflammation, and myocardial infarction. SenePy characterizes cell-type-specific in vivo cellular senescence and could lead to the identification of genes that serve as mediators of cellular senescence and disease progression.

The advance in microbiome and metabolome studies has generated rich omics data revealing the involvement of the microbial community in host disease pathogenesis through interactions with their host at a metabolic level. However, the computational tools to uncover these relationships are just emerging. Here, we present MiMeNet, a neural network framework for modeling microbe-metabolite relationships. Using ten iterations of 10-fold cross-validation on three paired microbiome-metabolome datasets, we show that MiMeNet more accurately predicts metabolite abundances (mean Spearman correlation coefficients increase from 0.108 to 0.309, 0.276 to 0.457, and -0.272 to 0.264) and identifies more well-predicted metabolites (increase in the number of well-predicted metabolites from 198 to 366, 104 to 143, and 4 to 29) compared to state-of-art linear models for individual metabolite predictions. Additionally, we demonstrate that MiMeNet can group microbes and metabolites with similar interaction patterns and functions to illuminate the underlying structure of the microbe-metabolite interaction network, which could potentially shed light on uncharacterized metabolites through “Guilt by Association”. Our results demonstrated that MiMeNet is a powerful tool to provide insights into the causes of metabolic dysregulation in disease, facilitating future hypothesis generation at the interface of the microbiome and metabolomics.

Progesterone, via its nuclear receptor (PR), exerts an overall tumorigenic effect on both uterine fibroid (leiomyoma) and breast cancer tissues, whereas the antiprogestin RU486 inhibits growth of these tissues through an unknown mechanism. Here, we determined the interaction between common or cell-specific genome-wide binding sites of PR and mRNA expression in RU486-treated uterine leiomyoma and breast cancer cells. ChIP-sequencing revealed 31,457 and 7,034 PR-binding sites in breast cancer and uterine leiomyoma cells, respectively; 1,035 sites overlapped in both cell types. Based on the chromatin-PR interaction in both cell types, we statistically refined the consensus progesterone response element to G•ACA• • •TGT•C. We identified two striking differences between uterine leiomyoma and breast cancer cells. First, the cis-regulatory elements for HSF, TEF-1, and C/EBPα and β were statistically enriched at genomic RU486/PR-targets in uterine leiomyoma, whereas E2F, FOXO1, FOXA1, and FOXF sites were preferentially enriched in breast cancer cells. Second, 51.5% of RU486-regulated genes in breast cancer cells but only 6.6% of RU486-regulated genes in uterine leiomyoma cells contained a PR-binding site within 5 kb from their transcription start sites (TSSs), whereas 75.4% of RU486-regulated genes contained a PR-binding site farther than 50 kb from their TSSs in uterine leiomyoma cells. RU486 regulated only seven mRNAs in both cell types. Among these, adipophilin (PLIN2), a pro-differentiation gene, was induced via RU486 and PR via the same regulatory region in both cell types. Our studies have identified molecular components in a RU486/PR-controlled gene network involved in the regulation of cell growth, cell migration, and extracellular matrix function. Tissue-specific and common patterns of genome-wide PR binding and gene regulation may determine the therapeutic effects of antiprogestins in uterine fibroids and breast cancer.

Studies have demonstrated the phenotypic heterogeneity of vascular endothelial cells (ECs) within a vascular bed; however, little is known about how distinct endothelial subpopulations in a particular organ respond to an inflammatory stimulus. We performed single-cell RNA-Seq of 35,973 lung ECs obtained during baseline as well as postinjury time points after inflammatory lung injury induced by LPS. Seurat clustering and gene expression pathway analysis identified 2 major subpopulations in the lung microvascular endothelium, a subpopulation enriched for expression of immune response genes such as MHC genes (immuneEC) and another defined by increased expression of vascular development genes such as Sox17 (devEC). The presence of immuneEC and devEC subpopulations was also observed in nonhuman primate lungs infected with SARS-CoV-2 and murine lungs infected with H1N1 influenza virus. After the peak of inflammatory injury, we observed the emergence of a proliferative lung EC subpopulation. Overexpression of Sox17 prevented inflammatory activation in ECs. Thus, there appeared to be a “division of labor” within the lung microvascular endothelium in which some ECs showed propensity for inflammatory signaling and others for endothelial regeneration. These results provide underpinnings for the development of targeted therapies to limit inflammatory lung injury and promote regeneration.