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.

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.

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.

Genome-Scale CRISPR-Mediated Control of Gene Repression and Activation.

While the catalog of mammalian transcripts and their expression levels in different cell types and disease states is rapidly expanding, our understanding of transcript function lags behind. We present a robust technology enabling systematic investigation of the cellular consequences of repressing or inducing individual transcripts. We identify rules for specific targeting of transcriptional repressors (CRISPRi), typically achieving 90%-99% knockdown with minimal off-target effects, and activators (CRISPRa) to endogenous genes via endonuclease-deficient Cas9. Together they enable modulation of gene expression over a ∼1,000-fold range. Using these rules, we construct genome-scale CRISPRi and CRISPRa libraries, each of which we validate with two pooled screens. Growth-based screens identify essential genes, tumor suppressors, and regulators of differentiation. Screens for sensitivity to a cholera-diphtheria toxin provide broad insights into the mechanisms of pathogen entry, retrotranslocation and toxicity. Our results establish CRISPRi and CRISPRa as powerful tools that provide rich and complementary information for mapping complex pathways.

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.

Polyubiquitinated PCNA recruits the ZRANB3 translocase to maintain genomic integrity after replication stress.

Completion of DNA replication after replication stress depends on PCNA, which undergoes monoubiquitination to stimulate direct bypass of DNA lesions by specialized DNA polymerases or is polyubiquitinated to promote recombination-dependent DNA synthesis across DNA lesions by template switching mechanisms. Here we report that the ZRANB3 translocase, a SNF2 family member related to the SIOD disorder SMARCAL1 protein, is recruited by polyubiquitinated PCNA to promote fork restart following replication arrest. ZRANB3 depletion in mammalian cells results in an increased frequency of sister chromatid exchange and DNA damage sensitivity after treatment with agents that cause replication stress. Using in vitro biochemical assays, we show that recombinant ZRANB3 remodels DNA structures mimicking stalled replication forks and disassembles recombination intermediates. We therefore propose that ZRANB3 maintains genomic stability at stalled or collapsed replication forks by facilitating fork restart and limiting inappropriate recombination that could occur during template switching events.

A genome-wide camptothecin sensitivity screen identifies a mammalian MMS22L-NFKBIL2 complex required for genomic stability.

Replication stress involving collision of replisomes with camptothecin (CPT)-stabilized DNA-Topoisomerase I adducts activates an ATR-dependent pathway to promote repair by homologous recombination. To identify human genes that protect cells from such replication stress, we performed a genome-wide CPT sensitivity screen. Among numerous candidate genes are two previously unstudied proteins: the ankyrin repeat protein NFKBIL2 and C6ORF167 (MMS22L), distantly related to yeast replication stress regulator Mms22p. MMS22L and NFKBIL2 interact with each other and with FACT (facilitator of chromatin transcription) and MCM (minichromosome maintenance) complexes. Cells depleted of NFKBIL2 or MMS22L are sensitive to DNA-damaging agents, load phosphorylated RPA onto chromatin in a CTIP-dependent manner, activate the ATR/ATRIP-CHK1 and double-strand break repair signaling pathways, and are defective in HR. This study identifies MMS22L-NFKBIL2 as components of the replication stress control pathway and provides a resource for discovery of additional components of this pathway.

A bioinformatics method identifies prominent off-targeted transcripts in RNAi screens.

Because off-target effects hamper interpretation and validation of RNAi screen data, we developed a bioinformatics method, genome-wide enrichment of seed sequence matches (GESS), to identify candidate off-targeted transcripts in primary screening data. GESS analysis revealed a prominent off-targeted transcript in several screens, including MAD2 (MAD2L1) in a screen for genes required for the spindle assembly checkpoint. GESS analysis results can enhance the validation rate in RNAi screens.

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.

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.

A chromatin localization screen reveals poly (ADP ribose)-regulated recruitment of the repressive polycomb and NuRD complexes to sites of DNA damage.

Many proteins that respond to DNA damage are recruited to DNA lesions. We used a proteomics approach that coupled isotopic labeling with chromatin fractionation and mass spectrometry to uncover proteins that associate with damaged DNA, many of which are involved in DNA repair or nucleolar function. We show that polycomb group members are recruited by poly(ADP ribose) polymerase (PARP) to DNA lesions following UV laser microirradiation. Loss of polycomb components results in IR sensitivity of mammalian cells and Caenorhabditis elegans. PARP also recruits two components of the repressive nucleosome remodeling and deacetylase (NuRD) complex, chromodomain helicase DNA-binding protein 4 (CHD4) and metastasis associated 1 (MTA1), to DNA lesions. PARP plays a role in removing nascent RNA and elongating RNA polymerase II from sites of DNA damage. We propose that PARP sets up a transient repressive chromatin structure at sites of DNA damage to block transcription and facilitate DNA repair.

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.

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.

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.

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.

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.