A CTCF-binding site in the Mdm1-Il22-Ifng locus shapes cytokine expression profiles and plays a critical role in early Th1 cell fate specification.

Cytokine expression during T cell differentiation is a highly regulated process that involves long-range promoter-enhancer and CTCF-CTCF contacts at cytokine loci. Here, we investigated the impact of dynamic chromatin loop formation within the topologically associating domain (TAD) in regulating the expression of interferon gamma (IFN-γ) and interleukin-22 (IL-22); these cytokine loci are closely located in the genome and are associated with complex enhancer landscapes, which are selectively active in type 1 and type 3 lymphocytes. In situ Hi-C analyses revealed inducible TADs that insulated Ifng and Il22 enhancers during Th1 cell differentiation. Targeted deletion of a 17 bp boundary motif of these TADs imbalanced Th1- and Th17-associated immunity, both in vitro and in vivo, upon Toxoplasma gondii infection. In contrast, this boundary element was dispensable for cytokine regulation in natural killer cells. Our findings suggest that precise cytokine regulation relies on lineage- and developmental stage-specific interactions of 3D chromatin architectures and enhancer landscapes.

Functional dissection of complex and molecular trait variants at single nucleotide resolution.

Identifying the causal variants and mechanisms that drive complex traits and diseases remains a core problem in human genetics. The majority of these variants have individually weak effects and lie in non-coding gene-regulatory elements where we lack a complete understanding of how single nucleotide alterations modulate transcriptional processes to affect human phenotypes. To address this, we measured the activity of 221,412 trait-associated variants that had been statistically fine-mapped using a Massively Parallel Reporter Assay (MPRA) in 5 diverse cell-types. We show that MPRA is able to discriminate between likely causal variants and controls, identifying 12,025 regulatory variants with high precision. Although the effects of these variants largely agree with orthogonal measures of function, only 69% can plausibly be explained by the disruption of a known transcription factor (TF) binding motif. We dissect the mechanisms of 136 variants using saturation mutagenesis and assign impacted TFs for 91% of variants without a clear canonical mechanism. Finally, we provide evidence that epistasis is prevalent for variants in close proximity and identify multiple functional variants on the same haplotype at a small, but important, subset of trait-associated loci. Overall, our study provides a systematic functional characterization of likely causal common variants underlying complex and molecular human traits, enabling new insights into the regulatory grammar underlying disease risk.

FOXO1 is a master regulator of memory programming in CAR T cells.

A major limitation of chimeric antigen receptor (CAR) T cell therapies is the poor persistence of these cells in vivo1. The expression of memory-associated genes in CAR T cells is linked to their long-term persistence in patients and clinical efficacy2-6, suggesting that memory programs may underpin durable CAR T cell function. Here we show that the transcription factor FOXO1 is responsible for promoting memory and restraining exhaustion in human CAR T cells. Pharmacological inhibition or gene editing of endogenous FOXO1 diminished the expression of memory-associated genes, promoted an exhaustion-like phenotype and impaired the antitumour activity of CAR T cells. Overexpression of FOXO1 induced a gene-expression program consistent with T cell memory and increased chromatin accessibility at FOXO1-binding motifs. CAR T cells that overexpressed FOXO1 retained their function, memory potential and metabolic fitness in settings of chronic stimulation, and exhibited enhanced persistence and tumour control in vivo. By contrast, overexpression of TCF1 (encoded by TCF7) did not enforce canonical memory programs or enhance the potency of CAR T cells. Notably, FOXO1 activity correlated with positive clinical outcomes of patients treated with CAR T cells or tumour-infiltrating lymphocytes, underscoring the clinical relevance of FOXO1 in cancer immunotherapy. Our results show that overexpressing FOXO1 can increase the antitumour activity of human CAR T cells, and highlight memory reprogramming as a broadly applicable approach for optimizing therapeutic T cell states.

Convergent epigenetic evolution drives relapse in acute myeloid leukemia.

Relapse of acute myeloid leukemia (AML) is highly aggressive and often treatment refractory. We analyzed previously published AML relapse cohorts and found that 40% of relapses occur without changes in driver mutations, suggesting that non-genetic mechanisms drive relapse in a large proportion of cases. We therefore characterized epigenetic patterns of AML relapse using 26 matched diagnosis-relapse samples with ATAC-seq. This analysis identified a relapse-specific chromatin accessibility signature for mutationally stable AML, suggesting that AML undergoes epigenetic evolution at relapse independent of mutational changes. Analysis of leukemia stem cell (LSC) chromatin changes at relapse indicated that this leukemic compartment underwent significantly less epigenetic evolution than non-LSCs, while epigenetic changes in non-LSCs reflected overall evolution of the bulk leukemia. Finally, we used single-cell ATAC-seq paired with mitochondrial sequencing (mtscATAC) to map clones from diagnosis into relapse along with their epigenetic features. We found that distinct mitochondrially-defined clones exhibit more similar chromatin accessibility at relapse relative to diagnosis, demonstrating convergent epigenetic evolution in relapsed AML. These results demonstrate that epigenetic evolution is a feature of relapsed AML and that convergent epigenetic evolution can occur following treatment with induction chemotherapy.

Acute myeloid leukemia (or AML for short) is a type of blood cancer characterized by abnormally high production of immature white blood cells. Despite advances in AML treatment, many patients relapse after an initially successful first round of treatment. As a result, understanding the factors contributing to relapse is essential for developing effective treatments for the disease. Like most cancers, AML can evolve because of changes to the DNA sequence in cells that cause them to grow uncontrollably or resist treatment. Alongside these genetic mutations, AML cells also undergo 'epigenetic' changes, where regions of the DNA are modified and genes can be switched on or off without altering the DNA sequence. Previous research has demonstrated that epigenetic changes contribute to the development of AML, however, it was not clear if these changes could also make cells resistant to treatment without acquiring new DNA mutations. Nuno, Azizi et al. addressed this question by analyzing the epigenetic states of AML cells from 26 patients at the time of their diagnosis and after treatment when the disease had relapsed. Analysis revealed that almost half of the patients with AML experienced a relapse without acquiring new DNA mutations. Instead, these AML cells developed specific epigenetic changes that helped them to resist cancer treatment. Moreover, studying individual AML cells from different patients showed that the cells became more epigenetically similar at relapse, suggesting that they converge towards a more treatment-resistant disease. Future experiments will determine exactly how these epigenetic changes lead to treatment resistance. Currently, most of the drugs used to treat AML are either chemotherapies or ones that target specific DNA mutations. The findings of Nuno, Azizi et al. suggest that drugs targeting specific epigenetic changes may be more effective for some patients. Further studies will be needed to determine which patients may benefit and which epigenetic drugs could be useful.

Transcript-specific enrichment enables profiling rare cell states via scRNA-seq.

Single-cell genomics technologies have accelerated our understanding of cell-state heterogeneity in diverse contexts. Although single-cell RNA sequencing (scRNA-seq) identifies many rare populations of interest that express specific marker transcript combinations, traditional flow sorting limits our ability to enrich these populations for further profiling, including requiring cell surface markers with high-fidelity antibodies. Additionally, many single-cell studies require the isolation of nuclei from tissue, eliminating the ability to enrich learned rare cell states based on extranuclear protein markers. To address these limitations, we describe Programmable Enrichment via RNA Flow-FISH by sequencing (PERFF-seq), a scalable assay that enables scRNA-seq profiling of subpopulations from complex cellular mixtures defined by the presence or absence of specific RNA transcripts. Across immune populations (n = 141,227 cells) and fresh-frozen and formalin-fixed paraffin-embedded brain tissue (n = 29,522 nuclei), we demonstrate the sorting logic that can be used to enrich for cell populations via RNA-based cytometry followed by high-throughput scRNA-seq. Our approach provides a rational, programmable method for studying rare populations identified by one or more marker transcripts.

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.

Multi-modal skin atlas identifies a multicellular immune-stromal community associated with altered cornification and specific T cell expansion in atopic dermatitis.

In healthy skin, a cutaneous immune system maintains the balance between tolerance towards innocuous environmental antigens and immune responses against pathological agents. In atopic dermatitis (AD), barrier and immune dysfunction result in chronic tissue inflammation. Our understanding of the skin tissue ecosystem in AD remains incomplete with regard to the hallmarks of pathological barrier formation, and cellular state and clonal composition of disease-promoting cells. Here, we generated a multi-modal cell census of 310,691 cells spanning 86 cell subsets from whole skin tissue of 19 adult individuals, including non-lesional and lesional skin from 11 AD patients, and integrated it with 396,321 cells from four studies into a comprehensive human skin cell atlas in health and disease. Reconstruction of human keratinocyte differentiation from basal to cornified layers revealed a disrupted cornification trajectory in AD. This disrupted epithelial differentiation was associated with signals from a unique immune and stromal multicellular community comprised of MMP12 + dendritic cells (DCs), mature migratory DCs, cycling ILCs, NK cells, inflammatory CCL19 + IL4I1 + fibroblasts, and clonally expanded IL13 + IL22 + IL26 + T cells with overlapping type 2 and type 17 characteristics. Cell subsets within this immune and stromal multicellular community were connected by multiple inter-cellular positive feedback loops predicted to impact community assembly and maintenance. AD GWAS gene expression was enriched both in disrupted cornified keratinocytes and in cell subsets from the lesional immune and stromal multicellular community including IL13 + IL22 + IL26 + T cells and ILCs, suggesting that epithelial or immune dysfunction in the context of the observed cellular communication network can initiate and then converge towards AD. Our work highlights specific, disease-associated cell subsets and interactions as potential targets in progression and resolution of chronic inflammation.

FOXO1 is a master regulator of CAR T memory programming.

Poor CAR T persistence limits CAR T cell therapies for B cell malignancies and solid tumors1,2. The expression of memory-associated genes such as TCF7 (protein name TCF1) is linked to response and long-term persistence in patients3-7, thereby implicating memory programs in therapeutic efficacy. Here, we demonstrate that the pioneer transcription factor, FOXO1, is responsible for promoting memory programs and restraining exhaustion in human CAR T cells. Pharmacologic inhibition or gene editing of endogenous FOXO1 in human CAR T cells diminished the expression of memory-associated genes, promoted an exhaustion-like phenotype, and impaired antitumor activity in vitro and in vivo. FOXO1 overexpression induced a gene expression program consistent with T cell memory and increased chromatin accessibility at FOXO1 binding motifs. FOXO1-overexpressing cells retained function, memory potential, and metabolic fitness during settings of chronic stimulation and exhibited enhanced persistence and antitumor activity in vivo. In contrast, TCF1 overexpression failed to enforce canonical memory programs or enhance CAR T cell potency. Importantly, endogenous FOXO1 activity correlated with CAR T and TIL responses in patients, underscoring its clinical relevance in cancer immunotherapy. Our results demonstrate that memory reprogramming through FOXO1 can enhance the persistence and potency of human CAR T cells and highlights the utility of pioneer factors, which bind condensed chromatin and induce local epigenetic remodeling, for optimizing therapeutic T cell states.

Latent human herpesvirus 6 is reactivated in CAR T cells.

Cell therapies have yielded durable clinical benefits for patients with cancer, but the risks associated with the development of therapies from manipulated human cells are understudied. For example, we lack a comprehensive understanding of the mechanisms of toxicities observed in patients receiving T cell therapies, including recent reports of encephalitis caused by reactivation of human herpesvirus 6 (HHV-6)1. Here, through petabase-scale viral genomics mining, we examine the landscape of human latent viral reactivation and demonstrate that HHV-6B can become reactivated in cultures of human CD4+ T cells. Using single-cell sequencing, we identify a rare population of HHV-6 'super-expressors' (about 1 in 300-10,000 cells) that possess high viral transcriptional activity, among research-grade allogeneic chimeric antigen receptor (CAR) T cells. By analysing single-cell sequencing data from patients receiving cell therapy products that are approved by the US Food and Drug Administration2 or are in clinical studies3-5, we identify the presence of HHV-6-super-expressor CAR T cells in patients in vivo. Together, the findings of our study demonstrate the utility of comprehensive genomics analyses in implicating cell therapy products as a potential source contributing to the lytic HHV-6 infection that has been reported in clinical trials1,6-8 and may influence the design and production of autologous and allogeneic cell therapies.

SND1 binds SARS-CoV-2 negative-sense RNA and promotes viral RNA synthesis through NSP9.

Regulation of viral RNA biogenesis is fundamental to productive SARS-CoV-2 infection. To characterize host RNA-binding proteins (RBPs) involved in this process, we biochemically identified proteins bound to genomic and subgenomic SARS-CoV-2 RNAs. We find that the host protein SND1 binds the 5' end of negative-sense viral RNA and is required for SARS-CoV-2 RNA synthesis. SND1-depleted cells form smaller replication organelles and display diminished virus growth kinetics. We discover that NSP9, a viral RBP and direct SND1 interaction partner, is covalently linked to the 5' ends of positive- and negative-sense RNAs produced during infection. These linkages occur at replication-transcription initiation sites, consistent with NSP9 priming viral RNA synthesis. Mechanistically, SND1 remodels NSP9 occupancy and alters the covalent linkage of NSP9 to initiating nucleotides in viral RNA. Our findings implicate NSP9 in the initiation of SARS-CoV-2 RNA synthesis and unravel an unsuspected role of a cellular protein in orchestrating viral RNA production.

Convergent Epigenetic Evolution Drives Relapse in Acute Myeloid Leukemia.

Relapse of acute myeloid leukemia (AML) is highly aggressive and often treatment refractory. We analyzed previously published AML relapse cohorts and found that 40% of relapses occur without changes in driver mutations, suggesting that non-genetic mechanisms drive relapse in a large proportion of cases. We therefore characterized epigenetic patterns of AML relapse using 26 matched diagnosis-relapse samples with ATAC-seq. This analysis identified a relapse-specific chromatin accessibility signature for mutationally stable AML, suggesting that AML undergoes epigenetic evolution at relapse independent of mutational changes. Analysis of leukemia stem cell (LSC) chromatin changes at relapse indicated that this leukemic compartment underwent significantly less epigenetic evolution than non-LSCs, while epigenetic changes in non-LSCs reflected overall evolution of the bulk leukemia. Finally, we used single-cell ATAC-seq paired with mitochondrial sequencing (mtscATAC) to map clones from diagnosis into relapse along with their epigenetic features. We found that distinct mitochondrially-defined clones exhibit more similar chromatin accessibility at relapse relative to diagnosis, demonstrating convergent epigenetic evolution in relapsed AML. These results demonstrate that epigenetic evolution is a feature of relapsed AML and that convergent epigenetic evolution can occur following treatment with induction chemotherapy.

Single-cell multi-omics of mitochondrial DNA disorders reveals dynamics of purifying selection across human immune cells.

Pathogenic mutations in mitochondrial DNA (mtDNA) compromise cellular metabolism, contributing to cellular heterogeneity and disease. Diverse mutations are associated with diverse clinical phenotypes, suggesting distinct organ- and cell-type-specific metabolic vulnerabilities. Here we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from six patients across the phenotypic spectrum of single large-scale mtDNA deletions (SLSMDs). By profiling 206,663 cells, we reveal the dynamics of pathogenic mtDNA deletion heteroplasmy consistent with purifying selection and distinct metabolic vulnerabilities across T-cell states in vivo and validate these observations in vitro. By extending analyses to hematopoietic and erythroid progenitors, we reveal mtDNA dynamics and cell-type-specific gene regulatory adaptations, demonstrating the context-dependence of perturbing mitochondrial genomic integrity. Collectively, we report pathogenic mtDNA heteroplasmy dynamics of individual blood and immune cells across lineages, demonstrating the power of single-cell multi-omics for revealing fundamental properties of mitochondrial genetics.

Co-opting signalling molecules enables logic-gated control of CAR T cells.

Although chimeric antigen receptor (CAR) T cells have altered the treatment landscape for B cell malignancies, the risk of on-target, off-tumour toxicity has hampered their development for solid tumours because most target antigens are shared with normal cells1,2. Researchers have attempted to apply Boolean-logic gating to CAR T cells to prevent toxicity3-5; however, a truly safe and effective logic-gated CAR has remained elusive6. Here we describe an approach to CAR engineering in which we replace traditional CD3ζ domains with intracellular proximal T cell signalling molecules. We show that certain proximal signalling CARs, such as a ZAP-70 CAR, can activate T cells and eradicate tumours in vivo while bypassing upstream signalling proteins, including CD3ζ. The primary role of ZAP-70 is to phosphorylate LAT and SLP-76, which form a scaffold for signal propagation. We exploited the cooperative role of LAT and SLP-76 to engineer logic-gated intracellular network (LINK) CAR, a rapid and reversible Boolean-logic AND-gated CAR T cell platform that outperforms other systems in both efficacy and prevention of on-target, off-tumour toxicity. LINK CAR will expand the range of molecules that can be targeted with CAR T cells, and will enable these powerful therapeutic agents to be used for solid tumours and diverse diseases such as autoimmunity7 and fibrosis8. In addition, this work shows that the internal signalling machinery of cells can be repurposed into surface receptors, which could open new avenues for cellular engineering.

Mitochondrial single-cell ATAC-seq for high-throughput multi-omic detection of mitochondrial genotypes and chromatin accessibility.

Natural sequence variation within mitochondrial DNA (mtDNA) contributes to human phenotypes and may serve as natural genetic markers in human cells for clonal and lineage tracing. We recently developed a single-cell multi-omic approach, called 'mitochondrial single-cell assay for transposase-accessible chromatin with sequencing' (mtscATAC-seq), enabling concomitant high-throughput mtDNA genotyping and accessible chromatin profiling. Specifically, our technique allows the mitochondrial genome-wide inference of mtDNA variant heteroplasmy along with information on cell state and accessible chromatin variation in individual cells. Leveraging somatic mtDNA mutations, our method further enables inference of clonal relationships among native ex vivo-derived human cells not amenable to genetic engineering-based clonal tracing approaches. Here, we provide a step-by-step protocol for the use of mtscATAC-seq, including various cell-processing and flow cytometry workflows, by using primary hematopoietic cells, subsequent single-cell genomic library preparation and sequencing that collectively take ~3-4 days to complete. We discuss experimental and computational data quality control metrics and considerations for the extension to other mammalian tissues. Overall, mtscATAC-seq provides a broadly applicable platform to map clonal relationships between cells in human tissues, investigate fundamental aspects of mitochondrial genetics and enable additional modes of multi-omic discovery.

Massively Parallel Profiling of Accessible Chromatin and Proteins with ASAP-Seq.

While methods such as the Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) enable a comprehensive characterization of regulatory DNA, additional measurements are required to characterize the multifaceted nature of eukaryotic cells. Here, we delineate the ATAC with Select Antigen Profiling by sequencing (ASAP-seq) protocol, a scalable approach to quantifying proteins via oligo-tagged antibodies alongside accessible DNA in thousands of single cells. Critically, our method utilizes a custom bridge oligo that enables the utilization of a variety of oligo-conjugated antibodies, enabling the utilization and repurposing of other commercial products. The ASAP-seq method can be completed with straightforward experimental and computational modifications existing single-cell ATAC-seq workflows but yields distinct modalities underlying complex cellular states, including estimation of protein abundance on the cell surface as well as intracellular and intranuclear factors.

Concomitant Sequencing of Accessible Chromatin and Mitochondrial Genomes in Single Cells Using mtscATAC-Seq.

Mitochondria are unique organelles of eukaryotic cells that carry their own multicopy number and circular genome. In most mammals, including humans and mice, the size of the chromosome is ~16,000 base pairs and unlike nuclear DNA, mitochondrial DNA (mtDNA) is not densely compacted. This results in mtDNA to be highly accessible for enzymes such as the Tn5 transposase, commonly used for accessible chromatin profiling of nuclear chromatinized DNA. Here, we describe a method for the concomitant sequencing of mtDNA and accessible chromatin in thousands of individual cells via the mitochondrial single-cell assay for transposase accessible chromatin by sequencing (mtscATAC-seq). Our approach extends the utility of existing scATAC-seq products and protocols as we (Nam et al, Nat Rev Genet 22:3-18, 2021) fix cells using formaldehyde to retain mitochondria and its mtDNA within its originating cell, (Buenrostro et al, Nat Methods 10:1213-1218, 2013) modify lysis conditions to permeabilize cells and mitochondria, and (Corces et al, Nat Methods 14:959-962, 2017) optimize bioinformatic processing protocols to collectively increase mitochondrial genome coverage for downstream analysis. Here, we discuss the essentials for the experimental and computational methodologies to generate and analyze thousands of multiomic profiles of single cells over the course of a few days, enabling the profiling of accessible chromatin and mtDNA genotypes to reconstruct clonal relationships and studies of mitochondrial genetics and disease.