Engineered cell entry links receptor biology with single-cell genomics.

Cells communicate with each other via receptor-ligand interactions. Here, we describe lentiviral-mediated cell entry by engineered receptor-ligand interaction (ENTER) to display ligand proteins, deliver payloads, and record receptor specificity. We optimize ENTER to decode interactions between T cell receptor (TCR)-MHC peptides, antibody-antigen, and other receptor-ligand pairs. A viral presentation strategy allows ENTER to capture interactions between B cell receptor and any antigen. We engineer ENTER to deliver genetic payloads to antigen-specific T or B cells to selectively modulate cellular behavior in mixed populations. Single-cell readout of ENTER by RNA sequencing (ENTER-seq) enables multiplexed enumeration of antigen specificities, TCR clonality, cell type, and states of individual T cells. ENTER-seq of CMV-seropositive patient blood samples reveals the viral epitopes that drive effector memory T cell differentiation and inter-clonal vs. intra-clonal phenotypic diversity targeting the same epitope. ENTER technology enables systematic discovery of receptor specificity, linkage to cell fates, and antigen-specific cargo delivery.

Single-Cell Analyses Identify Brain Mural Cells Expressing CD19 as Potential Off-Tumor Targets for CAR-T Immunotherapies.

CD19-directed immunotherapies are clinically effective for treating B cell malignancies but also cause a high incidence of neurotoxicity. A subset of patients treated with chimeric antigen receptor (CAR) T cells or bispecific T cell engager (BiTE) antibodies display severe neurotoxicity, including fatal cerebral edema associated with T cell infiltration into the brain. Here, we report that mural cells, which surround the endothelium and are critical for blood-brain-barrier integrity, express CD19. We identify CD19 expression in brain mural cells using single-cell RNA sequencing data and confirm perivascular staining at the protein level. CD19 expression in the brain begins early in development alongside the emergence of mural cell lineages and persists throughout adulthood across brain regions. Mouse mural cells demonstrate lower levels of Cd19 expression, suggesting limitations in preclinical animal models of neurotoxicity. These data suggest an on-target mechanism for neurotoxicity in CD19-directed therapies and highlight the utility of human single-cell atlases for designing immunotherapies.

Divergent clonal differentiation trajectories of T cell exhaustion.

Chronic antigen exposure during viral infection or cancer promotes an exhausted T cell (Tex) state with reduced effector function. However, whether all antigen-specific T cell clones follow the same Tex differentiation trajectory remains unclear. Here, we generate a single-cell multiomic atlas of T cell exhaustion in murine chronic viral infection that redefines Tex phenotypic diversity, including two late-stage Tex subsets with either a terminal exhaustion (Texterm) or a killer cell lectin-like receptor-expressing cytotoxic (TexKLR) phenotype. We use paired single-cell RNA and T cell receptor sequencing to uncover clonal differentiation trajectories of Texterm-biased, TexKLR-biased or divergent clones that acquire both phenotypes. We show that high T cell receptor signaling avidity correlates with Texterm, whereas low avidity correlates with effector-like TexKLR fate. Finally, we identify similar clonal differentiation trajectories in human tumor-infiltrating lymphocytes. These findings reveal clonal heterogeneity in the T cell response to chronic antigen that influences Tex fates and persistence.

Distinct epigenomic landscapes underlie tissue-specific memory T cell differentiation.

The memory CD8+ T cell pool contains phenotypically and transcriptionally heterogeneous subsets with specialized functions and recirculation patterns. Here, we examined the epigenetic landscape of CD8+ T cells isolated from seven non-lymphoid organs across four distinct infection models, alongside their circulating T cell counterparts. Using single-cell transposase-accessible chromatin sequencing (scATAC-seq), we found that tissue-resident memory T (TRM) cells and circulating memory T (TCIRC) cells develop along distinct epigenetic trajectories. We identified organ-specific transcriptional regulators of TRM cell development, including FOSB, FOS, FOSL1, and BACH2, and defined an epigenetic signature common to TRM cells across organs. Finally, we found that although terminal TEX cells share accessible regulatory elements with TRM cells, they are defined by TEX-specific epigenetic features absent from TRM cells. Together, this comprehensive data resource shows that TRM cell development is accompanied by dynamic transcriptome alterations and chromatin accessibility changes that direct tissue-adapted and functionally distinct T cell states.

Promoter of lncRNA Gene PVT1 Is a Tumor-Suppressor DNA Boundary Element.

Noncoding mutations in cancer genomes are frequent but challenging to interpret. PVT1 encodes an oncogenic lncRNA, but recurrent translocations and deletions in human cancers suggest alternative mechanisms. Here, we show that the PVT1 promoter has a tumor-suppressor function that is independent of PVT1 lncRNA. CRISPR interference of PVT1 promoter enhances breast cancer cell competition and growth in vivo. The promoters of the PVT1 and the MYC oncogenes, located 55 kb apart on chromosome 8q24, compete for engagement with four intragenic enhancers in the PVT1 locus, thereby allowing the PVT1 promoter to regulate pause release of MYC transcription. PVT1 undergoes developmentally regulated monoallelic expression, and the PVT1 promoter inhibits MYC expression only from the same chromosome via promoter competition. Cancer genome sequencing identifies recurrent mutations encompassing the human PVT1 promoter, and genome editing verified that PVT1 promoter mutation promotes cancer cell growth. These results highlight regulatory sequences of lncRNA genes as potential disease-associated DNA elements.

Impaired mitochondrial oxidative phosphorylation limits the self-renewal of T cells exposed to persistent antigen.

The majority of tumor-infiltrating T cells exhibit a terminally exhausted phenotype, marked by a loss of self-renewal capacity. How repetitive antigenic stimulation impairs T cell self-renewal remains poorly defined. Here, we show that persistent antigenic stimulation impaired ADP-coupled oxidative phosphorylation. The resultant bioenergetic compromise blocked proliferation by limiting nucleotide triphosphate synthesis. Inhibition of mitochondrial oxidative phosphorylation in activated T cells was sufficient to suppress proliferation and upregulate genes linked to T cell exhaustion. Conversely, prevention of mitochondrial oxidative stress during chronic T cell stimulation allowed sustained T cell proliferation and induced genes associated with stem-like progenitor T cells. As a result, antioxidant treatment enhanced the anti-tumor efficacy of chronically stimulated T cells. These data reveal that loss of ATP production through oxidative phosphorylation limits T cell proliferation and effector function during chronic antigenic stimulation. Furthermore, treatments that maintain redox balance promote T cell self-renewal and enhance anti-tumor immunity.

Deciphering cell states and genealogies of human haematopoiesis.

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

A human autoimmune organoid model reveals IL-7 function in coeliac disease.

In vitro models of autoimmunity are constrained by an inability to culture affected epithelium alongside the complex tissue-resident immune microenvironment. Coeliac disease (CeD) is an autoimmune disease in which dietary gluten-derived peptides bind to the major histocompatibility complex (MHC) class II human leukocyte antigen molecules (HLA)-DQ2 or HLA-DQ8 to initiate immune-mediated duodenal mucosal injury1-4. Here, we generated air-liquid interface (ALI) duodenal organoids from intact fragments of endoscopic biopsies that preserve epithelium alongside native mesenchyme and tissue-resident immune cells as a unit without requiring reconstitution. The immune diversity of ALI organoids spanned T cells, B and plasma cells, natural killer (NK) cells and myeloid cells, with extensive T-cell and B-cell receptor repertoires. HLA-DQ2.5-restricted gluten peptides selectively instigated epithelial destruction in HLA-DQ2.5-expressing organoids derived from CeD patients, and this was antagonized by blocking MHC-II or NKG2C/D. Gluten epitopes stimulated a CeD organoid immune network response in lymphoid and myeloid subsets alongside anti-transglutaminase 2 (TG2) autoantibody production. Functional studies in CeD organoids revealed that interleukin-7 (IL-7) is a gluten-inducible pathogenic modulator that regulates CD8+ T-cell NKG2C/D expression and is necessary and sufficient for epithelial destruction. Furthermore, endogenous IL-7 was markedly upregulated in patient biopsies from active CeD compared with remission disease from gluten-free diets, predominantly in lamina propria mesenchyme. By preserving the epithelium alongside diverse immune populations, this human in vitro CeD model recapitulates gluten-dependent pathology, enables mechanistic investigation and establishes a proof of principle for the organoid modelling of autoimmunity.

ecDNA hubs drive cooperative intermolecular oncogene expression.

Extrachromosomal DNA (ecDNA) is prevalent in human cancers and mediates high expression of oncogenes through gene amplification and altered gene regulation1. Gene induction typically involves cis-regulatory elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs-clusters of around 10-100 ecDNAs within the nucleus-enable intermolecular enhancer-gene interactions to promote oncogene overexpression. ecDNAs that encode multiple distinct oncogenes form hubs in diverse cancer cell types and primary tumours. Each ecDNA is more likely to transcribe the oncogene when spatially clustered with additional ecDNAs. ecDNA hubs are tethered by the bromodomain and extraterminal domain (BET) protein BRD4 in a MYC-amplified colorectal cancer cell line. The BET inhibitor JQ1 disperses ecDNA hubs and preferentially inhibits ecDNA-derived-oncogene transcription. The BRD4-bound PVT1 promoter is ectopically fused to MYC and duplicated in ecDNA, receiving promiscuous enhancer input to drive potent expression of MYC. Furthermore, the PVT1 promoter on an exogenous episome suffices to mediate gene activation in trans by ecDNA hubs in a JQ1-sensitive manner. Systematic silencing of ecDNA enhancers by CRISPR interference reveals intermolecular enhancer-gene activation among multiple oncogene loci that are amplified on distinct ecDNAs. Thus, protein-tethered ecDNA hubs enable intermolecular transcriptional regulation and may serve as units of oncogene function and cooperative evolution and as potential targets for cancer therapy.

Polycomb-mediated genome architecture enables long-range spreading of H3K27 methylation.

Polycomb-group proteins play critical roles in gene silencing through the deposition of histone H3 lysine 27 trimethylation (H3K27me3) and chromatin compaction. This process is essential for embryonic stem cell (ESC) pluripotency, differentiation, and development. Polycomb repressive complex 2 (PRC2) can both read and write H3K27me3, enabling progressive spreading of H3K27me3 on the linear genome. Long-range Polycomb-associated DNA contacts have also been described, but their regulation and role in gene silencing remain unclear. Here, we apply H3K27me3 HiChIP, a protein-directed chromosome conformation method, and optical reconstruction of chromatin architecture to profile long-range Polycomb-associated DNA loops that span tens to hundreds of megabases across multiple topological associated domains in mouse ESCs and human induced pluripotent stem cells. We find that H3K27me3 loop anchors are enriched for Polycomb nucleation points and coincide with key developmental genes. Genetic deletion of H3K27me3 loop anchors results in disruption of spatial contact between distant loci and altered H3K27me3 in cis, both locally and megabases away on the same chromosome. In mouse embryos, loop anchor deletion leads to ectopic activation of the partner gene, suggesting that Polycomb-associated loops control gene silencing during development. Further, we find that alterations in PRC2 occupancy resulting from an RNA binding­deficient EZH2 mutant are accompanied by loss of Polycomb-associated DNA looping. Together, these results suggest PRC2 uses RNA binding to enhance long-range chromosome folding and H3K27me3 spreading. Developmental gene loci have unique roles in Polycomb spreading, emerging as important architectural elements of the epigenome.

BABEL enables cross-modality translation between multiomic profiles at single-cell resolution.

Simultaneous profiling of multiomic modalities within a single cell is a grand challenge for single-cell biology. While there have been impressive technical innovations demonstrating feasibility-for example, generating paired measurements of single-cell transcriptome (single-cell RNA sequencing [scRNA-seq]) and chromatin accessibility (single-cell assay for transposase-accessible chromatin using sequencing [scATAC-seq])-widespread application of joint profiling is challenging due to its experimental complexity, noise, and cost. Here, we introduce BABEL, a deep learning method that translates between the transcriptome and chromatin profiles of a single cell. Leveraging an interoperable neural network model, BABEL can predict single-cell expression directly from a cell's scATAC-seq and vice versa after training on relevant data. This makes it possible to computationally synthesize paired multiomic measurements when only one modality is experimentally available. Across several paired single-cell ATAC and gene expression datasets in human and mouse, we validate that BABEL accurately translates between these modalities for individual cells. BABEL also generalizes well to cell types within new biological contexts not seen during training. Starting from scATAC-seq of patient-derived basal cell carcinoma (BCC), BABEL generated single-cell expression that enabled fine-grained classification of complex cell states, despite having never seen BCC data. These predictions are comparable to analyses of experimental BCC scRNA-seq data for diverse cell types related to BABEL's training data. We further show that BABEL can incorporate additional single-cell data modalities, such as protein epitope profiling, thus enabling translation across chromatin, RNA, and protein. BABEL offers a powerful approach for data exploration and hypothesis generation.

Integrated spatial multiomics reveals fibroblast fate during tissue repair.

In the skin, tissue injury results in fibrosis in the form of scars composed of dense extracellular matrix deposited by fibroblasts. The therapeutic goal of regenerative wound healing has remained elusive, in part because principles of fibroblast programming and adaptive response to injury remain incompletely understood. Here, we present a multimodal -omics platform for the comprehensive study of cell populations in complex tissue, which has allowed us to characterize the cells involved in wound healing across both time and space. We employ a stented wound model that recapitulates human tissue repair kinetics and multiple Rainbow transgenic lines to precisely track fibroblast fate during the physiologic response to skin injury. Through integrated analysis of single cell chromatin landscapes and gene expression states, coupled with spatial transcriptomic profiling, we are able to impute fibroblast epigenomes with temporospatial resolution. This has allowed us to reveal potential mechanisms controlling fibroblast fate during migration, proliferation, and differentiation following skin injury, and thereby reexamine the canonical phases of wound healing. These findings have broad implications for the study of tissue repair in complex organ systems.

Invariant natural killer T-cell subsets have diverse graft-versus-host-disease-preventing and antitumor effects.

Invariant natural killer T (iNKT) cells are a T-cell subset with potent immunomodulatory properties. Experimental evidence in mice and observational studies in humans indicate that iNKT cells have antitumor potential as well as the ability to suppress acute and chronic graft-versus-host-disease (GVHD). Murine iNKT cells differentiate during thymic development into iNKT1, iNKT2, and iNKT17 sublineages, which differ transcriptomically and epigenomically and have subset-specific developmental requirements. Whether distinct iNKT sublineages also differ in their antitumor effect and their ability to suppress GVHD is currently unknown. In this work, we generated highly purified murine iNKT sublineages, characterized their transcriptomic and epigenomic landscape, and assessed specific functions. We show that iNKT2 and iNKT17, but not iNKT1, cells efficiently suppress T-cell activation in vitro and mitigate murine acute GVHD in vivo. Conversely, we show that iNKT1 cells display the highest antitumor activity against murine B-cell lymphoma cells both in vitro and in vivo. Thus, we report for the first time that iNKT sublineages have distinct and different functions, with iNKT1 cells having the highest antitumor activity and iNKT2 and iNKT17 cells having immune-regulatory properties. These results have important implications for the translation of iNKT cell therapies to the clinic for cancer immunotherapy as well as for the prevention and treatment of GVHD.

Identification of a T-bethi Quiescent Exhausted CD8 T Cell Subpopulation That Can Differentiate into TIM3+CX3CR1+ Effectors and Memory-like Cells.

Persistent Ag induces a dysfunctional CD8 T cell state known as "exhaustion" characterized by PD-1 expression. Nevertheless, exhausted CD8 T cells retain functionality through continued differentiation of progenitor into effector cells. However, it remains ill-defined how CD8 T cell effector responses are sustained in situ. In this study, we show using the mouse chronic lymphocytic choriomeningitis virus infection model that CX3CR1+ CD8 T cells contain a T-bet-dependent TIM3-PD-1lo subpopulation that is distinct from the TIM3+CX3CR1+PD-1+ proliferative effector subset. The TIM3-CX3CR1+ cells are quiescent and express a low but significant level of the transcription factor TCF-1, demonstrating similarity to TCF-1hi progenitor CD8 T cells. Furthermore, following the resolution of lymphocytic choriomeningitis virus viremia, a substantial proportion of TCF-1+ memory-like CD8 T cells show evidence of CX3CR1 expression during the chronic phase of the infection. Our results suggest a subset of the CX3CR1+ exhausted population demonstrates progenitor-like features that support the generation of the CX3CR1+ effector pool from the TCF-1hi progenitors and contribute to the memory-like pool following the resolution of viremia.

Regulation of immune signal integration and memory by inflammation-induced chromosome conformation.

3-dimensional (3D) genome conformation is central to gene expression regulation, yet our understanding of its contribution to rapid transcriptional responses, signal integration, and memory in immune cells is limited. Here, we study the molecular regulation of the inflammatory response in primary macrophages using integrated transcriptomic, epigenomic, and chromosome conformation data, including base pair-resolution Micro-Capture C. We demonstrate that interleukin-4 (IL-4) primes the inflammatory response in macrophages by stably rewiring 3D genome conformation, juxtaposing endotoxin-, interferon-gamma-, and dexamethasone-responsive enhancers in close proximity to their cognate gene promoters. CRISPR-based perturbations of enhancer-promoter contacts or CCCTC-binding factor (CTCF) boundary elements demonstrated that IL-4-driven conformation changes are indispensable for enhanced and synergistic endotoxin-induced transcriptional responses, as well as transcriptional memory following stimulus removal. Moreover, transcriptional memory mediated by changes in chromosome conformation often occurred in the absence of changes in chromatin accessibility or histone modifications. Collectively, these findings demonstrate that rapid and memory transcriptional responses to immunological stimuli are encoded in the 3D genome.

Disparate pathways for extrachromosomal DNA biogenesis and genomic DNA repair.

Oncogene amplification on extrachromosomal DNA (ecDNA) is a pervasive driver event in cancer, yet our understanding of how ecDNA forms is limited. Here, we couple a CRISPR-based method for ecDNA induction with extensive characterization of newly formed ecDNA to examine their biogenesis. We find that DNA circularization is efficient, irrespective of 3D genome context, with formation of 800kb, 1 Mb, and 1.8 Mb ecDNAs reaching or exceeding 15%. We show non-homologous end joining and microhomology-mediated end joining both contribute to ecDNA formation, while inhibition of DNA-PKcs and ATM have opposing impacts on ecDNA formation. EcDNA and the corresponding chromosomal excision scar can form at significantly different rates and respond differently to DNA-PKcs and ATM inhibition. Taken together, our results support a model of ecDNA formation in which double strand break ends dissociate from their legitimate ligation partners prior to joining of illegitimate ends to form the ecDNA and excision scar.

Disparate pathways for extrachromosomal DNA biogenesis and genomic DNA repair.

Oncogene amplification on extrachromosomal DNA (ecDNA) is a pervasive driver event in cancer, yet our understanding of how ecDNA forms is limited. Here, we couple a CRISPR-based method for induction of ecDNA with extensive characterization of newly formed ecDNA to examine ecDNA biogenesis. We find that DNA circularization is efficient, irrespective of 3D genome context, with formation of a 1 Mb and 1.8 Mb ecDNA both reaching 15%. We show non-homologous end joining and microhomology mediated end joining both contribute to ecDNA formation, while inhibition of DNA-PKcs and ATM have opposing impacts on ecDNA formation. EcDNA and the corresponding chromosomal excision scar form at significantly different rates and respond differently to DNA-PKcs and ATM inhibition. Taken together, our results support a model of ecDNA formation in which double strand break ends dissociate from their legitimate ligation partners prior to joining of illegitimate ends to form the ecDNA and excision scar.

RNA-binding proteins direct myogenic cell fate decisions.

RNA-binding proteins (RBPs), essential for skeletal muscle regeneration, cause muscle degeneration and neuromuscular disease when mutated. Why mutations in these ubiquitously expressed RBPs orchestrate complex tissue regeneration and direct cell fate decisions in skeletal muscle remains poorly understood. Single-cell RNA-sequencing of regenerating Mus musculus skeletal muscle reveals that RBP expression, including the expression of many neuromuscular disease-associated RBPs, is temporally regulated in skeletal muscle stem cells and correlates with specific stages of myogenic differentiation. By combining machine learning with RBP engagement scoring, we discovered that the neuromuscular disease-associated RBP Hnrnpa2b1 is a differentiation-specifying regulator of myogenesis that controls myogenic cell fate transitions during terminal differentiation in mice. The timing of RBP expression specifies cell fate transitions by providing post-transcriptional regulation of messenger RNAs that coordinate stem cell fate decisions during tissue regeneration.