Accurate classification of BRCA1 variants with saturation genome editing

Structural and functional prediction, evaluation, and validation in the post-sequencing era

The surge of genome sequencing data has underlined substantial genetic variants of uncertain significance (VUS). The decryption of VUS discovered by sequencing poses a major challenge in the post-sequencing era. Although experimental assays have progressed in classifying VUS, only a tiny fraction of the human genes have been explored experimentally. Thus, it is urgently needed to generate state-of-the-art functional predictors of VUS in silico. Artificial intelligence (AI) is an invaluable tool to assist in the identification of VUS with high efficiency and accuracy. An increasing number of studies indicate that AI has brought an exciting acceleration in the interpretation of VUS, and our group has already used AI to develop protein structure-based prediction models. In this review, we provide an overview of the previous research on AI-based prediction of missense variants, and elucidate the challenges and opportunities for protein structure-based variant prediction in the post-sequencing era.

Using AI-predicted protein structures as a reference to predict loss-of-function activity in tumor suppressor breast cancer genes

Background

The loss-of-function (LOF) classification of most missense variants in tumor suppressor breast cancer genes BRCA1, BRCA2, PALB2, and RAD51C remains unclassified and confounds clinical actionability. Classifying these variants is challenging due to their rarity, leading clinicians to rely on in silico predictive methods. Protein stability changes are associated with function, making stability predictors valuable. Stability predictions upon missense variant perturbations require high-resolution protein structures. However, the availability of these high-resolution structures is lacking. This study explores using generative AI to predict high-resolution protein structures, which can then be analyzed with in silico protein stability prediction methods to assess LOF activity in ordered regions of the protein. This study also determines the appropriate in silico protein stability and dedicated in silico missense prediction methods in dbNSFP v4.7 database to predict LOF activity in ordered regions of these four genes. Functional classifications from homology recombination DNA repair (HDR) assays and variant classifications from the ClinVar database provide a reliable dataset for evaluating the performance of these in silico prediction methods.

Results

Complex AlphaFold2 structures of the BRCA1-C terminal (BRCT) domain and the DNA-binding (DB) domain of BRCA2, analyzed using protein stability tool FoldX predicts LOF activity from missense variants significantly better than experimentally-derived structures in ordered regions. The BRCT domain achieved an Area Under the Curve (AUC)= 0.861 (95 % CI:0.858-0.863) and AUC= 0.842 (95 % CI:0.840-0.845), while the DB domain achieved an AUC= 0.836 (95 % CI:0.8322-0.841), compared to AUC= 0.847 (95 % CI:0.844-0.850) and AUC= 0.835 (95 % CI:0.832-0.837) from the BRCT domain, and AUC= 0.830 (95 % CI:0.821-0.8320) from the DB domain from experimentally-derived structures. Protein stability does not predict LOF activity from missense variants better than dedicated in silico missense predictors. Overall, we find that AlphaMissense ranks highly, with an average AUC= 0.890 (95 % CI 0.886-0.895) from ordered regions across these four cancer genes, compared to all other in silico missense predictors present in the dbNSFP database.

Conclusions

The study reveals that generative AI protein predicted structures can outperform experimentally-derived structures in evaluating LOF activity from predicted protein stability in ordered regions of genes BRCA1, BRCA2, PALB2 and RAD51C. The study also highlights the predictive performance of AlphaMissense as the premier in silico missense prediction method to predict LOF activity from missense variants in these four tumor suppressor breast cancer genes. The code for this study can be downloaded for free on GitHub (https://github.com/rohandavidg/CarePred).

ResMAP-a saturation mutagenesis platform enabling parallel profiling of target-specific resistance-conferring mutations in Plasmodium

New and improved drugs are required for the treatment and ultimate eradication of malaria. The efficacy of front-line therapies is now threatened by emerging drug resistance; thus, new tools to support the development of drugs with a lower propensity for resistance are needed. Here, we describe the development of a RESistance Mapping And Profiling (ResMAP) platform for the identification of resistance-conferring mutations in Plasmodium drug targets. Proof-of-concept studies focused on interrogating the antimalarial drug target, Plasmodium falciparum lysyl tRNA synthetase (PfKRS). Saturation mutagenesis was used to construct a plasmid library encoding all conceivable mutations within a 20-residue span at the base of the PfKRS ATP-binding pocket. The superior transfection efficiency of Plasmodium knowlesi was exploited to generate a high coverage parasite library expressing PfKRS bearing all possible amino acid changes within this region of the enzyme. The selection of the library with PfKRS inhibitors, cladosporin and DDD01510706, successfully identified multiple resistance-conferring substitutions. Genetic validation of a subset of these mutations confirmed their direct role in resistance, with computational modeling used to dissect the structural basis of resistance. The application of ResMAP to inform the development of resistance-resilient antimalarials of the future is discussed.

IMPORTANCE

An increase in treatment failures for malaria highlights an urgent need for new tools to understand and minimize the spread of drug resistance. We describe the development of a RESistance Mapping And Profiling (ResMAP) platform for the identification of resistance-conferring mutations in Plasmodium spp, the causative agent of malaria. Saturation mutagenesis was used to generate a mutation library containing all conceivable mutations for a region of the antimalarial-binding site of a promising drug target, Plasmodium falciparum lysyl tRNA synthetase (PfKRS). Screening of this high-coverage library with characterized PfKRS inhibitors revealed multiple resistance-conferring substitutions including several clinically relevant mutations. Genetic validation of these mutations confirmed resistance of up to 100-fold and computational modeling dissected their role in drug resistance. We discuss potential applications of this data including the potential to design compounds that can bypass the most serious resistance mutations and future resistance surveillance.

FUSE: Improving the estimation and imputation of variant impacts in functional screening

Deep mutational scanning enables high-throughput functional assessment of genetic variants. While phenotypic measurements from screening assays generally align with clinical outcomes, experimental noise may affect the accuracy of individual variant estimates. We developed the FUSE (functional substitution estimation) pipeline, which leverages measurements collectively within screening assays to improve the estimation of variant impacts. Drawing data from 115 published functional assays, FUSE assesses the mean functional effect per amino acid position and makes estimates for individual allelic variants. It enhances the correlation of variant functional effects from different assay platforms and increases the classification accuracy of missense variants in ClinVar across 29 genes (area under the receiver operating characteristic [ROC] curve [AUC] from 0.83 to 0.90). In UK Biobank patients with rare missense variants in BRCA1, LDLR, or TP53, FUSE improves the classification accuracy of associated phenotypes. FUSE can also impute variant effects for substitutions not experimentally screened. This approach improves accuracy and broadens the utility of data from functional screening.

High-resolution functional mapping of RAD51C by saturation genome editing

Pathogenic variants in RAD51C confer an elevated risk of breast and ovarian cancer, while individuals homozygous for specific RAD51C alleles may develop Fanconi anemia. Using saturation genome editing (SGE), we functionally assess 9,188 unique variants, including >99.5% of all possible coding sequence single-nucleotide alterations. By computing changes in variant abundance and Gaussian mixture modeling (GMM), we functionally classify 3,094 variants to be disruptive and use clinical truth sets to reveal an accuracy/concordance of variant classification >99.9%. Cell fitness was the primary assay readout allowing us to observe a phenomenon where specific missense variants exhibit distinct depletion kinetics potentially suggesting that they represent hypomorphic alleles. We further explored our exhaustive functional map, revealing critical residues on the RAD51C structure and resolving variants found in cancer-segregating kindred. Furthermore, through interrogation of UK Biobank and a large multi-center ovarian cancer cohort, we find significant associations between SGE-depleted variants and cancer diagnoses.

Pancreatic Cancer Screening among High-risk Individuals

Pancreatic ductal adenocarcinoma (PDAC) continues to remain one of the leading causes of cancer-related death. Unlike other malignancies where universal screening is recommended, the same cannot be said for PDAC. The purpose of this study is to review which patients are at high risk of developing PDAC and therefore candidates for screening, methods/frequency of screening, and risk for these groups of patients.

Federated analysis of autosomal recessive coding variants in 29,745 developmental disorder patients from diverse populations

Autosomal recessive coding variants are well-known causes of rare disorders. We quantified the contribution of these variants to developmental disorders in a large, ancestrally diverse cohort comprising 29,745 trios, of whom 20.4% had genetically inferred non-European ancestries. The estimated fraction of patients attributable to exome-wide autosomal recessive coding variants ranged from ~2-19% across genetically inferred ancestry groups and was significantly correlated with average autozygosity. Established autosomal recessive developmental disorder-associated (ARDD) genes explained 84.0% of the total autosomal recessive coding burden, and 34.4% of the burden in these established genes was explained by variants not already reported as pathogenic in ClinVar. Statistical analyses identified two novel ARDD genes: KBTBD2 and ZDHHC16. This study expands our understanding of the genetic architecture of developmental disorders across diverse genetically inferred ancestry groups and suggests that improving strategies for interpreting missense variants in known ARDD genes may help diagnose more patients than discovering the remaining genes.

Leaving no patient behind! Expert recommendation in the use of innovative technologies for diagnosing rare diseases

Genetic diagnosis plays a crucial role in rare diseases, particularly with the increasing availability of emerging and accessible treatments. The International Rare Diseases Research Consortium (IRDiRC) has set its primary goal as: "Ensuring that all patients who present with a suspected rare disease receive a diagnosis within one year if their disorder is documented in the medical literature". Despite significant advances in genomic sequencing technologies, more than half of the patients with suspected Mendelian disorders remain undiagnosed. In response, IRDiRC proposes the establishment of "a globally coordinated diagnostic and research pipeline". To help facilitate this, IRDiRC formed the Task Force on Integrating New Technologies for Rare Disease Diagnosis. This multi-stakeholder Task Force aims to provide an overview of the current state of innovative diagnostic technologies for clinicians and researchers, focusing on the patient's diagnostic journey. Herein, we provide an overview of a broad spectrum of emerging diagnostic technologies involving genomics, epigenomics and multi-omics, functional testing and model systems, data sharing, bioinformatics, and Artificial Intelligence (AI), highlighting their advantages, limitations, and the current state of clinical adaption. We provide expert recommendations outlining the stepwise application of these innovative technologies in the diagnostic pathways while considering global differences in accessibility. The importance of FAIR (Findability, Accessibility, Interoperability, and Reusability) and CARE (Collective benefit, Authority to control, Responsibility, and Ethics) data management is emphasized, along with the need for enhanced and continuing education in medical genomics. We provide a perspective on future technological developments in genome diagnostics and their integration into clinical practice. Lastly, we summarize the challenges related to genomic diversity and accessibility, highlighting the significance of innovative diagnostic technologies, global collaboration, and equitable access to diagnosis and treatment for people living with rare disease.

Variability in proliferative and migratory defects in Hirschsprung disease-associated RET pathogenic variants

Despite the extensive genetic heterogeneity of Hirschsprung disease (HSCR; congenital colonic aganglionosis) 72% of patients harbor pathogenic variants in 10 genes that form a gene regulatory network (GRN) controlling the development of the enteric nervous system (ENS). Among these genes, the receptor tyrosine kinase gene RET is the most significant contributor, accounting for pathogenic variants in 12%-50% of patients depending on phenotype. RET plays a critical role in the proliferation and migration of ENS precursors, and defects in these processes lead to HSCR. However, despite the gene's importance in HSCR, the functional consequences of RET pathogenic variants and their mechanism of disease remain poorly understood. To address this, we investigated the proliferative and migratory phenotypes in a RET-dependent neural crest-derived cell line harboring one of five missense (L56M, E178Q, Y791F, S922Y, F998L) or three nonsense (Y204X, R770X, Y981X) pathogenic heterozygous variants. Using a combination of cDNA-based and CRISPR-based PRIME editing coupled with quantitative proliferation and migration assays, we detected significant losses in cell proliferation and migration in three missense (E178Q, S922Y, F998L) and all nonsense variants. Our data suggests that the Y791F variant, whose pathogenicity has been debated, is likely not pathogenic. Importantly, the severity of migration loss did not consistently correlate with proliferation defects, and the phenotypic severity of nonsense variants was independent of their position within the RET protein. This study highlights the necessity and feasibility of targeted functional assays to accurately assess the pathogenicity of HSCR-associated variants, rather than relying solely on machine learning predictions, which could themselves be refined by incorporating such functional data.

Genomic Language Models: Opportunities and Challenges

Large language models (LLMs) are having transformative impacts across a wide range of scientific fields, particularly in the biomedical sciences. Just as the goal of Natural Language Processing is to understand sequences of words, a major objective in biology is to understand biological sequences. Genomic Language Models (gLMs), which are LLMs trained on DNA sequences, have the potential to significantly advance our understanding of genomes and how DNA elements at various scales interact to give rise to complex functions. To showcase this potential, we highlight key applications of gLMs, including functional constraint prediction, sequence design, and transfer learning. Despite notable recent progress, however, developing effective and efficient gLMs presents numerous challenges, especially for species with large, complex genomes. Here, we discuss major considerations for developing and evaluating gLMs.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - in the onset and phenotypic presentation of rare diseases. Here, we comprehensively map individual polygenic liability for 1102 open-source PGS in a cohort of 3059 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. Using this resource, we demonstrate increased polygenic liability in probands with an inherited candidate disease variant (VUS) compared to unaffected carrier parents. Further, we show an enrichment for large-effect rare variants in putative core PGS genes for associated complex traits. Overall, our study supports and expands on previous findings of complex trait associations in rare diseases, implicates polygenic liability as a potential mechanism underlying variable penetrance of candidate causal variants, and provides a framework for identifying novel candidate rare disease genes.

Assigning credit where it is due: an information content score to capture the clinical value of multiplexed assays of variant effect

BACKGROUND

A variant can be pathogenic or benign with relation to a human disease. Current classification categories from benign to pathogenic reflect a probabilistic summary of the current understanding. A primary metric of clinical utility for multiplexed assays of variant effect (MAVE) is the number of variants that can be reclassified from uncertain significance (VUS). However, a gap in this measure of utility is that it underrepresents the information gained from MAVEs. The aim of this study was to develop an improved quantification metric for MAVE utility. We propose adopting an information content approach that includes data that does not reclassify variants will better reflect true information gain. We adopted an information content approach to evaluate the information gain, in bits, for MAVEs of BRCA1, PTEN, and TP53. Here, one bit represents the amount of information required to completely classify a single variant starting from no information.

RESULTS

BRCA1 MAVEs produced a total of 831.2 bits of information, 6.58% of the total missense information in BRCA1 and a 22-fold increase over the information that only contributed to VUS reclassification. PTEN MAVEs produced 2059.6 bits of information which represents 32.8% of the total missense information in PTEN and an 85-fold increase over the information that contributed to VUS reclassification. TP53 MAVEs produced 277.8 bits of information which represents 6.22% of the total missense information in TP53 and a 3.5-fold increase over the information that contributed to VUS reclassification.

CONCLUSIONS

An information content approach will more accurately portray information gained through MAVE mapping efforts than by counting the number of variants reclassified. This information content approach may also help define the impact of guideline changes that modify the information definitions used to classify groups of variants.

Co-observation of germline pathogenic variants in breast cancer predisposition genes: Results from analysis of the BRIDGES sequencing dataset

Co-observation of a gene variant with a pathogenic variant in another gene that explains the disease presentation has been designated as evidence against pathogenicity for commonly used variant classification guidelines. Multiple variant curation expert panels have specified, from consensus opinion, that this evidence type is not applicable for the classification of breast cancer predisposition gene variants. Statistical analysis of sequence data for 55,815 individuals diagnosed with breast cancer from the BRIDGES sequencing project was undertaken to formally assess the utility of co-observation data for germline variant classification. Our analysis included expected loss-of-function variants in 11 breast cancer predisposition genes and pathogenic missense variants in BRCA1, BRCA2, and TP53. We assessed whether co-observation of pathogenic variants in two different genes occurred more or less often than expected under the assumption of independence. Co-observation of pathogenic variants in each of BRCA1, BRCA2, and PALB2 with the remaining genes was less frequent than expected. This evidence for depletion remained after adjustment for age at diagnosis, study design (familial versus population-based), and country. Co-observation of a variant of uncertain significance in BRCA1, BRCA2, or PALB2 with a pathogenic variant in another breast cancer gene equated to supporting evidence against pathogenicity following criterion strength assignment based on the likelihood ratio and showed utility in reclassification of missense BRCA1 and BRCA2 variants identified in BRIDGES. Our approach has applicability for assessing the value of co-observation as a predictor of variant pathogenicity in other clinical contexts, including for gene-specific guidelines developed by ClinGen Variant Curation Expert Panels.

Analysis of more than 400,000 women provides case-control evidence for BRCA1 and BRCA2 variant classification

Clinical genetic testing identifies variants causal for hereditary cancer, information that is used for risk assessment and clinical management. Unfortunately, some variants identified are of uncertain clinical significance (VUS), complicating patient management. Case-control data is one evidence type used to classify VUS, and previous findings indicate that case-control likelihood ratios (LRs) outperform odds ratios for variant classification. As an initiative of the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) Analytical Working Group we analyzed germline sequencing data of BRCA1 and BRCA2 from 96,691 female breast cancer cases and 303,925 unaffected controls from three studies: the BRIDGES study of the Breast Cancer Association Consortium, the Cancer Risk Estimates Related to Susceptibility consortium, and the UK Biobank. We observed 11,227 BRCA1 and BRCA2 variants, with 6,921 being coding, covering 23.4% of BRCA1 and BRCA2 VUS in ClinVar and 19.2% of ClinVar curated (likely) benign or pathogenic variants. Case-control LR evidence was highly consistent with ClinVar assertions for (likely) benign or pathogenic variants; exhibiting 99.1% sensitivity and 95.4% specificity for BRCA1 and 92.2% sensitivity and 86.6% specificity for BRCA2. This approach provides case-control evidence for 785 unclassified variants, that can serve as a valuable element for clinical classification.

Venous thromboembolic disease genetics: from variants to function

Venous thromboembolic disease (VTE) is a prevalent and potentially life-threatening vascular disease, including both deep vein thrombosis and pulmonary embolism. This review will focus on recent insights into the heritable factors that influence an individual's risk for VTE. Here, we will explore not only the discovery of new genetic risk variants but also the importance of functional characterization of these variants. These genome-wide studies should lead to a better understanding of the biological role of genes inside and outside of the canonical coagulation system in thrombus formation and lead to an improved ability to predict an individual's risk of VTE. Further understanding of the molecular mechanisms altered by genetic variation in VTE risk will be accelerated by further human genome sequencing efforts and the use of functional genetic screens.

Deciphering the impact of genomic variation on function

Our genomes influence nearly every aspect of human biology-from molecular and cellular functions to phenotypes in health and disease. Studying the differences in DNA sequence between individuals (genomic variation) could reveal previously unknown mechanisms of human biology, uncover the basis of genetic predispositions to diseases, and guide the development of new diagnostic tools and therapeutic agents. Yet, understanding how genomic variation alters genome function to influence phenotype has proved challenging. To unlock these insights, we need a systematic and comprehensive catalogue of genome function and the molecular and cellular effects of genomic variants. Towards this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations and predictive modelling to investigate the relationships among genomic variation, genome function and phenotypes. IGVF will create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how such effects connect through gene-regulatory and protein-interaction networks. These experimental data, computational predictions and accompanying standards and pipelines will be integrated into an open resource that will catalyse community efforts to explore how our genomes influence biology and disease across populations.

Custom target-sequencing in triple-negative and luminal breast cancer from young Brazilian patients

OBJECTIVES

To identify somatic mutations in tumors from young women with triple-negative or luminal breast cancer, through targeted sequencing and to explore the cancer driver potential of these gene variants.

METHODS

A customized gene panel was assembled based on data from previous sequencing studies of breast cancer from young women. Triple-negative and luminal tumors and paired blood samples from young breast cancer patients were sequenced, and identified gene variants were searched for their driver potential, in databases and literature. Additionally, the authors performed an exploratory analysis using large, curated databases to evaluate the frequency of somatic mutations in this gene panel in tumors stratified by age groups (every 10 years).

RESULTS

A total of 28 young women had their tumoral tissue and blood samples sequenced. Using a customized panel of 64 genes, the authors could detect cancer drivers in 11/12 (91.7 %) TNBC samples and 11/16 (68.7 %) luminal samples. Among TNBC patients, the most frequent cancer driver was TP53, followed by NF1, NOTCH1 and PTPN13. In luminal samples, PIK3CA and GATA3 were the main cancer drivers, and other drivers were GRHL2 and SMURF2. CACNA1E was involved in both TN and luminal BC. The exploratory analysis also indicated a role for SMURF2 in luminal BC development in young patients.

CONCLUSIONS

The data further indicates that some cancer drivers are more common in a specific breast cancer subtype from young patients, such as TP53 in TNBC and PIK3CA and GATA3 in luminal samples. These results also provide additional evidence that some genes not considered classical cancer-causing genes, such as CACNA1E, GRHL2 and SMURF2 might be cancer drivers in this age group.

LAT encodes T cell activation pathway balance

Immune cells transduce environmental stimuli into responses essential for host health via complex signaling cascades. T cells, in particular, leverage their unique T cell receptors (TCRs) to detect specific Human Leukocyte Antigen (HLA)-presented peptides. TCR activation is then relayed via linker for activation of T cells (LAT), a TCR-proximal disordered adapter protein, which organizes protein partners and mediates the propagation of signals down diverse pathways including NFAT and AP-1. Here, we studied how balanced downstream pathway activation is encoded in the amino acid sequence of LAT. To comprehensively profile the sequence-function relationship of LAT, we developed a pooled, single-cell, high-content screening approach in which a large series of mutants in the LAT protein were analyzed to characterize their effects on T cell activation. Measuring epigenetic, transcriptomic, and cell surface protein dynamics of single cells harboring distinct LAT mutants, we found functional regions spanning over 40% of the LAT amino acid sequence. Conserved sequence motifs for protein interactions along with charge distribution are critical sequence features, and contribute to interpretation of human genetic variation in LAT. While mutant defect severity spans from moderate to complete loss of function, nearly all defective mutants, irrespective of their position in LAT, confer balanced defects across all downstream pathways. To understand the molecular basis for this observation, we performed proximal protein labeling which demonstrated that disruption of LAT interaction with a single partner protein indirectly disrupts other partner interactions, likely through the dual roles of these proteins as effectors of downstream pathways and bridging factors between LAT molecules. Overall, we report widely distributed functional regions throughout a disordered adapter and a precise physical organization of LAT and interacting molecules which constrains signaling outputs. More broadly, we describe an approach for interrogating sequence-function relationships for proteins with complex activities across regulatory layers of the cell.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Concurrent RB1 Loss and BRCA Deficiency Predicts Enhanced Immunologic Response and Long-term Survival in Tubo-ovarian High-grade Serous Carcinoma

PURPOSE

The purpose of this study was to evaluate RB1 expression and survival across ovarian carcinoma histotypes and how co-occurrence of BRCA1 or BRCA2 (BRCA) alterations and RB1 loss influences survival in tubo-ovarian high-grade serous carcinoma (HGSC).

EXPERIMENTAL DESIGN

RB1 protein expression was classified by immunohistochemistry in ovarian carcinomas of 7,436 patients from the Ovarian Tumor Tissue Analysis consortium. We examined RB1 expression and germline BRCA status in a subset of 1,134 HGSC, and related genotype to overall survival (OS), tumor-infiltrating CD8+ lymphocytes, and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cells with and without BRCA1 alterations to model co-loss with treatment response. We performed whole-genome and transcriptome data analyses on 126 patients with primary HGSC to characterize tumors with concurrent BRCA deficiency and RB1 loss.

RESULTS

RB1 loss was associated with longer OS in HGSC but with poorer prognosis in endometrioid ovarian carcinoma. Patients with HGSC harboring both RB1 loss and pathogenic germline BRCA variants had superior OS compared with patients with either alteration alone, and their median OS was three times longer than those without pathogenic BRCA variants and retained RB1 expression (9.3 vs. 3.1 years). Enhanced sensitivity to cisplatin and paclitaxel was seen in BRCA1-altered cells with RB1 knockout. Combined RB1 loss and BRCA deficiency correlated with transcriptional markers of enhanced IFN response, cell-cycle deregulation, and reduced epithelial-mesenchymal transition. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1.

CONCLUSIONS

Co-occurrence of RB1 loss and BRCA deficiency was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.

Understanding variants of unknown significance: the computational frontier

The rapid advancement of sequencing technologies has led to the identification of numerous mutations in cancer genomes, many of which are variants of unknown significance (VUS). Computational models are increasingly being used to predict the functional impact of these mutations, in both coding and noncoding regions. Integration of these models with emerging genomic datasets will refine our understanding of mutation effects and guide clinical decision making. Future advancements in modeling protein interactions and transcriptional regulation will further enhance our ability to interpret VUS. Periodic incorporation of these developments into VUS reclassification practice has the potential to significantly improve personalized cancer care.

Germline Genetic Testing for Hereditary Breast and Ovarian Cancer: Current Concepts in Risk Evaluation

Our understanding of hereditary breast and ovarian cancer has significantly improved over the past two decades. In addition to BRCA1/2, pathogenic variants in several other DNA-repair genes have been shown to increase the risks of breast and ovarian cancer. The magnitude of cancer risk is impacted not only by the gene involved, but also by family history of cancer, polygenic risk scores, and, in certain genes, pathogenic variant type or location. While estimates of breast and ovarian cancer risk associated with pathogenic variants are available, these are predominantly based on studies of high-risk populations with young age at diagnosis of cancer, multiple primary cancers, or family history of cancer. More recently, breast cancer risk for germline pathogenic variant carriers has been estimated from population-based studies. Here, we provide a review of the field of germline genetic testing and risk evaluation for hereditary breast and ovarian cancers in high-risk and population-based settings.

Scalable approaches for generating, validating and incorporating data from high-throughput functional assays to improve clinical variant classification

As the adoption and scope of genetic testing continue to expand, interpreting the clinical significance of DNA sequence variants at scale remains a formidable challenge, with a high proportion classified as variants of uncertain significance (VUSs). Genetic testing laboratories have historically relied, in part, on functional data from academic literature to support variant classification. High-throughput functional assays or multiplex assays of variant effect (MAVEs), designed to assess the effects of DNA variants on protein stability and function, represent an important and increasingly available source of evidence for variant classification, but their potential is just beginning to be realized in clinical lab settings. Here, we describe a framework for generating, validating and incorporating data from MAVEs into a semi-quantitative variant classification method applied to clinical genetic testing. Using single-cell gene expression measurements, cellular evidence models were built to assess the effects of DNA variation in 44 genes of clinical interest. This framework was also applied to models for an additional 22 genes with previously published MAVE datasets. In total, modeling data was incorporated from 24 genes into our variant classification method. These data contributed evidence for classifying 4043 observed variants in over 57,000 individuals. Genetic testing laboratories are uniquely positioned to generate, analyze, validate, and incorporate evidence from high-throughput functional data and ultimately enable the use of these data to provide definitive clinical variant classifications for more patients.

Population Diversity at the Single-Cell Level

Population-scale single-cell genomics is a transformative approach for unraveling the intricate links between genetic and cellular variation. This approach is facilitated by cutting-edge experimental methodologies, including the development of high-throughput single-cell multiomics and advances in multiplexed environmental and genetic perturbations. Examining the effects of natural or synthetic genetic variants across cellular contexts provides insights into the mutual influence of genetics and the environment in shaping cellular heterogeneity. The development of computational methodologies further enables detailed quantitative analysis of molecular variation, offering an opportunity to examine the respective roles of stochastic, intercellular, and interindividual variation. Future opportunities lie in leveraging long-read sequencing, refining disease-relevant cellular models, and embracing predictive and generative machine learning models. These advancements hold the potential for a deeper understanding of the genetic architecture of human molecular traits, which in turn has important implications for understanding the genetic causes of human disease.

Large-scale discovery of chromatin dysregulation induced by oncofusions and other protein-coding variants

Population-scale databases have expanded to millions of protein-coding variants, yet insight into their mechanistic consequences has lagged. Here we present PROD-ATAC, a high-throughput method for discovering the effects of protein-coding variants on chromatin regulation. A pooled variant library is expressed in a disease-agnostic cell line, and single-cell assay for transposase-accessible chromatin resolves each variant's effect on the chromatin landscape. Using PROD-ATAC, we characterized the effects of more than 100 oncofusions (cancer-causing chimeric proteins) and controls and revealed that chromatin remodeling is common to fusions spanning an enormous range of fusion frequencies. Furthermore, fusion-induced dysregulation can be context agnostic, as observed mechanisms often overlapped with cancer and cell-type-specific prior knowledge. We also showed that gain-of-function activity is common among oncofusions. This work begins to outline a global map of fusion-induced chromatin alterations. We suggest that there might be convergent mechanisms among disparate oncofusions and shared modes of dysregulation among fusions present in tumors at different frequencies. PROD-ATAC is generalizable to any set of protein-coding variants.

In Silico Prediction of BRCA1 and BRCA2 Variants with Conflicting Clinical Interpretation in a Cohort of Breast Cancer Patients

Germline BRCA1/2 alteration has been linked to an increased risk of hereditary breast and ovarian cancer syndromes. As a result, genetic testing, based on NGS, allows us to identify a high number of variants of uncertain significance (VUS) or conflicting interpretation of pathogenicity (CIP) variants. The identification of CIP/VUS is often considered inconclusive and clinically not actionable for the patients' and unaffected carriers' management. In this context, their assessment and classification remain a significant challenge. The aim of the study was to investigate whether the in silico prediction tools (PolyPhen-2, SIFT, Mutation Taster and PROVEAN) could predict the potential clinical impact and significance of BRCA1/2 CIP/VUS alterations, eventually impacting the clinical management of Breast Cancer subjects. In a cohort of 860 BC patients, 10.6% harbored BRCA1 or BRCA2 CIP/VUS alterations, mostly observed in BRCA2 sequences (85%). Among them, forty-two out of fifty-five alterations were predicted as damaging, with at least one in silico that used tools. Prediction agreement of the four tools was achieved in 45.5% of patients. Moreover, the highest consensus was obtained in twelve out of forty-two (28.6%) mutations by considering three out of four in silico algorithms. The use of prediction tools may help to identify variants with a potentially damaging effect. The lack of substantial agreement between the different algorithms suggests that the bioinformatic approaches should be combined with the personal and family history of the cancer patients.

Comprehensive deletion scan of anti-CRISPR AcrIIA4 reveals essential and dispensable domains for Cas9 inhibition

Delineating a protein's essential and dispensable domains provides critical insight into how it carries out its function. Here, we developed a high-throughput method to synthesize and test the functionality of all possible in-frame and continuous deletions in a gene of interest, enabling rapid and unbiased determination of protein domain importance. Our approach generates precise deletions using a CRISPR library framework that is free from constraints of gRNA target site availability and efficacy. We applied our method to AcrIIA4, a phage-encoded anti-CRISPR protein that robustly inhibits SpCas9. Extensive structural characterization has shown that AcrIIA4 physically occupies the DNA-binding interfaces of several SpCas9 domains; nonetheless, the importance of each AcrIIA4 interaction for SpCas9 inhibition is unknown. We used our approach to determine the essential and dispensable regions of AcrIIA4. Surprisingly, not all contacts with SpCas9 were required, and in particular, we found that the AcrIIA4 loop that inserts into SpCas9's RuvC catalytic domain can be deleted. Our results show that AcrIIA4 inhibits SpCas9 primarily by blocking PAM binding, and that its interaction with the SpCas9 catalytic domain is inessential.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Saturation genome editing maps the functional spectrum of pathogenic VHL alleles

To maximize the impact of precision medicine approaches, it is critical to identify genetic variants underlying disease and to accurately quantify their functional effects. A gene exemplifying the challenge of variant interpretation is the von Hippel-Lindautumor suppressor (VHL). VHL encodes an E3 ubiquitin ligase that regulates the cellular response to hypoxia. Germline pathogenic variants in VHL predispose patients to tumors including clear cell renal cell carcinoma (ccRCC) and pheochromocytoma, and somatic VHL mutations are frequently observed in sporadic renal cancer. Here we optimize and apply saturation genome editing to assay nearly all possible single-nucleotide variants (SNVs) across VHL's coding sequence. To delineate mechanisms, we quantify mRNA dosage effects and compare functional effects in isogenic cell lines. Function scores for 2,268 VHL SNVs identify a core set of pathogenic alleles driving ccRCC with perfect accuracy, inform differential risk across tumor types and reveal new mechanisms by which variants impact function. These results have immediate utility for classifying VHL variants encountered clinically and illustrate how precise functional measurements can resolve pleiotropic and dosage-dependent genotype-phenotype relationships across complete genes.

Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants

Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.

Saturation genome editing of BAP1 functionally classifies somatic and germline variants

Many variants that we inherit from our parents or acquire de novo or somatically are rare, limiting the precision with which we can associate them with disease. We performed exhaustive saturation genome editing (SGE) of BAP1, the disruption of which is linked to tumorigenesis and altered neurodevelopment. We experimentally characterized 18,108 unique variants, of which 6,196 were found to have abnormal functions, and then used these data to evaluate phenotypic associations in the UK Biobank. We also characterized variants in a large population-ascertained tumor collection, in cancer pedigrees and ClinVar, and explored the behavior of cancer-associated variants compared to that of variants linked to neurodevelopmental phenotypes. Our analyses demonstrated that disruptive germline BAP1 variants were significantly associated with higher circulating levels of the mitogen IGF-1, suggesting a possible pathological mechanism and therapeutic target. Furthermore, we built a variant classifier with >98% sensitivity and specificity and quantify evidence strengths to aid precision variant interpretation.

Discovering mechanisms of human genetic variation and controlling cell states at scale

Population-scale sequencing efforts have catalogued substantial genetic variation in humans such that variant discovery dramatically outpaces interpretation. We discuss how single-cell sequencing is poised to reveal genetic mechanisms at a rate that may soon approach that of variant discovery. The functional genomics toolkit is sufficiently modular to systematically profile almost any type of variation within increasingly diverse contexts and with molecularly comprehensive and unbiased readouts. As a result, we can construct deep phenotypic atlases of variant effects that span the entire regulatory cascade. The same conceptual approach to interpreting genetic variation should be applied to engineering therapeutic cell states. In this way, variant mechanism discovery and cell state engineering will become reciprocating and iterative processes towards genomic medicine.

Accelerated drug-resistant variant discovery with an enhanced, scalable mutagenic base editor platform

Personalized cancer therapeutics bring directed treatment options to patients based on their tumor's genetic signature. Unfortunately, tumor genomes are remarkably adaptable, and acquired resistance through gene mutation frequently occurs. Identifying mutations that promote resistance within drug-treated patient populations can be cost, resource, and time intensive. Accordingly, base editing, enabled by Cas9-deaminase domain fusions, has emerged as a promising approach for rapid, large-scale gene variant screening in situ. Here, we adapt and optimize a conditional activation-induced cytidine deaminase (AID)-dead Cas9 (dCas9) system, which demonstrates greater heterogeneity of edits with an expanded footprint compared to the most commonly utilized cytosine base editor, BE4. In combination with a custom single guide RNA (sgRNA) library, we identify individual and compound variants in epidermal growth factor receptor (EGFR) and v-raf murine sarcoma viral oncogene homolog B1 (BRAF) that confer resistance to established EGFR inhibitors. This system and analytical pipeline provide a simple, highly scalable platform for cis or trans drug-modifying variant discovery and for uncovering valuable insights into protein structure-function relationships.

In vitro evolution and whole genome analysis to study chemotherapy drug resistance in haploid human cells

In vitro evolution and whole genome analysis has proven to be a powerful method for studying the mechanism of action of small molecules in many haploid microbes but has generally not been applied to human cell lines in part because their diploid state complicates the identification of variants that confer drug resistance. To determine if haploid human cells could be used in MOA studies, we evolved resistance to five different anticancer drugs (doxorubicin, gemcitabine, etoposide, topotecan, and paclitaxel) using a near-haploid cell line (HAP1) and then analyzed the genomes of the drug resistant clones, developing a bioinformatic pipeline that involved filtering for high frequency alleles predicted to change protein sequence, or alleles which appeared in the same gene for multiple independent selections with the same compound. Applying the filter to sequences from 28 drug resistant clones identified a set of 21 genes which was strongly enriched for known resistance genes or known drug targets (TOP1, TOP2A, DCK, WDR33, SLCO3A1). In addition, some lines carried structural variants that encompassed additional known resistance genes (ABCB1, WWOX and RRM1). Gene expression knockdown and knockout experiments of 10 validation targets showed a high degree of specificity and accuracy in our calls and demonstrates that the same drug resistance mechanisms found in diverse clinical samples can be evolved, discovered and studied in an isogenic background.

A Novel Germline Mutation of BRCA1 and Integrated Analysis With Somatic Mutation in a Chinese Multi-Cancer Family

The presence of mutations in the BRCA1 gene (MIM: 113705) is widely recognized as a significant genetic predisposition for ovarian cancer. This study investigated the genomic mutations in a Chinese family with a history of ovarian, breast, and rectal adenocarcinoma. A novel germline mutation (Phe1695Val) in BRCA1 was identified through whole-exome sequencing. Subsequently, we performed whole-genome sequencing to identify somatic mutations and analyze mutational signatures in individuals carrying the novel germline mutation. Our findings revealed a correlation between somatic mutational signatures and the BRCA1 germline mutation in the proband with ovarian cancer, while no such association was observed in the tumor tissue from the patient with breast cancer. Furthermore, distinct somatic driver mutations were identified, a truncated mutation in the TP53 gene in the ovarian tumor tissue, and a hotspot mutation in the PIK3CA gene in the breast cancer. According to our findings, the BRCA1 F1695V mutation is linked to ovarian cancer susceptibility in the family and causes specific somatic mutational profiles.

High-throughput assays to assess variant effects on disease

Interpreting the wealth of rare genetic variants discovered in population-scale sequencing efforts and deciphering their associations with human health and disease present a critical challenge due to the lack of sufficient clinical case reports. One promising avenue to overcome this problem is deep mutational scanning (DMS), a method of introducing and evaluating large-scale genetic variants in model cell lines. DMS allows unbiased investigation of variants, including those that are not found in clinical reports, thus improving rare disease diagnostics. Currently, the main obstacle limiting the full potential of DMS is the availability of functional assays that are specific to disease mechanisms. Thus, we explore high-throughput functional methodologies suitable to examine broad disease mechanisms. We specifically focus on methods that do not require robotics or automation but instead use well-designed molecular tools to transform biological mechanisms into easily detectable signals, such as cell survival rate, fluorescence or drug resistance. Here, we aim to bridge the gap between disease-relevant assays and their integration into the DMS framework.

Structure-based network analysis predicts pathogenic variants in human proteins associated with inherited retinal disease

Advances in gene sequencing technologies have accelerated the identification of genetic variants, but better tools are needed to understand which are causal of disease. This would be particularly useful in fields where gene therapy is a potential therapeutic modality for a disease-causing variant such as inherited retinal disease (IRD). Here, we apply structure-based network analysis (SBNA), which has been successfully utilized to identify variant-constrained amino acid residues in viral proteins, to identify residues that may cause IRD if subject to missense mutation. SBNA is based entirely on structural first principles and is not fit to specific outcome data, which makes it distinct from other contemporary missense prediction tools. In 4 well-studied human disease-associated proteins (BRCA1, HRAS, PTEN, and ERK2) with high-quality structural data, we find that SBNA scores correlate strongly with deep mutagenesis data. When applied to 47 IRD genes with available high-quality crystal structure data, SBNA scores reliably identified disease-causing variants according to phenotype definitions from the ClinVar database. Finally, we applied this approach to 63 patients at Massachusetts Eye and Ear (MEE) with IRD but for whom no genetic cause had been identified. Untrained models built using SBNA scores and BLOSUM62 scores for IRD-associated genes successfully predicted the pathogenicity of novel variants (AUC = 0.851), allowing us to identify likely causative disease variants in 40 IRD patients. Model performance was further augmented by incorporating orthogonal data from EVE scores (AUC = 0.927), which are based on evolutionary multiple sequence alignments. In conclusion, SBNA can used to successfully identify variants as causal of disease in human proteins and may help predict variants causative of IRD in an unbiased fashion.

Rosace: a robust deep mutational scanning analysis framework employing position and mean-variance shrinkage

Deep mutational scanning (DMS) measures the effects of thousands of genetic variants in a protein simultaneously. The small sample size renders classical statistical methods ineffective. For example, p-values cannot be correctly calibrated when treating variants independently. We propose Rosace, a Bayesian framework for analyzing growth-based DMS data. Rosace leverages amino acid position information to increase power and control the false discovery rate by sharing information across parameters via shrinkage. We also developed Rosette for simulating the distributional properties of DMS. We show that Rosace is robust to the violation of model assumptions and is more powerful than existing tools.

Development of multiplexed orthogonal base editor (MOBE) systems

Base editors (BEs) enable efficient, programmable installation of point mutations while avoiding the use of double-strand breaks. Simultaneous application of two or more different BEs, such as an adenine BE (which converts A·T base pairs to G·C) and a cytosine BE (which converts C·G base pairs to T·A), is not feasible because guide RNA crosstalk results in non-orthogonal editing, with all BEs modifying all target loci. Here we engineer both adenine BEs and cytosine BEs that can be orthogonally multiplexed by using RNA aptamer-coat protein systems to recruit the DNA-modifying enzymes directly to the guide RNAs. We generate four multiplexed orthogonal BE systems that enable rates of precise co-occurring edits of up to 7.1% in the same DNA strand without enrichment or selection strategies. The addition of a fluorescent enrichment strategy increases co-occurring edit rates up to 24.8% in human cells. These systems are compatible with expanded protospacer adjacent motif and high-fidelity Cas9 variants, function well in multiple cell types, have equivalent or reduced off-target propensities compared with their parental systems and can model disease-relevant point mutation combinations.

Analyzing the functional effects of DNA variants with gene editing

Continual advancements in genomics have led to an ever-widening disparity between the rate of discovery of genetic variants and our current understanding of their functions and potential roles in disease. Systematic methods for phenotyping DNA variants are required to effectively translate genomics data into improved outcomes for patients with genetic diseases. To make the biggest impact, these approaches must be scalable and accurate, faithfully reflect disease biology, and define complex disease mechanisms. We compare current methods to analyze the function of variants in their endogenous DNA context using genome editing strategies, such as saturation genome editing, base editing and prime editing. We discuss how these technologies can be linked to high-content readouts to gain deep mechanistic insights into variant effects. Finally, we highlight key challenges that need to be addressed to bridge the genotype to phenotype gap, and ultimately improve the diagnosis and treatment of genetic diseases.

A systematic analysis of regression models for protein engineering

To optimize proteins for particular traits holds great promise for industrial and pharmaceutical purposes. Machine Learning is increasingly applied in this field to predict properties of proteins, thereby guiding the experimental optimization process. A natural question is: How much progress are we making with such predictions, and how important is the choice of regressor and representation? In this paper, we demonstrate that different assessment criteria for regressor performance can lead to dramatically different conclusions, depending on the choice of metric, and how one defines generalization. We highlight the fundamental issues of sample bias in typical regression scenarios and how this can lead to misleading conclusions about regressor performance. Finally, we make the case for the importance of calibrated uncertainty in this domain.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

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.

Minimum information and guidelines for reporting a multiplexed assay of variant effect

Multiplexed assays of variant effect (MAVEs) have emerged as a powerful approach for interrogating thousands of genetic variants in a single experiment. The flexibility and widespread adoption of these techniques across diverse disciplines have led to a heterogeneous mix of data formats and descriptions, which complicates the downstream use of the resulting datasets. To address these issues and promote reproducibility and reuse of MAVE data, we define a set of minimum information standards for MAVE data and metadata and outline a controlled vocabulary aligned with established biomedical ontologies for describing these experimental designs.

Clinical implications of VUS reclassification in a single-centre series from application of ACMG/AMP classification rules specified for BRCA1/2

BACKGROUND

BRCA1/2 testing is crucial to guide clinical decisions in patients with hereditary breast/ovarian cancer, but detection of variants of uncertain significance (VUSs) prevents proper management of carriers. The ENIGMA (Evidence-based Network for the Interpretation of Germline Mutant Alleles) BRCA1/2 Variant Curation Expert Panel (VCEP) has recently developed BRCA1/2 variant classification guidelines consistent with ClinGen processes, specified against the ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular-Pathology) classification framework.

METHODS

The ClinGen-approved BRCA1/2-specified ACMG/AMP classification guidelines were applied to BRCA1/2 VUSs identified from 2011 to 2022 in a series of patients, retrieving information from the VCEP documentation, public databases, literature and ENIGMA unpublished data. Then, we critically re-evaluated carrier families based on new results and checked consistency of updated classification with main sources for clinical interpretation of BRCA1/2 variants.

RESULTS

Among 166 VUSs detected in 231 index cases, 135 (81.3%) found in 197 index cases were classified by applying BRCA1/2-specified ACMG/AMP criteria: 128 (94.8%) as Benign/Likely Benign and 7 (5.2%) as Pathogenic/Likely Pathogenic. The average time from the first report as 'VUS' to classification using this approach was 49.4 months. Considering that 15 of these variants found in 64 families had already been internally reclassified prior to this work, this study provided 121 new reclassifications among the 151 (80.1%) remaining VUSs, relevant to 133/167 (79.6%) families.

CONCLUSIONS

These results demonstrated the effectiveness of new BRCA1/2 ACMG/AMP classification guidelines for VUS classification within a clinical cohort, and their important clinical impact. Furthermore, they suggested a cadence of no more than 3 years for regular review of VUSs, which however requires time, expertise and resources.

Defining and Reducing Variant Classification Disparities

Background

Multiplexed Assays of Variant Effects (MAVEs) can test all possible single variants in a gene of interest. The resulting saturation-style data may help resolve variant classification disparities between populations, especially for variants of uncertain significance (VUS).

Methods

We analyzed clinical significance classifications in 213,663 individuals of European-like genetic ancestry versus 206,975 individuals of non-European-like genetic ancestry from All of Us and the Genome Aggregation Database. Then, we incorporated clinically calibrated MAVE data into the Clinical Genome Resource's Variant Curation Expert Panel rules to automate VUS reclassification for BRCA1, TP53, and PTEN .

Results

Using two orthogonal statistical approaches, we show a higher prevalence ( p ≤5.95e-06) of VUS in individuals of non-European-like genetic ancestry across all medical specialties assessed in all three databases. Further, in the non-European-like genetic ancestry group, higher rates of Benign or Likely Benign and variants with no clinical designation ( p ≤2.5e-05) were found across many medical specialties, whereas Pathogenic or Likely Pathogenic assignments were higher in individuals of European-like genetic ancestry ( p ≤2.5e-05). Using MAVE data, we reclassified VUS in individuals of non-European-like genetic ancestry at a significantly higher rate in comparison to reclassified VUS from European-like genetic ancestry ( p =9.1e-03) effectively compensating for the VUS disparity. Further, essential code analysis showed equitable impact of MAVE evidence codes but inequitable impact of allele frequency ( p =7.47e-06) and computational predictor ( p =6.92e-05) evidence codes for individuals of non-European-like genetic ancestry.

Conclusions

Generation of saturation-style MAVE data should be a priority to reduce VUS disparities and produce equitable training data for future computational predictors.

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.

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.

Precise large-fragment deletions in mammalian cells and mice generated by dCas9-controlled CRISPR/Cas3

Currently, the Cas9 and Cas12a systems are widely used for genome editing, but their ability to precisely generate large chromosome fragment deletions is limited. Type I-E CRISPR mediates broad and unidirectional DNA degradation, but controlling the size of Cas3-mediated DNA deletions has proven elusive thus far. Here, we demonstrate that the endonuclease deactivation of Cas9 (dCas9) can precisely control Cas3-mediated large-fragment deletions in mammalian cells. In addition, we report the elimination of the Y chromosome and precise retention of the Sry gene in mice using CRISPR/Cas3 and dCas9-controlled CRISPR/Cas3, respectively. In conclusion, dCas9-controlled CRISPR/Cas3-mediated precise large-fragment deletion provides an approach for establishing animal models by chromosome elimination. This method also holds promise as a potential therapeutic strategy for treating fragment mutations or human aneuploidy diseases that involve additional chromosomes.