Identification of Eight Spliceogenic Variants in BRCA2 Exon 16 by Minigene Assays

Real-world evidence: Risdiplam in a patient with spinal muscular atrophy type I with a novel splicing mutation and one SMN2 copy

Spinal muscular atrophy (SMA), which results from the deletion or/and mutation in the SMN1 gene, is an autosomal recessive neuromuscular disorder that leads to weakness and muscle atrophy. SMN2 is a paralogous gene of SMN1. SMN2 copy number affects the severity of SMA, but its role in patients treated with disease modifying therapies is unclear. The most appropriate individualized treatment for SMA has not yet been determined. Here, we reported a case of SMA type I with normal breathing and swallowing function. We genetically confirmed that this patient had a compound heterozygous variant: one deleted SMN1 allele and a novel splice mutation c.628-3T>G in the retained allele, with one SMN2 copy. Patient-derived sequencing of 4 SMN1 cDNA clones showed that this intronic single transversion mutation results in an alternative exon (e)5 3' splice site, which leads to an additional 2 nucleotides (AG) at the 5' end of e5, thereby explaining why the patient with only one copy of SMN2 had a mild clinical phenotype. Additionally, a minigene assay of wild type and mutant SMN1 in HEK293T cells also demonstrated that this transversion mutation induced e5 skipping. Considering treatment cost and goals of avoiding pain caused by injections and starting treatment as early as possible, risdiplam was prescribed for this patient. However, the patient showed remarkable clinical improvements after treatment with risdiplam for 7 months despite carrying only one copy of SMN2. This study is the first report on the treatment of risdiplam in a patient with one SMN2 copy in a real-world setting. These findings expand the mutation spectrum of SMA and provide accurate genetic counseling information, as well as clarify the molecular mechanism of careful genotype-phenotype correlation of the patient.

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.

Identified eleven exon variants in PKD1 and PKD2 genes that altered RNA splicing by minigene assay

BACKGROUND

Autosomal dominant polycystic kidney disease (ADPKD) is a common monogenic multisystem disease caused primarily by mutations in the PKD1 gene or PKD2 gene. There is increasing evidence that some of these variants, which are described as missense, synonymous or nonsense mutations in the literature or databases, may be deleterious by affecting the pre-mRNA splicing process.

RESULTS

This study aimed to determine the effect of these PKD1 and PKD2 variants on exon splicing combined with predictive bioinformatics tools and minigene assay. As a result, among the 19 candidate single nucleotide alterations, 11 variants distributed in PKD1 (c.7866C > A, c.7960A > G, c.7979A > T, c.7987C > T, c.11248C > G, c.11251C > T, c.11257C > G, c.11257C > T, c.11346C > T, and c.11393C > G) and PKD2 (c.1480G > T) were identified to result in exon skipping.

CONCLUSIONS

We confirmed that 11 variants in the gene of PKD1 and PKD2 affect normal splicing by interfering the recognition of classical splicing sites or by disrupting exon splicing enhancers and generating exon splicing silencers. This is the most comprehensive study to date on pre-mRNA splicing of exonic variants in ADPKD-associated disease-causing genes in consideration of the increasing number of identified variants in PKD1 and PKD2 gene in recent years. These results emphasize the significance of assessing the effect of exon single nucleotide variants in ADPKD at the mRNA level.

Molecular Basis of Unequal Alternative Splicing of Human SCD5 and Its Alteration by Natural Genetic Variations

Alternative splicing (AS) is a major means of post-transcriptional control of gene expression, and provides a dynamic versatility of protein isoforms. Cancer-related AS disorders have diagnostic, prognostic and therapeutic values. Changes in the expression and AS of human stearoyl-CoA desaturase-5 (SCD5) are promising specific tumor markers, although the transcript variants (TVs) of the gene have not yet been confirmed. Our in silico, in vitro and in vivo study focuses on the distribution of SCD5 TVs (A and B) in human tissues, the functionality of the relevant splice sites, and their modulation by certain single-nucleotide variations (SNVs). An order of magnitude higher SCD5A expression was found compared with SCD5B. This unequal splicing is attributed to a weaker recognition of the SCD5B-specific splicing acceptor site, based on predictions confirmed by an optimized minigene assay. The pronounced dominance of SCD5A was largely modified (rs1430176385_A, rs1011850309_A) or even inverted (rs1011850309_C) by natural SNVs at the TV-specific splice sites. Our results provide long missing data on the proportion of SCD5 TVs in human tissues and reveal mutation-driven changes in SCD5 AS, potentially affecting tumor-associated reprogramming of lipid metabolism, thus having prognostic significance, which may be utilized for novel and personalized therapeutic approaches.

Case Review: Whole-Exome Sequencing Analyses Identify Carriers of a Known Likely Pathogenic Intronic BRCA1 Variant in Ovarian Cancer Cases Clinically Negative for Pathogenic BRCA1 and BRCA2 Variants

Background: Detecting pathogenic intronic variants resulting in aberrant splicing remains a challenge in routine genetic testing. We describe germline whole-exome sequencing (WES) analyses and apply in silico predictive tools of familial ovarian cancer (OC) cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants. Methods: WES data from 27 familial OC cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants and 53 sporadic early-onset OC cases were analyzed for pathogenic variants in BRCA1 or BRCA2. WES data from carriers of pathogenic BRCA1 or BRCA2 variants were analyzed for pathogenic variants in 10 other OC predisposing genes. Loss of heterozygosity analysis was performed on tumor DNA from variant carriers. Results: BRCA1 c.5407-25T>A intronic variant, identified in two affected sisters and one sporadic OC case, is predicted to create a new splice effecting transcription of BRCA1. WES data from BRCA1 c.5407-25T>A carriers showed no evidence of pathogenic variants in other OC predisposing genes. Sequencing the tumor DNA from the variant carrier showed complete loss of the wild-type allele. Conclusions: The findings support BRCA1 c.5407-25T>A as a likely pathogenic variant and highlight the importance of investigating intronic sequences as causal variants in OC families where the involvement of BRCA1 is highly suggestive.

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

PALB2 loss-of-function variants confer high risk of developing breast cancer. Here we present a systematic functional analysis of PALB2 splice-site variants detected in approximately 113,000 women in the large-scale sequencing project Breast Cancer After Diagnostic Gene Sequencing (BRIDGES; https://bridges-research.eu/). Eighty-two PALB2 variants at the intron-exon boundaries were analyzed with MaxEntScan. Forty-two variants were selected for the subsequent splicing functional assays. For this purpose, three splicing reporter minigenes comprising exons 1-12 were constructed. The 42 potential spliceogenic variants were introduced into the minigenes by site-directed mutagenesis and assayed in MCF-7/MDA-MB-231 cells. Splicing anomalies were observed in 35 variants, 23 of which showed no traces or minimal amounts of the expected full-length transcripts of each minigene. More than 30 different variant-induced transcripts were characterized, 23 of which were predicted to truncate the PALB2 protein. The pathogenicity of all variants was interpreted according to an in-house adaptation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) variant classification scheme. Up to 23 variants were classified as pathogenic/likely pathogenic. Remarkably, three ±1,2 variants (c.49-2A>T, c.108+2T>C, and c.211+1G>A) were classified as variants of unknown significance, as they produced significant amounts of either in-frame transcripts of unknown impact on the PALB2 protein function or the minigene full-length transcripts. In conclusion, we have significantly contributed to the ongoing effort of identifying spliceogenic variants in the clinically relevant PALB2 cancer susceptibility gene. Moreover, we suggest some approaches to classify the findings in accordance with the ACMG-AMP rationale. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Mutation of c.244G>T in NR5A1 gene causing 46, XY DSD by affecting RNA splicing

Objective

To identify the pathogenic mechanism of the c.244G>T mutation in NR5A1 gene found in a Chinese patient with 46, XY disorders of sex development (DSD). Subjects and methods: Genomic DNA was extracted from a Chinese 46, XY DSD patient. Targeted next-generation and Sanger sequencing were performed to investigate and validate the gene mutation causing 46, XY DSD, respectively. In silico tools were used to predict the pathogenicity of the variant. Dual luciferase reporter gene assay and minigene splicing reporter assay were used to identify the pathogenicity of the variant.

Results

A novel heterozygous variant, c.244G>T (p.Ala82Ser), in NR5A1 gene was detected in the 46, XY DSD patient. Four of five silico tools predicting pathogenicity of missense variants indicated that the variant was pathogenic. However, in vitro functional study showed that p.Ala82Ser did not affect the transcriptional activity of NR5A1. In silico tools predicting the potential splicing loci revealed that c.244G>T led to aberrant splicing of NR5A1 RNA. Minigene splicing reporter assay confirmed that c.244G>T resulted in the deletion of exon2 or deletion of 19 nucleotides in 3′ end of exon2.

Conclusions

Mutation of c.244G>T in NR5A1 results in 46, XY DSD by inducing abnormal splicing of NR5A1 RNA instead of amino acid substitution of NR5A1.

Identification of seven exonic variants in the SLC4A1, ATP6V1B1, and ATP6V0A4 genes that alter RNA splicing by minigene assay

Primary distal renal tubular acidosis (dRTA) is a rare tubular disease associated with variants in SLC4A1, ATP6V0A4, ATP6V1B1, FOXⅠ1, or WDR72 genes. Currently, there is growing evidence that all types of exonic variants can alter splicing regulatory elements, affecting the precursor messenger RNA (pre-mRNA) splicing process. This study was to determine the consequences of variants associated with dRTA on pre-mRNA splicing combined with predictive bioinformatics tools and minigene assay. As a result, among the 15 candidate variants, 7 variants distributed in SLC4A1 (c.1765C>T, p.Arg589Cys), ATP6V1B1 (c.368G>T, p.Gly123Val; c.370C>T, p.Arg124Trp; c.484G>T, p.Glu162* and c.1102G>A, p.Glu368Lys) and ATP6V0A4 genes (c.322C>T, p.Gln108* and c.1572G>A, p.Pro524Pro) were identified to result in complete or incomplete exon skipping by either disruption of exonic splicing enhancers (ESEs) and generation of exonic splicing silencers, or interference with the recognition of the classic splicing site, or both. To our knowledge, this is the first study on pre-mRNA splicing of exonic variants in the dRTA-related genes. These results highlight the importance of assessing the effects of exonic variants at the mRNA level and suggest that minigene analysis is an effective tool for evaluating the effects of splicing on variants in vitro.

RAD51D Aberrant Splicing in Breast Cancer: Identification of Splicing Regulatory Elements and Minigene-Based Evaluation of 53 DNA Variants

RAD51D loss-of-function variants increase lifetime risk of breast and ovarian cancer. Splicing disruption is a frequent pathogenic mechanism associated with variants in susceptibility genes. Herein, we have assessed the splicing and clinical impact of splice-site and exonic splicing enhancer (ESE) variants identified through the study of ~113,000 women of the BRIDGES cohort. A RAD51D minigene with exons 2-9 was constructed in splicing vector pSAD. Eleven BRIDGES splice-site variants (selected by MaxEntScan) were introduced into the minigene by site-directed mutagenesis and tested in MCF-7 cells. The 11 variants disrupted splicing, collectively generating 25 different aberrant transcripts. All variants but one produced negligible levels (<3.4%) of the full-length (FL) transcript. In addition, ESE elements of the alternative exon 3 were mapped by testing four overlapping exonic microdeletions (≥30-bp), revealing an ESE-rich interval (c.202_235del) with critical sequences for exon 3 recognition that might have been affected by germline variants. Next, 26 BRIDGES variants and 16 artificial exon 3 single-nucleotide substitutions were also assayed. Thirty variants impaired splicing with variable amounts (0-65.1%) of the FL transcript, although only c.202G>A demonstrated a complete aberrant splicing pattern without the FL transcript. On the other hand, c.214T>C increased efficiency of exon 3 recognition, so only the FL transcript was detected (100%). In conclusion, 41 RAD51D spliceogenic variants (28 of which were from the BRIDGES cohort) were identified by minigene assays. We show that minigene-based mapping of ESEs is a powerful approach for identifying ESE hotspots and ESE-disrupting variants. Finally, we have classified nine variants as likely pathogenic according to ACMG/AMP-based guidelines, highlighting the complex relationship between splicing alterations and variant interpretation.

Imprecise Medicine: BRCA2 Variants of Uncertain Significance (VUS), the Challenges and Benefits to Integrate a Functional Assay Workflow with Clinical Decision Rules

Pathological mutations in homology-directed repair (HDR) genes impact both future cancer risk and therapeutic options for patients. HDR is a high-fidelity DNA repair pathway for resolving DNA double-strand breaks throughout the genome. BRCA2 is an essential protein that mediates the loading of RAD51 onto resected DNA breaks, a key step in HDR. Germline mutations in BRCA2 are associated with an increased risk for breast, ovarian, prostate, and pancreatic cancer. Clinical findings of germline or somatic BRCA2 mutations in tumors suggest treatment with platinum agents or PARP inhibitors. However, when genetic analysis reveals a variant of uncertain significance (VUS) in the BRCA2 gene, precision medicine-based decisions become complex. VUS are genetic changes with unknown pathological impact. Current statistics indicate that between 10–20% of BRCA sequencing results are VUS, and of these, more than 50% are missense mutations. Functional assays to determine the pathological outcome of VUS are urgently needed to provide clinical guidance regarding cancer risk and treatment options. In this review, we provide a brief overview of BRCA2 functions in HDR, describe how BRCA2 VUS are currently assessed in the clinic, and how genetic and biochemical functional assays could be integrated into the clinical decision process. We suggest a multi-step workflow composed of robust and accurate functional assays to correctly evaluate the potential pathogenic or benign nature of BRCA2 VUS. Success in this precision medicine endeavor will offer actionable information to patients and their physicians.

Functional Analysis of the PCCA and PCCB Gene Variants Predicted to Affect Splicing

It is estimated that up to one-third of all variants causing inherited diseases affect splicing; however, their deleterious effects and roles in disease pathogenesis are often not fully characterized. Given their prevalence and the development of various antisense-based splice-modulating approaches, pathogenic splicing variants have become an important object of genomic medicine. To improve the accuracy of variant interpretation in public mutation repositories, we applied the minigene splicing assay to study the effects of 24 variants that were predicted to affect normal splicing in the genes associated with propionic acidemia (PA)—PCCA and PCCB. As a result, 13 variants (including one missense and two synonymous variants) demonstrated a significant alteration of splicing with the predicted deleterious effect at the protein level and were characterized as spliceogenic loss-of-function variants. The analysis of the available data for the studied variants and application of the American College of Medical Genetics and the Association for Molecular Pathology (ACMG/AMP) guidelines allowed us to precisely classify five of the variants and change the pathogenic status of nine. Using the example of the PA genes, we demonstrated the utility of the minigene splicing assay in the fast and effective assessment of the spliceogenic effect for identified variants and highlight the necessity of their standardized classification.

Comprehensive Functional Characterization and Clinical Interpretation of 20 Splice-Site Variants of the RAD51C Gene

Hereditary breast and/or ovarian cancer is a highly heterogeneous disease with more than 10 known disease-associated genes. In the framework of the BRIDGES project (Breast Cancer Risk after Diagnostic Gene Sequencing), the RAD51C gene has been sequenced in 60,466 breast cancer patients and 53,461 controls. We aimed at functionally characterizing all the identified genetic variants that are predicted to disrupt the splicing process. Forty RAD51C variants of the intron-exon boundaries were bioinformatically analyzed, 20 of which were selected for splicing functional assays. To test them, a splicing reporter minigene with exons 2 to 8 was designed and constructed. This minigene generated a full-length transcript of the expected size (1062 nucleotides), sequence, and structure (Vector exon V1- RAD51C exons_2-8- Vector exon V2). The 20 candidate variants were genetically engineered into the wild type minigene and functionally assayed in MCF-7 cells. Nineteen variants (95%) impaired splicing, while 18 of them produced severe splicing anomalies. At least 35 transcripts were generated by the mutant minigenes: 16 protein-truncating, 6 in-frame, and 13 minor uncharacterized isoforms. According to ACMG/AMP-based standards, 15 variants could be classified as pathogenic or likely pathogenic variants: c.404G > A, c.405-6T > A, c.571 + 4A > G, c.571 + 5G > A, c.572-1G > T, c.705G > T, c.706-2A > C, c.706-2A > G, c.837 + 2T > C, c.905-3C > G, c.905-2A > C, c.905-2_905-1del, c.965 + 5G > A, c.1026 + 5_1026 + 7del, and c.1026 + 5G > T.

Variant effect on splicing regulatory elements, branchpoint usage, and pseudoexonization: Strategies to enhance bioinformatic prediction using hereditary cancer genes as exemplars

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre-existing large-scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.

UGT1A1 Variants c.864+5G>T and c.996+2_996+5del of a Crigler-Najjar Patient Induce Aberrant Splicing in Minigene Assays

A large fraction of DNA variants impairs pre-mRNA splicing in human hereditary disorders. Crigler-Najjar syndrome (CNS) is characterized by a severe unconjugated hyperbilirubinemia caused by variants in the UGT1A1 gene. We previously reported one CNS-type II patient with two splice-site variants in trans (c.864+5G>T and c.996+2_996+5del). According to MaxEntScan, both disrupt their corresponding donor sites (c.864+5G>T: 6.99 → 2.28; c.996+2_996+5del: 5.96 → -11.02), so they were selected for subsequent functional tests. Given the unavailability of patient RNA, we constructed an UGT1A1 splicing-reporter minigene with exons 1-4 to characterize the underlying splicing anomaly. The variant c.996+2_996+5del generated two aberrant transcripts, Δ(E2) (exon 2 skipping/64%) and ▼(E2q135) (intron retention of 135-nt/36%), which lead to the loss of 18 conserved amino-acids and the gain of 45 new ones of a critical functional domain, respectively. The c.864+5G>T variant mainly produced the aberrant transcript Δ(E1q141) (141-nt deletion/70.4%) and the full-length isoform (29.6%). Δ(E1q141) would provoke the loss of 47 amino-acids of the N-terminal domain that encodes for substrate specificity. Thus, the three anomalous transcripts are likely to inactivate UGT1A1. Moreover, this patient is also homozygous for the promoter variant A(TA)7TAA that decreases UGT1A1 expression by 70%, so the full-length transcript produced by c.864+5G>T would be even more reduced (<9%), thus supporting the diagnosis of CNS-type II. Therefore, minigenes represent valuable tools for the functional and clinical classifications of genetic variants.

Splicing profile by capture RNA-seq identifies pathogenic germline variants in tumor suppressor genes

Germline variants in tumor suppressor genes (TSGs) can result in RNA mis-splicing and predisposition to cancer. However, identification of variants that impact splicing remains a challenge, contributing to a substantial proportion of patients with suspected hereditary cancer syndromes remaining without a molecular diagnosis. To address this, we used capture RNA-sequencing (RNA-seq) to generate a splicing profile of 18 TSGs (APC, ATM, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, MLH1, MSH2, MSH6, MUTYH, NF1, PALB2, PMS2, PTEN, RAD51C, RAD51D, and TP53) in 345 whole-blood samples from healthy donors. We subsequently demonstrated that this approach can detect mis-splicing by comparing splicing profiles from the control dataset to profiles generated from whole blood of individuals previously identified with pathogenic germline splicing variants in these genes. To assess the utility of our TSG splicing profile to prospectively identify pathogenic splicing variants, we performed concurrent capture DNA and RNA-seq in a cohort of 1000 patients with suspected hereditary cancer syndromes. This approach improved the diagnostic yield in this cohort, resulting in a 9.1% relative increase in the detection of pathogenic variants, demonstrating the utility of performing simultaneous DNA and RNA genetic testing in a clinical context.

Minigene Splicing Assays Identify 12 Spliceogenic Variants of BRCA2 Exons 14 and 15

A relevant fraction of BRCA2 variants is associated with splicing alterations and with an increased risk of hereditary breast and ovarian cancer (HBOC). In this work, we have carried out a thorough study of variants from BRCA2 exons 14 and 15 reported at mutation databases. A total of 294 variants from exons 14 and 15 and flanking intronic sequences were analyzed with the online splicing tools NNSplice and Human Splicing Finder. Fifty-three out of these 294 variants were selected as candidate splicing variants. All variants but one, were introduced into the minigene MGBR2_ex14-20 (with exons 14–20) by site-directed mutagenesis and assayed in MCF-7 cells. Twelve of the remaining 52 variants (23.1%) impaired splicing at different degrees, yielding from 5 to 100% of aberrant transcripts. Nine variants affected the natural acceptor or donor sites of both exons and three affected putative enhancers or silencers. Fluorescent capillary electrophoresis revealed at least 10 different anomalous transcripts: (E14q5), Δ (E14p10), Δ(E14p246), Δ(E14q256), Δ(E14), Δ(E15p12), Δ(E15p13), Δ(E15p83), Δ(E15) and a 942-nt fragment of unknown structure. All transcripts, except for Δ(E14q256) and Δ(E15p12), are expected to truncate the BRCA2 protein. Nine variants induced severe splicing aberrations with more than 90% of abnormal transcripts. Thus, according to the guidelines of the American College of Medical Genetics and Genomics, eight variants should be classified as pathogenic (c.7008-2A > T, c.7008-1G > A, c.7435+1G > C, c.7436-2A > T, c.7436-2A > G, c.7617+1G > A, c.7617+1G > T, and c.7617+2T > G), one as likely pathogenic (c.7008-3C > G) and three remain as variants of uncertain clinical significance or VUS (c.7177A > G, c.7447A > G and c.7501C > T). In conclusion, functional assays by minigenes constitute a valuable strategy to primarily check the splicing impact of DNA variants and their clinical interpretation. While bioinformatics predictions of splice site variants were accurate, those of enhancer or silencer variants were poor (only 3/23 spliceogenic variants) which showed weak impacts on splicing (∼5–16% of aberrant isoforms). So, the Exonic Splicing Enhancer and Silencer (ESE and ESS, respectively) prediction algorithms require further improvement.

Mis-splicing in breast cancer: identification of pathogenic BRCA2 variants by systematic minigene assays

Splicing disruption is a common mechanism of gene inactivation associated with germline variants of susceptibility genes. To study the role of BRCA2 mis-splicing in hereditary breast/ovarian cancer (HBOC), we performed a comprehensive analysis of variants from BRCA2 exons 2-9, as well as the initial characterization of the regulatory mechanisms of such exons. A pSAD-based minigene with exons 2-9 was constructed and validated in MCF-7 cells, producing the expected transcript (1016-nt/V1-BRCA2_exons_2-9-V2). DNA variants from mutational databases were analyzed by NNSplice and Human Splicing Finder softwares. To refine ESE-variant prediction, we mapped the regulatory regions through a functional strategy whereby 26 exonic microdeletions were introduced into the minigene and tested in MCF-7 cells. Thus, we identified nine spliceogenic ESE-rich intervals where ESE-variants may be located. Combining bioinformatics and microdeletion assays, 83 variants were selected and genetically engineered in the minigene. Fifty-three changes impaired splicing: 28 variants disrupted the canonical sites, four created new ones, 10 abrogated enhancers, eight created silencers and three caused a double-effect. Notably, nine spliceogenic-ESE variants were located within ESE-containing intervals. Capillary electrophoresis and sequencing revealed more than 23 aberrant transcripts, where exon skipping was the most common event. Interestingly, variant c.67G>A triggered the usage of a noncanonical GC-donor 4-nt upstream. Thirty-six variants that induced severe anomalies (>60% aberrant transcripts) were analyzed according to the ACMG guidelines. Thus, 28 variants were classified as pathogenic, five as likely pathogenic and three as variants of uncertain significance. Interestingly, 13 VUS were reclassified as pathogenic or likely pathogenic variants. In conclusion, a large fraction of BRCA2 variants (∼64%) provoked splicing anomalies lending further support to the high prevalence of this disease-mechanism. The low accuracy of ESE-prediction algorithms may be circumvented by functional ESE-mapping that represents an optimal strategy to identify spliceogenic ESE-variants. Finally, systematic functional assays by minigenes depict a valuable tool for the initial characterization of splicing anomalies and the clinical interpretation of variants. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Identification of novel PROM1 mutations responsible for autosomal recessive maculopathy with rod-cone dystrophy

PURPOSE

To characterize two patients with macular and rod-cone dystrophy and identify the genetic basis for disease.

METHOD

Ophthalmic examinations were performed for the family and the peripheral blood samples were collected for whole exome sequencing. The mutated sequences of PROM1 gene were cloned and expressed in cultured cell lines after transient transfection followed by analysis with confocal microscopy and bridge-PCR.

RESULT

We reported that two patients, brothers in a family, were diagnosed with macular and rod-cone dystrophy. Phenotypically, both patients experience progressive visual impairment and nyctalopia. The fundus examination showed macular and choroid dystrophy with pigment deposits in the macular region. Functionally, photoreceptor response to electrophysiological stimulation was significantly compromised with more severe decline in rods. Genetic analysis by whole exome sequencing revealed two novel compound heterogeneous point mutations in PROM1 gene that co-segregate with patients in an autosomal recessive manner. Specifically, the c.C1902G(p.Y634X) nonsense mutation results in a truncated, labile, and mislocalized protein, while the c.C1682+3A>G intronic mutation disrupts messenger RNA splicing.

CONCLUSION

Our findings have identified two novel deleterious mutations in PROM1 gene that are associated with hereditary macular and rod-cone dystrophy in human.