Single-cell chromatin accessibility reveals malignant regulatory programs in primary human cancers.

To identify cancer-associated gene regulatory changes, we generated single-cell chromatin accessibility landscapes across eight tumor types as part of The Cancer Genome Atlas. Tumor chromatin accessibility is strongly influenced by copy number alterations that can be used to identify subclones, yet underlying cis-regulatory landscapes retain cancer type-specific features. Using organ-matched healthy tissues, we identified the "nearest healthy" cell types in diverse cancers, demonstrating that the chromatin signature of basal-like-subtype breast cancer is most similar to secretory-type luminal epithelial cells. Neural network models trained to learn regulatory programs in cancer revealed enrichment of model-prioritized somatic noncoding mutations near cancer-associated genes, suggesting that dispersed, nonrecurrent, noncoding mutations in cancer are functional. Overall, these data and interpretable gene regulatory models for cancer and healthy tissue provide a framework for understanding cancer-specific gene regulation.

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.

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.

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.

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.

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.

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.

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.

Reinfection with SARS-CoV-2 and Waning Humoral Immunity: A Case Report.

Recovery from COVID-19 is associated with production of anti-SARS-CoV-2 antibodies, but it is uncertain whether these confer immunity. We describe viral RNA shedding duration in hospitalized patients and identify patients with recurrent shedding. We sequenced viruses from two distinct episodes of symptomatic COVID-19 separated by 144 days in a single patient, to conclusively describe reinfection with a different strain harboring the spike variant D614G. This case of reinfection was one of the first cases of reinfection reported in 2020. With antibody, B cell and T cell analytics, we show correlates of adaptive immunity at reinfection, including a differential response in neutralizing antibodies to a D614G pseudovirus. Finally, we discuss implications for vaccine programs and begin to define benchmarks for protection against reinfection from SARS-CoV-2.

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.

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.

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.

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.