Toward a foundation model of causal cell and tissue biology with a Perturbation Cell and Tissue Atlas.

Comprehensively charting the biologically causal circuits that govern the phenotypic space of human cells has often been viewed as an insurmountable challenge. However, in the last decade, a suite of interleaved experimental and computational technologies has arisen that is making this fundamental goal increasingly tractable. Pooled CRISPR-based perturbation screens with high-content molecular and/or image-based readouts are now enabling researchers to probe, map, and decipher genetically causal circuits at increasing scale. This scale is now eminently suitable for the deployment of artificial intelligence and machine learning (AI/ML) to both direct further experiments and to predict or generate information that was not-and sometimes cannot-be gathered experimentally. By combining and iterating those through experiments that are designed for inference, we now envision a Perturbation Cell Atlas as a generative causal foundation model to unify human cell biology.

A transcription factor atlas of directed differentiation.

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

Proper acquisition of cell class identity in organoids allows definition of fate specification programs of the human cerebral cortex.

Realizing the full utility of brain organoids to study human development requires understanding whether organoids precisely replicate endogenous cellular and molecular events, particularly since acquisition of cell identity in organoids can be impaired by abnormal metabolic states. We present a comprehensive single-cell transcriptomic, epigenetic, and spatial atlas of human cortical organoid development, comprising over 610,000 cells, from generation of neural progenitors through production of differentiated neuronal and glial subtypes. We show that processes of cellular diversification correlate closely to endogenous ones, irrespective of metabolic state, empowering the use of this atlas to study human fate specification. We define longitudinal molecular trajectories of cortical cell types during organoid development, identify genes with predicted human-specific roles in lineage establishment, and uncover early transcriptional diversity of human callosal neurons. The findings validate this comprehensive atlas of human corticogenesis in vitro as a resource to prime investigation into the mechanisms of human cortical development.

Tissue-resident memory and circulating T cells are early responders to pre-surgical cancer immunotherapy.

Neoadjuvant immune checkpoint blockade has shown promising clinical activity. Here, we characterized early kinetics in tumor-infiltrating and circulating immune cells in oral cancer patients treated with neoadjuvant anti-PD-1 or anti-PD-1/CTLA-4 in a clinical trial (NCT02919683). Tumor-infiltrating CD8 T cells that clonally expanded during immunotherapy expressed elevated tissue-resident memory and cytotoxicity programs, which were already active prior to therapy, supporting the capacity for rapid response. Systematic target discovery revealed that treatment-expanded tumor T cell clones in responding patients recognized several self-antigens, including the cancer-specific antigen MAGEA1. Treatment also induced a systemic immune response characterized by expansion of activated T cells enriched for tumor-infiltrating T cell clonotypes, including both pre-existing and emergent clonotypes undetectable prior to therapy. The frequency of activated blood CD8 T cells, notably pre-treatment PD-1-positive KLRG1-negative T cells, was strongly associated with intra-tumoral pathological response. These results demonstrate how neoadjuvant checkpoint blockade induces local and systemic tumor immunity.

BACH2 regulates diversification of regulatory and proinflammatory chromatin states in TH17 cells.

Interleukin-17 (IL-17)-producing helper T (TH17) cells are heterogenous and consist of nonpathogenic TH17 (npTH17) cells that contribute to tissue homeostasis and pathogenic TH17 (pTH17) cells that mediate tissue inflammation. Here, we characterize regulatory pathways underlying TH17 heterogeneity and discover substantial differences in the chromatin landscape of npTH17 and pTH17 cells both in vitro and in vivo. Compared to other CD4+ T cell subsets, npTH17 cells share accessible chromatin configurations with regulatory T cells, whereas pTH17 cells exhibit features of both npTH17 cells and type 1 helper T (TH1) cells. Integrating single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) and single-cell RNA sequencing (scRNA-seq), we infer self-reinforcing and mutually exclusive regulatory networks controlling different cell states and predicted transcription factors regulating TH17 cell pathogenicity. We validate that BACH2 promotes immunomodulatory npTH17 programs and restrains proinflammatory TH1-like programs in TH17 cells in vitro and in vivo. Furthermore, human genetics implicate BACH2 in multiple sclerosis. Overall, our work identifies regulators of TH17 heterogeneity as potential targets to mitigate autoimmunity.

A cellular and spatial map of the choroid plexus across brain ventricles and ages.

The choroid plexus (ChP) in each brain ventricle produces cerebrospinal fluid (CSF) and forms the blood-CSF barrier. Here, we construct a single-cell and spatial atlas of each ChP in the developing, adult, and aged mouse brain. We delineate diverse cell types, subtypes, cell states, and expression programs in epithelial and mesenchymal cells across ages and ventricles. In the developing ChP, we predict a common progenitor pool for epithelial and neuronal cells, validated by lineage tracing. Epithelial and fibroblast cells show regionalized expression by ventricle, starting at embryonic stages and persisting with age, with a dramatic transcriptional shift with maturation, and a smaller shift in each aged cell type. With aging, epithelial cells upregulate host-defense programs, and resident macrophages upregulate interleukin-1ß (IL-1ß) signaling genes. Our atlas reveals cellular diversity, architecture and signaling across ventricles during development, maturation, and aging of the ChP-brain barrier.

Stem-like intestinal Th17 cells give rise to pathogenic effector T cells during autoimmunity.

While intestinal Th17 cells are critical for maintaining tissue homeostasis, recent studies have implicated their roles in the development of extra-intestinal autoimmune diseases including multiple sclerosis. However, the mechanisms by which tissue Th17 cells mediate these dichotomous functions remain unknown. Here, we characterized the heterogeneity, plasticity, and migratory phenotypes of tissue Th17 cells in vivo by combined fate mapping with profiling of the transcriptomes and TCR clonotypes of over 84,000 Th17 cells at homeostasis and during CNS autoimmune inflammation. Inter- and intra-organ single-cell analyses revealed a homeostatic, stem-like TCF1+ IL-17+ SLAMF6+ population that traffics to the intestine where it is maintained by the microbiota, providing a ready reservoir for the IL-23-driven generation of encephalitogenic GM-CSF+ IFN-γ+ CXCR6+ T cells. Our study defines a direct in vivo relationship between IL-17+ non-pathogenic and GM-CSF+ and IFN-γ+ pathogenic Th17 populations and provides a mechanism by which homeostatic intestinal Th17 cells direct extra-intestinal autoimmune disease.

Metabolic modeling of single Th17 cells reveals regulators of autoimmunity.

Metabolism is a major regulator of immune cell function, but it remains difficult to study the metabolic status of individual cells. Here, we present Compass, an algorithm to characterize cellular metabolic states based on single-cell RNA sequencing and flux balance analysis. We applied Compass to associate metabolic states with T helper 17 (Th17) functional variability (pathogenic potential) and recovered a metabolic switch between glycolysis and fatty acid oxidation, akin to known Th17/regulatory T cell (Treg) differences, which we validated by metabolic assays. Compass also predicted that Th17 pathogenicity was associated with arginine and downstream polyamine metabolism. Indeed, polyamine-related enzyme expression was enhanced in pathogenic Th17 and suppressed in Treg cells. Chemical and genetic perturbation of polyamine metabolism inhibited Th17 cytokines, promoted Foxp3 expression, and remodeled the transcriptome and epigenome of Th17 cells toward a Treg-like state. In vivo perturbations of the polyamine pathway altered the phenotype of encephalitogenic T cells and attenuated tissue inflammation in CNS autoimmunity.

Antigen dominance hierarchies shape TCF1+ progenitor CD8 T cell phenotypes in tumors.

CD8 T cell responses against different tumor neoantigens occur simultaneously, yet little is known about the interplay between responses and its impact on T cell function and tumor control. In mouse lung adenocarcinoma, we found that immunodominance is established in tumors, wherein CD8 T cell expansion is predominantly driven by the antigen that most stably binds MHC. T cells responding to subdominant antigens were enriched for a TCF1+ progenitor phenotype correlated with response to immune checkpoint blockade (ICB) therapy. However, the subdominant T cell response did not preferentially benefit from ICB due to a dysfunctional subset of TCF1+ cells marked by CCR6 and Tc17 differentiation. Analysis of human samples and sequencing datasets revealed that CCR6+ TCF1+ cells exist across human cancers and are not correlated with ICB response. Vaccination eliminated CCR6+ TCF1+ cells and dramatically improved the subdominant response, highlighting a strategy to optimally engage concurrent neoantigen responses against tumors.

B cell genomics behind cross-neutralization of SARS-CoV-2 variants and SARS-CoV.

Monoclonal antibodies (mAbs) are a focus in vaccine and therapeutic design to counteract severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants. Here, we combined B cell sorting with single-cell VDJ and RNA sequencing (RNA-seq) and mAb structures to characterize B cell responses against SARS-CoV-2. We show that the SARS-CoV-2-specific B cell repertoire consists of transcriptionally distinct B cell populations with cells producing potently neutralizing antibodies (nAbs) localized in two clusters that resemble memory and activated B cells. Cryo-electron microscopy structures of selected nAbs from these two clusters complexed with SARS-CoV-2 spike trimers show recognition of various receptor-binding domain (RBD) epitopes. One of these mAbs, BG10-19, locks the spike trimer in a closed conformation to potently neutralize SARS-CoV-2, the recently arising mutants B.1.1.7 and B.1.351, and SARS-CoV and cross-reacts with heterologous RBDs. Together, our results characterize transcriptional differences among SARS-CoV-2-specific B cells and uncover cross-neutralizing Ab targets that will inform immunogen and therapeutic design against coronaviruses.

Inhibitory CD161 receptor identified in glioma-infiltrating T cells by single-cell analysis.

T cells are critical effectors of cancer immunotherapies, but little is known about their gene expression programs in diffuse gliomas. Here, we leverage single-cell RNA sequencing (RNA-seq) to chart the gene expression and clonal landscape of tumor-infiltrating T cells across 31 patients with isocitrate dehydrogenase (IDH) wild-type glioblastoma and IDH mutant glioma. We identify potential effectors of anti-tumor immunity in subsets of T cells that co-express cytotoxic programs and several natural killer (NK) cell genes. Analysis of clonally expanded tumor-infiltrating T cells further identifies the NK gene KLRB1 (encoding CD161) as a candidate inhibitory receptor. Accordingly, genetic inactivation of KLRB1 or antibody-mediated CD161 blockade enhances T cell-mediated killing of glioma cells in vitro and their anti-tumor function in vivo. KLRB1 and its associated transcriptional program are also expressed by substantial T cell populations in other human cancers. Our work provides an atlas of T cells in gliomas and highlights CD161 and other NK cell receptors as immunotherapy targets.

Differential pre-malignant programs and microenvironment chart distinct paths to malignancy in human colorectal polyps.

Colorectal cancers (CRCs) arise from precursor polyps whose cellular origins, molecular heterogeneity, and immunogenic potential may reveal diagnostic and therapeutic insights when analyzed at high resolution. We present a single-cell transcriptomic and imaging atlas of the two most common human colorectal polyps, conventional adenomas and serrated polyps, and their resulting CRC counterparts. Integrative analysis of 128 datasets from 62 participants reveals adenomas arise from WNT-driven expansion of stem cells, while serrated polyps derive from differentiated cells through gastric metaplasia. Metaplasia-associated damage is coupled to a cytotoxic immune microenvironment preceding hypermutation, driven partly by antigen-presentation differences associated with tumor cell-differentiation status. Microsatellite unstable CRCs contain distinct non-metaplastic regions where tumor cells acquire stem cell properties and cytotoxic immune cells are depleted. Our multi-omic atlas provides insights into malignant progression of colorectal polyps and their microenvironment, serving as a framework for precision surveillance and prevention of CRC.

Targeting PGLYRP1 promotes antitumor immunity while inhibiting autoimmune neuroinflammation.

Co-inhibitory and checkpoint molecules suppress T cell function in the tumor microenvironment, thereby rendering T cells dysfunctional. Although immune checkpoint blockade is a successful treatment option for multiple human cancers, severe autoimmune-like adverse effects can limit its application. Here, we show that the gene encoding peptidoglycan recognition protein 1 (PGLYRP1) is highly coexpressed with genes encoding co-inhibitory molecules, indicating that it might be a promising target for cancer immunotherapy. Genetic deletion of Pglyrp1 in mice led to decreased tumor growth and an increased activation/effector phenotype in CD8+ T cells, suggesting an inhibitory function of PGLYRP1 in CD8+ T cells. Surprisingly, genetic deletion of Pglyrp1 protected against the development of experimental autoimmune encephalomyelitis, a model of autoimmune disease in the central nervous system. PGLYRP1-deficient myeloid cells had a defect in antigen presentation and T cell activation, indicating that PGLYRP1 might function as a proinflammatory molecule in myeloid cells during autoimmunity. These results highlight PGLYRP1 as a promising target for immunotherapy that, when targeted, elicits a potent antitumor immune response while protecting against some forms of tissue inflammation and autoimmunity.

A genetic disorder reveals a hematopoietic stem cell regulatory network co-opted in leukemia.

The molecular regulation of human hematopoietic stem cell (HSC) maintenance is therapeutically important, but limitations in experimental systems and interspecies variation have constrained our knowledge of this process. Here, we have studied a rare genetic disorder due to MECOM haploinsufficiency, characterized by an early-onset absence of HSCs in vivo. By generating a faithful model of this disorder in primary human HSCs and coupling functional studies with integrative single-cell genomic analyses, we uncover a key transcriptional network involving hundreds of genes that is required for HSC maintenance. Through our analyses, we nominate cooperating transcriptional regulators and identify how MECOM prevents the CTCF-dependent genome reorganization that occurs as HSCs differentiate. We show that this transcriptional network is co-opted in high-risk leukemias, thereby enabling these cancers to acquire stem cell properties. Collectively, we illuminate a regulatory network necessary for HSC self-renewal through the study of a rare experiment of nature.

Temporal single-cell atlas of non-neuronal retinal cells reveals dynamic, coordinated multicellular responses to central nervous system injury.

Non-neuronal cells are key to the complex cellular interplay that follows central nervous system insult. To understand this interplay, we generated a single-cell atlas of immune, glial and retinal pigment epithelial cells from adult mouse retina before and at multiple time points after axonal transection. We identified rare subsets in naive retina, including interferon (IFN)-response glia and border-associated macrophages, and delineated injury-induced changes in cell composition, expression programs and interactions. Computational analysis charted a three-phase multicellular inflammatory cascade after injury. In the early phase, retinal macroglia and microglia were reactivated, providing chemotactic signals concurrent with infiltration of CCR2+ monocytes from the circulation. These cells differentiated into macrophages in the intermediate phase, while an IFN-response program, likely driven by microglia-derived type I IFN, was activated across resident glia. The late phase indicated inflammatory resolution. Our findings provide a framework to decipher cellular circuitry, spatial relationships and molecular interactions following tissue injury.

Chromatin Potential Identified by Shared Single-Cell Profiling of RNA and Chromatin.

Cell differentiation and function are regulated across multiple layers of gene regulation, including modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of regulatory events leading to cell fate commitment. Here we developed simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq), a highly scalable approach for measurement of chromatin accessibility and gene expression in the same single cell, applicable to different tissues. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define domains of regulatory chromatin (DORCs) that significantly overlap with super-enhancers. During lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for lineage commitment. We computationally infer chromatin potential as a quantitative measure of chromatin lineage-priming and use it to predict cell fate outcomes. SHARE-seq is an extensible platform to study regulatory circuitry across diverse cells in tissues.

The Human and Mouse Enteric Nervous System at Single-Cell Resolution.

The enteric nervous system (ENS) coordinates diverse functions in the intestine but has eluded comprehensive molecular characterization because of the rarity and diversity of cells. Here we develop two methods to profile the ENS of adult mice and humans at single-cell resolution: RAISIN RNA-seq for profiling intact nuclei with ribosome-bound mRNA and MIRACL-seq for label-free enrichment of rare cell types by droplet-based profiling. The 1,187,535 nuclei in our mouse atlas include 5,068 neurons from the ileum and colon, revealing extraordinary neuron diversity. We highlight circadian expression changes in enteric neurons, show that disease-related genes are dysregulated with aging, and identify differences between the ileum and proximal/distal colon. In humans, we profile 436,202 nuclei, recovering 1,445 neurons, and identify conserved and species-specific transcriptional programs and putative neuro-epithelial, neuro-stromal, and neuro-immune interactions. The human ENS expresses risk genes for neuropathic, inflammatory, and extra-intestinal diseases, suggesting neuronal contributions to disease.

The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous large-scale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.

DNA Microscopy: Optics-free Spatio-genetic Imaging by a Stand-Alone Chemical Reaction.

Analyzing the spatial organization of molecules in cells and tissues is a cornerstone of biological research and clinical practice. However, despite enormous progress in molecular profiling of cellular constituents, spatially mapping them remains a disjointed and specialized machinery-intensive process, relying on either light microscopy or direct physical registration. Here, we demonstrate DNA microscopy, a distinct imaging modality for scalable, optics-free mapping of relative biomolecule positions. In DNA microscopy of transcripts, transcript molecules are tagged in situ with randomized nucleotides, labeling each molecule uniquely. A second in situ reaction then amplifies the tagged molecules, concatenates the resulting copies, and adds new randomized nucleotides to uniquely label each concatenation event. An algorithm decodes molecular proximities from these concatenated sequences and infers physical images of the original transcripts at cellular resolution with precise sequence information. Because its imaging power derives entirely from diffusive molecular dynamics, DNA microscopy constitutes a chemically encoded microscopy system.

Molecular Classification and Comparative Taxonomics of Foveal and Peripheral Cells in Primate Retina.

High-acuity vision in primates, including humans, is mediated by a small central retinal region called the fovea. As more accessible organisms lack a fovea, its specialized function and its dysfunction in ocular diseases remain poorly understood. We used 165,000 single-cell RNA-seq profiles to generate comprehensive cellular taxonomies of macaque fovea and peripheral retina. More than 80% of >60 cell types match between the two regions but exhibit substantial differences in proportions and gene expression, some of which we relate to functional differences. Comparison of macaque retinal types with those of mice reveals that interneuron types are tightly conserved. In contrast, projection neuron types and programs diverge, despite exhibiting conserved transcription factor codes. Key macaque types are conserved in humans, allowing mapping of cell-type and region-specific expression of >190 genes associated with 7 human retinal diseases. Our work provides a framework for comparative single-cell analysis across tissue regions and species.