Cancers infiltrated with T-cells are associated with a higher likelihood of response to PD-1/PD-L1 blockade. Counterintuitively, a correlation between epithelial–mesenchymal transition (EMT)-related gene expression and T-cell infiltration has been observed across tumor types. Here we demonstrate, using The Cancer Genome Atlas (TCGA) urothelial cancer dataset, that although a gene expression-based measure of infiltrating T-cell abundance and EMT-related gene expression are positively correlated, these signatures convey disparate prognostic information. We further demonstrate that non-hematopoietic stromal cells are a major source of EMT-related gene expression in bulk urothelial cancer transcriptomes. Finally, using a cohort of patients with metastatic urothelial cancer treated with a PD-1 inhibitor, nivolumab, we demonstrate that in patients with T-cell infiltrated tumors, higher EMT/stroma-related gene expression is associated with lower response rates and shorter progression-free and overall survival. Together, our findings suggest a stroma-mediated source of immune resistance in urothelial cancer and provide rationale for co-targeting PD-1 and stromal elements.

Using expression profiles from postmortem prefrontal cortex samples of 624 dementia patients and non‐demented controls, we investigated global disruptions in the co‐regulation of genes in two neurodegenerative diseases, late‐onset Alzheimer's disease (AD) and Huntington's disease (HD). We identified networks of differentially co‐expressed (DC) gene pairs that either gained or lost correlation in disease cases relative to the control group, with the former dominant for both AD and HD and both patterns replicating in independent human cohorts of AD and aging. When aligning networks of DC patterns and physical interactions, we identified a 242‐gene subnetwork enriched for independent AD/HD signatures. This subnetwork revealed a surprising dichotomy of gained/lost correlations among two inter‐connected processes, chromatin organization and neural differentiation, and included DNA methyltransferases, DNMT1 and DNMT3A , of which we predicted the former but not latter as a key regulator. To validate the inter‐connection of these two processes and our key regulator prediction, we generated two brain‐specific knockout (KO) mice and show that Dnmt1 KO signature significantly overlaps with the subnetwork ( P = 3.1 × 10 ⁻¹² ), while Dnmt3a KO signature does not ( P = 0.017).

Cells employ multiple levels of regulation, including transcriptional and translational regulation, that drive core biological processes and enable cells to respond to genetic and environmental changes. Small-molecule metabolites are one category of critical cellular intermediates that can influence as well as be a target of cellular regulations. Because metabolites represent the direct output of protein-mediated cellular processes, endogenous metabolite concentrations can closely reflect cellular physiological states, especially when integrated with other molecular-profiling data. Here we develop and apply a network reconstruction approach that simultaneously integrates six different types of data: endogenous metabolite concentration, RNA expression, DNA variation, DNA–protein binding, protein–metabolite interaction, and protein–protein interaction data, to construct probabilistic causal networks that elucidate the complexity of cell regulation in a segregating yeast population. Because many of the metabolites are found to be under strong genetic control, we were able to employ a causal regulator detection algorithm to identify causal regulators of the resulting network that elucidated the mechanisms by which variations in their sequence affect gene expression and metabolite concentrations. We examined all four expression quantitative trait loci (eQTL) hot spots with colocalized metabolite QTLs, two of which recapitulated known biological processes, while the other two elucidated novel putative biological mechanisms for the eQTL hot spots. It is now possible to score variations in DNA across whole genomes, RNA levels and alternative isoforms, metabolite levels, protein levels and protein state information, protein–protein interactions, and protein–DNA interactions, in a comprehensive fashion in populations of individuals. Interactions among these molecular entities define the complex web of biological processes that give rise to all higher order phenotypes, including disease. The development of analytical approaches that simultaneously integrate different dimensions of data is essential if we are to extract the meaning from large-scale data to elucidate the complexity of living systems. Here, we use a novel Bayesian network reconstruction algorithm that simultaneously integrates DNA variation, RNA levels, metabolite levels, protein–protein interaction data, protein–DNA binding data, and protein–small-molecule interaction data to construct molecular networks in yeast. We demonstrate that these networks can be used to infer causal relationships among genes, enabling the identification of novel genes that modulate cellular regulation. We show that our network predictions either recapitulate known biology or can be prospectively validated, demonstrating a high degree of accuracy in the predicted network.

Though discovered over 100 years ago, the molecular foundation of sporadic Alzheimer’s disease (AD) remains elusive. To better characterize the complex nature of AD, we constructed multiscale causal networks on a large human AD multi-omics dataset, integrating clinical features of AD, DNA variation, and gene- and protein-expression. These probabilistic causal models enabled detection, prioritization and replication of high-confidence master regulators of AD-associated networks, including the top predicted regulator, VGF. Overexpression of neuropeptide precursor VGF in 5xFAD mice partially rescued beta-amyloid-mediated memory impairment and neuropathology. Molecular validation of network predictions downstream of VGF was also achieved in this AD model, with significant enrichment for homologous genes identified as differentially expressed in 5xFAD brains overexpressing VGF. Our findings support a causal role for VGF in protecting against AD pathogenesis and progression.

To identify factors that regulate gut microbiota density and the impact of varied microbiota density on health, we assayed this fundamental ecosystem property in fecal samples across mammals, human disease, and therapeutic interventions. Physiologic features of the host (carrying capacity) and the fitness of the gut microbiota shape microbiota density. Therapeutic manipulation of microbiota density in mice altered host metabolic and immune homeostasis. In humans, gut microbiota density was reduced in Crohn’s disease, ulcerative colitis, and ileal pouch-anal anastomosis. The gut microbiota in recurrent Clostridium difficile infection had lower density and reduced fitness that were restored by fecal microbiota transplantation. Understanding the interplay between microbiota and disease in terms of microbiota density, host carrying capacity, and microbiota fitness provide new insights into microbiome structure and microbiome targeted therapeutics.