Computational identification of clonal cells in single-cell CRISPR screens

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.

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.

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Multiplex generation and single cell analysis of structural variants in a mammalian genome

The functional consequences of structural variants (SVs) in mammalian genomes are challenging to study. This is due to several factors, including: 1) their numerical paucity relative to other forms of standing genetic variation such as single nucleotide variants (SNVs) and short insertions or deletions (indels); 2) the fact that a single SV can involve and potentially impact the function of more than one gene and/or cis regulatory element; and 3) the relative immaturity of methods to generate and map SVs, either randomly or in targeted fashion, in in vitro or in vivo model systems. Towards addressing these challenges, we developed Genome-Shuffle-seq, a straightforward method that enables the multiplex generation and mapping of several major forms of SVs (deletions, inversions, translocations) throughout a mammalian genome. Genome-Shuffle-seq is based on the integration of "shuffle cassettes" to the genome, wherein each shuffle cassette contains components that facilitate its site-specific recombination (SSR) with other integrated shuffle cassettes (via Cre-loxP), its mapping to a specific genomic location (via T7-mediated in vitro transcription or IVT), and its identification in single-cell RNA-seq (scRNA-seq) data (via T7-mediated in situ transcription or IST). In this proof-of-concept, we apply Genome-Shuffle-seq to induce and map thousands of genomic SVs in mouse embryonic stem cells (mESCs) in a single experiment. Induced SVs are rapidly depleted from the cellular population over time, possibly due to Cre-mediated toxicity and/or negative selection on the rearrangements themselves. Leveraging T7 IST of barcodes whose positions are already mapped, we further demonstrate that we can efficiently genotype which SVs are present in association with each of many single cell transcriptomes in scRNA-seq data. Finally, preliminary evidence suggests our method may be a powerful means of generating extrachromosomal circular DNAs (ecDNAs). Looking forward, we anticipate that Genome-Shuffle-seq may be broadly useful for the systematic exploration of the functional consequences of SVs on gene expression, the chromatin landscape, and 3D nuclear architecture. We further anticipate potential uses for in vitro modeling of ecDNAs, as well as in paving the path to a minimal mammalian genome.

Technological Convergence: Highlighting the Power of CRISPR Single-Cell Perturbation Toolkit for Functional Interrogation of Enhancers

Higher eukaryotic enhancers, as a major class of regulatory elements, play a crucial role in the regulation of gene expression. Over the last decade, the development of sequencing technologies has flooded researchers with transcriptome-phenotype data alongside emerging candidate regulatory elements. Since most methods can only provide hints about enhancer function, there have been attempts to develop experimental and computational approaches that can bridge the gap in the causal relationship between regulatory regions and phenotypes. The coupling of two state-of-the-art technologies, also referred to as crisprQTL, has emerged as a promising high-throughput toolkit for addressing this question. This review provides an overview of the importance of studying enhancers, the core molecular foundation of crisprQTL, and recent studies utilizing crisprQTL to interrogate enhancer-phenotype correlations. Additionally, we discuss computational methods currently employed for crisprQTL data analysis. We conclude by pointing out common challenges, making recommendations, and looking at future prospects, with the aim of providing researchers with an overview of crisprQTL as an important toolkit for studying enhancers.

Context-Dependent Function of Long Noncoding RNA PURPL in Transcriptome Regulation during p53 Activation

PURPL is a p53-induced lncRNA that suppresses basal p53 levels. Here, we investigated PURPL upon p53 activation in liver cancer cells, where it is expressed at significantly higher levels than other cell types. Using isoform sequencing, we discovered novel PURPL transcripts that have a retained intron and/or previously unannotated exons. To determine PURPL function upon p53 activation, we performed transcriptome sequencing (RNA-Seq) after depleting PURPL using CRISPR interference (CRISPRi), followed by Nutlin treatment to induce p53. Strikingly, although loss of PURPL in untreated cells altered the expression of only 7 genes, loss of PURPL resulted in altered expression of ~800 genes upon p53 activation, revealing a context-dependent function of PURPL. Pathway analysis suggested that PURPL is important for fine-tuning the expression of specific genes required for mitosis. Consistent with these results, we observed a significant decrease in the percentage of mitotic cells upon PURPL depletion. Collectively, these data identify novel transcripts from the PURPL locus and suggest that PURPL delicately moderates the expression of mitotic genes in the context of p53 activation to control cell cycle arrest.