Joint genotypic and phenotypic outcome modeling improves base editing variant effect quantification.

CRISPR base editing screens enable analysis of disease-associated variants at scale; however, variable efficiency and precision confounds the assessment of variant-induced phenotypes. Here, we provide an integrated experimental and computational pipeline that improves estimation of variant effects in base editing screens. We use a reporter construct to measure guide RNA (gRNA) editing outcomes alongside their phenotypic consequences and introduce base editor screen analysis with activity normalization (BEAN), a Bayesian network that uses per-guide editing outcomes provided by the reporter and target site chromatin accessibility to estimate variant impacts. BEAN outperforms existing tools in variant effect quantification. We use BEAN to pinpoint common regulatory variants that alter low-density lipoprotein (LDL) uptake, implicating previously unreported genes. Additionally, through saturation base editing of LDLR, we accurately quantify missense variant pathogenicity that is consistent with measurements in UK Biobank patients and identify underlying structural mechanisms. This work provides a widely applicable approach to improve the power of base editing screens for disease-associated variant characterization.

Joint genotypic and phenotypic outcome modeling improves base editing variant effect quantification.

CRISPR base editing screens are powerful tools for studying disease-associated variants at scale. However, the efficiency and precision of base editing perturbations vary, confounding the assessment of variant-induced phenotypic effects. Here, we provide an integrated pipeline that improves the estimation of variant impact in base editing screens. We perform high-throughput ABE8e-SpRY base editing screens with an integrated reporter construct to measure the editing efficiency and outcomes of each gRNA alongside their phenotypic consequences. We introduce BEAN, a Bayesian network that accounts for per-guide editing outcomes and target site chromatin accessibility to estimate variant impacts. We show this pipeline attains superior performance compared to existing tools in variant classification and effect size quantification. We use BEAN to pinpoint common variants that alter LDL uptake, implicating novel genes. Additionally, through saturation base editing of LDLR, we enable accurate quantitative prediction of the effects of missense variants on LDL-C levels, which aligns with measurements in UK Biobank individuals, and identify structural mechanisms underlying variant pathogenicity. This work provides a widely applicable approach to improve the power of base editor screens for disease-associated variant characterization.

Systematic elucidation of genetic mechanisms underlying cholesterol uptake.

Genetic variation contributes greatly to LDL cholesterol (LDL-C) levels and coronary artery disease risk. By combining analysis of rare coding variants from the UK Biobank and genome-scale CRISPR-Cas9 knockout and activation screening, we substantially improve the identification of genes whose disruption alters serum LDL-C levels. We identify 21 genes in which rare coding variants significantly alter LDL-C levels at least partially through altered LDL-C uptake. We use co-essentiality-based gene module analysis to show that dysfunction of the RAB10 vesicle transport pathway leads to hypercholesterolemia in humans and mice by impairing surface LDL receptor levels. Further, we demonstrate that loss of function of OTX2 leads to robust reduction in serum LDL-C levels in mice and humans by increasing cellular LDL-C uptake. Altogether, we present an integrated approach that improves our understanding of the genetic regulators of LDL-C levels and provides a roadmap for further efforts to dissect complex human disease genetics.

A molecular glue approach to control the half-life of CRISPR-based technologies.

Cas9 is a programmable nuclease that has furnished transformative technologies, including base editors and transcription modulators (e.g., CRISPRi/a), but several applications of these technologies, including therapeutics, mandatorily require precision control of their half-life. For example, such control can help avert any potential immunological and adverse events in clinical trials. Current genome editing technologies to control the half-life of Cas9 are slow, have lower activity, involve fusion of large response elements (> 230 amino acids), utilize expensive controllers with poor pharmacological attributes, and cannot be implemented in vivo on several CRISPR-based technologies. We report a general platform for half-life control using the molecular glue, pomalidomide, that binds to a ubiquitin ligase complex and a response-element bearing CRISPR-based technology, thereby causing the latter's rapid ubiquitination and degradation. Using pomalidomide, we were able to control the half-life of large CRISPR-based technologies (e.g., base editors, CRISPRi) and small anti-CRISPRs that inhibit such technologies, allowing us to build the first examples of on-switch for base editors. The ability to switch on, fine-tune and switch-off CRISPR-based technologies with pomalidomide allowed complete control over their activity, specificity, and genome editing outcome. Importantly, the miniature size of the response element and favorable pharmacological attributes of the drug pomalidomide allowed control of activity of base editor in vivo using AAV as the delivery vehicle. These studies provide methods and reagents to precisely control the dosage and half-life of CRISPR-based technologies, propelling their therapeutic development.

Systematic elucidation of genetic mechanisms underlying cholesterol uptake.

Genetic variation contributes greatly to LDL cholesterol (LDL-C) levels and coronary artery disease risk. By combining analysis of rare coding variants from the UK Biobank and genome-scale CRISPR-Cas9 knockout and activation screening, we have substantially improved the identification of genes whose disruption alters serum LDL-C levels. We identify 21 genes in which rare coding variants significantly alter LDL-C levels at least partially through altered LDL-C uptake. We use co-essentiality-based gene module analysis to show that dysfunction of the RAB10 vesicle transport pathway leads to hypercholesterolemia in humans and mice by impairing surface LDL receptor levels. Further, we demonstrate that loss of function of OTX2 leads to robust reduction in serum LDL-C levels in mice and humans by increasing cellular LDL-C uptake. Altogether, we present an integrated approach that improves our understanding of genetic regulators of LDL-C levels and provides a roadmap for further efforts to dissect complex human disease genetics.

Peptide fusion improves prime editing efficiency.

Prime editing enables search-and-replace genome editing but is limited by low editing efficiency. We present a high-throughput approach, the Peptide Self-Editing sequencing assay (PepSEq), to measure how fusion of 12,000 85-amino acid peptides influences prime editing efficiency. We show that peptide fusion can enhance prime editing, prime-enhancing peptides combine productively, and a top dual peptide-prime editor increases prime editing significantly in multiple cell lines across dozens of target sites. Top prime-enhancing peptides function by increasing translation efficiency and serve as broadly useful tools to improve prime editing efficiency.

Small molecule inhibition of ATM kinase increases CRISPR-Cas9 1-bp insertion frequency.

Mutational outcomes following CRISPR-Cas9-nuclease cutting in mammalian cells have recently been shown to be predictable and, in certain cases, skewed toward single genotypes. However, the ability to control these outcomes remains limited, especially for 1-bp insertions, a common and therapeutically relevant class of repair outcomes. Here, through a small molecule screen, we identify the ATM kinase inhibitor KU-60019 as a compound capable of reproducibly increasing the fraction of 1-bp insertions relative to other Cas9 repair outcomes. Small molecule or genetic ATM inhibition increases 1-bp insertion outcome fraction across three human and mouse cell lines, two Cas9 species, and dozens of target sites, although concomitantly reducing the fraction of edited alleles. Notably, KU-60019 increases the relative frequency of 1-bp insertions to over 80% of edited alleles at several native human genomic loci and improves the efficiency of correction for pathogenic 1-bp deletion variants. The ability to increase 1-bp insertion frequency adds another dimension to precise template-free Cas9-nuclease genome editing.

Using CRISPR to understand and manipulate gene regulation.

Understanding how genes are expressed in the correct cell types and at the correct level is a key goal of developmental biology research. Gene regulation has traditionally been approached largely through observational methods, whereas perturbational approaches have lacked precision. CRISPR-Cas9 has begun to transform the study of gene regulation, allowing for precise manipulation of genomic sequences, epigenetic functionalization and gene expression. CRISPR-Cas9 technology has already led to the discovery of new paradigms in gene regulation and, as new CRISPR-based tools and methods continue to be developed, promises to transform our knowledge of the gene regulatory code and our ability to manipulate cell fate. Here, we discuss the current and future application of the emerging CRISPR toolbox toward predicting gene regulatory network behavior, improving stem cell disease modeling, dissecting the epigenetic code, reprogramming cell fate and treating diseases of gene dysregulation.

Detection of gene cis-regulatory element perturbations in single-cell transcriptomes.

We introduce poly-adenine CRISPR gRNA-based single-cell RNA-sequencing (pAC-Seq), a method that enables the direct observation of guide RNAs (gRNAs) in scRNA-seq. We use pAC-Seq to assess the phenotypic consequences of CRISPR/Cas9 based alterations of gene cis-regulatory regions. We show that pAC-Seq is able to detect cis-regulatory-induced alteration of target gene expression even when biallelic loss of target gene expression occurs in only ~5% of cells. This low rate of biallelic loss significantly increases the number of cells required to detect the consequences of changes to the regulatory genome, but can be ameliorated by transcript-targeted sequencing. Based on our experimental results we model the power to detect regulatory genome induced transcriptomic effects based on the rate of mono/biallelic loss, baseline gene expression, and the number of cells per target gRNA.

Machine learning based CRISPR gRNA design for therapeutic exon skipping.

Restoring gene function by the induced skipping of deleterious exons has been shown to be effective for treating genetic disorders. However, many of the clinically successful therapies for exon skipping are transient oligonucleotide-based treatments that require frequent dosing. CRISPR-Cas9 based genome editing that causes exon skipping is a promising therapeutic modality that may offer permanent alleviation of genetic disease. We show that machine learning can select Cas9 guide RNAs that disrupt splice acceptors and cause the skipping of targeted exons. We experimentally measured the exon skipping frequencies of a diverse genome-integrated library of 791 splice sequences targeted by 1,063 guide RNAs in mouse embryonic stem cells. We found that our method, SkipGuide, is able to identify effective guide RNAs with a precision of 0.68 (50% threshold predicted exon skipping frequency) and 0.93 (70% threshold predicted exon skipping frequency). We anticipate that SkipGuide will be useful for selecting guide RNA candidates for evaluation of CRISPR-Cas9-mediated exon skipping therapy.

Chemogenetic System Demonstrates That Cas9 Longevity Impacts Genome Editing Outcomes.

Prolonged Cas9 activity can hinder genome engineering as it causes off-target effects, genotoxicity, heterogeneous genome-editing outcomes, immunogenicity, and mosaicism in embryonic editing-issues which could be addressed by controlling the longevity of Cas9. Though some temporal controls of Cas9 activity have been developed, only cumbersome systems exist for modifying the lifetime. Here, we have developed a chemogenetic system that brings Cas9 in proximity to a ubiquitin ligase, enabling rapid ubiquitination and degradation of Cas9 by the proteasome. Despite the large size of Cas9, we were able to demonstrate efficient degradation in cells from multiple species. Furthermore, by controlling the Cas9 lifetime, we were able to bias the DNA repair pathways and the genotypic outcome for both templated and nontemplated genome editing. Finally, we were able to dosably control the Cas9 activity and specificity to ameliorate the off-target effects. The ability of this system to change the Cas9 lifetime and, therefore, bias repair pathways and specificity in the desired direction allows precision control of the genome editing outcome.

Comprehensive Mapping of Key Regulatory Networks that Drive Oncogene Expression.

Gene expression is controlled by the collective binding of transcription factors to cis-regulatory regions. Deciphering gene-centered regulatory networks is vital to understanding and controlling gene misexpression in human disease; however, systematic approaches to uncovering regulatory networks have been lacking. Here we present high-throughput interrogation of gene-centered activation networks (HIGAN), a pipeline that employs a suite of multifaceted genomic approaches to connect upstream signaling inputs, trans-acting TFs, and cis-regulatory elements. We apply HIGAN to understand the aberrant activation of the cytidine deaminase APOBEC3B, an intrinsic source of cancer hypermutation. We reveal that nuclear factor κB (NF-κB) and AP-1 pathways are the most salient trans-acting inputs, with minor roles for other inflammatory pathways. We identify a cis-regulatory architecture dominated by a major intronic enhancer that requires coordinated NF-κB and AP-1 activity with secondary inputs from distal regulatory regions. Our data demonstrate how integration of cis and trans genomic screening platforms provides a paradigm for building gene-centered regulatory networks.

Elucidation of remdesivir cytotoxicity pathways through genome-wide CRISPR-Cas9 screening and transcriptomics.

The adenosine analogue remdesivir has emerged as a front-line antiviral treatment for SARS-CoV-2, with preliminary evidence that it reduces the duration and severity of illness1.Prior clinical studies have identified adverse events1,2, and remdesivir has been shown to inhibit mitochondrial RNA polymerase in biochemical experiments7, yet little is known about the specific genetic pathways involved in cellular remdesivir metabolism and cytotoxicity. Through genome-wide CRISPR-Cas9 screening and RNA sequencing, we show that remdesivir treatment leads to a repression of mitochondrial respiratory activity, and we identify five genes whose loss significantly reduces remdesivir cytotoxicity. In particular, we show that loss of the mitochondrial nucleoside transporter SLC29A3 mitigates remdesivir toxicity without a commensurate decrease in SARS-CoV-2 antiviral potency and that the mitochondrial adenylate kinase AK2 is a remdesivir kinase required for remdesivir efficacy and toxicity. This work elucidates the cellular mechanisms of remdesivir metabolism and provides a candidate gene target to reduce remdesivir cytotoxicity.

Identification of determinants of differential chromatin accessibility through a massively parallel genome-integrated reporter assay.

A key mechanism in cellular regulation is the ability of the transcriptional machinery to physically access DNA. Transcription factors interact with DNA to alter the accessibility of chromatin, which enables changes to gene expression during development or disease or as a response to environmental stimuli. However, the regulation of DNA accessibility via the recruitment of transcription factors is difficult to study in the context of the native genome because every genomic site is distinct in multiple ways. Here we introduce the multiplexed integrated accessibility assay (MIAA), an assay that measures chromatin accessibility of synthetic oligonucleotide sequence libraries integrated into a controlled genomic context with low native accessibility. We apply MIAA to measure the effects of sequence motifs on cell type-specific accessibility between mouse embryonic stem cells and embryonic stem cell-derived definitive endoderm cells, screening 7905 distinct DNA sequences. MIAA recapitulates differential accessibility patterns of 100-nt sequences derived from natively differential genomic regions, identifying E-box motifs common to epithelial-mesenchymal transition driver transcription factors in stem cell-specific accessible regions that become repressed in endoderm. We show that a single binding motif for a key regulatory transcription factor is sufficient to open chromatin, and classify sets of stem cell-specific, endoderm-specific, and shared accessibility-modifying transcription factor motifs. We also show that overexpression of two definitive endoderm transcription factors, T and Foxa2, results in changes to accessibility in DNA sequences containing their respective DNA-binding motifs and identify preferential motif arrangements that influence accessibility.

A High-Throughput Genome-Integrated Assay Reveals Spatial Dependencies Governing Tcf7l2 Binding.

Predicting where transcription factors bind in the genome from their in vitro DNA-binding affinity is confounded by the large number of possible interactions with nearby transcription factors. To characterize the in vivo binding logic for the Wnt effector Tcf7l2, we developed a high-throughput screening platform in which thousands of synthesized DNA phrases are inserted into a specific genomic locus, followed by measurement of Tcf7l2 binding by DamID. Using this platform at two genomic loci in mouse embryonic stem cells, we show that while the binding of Tcf7l2 closely follows the in vitro motif-binding strength and is influenced by local chromatin accessibility, it is also strongly affected by the surrounding 99 bp of sequence. Through controlled sequence perturbation, we show that Oct4 and Klf4 motifs promote Tcf7l2 binding, particularly in the adjacent ∼50 bp and oscillating with a 10.8-bp phasing relative to these cofactor motifs, which matches the turn of a DNA helix.

A Multiplexed Barcodelet Single-Cell RNA-Seq Approach Elucidates Combinatorial Signaling Pathways that Drive ESC Differentiation.

Empirical optimization of stem cell differentiation protocols is time consuming, is laborintensive, and typically does not comprehensively interrogate all relevant signaling pathways. Here we describe barcodelet single-cell RNA sequencing (barRNA-seq), which enables systematic exploration of cellular perturbations by tagging individual cells with RNA "barcodelets" to identify them on the basis of the treatments they receive. We apply barRNA-seq to simultaneously manipulate up to seven developmental pathways and study effects on embryonic stem cell (ESC) germ layer specification and mesodermal specification, uncovering combinatorial effects of signaling pathway activation on gene expression. We further develop a data-driven framework for identifying combinatorial signaling perturbations that drive cells toward specific fates, including several annotated in an existing scRNA-seq gastrulation atlas, and use this approach to guide ESC differentiation into a notochord-like population. We expect that barRNA-seq will have broad utility for investigating and understanding how cooperative signaling pathways drive cell fate acquisition.

Author Correction: Predictable and precise template-free CRISPR editing of pathogenic variants.

In this Article, a data processing error affected Fig. 3e and Extended Data Table 2; these errors have been corrected online.

Predictable and precise template-free CRISPR editing of pathogenic variants.

Following Cas9 cleavage, DNA repair without a donor template is generally considered stochastic, heterogeneous and impractical beyond gene disruption. Here, we show that template-free Cas9 editing is predictable and capable of precise repair to a predicted genotype, enabling correction of disease-associated mutations in humans. We constructed a library of 2,000 Cas9 guide RNAs paired with DNA target sites and trained inDelphi, a machine learning model that predicts genotypes and frequencies of 1- to 60-base-pair deletions and 1-base-pair insertions with high accuracy (r = 0.87) in five human and mouse cell lines. inDelphi predicts that 5-11% of Cas9 guide RNAs targeting the human genome are 'precise-50', yielding a single genotype comprising greater than or equal to 50% of all major editing products. We experimentally confirmed precise-50 insertions and deletions in 195 human disease-relevant alleles, including correction in primary patient-derived fibroblasts of pathogenic alleles to wild-type genotype for Hermansky-Pudlak syndrome and Menkes disease. This study establishes an approach for precise, template-free genome editing.

A CRISPR view of gene regulation.

Due to plummeting costs, whole genome sequencing of patients and cancers will soon become routine medical practice; however, we cannot currently predict how non-coding genotype affects cellular gene expression. Gene regulation research has recently been dominated by observational approaches that correlate chromatin state with regulatory function. These approaches are limited to the available genotypes and cannot scratch the surface of possible sequence combinations, and thus there is a need for perturbation-based approaches to better understand how DNA encodes gene regulatory functions. CRISPR/Cas9 genome editing has revolutionized our ability to alter genome sequence, and CRISPR/Cas9-based assays have already begun to contribute to new paradigms of gene regulation. We discuss the variety of arenas in which current and future CRISPR-based technologies will aid in developing predictive understanding of how genome sequence leads to gene regulatory function.

A CD47-associated super-enhancer links pro-inflammatory signalling to CD47 upregulation in breast cancer.

CD47 is a cell surface molecule that inhibits phagocytosis of cells that express it by binding to its receptor, SIRPα, on macrophages and other immune cells. CD47 is expressed at different levels by neoplastic and normal cells. Here, to reveal mechanisms by which different neoplastic cells generate this dominant 'don't eat me' signal, we analyse the CD47 regulatory genomic landscape. We identify two distinct super-enhancers (SEs) associated with CD47 in certain cancer cell types. We show that a set of active constituent enhancers, located within the two CD47 SEs, regulate CD47 expression in different cancer cell types and that disruption of CD47 SEs reduces CD47 gene expression. Finally we report that the TNF-NFKB1 signalling pathway directly regulates CD47 by interacting with a constituent enhancer located within a CD47-associated SE specific to breast cancer. These results suggest that cancers can evolve SE to drive CD47 overexpression to escape immune surveillance.