Multiplexed Engineering and Analysis of Combinatorial Enhancer Activity in Single Cells

Deciphering the impact of genomic variation on function

Our genomes influence nearly every aspect of human biology-from molecular and cellular functions to phenotypes in health and disease. Studying the differences in DNA sequence between individuals (genomic variation) could reveal previously unknown mechanisms of human biology, uncover the basis of genetic predispositions to diseases, and guide the development of new diagnostic tools and therapeutic agents. Yet, understanding how genomic variation alters genome function to influence phenotype has proved challenging. To unlock these insights, we need a systematic and comprehensive catalogue of genome function and the molecular and cellular effects of genomic variants. Towards this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations and predictive modelling to investigate the relationships among genomic variation, genome function and phenotypes. IGVF will create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how such effects connect through gene-regulatory and protein-interaction networks. These experimental data, computational predictions and accompanying standards and pipelines will be integrated into an open resource that will catalyse community efforts to explore how our genomes influence biology and disease across populations.

Toward a foundation model of causal cell and tissue biology with a Perturbation Cell and Tissue Atlas

Comprehensively charting the biologically causal circuits that govern the phenotypic space of human cells has often been viewed as an insurmountable challenge. However, in the last decade, a suite of interleaved experimental and computational technologies has arisen that is making this fundamental goal increasingly tractable. Pooled CRISPR-based perturbation screens with high-content molecular and/or image-based readouts are now enabling researchers to probe, map, and decipher genetically causal circuits at increasing scale. This scale is now eminently suitable for the deployment of artificial intelligence and machine learning (AI/ML) to both direct further experiments and to predict or generate information that was not-and sometimes cannot-be gathered experimentally. By combining and iterating those through experiments that are designed for inference, we now envision a Perturbation Cell Atlas as a generative causal foundation model to unify human cell biology.

Oncogenic transcription factors instruct promoter-enhancer hubs in individual triple negative breast cancer cells

Sequencing-based mapping of ensemble pairwise interactions among regulatory elements support the existence of topological assemblies known as promoter-enhancer hubs or cliques in cancer. Yet, prevalence, regulators, and functions of promoter-enhancer hubs in individual cancer cells remain unclear. Here, we systematically integrated functional genomics, transcription factor screening, and optical mapping of promoter-enhancer interactions to identify key promoter-enhancer hubs, examine heterogeneity of their assembly, determine their regulators, and elucidate their role in gene expression control in individual triple negative breast cancer (TNBC) cells. Optical mapping of individual SOX9 and MYC alleles revealed the existence of frequent multiway interactions among promoters and enhancers within spatial hubs. Our single-allele studies further demonstrated that lineage-determining SOX9 and signaling-dependent NOTCH1 transcription factors compact MYC and SOX9 hubs. Together, our findings suggest that promoter-enhancer hubs are dynamic and heterogeneous topological assemblies, which are controlled by oncogenic transcription factors and facilitate subtype-restricted gene expression in cancer.

Comprehensive transcription factor perturbations recapitulate fibroblast transcriptional states

Cell atlas projects have nominated recurrent transcriptional states as drivers of biological processes and disease, but their origins, regulation, and properties remain unclear. To enable complementary functional studies, we developed a scalable approach for recapitulating cell states in vitro using CRISPR activation (CRISPRa) Perturb-seq. Aided by a novel multiplexing method, we activated 1,836 transcription factors in two cell types. Measuring 21,958 perturbations showed that CRISPRa activated targets within physiological ranges, that epigenetic features predicted activatable genes, and that the protospacer seed region drove an off-target effect. Perturbations recapitulated in vivo fibroblast states, including universal and inflammatory states, and identified KLF4 and KLF5 as key regulators of the universal state. Inducing the universal state suppressed disease-associated states, highlighting its therapeutic potential. Our findings cement CRISPRa as a tool for perturbing differentiated cells and indicate that in vivo states can be elicited via perturbation, enabling studies of clinically relevant states ex vivo.

Unlocking biological mechanisms with integrative functional genomics approaches

Reverse genetics offers precise functional insights into genes through the targeted manipulation of gene expression followed by phenotypic assessment. While these approaches have proven effective in model organisms such as Saccharomyces cerevisiae, large-scale genetic manipulations in human cells were historically unfeasible due to methodological limitations. However, recent advancements in functional genomics, particularly clustered regularly interspaced short palindromic repeats (CRISPR)-based screening technologies and next-generation sequencing platforms, have enabled pooled screening technologies that allow massively parallel, unbiased assessments of biological phenomena in human cells. This review provides a comprehensive overview of cutting-edge functional genomic screening technologies applicable to human cells, ranging from short hairpin RNA screens to modern CRISPR screens. Additionally, we explore the integration of CRISPR platforms with single-cell approaches to monitor gene expression, chromatin accessibility, epigenetic regulation, and chromatin architecture following genetic perturbations at the omics level. By offering an in-depth understanding of these genomic screening methods, this review aims to provide insights into more targeted and effective strategies for genomic research and personalized medicine.

Dissecting key regulators of transcriptome kinetics through scalable single-cell RNA profiling of pooled CRISPR screens

We present a combinatorial indexing method, PerturbSci-Kinetics, for capturing whole transcriptomes, nascent transcriptomes and single guide RNA (sgRNA) identities across hundreds of genetic perturbations at the single-cell level. Profiling a pooled CRISPR screen targeting various biological processes, we show the gene expression regulation during RNA synthesis, processing and degradation, miRNA biogenesis and mitochondrial mRNA processing, systematically decoding the genome-wide regulatory network that underlies RNA temporal dynamics at scale.

Multiome Perturb-seq unlocks scalable discovery of integrated perturbation effects on the transcriptome and epigenome

Single-cell CRISPR screens link genetic perturbations to transcriptional states, but high-throughput methods connecting these induced changes to their regulatory foundations are limited. Here we introduce Multiome Perturb-seq, extending single-cell CRISPR screens to simultaneously measure perturbation-induced changes in gene expression and chromatin accessibility. We apply Multiome Perturb-seq in a CRISPRi screen of 13 chromatin remodelers in human RPE-1 cells, achieving efficient assignment of sgRNA identities to single nuclei via an improved method for capturing barcode transcripts from nuclear RNA. We organize expression and accessibility measurements into coherent programs describing the integrated effects of perturbations on cell state, finding that ARID1A and SUZ12 knockdowns induce programs enriched for developmental features. Pseudotime analysis of perturbations connects accessibility changes to changes in gene expression, highlighting the value of multimodal profiling. Overall, our method provides a scalable and simply implemented system to dissect the regulatory logic underpinning cell state.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Cohesin-mediated 3D contacts tune enhancer-promoter regulation

Enhancers are key drivers of gene regulation thought to act via 3D physical interactions with the promoters of their target genes. However, genome-wide depletions of architectural proteins such as cohesin result in only limited changes in gene expression, despite a loss of contact domains and loops. Consequently, the role of cohesin and 3D contacts in enhancer function remains debated. Here, we developed CRISPRi of regulatory elements upon degron operation (CRUDO), a novel approach to measure how changes in contact frequency impact enhancer effects on target genes by perturbing enhancers with CRISPRi and measuring gene expression in the presence or absence of cohesin. We systematically perturbed all 1,039 candidate enhancers near five cohesin-dependent genes and identified 34 enhancer-gene regulatory interactions. Of 26 regulatory interactions with sufficient statistical power to evaluate cohesin dependence, 18 show cohesin-dependent effects. A decrease in enhancer-promoter contact frequency upon removal of cohesin is frequently accompanied by a decrease in the regulatory effect of the enhancer on gene expression, consistent with a contact-based model for enhancer function. However, changes in contact frequency and regulatory effects on gene expression vary as a function of distance, with distal enhancers (e.g., >50Kb) experiencing much larger changes than proximal ones (e.g., <50Kb). Because most enhancers are located close to their target genes, these observations can explain how only a small subset of genes - those with strong distal enhancers - are sensitive to cohesin. Together, our results illuminate how 3D contacts, influenced by both cohesin and genomic distance, tune enhancer effects on gene expression.

SpRY-mediated screens facilitate functional dissection of non-coding sequences at single-base resolution

CRISPR mutagenesis screens conducted with SpCas9 and other nucleases have identified certain cis-regulatory elements and genetic variants but at a limited resolution due to the absence of protospacer adjacent motif (PAM) sequences. Here, leveraging the broad targeting scope of the near-PAMless SpRY variant, we have demonstrated that saturated SpRY mutagenesis and base editing screens can faithfully identify functional regulatory elements and essential genetic variants for target gene expression at single-base resolution. We further extended this methodology to investigate a genome-wide association study (GWAS) locus at 10q22.1 associated with a red blood cell trait, where we identified potential enhancers regulating HK1 gene expression, despite not all of these enhancers exhibiting typical chromatin signatures. More importantly, our saturated base editing screens pinpoint multiple causal variants within this locus that would otherwise be missed by Bayesian statistical fine-mapping. Our approach is generally applicable to functional interrogation of all non-coding genomic elements while complementing other high-coverage CRISPR screens.

Analysis of single-cell CRISPR perturbations indicates that enhancers act multiplicatively and provides limited evidence for epistatic-like interactions

A single gene may have multiple enhancers, but how they work in concert to regulate transcription is poorly understood. To analyze enhancer interactions throughout the genome, we developed a generalized linear modeling framework, GLiMMIRS, for interrogating enhancer effects from single-cell CRISPR experiments. We applied GLiMMIRS to a published dataset and tested for interactions between 46,166 enhancer pairs and corresponding genes, including 264 'high-confidence' enhancer pairs. We found that enhancer effects combine multiplicatively but with limited evidence for further interactions. Only 31 enhancer pairs exhibited significant interactions (FDR < 0.1), of which none came from the high confidence subset and 20 were driven by outlier expression values. Additional analyses of a second CRISPR dataset and in silico enhancer perturbations with Enformer both support a multiplicative model of enhancer effects without interactions. Altogether, our results indicate that enhancer interactions are uncommon or have small effects that are difficult to detect.

Advances in computational and experimental approaches for deciphering transcriptional regulatory networks: Understanding the roles of cis-regulatory elements is essential, and recent research utilizing MPRAs, STARR-seq, CRISPR-Cas9, and machine learning has yielded valuable insights

Understanding the influence of cis-regulatory elements on gene regulation poses numerous challenges given complexities stemming from variations in transcription factor (TF) binding, chromatin accessibility, structural constraints, and cell-type differences. This review discusses the role of gene regulatory networks in enhancing understanding of transcriptional regulation and covers construction methods ranging from expression-based approaches to supervised machine learning. Additionally, key experimental methods, including MPRAs and CRISPR-Cas9-based screening, which have significantly contributed to understanding TF binding preferences and cis-regulatory element functions, are explored. Lastly, the potential of machine learning and artificial intelligence to unravel cis-regulatory logic is analyzed. These computational advances have far-reaching implications for precision medicine, therapeutic target discovery, and the study of genetic variations in health and disease.

Defining cis-regulatory elements and transcription factors that control human cortical interneuron development

Although human cortical interneurons (cINs) are a minority population in the cerebral cortex, disruption of interneuron development is a frequent contributor to neurodevelopmental disorders. Here, we utilized a model for deriving cINs from human embryonic stem cells to profile chromatin state changes and generate an atlas of cis-regulatory elements (CREs) controlling human cIN development. We used these data to define candidate transcription factors (TFs) that may bind these CREs to regulate interneuron progenitor specification. Among these were RFX3 and RFX4, risk genes for autism spectrum disorder (ASD) with uncharacterized roles in human neuronal development. Using RFX3 and RFX4 knockdown models, we demonstrated new requirements for both genes in interneuron progenitor specification, with RFX3 deficiency causing precocious neuronal differentiation while RFX4 deficiency instead resulted in cessation of progenitor cell proliferation. Together, this work both defined central features of cis-regulatory control and identified new TF requirements for human interneuron development.

Targeting T cell plasticity in kidney and gut inflammation by pooled single-cell CRISPR screening

Pro-inflammatory CD4+ T cells are major drivers of autoimmune diseases, yet therapies modulating T cell phenotypes to promote an anti-inflammatory state are lacking. Here, we identify T helper 17 (TH17) cell plasticity in the kidneys of patients with antineutrophil cytoplasmic antibody-associated glomerulonephritis on the basis of single-cell (sc) T cell receptor analysis and scRNA velocity. To uncover molecules driving T cell polarization and plasticity, we established an in vivo pooled scCRISPR droplet sequencing (iCROP-seq) screen and applied it to mouse models of glomerulonephritis and colitis. CRISPR-based gene targeting in TH17 cells could be ranked according to the resulting transcriptional perturbations, and polarization biases into T helper 1 (TH1) and regulatory T cells could be quantified. Furthermore, we show that iCROP-seq can facilitate the identification of therapeutic targets by efficient functional stratification of genes and pathways in a disease- and tissue-specific manner. These findings uncover TH17 to TH1 cell plasticity in the human kidney in the context of renal autoimmunity.

Decoding polygenic diseases: advances in noncoding variant prioritization and validation

Genome-wide association studies (GWASs) provide a key foundation for elucidating the genetic underpinnings of common polygenic diseases. However, these studies have limitations in their ability to assign causality to particular genetic variants, especially those residing in the noncoding genome. Over the past decade, technological and methodological advances in both analytical and empirical prioritization of noncoding variants have enabled the identification of causative variants by leveraging orthogonal functional evidence at increasing scale. In this review, we present an overview of these approaches and describe how this workflow provides the groundwork necessary to move beyond associations toward genetically informed studies on the molecular and cellular mechanisms of polygenic disease.

The ETO2 transcriptional cofactor maintains acute leukemia by driving a MYB/EP300-dependent stemness program

Transcriptional cofactors of the ETO family are recurrent fusion partners in acute leukemia. We characterized the ETO2 regulome by integrating transcriptomic and chromatin binding analyses in human erythroleukemia xenografts and controlled ETO2 depletion models. We demonstrate that beyond its well-established repressive activity, ETO2 directly activates transcription of MYB, among other genes. The ETO2-activated signature is associated with a poorer prognosis in erythroleukemia but also in other acute myeloid and lymphoid leukemia subtypes. Mechanistically, ETO2 colocalizes with EP300 and MYB at enhancers supporting the existence of an ETO2/MYB feedforward transcription activation loop (e.g., on MYB itself). Both small-molecule and PROTAC-mediated inhibition of EP300 acetyltransferases strongly reduced ETO2 protein, chromatin binding, and ETO2-activated transcripts. Taken together, our data show that ETO2 positively enforces a leukemia maintenance program that is mediated in part by the MYB transcription factor and that relies on acetyltransferase cofactors to stabilize ETO2 scaffolding activity.

Bidirectional epigenetic editing reveals hierarchies in gene regulation

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

Robust differential expression testing for single-cell CRISPR screens at low multiplicity of infection

Single-cell CRISPR screens (perturb-seq) link genetic perturbations to phenotypic changes in individual cells. The most fundamental task in perturb-seq analysis is to test for association between a perturbation and a count outcome, such as gene expression. We conduct the first-ever comprehensive benchmarking study of association testing methods for low multiplicity-of-infection (MOI) perturb-seq data, finding that existing methods produce excess false positives. We conduct an extensive empirical investigation of the data, identifying three core analysis challenges: sparsity, confounding, and model misspecification. Finally, we develop an association testing method - SCEPTRE low-MOI - that resolves these analysis challenges and demonstrates improved calibration and power.

Multiplex, single-cell CRISPRa screening for cell type specific regulatory elements

CRISPR-based gene activation (CRISPRa) is a promising therapeutic approach for gene therapy, upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis- regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to a library of 493 gRNAs targeting candidate cis- regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically upregulating intended target genes and no other neighboring genes within 1 Mb, including gRNAs yielding upregulation of six autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes in neurons. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans- acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate therapeutically relevant genes in a cell type-specific manner.

Robust differential expression testing for single-cell CRISPR screens at low multiplicity of infection

Single-cell CRISPR screens (perturb-seq) link genetic perturbations to phenotypic changes in individual cells. The most fundamental task in perturb-seq analysis is to test for association between a perturbation and a count outcome, such as gene expression. We conduct the first-ever comprehensive benchmarking study of association testing methods for low multiplicity-of-infection (MOI) perturb-seq data, finding that existing methods produce excess false positives. We conduct an extensive empirical investigation of the data, identifying three core analysis challenges: sparsity, confounding, and model misspecification. Finally, we develop an association testing method - SCEPTRE low-MOI - that resolves these analysis challenges and demonstrates improved calibration and power.

Recurrent mutations in tumor suppressor FBXW7 bypass Wnt/β-catenin addiction in cancer

Pathologic Wnt/β-catenin signaling drives various cancers, leading to multiple approaches to drug this pathway. Appropriate patient selection can maximize success of these interventions. Wnt ligand addiction is a druggable vulnerability in RNF43-mutant/RSPO-fusion cancers. However, pharmacologically targeting the biogenesis of Wnt ligands, e.g., with PORCN inhibitors, has shown mixed therapeutic responses, possibly due to tumor heterogeneity. Here, we show that the tumor suppressor FBXW7 is frequently mutated in RNF43-mutant/RSPO-fusion tumors, and FBXW7 mutations cause intrinsic resistance to anti-Wnt therapies. Mechanistically, FBXW7 inactivation stabilizes multiple oncoproteins including Cyclin E and MYC and antagonizes the cytostatic effect of Wnt inhibitors. Moreover, although FBXW7 mutations do not mitigate β-catenin degradation upon Wnt inhibition, FBXW7-mutant RNF43-mutant/RSPO-fusion cancers instead lose dependence on β-catenin signaling, accompanied by dedifferentiation and loss of lineage specificity. These FBXW7-mutant Wnt/β-catenin-independent tumors are susceptible to multi-cyclin-dependent kinase inhibition. An in-depth understanding of primary resistance to anti-Wnt/β-catenin therapies allows for more appropriate patient selection and use of alternative mechanism-based therapies.

The rise of single-cell transcriptomics in yeast

The field of single-cell omics has transformed our understanding of biological processes and is constantly advancing both experimentally and computationally. One of the most significant developments is the ability to measure the transcriptome of individual cells by single-cell RNA-seq (scRNA-seq), which was pioneered in higher eukaryotes. While yeast has served as a powerful model organism in which to test and develop transcriptomic technologies, the implementation of scRNA-seq has been significantly delayed in this organism, mainly because of technical constraints associated with its intrinsic characteristics, namely the presence of a cell wall, a small cell size and little amounts of RNA. In this review, we examine the current technologies for scRNA-seq in yeast and highlight their strengths and weaknesses. Additionally, we explore opportunities for developing novel technologies and the potential outcomes of implementing single-cell transcriptomics and extension to other modalities. Undoubtedly, scRNA-seq will be invaluable for both basic and applied yeast research, providing unique insights into fundamental biological processes.

Applications of single‑cell omics and spatial transcriptomics technologies in gastric cancer (Review)

Gastric cancer (GC) is a prominent contributor to global cancer-related mortalities, and a deeper understanding of its molecular characteristics and tumor heterogeneity is required. Single-cell omics and spatial transcriptomics (ST) technologies have revolutionized cancer research by enabling the exploration of cellular heterogeneity and molecular landscapes at the single-cell level. In the present review, an overview of the advancements in single-cell omics and ST technologies and their applications in GC research is provided. Firstly, multiple single-cell omics and ST methods are discussed, highlighting their ability to offer unique insights into gene expression, genetic alterations, epigenomic modifications, protein expression patterns and cellular location in tissues. Furthermore, a summary is provided of key findings from previous research on single-cell omics and ST methods used in GC, which have provided valuable insights into genetic alterations, tumor diagnosis and prognosis, tumor microenvironment analysis, and treatment response. In summary, the application of single-cell omics and ST technologies has revealed the levels of cellular heterogeneity and the molecular characteristics of GC, and holds promise for improving diagnostics, personalized treatments and patient outcomes in GC.

Enhancer-promoter specificity in gene transcription: molecular mechanisms and disease associations

Although often located at a distance from their target gene promoters, enhancers are the primary genomic determinants of temporal and spatial transcriptional specificity in metazoans. Since the discovery of the first enhancer element in simian virus 40, there has been substantial interest in unraveling the mechanism(s) by which enhancers communicate with their partner promoters to ensure proper gene expression. These research efforts have benefited considerably from the application of increasingly sophisticated sequencing- and imaging-based approaches in conjunction with innovative (epi)genome-editing technologies; however, despite various proposed models, the principles of enhancer-promoter interaction have still not been fully elucidated. In this review, we provide an overview of recent progress in the eukaryotic gene transcription field pertaining to enhancer-promoter specificity. A better understanding of the mechanistic basis of lineage- and context-dependent enhancer-promoter engagement, along with the continued identification of functional enhancers, will provide key insights into the spatiotemporal control of gene expression that can reveal therapeutic opportunities for a range of enhancer-related diseases.

GWAS-informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density

Over 1,100 independent signals have been identified with genome-wide association studies (GWAS) for bone mineral density (BMD), a key risk factor for mortality-increasing fragility fractures; however, the effector gene(s) for most remain unknown. Informed by a variant-to-gene mapping strategy implicating 89 non-coding elements predicted to regulate osteoblast gene expression at BMD GWAS loci, we executed a single-cell CRISPRi screen in human fetal osteoblast 1.19 cells (hFOBs). The BMD relevance of hFOBs was supported by heritability enrichment from cross-cell type stratified LD-score regression involving 98 cell types grouped into 15 tissues. 24 genes showed perturbation in the screen, with four (ARID5B, CC2D1B, EIF4G2, and NCOA3) exhibiting consistent effects upon siRNA knockdown on three measures of osteoblast maturation and mineralization. Lastly, additional heritability enrichments, genetic correlations, and multi-trait fine-mapping revealed that many BMD GWAS signals are pleiotropic and likely mediate their effects via non-bone tissues that warrant attention in future screens.

scPerturb: harmonized single-cell perturbation data

Analysis across a growing number of single-cell perturbation datasets is hampered by poor data interoperability. To facilitate development and benchmarking of computational methods, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform quality control pipelines and harmonize feature annotations. The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets. We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing, and demonstrate E-distance as a general distance measure between sets of single-cell expression profiles. We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations. The perturbation-response datasets and E-statistics computation software are publicly available at scperturb.org. This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.

Past, present, and future of CRISPR genome editing technologies

Genome editing has been a transformative force in the life sciences and human medicine, offering unprecedented opportunities to dissect complex biological processes and treat the underlying causes of many genetic diseases. CRISPR-based technologies, with their remarkable efficiency and easy programmability, stand at the forefront of this revolution. In this Review, we discuss the current state of CRISPR gene editing technologies in both research and therapy, highlighting limitations that constrain them and the technological innovations that have been developed in recent years to address them. Additionally, we examine and summarize the current landscape of gene editing applications in the context of human health and therapeutics. Finally, we outline potential future developments that could shape gene editing technologies and their applications in the coming years.

No country for old methods: New tools for studying microproteins

Microproteins encoded by small open reading frames (sORFs) have emerged as a fascinating frontier in genomics. Traditionally overlooked due to their small size, recent technological advancements such as ribosome profiling, mass spectrometry-based strategies and advanced computational approaches have led to the annotation of more than 7000 sORFs in the human genome. Despite the vast progress, only a tiny portion of these microproteins have been characterized and an important challenge in the field lies in identifying functionally relevant microproteins and understanding their role in different cellular contexts. In this review, we explore the recent advancements in sORF research, focusing on the new methodologies and computational approaches that have facilitated their identification and functional characterization. Leveraging these new tools hold great promise for dissecting the diverse cellular roles of microproteins and will ultimately pave the way for understanding their role in the pathogenesis of diseases and identifying new therapeutic targets.

Higher-order combinatorial chromatin perturbations by engineered CRISPR-Cas12a for functional genomics

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting 1-3 genomic sites per cell. To develop a tool for higher-order ( > 3) combinatorial targeting of genomic sites with CRISPRi in functional genomics screens, we engineered an Acidaminococcus Cas12a variant -- referred to as mul tiplexed transcriptional interference AsCas12a (multiAsCas12a). multiAsCas12a incorporates a key mutation, R1226A, motivated by the hypothesis of nicking-induced stabilization of the ribonucleoprotein:DNA complex for improving CRISPRi activity. multiAsCas12a significantly outperforms prior state-of-the-art Cas12a variants in combinatorial CRISPRi targeting using high-order multiplexed arrays of lentivirally transduced CRISPR RNAs (crRNA), including in high-throughput pooled screens using 6-plex crRNA array libraries. Using multiAsCas12a CRISPRi, we discover new enhancer elements and dissect the combinatorial function of cis-regulatory elements. These results instantiate a group testing framework for efficiently surveying potentially numerous combinations of chromatin perturbations for biological discovery and engineering.

Mapping the functional impact of non-coding regulatory elements in primary T cells through single-cell CRISPR screens

BACKGROUND

Drug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions.

RESULTS

To identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci.

CONCLUSIONS

We describe "primary T cell crisprQTL" - a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.

Next-generation forward genetic screens: uniting high-throughput perturbations with single-cell analysis

Programmable genome-engineering technologies, such as CRISPR (clustered regularly interspaced short palindromic repeats) nucleases and massively parallel CRISPR screens that capitalize on this programmability, have transformed biomedical science. These screens connect genes and noncoding genome elements to disease-relevant phenotypes, but until recently have been limited to individual phenotypes such as growth or fluorescent reporters of gene expression. By pairing massively parallel screens with high-dimensional profiling of single-cell types/states, we can now measure how individual genetic perturbations or combinations of perturbations impact the cellular transcriptome, proteome, and epigenome. We review technologies that pair CRISPR screens with single-cell multiomics and the unique opportunities afforded by extending pooled screens using deep multimodal phenotyping.

Emerging applications of single-cell profiling in precision medicine of atherosclerosis

Atherosclerosis is a chronic, progressive, inflammatory disease that occurs in the arterial wall. Despite recent advancements in treatment aimed at improving efficacy and prolonging survival, atherosclerosis remains largely incurable. In this review, we discuss emerging single-cell sequencing techniques and their novel insights into atherosclerosis. We provide examples of single-cell profiling studies that reveal phenotypic characteristics of atherosclerosis plaques, blood, liver, and the intestinal tract. Additionally, we highlight the potential clinical applications of single-cell analysis and propose that combining this approach with other techniques can facilitate early diagnosis and treatment, leading to more accurate medical interventions.

scDecouple: decoupling cellular response from infected proportion bias in scCRISPR-seq

Single-cell clustered regularly interspaced short palindromic repeats-sequencing (scCRISPR-seq) is an emerging high-throughput CRISPR screening technology where the true cellular response to perturbation is coupled with infected proportion bias of guide RNAs (gRNAs) across different cell clusters. The mixing of these effects introduces noise into scCRISPR-seq data analysis and thus obstacles to relevant studies. We developed scDecouple to decouple true cellular response of perturbation from the influence of infected proportion bias. scDecouple first models the distribution of gene expression profiles in perturbed cells and then iteratively finds the maximum likelihood of cell cluster proportions as well as the cellular response for each gRNA. We demonstrated its performance in a series of simulation experiments. By applying scDecouple to real scCRISPR-seq data, we found that scDecouple enhances the identification of biologically perturbation-related genes. scDecouple can benefit scCRISPR-seq data analysis, especially in the case of heterogeneous samples or complex gRNA libraries.

The Current Situation and Development Prospect of Whole-Genome Screening

High-throughput genetic screening is useful for discovering critical genes or gene sequences that trigger specific cell functions and/or phenotypes. Loss-of-function genetic screening is mainly achieved through RNA interference (RNAi), CRISPR knock-out (CRISPRko), and CRISPR interference (CRISPRi) technologies. Gain-of-function genetic screening mainly depends on the overexpression of a cDNA library and CRISPR activation (CRISPRa). Base editing can perform both gain- and loss-of-function genetic screening. This review discusses genetic screening techniques based on Cas9 nuclease, including Cas9-mediated genome knock-out and dCas9-based gene activation and interference. We compare these methods with previous genetic screening techniques based on RNAi and cDNA library overexpression and propose future prospects and applications for CRISPR screening.

Advancements in CRISPR screens for the development of cancer immunotherapy strategies

CRISPR screen technology enables systematic and scalable interrogation of gene function by using the CRISPR-Cas9 system to perturb gene expression. In the field of cancer immunotherapy, this technology has empowered the discovery of genes, biomarkers, and pathways that regulate tumor development and progression, immune reactivity, and the effectiveness of immunotherapeutic interventions. By conducting large-scale genetic screens, researchers have successfully identified novel targets to impede tumor growth, enhance anti-tumor immune responses, and surmount immunosuppression within the tumor microenvironment (TME). Here, we present an overview of CRISPR screens conducted in tumor cells for the purpose of identifying novel therapeutic targets. We also explore the application of CRISPR screens in immune cells to propel the advancement of cell-based therapies, encompassing T cells, natural killer cells, dendritic cells, and macrophages. Furthermore, we outline the crucial components necessary for the successful implementation of immune-specific CRISPR screens and explore potential directions for future research.

Decoding Heterogenous Single-cell Perturbation Responses

Understanding diverse responses of individual cells to the same perturbation is central to many biological and biomedical problems. Current methods, however, do not precisely quantify the strength of perturbation responses and, more importantly, reveal new biological insights from heterogeneity in responses. Here we introduce the perturbation-response score (PS), based on constrained quadratic optimization, to quantify diverse perturbation responses at a single-cell level. Applied to single-cell transcriptomes of large-scale genetic perturbation datasets (e.g., Perturb-seq), PS outperforms existing methods for quantifying partial gene perturbation responses. In addition, PS presents two major advances. First, PS enables large-scale, single-cell-resolution dosage analysis of perturbation, without the need to titrate perturbation strength. By analyzing the dose-response patterns of over 2,000 essential genes in Perturb-seq, we identify two distinct patterns, depending on whether a moderate reduction in their expression induces strong downstream expression alterations. Second, PS identifies intrinsic and extrinsic biological determinants of perturbation responses. We demonstrate the application of PS in contexts such as T cell stimulation, latent HIV-1 expression, and pancreatic cell differentiation. Notably, PS unveiled a previously unrecognized, cell-type-specific role of coiled-coil domain containing 6 (CCDC6) in guiding liver and pancreatic lineage decisions, where CCDC6 knockouts drive the endoderm cell differentiation towards liver lineage, rather than pancreatic lineage. The PS approach provides an innovative method for dose-to-function analysis and will enable new biological discoveries from single-cell perturbation datasets.

One sentence summary

We present a method to quantify diverse perturbation responses and discover novel biological insights in single-cell perturbation datasets.

Pooled Genome-Scale CRISPR Screens in Single Cells

Assigning functions to genes and learning how to control their expression are part of the foundation of cell biology and therapeutic development. An efficient and unbiased method to accomplish this is genetic screening, which historically required laborious clone generation and phenotyping and is still limited by scale today. The rapid technological progress on modulating gene function with CRISPR-Cas and measuring it in individual cells has now relaxed the major experimental constraints and enabled pooled screening with complex readouts from single cells. Here, we review the principles and practical considerations for pooled single-cell CRISPR screening. We discuss perturbation strategies, experimental model systems, matching the perturbation to the individual cells, reading out cell phenotypes, and data analysis. Our focus is on single-cell RNA sequencing and cell sorting-based readouts, including image-enabled cell sorting. We expect this transformative approach to fuel biomedical research for the next several decades.

Enhancer regulatory networks globally connect non-coding breast cancer loci to cancer genes

Genetic studies have associated thousands of enhancers with breast cancer. However, the vast majority have not been functionally characterized. Thus, it remains unclear how variant-associated enhancers contribute to cancer. Here, we perform single-cell CRISPRi screens of 3,512 regulatory elements associated with breast cancer to measure the impact of these regions on transcriptional phenotypes. Analysis of >500,000 single-cell transcriptomes in two breast cancer cell lines shows that perturbation of variant-associated enhancers disrupts breast cancer gene programs. We observe variant-associated enhancers that directly or indirectly regulate the expression of cancer genes. We also find one-to-multiple and multiple-to-one network motifs where enhancers indirectly regulate cancer genes. Notably, multiple variant-associated enhancers indirectly regulate TP53. Comparative studies illustrate sub-type specific functions between enhancers in ER+ and ER- cells. Finally, we developed the pySpade package to facilitate analysis of single-cell enhancer screens. Overall, we demonstrate that enhancers form regulatory networks that link cancer genes in the genome, providing a more comprehensive understanding of the contribution of enhancers to breast cancer development.

CRISPR-based epigenome editing: mechanisms and applications

Epigenomic anomalies contribute significantly to the development of numerous human disorders. The development of epigenetic research tools is essential for understanding how epigenetic marks contribute to gene expression. A gene-editing technique known as CRISPR (clustered regularly interspaced short palindromic repeats) typically targets a particular DNA sequence using a guide RNA (gRNA). CRISPR/Cas9 technology has been remodeled for epigenome editing by generating a 'dead' Cas9 protein (dCas9) that lacks nuclease activity and juxtaposing it with an epigenetic effector domain. Based on fusion partners of dCas9, a specific epigenetic state can be achieved. CRISPR-based epigenome editing has widespread application in drug screening, cancer treatment and regenerative medicine. This paper discusses the tools developed for CRISPR-based epigenome editing and their applications.

Applications of CRISPR technology in cellular immunotherapy

CRISPR technology has transformed multiple fields, including cancer and immunology. CRISPR-based gene editing and screening empowers direct genomic manipulation of immune cells, opening doors to unbiased functional genetic screens. These screens aid in the discovery of novel factors that regulate and reprogram immune responses, offering novel drug targets. The engineering of immune cells using CRISPR has sparked a transformation in the cellular immunotherapy field, resulting in a multitude of ongoing clinical trials. In this review, we discuss the development and applications of CRISPR and related gene editing technologies in immune cells, focusing on functional genomics screening, gene editing-based cell therapies, as well as future directions in this rapidly advancing field.

Gene regulatory network inference in the era of single-cell multi-omics

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

In-organoid single-cell CRISPR screening reveals determinants of hepatocyte differentiation and maturation

BACKGROUND

Harnessing hepatocytes for basic research and regenerative medicine demands a complete understanding of the genetic determinants underlying hepatocyte differentiation and maturation. Single-cell CRISPR screens in organoids could link genetic perturbations with parallel transcriptomic readout in single cells, providing a powerful method to delineate roles of cell fate regulators. However, a big challenge for identifying key regulators during data analysis is the low expression levels of transcription factors (TFs), which are difficult to accurately estimate due to noise and dropouts in single-cell sequencing. Also, it is often the changes in TF activities in the transcriptional cascade rather than the expression levels of TFs that are relevant to the cell fate transition.

RESULTS

Here, we develop Organoid-based Single-cell CRISPR screening Analyzed with Regulons (OSCAR), a framework using regulon activities as readouts to dissect gene knockout effects in organoids. In adult-stem-cell-derived liver organoids, we map transcriptomes in 80,576 cells upon 246 perturbations associated with transcriptional regulation of hepatocyte formation. Using OSCAR, we identify known and novel positive and negative regulators, among which Fos and Ubr5 are the top-ranked ones. Further single-gene loss-of-function assays demonstrate that Fos depletion in mouse and human liver organoids promote hepatocyte differentiation by specific upregulation of liver metabolic genes and pathways, and conditional knockout of Ubr5 in mouse liver delays hepatocyte maturation.

CONCLUSIONS

Altogether, we provide a framework to explore lineage specifiers in a rapid and systematic manner, and identify hepatocyte determinators with potential clinical applications.

CRISPR screening in hematology research: from bulk to single-cell level

The CRISPR genome editing technology has revolutionized the way gene function is studied. Genome editing can be achieved in single genes or for thousands of genes simultaneously in sensitive genetic screens. While conventional genetic screens are limited to bulk measurements of cell behavior, recent developments in single-cell technologies make it possible to combine CRISPR screening with single-cell profiling. In this way, cell behavior and gene expression can be monitored simultaneously, with the additional possibility of including data on chromatin accessibility and protein levels. Moreover, the availability of various Cas proteins leading to inactivation, activation, or other effects on gene function further broadens the scope of such screens. The integration of single-cell multi-omics approaches with CRISPR screening open the path to high-content information on the impact of genetic perturbations at single-cell resolution. Current limitations in cell throughput and data density need to be taken into consideration, but new technologies are rapidly evolving and are likely to easily overcome these limitations. In this review, we discuss the use of bulk CRISPR screening in hematology research, as well as the emergence of single-cell CRISPR screening and its added value to the field.

Gene regulatory network reconstruction: harnessing the power of single-cell multi-omic data

Inferring gene regulatory networks (GRNs) is a fundamental challenge in biology that aims to unravel the complex relationships between genes and their regulators. Deciphering these networks plays a critical role in understanding the underlying regulatory crosstalk that drives many cellular processes and diseases. Recent advances in sequencing technology have led to the development of state-of-the-art GRN inference methods that exploit matched single-cell multi-omic data. By employing diverse mathematical and statistical methodologies, these methods aim to reconstruct more comprehensive and precise gene regulatory networks. In this review, we give a brief overview on the statistical and methodological foundations commonly used in GRN inference methods. We then compare and contrast the latest state-of-the-art GRN inference methods for single-cell matched multi-omics data, and discuss their assumptions, limitations and opportunities. Finally, we discuss the challenges and future directions that hold promise for further advancements in this rapidly developing field.

The technological landscape and applications of single-cell multi-omics

Single-cell multi-omics technologies and methods characterize cell states and activities by simultaneously integrating various single-modality omics methods that profile the transcriptome, genome, epigenome, epitranscriptome, proteome, metabolome and other (emerging) omics. Collectively, these methods are revolutionizing molecular cell biology research. In this comprehensive Review, we discuss established multi-omics technologies as well as cutting-edge and state-of-the-art methods in the field. We discuss how multi-omics technologies have been adapted and improved over the past decade using a framework characterized by optimization of throughput and resolution, modality integration, uniqueness and accuracy, and we also discuss multi-omics limitations. We highlight the impact that single-cell multi-omics technologies have had in cell lineage tracing, tissue-specific and cell-specific atlas production, tumour immunology and cancer genetics, and in mapping of cellular spatial information in fundamental and translational research. Finally, we discuss bioinformatics tools that have been developed to link different omics modalities and elucidate functionality through the use of better mathematical modelling and computational methods.

Multiple genes in a single GWAS risk locus synergistically mediate aberrant synaptic development and function in human neurons

The mechanistic tie between genome-wide association study (GWAS)-implicated risk variants and disease-relevant cellular phenotypes remains largely unknown. Here, using human induced pluripotent stem cell (hiPSC)-derived neurons as a neurodevelopmental model, we identify multiple schizophrenia (SZ) risk variants that display allele-specific open chromatin (ASoC) and are likely to be functional. Editing the strongest ASoC SNP, rs2027349, near vacuolar protein sorting 45 homolog (VPS45) alters the expression of VPS45, lncRNA AC244033.2, and a distal gene, C1orf54. Notably, the transcriptomic changes in neurons are associated with SZ and other neuropsychiatric disorders. Neurons carrying the risk allele exhibit increased dendritic complexity and hyperactivity. Interestingly, individual/combinatorial gene knockdown shows that these genes alter cellular phenotypes in a non-additive synergistic manner. Our study reveals that multiple genes at a single GWAS risk locus mediate a compound effect on neural function, providing a mechanistic link between a non-coding risk variant and disease-related cellular phenotypes.

Manipulating and studying gene function in human pluripotent stem cell models

Human pluripotent stem cells (hPSCs) are uniquely suited to study human development and disease and promise to revolutionize regenerative medicine. These applications rely on robust methods to manipulate gene function in hPSC models. This comprehensive review aims to both empower scientists approaching the field and update experienced stem cell biologists. We begin by highlighting challenges with manipulating gene expression in hPSCs and their differentiated derivatives, and relevant solutions (transfection, transduction, transposition, and genomic safe harbor editing). We then outline how to perform robust constitutive or inducible loss-, gain-, and change-of-function experiments in hPSCs models, both using historical methods (RNA interference, transgenesis, and homologous recombination) and modern programmable nucleases (particularly CRISPR/Cas9 and its derivatives, i.e., CRISPR interference, activation, base editing, and prime editing). We further describe extension of these approaches for arrayed or pooled functional studies, including emerging single-cell genomic methods, and the related design and analytical bioinformatic tools. Finally, we suggest some directions for future advancements in all of these areas. Mastering the combination of these transformative technologies will empower unprecedented advances in human biology and medicine.

Dynamic network-guided CRISPRi screen identifies CTCF-loop-constrained nonlinear enhancer gene regulatory activity during cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those contributing to complex diseases, have weak effects on gene expression. Our gene regulatory network modeling identified that nonlinear enhancer gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Using human embryonic stem cell definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. We discovered a comprehensive set of enhancers for each of the core endoderm-specifying transcription factors. Many enhancers had strong effects mid-transition but weak effects post-transition, consistent with the nonlinear temporal responses to enhancer perturbation predicted by the modeling. Integrating three-dimensional genomic information, we were able to develop a CTCF-loop-constrained Interaction Activity model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and systematic enhancer discovery in both normal and pathological cell state transitions.

Breaking enhancers to gain insights into developmental defects

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.