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.

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.

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.

Multiomic analysis reveals conservation of cancer-associated fibroblast phenotypes across species and tissue of origin.

Cancer-associated fibroblasts (CAFs) are integral to the solid tumor microenvironment. CAFs were once thought to be a relatively uniform population of matrix-producing cells, but single-cell RNA sequencing has revealed diverse CAF phenotypes. Here, we further probed CAF heterogeneity with a comprehensive multiomics approach. Using paired, same-cell chromatin accessibility and transcriptome analysis, we provided an integrated analysis of CAF subpopulations over a complex spatial transcriptomic and proteomic landscape to identify three superclusters: steady state-like (SSL), mechanoresponsive (MR), and immunomodulatory (IM) CAFs. These superclusters are recapitulated across multiple tissue types and species. Selective disruption of underlying mechanical force or immune checkpoint inhibition therapy results in shifts in CAF subpopulation distributions and affected tumor growth. As such, the balance among CAF superclusters may have considerable translational implications. Collectively, this research expands our understanding of CAF biology, identifying regulatory pathways in CAF differentiation and elucidating therapeutic targets in a species- and tumor-agnostic manner.

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.

The evolutionary dynamics of extrachromosomal DNA in human cancers.

Oncogene amplification on extrachromosomal DNA (ecDNA) is a common event, driving aggressive tumor growth, drug resistance and shorter survival. Currently, the impact of nonchromosomal oncogene inheritance-random identity by descent-is poorly understood. Also unclear is the impact of ecDNA on somatic variation and selection. Here integrating theoretical models of random segregation, unbiased image analysis, CRISPR-based ecDNA tagging with live-cell imaging and CRISPR-C, we demonstrate that random ecDNA inheritance results in extensive intratumoral ecDNA copy number heterogeneity and rapid adaptation to metabolic stress and targeted treatment. Observed ecDNAs benefit host cell survival or growth and can change within a single cell cycle. ecDNA inheritance can predict, a priori, some of the aggressive features of ecDNA-containing cancers. These properties are facilitated by the ability of ecDNA to rapidly adapt genomes in a way that is not possible through chromosomal oncogene amplification. These results show how the nonchromosomal random inheritance pattern of ecDNA contributes to poor outcomes for patients with cancer.

Genome-wide CRISPR screens of T cell exhaustion identify chromatin remodeling factors that limit T cell persistence.

T cell exhaustion limits antitumor immunity, but the molecular determinants of this process remain poorly understood. Using a chronic stimulation assay, we performed genome-wide CRISPR-Cas9 screens to systematically discover regulators of T cell exhaustion, which identified an enrichment of epigenetic factors. In vivo CRISPR screens in murine and human tumor models demonstrated that perturbation of the INO80 and BAF chromatin remodeling complexes improved T cell persistence in tumors. In vivo Perturb-seq revealed distinct transcriptional roles of each complex and that depletion of canonical BAF complex members, including Arid1a, resulted in the maintenance of an effector program and downregulation of exhaustion-related genes in tumor-infiltrating T cells. Finally, Arid1a depletion limited the acquisition of exhaustion-associated chromatin accessibility and led to improved antitumor immunity. In summary, we provide an atlas of the genetic regulators of T cell exhaustion and demonstrate that modulation of epigenetic state can improve T cell responses in cancer immunotherapy.

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.

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.

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.

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.

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.

Tracking the immune response with single-cell genomics.

The immune system is composed of a diverse array of cell types, each with a specialized role in orchestrating the immune response to pathogens or cancer. Even within a single cell 'type,' individual cells can access a wide spectrum of differentiation and activation states, which reflect the physiological response of each cell to the tissue environment and immune stimuli. Thus, the cellular diversity of the immune system is inherently quite complex and understanding this complexity has greatly benefited from technologies that measure immune responses at single-cell resolution, in addition to the systems-level response as a whole. In this Commentary, we focus on recent work at the interface of immunology and single-cell genomics and highlight advances in technologies and their application to immune cells. In particular, we highlight recent single-cell genomic profiling studies of T cells, since somatic rearrangements in the T cell receptor (TCR) loci enable the tracking of clonal T cell responses through space and time. Finally, we discuss opportunities for future use of these technologies in understanding vaccination and the basis for effective vaccine-induced immunity.

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion.

Understanding complex tissues requires single-cell deconstruction of gene regulation with precision and scale. Here, we assess the performance of a massively parallel droplet-based method for mapping transposase-accessible chromatin in single cells using sequencing (scATAC-seq). We apply scATAC-seq to obtain chromatin profiles of more than 200,000 single cells in human blood and basal cell carcinoma. In blood, application of scATAC-seq enables marker-free identification of cell type-specific cis- and trans-regulatory elements, mapping of disease-associated enhancer activity and reconstruction of trajectories of cellular differentiation. In basal cell carcinoma, application of scATAC-seq reveals regulatory networks in malignant, stromal and immune cells in the tumor microenvironment. Analysis of scATAC-seq profiles from serial tumor biopsies before and after programmed cell death protein 1 blockade identifies chromatin regulators of therapy-responsive T cell subsets and reveals a shared regulatory program that governs intratumoral CD8+ T cell exhaustion and CD4+ T follicular helper cell development. We anticipate that scATAC-seq will enable the unbiased discovery of gene regulatory factors across diverse biological systems.

Clonal replacement of tumor-specific T cells following PD-1 blockade.

Immunotherapies that block inhibitory checkpoint receptors on T cells have transformed the clinical care of patients with cancer1. However, whether the T cell response to checkpoint blockade relies on reinvigoration of pre-existing tumor-infiltrating lymphocytes or on recruitment of novel T cells remains unclear2-4. Here we performed paired single-cell RNA and T cell receptor sequencing on 79,046 cells from site-matched tumors from patients with basal or squamous cell carcinoma before and after anti-PD-1 therapy. Tracking T cell receptor clones and transcriptional phenotypes revealed coupling of tumor recognition, clonal expansion and T cell dysfunction marked by clonal expansion of CD8+CD39+ T cells, which co-expressed markers of chronic T cell activation and exhaustion. However, the expansion of T cell clones did not derive from pre-existing tumor-infiltrating T lymphocytes; instead, the expanded clones consisted of novel clonotypes that had not previously been observed in the same tumor. Clonal replacement of T cells was preferentially observed in exhausted CD8+ T cells and evident in patients with basal or squamous cell carcinoma. These results demonstrate that pre-existing tumor-specific T cells may have limited reinvigoration capacity, and that the T cell response to checkpoint blockade derives from a distinct repertoire of T cell clones that may have just recently entered the tumor.

Rapid and reversible suppression of ALT by DAXX in osteosarcoma cells.

Many tumors maintain chromosome-ends through a telomerase-independent, DNA-templated mechanism called alternative lengthening of telomeres (ALT). While ALT occurs in only a subset of tumors, it is strongly associated with mutations in the genes ATRX and DAXX, which encode components of an H3.3 histone chaperone complex. The role of ATRX and DAXX mutations in potentiating the mechanism of ALT remains incompletely understood. Here we characterize an osteosarcoma cell line, G292, with wild-type ATRX but a unique chromosome translocation resulting in loss of DAXX function. While ATRX and DAXX form a complex in G292, this complex fails to localize to nuclear PML bodies. We demonstrate that introduction of wild type DAXX suppresses the ALT phenotype and restores the localization of ATRX/DAXX to PML bodies. Using an inducible system, we show that ALT-associated PML bodies are disrupted rapidly following DAXX induction and that ALT is again restored following withdrawal of DAXX.

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.