Mapping genotypes to chromatin accessibility profiles in single cells.

In somatic tissue differentiation, chromatin accessibility changes govern priming and precursor commitment towards cellular fates1-3. Therefore, somatic mutations are likely to alter chromatin accessibility patterns, as they disrupt differentiation topologies leading to abnormal clonal outgrowth. However, defining the impact of somatic mutations on the epigenome in human samples is challenging due to admixed mutated and wild-type cells. Here, to chart how somatic mutations disrupt epigenetic landscapes in human clonal outgrowths, we developed genotyping of targeted loci with single-cell chromatin accessibility (GoT-ChA). This high-throughput platform links genotypes to chromatin accessibility at single-cell resolution across thousands of cells within a single assay. We applied GoT-ChA to CD34+ cells from patients with myeloproliferative neoplasms with JAK2V617F-mutated haematopoiesis. Differential accessibility analysis between wild-type and JAK2V617F-mutant progenitors revealed both cell-intrinsic and cell-state-specific shifts within mutant haematopoietic precursors, including cell-intrinsic pro-inflammatory signatures in haematopoietic stem cells, and a distinct profibrotic inflammatory chromatin landscape in megakaryocytic progenitors. Integration of mitochondrial genome profiling and cell-surface protein expression measurement allowed expansion of genotyping onto DOGMA-seq through imputation, enabling single-cell capture of genotypes, chromatin accessibility, RNA expression and cell-surface protein expression. Collectively, we show that the JAK2V617F mutation leads to epigenetic rewiring in a cell-intrinsic and cell type-specific manner, influencing inflammation states and differentiation trajectories. We envision that GoT-ChA will empower broad future investigations of the critical link between somatic mutations and epigenetic alterations across clonal populations in malignant and non-malignant contexts.

Meiotic DNA break resection and recombination rely on chromatin remodeler Fun30.

DNA double-strand breaks (DSBs) are nucleolytically processed to generate single-stranded DNA tails for homologous recombination. In Saccharomyces cerevisiae meiosis, this 5'-to-3' resection involves initial nicking by the Mre11-Rad50-Xrs2 complex (MRX) plus Sae2, then exonucleolytic digestion by Exo1. Chromatin remodeling adjacent to meiotic DSBs is thought to be necessary for resection, but the relevant remodeling activity was unknown. Here we show that the SWI/SNF-like ATPase Fun30 plays a major, non-redundant role in resecting meiotic DSBs. A fun30 null mutation shortened resection tract lengths almost as severely as an exo1-nd (nuclease-dead) mutation, and resection was further shortened in the fun30 exo1-nd double mutant. Fun30 associates with chromatin in response to meiotic DSBs, and the constitutive positioning of nucleosomes governs resection endpoint locations in the absence of Fun30. We infer that Fun30 directly promotes both the MRX- and Exo1-dependent steps in resection, possibly by removing nucleosomes from broken chromatids. Moreover, we found that the extremely short resection in the fun30 exo1-nd double mutant is accompanied by compromised interhomolog recombination bias, leading to defects in recombination and chromosome segregation. Thus, this study also provides insight about the minimal resection lengths needed for robust recombination.

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

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

Discovery of target genes and pathways at GWAS loci by pooled single-cell CRISPR screens.

Most variants associated with complex traits and diseases identified by genome-wide association studies (GWAS) map to noncoding regions of the genome with unknown effects. Using ancestrally diverse, biobank-scale GWAS data, massively parallel CRISPR screens, and single-cell transcriptomic and proteomic sequencing, we discovered 124 cis-target genes of 91 noncoding blood trait GWAS loci. Using precise variant insertion through base editing, we connected specific variants with gene expression changes. We also identified trans-effect networks of noncoding loci when cis target genes encoded transcription factors or microRNAs. Networks were themselves enriched for GWAS variants and demonstrated polygenic contributions to complex traits. This platform enables massively parallel characterization of the target genes and mechanisms of human noncoding variants in both cis and trans.

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

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

mRNA COVID-19 vaccine elicits potent adaptive immune response without the persistent inflammation seen in SARS-CoV-2 infection.

SARS-CoV-2 infection and vaccination elicit potent immune responses. Our study presents a comprehensive multimodal single-cell dataset of peripheral blood of patients with acute COVID-19 and of healthy volunteers before and after receiving the SARS-CoV-2 mRNA vaccine and booster. We compared host immune responses to the virus and vaccine using transcriptional profiling, coupled with B/T cell receptor repertoire reconstruction. COVID-19 patients displayed an enhanced interferon signature and cytotoxic gene upregulation, absent in vaccine recipients. These findings were validated in an independent dataset. Analysis of B and T cell repertoires revealed that, while the majority of clonal lymphocytes in COVID-19 patients were effector cells, clonal expansion was more evident among circulating memory cells in vaccine recipients. Furthermore, while clonal αß T cell responses were observed in both COVID-19 patients and vaccine recipients, dramatic expansion of clonal γδT cells was found only in infected individuals. Our dataset enables comparative analyses of immune responses to infection versus vaccination, including clonal B and T cell responses. Integrating our data with publicly available datasets allowed us to validate our findings in larger cohorts. To our knowledge, this is the first dataset to include comprehensive profiling of longitudinal samples from healthy volunteers pre/post SARS-CoV-2 vaccine and booster.

Efficient combinatorial targeting of RNA transcripts in single cells with Cas13 RNA Perturb-seq.

Pooled CRISPR screens coupled with single-cell RNA-sequencing have enabled systematic interrogation of gene function and regulatory networks. Here, we introduce Cas13 RNA Perturb-seq (CaRPool-seq), which leverages the RNA-targeting CRISPR-Cas13d system and enables efficient combinatorial perturbations alongside multimodal single-cell profiling. CaRPool-seq encodes multiple perturbations on a cleavable CRISPR array that is associated with a detectable barcode sequence, allowing for the simultaneous targeting of multiple genes. We compared CaRPool-seq to existing Cas9-based methods, highlighting its unique strength to efficiently profile combinatorially perturbed cells. Finally, we apply CaRPool-seq to perform multiplexed combinatorial perturbations of myeloid differentiation regulators in an acute myeloid leukemia (AML) model system and identify extensive interactions between different chromatin regulators that can enhance or suppress AML differentiation phenotypes.

mRNA COVID-19 vaccine elicits potent adaptive immune response without the acute inflammation of SARS-CoV-2 infection.

SARS-CoV-2 infection and vaccination elicit potent immune responses. Our study presents a comprehensive multimodal single-cell analysis of blood from COVID-19 patients and healthy volunteers receiving the SARS-CoV-2 vaccine and booster. We profiled immune responses via transcriptional analysis and lymphocyte repertoire reconstruction. COVID-19 patients displayed an enhanced interferon signature and cytotoxic gene upregulation, absent in vaccine recipients. B and T cell repertoire analysis revealed clonal expansion among effector cells in COVID-19 patients and memory cells in vaccine recipients. Furthermore, while clonal αß T cell responses were observed in both COVID-19 patients and vaccine recipients, expansion of clonal γδ T cells was found only in infected individuals. Our dataset enables side-by-side comparison of immune responses to infection versus vaccination, including clonal B and T cell responses. Our comparative analysis shows that vaccination induces a robust, durable clonal B and T cell responses, without the severe inflammation associated with infection.

A genome-scale screen for synthetic drivers of T cell proliferation.

The engineering of autologous patient T cells for adoptive cell therapies has revolutionized the treatment of several types of cancer1. However, further improvements are needed to increase response and cure rates. CRISPR-based loss-of-function screens have been limited to negative regulators of T cell functions2-4 and raise safety concerns owing to the permanent modification of the genome. Here we identify positive regulators of T cell functions through overexpression of around 12,000 barcoded human open reading frames (ORFs). The top-ranked genes increased the proliferation and activation of primary human CD4+ and CD8+ T cells and their secretion of key cytokines such as interleukin-2 and interferon-γ. In addition, we developed the single-cell genomics method OverCITE-seq for high-throughput quantification of the transcriptome and surface antigens in ORF-engineered T cells. The top-ranked ORF-lymphotoxin-ß receptor (LTBR)-is typically expressed in myeloid cells but absent in lymphocytes. When overexpressed in T cells, LTBR induced profound transcriptional and epigenomic remodelling, leading to increased T cell effector functions and resistance to exhaustion in chronic stimulation settings through constitutive activation of the canonical NF-κB pathway. LTBR and other highly ranked genes improved the antigen-specific responses of chimeric antigen receptor T cells and γδ T cells, highlighting their potential for future cancer-agnostic therapies5. Our results provide several strategies for improving next-generation T cell therapies by the induction of synthetic cell programmes.

Multimodal single-cell analysis of cutaneous T-cell lymphoma reveals distinct subclonal tissue-dependent signatures.

Cutaneous T-cell lymphoma (CTCL) is a heterogeneous group of mature T-cell neoplasms characterized by the accumulation of clonal malignant CD4+ T cells in the skin. The most common variant of CTCL, mycosis fungoides (MF ), is confined to the skin in early stages but can be accompanied by extracutaneous dissemination of malignant T cells to the blood and lymph nodes in advanced stages of disease. Sézary syndrome (SS), a leukemic form of disease, is characterized by significant blood involvement. Little is known about the transcriptional and genomic relationship between skin- and blood-residing malignant T cells in CTCL. To identify and interrogate malignant clones in matched skin and blood from patients with leukemic MF and SS, we combine T-cell receptor clonotyping with quantification of gene expression and cell surface markers at the single cell level. Our data reveal clonal evolution at a transcriptional and genetic level within the malignant populations of individual patients. We highlight highly consistent transcriptional signatures delineating skin- and blood-derived malignant T cells. Analysis of these 2 populations suggests that environmental cues, along with genetic aberrations, contribute to transcriptional profiles of malignant T cells. Our findings indicate that the skin microenvironment in CTCL promotes a transcriptional response supporting rapid malignant expansion, as opposed to the quiescent state observed in the blood, potentially influencing efficacy of therapies. These results provide insight into tissue-specific characteristics of cancerous cells and underscore the need to address the patients' individual malignant profiles at the time of therapy to eliminate all subclones.

Integrated analysis of multimodal single-cell data.

The simultaneous measurement of multiple modalities represents an exciting frontier for single-cell genomics and necessitates computational methods that can define cellular states based on multimodal data. Here, we introduce "weighted-nearest neighbor" analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of 211,000 human peripheral blood mononuclear cells (PBMCs) with panels extending to 228 antibodies to construct a multimodal reference atlas of the circulating immune system. Multimodal analysis substantially improves our ability to resolve cell states, allowing us to identify and validate previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets and to interpret immune responses to vaccination and coronavirus disease 2019 (COVID-19). Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets and to look beyond the transcriptome toward a unified and multimodal definition of cellular identity.

Scalable, multimodal profiling of chromatin accessibility, gene expression and protein levels in single cells.

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 molecular regulation of inhibitory immune checkpoints with multimodal single-cell screens.

The expression of inhibitory immune checkpoint molecules, such as programmed death-ligand (PD-L)1, is frequently observed in human cancers and can lead to the suppression of T cell-mediated immune responses. Here, we apply expanded CRISPR-compatible (EC)CITE-seq, a technology that combines pooled CRISPR screens with single-cell mRNA and surface protein measurements, to explore the molecular networks that regulate PD-L1 expression. We also develop a computational framework, mixscape, that substantially improves the signal-to-noise ratio in single-cell perturbation screens by identifying and removing confounding sources of variation. Applying these tools, we identify and validate regulators of PD-L1 and leverage our multimodal data to identify both transcriptional and post-transcriptional modes of regulation. Specifically, we discover that the Kelch-like protein KEAP1 and the transcriptional activator NRF2 mediate the upregulation of PD-L1 after interferon (IFN)-γ stimulation. Our results identify a new mechanism for the regulation of immune checkpoints and present a powerful analytical framework for the analysis of multimodal single-cell perturbation screens.

Structural and functional characterization of the Spo11 core complex.

Spo11, which makes DNA double-strand breaks (DSBs) that are essential for meiotic recombination, has long been recalcitrant to biochemical study. We provide molecular analysis of Saccharomyces cerevisiae Spo11 purified with partners Rec102, Rec104 and Ski8. Rec102 and Rec104 jointly resemble the B subunit of archaeal topoisomerase VI, with Rec104 occupying a position similar to the Top6B GHKL-type ATPase domain. Unexpectedly, the Spo11 complex is monomeric (1:1:1:1 stoichiometry), consistent with dimerization controlling DSB formation. Reconstitution of DNA binding reveals topoisomerase-like preferences for duplex-duplex junctions and bent DNA. Spo11 also binds noncovalently but with high affinity to DNA ends mimicking cleavage products, suggesting a mechanism to cap DSB ends. Mutations that reduce DNA binding in vitro attenuate DSB formation, alter DSB processing and reshape the DSB landscape in vivo. Our data reveal structural and functional similarities between the Spo11 core complex and Topo VI, but also highlight differences reflecting their distinct biological roles.

Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells.

To better understand host-virus genetic dependencies and find potential therapeutic targets for COVID-19, we performed a genome-scale CRISPR loss-of-function screen to identify host factors required for SARS-CoV-2 viral infection of human alveolar epithelial cells. Top-ranked genes cluster into distinct pathways, including the vacuolar ATPase proton pump, Retromer, and Commander complexes. We validate these gene targets using several orthogonal methods such as CRISPR knockout, RNA interference knockdown, and small-molecule inhibitors. Using single-cell RNA-sequencing, we identify shared transcriptional changes in cholesterol biosynthesis upon loss of top-ranked genes. In addition, given the key role of the ACE2 receptor in the early stages of viral entry, we show that loss of RAB7A reduces viral entry by sequestering the ACE2 receptor inside cells. Overall, this work provides a genome-scale, quantitative resource of the impact of the loss of each host gene on fitness/response to viral infection.

Multiplexed detection of proteins, transcriptomes, clonotypes and CRISPR perturbations in single cells.

Multimodal single-cell assays provide high-resolution snapshots of complex cell populations, but are mostly limited to transcriptome plus an additional modality. Here, we describe expanded CRISPR-compatible cellular indexing of transcriptomes and epitopes by sequencing (ECCITE-seq) for the high-throughput characterization of at least five modalities of information from each single cell. We demonstrate application of ECCITE-seq to multimodal CRISPR screens with robust direct single-guide RNA capture and to clonotype-aware multimodal phenotyping of cancer samples.

S1-seq Assay for Mapping Processed DNA Ends.

During meiosis, the specialized cell division giving rise to gametes, numerous DNA double-strand breaks (DSBs) are introduced at multiple places throughout the genome by the topoisomerase-like protein Spo11. Homologous recombination, a highly conserved DSB repair pathway, is employed for their repair and ensures the formation of chiasmata and the proper segregation of homologous chromosomes. In the initial steps of recombination, end resection takes place, wherein Spo11 is endonucleolytically released and the 5'-terminal strands of each DSB are exonucleolytically processed, exposing the ssDNA necessary to identify a homologous repair template. DNA removed by DSB processing is reconstituted by DNA synthesis, which copies genetic information from the intact homologous template. We developed a next-generation sequencing assay, termed S1-seq, to study DNA end resection genome-wide at high-spatial resolution during yeast meiotic recombination. The assay relies on the fact that removal of the ssDNA tails of resected DSBs marks the position where resection stopped. Molecular features of resection are revealed by sequencing of these ssDNA-to-dsDNA junctions and comparison to high-resolution Spo11 DSB maps. We describe the experimental and computational methods for S1-seq as applied to meiosis in the SK1 strain of budding yeast Saccharomyces cerevisiae and discuss how it can also be applied to map DSBs and recombination intermediates.

A global view of meiotic double-strand break end resection.

DNA double-strand breaks that initiate meiotic recombination are exonucleolytically processed. This 5'→3' resection is a central, conserved feature of recombination but remains poorly understood. To address this lack, we mapped resection endpoints genome-wide at high resolution in Saccharomyces cerevisiae Full-length resection requires Exo1 exonuclease and the DSB-responsive kinase Tel1, but not Sgs1 helicase. Tel1 also promotes efficient and timely resection initiation. Resection endpoints display pronounced heterogeneity between genomic loci that reflects a tendency for nucleosomes to block Exo1, yet Exo1 also appears to digest chromatin with high processivity and at rates similar to naked DNA in vitro. This paradox points to nucleosome destabilization or eviction as a defining feature of the meiotic resection landscape.

Resection activity of the Sgs1 helicase alters the affinity of DNA ends for homologous recombination proteins in Saccharomyces cerevisiae.

The RecQ helicase family is critical during DNA damage repair, and mutations in these proteins are associated with Bloom, Werner, or Rothmund-Thompson syndromes in humans, leading to cancer predisposition and/or premature aging. In the budding yeast Saccharomyces cerevisiae, mutations in the RecQ homolog, SGS1, phenocopy many of the defects observed in the human syndromes. One challenge to studying RecQ helicases is that their disruption leads to a pleiotropic phenotype. Using yeast, we show that the separation-of-function allele of SGS1, sgs1-D664Δ, has impaired activity at DNA ends, resulting in a resection processivity defect. Compromising Sgs1 resection function in the absence of the Sae2 nuclease causes slow growth, which is alleviated by making the DNA ends accessible to Exo1 nuclease. Furthermore, fluorescent microscopy studies reveal that, when Sgs1 resection activity is compromised in sae2Δ cells, Mre11 repair foci persist. We suggest a model where the role of Sgs1 in end resection along with Sae2 is important for removing Mre11 from DNA ends during repair.