Splicing predictions, minigene analyses, and ACMG-AMP clinical classification of 42 germline PALB2 splice-site variants

Comprehensive splicing analysis of the alternatively spliced CHEK2 exons 8 and 10 reveals three enhancer/silencer-rich regions and 38 spliceogenic variants

Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Reclassifying BRCA1 c.4358-2A > G and BRCA2 c.475 + 5G > C variants from "Uncertain Significance" to "Pathogenic" based on minigene assays and clinical evidence

BACKGROUND

Pathogenic variants in BRCA genes play a crucial role in the pathogenesis of ovarian cancer. Intronic variants of uncertain significance (VUS) may contribute to pathogenicity by affecting splicing. Currently, the significance of many intronic variants in BRCA has not been clarified, impacting patient treatment strategies and the management of familial cases.

METHOD

A retrospective study was conducted to analyze BRCA intronic VUS in a cohort of 707 unrelated ovarian cancer patients at a single institution from 2018 to 2023. Three splicing predictors were employed to analyze detected intronic VUS. Variants predicted to have splicing alterations were selected for further validation through minigene assays. Patient and familial investigations were conducted to comprehend cancer incidence within pedigrees and the application of poly (ADP-ribose) polymerase inhibitors (PARPi) by the patients. In accordance with the guidelines of the American College of Medical Genetics and Genomics (ACMG), the intronic VUS were reclassified based on minigene assay results and clinical evidence.

RESULT

Approximately 9.8% (69/707) of patients were identified as carriers of 67 different VUS in BRCA1/2, with four intronic variants accounting for 6% (4/67) of all VUS. Splicing predictors indicated potential splicing alterations in splicing for BRCA1 c.4358-2A>G and BRCA2 c.475+5G>C variants. Minigene assays utilizing the pSPL3 exon trapping vector revealed that these variants induced changes in splicing sites and frameshift, resulting in premature termination of translation (p. Ala1453Glyfs and p. Pro143Glyfs). According to ACMG guidelines, BRCA1 c.4358-2A>G and BRCA2 c.475+5G>C were reclassified as pathogenic variants. Pedigree investigations were conducted on patients with BRCA1 c.4358-2A>G variant, and the detailed utilization of PARPi provided valuable insights into research on PARPi resistance.

CONCLUSION

Two intronic VUS were reclassified as pathogenic variants. A precise classification of variants is crucial for the effective treatment and management of both patients and healthy carriers.

Systematic Minigene-Based Splicing Analysis and Tentative Clinical Classification of 52 CHEK2 Splice-Site Variants

BACKGROUND

Disrupted pre-mRNA splicing is a frequent deleterious mechanism in hereditary cancer. We aimed to functionally analyze candidate spliceogenic variants of the breast cancer susceptibility gene CHEK2 by splicing reporter minigenes.

METHODS

A total of 128 CHEK2 splice-site variants identified in the Breast Cancer After Diagnostic Gene Sequencing (BRIDGES) project (https://cordis.europa.eu/project/id/634935) were analyzed with MaxEntScan and subsetted to 52 variants predicted to impact splicing. Three CHEK2 minigenes, which span all 15 exons, were constructed and validated. The 52 selected variants were then genetically engineered into the minigenes and assayed in MCF-7 (human breast adenocarcinoma) cells.

RESULTS

Of 52 variants, 46 (88.5%) impaired splicing. Some of them led to complex splicing patterns with up to 11 different transcripts. Thirty-four variants induced splicing anomalies without any trace or negligible amounts of the full-length transcript. A total of 89 different transcripts were annotated, which derived from different events: single- or multi-exon skipping, alternative site-usage, mutually exclusive exon inclusion, intron retention or combinations of the abovementioned events. Fifty-nine transcripts were predicted to introduce premature termination codons, 7 kept the original open-reading frame, 5 removed the translation start codon, 6 affected the 5'UTR (Untranslated Region), and 2 included missense variations. Analysis of variant c.684-2A > G revealed the activation of a non-canonical TG-acceptor site and exon 6 sequences critical for its recognition.

CONCLUSIONS

Incorporation of minigene read-outs into an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme allowed us to classify 32 CHEK2 variants (27 pathogenic/likely pathogenic and 5 likely benign). However, 20 variants (38%) remained of uncertain significance, reflecting in part the complex splicing patterns of this gene.

Finding significance: New perspectives in variant classification of the RAD51 regulators, BRCA2 and beyond

For many individuals harboring a variant of uncertain functional significance (VUS) in a homologous recombination (HR) gene, their risk of developing breast and ovarian cancer is unknown. Integral to the process of HR are BRCA1 and regulators of the central HR protein, RAD51, including BRCA2, PALB2, RAD51C and RAD51D. Due to advancements in sequencing technology and the continued expansion of cancer screening panels, the number of VUS identified in these genes has risen significantly. Standard practices for variant classification utilize different types of predictive, population, phenotypic, allelic and functional evidence. While variant analysis is improving, there remains a struggle to keep up with demand. Understanding the effects of an HR variant can aid in preventative care and is critical for developing an effective cancer treatment plan. In this review, we discuss current perspectives in the classification of variants in the breast and ovarian cancer genes BRCA1, BRCA2, PALB2, RAD51C and RAD51D.

ParSE-seq: A Calibrated Multiplexed Assay to Facilitate the Clinical Classification of Putative Splice-altering Variants

Background

Interpreting the clinical significance of putative splice-altering variants outside 2-base pair canonical splice sites remains difficult without functional studies.

Methods

We developed Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed minigene-based assay, to test variant effects on RNA splicing quantified by high-throughput sequencing. We studied variants in SCN5A, an arrhythmia-associated gene which encodes the major cardiac voltage-gated sodium channel. We used the computational tool SpliceAI to prioritize exonic and intronic candidate splice variants, and ClinVar to select benign and pathogenic control variants. We generated a pool of 284 barcoded minigene plasmids, transfected them into Human Embryonic Kidney (HEK293) cells and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), sequenced the resulting pools of splicing products, and calibrated the assay to the American College of Medical Genetics and Genomics scheme. Variants were interpreted using the calibrated functional data, and experimental data were compared to SpliceAI predictions. We further studied some splice-altering missense variants by cDNA-based automated patch clamping (APC) in HEK cells and assessed splicing and sodium channel function in CRISPR-edited iPSC-CMs.

Results

ParSE-seq revealed the splicing effect of 224 SCN5A variants in iPSC-CMs and 244 variants in HEK293 cells. The scores between the cell types were highly correlated (R2=0.84). In iPSCs, the assay had concordant scores for 21/22 benign/likely benign and 24/25 pathogenic/likely pathogenic control variants from ClinVar. 43/112 exonic variants and 35/70 intronic variants with determinate scores disrupted splicing. 11 of 42 variants of uncertain significance were reclassified, and 29 of 34 variants with conflicting interpretations were reclassified using the functional data. SpliceAI computational predictions correlated well with experimental data (AUC = 0.96). We identified 20 unique SCN5A missense variants that disrupted splicing, and 2 clinically observed splice-altering missense variants of uncertain significance had normal function when tested with the cDNA-based APC assay. A splice-altering intronic variant detected by ParSE-seq, c.1891-5C>G, also disrupted splicing and sodium current when introduced into iPSC-CMs at the endogenous locus by CRISPR editing.

Conclusions

ParSE-seq is a calibrated, multiplexed, high-throughput assay to facilitate the classification of candidate splice-altering variants.

Overview of the Genetic Causes of Hereditary Breast and Ovarian Cancer Syndrome in a Large French Patient Cohort

The use of multigene panel testing for patients with a predisposition to Hereditary Breast and Ovarian Cancer syndrome (HBOC) is increasing as the identification of mutations is useful for diagnosis and disease management. Here, we conducted a retrospective analysis of BRCA1/2 and non-BRCA gene sequencing in 4630 French HBOC suspected patients. Patients were investigated using a germline cancer panel including the 13 genes defined by The French Genetic and Cancer Group (GGC)-Unicancer. In the patients analyzed, 528 pathogenic and likely pathogenic variants (P/LP) were identified, including BRCA1 (n = 203, 38%), BRCA2 (n = 198, 37%), PALB2 (n = 46, 9%), RAD51C (n = 36, 7%), TP53 (n = 16, 3%), and RAD51D (n = 13, 2%). In addition, 35 novel (P/LP) variants, according to our knowledge, were identified, and double mutations in two distinct genes were found in five patients. Interestingly, retesting a subset of BRCA1/2-negative individuals with an expanded panel produced clinically relevant results in 5% of cases. Additionally, combining in silico (splicing impact prediction tools) and in vitro analyses (RT-PCR and Sanger sequencing) highlighted the deleterious impact of four candidate variants on splicing and translation. Our results present an overview of pathogenic variations of HBOC genes in the southeast of France, emphasizing the clinical relevance of cDNA analysis and the importance of retesting BRCA-negative individuals with an expanded panel.

Minigene Assay as an Effective Molecular Diagnostic Strategy in Determining the Pathogenicity of Noncanonical Splice-Site Variants in FLCN

Primary spontaneous pneumothorax (PSP) or pulmonary cyst is one of the manifestations of Birt-Hogg-Dubé syndrome, which is caused by pathogenic variants in FLCN gene. Genetic testing in patients with PSP identifies a certain number of missense or intronic variants. These variants are usually considered as variants of uncertain significance, whose functional interpretations pose a challenge in clinical genetics. To improve recognition of pathogenic splice-altering variants in FLCN gene, computational tools are used to prioritize potential splice-altering variants and then a hybrid minigene assay is performed to verify the RNA splicing pattern. Herein, variants in FLCN exon 11 and its flanking sequence are focused. Eight variants detected in 11 patients with PSP are evaluated, and six variants are prioritized by in silico tools as potential splice-altering variants of uncertain significance. Four variants (c.1177-5_1177-3delCTC, c.1292_1300+4del, c.1300+4C>T, and c.1300+5G>A) are demonstrated by minigene assay to alter RNA splicing of FLCN, and the last three of them are novel. RT-PCR of patient-derived RNA gives consistent results. Genotype-phenotype correlation analysis in patients with PSP with these variants demonstrates good concordance. Our results underline the importance of RNA analysis, which could provide molecular evidence for pathogenicity of a variant, and provide essential information for the clinical interpretation of variants. Combining the clinical information, a definitive diagnosis could be made.

Functional Assays Reclassify Suspected Splice-Altering Variants of Uncertain Significance in Mendelian Channelopathies

BACKGROUND

Rare protein-altering variants in SCN5A, KCNQ1, and KCNH2 are major causes of Brugada syndrome and the congenital long QT syndrome. While splice-altering variants lying outside 2-bp canonical splice sites can cause these diseases, their role remains poorly described. We implemented 2 functional assays to assess 12 recently reported putative splice-altering variants of uncertain significance and 1 likely pathogenic variant without functional data observed in Brugada syndrome and long QT syndrome probands.

METHODS

We deployed minigene assays to assess the splicing consequences of 10 variants. Three variants incompatible with the minigene approach were introduced into control induced pluripotent stem cells by CRISPR genome editing. We differentiated cells into induced pluripotent stem cell-derived cardiomyocytes and studied splicing outcomes by reverse transcription-polymerase chain reaction. We used the American College of Medical Genetics and Genomics functional assay criteria (PS3/BS3) to reclassify variants.

RESULTS

We identified aberrant splicing, with presumed disruption of protein sequence, in 8/10 variants studied using the minigene assay and 1/3 studied in induced pluripotent stem cell-derived cardiomyocytes. We reclassified 8 variants of uncertain significance to likely pathogenic, 1 variant of uncertain significance to likely benign, and 1 likely pathogenic variant to pathogenic.

CONCLUSIONS

Functional assays reclassified splice-altering variants outside canonical splice sites in Brugada Syndrome- and long QT syndrome-associated genes.

Clinical, splicing, and functional analysis to classify BRCA2 exon 3 variants: Application of a points-based ACMG/AMP approach

Skipping of BRCA2 exon 3 (∆E3) is a naturally occurring splicing event, complicating clinical classification of variants that may alter ∆E3 expression. This study used multiple evidence types to assess pathogenicity of 85 variants in/near BRCA2 exon 3. Bioinformatically predicted spliceogenic variants underwent mRNA splicing analysis using minigenes and/or patient samples. ∆E3 was measured using quantitative analysis. A mouse embryonic stem cell (mESC) based assay was used to determine the impact of 18 variants on mRNA splicing and protein function. For each variant, population frequency, bioinformatic predictions, clinical data, and existing mRNA splicing and functional results were collated. Variant class was assigned using a gene-specific adaptation of ACMG/AMP guidelines, following a recently proposed points-based system. mRNA and mESC analysis combined identified six variants with transcript and/or functional profiles interpreted as loss of function. Cryptic splice site use for acceptor site variants generated a transcript encoding a shorter protein that retains activity. Overall, 69/85 (81%) variants were classified using the points-based approach. Our analysis shows the value of applying gene-specific ACMG/AMP guidelines using a points-based approach and highlights the consideration of cryptic splice site usage to appropriately assign PVS1 code strength.

Splicing Analysis of 16 PALB2 ClinVar Variants by Minigene Assays: Identification of Six Likely Pathogenic Variants

PALB2 loss-of-function variants are associated with significant increased risk of breast cancer as well as other types of tumors. Likewise, splicing disruptions are a common mechanism of disease susceptibility. Indeed, we previously showed, by minigene assays, that 35 out of 42 PALB2 variants impaired splicing. Taking advantage of one of these constructs (mgPALB2_ex1-3), we proceeded to analyze other variants at exons 1 to 3 reported at the ClinVar database. Thirty-one variants were bioinformatically analyzed with MaxEntScan and SpliceAI. Then, 16 variants were selected for subsequent RNA assays. We identified a total of 12 spliceogenic variants, 11 of which did not produce any trace of the expected minigene full-length transcript. Interestingly, variant c.49-1G > A mimicked previous outcomes in patient RNA (transcript ∆(E2p6)), supporting the reproducibility of the minigene approach. A total of eight variant-induced transcripts were characterized, three of which (∆(E1q17), ∆(E3p11), and ∆(E3)) were predicted to introduce a premature termination codon and to undergo nonsense-mediated decay, and five (▼(E1q9), ∆(E2p6), ∆(E2), ▼(E3q48)-a, and ▼(E3q48)-b) maintained the reading frame. According to an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme, which integrates mgPALB2 data, six PALB2 variants were classified as pathogenic/likely pathogenic, five as VUS, and five as likely benign. Furthermore, five ±1,2 variants were catalogued as VUS because they produced significant proportions of in-frame transcripts of unknown impact on protein function.

Minigene-based splicing analysis and ACMG/AMP-based tentative classification of 56 ATM variants

The ataxia telangiectasia-mutated (ATM) protein is a major coordinator of the DNA damage response pathway. ATM loss-of-function variants are associated with 2-fold increased breast cancer risk. We aimed at identifying and classifying spliceogenic ATM variants detected in subjects of the large-scale sequencing project BRIDGES. A total of 381 variants at the intron-exon boundaries were identified, 128 of which were predicted to be spliceogenic. After further filtering, we ended up selecting 56 variants for splicing analysis. Four functional minigenes (mgATM) spanning exons 4-9, 11-17, 25-29, and 49-52 were constructed in the splicing plasmid pSAD. Selected variants were genetically engineered into the four constructs and assayed in MCF-7/HeLa cells. Forty-eight variants (85.7%) impaired splicing, 32 of which did not show any trace of the full-length (FL) transcript. A total of 43 transcripts were identified where the most prevalent event was exon/multi-exon skipping. Twenty-seven transcripts were predicted to truncate the ATM protein. A tentative ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme that integrates mgATM data allowed us to classify 29 ATM variants as pathogenic/likely pathogenic and seven variants as likely benign. Interestingly, the likely pathogenic variant c.1898+2T>G generated 13% of the minigene FL-transcript due to the use of a noncanonical GG-5'-splice-site (0.014% of human donor sites). Circumstantial evidence in three ATM variants (leakiness uncovered by our mgATM analysis together with clinical data) provides some support for a dosage-sensitive expression model in which variants producing ≥30% of FL-transcripts would be predicted benign, while variants producing ≤13% of FL-transcripts might be pathogenic. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Cutting-Edge AI Technologies Meet Precision Medicine to Improve Cancer Care

To provide precision medicine for better cancer care, researchers must work on clinical patient data, such as electronic medical records, physiological measurements, biochemistry, computerized tomography scans, digital pathology, and the genetic landscape of cancer tissue. To interpret big biodata in cancer genomics, an operational flow based on artificial intelligence (AI) models and medical management platforms with high-performance computing must be set up for precision cancer genomics in clinical practice. To work in the fast-evolving fields of patient care, clinical diagnostics, and therapeutic services, clinicians must understand the fundamentals of the AI tool approach. Therefore, the present article covers the following four themes: (i) computational prediction of pathogenic variants of cancer susceptibility genes; (ii) AI model for mutational analysis; (iii) single-cell genomics and computational biology; (iv) text mining for identifying gene targets in cancer; and (v) the NVIDIA graphics processing units, DRAGEN field programmable gate arrays systems and AI medical cloud platforms in clinical next-generation sequencing laboratories. Based on AI medical platforms and visualization, large amounts of clinical biodata can be rapidly copied and understood using an AI pipeline. The use of innovative AI technologies can deliver more accurate and rapid cancer therapy targets.

Minigene Splicing Assays Identify 20 Spliceogenic Variants of the Breast/Ovarian Cancer Susceptibility Gene RAD51C

RAD51C loss-of-function variants are associated with an increased risk of breast and ovarian cancers. Likewise, splicing disruptions are a frequent mechanism of gene inactivation. Taking advantage of a previous splicing-reporter minigene with exons 2-8 (mgR51C_ex2-8), we proceeded to check its impact on the splicing of candidate ClinVar variants. A total of 141 RAD51C variants at the intron/exon boundaries were analyzed with MaxEntScan. Twenty variants were selected and genetically engineered into the wild-type minigene. All the variants disrupted splicing, and 18 induced major splicing anomalies without any trace or minimal amounts (<2.4%) of the minigene full-length (FL) transcript. Twenty-seven transcripts (including the wild-type and r.904A FL transcripts) were identified by fluorescent fragment electrophoresis; of these, 14 were predicted to truncate the RAD51C protein, 3 kept the reading frame, and 8 minor isoforms (1.1−4.7% of the overall expression) could not be characterized. Finally, we performed a tentative interpretation of the variants according to an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme, classifying 16 variants as likely pathogenic. Minigene assays have been proven as valuable tools for the initial characterization of potential spliceogenic variants. Hence, minigene mgR51C_ex2-8 provided useful splicing data for 40 RAD51C variants.

El-Hattab-Alkuraya syndrome caused by biallelic WDR45B pathogenic variants: further delineation of the phenotype and genotype

Homozygous pathogenic variants in WDR45B were first identified in six subjects from three unrelated families with global development delay, refractory seizures, spastic quadriplegia, and brain malformations. Since the initial report in 2018, no further cases have been described. In this report, we present 12 additional individuals from seven unrelated families and their clinical, radiological, and molecular findings. Six different variants in WDR45B were identified, five of which are novel. Microcephaly and global developmental delay were observed in all subjects, and seizures and spastic quadriplegia in most. Common findings on brain imaging include cerebral atrophy, ex-vacuo ventricular dilatation, brainstem volume loss, and symmetric under-opercularization. El-Hattab-Alkuraya syndrome is associated with a consistent phenotype characterized by early onset cerebral atrophy resulting in microcephaly, developmental delay, spastic quadriplegia, and seizures. The phenotype appears to be more severe among individuals with loss-of-function variants whereas those with missense variants were less severely affected suggesting a potential genotype-phenotype correlation in this disorder. A brain imaging pattern emerges which is consistent among individuals with loss-of-function variants and could potentially alert the neuroradiologists or clinician to consider WDR45B-related El-Hattab-Alkuraya syndrome. This article is protected by copyright. All rights reserved.