Identification of pathogenic variants in cancer genes using base editing screens with editing efficiency correction

Predicting functional UTR variants by integrating region-specific features

The untranslated region (UTR) of messenger ribonucleic acid (mRNA), including the 5'UTR and 3'UTR, plays a critical role in regulating gene expression and translation. Variants within the UTR can lead to changes associated with human traits and diseases; however, computational prediction of UTR variant effect is challenging. Current noncoding variant prediction mainly focuses on the promoters and enhancers, neglecting the unique sequence of the UTR and thereby limiting their predictive accuracy. In this study, using consolidated datasets of UTR variants from disease databases and large-scale experimental data, we systematically analyzed more than 50 region-specific features of UTR, including functional elements, secondary structure, sequence composition and site conservation. Our analysis reveals that certain features, such as C/G-related sequence composition in 5'UTR and A/T-related sequence composition in 3'UTR, effectively differentiate between nonfunctional and functional variant sets, unveiling potential sequence determinants of functional UTR variants. Leveraging these insights, we developed two classification models to predict functional UTR variants using machine learning, achieving an area under the curve (AUC) value of 0.94 for 5'UTR and 0.85 for 3'UTR, outperforming all existing methods. Our models will be valuable for enhancing clinical interpretation of genetic variants, facilitating the prediction and management of disease risk.

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.

Deconstructing cancer with precision genome editing

Recent advances in genome editing technologies are allowing investigators to engineer and study cancer-associated mutations in their endogenous genetic contexts with high precision and efficiency. Of these, base editing and prime editing are quickly becoming gold-standards in the field due to their versatility and scalability. Here, we review the merits and limitations of these precision genome editing technologies, their application to modern cancer research, and speculate how these could be integrated to address future directions in the field.

Recent advances in CRISPR-based functional genomics for the study of disease-associated genetic variants

Advances in sequencing technology have greatly increased our ability to gather genomic data, yet understanding the impact of genetic mutations, particularly variants of uncertain significance (VUSs), remains a challenge in precision medicine. The CRISPR‒Cas system has emerged as a pivotal tool for genome engineering, enabling the precise incorporation of specific genetic variations, including VUSs, into DNA to facilitate their functional characterization. Additionally, the integration of CRISPR‒Cas technology with sequencing tools allows the high-throughput evaluation of mutations, transforming uncertain genetic data into actionable insights. This allows researchers to comprehensively study the functional consequences of point mutations, paving the way for enhanced understanding and increasing application to precision medicine. This review summarizes the current genome editing tools utilizing CRISPR‒Cas systems and their combination with sequencing tools for functional genomics, with a focus on point mutations.

The Current Situation and Development Prospect of Whole-Genome Screening

High-throughput genetic screening is useful for discovering critical genes or gene sequences that trigger specific cell functions and/or phenotypes. Loss-of-function genetic screening is mainly achieved through RNA interference (RNAi), CRISPR knock-out (CRISPRko), and CRISPR interference (CRISPRi) technologies. Gain-of-function genetic screening mainly depends on the overexpression of a cDNA library and CRISPR activation (CRISPRa). Base editing can perform both gain- and loss-of-function genetic screening. This review discusses genetic screening techniques based on Cas9 nuclease, including Cas9-mediated genome knock-out and dCas9-based gene activation and interference. We compare these methods with previous genetic screening techniques based on RNAi and cDNA library overexpression and propose future prospects and applications for CRISPR screening.

Deep mutational scanning of proteins in mammalian cells

Protein mutagenesis is essential for unveiling the molecular mechanisms underlying protein function in health, disease, and evolution. In the past decade, deep mutational scanning methods have evolved to support the functional analysis of nearly all possible single-amino acid changes in a protein of interest. While historically these methods were developed in lower organisms such as E. coli and yeast, recent technological advancements have resulted in the increased use of mammalian cells, particularly for studying proteins involved in human disease. These advancements will aid significantly in the classification and interpretation of variants of unknown significance, which are being discovered at large scale due to the current surge in the use of whole-genome sequencing in clinical contexts. Here, we explore the experimental aspects of deep mutational scanning studies in mammalian cells and report the different methods used in each step of the workflow, ultimately providing a useful guide toward the design of such studies.

A split and inducible adenine base editor for precise in vivo base editing

DNA base editors use deaminases fused to a programmable DNA-binding protein for targeted nucleotide conversion. However, the most widely used TadA deaminases lack post-translational control in living cells. Here, we present a split adenine base editor (sABE) that utilizes chemically induced dimerization (CID) to control the catalytic activity of the deoxyadenosine deaminase TadA-8e. sABE shows high on-target editing activity comparable to the original ABE with TadA-8e (ABE8e) upon rapamycin induction while maintaining low background activity without induction. Importantly, sABE exhibits a narrower activity window on DNA and higher precision than ABE8e, with an improved single-to-double ratio of adenine editing and reduced genomic and transcriptomic off-target effects. sABE can achieve gene knockout through multiplex splice donor disruption in human cells. Furthermore, when delivered via dual adeno-associated virus vectors, sABE can efficiently convert a single A•T base pair to a G•C base pair on the PCSK9 gene in mouse liver, demonstrating in vivo CID-controlled DNA base editing. Thus, sABE enables precise control of base editing, which will have broad implications for basic research and in vivo therapeutic applications.

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.

Saturation genome editing of 11 codons and exon 13 of BRCA2 coupled with chemotherapeutic drug response accurately determines pathogenicity of variants

The unknown pathogenicity of a significant number of variants found in cancer-related genes is attributed to limited epidemiological data, resulting in their classification as variant of uncertain significance (VUS). To date, Breast Cancer gene-2 (BRCA2) has the highest number of VUSs, which has necessitated the development of several robust functional assays to determine their functional significance. Here we report the use of a humanized-mouse embryonic stem cell (mESC) line expressing a single copy of the human BRCA2 for a CRISPR-Cas9-based high-throughput functional assay. As a proof-of-principle, we have saturated 11 codons encoded by BRCA2 exons 3, 18, 19 and all possible single-nucleotide variants in exon 13 and multiplexed these variants for their functional categorization. Specifically, we used a pool of 180-mer single-stranded donor DNA to generate all possible combination of variants. Using a high throughput sequencing-based approach, we show a significant drop in the frequency of non-functional variants, whereas functional variants are enriched in the pool of the cells. We further demonstrate the response of these variants to the DNA-damaging agents, cisplatin and olaparib, allowing us to use cellular survival and drug response as parameters for variant classification. Using this approach, we have categorized 599 BRCA2 variants including 93-single nucleotide variants (SNVs) across the 11 codons, of which 28 are reported in ClinVar. We also functionally categorized 252 SNVs from exon 13 into 188 functional and 60 non-functional variants, demonstrating that saturation genome editing (SGE) coupled with drug sensitivity assays can enhance functional annotation of BRCA2 VUS.

Base editor screens for in situ mutational scanning at scale

A fundamental challenge in biology is understanding the molecular details of protein function. How mutations alter protein activity, regulation, and response to drugs is of critical importance to human health. Recent years have seen the emergence of pooled base editor screens for in situ mutational scanning: the interrogation of protein sequence-function relationships by directly perturbing endogenous proteins in live cells. These studies have revealed the effects of disease-associated mutations, discovered novel drug resistance mechanisms, and generated biochemical insights into protein function. Here, we discuss how this "base editor scanning" approach has been applied to diverse biological questions, compare it with alternative techniques, and describe the emerging challenges that must be addressed to maximize its utility. Given its broad applicability toward profiling mutations across the proteome, base editor scanning promises to revolutionize the investigation of proteins in their native contexts.

Revolutionizing DNA repair research and cancer therapy with CRISPR-Cas screens

All organisms possess molecular mechanisms that govern DNA repair and associated DNA damage response (DDR) processes. Owing to their relevance to human disease, most notably cancer, these mechanisms have been studied extensively, yet new DNA repair and/or DDR factors and functional interactions between them are still being uncovered. The emergence of CRISPR technologies and CRISPR-based genetic screens has enabled genome-scale analyses of gene-gene and gene-drug interactions, thereby providing new insights into cellular processes in distinct DDR-deficiency genetic backgrounds and conditions. In this Review, we discuss the mechanistic basis of CRISPR-Cas genetic screening approaches and describe how they have contributed to our understanding of DNA repair and DDR pathways. We discuss how DNA repair pathways are regulated, and identify and characterize crosstalk between them. We also highlight the impacts of CRISPR-based studies in identifying novel strategies for cancer therapy, and in understanding, overcoming and even exploiting cancer-drug resistance, for example in the contexts of PARP inhibition, homologous recombination deficiencies and/or replication stress. Lastly, we present the DDR CRISPR screen (DDRcs) portal , in which we have collected and reanalysed data from CRISPR screen studies and provide a tool for systematically exploring them.

Base editors: development and applications in biomedicine

Base editor (BE) is a gene-editing tool developed by combining the CRISPR/Cas system with an individual deaminase, enabling precise single-base substitution in DNA or RNA without generating a DNA double-strand break (DSB) or requiring donor DNA templates in living cells. Base editors offer more precise and secure genome-editing effects than other conventional artificial nuclease systems, such as CRISPR/Cas9, as the DSB induced by Cas9 will cause severe damage to the genome. Thus, base editors have important applications in the field of biomedicine, including gene function investigation, directed protein evolution, genetic lineage tracing, disease modeling, and gene therapy. Since the development of the two main base editors, cytosine base editors (CBEs) and adenine base editors (ABEs), scientists have developed more than 100 optimized base editors with improved editing efficiency, precision, specificity, targeting scope, and capacity to be delivered in vivo, greatly enhancing their application potential in biomedicine. Here, we review the recent development of base editors, summarize their applications in the biomedical field, and discuss future perspectives and challenges for therapeutic applications.

Joint estimation and imputation of variant functional effects using high throughput assay data

Deep mutational scanning assays enable the functional assessment of variants in high throughput. Phenotypic measurements from these assays are broadly concordant with clinical outcomes but are prone to noise at the individual variant level. We develop a framework to exploit related measurements within and across experimental assays to jointly estimate variant impact. Drawing from a large corpus of deep mutational scanning data, we collectively estimate the mean functional effect per AA residue position within each gene, normalize observed functional effects by substitution type, and make estimates for individual allelic variants with a pipeline called FUSE (Functional Substitution Estimation). FUSE improves the correlation of functional screening datasets covering the same variants, better separates estimated functional impacts for known pathogenic and benign variants (ClinVar BRCA1, p=2.24×10-51), and increases the number of variants for which predictions can be made (2,741 to 10,347) by inferring additional variant effects for substitutions not experimentally screened. For UK Biobank patients who carry a rare variant in TP53, FUSE significantly improves the separation of patients who develop cancer syndromes from those without cancer (p=1.77×10-6). These approaches promise to improve estimates of variant impact and broaden the utility of screening data generated from functional assays.

The CRISPR Revolution in the Drug Discovery Workflow: An Industry Perspective

In a relatively short time, the pharmaceutical industry has witnessed a rapid integration of the CRISPR technology in multiple areas of research, development, therapy, and diagnostics. A unique feature to this system compared with other technologies is the exceptional versatility in adapting to the broad range of needs across the drug discovery pipeline, such as target identification, cell engineering, and in vivo modeling. As a consequence, the CRISPR toolbox has been evolving to address key questions around preclinical and clinical drug development. In this review, we provide a high-level perspective of how CRISPR has impacted several aspects of the drug discovery workflow and the future ahead for this exciting technology.

CRISPR screening in cancer stem cells

Cancer stem cells (CSCs) are a subpopulation of tumor cells with self-renewal ability. Increasing evidence points to the critical roles of CSCs in tumorigenesis, metastasis, therapy resistance, and cancer relapse. As such, the elimination of CSCs improves cancer treatment outcomes. However, challenges remain due to limited understanding of the molecular mechanisms governing self-renewal and survival of CSCs. Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 screening has been increasingly used to identify genetic determinants in cancers. In this primer, we discuss the progress made and emerging opportunities of coupling advanced CRISPR screening systems with CSC models to reveal the understudied vulnerabilities of CSCs.

CRISPR-based genome editing through the lens of DNA repair

Genome editing technologies operate by inducing site-specific DNA perturbations that are resolved by cellular DNA repair pathways. Products of genome editors include DNA breaks generated by CRISPR-associated nucleases, base modifications induced by base editors, DNA flaps created by prime editors, and integration intermediates formed by site-specific recombinases and transposases associated with CRISPR systems. Here, we discuss the cellular processes that repair CRISPR-generated DNA lesions and describe strategies to obtain desirable genomic changes through modulation of DNA repair pathways. Advances in our understanding of the DNA repair circuitry, in conjunction with the rapid development of innovative genome editing technologies, promise to greatly enhance our ability to improve food production, combat environmental pollution, develop cell-based therapies, and cure genetic and infectious diseases.

OUP accepted manuscript

CRISPR/Cas base editors promise nucleotide-level control over DNA sequences, but the determinants of their activity remain incompletely understood. We measured base editing frequencies in two human cell lines for two cytosine and two adenine base editors at ∼14 000 target sequences and find that base editing activity is sequence-biased, with largest effects from nucleotides flanking the target base. Whether a base is edited depends strongly on the combination of its position in the target and the preceding base, acting to widen or narrow the effective editing window. The impact of features on editing rate depends on the position, with sequence bias efficacy mainly influencing bases away from the center of the window. We use these observations to train a machine learning model to predict editing activity per position, with accuracy ranging from 0.49 to 0.72 between editors, and with better generalization across datasets than existing tools. We demonstrate the usefulness of our model by predicting the efficacy of disease mutation correcting guides, and find that most of them suffer from more unwanted editing than pure outcomes. This work unravels the position-specificity of base editing biases and allows more efficient planning of editing campaigns in experimental and therapeutic contexts.

Mechanistic insights into genetic susceptibility to prostate cancer

Prostate cancer (PCa) is the second most common cancer in men and is a highly heritable disease that affects millions of individuals worldwide. Genome-wide association studies have to date discovered nearly 270 genetic loci harboring hundreds of single nucleotide polymorphisms (SNPs) that are associated with PCa susceptibility. In contrast, the functional characterization of the mechanisms underlying PCa risk association is still growing. Given that PCa risk-associated SNPs are highly enriched in noncoding cis-regulatory genomic regions, accumulating evidence suggests a widespread modulation of transcription factor chromatin binding and allelic enhancer activity by these noncoding SNPs, thereby dysregulating gene expression. Emerging studies have shown that a proportion of noncoding variants can modulate the formation of transcription factor complexes at enhancers and CTCF-mediated 3D genome architecture. Interestingly, DNA methylation-regulated CTCF binding could orchestrate a long-range chromatin interaction between PCa risk enhancer and causative genes. Additionally, one-causal-variant-two-risk genes or multiple-risk-variant-multiple-genes are prevalent in some PCa risk-associated loci. In this review, we will discuss the current understanding of the general principles of SNP-mediated gene regulation, experimental advances, and functional evidence supporting the mechanistic roles of several PCa genetic loci with potential clinical impact on disease prevention and treatment.

Towards a CRISPeR understanding of homologous recombination with high-throughput functional genomics

CRISPR-dependent genome editing enables the study of genes and mutations on a large scale. Here we review CRISPR-based functional genomics technologies that generate gene knockouts and single nucleotide variants (SNVs) and discuss how their use has provided new important insights into the function of homologous recombination (HR) genes. In particular, we highlight discoveries from CRISPR screens that have contributed to define the response to PARP inhibition in cells deficient for the HR genes BRCA1 and BRCA2, uncover genes whose loss causes synthetic lethality in combination with BRCA1/2 deficiency, and characterize the function of BRCA1/2 SNVs of uncertain clinical significance. Further use of these approaches, combined with next-generation CRISPR-based technologies, will aid to dissect the genetic network of the HR pathway, define the impact of HR mutations on cancer etiology and treatment, and develop novel targeted therapies for HR-deficient tumors.

Comprehensive Identification of Deleterious TP53 Missense VUS Variants Based on Their Impact on TP53 Structural Stability

TP53 plays critical roles in maintaining genome stability. Deleterious genetic variants damage the function of TP53, causing genome instability and increased cancer risk. Of the large quantity of genetic variants identified in TP53, however, many remain functionally unclassified as variants of unknown significance (VUS) due to the lack of evidence. This is reflected by the presence of 749 (42%) VUS of the 1785 germline variants collected in the ClinVar database. In this study, we addressed the deleteriousness of TP53 missense VUS. Utilizing the protein structure-based Ramachandran Plot-Molecular Dynamics Simulation (RPMDS) method that we developed, we measured the effects of missense VUS on TP53 structural stability. Of the 340 missense VUS tested, we observed deleterious evidence for 193 VUS, as reflected by the TP53 structural changes caused by the VUS-substituted residues. We compared the results from RPMDS with those from other in silico methods and observed higher specificity of RPMDS in classification of TP53 missense VUS than these methods. Data from our current study address a long-standing challenge in classifying the missense VUS in TP53, one of the most important tumor suppressor genes.

Identification of pathogenic variants in cancer genes using base editing screens with editing efficiency correction

BACKGROUND

Millions of nucleotide variants are identified through cancer genome sequencing and it is clinically important to identify the pathogenic variants among them. By introducing base substitutions at guide RNA target regions in the genome, CRISPR-Cas9-based base editors provide the possibility for evaluating a large number of variants in their genomic context. However, the variability in editing efficiency and the complexity of outcome mapping are two existing problems for assigning guide RNA effects to variants in base editing screens.

RESULTS

To improve the identification of pathogenic variants, we develop a framework to combine base editing screens with sgRNA efficiency and outcome mapping. We apply the method to evaluate more than 9000 variants across all the exons of BRCA1 and BRCA2 genes. Our efficiency-corrected scoring model identifies 910 loss-of-function variants for BRCA1/2, including 151 variants in the noncoding part of the genes such as the 5' untranslated regions. Many of them are identified in cancer patients and are reported as "benign/likely benign" or "variants of uncertain significance" by clinicians. Our data suggest a need to re-evaluate their clinical significance, which may be helpful for risk assessment and treatment of breast and ovarian cancer.

CONCLUSIONS

Our results suggest that base editing screens with efficiency correction is a powerful strategy to identify pathogenic variants in a high-throughput manner. Applying this strategy to assess variants in both coding and noncoding regions of the genome could have a direct impact on the interpretation of cancer variants.