Sources of Error in Mammalian Genetic Screens

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.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Engineered CRISPR-Cas12a for higher-order combinatorial chromatin perturbations

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 one to three genomic sites per cell. We engineer an Acidaminococcus Cas12a (AsCas12a) variant, multiplexed transcriptional interference AsCas12a (multiAsCas12a), that incorporates R1226A, a mutation that stabilizes the ribonucleoprotein-DNA complex via DNA nicking. The multiAsCas12a-KRAB fusion improves CRISPRi activity over DNase-dead AsCas12a-KRAB fusions, often rescuing the activities of lentivirally delivered CRISPR RNAs (crRNA) that are inactive when used with the latter. multiAsCas12a-KRAB supports CRISPRi using 6-plex crRNA arrays in high-throughput pooled screens. Using multiAsCas12a-KRAB, we discover enhancer elements and dissect the combinatorial function of cis-regulatory elements in human cells. These results instantiate a group testing framework for efficiently surveying numerous combinations of chromatin perturbations for biological discovery and engineering.

Prediction of on-target and off-target activity of CRISPR-Cas13d guide RNAs using deep learning

Transcriptome engineering applications in living cells with RNA-targeting CRISPR effectors depend on accurate prediction of on-target activity and off-target avoidance. Here we design and test ~200,000 RfxCas13d guide RNAs targeting essential genes in human cells with systematically designed mismatches and insertions and deletions (indels). We find that mismatches and indels have a position- and context-dependent impact on Cas13d activity, and mismatches that result in G-U wobble pairings are better tolerated than other single-base mismatches. Using this large-scale dataset, we train a convolutional neural network that we term targeted inhibition of gene expression via gRNA design (TIGER) to predict efficacy from guide sequence and context. TIGER outperforms the existing models at predicting on-target and off-target activity on our dataset and published datasets. We show that TIGER scoring combined with specific mismatches yields the first general framework to modulate transcript expression, enabling the use of RNA-targeting CRISPRs to precisely control gene dosage.

CROPseq-multi: a versatile solution for multiplexed perturbation and decoding in pooled CRISPR screens

Forward genetic screens seek to dissect complex biological systems by systematically perturbing genetic elements and observing the resulting phenotypes. While standard screening methodologies introduce individual perturbations, multiplexing perturbations improves the performance of single-target screens and enables combinatorial screens for the study of genetic interactions. Current tools for multiplexing perturbations are incompatible with pooled screening methodologies that require mRNA-embedded barcodes, including some microscopy and single cell sequencing approaches. Here, we report the development of CROPseq-multi, a CROPseq1-inspired lentiviral system to multiplex Streptococcus pyogenes (Sp) Cas9-based perturbations with mRNA-embedded barcodes. CROPseq-multi has equivalent per-guide activity to CROPseq and low lentiviral recombination frequencies. CROPseq-multi is compatible with enrichment screening methodologies and optical pooled screens, and is extensible to screens with single-cell sequencing readouts. For optical pooled screens, an optimized and multiplexed in situ detection protocol improves barcode detection efficiency 10-fold, enables detection of recombination events, and increases decoding efficiency 3-fold relative to CROPseq. CROPseq-multi is a widely applicable multiplexing solution for diverse SpCas9-based genetic screening approaches.

CRISPR activation screens: navigating technologies and applications

Clustered regularly interspaced short palindromic repeats (CRISPR) activation (CRISPRa) has become an integral part of the molecular biology toolkit. CRISPRa genetic screens are an exciting high-throughput means of identifying genes the upregulation of which is sufficient to elicit a given phenotype. Activation machinery is continually under development to achieve greater, more robust, and more consistent activation. In this review, we offer a succinct technological overview of available CRISPRa architectures and a comprehensive summary of pooled CRISPRa screens. Furthermore, we discuss contemporary applications of CRISPRa across broad fields of research, with the aim of presenting a view of exciting emerging applications for CRISPRa screening.

Modular vector assembly enables rapid assessment of emerging CRISPR technologies

The diversity of CRISPR systems, coupled with scientific ingenuity, has led to an explosion of applications; however, to test newly described innovations in their model systems, researchers typically embark on cumbersome, one-off cloning projects to generate custom reagents that are optimized for their biological questions. Here, we leverage Golden Gate cloning to create the Fragmid toolkit, a modular set of CRISPR cassettes and delivery technologies, along with a web portal, resulting in a combinatorial platform that enables scalable vector assembly within days. We further demonstrate that multiple CRISPR technologies can be assessed in parallel in a pooled screening format using this resource, enabling the rapid optimization of both novel technologies and cellular models. These results establish Fragmid as a robust system for the rapid design of CRISPR vectors, and we anticipate that this assembly approach will be broadly useful for systematic development, comparison, and dissemination of CRISPR technologies.

Deep learning models to predict the editing efficiencies and outcomes of diverse base editors

Applications of base editing are frequently restricted by the requirement for a protospacer adjacent motif (PAM), and selecting the optimal base editor (BE) and single-guide RNA pair (sgRNA) for a given target can be difficult. To select for BEs and sgRNAs without extensive experimental work, we systematically compared the editing windows, outcomes and preferred motifs for seven BEs, including two cytosine BEs, two adenine BEs and three C•G to G•C BEs at thousands of target sequences. We also evaluated nine Cas9 variants that recognize different PAM sequences and developed a deep learning model, DeepCas9variants, for predicting which variants function most efficiently at sites with a given target sequence. We then develop a computational model, DeepBE, that predicts editing efficiencies and outcomes of 63 BEs that were generated by incorporating nine Cas9 variants as nickase domains into the seven BE variants. The predicted median efficiencies of BEs with DeepBE-based design were 2.9- to 20-fold higher than those of rationally designed SpCas9-containing BEs.

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.

Refining chimeric antigen receptors via barcoded protein domain combination pooled screening

Chimeric antigen receptor (CAR)-T cells represent a promising frontier in cancer immunotherapy. However, the current process for developing new CAR constructs is time consuming and inefficient. To address this challenge and expedite the evaluation and comparison of full-length CAR designs, we have devised a novel cloning strategy. This strategy involves the sequential assembly of individual CAR domains using blunt ligation, with each domain being assigned a unique DNA barcode. Applying this method, we successfully generated 360 CAR constructs that specifically target clinically validated tumor antigens CD19 and GD2. By quantifying changes in barcode frequencies through next-generation sequencing, we characterize CARs that best mediate proliferation and expansion of transduced T cells. The screening revealed a crucial role for the hinge domain in CAR functionality, with CD8a and IgG4 hinges having opposite effects in the surface expression, cytokine production, and antitumor activity in CD19- versus GD2-based CARs. Importantly, we discovered two novel CD19-CAR architectures containing the IgG4 hinge domain that mediate superior in vivo antitumor activity compared with the construct used in Kymriah, a U.S. Food and Drug Administration (FDA)-approved therapy. This novel screening approach represents a major advance in CAR engineering, enabling accelerated development of cell-based cancer immunotherapies.

Modular vector assembly enables rapid assessment of emerging CRISPR technologies

The diversity of CRISPR systems, coupled with scientific ingenuity, has led to an explosion of applications; however, to test newly-described innovations in their model systems, researchers typically embark on cumbersome, one-off cloning projects to generate custom reagents that are optimized for their biological questions. Here, we leverage Golden Gate cloning to create the Fragmid toolkit, a modular set of CRISPR cassettes and delivery technologies, along with a web portal, resulting in a combinatorial platform that enables scalable vector assembly within days. We further demonstrate that multiple CRISPR technologies can be assessed in parallel in a pooled screening format using this resource, enabling the rapid optimization of both novel technologies and cellular models. These results establish Fragmid as a robust system for the rapid design of CRISPR vectors, and we anticipate that this assembly approach will be broadly useful for systematic development, comparison, and dissemination of CRISPR technologies.

Modular pooled discovery of synthetic knockin sequences to program durable cell therapies

Chronic stimulation can cause T cell dysfunction and limit the efficacy of cellular immunotherapies. Improved methods are required to compare large numbers of synthetic knockin (KI) sequences to reprogram cell functions. Here, we developed modular pooled KI screening (ModPoKI), an adaptable platform for modular construction of DNA KI libraries using barcoded multicistronic adaptors. We built two ModPoKI libraries of 100 transcription factors (TFs) and 129 natural and synthetic surface receptors (SRs). Over 30 ModPoKI screens across human TCR- and CAR-T cells in diverse conditions identified a transcription factor AP4 (TFAP4) construct that enhanced fitness of chronically stimulated CAR-T cells and anti-cancer function in vitro and in vivo. ModPoKI's modularity allowed us to generate an ∼10,000-member library of TF combinations. Non-viral KI of a combined BATF-TFAP4 polycistronic construct enhanced fitness. Overexpressed BATF and TFAP4 co-occupy and regulate key gene targets to reprogram T cell function. ModPoKI facilitates the discovery of complex gene constructs to program cellular functions.

Optimization of Cas12a for multiplexed genome-scale transcriptional activation

Cas12a CRISPR technology, unlike Cas9, allows for facile multiplexing of guide RNAs from a single transcript, simplifying combinatorial perturbations. While Cas12a has been implemented for multiplexed knockout genetic screens, it has yet to be optimized for CRISPR activation (CRISPRa) screens in human cells. Here, we develop a new Cas12a-based transactivation domain (TAD) recruitment system using the ALFA nanobody and demonstrate simultaneous activation of up to four genes. We screen a genome-wide library to identify modulators of growth and MEK inhibition, and we compare these results with those obtained with open reading frame (ORF) overexpression and Cas9-based CRISPRa. We find that the activity of multiplexed arrays is largely predictable from the best-performing guide and provide criteria for selecting active guides. We anticipate that these results will greatly accelerate the exploration of gene function and combinatorial phenotypes at scale.

Massively parallel evaluation and computational prediction of the activities and specificities of 17 small Cas9s

Recently, various small Cas9 orthologs and variants have been reported for use in in vivo delivery applications. Although small Cas9s are particularly suited for this purpose, selecting the most optimal small Cas9 for use at a specific target sequence continues to be challenging. Here, to this end, we have systematically compared the activities of 17 small Cas9s for thousands of target sequences. For each small Cas9, we have characterized the protospacer adjacent motif and determined optimal single guide RNA expression formats and scaffold sequence. High-throughput comparative analyses revealed distinct high- and low-activity groups of small Cas9s. We also developed DeepSmallCas9, a set of computational models predicting the activities of the small Cas9s at matched and mismatched target sequences. Together, this analysis and these computational models provide a useful guide for researchers to select the most suitable small Cas9 for specific applications.

Salvage of ribose from uridine or RNA supports glycolysis in nutrient-limited conditions

Glucose is vital for life, serving as both a source of energy and carbon building block for growth. When glucose is limiting, alternative nutrients must be harnessed. To identify mechanisms by which cells can tolerate complete loss of glucose, we performed nutrient-sensitized genome-wide genetic screens and a PRISM growth assay across 482 cancer cell lines. We report that catabolism of uridine from the medium enables the growth of cells in the complete absence of glucose. While previous studies have shown that uridine can be salvaged to support pyrimidine synthesis in the setting of mitochondrial oxidative phosphorylation deficiency1, our work demonstrates that the ribose moiety of uridine or RNA can be salvaged to fulfil energy requirements via a pathway based on: (1) the phosphorylytic cleavage of uridine by uridine phosphorylase UPP1/UPP2 into uracil and ribose-1-phosphate (R1P), (2) the conversion of uridine-derived R1P into fructose-6-P and glyceraldehyde-3-P by the non-oxidative branch of the pentose phosphate pathway and (3) their glycolytic utilization to fuel ATP production, biosynthesis and gluconeogenesis. Capacity for glycolysis from uridine-derived ribose appears widespread, and we confirm its activity in cancer lineages, primary macrophages and mice in vivo. An interesting property of this pathway is that R1P enters downstream of the initial, highly regulated steps of glucose transport and upper glycolysis. We anticipate that 'uridine bypass' of upper glycolysis could be important in the context of disease and even exploited for therapeutic purposes.

Massively Parallel Reporter Assays for High-Throughput In Vivo Analysis of Cis-Regulatory Elements

The rapid improvement of descriptive genomic technologies has fueled a dramatic increase in hypothesized connections between cardiovascular gene expression and phenotypes. However, in vivo testing of these hypotheses has predominantly been relegated to slow, expensive, and linear generation of genetically modified mice. In the study of genomic cis-regulatory elements, generation of mice featuring transgenic reporters or cis-regulatory element knockout remains the standard approach. While the data obtained is of high quality, the approach is insufficient to keep pace with candidate identification and therefore results in biases introduced during the selection of candidates for validation. However, recent advances across a range of disciplines are converging to enable functional genomic assays that can be conducted in a high-throughput manner. Here, we review one such method, massively parallel reporter assays (MPRAs), in which the activities of thousands of candidate genomic regulatory elements are simultaneously assessed via the next-generation sequencing of a barcoded reporter transcript. We discuss best practices for MPRA design and use, with a focus on practical considerations, and review how this emerging technology has been successfully deployed in vivo. Finally, we discuss how MPRAs are likely to evolve and be used in future cardiovascular research.

Sniper2L is a high-fidelity Cas9 variant with high activity

Although several high-fidelity SpCas9 variants have been reported, it has been observed that this increased specificity is associated with reduced on-target activity, limiting the applications of the high-fidelity variants when efficient genome editing is required. Here, we developed an improved version of Sniper-Cas9, Sniper2L, which represents an exception to this trade-off trend as it showed higher specificity with retained high activity. We evaluated Sniper2L activities at a large number of target sequences and developed DeepSniper, a deep learning model that can predict the activity of Sniper2L. We also confirmed that Sniper2L can induce highly efficient and specific editing at a large number of target sequences when it is delivered as a ribonucleoprotein complex. Mechanically, the high specificity of Sniper2L originates from its superior ability to avoid unwinding a target DNA containing even a single mismatch. We envision that Sniper2L will be useful when efficient and specific genome editing is required.

Massively Parallel CRISPR-Based Genetic Perturbation Screening at Single-Cell Resolution

The clustered regularly interspaced short palindromic repeats (CRISPR)-based genetic screening has been demonstrated as a powerful approach for unbiased functional genomics research. Single-cell CRISPR screening (scCRISPR) techniques, which result from the combination of single-cell toolkits and CRISPR screening, allow dissecting regulatory networks in complex biological systems at unprecedented resolution. These methods allow cells to be perturbed en masse using a pooled CRISPR library, followed by high-content phenotyping. This is technically accomplished by annotating cells with sgRNA-specific barcodes or directly detectable sgRNAs. According to the integration of distinct single-cell technologies, these methods principally fall into four categories: scCRISPR with RNA-seq, scCRISPR with ATAC-seq, scCRISPR with proteome probing, and imaging-based scCRISPR. scCRISPR has deciphered genotype-phenotype relationships, genetic regulations, tumor biological issues, and neuropathological mechanisms. This review provides insight into the technical breakthrough of scCRISPR by systematically summarizing the advancements of various scCRISPR methodologies and analyzing their merits and limitations. In addition, an application-oriented approach guide is offered to meet researchers' individualized demands.

CRISPR in cancer biology and therapy

Over the past decade, CRISPR has become as much a verb as it is an acronym, transforming biomedical research and providing entirely new approaches for dissecting all facets of cell biology. In cancer research, CRISPR and related tools have offered a window into previously intractable problems in our understanding of cancer genetics, the noncoding genome and tumour heterogeneity, and provided new insights into therapeutic vulnerabilities. Here, we review the progress made in the development of CRISPR systems as a tool to study cancer, and the emerging adaptation of these technologies to improve diagnosis and treatment.

A multiplexed epitope barcoding strategy that enables dynamic cellular phenotypic screens

Pooled genetic libraries have improved screening throughput for mapping genotypes to phenotypes. However, selectable phenotypes are limited, restricting screening to outcomes with a low spatiotemporal resolution. Here, we integrated live-cell imaging with pooled library-based screening. To enable intracellular multiplexing, we developed a method called EPICode that uses a combination of short epitopes, which can also appear in various subcellular locations. EPICode thus enables the use of live-cell microscopy to characterize a phenotype of interest over time, including after sequential stimulatory/inhibitory manipulations, and directly connects behavior to the cellular genotype. To test EPICode's capacity against an important milestone-engineering and optimizing dynamic, live-cell reporters-we developed a live-cell PKA kinase translocation reporter with improved sensitivity and specificity. The use of epitopes as fluorescent barcodes introduces a scalable strategy for high-throughput screening broadly applicable to protein engineering and drug discovery settings where image-based phenotyping is desired.

Gain-of-function genetic screening identifies the antiviral function of TMEM120A via STING activation

Zika virus (ZIKV) infection can be associated with neurological pathologies, such as microcephaly in newborns and Guillain-Barre syndrome in adults. Effective therapeutics are currently not available. As such, a comprehensive understanding of virus-host interactions may guide the development of medications for ZIKV. Here we report a human genome-wide overexpression screen to identify host factors that regulate ZIKV infection and find TMEM120A as a ZIKV restriction factor. TMEM120A overexpression significantly inhibits ZIKV replication, while TMEM120A knockdown increases ZIKV infection in cell lines. Moreover, Tmem120a knockout in mice facilitates ZIKV infection in primary mouse embryonic fibroblasts (MEF) cells. Mechanistically, the antiviral activity of TMEM120A is dependent on STING, as TMEM120A interacts with STING, promotes the translocation of STING from the endoplasmic reticulum (ER) to ER-Golgi intermediate compartment (ERGIC) and enhances the phosphorylation of downstream TBK1 and IRF3, resulting in the expression of multiple antiviral cytokines and interferon-stimulated genes. In summary, our gain-of-function screening identifies TMEM120A as a key activator of the antiviral signaling of STING.

Understanding the interplay between host and viral factors during infection is essential for the interactome of infection. Here the authors perform a gain-of-function screen to identify factors involved during Zika virus infection and identify TMEM120A as a key factor in the STING mediated immune responses.

Charting oncogenicity of genes and variants across lineages via multiplexed screens in teratomas

Deconstructing tissue-specific effects of genes and variants on proliferation is critical to understanding cellular transformation and systematically selecting cancer therapeutics. This requires scalable methods for multiplexed genetic screens tracking fitness across time, across lineages, and in a suitable niche, since physiological cues influence functional differences. Towards this, we present an approach, coupling single-cell cancer driver screens in teratomas with hit enrichment by serial teratoma reinjection, to simultaneously screen drivers across multiple lineages in vivo. Using this system, we analyzed population shifts and lineage-specific enrichment for 51 cancer associated genes and variants, profiling over 100,000 cells spanning over 20 lineages, across two rounds of serial reinjection. We confirmed that c-MYC alone or combined with myristoylated AKT1 potently drives proliferation in progenitor neural lineages, demonstrating signatures of malignancy. Additionally, mutant MEK1 S218D/S222D provides a proliferative advantage in mesenchymal lineages like fibroblasts. Our method provides a powerful platform for multi-lineage longitudinal study of oncogenesis.

Application of CHyMErA Cas9-Cas12a combinatorial genome-editing platform for genetic interaction mapping and gene fragment deletion screening

CRISPR-based forward genetic screening represents a powerful approach for the systematic characterization of gene function. Recent efforts have been directed toward establishing CRISPR-based tools for the programmable delivery of combinatorial genetic perturbations, most of which are mediated by a single nuclease and the expression of structurally identical guide backbones from two promoters. In contrast, we have developed CHyMErA (Cas hybrid for multiplexed editing and screening applications), which is based on the co-expression of Cas9 and Cas12a nucleases in conjunction with a hybrid guide RNA (hgRNA) engineered by the fusion of Cas9 and Cas12a guides and expressed from a single U6 promoter. CHyMErA is suitable for the high-throughput deletion of genetic segments including the excision of individual exons. Furthermore, CHyMErA enables the concomitant targeting of two or more genes and can thus be used for the systematic mapping of genetic interactions in mammalian cells. CHyMErA can also be applied for the perturbation of paralogous gene pairs, thereby allowing the capturing of phenotypic roles that would otherwise be masked because of genetic redundancy. Here, we provide instructions for the cloning of hgRNA screening libraries and individual hgRNA constructs and offer guidelines for designing and performing combinatorial pooled genetic screens using CHyMErA. Starting with the generation of Cas9- and Cas12a-expressing cell lines, CHyMErA screening can be implemented within 15-20 weeks.

Predicting the efficiency of prime editing guide RNAs in human cells

Prime editing enables the introduction of virtually any small-sized genetic change without requiring donor DNA or double-strand breaks. However, evaluation of prime editing efficiency requires time-consuming experiments, and the factors that affect efficiency have not been extensively investigated. In this study, we performed high-throughput evaluation of prime editor 2 (PE2) activities in human cells using 54,836 pairs of prime editing guide RNAs (pegRNAs) and their target sequences. The resulting data sets allowed us to identify factors affecting PE2 efficiency and to develop three computational models to predict pegRNA efficiency. For a given target sequence, the computational models predict efficiencies of pegRNAs with different lengths of primer binding sites and reverse transcriptase templates for edits of various types and positions. Testing the accuracy of the predictions using test data sets that were not used for training, we found Spearman's correlations between 0.47 and 0.81. Our computational models and information about factors affecting PE2 efficiency will facilitate practical application of prime editing.

Massively parallel functional testing of MSH2 missense variants conferring Lynch syndrome risk

The lack of functional evidence for the majority of missense variants limits their clinical interpretability and poses a key barrier to the broad utility of carrier screening. In Lynch syndrome (LS), one of the most highly prevalent cancer syndromes, nearly 90% of clinically observed missense variants are deemed “variants of uncertain significance” (VUS). To systematically resolve their functional status, we performed a massively parallel screen in human cells to identify loss-of-function missense variants in the key DNA mismatch repair factor MSH2. The resulting functional effect map is substantially complete, covering 94% of the 17,746 possible variants, and is highly concordant (96%) with existing functional data and expert clinicians’ interpretations. The large majority (89%) of missense variants were functionally neutral, perhaps unexpectedly in light of its evolutionary conservation. These data provide ready-to-use functional evidence to resolve the ∼1,300 extant missense VUSs in MSH2 and may facilitate the prospective classification of newly discovered variants in the clinic.

A systematic evaluation of the design and context dependencies of massively parallel reporter assays

Massively parallel reporter assays (MPRAs) functionally screen thousands of sequences for regulatory activity in parallel. To date, there are limited studies that systematically compare differences in MPRA design. Here, we screen a library of 2,440 candidate liver enhancers and controls for regulatory activity in HepG2 cells using nine different MPRA designs. We identify subtle but significant differences that correlate with epigenetic and sequence-level features, as well as differences in dynamic range and reproducibility. We also validate that enhancer activity is largely independent of orientation, at least for our library and designs. Finally, we assemble and test the same enhancers as 192-mers, 354-mers and 678-mers and observe sizable differences. This work provides a framework for the experimental design of high-throughput reporter assays, suggesting that the extended sequence context of tested elements and to a lesser degree the precise assay, influence MPRA results.

Prediction of the sequence-specific cleavage activity of Cas9 variants

Several Streptococcus pyogenes Cas9 (SpCas9) variants have been developed to improve an enzyme's specificity or to alter or broaden its protospacer-adjacent motif (PAM) compatibility, but selecting the optimal variant for a given target sequence and application remains difficult. To build computational models to predict the sequence-specific activity of 13 SpCas9 variants, we first assessed their cleavage efficiency at 26,891 target sequences. We found that, of the 256 possible four-nucleotide NNNN sequences, 156 can be used as a PAM by at least one of the SpCas9 variants. For the high-fidelity variants, overall activity could be ranked as SpCas9 ≥ Sniper-Cas9 > eSpCas9(1.1) > SpCas9-HF1 > HypaCas9 ≈ xCas9 >> evoCas9, whereas their overall specificities could be ranked as evoCas9 >> HypaCas9 ≥ SpCas9-HF1 ≈ eSpCas9(1.1) > xCas9 > Sniper-Cas9 > SpCas9. Using these data, we developed 16 deep-learning-based computational models that accurately predict the activity of these variants at any target sequence.

Rational Reprogramming of Cellular States by Combinatorial Perturbation

Ectopic expression of transcription factors (TFs) can reprogram cell state. However, because of the large combinatorial space of possible TF cocktails, it remains difficult to identify TFs that reprogram specific cell types. Here, we develop Reprogram-Seq to experimentally screen thousands of TF cocktails for reprogramming performance. Reprogram-Seq leverages organ-specific cell-atlas data with single-cell perturbation and computational analysis to predict, evaluate, and optimize TF combinations that reprogram a cell type of interest. Focusing on the cardiac system, we perform Reprogram-Seq on MEFs using an undirected library of 48 cardiac factors and, separately, a directed library of 10 epicardial-related TFs. We identify a combination of three TFs, which efficiently reprogram MEFs to epicardial-like cells that are transcriptionally, molecularly, morphologically, and functionally similar to primary epicardial cells. Reprogram-Seq holds promise to accelerate the generation of specific cell types for regenerative medicine.

Dissection of c-AMP Response Element Architecture by Using Genomic and Episomal Massively Parallel Reporter Assays

In eukaryotes, transcription factors (TFs) orchestrate gene expression by binding to TF-binding sites (TFBSs) and localizing transcriptional co-regulators and RNA polymerase II to cis-regulatory elements. However, we lack a basic understanding of the relationship between TFBS composition and their quantitative transcriptional responses. Here, we measured expression driven by 17,406 synthetic cis-regulatory elements with varied compositions of a model TFBS, the c-AMP response element (CRE) by using massively parallel reporter assays (MPRAs). We find CRE number, affinity, and promoter proximity largely determines expression. In addition, we observe expression modulation based on the spacing between CREs and CRE distance to the promoter, where expression follows a helical periodicity. Finally, we compare library expression between an episomal MPRA and a genomically integrated MPRA, where a single cis-regulatory element is assayed per cell at a defined locus. These assays largely recapitulate each other, although weaker, non-canonical CREs exhibit greater activity in a genomic context.

Pooled Knockin Targeting for Genome Engineering of Cellular Immunotherapies

Adoptive transfer of genetically modified immune cells holds great promise for cancer immunotherapy. CRISPR knockin targeting can improve cell therapies, but more high-throughput methods are needed to test which knockin gene constructs most potently enhance primary cell functions in vivo. We developed a widely adaptable technology to barcode and track targeted integrations of large non-viral DNA templates and applied it to perform pooled knockin screens in primary human T cells. Pooled knockin of dozens of unique barcoded templates into the T cell receptor (TCR)-locus revealed gene constructs that enhanced fitness in vitro and in vivo. We further developed pooled knockin sequencing (PoKI-seq), combining single-cell transcriptome analysis and pooled knockin screening to measure cell abundance and cell state ex vivo and in vivo. This platform nominated a novel transforming growth factor β (TGF-β) R2-41BB chimeric receptor that improved solid tumor clearance. Pooled knockin screening enables parallelized re-writing of endogenous genetic sequences to accelerate discovery of knockin programs for cell therapies.

Genetic interaction mapping and exon-resolution functional genomics with a hybrid Cas9-Cas12a platform

Systematic mapping of genetic interactions (GIs) and interrogation of the functions of sizable genomic segments in mammalian cells represent important goals of biomedical research. To advance these goals, we present a CRISPR (clustered regularly interspaced short palindromic repeats)-based screening system for combinatorial genetic manipulation that employs coexpression of CRISPR-associated nucleases 9 and 12a (Cas9 and Cas12a) and machine-learning-optimized libraries of hybrid Cas9-Cas12a guide RNAs. This system, named Cas Hybrid for Multiplexed Editing and screening Applications (CHyMErA), outperforms genetic screens using Cas9 or Cas12a editing alone. Application of CHyMErA to the ablation of mammalian paralog gene pairs reveals extensive GIs and uncovers phenotypes normally masked by functional redundancy. Application of CHyMErA in a chemogenetic interaction screen identifies genes that impact cell growth in response to mTOR pathway inhibition. Moreover, by systematically targeting thousands of alternative splicing events, CHyMErA identifies exons underlying human cell line fitness. CHyMErA thus represents an effective screening approach for GI mapping and the functional analysis of sizable genomic regions, such as alternative exons.

An improved platform for functional assessment of large protein libraries in mammalian cells

Multiplex genetic assays can simultaneously test thousands of genetic variants for a property of interest. However, limitations of existing multiplex assay methods in cultured mammalian cells hinder the breadth, speed and scale of these experiments. Here, we describe a series of improvements that greatly enhance the capabilities of a Bxb1 recombinase-based landing pad system for conducting different types of multiplex genetic assays in various mammalian cell lines. We incorporate the landing pad into a lentiviral vector, easing the process of generating new landing pad cell lines. We also develop several new landing pad versions, including one where the Bxb1 recombinase is expressed from the landing pad itself, improving recombination efficiency more than 2-fold and permitting rapid prototyping of transgenic constructs. Other versions incorporate positive and negative selection markers that enable drug-based enrichment of recombinant cells, enabling the use of larger libraries and reducing costs. A version with dual convergent promoters allows enrichment of recombinant cells independent of transgene expression, permitting the assessment of libraries of transgenes that perturb cell growth and survival. Lastly, we demonstrate these improvements by assessing the effects of a combinatorial library of oncogenes and tumor suppressors on cell growth. Collectively, these advancements make multiplex genetic assays in diverse cultured cell lines easier, cheaper and more effective, facilitating future studies probing how proteins impact cell function, using transgenic variant libraries tested individually or in combination.

Massively parallel profiling and predictive modeling of the outcomes of CRISPR/Cas9-mediated double-strand break repair

Non-homologous end-joining (NHEJ) plays an important role in double-strand break (DSB) repair of DNA. Recent studies have shown that the error patterns of NHEJ are strongly biased by sequence context, but these studies were based on relatively few templates. To investigate this more thoroughly, we systematically profiled ∼1.16 million independent mutational events resulting from CRISPR/Cas9-mediated cleavage and NHEJ-mediated DSB repair of 6872 synthetic target sequences, introduced into a human cell line via lentiviral infection. We find that: (i) insertions are dominated by 1 bp events templated by sequence immediately upstream of the cleavage site, (ii) deletions are predominantly associated with microhomology and (iii) targets exhibit variable but reproducible diversity with respect to the number and relative frequency of the mutational outcomes to which they give rise. From these data, we trained a model that uses local sequence context to predict the distribution of mutational outcomes. Exploiting the bias of NHEJ outcomes towards microhomology mediated events, we demonstrate the programming of deletion patterns by introducing microhomology to specific locations in the vicinity of the DSB site. We anticipate that our results will inform investigations of DSB repair mechanisms as well as the design of CRISPR/Cas9 experiments for diverse applications including genome-wide screens, gene therapy, lineage tracing and molecular recording.

An LKB1-SIK Axis Suppresses Lung Tumor Growth and Controls Differentiation

The kinase LKB1 is a critical tumor suppressor in sporadic and familial human cancers, yet the mechanisms by which it suppresses tumor growth remain poorly understood. To investigate the tumor-suppressive capacity of four canonical families of LKB1 substrates in vivo, we used CRISPR/Cas9-mediated combinatorial genome editing in a mouse model of oncogenic KRAS-driven lung adenocarcinoma. We demonstrate that members of the SIK family are critical for constraining tumor development. Histologic and gene-expression similarities between LKB1- and SIK-deficient tumors suggest that SIKs and LKB1 operate within the same axis. Furthermore, a gene-expression signature reflecting SIK deficiency is enriched in LKB1-mutant human lung adenocarcinomas and is regulated by LKB1 in human cancer cell lines. Together, these findings reveal a key LKB1-SIK tumor-suppressive axis and underscore the need to redirect efforts to elucidate the mechanisms through which LKB1 mediates tumor suppression. SIGNIFICANCE: Uncovering the effectors of frequently altered tumor suppressor genes is critical for understanding the fundamental driving forces of cancer growth. Our identification of the SIK family of kinases as effectors of LKB1-mediated tumor suppression will refocus future mechanistic studies and may lead to new avenues for genotype-specific therapeutic interventions.This article is highlighted in the In This Issue feature, p. 1469.

Optical Pooled Screens in Human Cells

Genetic screens are critical for the systematic identification of genes underlying cellular phenotypes. Pooling gene perturbations greatly improves scalability but is not compatible with imaging of complex and dynamic cellular phenotypes. Here, we introduce a pooled approach for optical genetic screens in mammalian cells. We use targeted in situ sequencing to demultiplex a library of genetic perturbations following image-based phenotyping. We screened a set of 952 genes across millions of cells for involvement in nuclear factor κB (NF-κB) signaling by imaging the translocation of RelA (p65) to the nucleus. Screening at a single time point across 3 cell lines recovered 15 known pathway components, while repeating the screen with live-cell imaging revealed a role for Mediator complex subunits in regulating the duration of p65 nuclear retention. These results establish a highly multiplexed approach to image-based screens of spatially and temporally defined phenotypes with pooled libraries.

Optimized CRISPR guide RNA design for two high-fidelity Cas9 variants by deep learning

Highly specific Cas9 nucleases derived from SpCas9 are valuable tools for genome editing, but their wide applications are hampered by a lack of knowledge governing guide RNA (gRNA) activity. Here, we perform a genome-scale screen to measure gRNA activity for two highly specific SpCas9 variants (eSpCas9(1.1) and SpCas9-HF1) and wild-type SpCas9 (WT-SpCas9) in human cells, and obtain indel rates of over 50,000 gRNAs for each nuclease, covering ~20,000 genes. We evaluate the contribution of 1,031 features to gRNA activity and develope models for activity prediction. Our data reveals that a combination of RNN with important biological features outperforms other models for activity prediction. We further demonstrate that our model outperforms other popular gRNA design tools. Finally, we develop an online design tool DeepHF for the three Cas9 nucleases. The database, as well as the designer tool, is freely accessible via a web server, http://www.DeepHF.com/ .

Interspecies analysis of MYC targets identifies tRNA synthetases as mediators of growth and survival in MYC-overexpressing cells

Aberrant MYC oncogene activation is one of the most prevalent characteristics of cancer. By overlapping datasets of Drosophila genes that are insulin-responsive and also regulate nucleolus size, we enriched for Myc target genes required for cellular biosynthesis. Among these, we identified the aminoacyl tRNA synthetases (aaRSs) as essential mediators of Myc growth control in Drosophila and found that their pharmacologic inhibition is sufficient to kill MYC-overexpressing human cells, indicating that aaRS inhibitors might be used to selectively target MYC-driven cancers. We suggest a general principle in which oncogenic increases in cellular biosynthesis sensitize cells to disruption of protein homeostasis.

Imaging-based pooled CRISPR screening reveals regulators of lncRNA localization

Pooled-library CRISPR screening provides a powerful means to discover genetic factors involved in cellular processes in a high-throughput manner. However, the phenotypes accessible to pooled-library screening are limited. Complex phenotypes, such as cellular morphology and subcellular molecular organization, as well as their dynamics, require imaging-based readout and are currently beyond the reach of pooled-library CRISPR screening. Here we report an all imaging-based pooled-library CRISPR screening approach that combines high-content phenotype imaging with high-throughput single guide RNA (sgRNA) identification in individual cells. In this approach, sgRNAs are codelivered to cells with corresponding barcodes placed at the 3' untranslated region of a reporter gene using a lentiviral delivery system with reduced recombination-induced sgRNA-barcode mispairing. Multiplexed error-robust fluorescence in situ hybridization (MERFISH) is used to read out the barcodes and hence identify the sgRNAs with high accuracy. We used this approach to screen 162 sgRNAs targeting 54 RNA-binding proteins for their effects on RNA localization to nuclear compartments and uncovered previously unknown regulatory factors for nuclear RNA localization. Notably, our screen revealed both positive and negative regulators for the nuclear speckle localization of a long noncoding RNA, MALAT1, suggesting a dynamic regulation of lncRNA localization in subcellular compartments.

Mapping Cellular Reprogramming via Pooled Overexpression Screens with Paired Fitness and Single-Cell RNA-Sequencing Readout

Understanding the effects of genetic perturbations on the cellular state has been challenging using traditional pooled screens, which typically rely on the delivery of a single perturbation per cell and unidimensional phenotypic readouts. Here, we use barcoded open reading frame overexpression libraries coupled with single-cell RNA sequencing to assay cell state and fitness, a technique we call SEUSS (scalable functional screening by sequencing). Using SEUSS, we perturbed hPSCs with a library of developmentally critical transcription factors (TFs) and assayed the impact of TF overexpression on fitness and transcriptomic states. We further leveraged the versatility of the ORF library approach to assay mutant genes and whole gene families. From the transcriptomic responses, we built genetic co-regulatory networks to identify altered gene modules and found that KLF4 and SNAI2 drive opposing effects along the epithelial-mesenchymal transition axis. From the fitness responses, we identified ETV2 as a driver of reprogramming toward an endothelial-like state.

Protein Barcodes Enable High-Dimensional Single-Cell CRISPR Screens

CRISPR pools are being widely employed to identify gene functions. However, current technology, which utilizes DNA as barcodes, permits limited phenotyping and bulk-cell resolution. To enable novel screening capabilities, we developed a barcoding system operating at the protein level. We synthesized modules encoding triplet combinations of linear epitopes to generate >100 unique protein barcodes (Pro-Codes). Pro-Code-expressing vectors were introduced into cells and analyzed by CyTOF mass cytometry. Using just 14 antibodies, we detected 364 Pro-Code populations; establishing the largest set of protein-based reporters. By pairing each Pro-Code with a different CRISPR, we simultaneously analyzed multiple phenotypic markers, including phospho-signaling, on dozens of knockouts. Pro-Code/CRISPR screens found two interferon-stimulated genes, the immunoproteasome component Psmb8 and a chaperone Rtp4, are important for antigen-dependent immune editing of cancer cells and identified Socs1 as a negative regulator of Pd-l1. The Pro-Code technology enables simultaneous high-dimensional protein-level phenotyping of 100s of genes with single-cell resolution.

A Multiplex Homology-Directed DNA Repair Assay Reveals the Impact of More Than 1,000 BRCA1 Missense Substitution Variants on Protein Function

Loss-of-function pathogenic variants in BRCA1 confer a predisposition to breast and ovarian cancer. Genetic testing for sequence changes in BRCA1 frequently reveals a missense variant for which the impact on cancer risk and on the molecular function of BRCA1 is unknown. Functional BRCA1 is required for the homology-directed repair (HDR) of double-strand DNA breaks, a critical activity for maintaining genome integrity and tumor suppression. Here, we describe a multiplex HDR reporter assay for concurrently measuring the effects of hundreds of variants of BRCA1 for their role in DNA repair. Using this assay, we characterized the effects of 1,056 amino acid substitutions in the first 192 residues of BRCA1. Benchmarking these results against variants with known effects on DNA repair function or on cancer predisposition, we demonstrate accurate discrimination of loss-of-function versus benign missense variants. We anticipate that this assay can be used to functionally characterize BRCA1 missense variants at scale, even before the variants are observed in results from genetic testing.

Mapping the Genetic Landscape of Human Cells

Seminal yeast studies have established the value of comprehensively mapping genetic interactions (GIs) for inferring gene function. Efforts in human cells using focused gene sets underscore the utility of this approach, but the feasibility of generating large-scale, diverse human GI maps remains unresolved. We developed a CRISPR interference platform for large-scale quantitative mapping of human GIs. We systematically perturbed 222,784 gene pairs in two cancer cell lines. The resultant maps cluster functionally related genes, assigning function to poorly characterized genes, including TMEM261, a new electron transport chain component. Individual GIs pinpoint unexpected relationships between pathways, exemplified by a specific cholesterol biosynthesis intermediate whose accumulation induces deoxynucleotide depletion, causing replicative DNA damage and a synthetic-lethal interaction with the ATR/9-1-1 DNA repair pathway. Our map provides a broad resource, establishes GI maps as a high-resolution tool for dissecting gene function, and serves as a blueprint for mapping the genetic landscape of human cells.

Frequent sgRNA-barcode recombination in single-cell perturbation assays

Simultaneously detecting CRISPR-based perturbations and induced transcriptional changes in the same cell is a powerful approach to unraveling genome function. Several lentiviral approaches have been developed, some of which rely on the detection of distally located genetic barcodes as an indirect proxy of sgRNA identity. Since barcodes are often several kilobases from their corresponding sgRNAs, viral recombination-mediated swapping of barcodes and sgRNAs is feasible. Using a self-circularization-based sgRNA-barcode library preparation protocol, we estimate the recombination rate to be ~50% and we trace this phenomenon to the pooled viral packaging step. Recombination is random, and decreases the signal-to-noise ratio of the assay. Our results suggest that alternative approaches can increase the throughput and sensitivity of single-cell perturbation assays.

Uncoupling of sgRNAs from their associated barcodes during PCR amplification of combinatorial CRISPR screens

Many implementations of pooled screens in mammalian cells rely on linking an element of interest to a barcode, with the latter subsequently quantitated by next generation sequencing. However, substantial uncoupling between these paired elements during lentiviral production has been reported, especially as the distance between elements increases. We detail that PCR amplification is another major source of uncoupling, and becomes more pronounced with increased amounts of DNA template molecules and PCR cycles. To lessen uncoupling in systems that use paired elements for detection, we recommend minimizing the distance between elements, using low and equal template DNA inputs for plasmid and genomic DNA during PCR, and minimizing the number of PCR cycles. We also present a vector design for conducting combinatorial CRISPR screens that enables accurate barcode-based detection with a single short sequencing read and minimal uncoupling.

Profound Tissue Specificity in Proliferation Control Underlies Cancer Drivers and Aneuploidy Patterns

Genomics has provided a detailed structural description of the cancer genome. Identifying oncogenic drivers that work primarily through dosage changes is a current challenge. Unrestrained proliferation is a critical hallmark of cancer. We constructed modular, barcoded libraries of human open reading frames (ORFs) and performed screens for proliferation regulators in multiple cell types. Approximately 10% of genes regulate proliferation, with most performing in an unexpectedly highly tissue-specific manner. Proliferation drivers in a given cell type showed specific enrichment in somatic copy number changes (SCNAs) from cognate tumors and helped predict aneuploidy patterns in those tumors, implying that tissue-type-specific genetic network architectures underlie SCNA and driver selection in different cancers. In vivo screening confirmed these results. We report a substantial contribution to the catalog of SCNA-associated cancer drivers, identifying 147 amplified and 107 deleted genes as potential drivers, and derive insights about the genetic network architecture of aneuploidy in tumors.

On the design of CRISPR-based single-cell molecular screens

Several groups recently coupled CRISPR perturbations and single-cell RNA-seq for pooled genetic screens. We demonstrate that vector designs of these studies are susceptible to ∼50% swapping of guide RNA-barcode associations because of lentiviral template switching. We optimized a published alternative, CROP-seq, in which the guide RNA also serves as the barcode, and here confirm that this strategy performs robustly and doubled the rate at which guides are assigned to cells to 94%.

Synergistic drug combinations for cancer identified in a CRISPR screen for pairwise genetic interactions

Identification of effective combination therapies is critical to address the emergence of drug-resistant cancers, but direct screening of all possible drug combinations is infeasible. Here we introduce a CRISPR-based double knockout (CDKO) system that improves the efficiency of combinatorial genetic screening using an effective strategy for cloning and sequencing paired single guide RNA (sgRNA) libraries and a robust statistical scoring method for calculating genetic interactions (GIs) from CRISPR-deleted gene pairs. We applied CDKO to generate a large-scale human GI map, comprising 490,000 double-sgRNAs directed against 21,321 pairs of drug targets in K562 leukemia cells and identified synthetic lethal drug target pairs for which corresponding drugs exhibit synergistic killing. These included the BCL2L1 and MCL1 combination, which was also effective in imatinib-resistant cells. We further validated this system by identifying known and previously unidentified GIs between modifiers of ricin toxicity. This work provides an effective strategy to screen synergistic drug combinations in high-throughput and a CRISPR-based tool to dissect functional GI networks.

A Multiplexed Single-Cell CRISPR Screening Platform Enables Systematic Dissection of the Unfolded Protein Response

Functional genomics efforts face tradeoffs between number of perturbations examined and complexity of phenotypes measured. We bridge this gap with Perturb-seq, which combines droplet-based single-cell RNA-seq with a strategy for barcoding CRISPR-mediated perturbations, allowing many perturbations to be profiled in pooled format. We applied Perturb-seq to dissect the mammalian unfolded protein response (UPR) using single and combinatorial CRISPR perturbations. Two genome-scale CRISPR interference (CRISPRi) screens identified genes whose repression perturbs ER homeostasis. Subjecting ∼100 hits to Perturb-seq enabled high-precision functional clustering of genes. Single-cell analyses decoupled the three UPR branches, revealed bifurcated UPR branch activation among cells subject to the same perturbation, and uncovered differential activation of the branches across hits, including an isolated feedback loop between the translocon and IRE1α. These studies provide insight into how the three sensors of ER homeostasis monitor distinct types of stress and highlight the ability of Perturb-seq to dissect complex cellular responses.