Search

This is the accepted manuscript version of the work published in its final form as Fonseca, R., Burn, T. N., Gandolfo, L. C., Devi, S., Park, S. L., Obers, A., Evrard, M., Christo, S. N., Buquicchio, F. A., Lareau, C. A., McDonald, K. M., Sandford, S. K., Zamudio, N. M., Zanluqui, N. G., Zaid, A., Speed, T. P., Satpathy, A. T., Mueller, S. N., Carbone, F. R., & Mackay, L. K. (2022). Runx3 drives a CD8+ T cell tissue residency program that is absent in CD4+ T cells.Nature Immunology,23(8), 1236-1245. https://doi.org/10.1038/s41590-022-01273-4 Deposited by shareyourpaper.org and openaccessbutton.org. We've taken reasonable steps to ensure this content doesn't violate copyright. However, if you think it does you can request a takedown by emailing help@openaccessbutton.org.

Supplementary Data from Mitochondrial DNA Mutations as Natural Barcodes for Lineage Tracing of Murine Tumor Models

Murine models are indispensable tools for functional genomic studies and preclinical testing of novel therapeutic approaches. Mitochondrial single-cell assay for transposase-accessible chromatin using sequencing (mtscATAC-seq) enables the dissection of cellular heterogeneity and clonal dynamics by capturing chromatin accessibility, copy-number variations (CNV), and mitochondrial DNA (mtDNA) mutations, yet its applicability to murine studies remains unexplored. By leveraging mtscATAC-seq in novel chronic lymphocytic leukemia and Richter syndrome mouse models, we report the detection of mtDNA mutations, particularly in highly proliferative murine cells, alongside CNV and chromatin state changes indicative of clonal evolution upon secondary transplant. This study thus demonstrates the feasibility and utility of multi-modal single-cell and natural barcoding approaches to characterize murine cancer models.Significance:mtDNA mutations can serve as natural barcodes to enable lineage tracing in murine cancer models, which can be used to provide new insights into disease biology and to identify therapeutic vulnerabilities.

Master regulators, such as the hematopoietic transcription factor (TF) GATA1, play an essential role in orchestrating lineage commitment and differentiation. However, the precise mechanisms by which such TFs regulate transcription through interactions with specific cis-regulatory elements remain incompletely understood. Here, we describe a form of congenital hemolytic anemia caused by missense mutations in an intrinsically disordered region of GATA1, with a poorly understood role in transcriptional regulation. Through integrative functional approaches, we demonstrate that these mutations perturb GATA1 transcriptional activity by partially impairing nuclear localization and selectively altering precise chromatin occupancy by GATA1. These alterations in chromatin occupancy and concordant chromatin accessibility changes alter faithful gene expression, with failure to both effectively silence and activate select genes necessary for effective terminal red cell production. We demonstrate how disease-causing mutations can reveal regulatory mechanisms that enable the faithful genomic targeting of master TFs during cellular differentiation.

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.

The state and behavior of a cell can be influenced by both genetic and environmental factors. In particular, tumor progression is determined by underlying genetic aberrations(1–4) as well as the makeup of the tumor microenvironment(5,6). Quantifying the contributions of these factors requires new technologies that can accurately measure the spatial location of genomic sequence together with phenotypic readouts. Here we developed slide-DNA-seq, a method for capturing spatially-resolved DNA sequences from intact tissue sections. We demonstrate that this method accurately preserves local tumor architecture and enables de novo discovery of distinct tumor clones and their copy number alterations. We then apply slide-DNA-seq to a mouse metastasis model and a primary human cancer, revealing that clonal populations are confined to distinct spatial regions. Moreover, through integration with spatial transcriptomics, we uncover distinct sets of genes that are associated with clone-specific genetic aberrations, the local tumor microenvironment, or both. Together, this multi-modal spatial genomics approach provides a versatile platform for quantifying how cell-intrinsic and extrinsic factors contribute to gene expression, protein abundance, and other cellular phenotypes.

The remarkable specificity of antibodies has enabled precision cancer immunotherapies, including chimeric antigen receptor T cells and antibody-drug conjugates. In parallel, single-cell genomics technologies present the possibility of a comprehensive annotation of antigen expression throughout tissues of the human body and on cancer cells. We reflect on the rationale for antigen targets currently used in immunotherapies, their adverse effects revealed in the clinic, and the opportunity to utilize large genomics datasets to de-risk potential targets and nominate optimal antigens for therapy.

While cancers evolve during disease progression and in response to therapy, temporal dynamics remain difficult to study in humans due to the lack of consistent barcodes marking individual clones in vivo. We employ mitochondrial single-cell assay for transposase-accessible chromatin with sequencing to profile 163,279 cells from 9 patients with chronic lymphocytic leukemia (CLL) collected across disease course and utilize mitochondrial DNA (mtDNA) mutations as natural genetic markers of cancer clones. We observe stable propagation of mtDNA mutations over years in the absence of strong selective pressure, indicating clonal persistence, but dramatic changes following tight bottlenecks, including disease transformation and relapse posttherapy, paralleled by acquisition of copy-number variants and changes in chromatin accessibility and gene expression. Furthermore, we link CLL subclones to distinct chromatin states, providing insight into nongenetic sources of relapse. mtDNA mutations thus mirror disease history and provide naturally occurring genetic barcodes to enable patient-specific study of cancer subclonal dynamics.Significance:Single-cell multi-omic profiling of CLL reveals the utility of somatic mtDNA mutations as in vivo barcodes, which mark subclones that can evolve over time along with changes in accessible chromatin and gene expression profiles to capture dynamics of disease evolution.See related commentary by Hilton and Scott, p. 2965.This article is highlighted in the In This Issue feature, p. 2945

The recent development of experimental methods for measuring chromatin state at single-cell resolution has created a need for computational tools capable of analyzing these datasets. Here we developed Signac, a comprehensive toolkit for the analysis of single-cell chromatin data. Signac enables an end-to-end analysis of single-cell chromatin data, including peak calling, quantification, quality control, dimension reduction, clustering, integration with single-cell gene expression datasets, DNA motif analysis and interactive visualization. Through its seamless compatibility with the Seurat package, Signac facilitates the analysis of diverse multimodal single-cell chromatin data, including datasets that co-assay DNA accessibility with gene expression, protein abundance and mitochondrial genotype. We demonstrate scaling of the Signac framework to analyze datasets containing over 700,000 cells.

Multimodal measurements of single-cell profiles are proving increasingly useful for characterizing cell states and regulatory mechanisms. In the present study, we developed PHAGE-ATAC (Assay for Transposase-Accessible Chromatin), a massively parallel droplet-based method that uses phage displaying, engineered, camelid single-domain antibodies ('nanobodies') for simultaneous single-cell measurements of protein levels and chromatin accessibility profiles, and mitochondrial DNA-based clonal tracing. We use PHAGE-ATAC for multimodal analysis in primary human immune cells, sample multiplexing, intracellular protein analysis and the detection of SARS-CoV-2 spike protein in human cell populations. Finally, we construct a synthetic high-complexity phage library for selection of antigen-specific nanobodies that bind cells of particular molecular profiles, opening an avenue for protein detection, cell characterization and screening with single-cell genomics.

Cells experience intrinsic and extrinsic pressures that affect their proclivity to expand and persistin vivo. In congenital disorders caused by loss-of-function mutations in mitochondrial DNA (mtDNA), metabolic vulnerabilities may result in cell-type specific phenotypes and depletion of pathogenic alleles, contributing to purifying selection. However, the impact of pathogenic mtDNA mutations on the cellular hematopoietic landscape is not well understood. Here, we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from Pearson syndrome patients. We resolve the interdependence between pathogenic mtDNA and lineage, including purifying selection against deletions in effector/memory CD8 T-cell populations and recent thymic emigrants and dynamics in other hematopoietic populations. Our mapping of lineage-specific purifying selection dynamics in primary cells from patients carrying pathogenic heteroplasmy provides a new perspective on recurrent clinical phenotypes in mitochondrial disorders, including cancer and infection, with potential broader relevance to age-related immune dysfunction.

The mutualistic relationship of gut-resident microbiota and the host immune system promotes homeostasis that ensures maintenance of the microbial community and of a largely non-aggressive immune cell compartment1,2. The consequences of disturbing this balance include proximal inflammatory conditions, such as Crohn's disease, and systemic illnesses. This equilibrium is achieved in part through the induction of both effector and suppressor arms of the adaptive immune system. Helicobacter species induce T regulatory (Treg) and T follicular helper (TFH) cells under homeostatic conditions, but induce inflammatory T helper 17 (TH17) cells when induced Treg (iTreg) cells are compromised3,4. How Helicobacter and other gut bacteria direct T cells to adopt distinct functions remains poorly understood. Here we investigated the cells and molecular components required for iTreg cell differentiation. We found that antigen presentation by cells expressing RORγt, rather than by classical dendritic cells, was required and sufficient for induction of Treg cells. These RORγt+ cells-probably type 3 innate lymphoid cells and/or Janus cells5-require the antigen-presentation machinery, the chemokine receptor CCR7 and the TGFβ activator αv integrin. In the absence of any of these factors, there was expansion of pathogenic TH17 cells instead of iTreg cells, induced by CCR7-independent antigen-presenting cells. Thus, intestinal commensal microbes and their products target multiple antigen-presenting cells with pre-determined features suited to directing appropriate T cell differentiation programmes, rather than a common antigen-presenting cell that they endow with appropriate functions.

Recent technological advances have enabled massively parallel chromatin profiling with scATAC-seq (single-cell assay for transposase accessible chromatin by sequencing). Here we present ATAC with select antigen profiling by sequencing (ASAP-seq), a tool to simultaneously profile accessible chromatin and protein levels. Our approach pairs sparse scATAC-seq data with robust detection of hundreds of cell surface and intracellular protein markers and optional capture of mitochondrial DNA for clonal tracking, capturing three distinct modalities in single cells. ASAP-seq uses a bridging approach that repurposes antibody:oligonucleotide conjugates designed for existing technologies that pair protein measurements with single-cell RNA sequencing. Together with DOGMA-seq, an adaptation of CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) for measuring gene activity across the central dogma of gene regulation, we demonstrate the utility of systematic multi-omic profiling by revealing coordinated and distinct changes in chromatin, RNA and surface proteins during native hematopoietic differentiation and peripheral blood mononuclear cell stimulation and as a combinatorial decoder and reporter of multiplexed perturbations in primary 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 still being understood. For example, we currently lack a comprehensive understanding of the mechanisms of neurotoxicity observed in patients receiving T cell therapies, including recent reports of encephalitis caused by human herpesvirus 6 (HHV-6) reactivation1. Here, via petabase-scale viral RNA data mining, we examine the landscape of human latent viral reactivation and demonstrate that HHV-6B can become reactivated in human CD4+ T cells in standard in vitro cultures. Using single-cell sequencing, we identify a rare population of HHV-6 ‘super-expressors’ (~1 in 300-10,000 cells) that possess high viral transcriptional activity in chimeric antigen receptor (CAR) T cell culture before spreading to infect other cells in vitro. Through the analysis of single-cell sequencing data from patients receiving cell therapy products that are FDA-approved2 or used in clinical studies3,4, we identify the presence of CAR+, HHV-6 super-expressor T cells in vivo. Together, our study implicates cell therapy products as a source of lytic HHV-6 reported in clinical trials1,5–7 and has broad implications for the design, production, and monitoring of cell therapies.

Murine models are indispensable tools for functional genomic studies and preclinical testing of novel therapeutic approaches. Mitochondrial single-cell assay for transposase-accessible chromatin using sequencing (mtscATAC-seq) enables the dissection of cellular heterogeneity and clonal dynamics by capturing chromatin accessibility, copy-number variations (CNV), and mitochondrial DNA (mtDNA) mutations, yet its applicability to murine studies remains unexplored. By leveraging mtscATAC-seq in novel chronic lymphocytic leukemia and Richter syndrome mouse models, we report the detection of mtDNA mutations, particularly in highly proliferative murine cells, alongside CNV and chromatin state changes indicative of clonal evolution upon secondary transplant. This study thus demonstrates the feasibility and utility of multi-modal single-cell and natural barcoding approaches to characterize murine cancer models.Significance:mtDNA mutations can serve as natural barcodes to enable lineage tracing in murine cancer models, which can be used to provide new insights into disease biology and to identify therapeutic vulnerabilities.

Characterizing the interactions that SARS-CoV-2 viral RNAs make with host cell proteins during infection can improve our understanding of viral RNA functions and the host innate immune response. Using RNA antisense purification and mass spectrometry, we identified up to 104 human proteins that directly and specifically bind to SARS-CoV-2 RNAs in infected human cells. We integrated the SARS-CoV-2 RNA interactome with changes in proteome abundance induced by viral infection and linked interactome proteins to cellular pathways relevant to SARS-CoV-2 infections. We demonstrated by genetic perturbation that cellular nucleic acid-binding protein (CNBP) and La-related protein 1 (LARP1), two of the most strongly enriched viral RNA binders, restrict SARS-CoV-2 replication in infected cells and provide a global map of their direct RNA contact sites. Pharmacological inhibition of three other RNA interactome members, PPIA, ATP1A1, and the ARP2/3 complex, reduced viral replication in two human cell lines. The identification of host dependency factors and defence strategies as presented in this work will improve the design of targeted therapeutics against SARS-CoV-2., Interactions between SARS-CoV-2 viral RNAs and host cell proteins during infection are evaluated to improve our understanding of viral RNA functions and the host innate immune response.

Many mitochondrial diseases are caused by mutations in mitochondrial DNA (mtDNA). Patients' cells contain a mixture of mutant and nonmutant mtDNA (a phenomenon called heteroplasmy). The proportion of mutant mtDNA varies across patients and among tissues within a patient. We simultaneously assayed single-cell heteroplasmy and cell state in thousands of blood cells obtained from three unrelated patients who had A3243G-associated mitochondrial encephalomyopathy, lactic acidosis, and strokelike episodes. We observed a broad range of heteroplasmy across all cell types but also found markedly reduced heteroplasmy in T cells, a finding consistent with purifying selection within this lineage. We observed this pattern in six additional patients who had heteroplasmic A3243G without strokelike episodes. (Funded by the Marriott Foundation and others.).

Existing adenine and cytosine base editors induce only a single type of modification, limiting the range of DNA alterations that can be created. Here we describe a CRISPR-Cas9-based synchronous programmable adenine and cytosine editor (SPACE) that can concurrently introduce A-to-G and C-to-T substitutions with minimal RNA off-target edits. SPACE expands the range of possible DNA sequence alterations, broadening the research applications of CRISPR base editors., Editorial summary Base editor that concurrently modifies both adenine and cytosine broadens the potential applications of base editing.

Genome-wide association studies have associated thousands of genetic variants with complex traits and diseases, but pinpointing the causal variant(s) among those in tight linkage disequilibrium with each associated variant remains a major challenge. Here, we use seven experimental assays to characterize all common variants at the multiple disease-associated TNFAIP3 locus in five disease-relevant immune cell lines, based on a set of features related to regulatory potential. Trait/disease-associated variants are enriched among SNPs prioritized based on either: (1) residing within CRISPRi-sensitive regulatory regions, or (2) localizing in a chromatin accessible region while displaying allele-specific reporter activity. Of the 15 trait/disease-associated haplotypes at TNFAIP3, 9 have at least one variant meeting one or both of these criteria, 5 of which are further supported by genetic fine-mapping. Our work provides a comprehensive strategy to characterize genetic variation at important disease-associated loci, and aids in the effort to identify trait causal genetic variants., While genome-wide association studies have yielded thousands of trait-associated loci, identifying causal variants remains challenging. Here, the authors perform seven genomics assays in various cell types to prioritize genetic variants in the TNFAIP3 locus, and report high-priority variants within disease-associated haplotypes.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

The mutualistic relationship of gut-resident microbiota and cells of the host immune system promotes homeostasis that ensures maintenance of the microbial community and of a poised, but largely non-aggressive, immune cell compartment1, 2. Consequences of disturbing this balance, by environmental or genetic factors, include proximal inflammatory conditions, like Crohn’s disease, and systemic illnesses, both metabolic and autoimmune. One of the means by which this equilibrium is achieved is through induction of both effector and suppressor or regulatory arms of the adaptive immune system. In mice, Helicobacter species induce regulatory (iTreg) and follicular helper (Tfh) T cells in the colon-draining mesenteric lymph nodes under homeostatic conditions, but can instead induce inflammatory Th17 cells when iTreg cells are compromised3, 4. How Helicobacter hepaticus and other gut bacteria direct T cells to adopt distinct functions remains poorly understood. Here, we investigated which cells and molecular components are required to convey the microbial instruction for the iTreg differentiation program. We found that antigen presentation by cells expressing RORγt, rather than by classical dendritic cells, was both required and sufficient for iTreg induction. These RORγt+ cells, likely to be type 3 innate lymphoid cells (ILC3) and/or a recently-described population of Aire+ cells termed Janus cells5, require the MHC class II antigen presentation machinery, the chemokine receptor CCR7, and αv integrin, which activates TGF-β, for iTreg cell differentiation. In the absence of any of these, instead of iTreg cells there was expansion of microbiota-specific pathogenic Th17 cells, which were induced by other antigen presenting cells (APCs) that did not require CCR7. Thus, intestinal commensal microbes and their products target multiple APCs with pre-determined features suited to directing appropriate T cell differentiation programs, rather than a common APC that they endow with appropriate functions. Our results illustrate the ability of microbiota to exploit specialized functions of distinct innate immune system cells, targeting them to achieve the desired composition of equipoised T cells, thus maintaining tolerance.

The precise timing of T cell priming during infection remains unclear. Here, we mapped the cellular dynamics of all immune lineages during acute infection with Listeria monocytogenes (Lm). We identified highly transient DC activation 2 days post-infection that functions as a critical time window for priming effector T cells. Regulation of this transient state was mediated by DC extrinsic IFN{gamma} provided by lymphocytes. Furthermore, antigen-specific T cells that arrive late to the site of priming and miss peak DC activation acquire only memory T cell fates. This temporal regulation of fate is recapitulated by CD8+ DCs ex vivo, suggesting that shifts in activation state of a single antigen presenting cell population alter T cell fates. These results uncover a novel mechanism for temporal regulation of T cell differentiation during a dynamic immune response to acute infection. One-Sentence SummaryThe immune system generates a balanced protective response through rapidly shifting activation states and recruitment of new cells into an ongoing inflammatory landscape.

The germinal center (GC) response is critical for both effective adaptive immunity and establishing peripheral tolerance by limiting autoreactive B cells. Dysfunction in these processes can lead to defective immune responses to infection or contribute to autoimmune disease. To understand the gene regulatory principles underlying the GC response, we generated a single-cell transcriptomic and epigenomic atlas of the human tonsil, a widely studied and representative lymphoid tissue. We characterize diverse immune cell subsets and build a trajectory of dynamic gene expression and transcription factor activity during B cell activation, GC formation, and plasma cell differentiation. We subsequently leverage cell type-specific transcriptomic and epigenomic maps to interpret potential regulatory impact of genetic variants implicated in autoimmunity, revealing that many exhibit their greatest regulatory potential in GC-associated cellular populations. These included gene loci linked with known roles in GC biology (IL21, IL21R, IL4R, BCL6) and transcription factors regulating B cell differentiation (POU2AF1, HHEX). Together, these analyses provide a powerful new cell type-resolved resource for the interpretation of cellular and genetic causes underpinning autoimmune disease.

BACKGROUND: Idiopathic aplastic anemia is a potentially lethal disease, characterized by T-cell mediated autoimmune attack of bone marrow hematopoietic stem cells. Standard of care therapies (stem cell transplantation or immunosuppression) are effective but associated with a risk of serious toxicities. METHODS: An 18-year-old man presented with aplastic anemia in the context of a germline gain-of-function variant in STAT1. Treatment with the JAK1 inhibitor itacitinib resulted in a rapid resolution of aplastic anemia and a sustained recovery of hematopoiesis. Peripheral blood and bone marrow samples were compared before and after JAK1 inhibitor therapy. FINDINGS: Following therapy, samples showed a decrease in the plasma concentration of interferon-γ, a decrease in PD1-positive exhausted CD8+ T cells population, and a decrease in an interferon responsive myeloid population. Single cell analysis of chromatin accessibility showed decreased accessibility of STAT1 across CD4+ and CD8+ T cells, as well as CD14+ monocytes. To query whether other cases of aplastic anemia share a similar STAT1-mediated pathophysiology, we examined a cohort of patients with idiopathic aplastic anemia. Bone marrow from six of nine patients also displayed abnormal STAT1 hyper-activation. CONCLUSIONS: These findings raise the possibility that STAT1-hyperactivition defines a subset of idiopathic aplastic anemia patients for whom JAK inhibition may be an efficacious therapy.

The precise timing of T cell priming during infection remains unclear. Here, we mapped the cellular dynamics of all immune lineages during acute infection with Listeria monocytogenes (Lm). We identified highly transient DC activation 2 days post-infection that functions as a critical time window for priming effector T cells. Regulation of this transient state was mediated by DC extrinsic IFNγ provided by lymphocytes. Furthermore, antigen-specific T cells that arrive late to the site of priming and miss peak DC activation acquire only memory T cell fates. This temporal regulation of fate is recapitulated by CD8+ DCs ex vivo, suggesting that shifts in activation state of a single antigen presenting cell population alter T cell fates. These results uncover a novel mechanism for temporal regulation of T cell differentiation during a dynamic immune response to acute infection.One-Sentence SummaryThe immune system generates a balanced protective response through rapidly shifting activation states and recruitment of new cells into an ongoing inflammatory landscape.

Cells require coordinated control over gene expression when responding to environmental stimuli. Here we apply scATAC-seq and single-cell RNA sequencing (scRNA-seq) in resting and stimulated human blood cells. Collectively, we generate ~91,000 single-cell profiles, allowing us to probe the cis-regulatory landscape of the immunological response across cell types, stimuli, and time. Advancing tools to integrate multi-omics data, we develop functional inference of gene regulation (FigR), a framework to computationally pair scA-TAC-seq with scRNA-seq cells, connect distal cis-regulatory elements to genes, and infer gene-regulatory networks (GRNs) to identify candidate transcription factor (TF) regulators. Utilizing these paired multi-omics data, we define domains of regulatory chromatin (DORCs) of immune stimulation and find that cells alter chromatin accessibility and gene expression at timescales of minutes. Construction of the stimulation GRN elucidates TF activity at disease-associated DORCs. Overall, FigR enables elucidation of regulatory interactions across single-cell data, providing new opportunities to understand the function of cells within tissues.

Cells require coordinated control over gene expression when responding to environmental stimuli. Here, we apply scATAC-seq and scRNA-seq in resting and stimulated human blood cells. Collectively, we generate ∼91,000 single-cell profiles, allowing us to probe the cis-regulatory landscape of immunological response across cell types, stimuli and time. Advancing tools to integrate multi-omic data, we develop FigR - a framework to computationally pair scATAC-seq with scRNA-seq cells, connect distal cis-regulatory elements to genes, and infer gene regulatory networks (GRNs) to identify candidate TF regulators. Utilizing these paired multi-omic data, we define Domains of Regulatory Chromatin (DORCs) of immune stimulation and find that cells alter chromatin accessibility prior to production of gene expression at time scales of minutes. Further, the construction of the stimulation GRN elucidates TF activity at disease-associated DORCs. Overall, FigR enables the elucidation of regulatory interactions across single-cell data, providing new opportunities to understand the function of cells within tissues.

Cytosine or adenine base editors (CBEs or ABEs) can introduce specific DNA C-to-T or A-to-G alterations1-4. However, we recently demonstrated that they can also induce transcriptome-wide guide-RNA-independent editing of RNA bases5, and created selective curbing of unwanted RNA editing (SECURE)-BE3 variants that have reduced unwanted RNA-editing activity5. Here we describe structure-guided engineering of SECURE-ABE variants with reduced off-target RNA-editing activity and comparable on-target DNA-editing activity that are also among the smallest Streptococcus pyogenes Cas9 base editors described to date. We also tested CBEs with cytidine deaminases other than APOBEC1 and found that the human APOBEC3A-based CBE induces substantial editing of RNA bases, whereas an enhanced APOBEC3A-based CBE6, human activation-induced cytidine deaminase-based CBE7, and the Petromyzon marinus cytidine deaminase-based CBE Target-AID4 induce less editing of RNA. Finally, we found that CBEs and ABEs that exhibit RNA off-target editing activity can also self-edit their own transcripts, thereby leading to heterogeneity in base-editor coding sequences.

Recent technical advancements have facilitated the mapping of epigenomes at single-cell resolution; however, the throughput and quality of these methods have limited their widespread adoption. Here we describe a high-quality (10(5) nuclear fragments per cell) droplet-microfluidics-based method for single-cell profiling of chromatin accessibility. We use this approach, named ‘droplet single-cell assay for transposase-accessible chromatin using sequencing’ (dscATAC-seq), to assay 46,653 cells for the unbiased discovery of cell types and regulatory elements in adult mouse brain. We further increase the throughput of this platform by combining it with combinatorial indexing (dsciATAC-seq), enabling single-cell studies at a massive scale. We demonstrate the utility of this approach by measuring chromatin accessibility across 136,463 resting and stimulated human bone marrow-derived cells to reveal changes in the cis- and trans-regulatory landscape across cell types and under stimulatory conditions at single-cell resolution. Altogether, we describe a total of 510,123 single-cell profiles, demonstrating the scalability and flexibility of this droplet-based platform.

Single-cell transcriptomic assays have enabled the de novo reconstruction of lineage differentiation trajectories, along with the characterization of cellular heterogeneity and state transitions. Several methods have been developed for reconstructing developmental trajectories from single-cell transcriptomic data, but efforts on analyzing single-cell epigenomic data and on trajectory visualization remain limited. Here we present STREAM, an interactive pipeline capable of disentangling and visualizing complex branching trajectories from both single-cell transcriptomic and epigenomic data. We have tested STREAM on several synthetic and real datasets generated with different single-cell technologies. We further demonstrate its utility for understanding myoblast differentiation and disentangling known heterogeneity in hematopoiesis for different organisms. STREAM is an open-source software package., The increasing accessibility of single cell omics technologies beyond transcriptomics demands parallel advances in analysis. Here, the authors introduce STREAM, a pipeline for reconstruction and visualization of differentiation trajectories from both single-cell RNA-seq and ATAC-seq data.

Perturbations to mammalian switch/sucrose non-fermentable (mSWI/SNF) chromatin remodeling complexes have been widely implicated as driving events in cancer1. One such perturbation is the dual loss of the SMARCA4 and SMARCA2 ATPase subunits in small cell carcinoma of the ovary, hypercalcemic type (SCCOHT)2-5, SMARCA4-deficient thoracic sarcomas6 and dedifferentiated endometrial carcinomas7. However, the consequences of dual ATPase subunit loss on mSWI/SNF complex subunit composition, chromatin targeting, DNA accessibility and gene expression remain unknown. Here we identify an ATPase module of subunits that is required for functional specification of the Brahma-related gene-associated factor (BAF) and polybromo-associated BAF (PBAF) mSWI/SNF family subcomplexes. Using SMARCA4/2 ATPase mutant variants, we define the catalytic activity-dependent and catalytic activity-independent contributions of the ATPase module to the targeting of BAF and PBAF complexes on chromatin genome-wide. Finally, by linking distinct mSWI/SNF complex target sites to tumor-suppressive gene expression programs, we clarify the transcriptional consequences of SMARCA4/2 dual loss in SCCOHT.

Widespread linkage disequilibrium and incomplete annotation of cell-to-cell state variation represent substantial challenges to elucidating mechanisms of trait-associated genetic variation. Here we perform genetic fine-mapping for blood cell traits in the UK Biobank to identify putative causal variants. These variants are enriched in genes encoding proteins in trait-relevant biological pathways and in accessible chromatin of hematopoietic progenitors. For regulatory variants, we explore patterns of developmental enhancer activity, predict molecular mechanisms, and identify likely target genes. In several instances, we localize multiple independent variants to the same regulatory element or gene. We further observe that variants with pleiotropic effects preferentially act in common progenitor populations to direct the production of distinct lineages. Finally, we leverage fine-mapped variants in conjunction with continuous epigenomic annotations to identify trait-cell type enrichments within closely related populations and in single cells. Our study provides a comprehensive framework for single-variant and single-cell analyses of genetic associations.