Search

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

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.

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.

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.

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.

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.

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.

The germinal center (GC) response is critical for both effective adaptive immunity and establishing peripheral tolerance by limiting auto-reactive B cells. Dysfunction in these processes can lead to defects in immune response to pathogens 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 cell populations. Together, these analyses provide a powerful new cell type-resolved resource for the interpretation of cellular and genetic causes underpinning autoimmune disease.One sentence summarySingle-cell chromatin accessibility landscapes of immune cell subsets reveal regulatory potential of autoimmune-associated genetic variants during the germinal center response.

Reconstructing lineage relationships in complex tissues can reveal mechanisms underlying development and disease. Recent methods combine single-cell transcriptomics with mitochondrial DNA variant detection to establish lineage relationships in primary human cells, but are not scalable to interrogate complex tissues. To overcome this limitation, here we develop a technology for high-confidence detection of mitochondrial mutations from high-throughput single-cell RNA-sequencing. We use the new method to identify skewed immune cell expansions in primary human clonal hematopoiesis.

Sarcoidosis is a systemic inflammatory disease characterized by infiltration of immune cells into granulomas. Previous gene expression studies using heterogeneous cell mixtures lack insight into cell-type-specific immune dysregulation. We performed the first single-cell RNA-sequencing study of sarcoidosis in peripheral immune cells in 48 patients and controls. Following unbiased clustering, differentially expressed genes were identified for 18 cell types and bioinformatically assessed for function and pathway enrichment. Our results reveal persistent activation of circulating classical monocytes with subsequent upregulation of trafficking molecules. Specifically, classical monocytes upregulated distinct markers of activation including adhesion molecules, pattern recognition receptors, and chemokine receptors, as well as enrichment of immunoregulatory pathways HMGB1, mTOR, and ephrin receptor signaling. Predictive modeling implicated TGFβ and mTOR signaling as drivers of persistent monocyte activation. Additionally, sarcoidosis T cell subsets displayed patterns of dysregulation. CD4 naïve T cells were enriched for markers of apoptosis and Th17/Treg differentiation, while effector T cells showed enrichment of anergy-related pathways. Differentially expressed genes in regulatory T cells suggested dysfunctional p53, cell death, and TNFR2 signaling. Using more sensitive technology and more precise units of measure, we identify cell-type specific, novel inflammatory and regulatory pathways. Based on our findings, we suggest a novel model involving four convergent arms of dysregulation: persistent hyperactivation of innate and adaptive immunity via classical monocytes and CD4 naïve T cells, regulatory T cell dysfunction, and effector T cell anergy. We further our understanding of the immunopathology of sarcoidosis and point to novel therapeutic targets.

© 2020 Elsevier Inc. Cell differentiation and function are regulated across multiple layers of gene regulation, including modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of regulatory events leading to cell fate commitment. Here we developed simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq), a highly scalable approach for measurement of chromatin accessibility and gene expression in the same single cell, applicable to different tissues. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define domains of regulatory chromatin (DORCs) that significantly overlap with super-enhancers. During lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for lineage commitment. We computationally infer chromatin potential as a quantitative measure of chromatin lineage-priming and use it to predict cell fate outcomes. SHARE-seq is an extensible platform to study regulatory circuitry across diverse cells in tissues.

Multi-modal measurements of single cell profiles are a powerful tool for characterizing cell states and regulatory mechanisms. While current methods allow profiling of RNA along with either chromatin or protein levels, connecting chromatin state to protein levels remains a barrier. Here, we developed PHAGE-ATAC, a method that uses engineered camelid single-domain antibody (‘nanobody’)-displaying phages for simultaneous single-cell measurement of surface proteins, chromatin accessibility profiles, and mtDNA-based clonal tracing through a massively parallel droplet-based assay of single-cell transposase-accessible chromatin with sequencing (ATAC-seq). We demonstrate PHAGE-ATAC for multimodal analysis in primary human immune cells and for sample multiplexing. Finally, we construct a synthetic high-complexity phage library for selection of novel antigen-specific nanobodies that bind cells of particular molecular profiles, opening a new avenue for protein detection, cell characterization and screening with single-cell genomics.

The transcription factor Foxp3 plays crucial roles for Treg cell development and function. Conserved non-coding sequences (CNSs) at the Foxp3 locus control Foxp3 transcription, but how they developmentally contribute to Treg cell lineage specification remains obscure. Here, we show that among Foxp3 CNSs, the promoter-upstream CNS0 and the intergenic CNS3, which bind distinct transcription factors, were activated at early stages of thymocyte differentiation prior to Foxp3 promoter activation, with sequential genomic looping bridging these regions and the promoter. While deletion of either CNS0 or CNS3 partially compromised thymic Treg cell generation, deletion of both completely abrogated the generation and impaired the stability of Foxp3 expression in residual Treg cells. As a result, CNS0 and CNS3 double-deleted mice succumbed to lethal systemic autoimmunity and inflammation. Thus, hierarchical and coordinated activation of Foxp3 CNS0 and CNS3 initiates and stabilizes Foxp3 gene expression, thereby crucially controlling Treg cell development, maintenance, and consequently immunological self-tolerance.

SARS-CoV-2 infections pose a global threat to human health and an unprecedented research challenge. Among the most urgent tasks is obtaining a detailed understanding of the molecular interactions that facilitate viral replication or contribute to host defense mechanisms in infected cells. While SARS-CoV-2 co-opts cellular factors for viral translation and genome replication, a comprehensive map of the host cell proteome in direct contact with viral RNA has not been elucidated. Here, we use RNA antisense purification and mass spectrometry (RAP-MS) to obtain an unbiased and quantitative picture of the human proteome that directly binds the SARS-CoV-2 RNA in infected human cells. We discover known host factors required for coronavirus replication, regulators of RNA metabolism and host defense pathways, along with dozens of potential drug targets among direct SARS-CoV-2 binders. We further integrate the SARS-CoV-2 RNA interactome with proteome dynamics induced by viral infection, linking interactome proteins to the emerging biology of SARS-CoV-2 infections. Validating RAP-MS, we show that CNBP, a regulator of proinflammatory cytokines, directly engages the SARS-CoV-2 RNA. Supporting the functional relevance of identified interactors, we show that the interferon-induced protein RYDEN suppresses SARS-CoV-2 ribosomal frameshifting and demonstrate that inhibition of SARS-CoV-2-bound proteins is sufficient to manipulate viral replication. The SARS-CoV-2 RNA interactome provides an unprecedented molecular perspective on SARS-CoV-2 infections and enables the systematic dissection of host dependency factors and host defense strategies, a crucial prerequisite for designing novel therapeutic strategies.

SummaryCell differentiation and function are regulated across multiple layers of gene regulation, including the modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of the regulatory events leading to cell fate commitment. Here, we developed SHARE-seq, a highly scalable approach for measurement of chromatin accessibility and gene expression within the same single cell. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define Domains of Regulatory Chromatin (DORCs), which significantly overlap with super-enhancers. We show that during lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting changes in chromatin accessibility may prime cells for lineage commitment. We therefore develop a computational strategy (chromatin potential) to quantify chromatin lineage-priming and predict cell fate outcomes. Together, SHARE-seq provides an extensible platform to study regulatory circuitry across diverse cells within tissues.

SummaryBlood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including 563,946 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering the full allele frequency spectrum of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood cell traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell GWAS to interrogate clinically meaningful variants across the full allelic spectrum of human variation.

SUMMARYMost loci identified by GWAS have been found in populations of European ancestry (EA). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EA individuals, we identified 5,552 trait-variant associations at P−9, including 71 novel loci not found in EA populations. We also identified novel ancestry-specific variants not found in EA, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL7 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 EA-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations, and compared genetic architecture and the impact 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.

Regulatory networks that maintain functional, differentiated cell states are often dysregulated in tumor development. Here, we use single-cell epigenomics to profile chromatin state transitions in a mouse model of lung adenocarcinoma (LUAD). We identify an epigenomic continuum representing loss of cellular identity and progression toward a metastatic state. We define co-accessible regulatory programs and infer key activating and repressive chromatin regulators of these cell states. Among these co-accessibility programs, we identify a pre-metastatic transition, characterized by activation of RUNX transcription factors, which mediates extracellular matrix remodeling to promote metastasis and is predictive of survival across human LUAD patients. Together, these results demonstrate the power of single-cell epigenomics to identify regulatory programs to uncover mechanisms and key biomarkers of tumor progression.

Summary Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation., Graphical Abstract, Highlights • Largest genome-wide association study of blood cell traits to date • Empiric assessments of omnigenic and infinitesimal models of polygenic variation • Functional insights into how genetic variants impact human hematopoiesis • Assessment of the effect of polygenic trait scores upon blood diseases, Analysis of blood cell traits in the UK Biobank and other cohorts illuminates the full genetic architecture of hematopoietic phenotypes, with evidence supporting the omnigenic model for complex traits and linking polygenic burden with monogenic blood diseases.

Reconstitution of donor hematopoiesis after allogeneic hematopoietic stem cell transplantation (HSCT) forms the basis for effective graft-versus-leukemia responses, but mixed chimerism is not an infrequent outcome. How the donor and host hematopoietic system interact under conditions of mixed chimerism remains incompletely understood. Multi-modal single cell sequencing platforms are increasingly available and provide information regarding cell identities and interactions at high resolution. However, the analysis of post-transplant immune reconstitution requires consistent distinguishing of donor- and recipient-derived cells, which for sparse single cell sequencing data until now has remained a challenge. Recently, mitochondrial DNA (mtDNA) mutations have been recognized for their potential as personal genetic barcodes that can be detected with mitochondrial single cell assay for transposase-accessible chromatin with sequencing (mtscATAC-seq). We hypothesized that individual-specific mtDNA mutations could provide a sensitive and robust approach for distinguishing donor- from recipient-derived cells, and therefore tested this approach on bone marrow (BM) samples from patients with relapsed acute myeloid leukemia (AML) post-HSCT. We employed ATAC with select cell surface antigen profiling by sequencing (ASAP-seq), which enables the detection of mtDNA mutations within distinct surface marker-defined cell populations alongside chromatin accessibility. We selected serial samples collected from the ETCTN 10026 study, which tested combined decitabine (days 1-5, every 4 weeks, start cycle 0) and ipilimumab (day 1, every 4 weeks, start cycle 1) in relapsed AML post-HSCT. We focused on 13 samples (study entry, on treatment and disease progression) from 3 patients: AML1012 (HSCT from a matched related donor [MRD]), and AML1010 and AML1026 (matched unrelated donor [MUD]-HSCT). In total, we obtained 33,943 ASAP-seq profiles, including 3,283 single T cells. While clustering using single cell chromatin profiles alone only allowed identification of either CD4 + or CD8 + T cells, integration with surface marker expression enabled more detailed annotations of 8 T cell subpopulations and NK cells. Further, phenotypically distinct subpopulations such as CD57 + CD4 + and CD8 + T cells shared highly similar chromatin profiles, and 11.1% of CD4 + and 33.7% of CD8 + T cells would have been mislabeled based on clustering of chromatin profiles alone. Thus, ASAP-seq identified T cell subsets with markedly improved accuracy and resolution than scATAC-seq alone. Upon evaluation of mtDNA mutations to discriminate donor- and recipient-derived single T cells, we found that this was unreliable for MRD-HSCT (AML1012), but highly robust in the setting of MUD-HSCT (AML1010, AML1026), consistent with maternal inheritance of mitochondrial genomes. For the latter two patients, we identified 48 donor- and 26 recipient-specific mtDNA mutations, all with high heteroplasmy (range 82 - 99%). Presence of donor- and recipient-derived mtDNA mutations was mutually exclusive, and recipient-specific mtDNA mutations were also detectable in AML cells. Clinical bulk and mtDNA mutation-based single T cell chimerisms were highly correlated (r = 0.97). AML1010 had sustained complete T cell chimerism (>97%) during study treatment. In AML1026, the mtDNA mutation-based T cell chimerism rose from 55% to 71% after 1 cycle of decitabine and then remained stable until disease progression 3 months later. This was associated with increased percentage of donor-derived CD4 + T cells (45% [study entry] vs. 71% [after 1 cycle of decitabine], p < 0.01), while donor-derived CD8 + T cells remained unchanged at 76%. Across all studied timepoints in AML1026, donor versus recipient skewing was also highest in CD4 + T cell subsets, with fewer naïve (20% vs. 31%, p < 0.01) but more donor-derived CD57 + CD4 + T cells (13% vs. 3%, p < 0.01). We demonstrate that mtDNA mutations can discriminate between donor- and recipient-derived single cells, enabling detection and in-depth characterization of chimeric immune cell dynamics after MUD HSCT. This approach will allow to systematically dissect conditions of mixed chimerism in the post-transplant setting with larger studies. Disclosures DeAngelo: Abbvie: Research Funding; Takeda: Consultancy; Servier: Consultancy; Pfizer: Consultancy; Novartis: Consultancy, Research Funding; Jazz: Consultancy; Incyte: Consultancy; Forty-Seven: Consultancy; Autolus: Consultancy; Amgen: Consultancy; Agios: Consultancy; Blueprint: Research Funding; Glycomimetrics: Research Funding. Neuberg: Madrigal Pharmaceuticals: Other: Stock ownership; Pharmacyclics: Research Funding. Sankaran: Cellarity: Consultancy; Forma: Consultancy; Novartis: Consultancy; Branch Biosciences: Consultancy; Ensoma: Consultancy. Soiffer: Jasper: Consultancy; Jazz Pharmaceuticals, USA: Consultancy; Precision Biosciences, USA: Consultancy; Juno Therapeutics, USA: Other: Data Safety Monitoring Board; Kiadis, Netherlands: Membership on an entity's Board of Directors or advisory committees; Rheos Therapeutics, USA: Consultancy; Gilead, USA: Other: Career Development Award Committee; NMPD - Be the Match, USA: Membership on an entity's Board of Directors or advisory committees; Takeda: Consultancy. Garcia: Genentech: Research Funding; Prelude: Research Funding; Pfizer: Research Funding; AstraZeneca: Research Funding; Takeda: Consultancy, Membership on an entity's Board of Directors or advisory committees; Astellas: Consultancy, Membership on an entity's Board of Directors or advisory committees; AbbVie: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding. Wu: Pharmacyclics: Research Funding; BioNTech: Current equity holder in publicly-traded company. OffLabel Disclosure: ipilimumab to modulate anti-leukemia immunity in the post-transplant and transplant-naive context

Summary Gain-of-function Notch mutations are recurrent in mature small B cell lymphomas such as mantle cell lymphoma (MCL) and chronic lymphocytic leukemia (CLL), but the Notch target genes that contribute to B cell oncogenesis are largely unknown. We performed integrative analysis of Notch-regulated transcripts, genomic binding of Notch transcription complexes, and genome conformation data to identify direct Notch target genes in MCL cell lines. This B cell Notch regulome is largely controlled through Notch-bound distal enhancers and includes genes involved in B cell receptor and cytokine signaling and the oncogene MYC , which sustains proliferation of Notch-dependent MCL cell lines via a Notch-regulated lineage-restricted enhancer complex. Expression of direct Notch target genes is associated with Notch activity in an MCL xenograft model and in CLL lymph node biopsies. Our findings provide key insights into the role of Notch in MCL and other B cell malignancies and have important implications for therapeutic targeting of Notch-dependent oncogenic pathways.

Sarcoidosis is a complex disease of unknown etiology characterized by the presence of granulomatous inflammation. Though various immune system pathways have been implicated in disease, the relationship between the genetic determinants of sarcoidosis and other inflammatory disorders has not been characterized. Herein, we examined the degree of genetic pleiotropy common to sarcoidosis and other inflammatory disorders to identify shared pathways and disease systems pertinent to sarcoidosis onset. To achieve this, we quantify the association of common variant polygenic risk scores from nine complex inflammatory disorders with sarcoidosis risk. Enrichment analyses of genes implicated in pleiotropic associations were further used to elucidate candidate pathways. In European-Americans, we identify significant pleiotropy between risk of sarcoidosis and risk of asthma (R2=2.03%; P=8.89 × 10−9), celiac disease (R2=2.03%; P=8.21 × 10−9), primary biliary cirrhosis (R2=2.43%; P=2.01 × 10−10) and rheumatoid arthritis (R2=4.32%; P=2.50 × 10−17). These associations validate in African Americans only after accounting for the proportion of genome-wide European ancestry, where we demonstrate similar effects of polygenic risk for African-Americans with the highest levels of European ancestry. Variants and genes implicated in European-American pleiotropic associations were enriched for pathways involving interleukin-12, interleukin-27 and cell adhesion molecules, corroborating the hypothesized immunopathogenesis of disease.

The work by Tang et al. provides a comprehensive, single-cell, transcriptomic analysis of kidney and blood cells from the adult zebrafish, identifying novel cell types, including two classes of NK immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells., Recent advances in single-cell, transcriptomic profiling have provided unprecedented access to investigate cell heterogeneity during tissue and organ development. In this study, we used massively parallel, single-cell RNA sequencing to define cell heterogeneity within the zebrafish kidney marrow, constructing a comprehensive molecular atlas of definitive hematopoiesis and functionally distinct renal cells found in adult zebrafish. Because our method analyzed blood and kidney cells in an unbiased manner, our approach was useful in characterizing immune-cell deficiencies within DNA–protein kinase catalytic subunit (prkdc), interleukin-2 receptor γ a (il2rga), and double-homozygous–mutant fish, identifying blood cell losses in T, B, and natural killer cells within specific genetic mutants. Our analysis also uncovered novel cell types, including two classes of natural killer immune cells, classically defined and erythroid-primed hematopoietic stem and progenitor cells, mucin-secreting kidney cells, and kidney stem/progenitor cells. In total, our work provides the first, comprehensive, single-cell, transcriptomic analysis of kidney and marrow cells in the adult zebrafish.

Two studies reveal a role for the transcriptional regulator BATF3 as a T cell–intrinsic factor mediating effective memory responses. This finding opens future avenues of investigation and opportunities to enhance cellular immunotherapy.

Our ability to track cellular dynamics in humans over time in vivo has been limited. Here, we demonstrate how somatic mutations in mitochondrial DNA (mtDNA) can be used to longitudinally track the dynamic output of hematopoietic stem and progenitor cells in humans. Over the course of 3 years of blood sampling in a single individual, our analyses reveal somatic mtDNA sequence variation and evolution reminiscent of models of hematopoiesis established by genetic labeling approaches. Furthermore, we observe fluctuations in mutation heteroplasmy, coinciding with specific clinical events, such as infections, and further identify lineage-specific somatic mtDNA mutations in longitudinally sampled circulating blood cell subsets in individuals with leukemia. Collectively, these observations indicate the significant potential of using tracking of somatic mtDNA sequence variation as a broadly applicable approach to systematically assess hematopoietic clonal dynamics in human health and disease.

Myeloproliferative neoplasms (MPNs) are blood cancers characterized by excessive production of mature myeloid cells that result from the acquisition of somatic driver mutations in hematopoietic stem cells (HSCs)1. While substantial progress has been made to define the causal somatic mutation profile for MPNs2, epidemiologic studies indicate a significant heritable component for the disease that is among the highest known for all cancers3. However, only a limited set of genetic risk loci have been identified, and the underlying biological mechanisms leading to MPN acquisition remain unexplained. Here, to define the inherited risk profile, we conducted the largest genome-wide association study of MPNs to date (978,913 individuals with 3,224 cases) and identified 14 genome-wide significant loci, as well as a polygenic signature that increases the odds for disease acquisition by nearly 3-fold between the top and median deciles. Interestingly, we find a shared genetic architecture between MPN risk and several hematopoietic traits spanning distinct lineages, as well as an association between increased MPN risk and longer leukocyte telomere length, collectively implicating HSC function and self-renewal. Strikingly, we find a significant enrichment for risk variants mapping to accessible chromatin in HSCs compared with other hematopoietic populations. Finally, gene mapping identifies modulators of HSC biology and targeted variant-to-function analyses suggest likely roles for CHEK2 and GFI1B in altering HSC function to confer disease risk. Overall, we demonstrate the power of human genetic studies to illuminate a previously unappreciated mechanism for MPN risk through modulation of HSC function.

Myeloproliferative neoplasms (MPNs) are blood cancers that are characterized by the excessive production of mature myeloid cells and arise from the acquisition of somatic driver mutations in haematopoietic stem cells (HSCs). Epidemiological studies indicate a substantial heritable component of MPNs that is among the highest known for cancers1. However, only a limited number of genetic risk loci have been identified, and the underlying biological mechanisms that lead to the acquisition of MPNs remain unclear. Here, by conducting a large-scale genome-wide association study (3,797 cases and 1,152,977 controls), we identify 17 MPN risk loci (P

Background Recent innovations in single-cell Assay for Transposase Accessible Chromatin using sequencing (scATAC-seq) enable profiling of the epigenetic landscape of thousands of individual cells. scATAC-seq data analysis presents unique methodological challenges. scATAC-seq experiments sample DNA, which, due to low copy numbers (diploid in humans), lead to inherent data sparsity (1–10% of peaks detected per cell) compared to transcriptomic (scRNA-seq) data (10–45% of expressed genes detected per cell). Such challenges in data generation emphasize the need for informative features to assess cell heterogeneity at the chromatin level. Results We present a benchmarking framework that is applied to 10 computational methods for scATAC-seq on 13 synthetic and real datasets from different assays, profiling cell types from diverse tissues and organisms. Methods for processing and featurizing scATAC-seq data were compared by their ability to discriminate cell types when combined with common unsupervised clustering approaches. We rank evaluated methods and discuss computational challenges associated with scATAC-seq analysis including inherently sparse data, determination of features, peak calling, the effects of sequencing coverage and noise, and clustering performance. Running times and memory requirements are also discussed. Conclusions This reference summary of scATAC-seq methods offers recommendations for best practices with consideration for both the non-expert user and the methods developer. Despite variation across methods and datasets, SnapATAC, Cusanovich2018, and cisTopic outperform other methods in separating cell populations of different coverages and noise levels in both synthetic and real datasets. Notably, SnapATAC is the only method able to analyze a large dataset (> 80,000 cells).

Studies of allelic variation underlying genetic blood disorders have provided important insights into human hematopoiesis. Most often, the identified pathogenic mutations result in loss-of-function or missense changes. However, assessing the pathogenicity of noncoding variants can be challenging. Here, we characterize two unrelated patients with a distinct presentation of dyserythropoietic anemia and other impairments in hematopoiesis associated with an intronic mutation in GATA1 that is 24 nucleotides upstream of the canonical splice acceptor site. Functional studies demonstrate that this single-nucleotide alteration leads to reduced canonical splicing and increased use of an alternative splice acceptor site that causes a partial intron retention event. The resultant altered GATA1 contains a five–amino acid insertion at the C-terminus of the C-terminal zinc finger and has no observable activity. Collectively, our results demonstrate how altered splicing of GATA1, which reduces levels of the normal form of this master transcription factor, can result in distinct changes in human hematopoiesis.

Lineage tracing provides key insights into the fate of individual cells in complex organisms. Although effective genetic labeling approaches are available in model systems, in humans, most approaches require detection of nuclear somatic mutations, which have high error rates, limited scale, and do not capture cell state information. Here, we show that somatic mutations in mtDNA can be tracked by single-cell RNA or assay for transposase accessible chromatin (ATAC) sequencing. We leverage somatic mtDNA mutations as natural genetic barcodes and demonstrate their utility as highly accurate clonal markers to infer cellular relationships. We track native human cells both in vitro and in vivo and relate clonal dynamics to gene expression and chromatin accessibility. Our approach should allow clonal tracking at a 1,000-fold greater scale than with nuclear genome sequencing, with simultaneous information on cell state, opening the way to chart cellular dynamics in human health and disease. Using single-cell sequencing technologies, somatic mutations in mtDNA can be used as natural genetic barcodes to study cellular states and clonal dynamics., National Institutes of Health (U.S.) (Grant R33 CA202820)

Enhancer elements in the human genome control how genes are expressed in specific cell types and harbor thousands of genetic variants that influence risk for common diseases1-4. Yet, we still do not know how enhancers regulate specific genes, and we lack general rules to predict enhancer-gene connections across cell types5,6. We developed an experimental approach, CRISPRi-FlowFISH, to perturb enhancers in the genome, and we applied it to test >3,500 potential enhancer-gene connections for 30 genes. We found that a simple activity-by-contact model substantially outperformed previous methods at predicting the complex connections in our CRISPR dataset. This activity-by-contact model allows us to construct genome-wide maps of enhancer-gene connections in a given cell type, on the basis of chromatin state measurements. Together, CRISPRi-FlowFISH and the activity-by-contact model provide a systematic approach to map and predict which enhancers regulate which genes, and will help to interpret the functions of the thousands of disease risk variants in the noncoding genome.

The core cohesin subunit STAG2 is recurrently mutated in Ewing sarcoma but its biological role is less clear. Herein, we demonstrate that cohesin complexes containing STAG2 occupy enhancer and polycomb repressive complex (PRC2) marked regulatory regions. Genetic suppression of STAG2 leads to a compensatory increase in cohesin-STAG1 complexes, but not in enhancer rich regions, and results in reprogramming of cis-chromatin interactions. Strikingly, in STAG2 knockout cells, the oncogenic genetic program driven by the fusion transcription factor EWS/FLI1 was highly perturbed, in part due to altered enhancer-promoter contacts. Moreover, loss of STAG2 also disrupted PRC2-mediated regulation of gene expression. Combined, these transcriptional changes converged to modulate EWS/FLI1, migratory and neurodevelopmental programs. Finally, consistent with clinical observations, functional studies revealed that loss of STAG2 enhances the metastatic potential of Ewing sarcoma xenografts. Our findings demonstrate that STAG2 mutations can alter chromatin architecture and transcriptional programs to promote an aggressive cancer phenotype.

Children diagnosed with diffuse intrinsic pontine glioma (DIPG), a type of high grade glioma in the brainstem, currently have a dismal 5-year overall survival of only 2%. The majority of DIPG patients harbor a K27M mutation in histone 3.3 encoding genes (H3.3K27M). To understand if the aberrant epigenetic landscape induced by H3.3K27M provides an opportunity for novel targeted therapies, we conducted the first CRISPR/Cas9 screen using a focused library of 1,350 epigenetic regulatory and cancer related genes in six H3.3K27M DIPG patient-derived primary neurosphere cell lines. We identified gene dependencies in chromatin regulators, polycomb repressive complexes 1 and 2 (PRC1 and PRC2), histone demethylases, acetyltransferases and deacetylators as novel tumor cell dependencies in DIPG. We hypothesized that targeting dysregulated functions of chromatin regulators by genetically deleting and chemically targeting these epigenetically induced vulnerabilities, we could ameliorate, or even reverse the downstream oncogenic effects of the aberrant epigenetic landscape of DIPG. In our secondary CRISPR nanoscreen, we first used six single guide RNAs (sgRNA) to knockout each gene using CRISPR/Cas9 ribonucleoprotein nucleofections, followed by use of three best sgRNAs combined with homology directed repair templates. Compared to lentiviral delivery, nucleofection is a rapid method, with reduced off-target toxicity, suitable for single gene knockouts in DIPG neurospheres. Secondary CRISPR validations confirmed dependencies in BMI1, CBX4, KDM1A, EZH2, EED, SUZ12, HDAC2, and EP300. Next, we conducted a chemical screen using 20 inhibitors and degraders to target the aberrant activity of HDAC, KDM1A, P300/CBP, PRC1 and PRC2. We identified eight chemical compounds that were effective in H3.3K27M DIPG neurosphere cell lines at low drug concentrations. Among these, an inhibitor and degrader targeting P300/CBP activity indicates a novel strategy of epigenetic therapy in DIPG. Through our combinatorial testing, we will identify a synergistic combination of epigenetic therapy for treating children diagnosed with H3.3K27M DIPG.

Widespread changes to DNA methylation and chromatin are well documented in cancer, but the fate of higher-order chromosomal structure remains obscure. Here we integrated topological maps for colon tumors and normal colons with epigenetic, transcriptional, and imaging data to characterize alterations to chromatin loops, topologically associated domains, and large-scale compartments. We found that spatial partitioning of the open and closed genome compartments is profoundly compromised in tumors. This reorganization is accompanied by compartment-specific hypomethylation and chromatin changes. Additionally, we identify a compartment at the interface between the canonical A and B compartments that is reorganized in tumors. Remarkably, similar shifts were evident in non-malignant cells that have accumulated excess divisions. Our analyses suggest that these topological changes repress stemness and invasion programs while inducing anti-tumor immunity genes and may therefore restrain malignant progression. Our findings call into question the conventional view that tumor-associated epigenomic alterations are primarily oncogenic.

A hallmark of cancer development is the loss of cellular identity due to the dysregulation of regulatory networks that maintain functional, differentiated states. Chromatin state has been linked to this control of cellular identity and developmental processes; however, the mechanisms by which these regulatory landscapes are disrupted during cancer progression are not well understood. To this end, we leveraged an optimized single-cell ATAC-sequencing (scATAC-seq) strategy using a combinatorial indexing-based technology to assess chromatin state changes that occur in a mouse model of lung adenocarcinoma (LUAD). We utilized a well-established LUAD mouse model in which mice conditionally express the KrasG12D mutation and lose p53 expression in alveolar type II cells following Cre-mediated recombination, termed the KP model. In this model, KP tumors progress from early-stage hyperplasias to metastatic, advanced adenocarcinomas without acquiring additional somatic mutations, suggesting that these phenotypic transitions may be driven in part by epigenetic mechanisms. Using scATAC-seq, we profiled 16,044 normal lung, primary KP tumor, and metastatic cells to assess chromatin state changes across LUAD progression. We found that LUAD tumor evolution was characterized by a heterogeneous epigenomic continuum that was associated with loss of alveolar identity and progression to metastasis. Interestingly, we identified a reproducible metastatic state that was homogeneous, suggesting that cancer cells funnel towards a stable epigenetic state. In addition, we utilized novel computational tools to assess combinatorial transcription factor (TF)-driven programs that delineate stages of cancer progression in the KP model. From these analyses, we identified a clear transition point that corresponded to progressive gain of RUNX2 transcription factor activity. Using CRISPR activation and knockout strategies, we found that RUNX2 regulated the chromatin accessibility of several regulatory programs involved in extracellular matrix remodeling, suggesting that RUNX2 signaling is a critical event in LUAD progression and metastasis. In addition, gene signatures derived from this RUNX-activated state were highly predictive of survival in human LUAD. Together, these results demonstrated the power of single-cell epigenomics to identify TF-driven regulatory programs as key biomarkers of tumor progression. This abstract is also being presented as Poster A40. Citation Format: Lindsay M. LaFave, Vinay Kartha, Sai Ma, Kevin Meli, Isabella Del Priore, Caleb Lareau, Venkat Sanker, Santiago Naranjo, Peter Westcott, Zachary Chiang, Alison Brack, Travis Law, Aviv Regev, Jason D Buenrostro, Tyler Jacks. Leveraging single-cell epigenomics to uncover regulatory programs in lung adenocarcinoma [abstract]. In: Proceedings of the AACR Special Conference on the Evolving Landscape of Cancer Modeling; 2020 Mar 2-5; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2020;80(11 Suppl):Abstract nr PR08.

Genome-wide association studies (GWAS) have identified thousands of variants associated with human diseases and traits. However, the majority of GWAS-implicated variants are in non-coding regions of the genome and require in depth follow-up to identify target genes and decipher biological mechanisms. Here, rather than focusing on causal variants, we have undertaken a pooled loss-of-function screen in primary hematopoietic cells to interrogate 389 candidate genes contained in 75 loci associated with red blood cell traits. Using this approach, we identify 77 genes at 38 GWAS loci, with most loci harboring 1–2 candidate genes. Importantly, the hit set was strongly enriched for genes validated through orthogonal genetic approaches. Genes identified by this approach are enriched in specific and relevant biological pathways, allowing regulators of human erythropoiesis and modifiers of blood diseases to be defined. More generally, this functional screen provides a paradigm for gene-centric follow up of GWAS for a variety of human diseases and traits.

CRISPR-Cas base-editor technology enables targeted nucleotide alterations, and is being increasingly used for research and potential therapeutic applications1,2. The most widely used cytosine base editors (CBEs) induce deamination of DNA cytosines using the rat APOBEC1 enzyme, which is targeted by a linked Cas protein-guide RNA complex3,4. Previous studies of the specificity of CBEs have identified off-target DNA edits in mammalian cells5,6. Here we show that a CBE with rat APOBEC1 can cause extensive transcriptome-wide deamination of RNA cytosines in human cells, inducing tens of thousands of C-to-U edits with frequencies ranging from 0.07% to 100% in 38-58% of expressed genes. CBE-induced RNA edits occur in both protein-coding and non-protein-coding sequences and generate missense, nonsense, splice site, and 5' and 3' untranslated region mutations. We engineered two CBE variants bearing mutations in rat APOBEC1 that substantially decreased the number of RNA edits (by more than 390-fold and more than 3,800-fold) in human cells. These variants also showed more precise on-target DNA editing than the wild-type CBE and, for most guide RNAs tested, no substantial reduction in editing efficiency. Finally, we show that an adenine base editor7 can also induce transcriptome-wide RNA edits. These results have implications for the use of base editors in both research and clinical settings, illustrate the feasibility of engineering improved variants with reduced RNA editing activities, and suggest the need to more fully define and characterize the RNA off-target effects of deaminase enzymes in base editor platforms.

Summary A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

Schaefer et al. recently advanced the provocative conclusion that CRISPR-Cas9 nuclease can induce off-target alterations at genomic loci that do not resemble the intended on-target site. Using high-coverage whole genome sequencing (WGS), these authors reported finding SNPs and indels in two CRISPR-Cas9-treated mice that were not present in a single untreated control mouse. On the basis of this association, Schaefer et al. concluded that these sequence variants were caused by CRISPR-Cas9. This new proposed CRISPR-Cas9 off-target activity runs contrary to previously published work and, if the authors are correct, could have profound implications for research and therapeutic applications. Here we demonstrate that the simplest interpretation of Schaefer et al.9s data is that the two CRISPR-Cas9-treated mice are actually more closely related genetically to each other than to the control mouse. This strongly suggests that the so-called “unexpected mutations” simply represent SNPs and indels shared in common by these mice prior to nuclease treatment. In addition, given the genomic and sequence distribution profiles of these variants, we show that it is challenging to explain how CRISPR-Cas9 might be expected to induce such changes. Finally, we argue that the lack of appropriate controls in Schaefer et al.9s experimental design precludes assignment of causality to CRISPR-Cas9. Given these substantial issues, we urge Schaefer et al. to revise or re-state the original conclusions of their published work so as to avoid leaving misleading and unsupported statements to persist in the literature.

Incomplete annotation of cell-to-cell state variance and widespread linkage disequilibrium in the human genome represent significant challenges to elucidating mechanisms of trait-associated genetic variation. Here, using data from the UK Biobank, we perform genetic fine-mapping for 16 blood cell traits to quantify posterior probabilities of association while allowing for multiple independent signals per region. We observe an enrichment of fine-mapped variants in accessible chromatin of lineage-committed hematopoietic progenitor cells. Further, we develop a novel analytic framework that identifies “core gene” cell type enrichments and show that this approach uniquely resolves relevant cell types within closely related populations. Applying our approach to single cell chromatin accessibility data, we discover significant heterogeneity within classically defined multipotential progenitor populations. Finally, using several lines of empirical evidence, we identify relevant cell types, predict target genes, and propose putative causal mechanisms for fine-mapped variants. In total, our study provides an analytic framework for single-variant and single-cell analyses to elucidate putative causal variants and cell types from GWAS and high-resolution epigenomic assays.

Human hematopoiesis involves cellular differentiation of multipotent cells into progressively more lineage-restricted states. While the chromatin accessibility landscape of this process has been explored in defined populations, single-cell regulatory variation has been hidden by ensemble averaging. We collected single-cell chromatin accessibility profiles across 10 populations of immunophenotypically defined human hematopoietic cell types and constructed a chromatin accessibility landscape of human hematopoiesis to characterize differentiation trajectories. We find variation consistent with lineage bias toward different developmental branches in multipotent cell types. We observe heterogeneity within common myeloid progenitors (CMPs) and granulocyte-macrophage progenitors (GMPs) and develop a strategy to partition GMPs along their differentiation trajectory. Furthermore, we integrated single-cell RNA sequencing (scRNA-seq) data to associate transcription factors to chromatin accessibility changes and regulatory elements to target genes through correlations of expression and regulatory element accessibility. Overall, this work provides a framework for integrative exploration of complex regulatory dynamics in a primary human tissue at single-cell resolution.

Schaefer et al. recently advanced the provocative conclusion that CRISPR-Cas9 nuclease can induce off-target alterations at genomic loci that do not resemble the intended on-target site.1 Using high-coverage whole genome sequencing (WGS), these authors reported finding SNPs and indels in two CRISPR-Cas9-treated mice that were not present in a single untreated control mouse. On the basis of this association, Schaefer et al. concluded that these sequence variants were caused by CRISPR-Cas9. This new proposed CRISPR-Cas9 off-target activity runs contrary to previously published work2–8 and, if the authors are correct, could have profound implications for research and therapeutic applications. Here, we demonstrate that the simplest interpretation of Schaefer et al.’s data is that the two CRISPR-Cas9-treated mice are actually more closely related genetically to each other than to the control mouse. This strongly suggests that the so-called “unexpected mutations” simply represent SNPs and indels shared in common by these mice prior to nuclease treatment. In addition, given the genomic and sequence distribution profiles of these variants, we show that it is challenging to explain how CRISPR-Cas9 might be expected to induce such changes. Finally, we argue that the lack of appropriate controls in Schaefer et al.’s experimental design precludes assignment of causality to CRISPR-Cas9. Given these substantial issues, we urge Schaefer et al. to revise or re-state the original conclusions of their published work so as to avoid leaving misleading and unsupported statements to persist in the literature.