A benchmarked, high-efficiency prime editing platform for multiplexed dropout screening.

Prime editing installs precise edits into the genome with minimal unwanted byproducts, but low and variable editing efficiencies have complicated application of the approach to high-throughput functional genomics. Leveraging several recent advances, we assembled a prime editing platform capable of high-efficiency substitution editing across a set of engineered prime editing guide RNAs (epegRNAs) and corresponding target sequences (80% median intended editing). Then, using a custom library of 240,000 epegRNAs targeting >17,000 codons with 175 different substitution types, we benchmarked our platform for functional interrogation of small substitution variants (1-3 nucleotides) targeted to essential genes. Resulting data identified negative growth phenotypes for nonsense mutations targeted to ~8,000 codons, and comparing those phenotypes to results from controls demonstrated high specificity. We also observed phenotypes for synonymous mutations that disrupted splice site motifs at 3' exon boundaries. Altogether, we establish and benchmark a high-throughput prime editing approach for functional characterization of genetic variants with simple readouts from multiplexed experiments.

Improving prime editing with an endogenous small RNA-binding protein.

Prime editing enables the precise modification of genomes through reverse transcription of template sequences appended to the 3' ends of CRISPR-Cas guide RNAs1. To identify cellular determinants of prime editing, we developed scalable prime editing reporters and performed genome-scale CRISPR-interference screens. From these screens, a single factor emerged as the strongest mediator of prime editing: the small RNA-binding exonuclease protection factor La. Further investigation revealed that La promotes prime editing across approaches (PE2, PE3, PE4 and PE5), edit types (substitutions, insertions and deletions), endogenous loci and cell types but has no consistent effect on genome-editing approaches that rely on standard, unextended guide RNAs. Previous work has shown that La binds polyuridine tracts at the 3' ends of RNA polymerase III transcripts2. We found that La functionally interacts with the 3' ends of polyuridylated prime editing guide RNAs (pegRNAs). Guided by these results, we developed a prime editor protein (PE7) fused to the RNA-binding, N-terminal domain of La. This editor improved prime editing with expressed pegRNAs and engineered pegRNAs (epegRNAs), as well as with synthetic pegRNAs optimized for La binding. Together, our results provide key insights into how prime editing components interact with the cellular environment and suggest general strategies for stabilizing exogenous small RNAs therein.

Efficient prime editing in two-cell mouse embryos using PEmbryo.

Using transient inhibition of DNA mismatch repair during a permissive stage of development, we demonstrate highly efficient prime editing of mouse embryos with few unwanted, local byproducts (average 58% precise edit frequency, 0.5% on-target error frequency across 13 substitution edits at 8 sites), enabling same-generation phenotyping of founders. Whole-genome sequencing reveals that mismatch repair inhibition increases off-target indels at low-complexity regions in the genome without any obvious phenotype in mice.

Membrane-proximal motifs encode differences in signaling strength between type I and III interferon receptors.

Interferons (IFNs) play crucial roles in antiviral defenses. Despite using the same Janus-activated kinase (JAK)-signal transducer and activator of transcription (STAT) signaling cascade, type I and III IFN receptors differ in the magnitude and dynamics of their signaling in terms of STAT phosphorylation, gene transcription, and antiviral responses. These differences are not due to ligand-binding affinity and receptor abundance. Here, we investigated the ability of the intracellular domains (ICDs) of IFN receptors to differentiate between type I and III IFN signaling. We engineered synthetic, heterodimeric type I and III IFN receptors that were stably expressed at similar amounts in human cells and responded to a common ligand. We found that our synthetic type I IFN receptors stimulated STAT phosphorylation and gene expression to greater extents than did the corresponding type III IFN receptors. Furthermore, we identified short "box motifs" within ICDs that bind to JAK1 that were sufficient to encode differences between the type I and III IFN receptors. Together, our results indicate that specific regions within the ICDs of IFN receptor subunits encode different downstream signaling strengths that enable type I and III IFN receptors to produce distinct signaling outcomes.

Mapping the Genetic Interaction Network of PARP inhibitor Response.

Genetic interactions have long informed our understanding of the coordinated proteins and pathways that respond to DNA damage in mammalian cells, but systematic interrogation of the genetic network underlying that system has yet to be achieved. Towards this goal, we measured 147,153 pairwise interactions among genes implicated in PARP inhibitor (PARPi) response. Evaluating genetic interactions at this scale, with and without exposure to PARPi, revealed hierarchical organization of the pathways and complexes that maintain genome stability during normal growth and defined changes that occur upon accumulation of DNA lesions due to cytotoxic doses of PARPi. We uncovered unexpected relationships among DNA repair genes, including context-specific buffering interactions between the minimally characterized AUNIP and BRCA1-A complex genes. Our work thus establishes a foundation for mapping differential genetic interactions in mammalian cells and provides a comprehensive resource for future studies of DNA repair and PARP inhibitors.

Single-cell massively-parallel multiplexed microbial sequencing (M3-seq) identifies rare bacterial populations and profiles phage infection.

Bacterial populations are highly adaptive. They can respond to stress and survive in shifting environments. How the behaviours of individual bacteria vary during stress, however, is poorly understood. To identify and characterize rare bacterial subpopulations, technologies for single-cell transcriptional profiling have been developed. Existing approaches show some degree of limitation, for example, in terms of number of cells or transcripts that can be profiled. Due in part to these limitations, few conditions have been studied with these tools. Here we develop massively-parallel, multiplexed, microbial sequencing (M3-seq)-a single-cell RNA-sequencing platform for bacteria that pairs combinatorial cell indexing with post hoc rRNA depletion. We show that M3-seq can profile bacterial cells from different species under a range of conditions in single experiments. We then apply M3-seq to hundreds of thousands of cells, revealing rare populations and insights into bet-hedging associated with stress responses and characterizing phage infection.

Recording morphogen signals reveals origins of gastruloid symmetry breaking.

When cultured in three dimensional spheroids, mammalian stem cells can reproducibly self-organize a single anterior-posterior axis and sequentially differentiate into structures resembling the primitive streak and tailbud. Whereas the embryo's body axes are instructed by spatially patterned extra-embryonic cues, it is unknown how these stem cell gastruloids break symmetry to reproducibly define a single anterior-posterior (A-P) axis. Here, we use synthetic gene circuits to trace how early intracellular signals predict cells' future anterior-posterior position in the gastruloid. We show that Wnt signaling evolves from a homogeneous state to a polarized state, and identify a critical 6-hour time period when single-cell Wnt activity predicts future cellular position, prior to the appearance of polarized signaling patterns or morphology. Single-cell RNA sequencing and live-imaging reveal that early Wnt-high and Wnt-low cells contribute to distinct cell types and suggest that axial symmetry breaking is driven by sorting rearrangements involving differential cell adhesion. We further extend our approach to other canonical embryonic signaling pathways, revealing that even earlier heterogeneity in TGFß signaling predicts A-P position and modulates Wnt signaling during the critical time period. Our study reveals a sequence of dynamic cellular processes that transform a uniform cell aggregate into a polarized structure and demonstrates that a morphological axis can emerge out of signaling heterogeneity and cell movements even in the absence of exogenous patterning cues.

Enhanced prime editing systems by manipulating cellular determinants of editing outcomes.

While prime editing enables precise sequence changes in DNA, cellular determinants of prime editing remain poorly understood. Using pooled CRISPRi screens, we discovered that DNA mismatch repair (MMR) impedes prime editing and promotes undesired indel byproducts. We developed PE4 and PE5 prime editing systems in which transient expression of an engineered MMR-inhibiting protein enhances the efficiency of substitution, small insertion, and small deletion prime edits by an average 7.7-fold and 2.0-fold compared to PE2 and PE3 systems, respectively, while improving edit/indel ratios by 3.4-fold in MMR-proficient cell types. Strategic installation of silent mutations near the intended edit can enhance prime editing outcomes by evading MMR. Prime editor protein optimization resulted in a PEmax architecture that enhances editing efficacy by 2.8-fold on average in HeLa cells. These findings enrich our understanding of prime editing and establish prime editing systems that show substantial improvement across 191 edits in seven mammalian cell types.

Mapping the genetic landscape of DNA double-strand break repair.

Cells repair DNA double-strand breaks (DSBs) through a complex set of pathways critical for maintaining genomic integrity. To systematically map these pathways, we developed a high-throughput screening approach called Repair-seq that measures the effects of thousands of genetic perturbations on mutations introduced at targeted DNA lesions. Using Repair-seq, we profiled DSB repair products induced by two programmable nucleases (Cas9 and Cas12a) in the presence or absence of oligonucleotides for homology-directed repair (HDR) after knockdown of 476 genes involved in DSB repair or associated processes. The resulting data enabled principled, data-driven inference of DSB end joining and HDR pathways. Systematic interrogation of this data uncovered unexpected relationships among DSB repair genes and demonstrated that repair outcomes with superficially similar sequence architectures can have markedly different genetic dependencies. This work provides a foundation for mapping DNA repair pathways and for optimizing genome editing across diverse modalities.

Efficient C•G-to-G•C base editors developed using CRISPRi screens, target-library analysis, and machine learning.

Programmable C•G-to-G•C base editors (CGBEs) have broad scientific and therapeutic potential, but their editing outcomes have proved difficult to predict and their editing efficiency and product purity are often low. We describe a suite of engineered CGBEs paired with machine learning models to enable efficient, high-purity C•G-to-G•C base editing. We performed a CRISPR interference (CRISPRi) screen targeting DNA repair genes to identify factors that affect C•G-to-G•C editing outcomes and used these insights to develop CGBEs with diverse editing profiles. We characterized ten promising CGBEs on a library of 10,638 genomically integrated target sites in mammalian cells and trained machine learning models that accurately predict the purity and yield of editing outcomes (R = 0.90) using these data. These CGBEs enable correction to the wild-type coding sequence of 546 disease-related transversion single-nucleotide variants (SNVs) with >90% precision (mean 96%) and up to 70% efficiency (mean 14%). Computational prediction of optimal CGBE-single-guide RNA pairs enables high-purity transversion base editing at over fourfold more target sites than achieved using any single CGBE variant.

Combinatorial single-cell CRISPR screens by direct guide RNA capture and targeted sequencing.

Single-cell CRISPR screens enable the exploration of mammalian gene function and genetic regulatory networks. However, use of this technology has been limited by reliance on indirect indexing of single-guide RNAs (sgRNAs). Here we present direct-capture Perturb-seq, a versatile screening approach in which expressed sgRNAs are sequenced alongside single-cell transcriptomes. Direct-capture Perturb-seq enables detection of multiple distinct sgRNA sequences from individual cells and thus allows pooled single-cell CRISPR screens to be easily paired with combinatorial perturbation libraries that contain dual-guide expression vectors. We demonstrate the utility of this approach for high-throughput investigations of genetic interactions and, leveraging this ability, dissect epistatic interactions between cholesterol biogenesis and DNA repair. Using direct capture Perturb-seq, we also show that targeting individual genes with multiple sgRNAs per cell improves efficacy of CRISPR interference and activation, facilitating the use of compact, highly active CRISPR libraries for single-cell screens. Last, we show that hybridization-based target enrichment permits sensitive, specific sequencing of informative transcripts from single-cell RNA-seq experiments.

Pervasive functional translation of noncanonical human open reading frames.

Ribosome profiling has revealed pervasive but largely uncharacterized translation outside of canonical coding sequences (CDSs). In this work, we exploit a systematic CRISPR-based screening strategy to identify hundreds of noncanonical CDSs that are essential for cellular growth and whose disruption elicits specific, robust transcriptomic and phenotypic changes in human cells. Functional characterization of the encoded microproteins reveals distinct cellular localizations, specific protein binding partners, and hundreds of microproteins that are presented by the human leukocyte antigen system. We find multiple microproteins encoded in upstream open reading frames, which form stable complexes with the main, canonical protein encoded on the same messenger RNA, thereby revealing the use of functional bicistronic operons in mammals. Together, our results point to a family of functional human microproteins that play critical and diverse cellular roles.

Prime Editing: Precision Genome Editing by Reverse Transcription.

Genome editing is a method for making targeted sequence changes to the genomes of living cells. Prime editing, recently reported by Anzalone et al. (2019), is a new technology that uses reverse transcription to "write" programmed sequence changes into genomic DNA and thus promises significant technical advances.

Molecular recording of mammalian embryogenesis.

Ontogeny describes the emergence of complex multicellular organisms from single totipotent cells. This field is particularly challenging in mammals, owing to the indeterminate relationship between self-renewal and differentiation, variation in progenitor field sizes, and internal gestation in these animals. Here we present a flexible, high-information, multi-channel molecular recorder with a single-cell readout and apply it as an evolving lineage tracer to assemble mouse cell-fate maps from fertilization through gastrulation. By combining lineage information with single-cell RNA sequencing profiles, we recapitulate canonical developmental relationships between different tissue types and reveal the nearly complete transcriptional convergence of endodermal cells of extra-embryonic and embryonic origins. Finally, we apply our cell-fate maps to estimate the number of embryonic progenitor cells and their degree of asymmetric partitioning during specification. Our approach enables massively parallel, high-resolution recording of lineage and other information in mammalian systems, which will facilitate the construction of a quantitative framework for understanding developmental processes.

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.

Perturb-Seq: Dissecting Molecular Circuits with Scalable Single-Cell RNA Profiling of Pooled Genetic Screens.

Genetic screens help infer gene function in mammalian cells, but it has remained difficult to assay complex phenotypes-such as transcriptional profiles-at scale. Here, we develop Perturb-seq, combining single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-based perturbations to perform many such assays in a pool. We demonstrate Perturb-seq by analyzing 200,000 cells in immune cells and cell lines, focusing on transcription factors regulating the response of dendritic cells to lipopolysaccharide (LPS). Perturb-seq accurately identifies individual gene targets, gene signatures, and cell states affected by individual perturbations and their genetic interactions. We posit new functions for regulators of differentiation, the anti-viral response, and mitochondrial function during immune activation. By decomposing many high content measurements into the effects of perturbations, their interactions, and diverse cell metadata, Perturb-seq dramatically increases the scope of pooled genomic assays.

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.

Compact and highly active next-generation libraries for CRISPR-mediated gene repression and activation.

We recently found that nucleosomes directly block access of CRISPR/Cas9 to DNA (Horlbeck et al., 2016). Here, we build on this observation with a comprehensive algorithm that incorporates chromatin, position, and sequence features to accurately predict highly effective single guide RNAs (sgRNAs) for targeting nuclease-dead Cas9-mediated transcriptional repression (CRISPRi) and activation (CRISPRa). We use this algorithm to design next-generation genome-scale CRISPRi and CRISPRa libraries targeting human and mouse genomes. A CRISPRi screen for essential genes in K562 cells demonstrates that the large majority of sgRNAs are highly active. We also find CRISPRi does not exhibit any detectable non-specific toxicity recently observed with CRISPR nuclease approaches. Precision-recall analysis shows that we detect over 90% of essential genes with minimal false positives using a compact 5 sgRNA/gene library. Our results establish CRISPRi and CRISPRa as premier tools for loss- or gain-of-function studies and provide a general strategy for identifying Cas9 target sites.

A Systematic Analysis of Factors Localized to Damaged Chromatin Reveals PARP-Dependent Recruitment of Transcription Factors.

Localization to sites of DNA damage is a hallmark of DNA damage response (DDR) proteins. To identify DDR factors, we screened epitope-tagged proteins for localization to sites of chromatin damaged by UV laser microirradiation and found >120 proteins that localize to damaged chromatin. These include the BAF tumor suppressor complex and the amyotrophic lateral sclerosis (ALS) candidate protein TAF15. TAF15 contains multiple domains that bind damaged chromatin in a poly-(ADP-ribose) polymerase (PARP)-dependent manner, suggesting a possible role as glue that tethers multiple PAR chains together. Many positives were transcription factors; > 70% of randomly tested transcription factors localized to sites of DNA damage, and of these, ∼90% were PARP dependent for localization. Mutational analyses showed that localization to damaged chromatin is DNA-binding-domain dependent. By examining Hoechst staining patterns at damage sites, we see evidence of chromatin decompaction that is PARP dependent. We propose that PARP-regulated chromatin remodeling at sites of damage allows transient accessibility of DNA-binding proteins.