Search

Mapping enhancers and their target genes in specific cell types is crucial for understanding gene regulation and human disease genetics. However, accurately predicting enhancer-gene regulatory interactions from single-cell datasets has been challenging. Here, we introduce a new family of classification models, scE2G, to predict enhancer-gene regulation. These models use features from single-cell ATAC-seq or multiomic RNA and ATAC-seq data and are trained on a CRISPR perturbation dataset including >10,000 evaluated element-gene pairs. We benchmark scE2G models against CRISPR perturbations, fine-mapped eQTLs, and GWAS variant-gene associations and demonstrate state-of-the-art performance at prediction tasks across multiple cell types and categories of perturbations. We apply scE2G to build maps of enhancer-gene regulatory interactions in heterogeneous tissues and interpret noncoding variants associated with complex traits, nominating regulatory interactions linking INPP4B and IL15 to lymphocyte counts. The scE2G models will enable accurate mapping of enhancer-gene regulatory interactions across thousands of diverse human cell types.

Metastasis is the leading cause of cancer-related deaths. It is unclear how intratumor heterogeneity (ITH) contributes to metastasis and how metastatic cells adapt to distant tissue environments. The study of these adaptations is challenged by the limited access to patient material and a lack of experimental models that appropriately recapitulate ITH. To investigate metastatic cell adaptations and the contribution of ITH to metastasis, we analyzed single-cell transcriptomes of matched primary tumors and metastases from patient-derived xenograft models of breast cancer. We found profound transcriptional differences between the primary tumor and metastatic cells. Primary tumors upregulated several metabolic genes, whereas motility pathway genes were upregulated in micrometastases, and stress response signaling was upregulated during progression. Additionally, we identified primary tumor gene signatures that were associated with increased metastatic potential and correlated with patient outcomes. Immune-regulatory control pathways were enriched in poorly metastatic primary tumors, whereas genes involved in epithelial-mesenchymal transition were upregulated in highly metastatic tumors. We found that ITH was dominated by epithelial-mesenchymal plasticity (EMP), which presented as a dynamic continuum with intermediate EMP cell states characterized by specific genes such as CRYAB and S100A2. Elevated expression of an intermediate EMP signature correlated with worse patient outcomes. Our findings identified inhibition of the intermediate EMP cell state as a potential therapeutic target to block metastasis.

Tumor metastasis requires systemic remodeling of distant organ microenvironments that impacts immune cell phenotypes, population structure, and intercellular communication. However, our understanding of immune phenotypic dynamics in the metastatic niche remains incomplete. Here, we longitudinally assayed lung immune transcriptional profiles in the polyomavirus middle T antigen (PyMT) and 4T1 metastatic breast cancer models from primary tumorigenesis, through pre-metastatic niche formation, to the final stages of metastatic outgrowth at single-cell resolution. Computational analyses of these data revealed a TLR-NFκB inflammatory program enacted by both peripherally derived and tissue-resident myeloid cells that correlated with pre-metastatic niche formation and mirrored CD14 + "activated" myeloid cells in the primary tumor. Moreover, we observed that primary tumor and metastatic niche natural killer (NK) cells are differentially regulated in mice and human patient samples, with the metastatic niche featuring elevated cytotoxic NK cell proportions. Finally, we identified cell-type-specific dynamic regulation of IGF1 and CCL6 signaling during metastatic progression that represents anti-metastatic immunotherapy candidate pathways., Competing Interests: Declaration of interests C.S.M. holds patents related to MULTI-seq. A.T.S. is a founder of Immunai, Cartography Biosciences, and Prox Biosciences, an advisor to Zafrens and Wing Venture Capital, and receives research funding from Merck Research Laboratories., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Tissue folds are structural motifs critical to organ function. In the intestine, bending of a flat epithelium into a periodic pattern of folds gives rise to villi, finger-like protrusions that enable nutrient absorption. However, the molecular and mechanical processes driving villus morphogenesis remain unclear. Here, we identify an active mechanical mechanism that simultaneously patterns and folds the intestinal epithelium to initiate villus formation. At the cellular level, we find that PDGFRA+ subepithelial mesenchymal cells generate myosin II-dependent forces sufficient to produce patterned curvature in neighboring tissue interfaces. This symmetry-breaking process requires altered cell and extracellular matrix interactions that are enabled by matrix metalloproteinase-mediated tissue fluidization. Computational models, together with in vitro and in vivo experiments, revealed that these cellular features manifest at the tissue level as differences in interfacial tensions that promote mesenchymal aggregation and interface bending through a process analogous to the active dewetting of a thin liquid film., Competing Interests: Declaration of interests Z.J.G. and C.S.M. hold patents related to the MULTI-seq barcoding method., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Gene regulatory networks within cells modulate the expression of the genome in response to signals and changing environmental conditions. Reconstructions of gene regulatory networks can reveal the information processing and control principles used by cells to maintain homeostasis and execute cell-state transitions. Here, we introduce a computational framework, D-SPIN, that generates quantitative models of gene regulatory networks from single-cell mRNA-seq datasets collected across thousands of distinct perturbation conditions. D-SPIN models the cell as a collection of interacting gene-expression programs, and constructs a probabilistic model to infer regulatory interactions between gene-expression programs and external perturbations. Using large Perturb-seq and drug-response datasets, we demonstrate that D-SPIN models reveal the organization of cellular pathways, sub-functions of macromolecular complexes, and the logic of cellular regulation of transcription, translation, metabolism, and protein degradation in response to gene knockdown perturbations. D-SPIN can also be applied to dissect drug response mechanisms in heterogeneous cell populations, elucidating how combinations of immunomodulatory drugs can induce novel cell states through additive recruitment of gene expression programs. D-SPIN provides a computational framework for constructing interpretable models of gene-regulatory networks to reveal principles of cellular information processing and physiological control., Competing Interests: Declaration of Interests The work was supported in part by funds from a Caltech Amgen Research Collaboration Award and reagent gifts from 10x Genomics. MWT received funding from Adaptive Biotechnologies for unrelated work. MWT is a member of the advisory board of Cell Systems. MWT is a co-founder of CognitiveAI and YurtsAI. ZJG is a co-founder of Scribe Biosciences. The Regents of the University of California with ZJG and CSM as inventors have filed patent applications related to MULTI-seq.

A key aspect of nutrient absorption is the exquisite division of labour across the length of the small intestine, with individual nutrients taken up at different proximal:distal positions. For millennia, the small intestine was thought to comprise three segments with indefinite borders: the duodenum, jejunum and ileum. By examining the fine-scale longitudinal transcriptional patterns that span the mouse and human small intestine, we instead identified five domains of nutrient absorption that mount distinct responses to dietary changes, and three regional stem cell populations. Molecular domain identity can be detected with machine learning, which provides a systematic method to computationally identify intestinal domains in mice. We generated a predictive model of transcriptional control of domain identity and validated the roles of Ppar-δ and Cdx1 in patterning lipid metabolism-associated genes. These findings represent a foundational framework for the zonation of absorption across the mammalian small intestine., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

The rise and fall of estrogen and progesterone across menstrual cycles and during pregnancy regulates breast development and modifies cancer risk. How these hormones impact each cell type in the breast remains poorly understood because they act indirectly through paracrine networks. Using single-cell analysis of premenopausal breast tissue, we reveal a network of coordinated transcriptional programs representing the tissue-level response to changing hormone levels. Our computational approach, DECIPHER-seq, leverages person-to-person variability in breast composition and cell state to uncover programs that co-vary across individuals. We use differences in cell-type proportions to infer a subset of programs that arise from direct cell-cell interactions regulated by hormones. Further, we demonstrate that prior pregnancy and obesity modify hormone responsiveness through distinct mechanisms: obesity reduces the proportion of hormone-responsive cells, whereas pregnancy dampens the direct response of these cells to hormones. Together, these results provide a comprehensive map of the cycling human breast., Competing Interests: Declaration of interests Z.J.G. and C.S.M. hold patents related to the MULTI-seq barcoding method. Z.J.G. is an equity holder in Scribe Biosciences and Provenance bio and a member of the scientific advisory board of Serotiny Bio. C.S.M. is a consultant for ImYoo. Since January 10, 2022, L.M.M. is an employee of Genentech, a member of the Roche group., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Tumor evolution is driven by the progressive acquisition of genetic and epigenetic alterations that enable uncontrolled growth and expansion to neighboring and distal tissues. The study of phylogenetic relationships between cancer cells provides key insights into these processes. Here, we introduced an evolving lineage-tracing system with a single-cell RNA-seq readout into a mouse model of Kras;Trp53(KP)-driven lung adenocarcinoma and tracked tumor evolution from single-transformed cells to metastatic tumors at unprecedented resolution. We found that the loss of the initial, stable alveolar-type2-like state was accompanied by a transient increase in plasticity. This was followed by the adoption of distinct transcriptional programs that enable rapid expansion and, ultimately, clonal sweep of stable subclones capable of metastasizing. Finally, tumors develop through stereotypical evolutionary trajectories, and perturbing additional tumor suppressors accelerates progression by creating novel trajectories. Our study elucidates the hierarchical nature of tumor evolution and, more broadly, enables in-depth studies of tumor progression., Competing Interests: Declaration of interests J.S.W. declares outside interest in 5 AM Venture, Amgen, Chroma Medicine, KSQ Therapeutics, Maze Therapeutics, Tenaya Therapeutics, and Tessera Therapeutics. T.J. is a member of the Board of Directors of Amgen and Thermo Fisher Scientific, is a co-founder of Dragonfly Therapeutics and T2 Biosystems, and is the president of Break Through Cancer. T.J. serves on the Scientific Advisory Board of Dragonfly Therapeutics, SQZ Biotech, and Skyhawk Therapeutics. None of these affiliations represent a conflict of interest with respect to this study. T.G.B. is an advisor to Array BioPharma, Revolution Medicines, Novartis, AstraZeneca, Takeda, Springworks, Jazz Pharmaceuticals, Relay Therapeutics, Rain Therapeutics, and Engine Biosciences and receives research funding from Novartis, Strategia, Kinnate, and Revolution Medicines. J.M.R. consults for Maze Therapeutics and Waypoint Bio. Z.J.G. is an equity holder in Scribe Biosciences and Provenance bio and a member of the SAB of Serotiny Bio., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Enhancers integrate transcription factor signaling pathways that drive cell fate specification in the developing brain. We paired enhancer labeling and single-cell RNA-sequencing (scRNA-seq) to delineate and distinguish specification of neuronal lineages in mouse medial, lateral, and caudal ganglionic eminences (MGE, LGE, and CGE) at embryonic day (E)11.5. We show that scRNA-seq clustering using transcription factors improves resolution of regional and developmental populations, and that enhancer activities identify specific and overlapping GE-derived neuronal populations. First, we mapped the activities of seven evolutionarily conserved brain enhancers at single-cell resolution in vivo, finding that the selected enhancers had diverse activities in specific progenitor and neuronal populations across the GEs. We then applied enhancer-based labeling, scRNA-seq, and analysis of in situ hybridization data to distinguish transcriptionally distinct and spatially defined subtypes of MGE-derived GABAergic and cholinergic projection neurons and interneurons. Our results map developmental origins and specification paths underlying neurogenesis in the embryonic basal ganglia and showcase the power of scRNA-seq combined with enhancer-based labeling to resolve the complex paths of neuronal specification underlying mouse brain development.

Microglia are resident macrophages in the brain that emerge in early development and respond to the local environment by altering their molecular and phenotypic states. Fundamental questions about microglia diversity and function during development remain unanswered because we lack experimental strategies to interrogate their interactions with other cell types and responses to perturbations ex vivo. We compared human microglia states across culture models, including cultured primary and pluripotent stem cell-derived microglia. We developed a "report card" of gene expression signatures across these distinct models to facilitate characterization of their responses across experimental models, perturbations, and disease conditions. Xenotransplantation of human microglia into cerebral organoids allowed us to characterize key transcriptional programs of developing microglia in vitro and reveal that microglia induce transcriptional changes in neural stem cells and decrease interferon signaling response genes. Microglia additionally accelerate the emergence of synchronized oscillatory network activity in brain organoids by modulating synaptic density., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)