Splicing and multifactorial analysis of intronic BRCA1 and BRCA2 sequence variants identifies clinically significant splicing aberrations up to 12 nucleotides from the intron/exon boundary

HRness in Breast and Ovarian Cancers

Ovarian and breast cancers are currently defined by the main pathways involved in the tumorigenesis. The majority are carcinomas, originating from epithelial cells that are in constant division and subjected to cyclical variations of the estrogen stimulus during the female hormonal cycle, therefore being vulnerable to DNA damage. A portion of breast and ovarian carcinomas arises in the context of DNA repair defects, in which genetic instability is the backdrop for cancer initiation and progression. For these tumors, DNA repair deficiency is now increasingly recognized as a target for therapeutics. In hereditary breast/ovarian cancers (HBOC), tumors with BRCA1/2 mutations present an impairment of DNA repair by homologous recombination (HR). For many years, BRCA1/2 mutations were only screened on germline DNA, but now they are also searched at the tumor level to personalize treatment. The reason of the inactivation of this pathway remains uncertain for most cases, even in the presence of a HR-deficient signature. Evidence indicates that identifying the mechanism of HR inactivation should improve both genetic counseling and therapeutic response, since they can be useful as new biomarkers of response.

Characterization of splice-altering mutations in inherited predisposition to cancer

Mutations responsible for inherited disease may act by disrupting normal transcriptional splicing. Such mutations can be difficult to detect, and their effects difficult to characterize, because many lie deep within exons or introns where they may alter splice enhancers or silencers or introduce new splice acceptors or donors. Multiple mutation-specific and genome-wide approaches have been developed to evaluate these classes of mutations. We introduce a complementary experimental approach, cBROCA, which yields qualitative and quantitative assessments of the effects of genomic mutations on transcriptional splicing of tumor suppressor genes. cBROCA analysis is undertaken by deriving complementary DNA (cDNA) from puromycin-treated patient lymphoblasts, hybridizing the cDNA to the BROCA panel of tumor suppressor genes, and then multiplex sequencing to very high coverage. At each splice junction suggested by split sequencing reads, read depths of test and control samples are compared. Significant Z scores indicate altered transcripts, over and above naturally occurring minor transcripts, and comparisons of read depths indicate relative abundances of mutant and normal transcripts. BROCA analysis of genomic DNA suggested 120 rare mutations from 150 families with cancers of the breast, ovary, uterus, or colon, in >600 informative genotyped relatives. cBROCA analysis of their transcripts revealed a wide variety of consequences of abnormal splicing in tumor suppressor genes, including whole or partial exon skipping, exonification of intronic sequence, loss or gain of exonic and intronic splicing enhancers and silencers, complete intron retention, hypomorphic alleles, and combinations of these alterations. Combined with pedigree analysis, cBROCA sequencing contributes to understanding the clinical consequences of rare inherited mutations.

Comprehensive Assessment of BARD1 Messenger Ribonucleic Acid Splicing With Implications for Variant Classification

Introduction: Case-control analyses have shown BARD1 variants to be associated with up to >2-fold increase in risk of breast cancer, and potentially greater risk of triple negative breast cancer. BARD1 is included in several gene sequencing panels currently marketed for the prediction of risk of cancer, however there are no gene-specific guidelines for the classification of BARD1 variants. We present the most comprehensive assessment of BARD1 messenger RNA splicing, and demonstrate the application of these data for the classification of truncating and splice site variants according to American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines. Methods: Nanopore sequencing, short-read RNA-seq (whole transcriptome and targeted), and capillary electrophoresis analysis were performed by four laboratories to investigate alternative BARD1 splicing in blood, breast, and fimbriae/ovary related specimens from non-cancer affected tissues. Splicing data were also collated from published studies of nine different tissues. The impact of the findings for PVS1 annotation was assessed for truncating and splice site variants. Results: We identified 62 naturally occurring alternative spliced BARD1 splicing events, including 19 novel events found by next generation sequencing and/or reverse transcription PCR analysis performed for this study. Quantitative analysis showed that naturally occurring splicing events causing loss of clinically relevant domains or nonsense mediated decay can constitute up to 11.9% of overlapping natural junctions, suggesting that aberrant splicing can be tolerated up to this level. Nanopore sequencing of whole BARD1 transcripts characterized 16 alternative isoforms from healthy controls, revealing that the most complex transcripts combined only two alternative splicing events. Bioinformatic analysis of ClinVar submitted variants at or near BARD1 splice sites suggest that all consensus splice site variants in BARD1 should be considered likely pathogenic, with the possible exception of variants at the donor site of exon 5. Conclusions: No BARD1 candidate rescue transcripts were identified in this study, indicating that all premature translation-termination codons variants can be annotated as PVS1. Furthermore, our analysis suggests that all donor and acceptor (IVS+/-1,2) variants can be considered PVS1 or PVS1_strong, with the exception of variants targeting the exon 5 donor site, that we recommend considering as PVS1_moderate.

Large scale multifactorial likelihood quantitative analysis of BRCA1 and BRCA2 variants: An ENIGMA resource to support clinical variant classification

The multifactorial likelihood analysis method has demonstrated utility for quantitative assessment of variant pathogenicity for multiple cancer syndrome genes. Independent data types currently incorporated in the model for assessing BRCA1 and BRCA2 variants include clinically calibrated prior probability of pathogenicity based on variant location and bioinformatic prediction of variant effect, co-segregation, family cancer history profile, co-occurrence with a pathogenic variant in the same gene, breast tumor pathology, and case-control information. Research and clinical data for multifactorial likelihood analysis were collated for 1,395 BRCA1/2 predominantly intronic and missense variants, enabling classification based on posterior probability of pathogenicity for 734 variants: 447 variants were classified as (likely) benign, and 94 as (likely) pathogenic; and 248 classifications were new or considerably altered relative to ClinVar submissions. Classifications were compared with information not yet included in the likelihood model, and evidence strengths aligned to those recommended for ACMG/AMP classification codes. Altered mRNA splicing or function relative to known nonpathogenic variant controls were moderately to strongly predictive of variant pathogenicity. Variant absence in population datasets provided supporting evidence for variant pathogenicity. These findings have direct relevance for BRCA1 and BRCA2 variant evaluation, and justify the need for gene-specific calibration of evidence types used for variant classification.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5′ and 3′ consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Incorporation of semi-quantitative analysis of splicing alterations for the clinical interpretation of variants in BRCA1 and BRCA2 genes

BRCA1 and BRCA2 (BRCA1/2) genetic variants that disrupt messenger RNA splicing are commonly associated with increased risks of developing breast/ovarian cancer. The majority of splicing studies published to date rely on qualitative methodologies (i.e., Sanger sequencing), but it is necessary to incorporate semi-quantitative or quantitative approaches to accurately interpret the clinical significance of spliceogenic variants. Here, we characterize the splicing impact of 31 BRCA1/2 variants using semi-quantitative capillary electrophoresis of fluorescent amplicons (CE), Sanger sequencing and allele-specific assays. A total of 14 variants were found to disrupt splicing. Allelic-specific assays could be performed for BRCA1 c.302-1G>A and BRCA2 c.516+2T>A, c.1909+1G>A, c.8332-13T>G, c.8332-2A>G, c.8954-2A>T variants, showing a monoallelic contribution to full-length transcript expression that was concordant with semi-quantitative data. The splicing fraction of alternative and aberrant transcripts was also measured by CE, facilitating variant interpretation. Following Evidence-based Network for the Interpretation of Germline Mutant Alleles criteria, we successfully classified eight variants as pathogenic (Class 5), five variants as likely pathogenic (Class 4), and 14 variants as benign (Class 1). We also provide splicing data for four variants classified as uncertain (Class 3), which produced a "leaky" splicing effect or introduced a missense change in the protein sequence, that will require further assessment to determine their clinical significance.

BRCA1-No Matter How You Splice It

BRCA1 (breast cancer 1, early onset), a well-known breast cancer susceptibility gene, is a highly alternatively spliced gene. BRCA1 alternative splicing may serve as an alternative regulatory mechanism for the inactivation of the BRCA1 gene in both hereditary and sporadic breast cancers, and other BRCA1-associated cancers. The alternative transcripts of BRCA1 can mimic known functions, possess unique functions compared with the full-length BRCA1 transcript, and in some cases, appear to function in opposition to full-length BRCA1 In this review, we will summarize the functional "naturally occurring" alternative splicing transcripts of BRCA1 and then discuss the latest next-generation sequencing-based detection methods and techniques to detect alternative BRCA1 splicing patterns and their potential use in cancer diagnosis, prognosis, and therapy.

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.

Usefulness and Limitations of Comprehensive Characterization of mRNA Splicing Profiles in the Definition of the Clinical Relevance of BRCA1/2 Variants of Uncertain Significance

Highly penetrant variants of BRCA1/2 genes are involved in hereditary predisposition to breast and ovarian cancer. The detection of pathogenic BRCA variants has a considerable clinical impact, allowing appropriate cancer-risk management. However, a major drawback is represented by the identification of variants of uncertain significance (VUS). Many VUS potentially affect mRNA splicing, making transcript analysis an essential step for the definition of their pathogenicity. Here, we characterize the impact on splicing of ten BRCA1/2 variants. Aberrant splicing patterns were demonstrated for eight variants whose alternative transcripts were fully characterized. Different events were observed, including exon skipping, intron retention, and usage of de novo and cryptic splice sites. Transcripts with premature stop codons or in-frame loss of functionally important residues were generated. Partial/complete splicing effect and quantitative contribution of different isoforms were assessed, leading to variant classification according to Evidence-based Network for the Interpretation of Mutant Alleles (ENIGMA) consortium guidelines. Two variants could be classified as pathogenic and two as likely benign, while due to a partial splicing effect, six variants remained of uncertain significance. The association with an undefined tumor risk justifies caution in recommending aggressive risk-reduction treatments, but prevents the possibility of receiving personalized therapies with potential beneficial effect. This indicates the need for applying additional approaches for the analysis of variants resistant to classification by gene transcript analyses.

Targeted RNA-seq successfully identifies normal and pathogenic splicing events in breast/ovarian cancer susceptibility and Lynch syndrome genes

A subset of genetic variants found through screening of patients with hereditary breast and ovarian cancer syndrome (HBOC) and Lynch syndrome impact RNA splicing. Through target enrichment of the transcriptome, it is possible to perform deep-sequencing and to identify the different and even rare mRNA isoforms. A targeted RNA-seq approach was used to analyse the naturally-occurring splicing events for a panel of 8 breast and/or ovarian cancer susceptibility genes (BRCA1, BRCA2, RAD51C, RAD51D, PTEN, STK11, CDH1, TP53), 3 Lynch syndrome genes (MLH1, MSH2, MSH6) and the fanconi anaemia SLX4 gene, in which monoallelic mutations were found in non-BRCA families. For BRCA1, BRCA2, RAD51C and RAD51D the results were validated by capillary electrophoresis and were compared to a non-targeted RNA-seq approach. We also compared splicing events from lymphoblastoid cell-lines with those from breast and ovarian fimbriae tissues. The potential of targeted RNA-seq to detect pathogenic changes in RNA-splicing was validated by the inclusion of samples with previously well characterized BRCA1/2 genetic variants. In our study, we update the catalogue of normal splicing events for BRCA1/2, provide an extensive catalogue of normal RAD51C and RAD51D alternative splicing, and list splicing events found for eight other genes. Additionally, we show that our approach allowed the identification of aberrant splicing events due to the presence of BRCA1/2 genetic variants and distinguished between complete and partial splicing events. In conclusion, targeted-RNA-seq can be very useful to classify variants based on their putative pathogenic impact on splicing.

A comparison of cosegregation analysis methods for the clinical setting

Quantitative cosegregation analysis can help evaluate the pathogenicity of genetic variants. However, genetics professionals without statistical training often use simple methods, reporting only qualitative findings. We evaluate the potential utility of quantitative cosegregation in the clinical setting by comparing three methods. One thousand pedigrees each were simulated for benign and pathogenic variants in BRCA1 and MLH1 using United States historical demographic data to produce pedigrees similar to those seen in the clinic. These pedigrees were analyzed using two robust methods, full likelihood Bayes factors (FLB) and cosegregation likelihood ratios (CSLR), and a simpler method, counting meioses. Both FLB and CSLR outperform counting meioses when dealing with pathogenic variants, though counting meioses is not far behind. For benign variants, FLB and CSLR greatly outperform as counting meioses is unable to generate evidence for benign variants. Comparing FLB and CSLR, we find that the two methods perform similarly, indicating that quantitative results from either of these methods could be combined in multifactorial calculations. Combining quantitative information will be important as isolated use of cosegregation in single families will yield classification for less than 1% of variants. To encourage wider use of robust cosegregation analysis, we present a website ( http://www.analyze.myvariant.org ) which implements the CSLR, FLB, and Counting Meioses methods for ATM, BRCA1, BRCA2, CHEK2, MEN1, MLH1, MSH2, MSH6, and PMS2. We also present an R package, CoSeg, which performs the CSLR analysis on any gene with user supplied parameters. Future variant classification guidelines should allow nuanced inclusion of cosegregation evidence against pathogenicity.

Prevalence of germline pathogenic BRCA1/2 variants in sequential epithelial ovarian cancer cases

Introduction

Poly(ADP-ribose) polymerase inhibitors significantly improve progression-free survival in platinum-sensitive high-grade serous and endometrioid ovarian carcinoma, with greatest benefits observed in women with a pathogenic BRCA1/2 variant. Consequently, the demand for germline BRCA1/2 testing in ovarian cancer has increased substantially, leading to the screening of unselected populations of patients. We aimed to determine the prevalence of pathogenic germline BRCA1/2 variants in women diagnosed with epithelial ovarian cancer, categorised according to the established risk factors for hereditary breast and ovarian cancer syndrome and the Manchester BRCA Score, to inform risk stratification.

Methods

A cohort of sequential epithelial ovarian cancer cases recruited between June 2013 and September 2018 underwent germline BRCA1/2 testing by next-generation sequencing and multiplex ligation-dependent probe amplification.

Results

Five hundred and fifty-seven patients were screened. Of these, 18% had inherited a pathogenic BRCA1/2 variant. The prevalence of pathogenic BRCA1/2 variants was >10% in women diagnosed with ovarian cancer earlier than 60 years of age (21%) and those diagnosed later than 60 years of age with a family history of breast and/or ovarian cancer (17%) or a medical history of breast cancer (34%). The prevalence of pathogenic BRCA1/2 variants was also >10% in women with a Manchester BRCA Score of ≥15 points (14%).

Discussion

Our study suggests that age at diagnosis, family history of breast and/or ovarian cancer, medical history of breast cancer or a Manchester BRCA Score of ≥15 points are associated with a >10% prevalence of germline pathogenic BRCA1/2 variants in epithelial ovarian cancer.

Genetic variants and haplotype combination in the bovine CRTC3 affected conformation traits in two Chinese native cattle breeds (Bos Taurus)

CREB-regulated transcription coactivator 3 (CRTC3) plays an extensive role in glucose and lipid metabolism. This study investigated the genetic variation and haplotype combination in CRTC3 and verified their contribution to bovine growth traits. Firstly, investigated the mRNA expression of CRTC3 in adult Qinchuan cattle and evaluated the effects that genetic variation of CRTC3 had on conformation and carcass traits in two Chinese cattle breeds (Qinchuan and Jiaxian). Four SNPs (single nucleotide polymorphisms) were identified including two in introns (SNP1: g.62652 A > G and SNP4: g.91297C > T) and two in exons (SNP2 g.62730C > T and SNP3: g.66478G > C). The association and haplotype combination results showed that there was an association with some growth and carcass traits(P < 0.05). Individuals with haplotype combination H1H1 (-AACCCCTT-) were associated with a conformation of a larger framed animal and an animal that produced a larger loin area. Variations in the CRTC3 genes and the haplotype combination H1H1 may be considered as molecular markers for carcass traits that are associated with more lean meat yield for use in cattle breeding programs in China.

Characterization of spliceogenic variants located in regions linked to high levels of alternative splicing: BRCA2 c.7976+5G > T as a case study

Many BRCA1 and BRCA2 (BRCA1/2) genetic variants have been studied at mRNA level and linked to hereditary breast and ovarian cancer due to splicing alteration. In silico tools are reliable when assessing variants located in consensus splice sites, but we may identify variants in complex genomic contexts for which bioinformatics is not precise enough. In this study, we characterize BRCA2 c.7976 + 5G > T variant located in intron 17 which has an atypical donor site (GC). This variant was identified in three unrelated Spanish families and we have detected exon 17 skipping as the predominant transcript occurring in carriers. We have also detected several isoforms (Δ16-18, Δ17,18, Δ18, and ▼17q²²⁴ ) at different expression levels among carriers and controls. This study remarks the challenge of interpreting genetic variants when multiple alternative isoforms are present, and that caution must be taken when using in silico tools to identify potential spliceogenic variants located in GC-AG introns.

Hereditary breast and ovarian cancer in Andalusian families: a genetic population study

Background

The BRCA1/2 mutation profile varies in Spain according to the geographical area studied. The mutational profile of BRCA1/2 in families at risk for hereditary breast and ovarian cancer has not so far been reported in Andalusia (southern Spain).

Methods

We analysed BRCA1/2 germline mutations in 562 high-risk cases with breast and/or ovarian cancer from Andalusian families from 2010 to 2015.

Results

Among the 562 cases, 120 (21.4%) carried a germline pathogenic mutation in BRCA1/2; 50 in BRCA1 (41.7%) and 70 in BRCA2 (58.3%). We detected 67 distinct mutations (29 in BRCA1 and 38 in BRCA2), of which 3 in BRCA1 (c.845C > A, c.1222_1223delAC, c.2527delA) and 5 in BRCA2 (c.293 T > G, c.5558_5559delGT, c.6034delT, c.6650_6654delAAGAT, c.6652delG) had not been previously described. The most frequent mutations in BRCA1 were c.5078_5080delCTG (10%) and c.5123C > A (10%), and in BRCA2 they were c.9018C > A (14%) and c.5720_5723delCTCT (8%). We identified 5 variants of unknown significance (VUS), all in BRCA2 (c.5836 T > C, c.6323G > T, c.9501 + 3A > T, c.8022_8030delGATAATGGA, c.10186A > C). We detected 76 polymorphisms (31 in BRCA1, 45 in BRCA2) not associated with breast cancer risk.

Conclusions

This is the first study reporting the mutational profile of BRCA1/2 in Andalusia. We identified 21.4% of patients harbouring BRCA1/2 mutations, 58.3% of them in BRCA2. We also characterized the clinical data, mutational profile, VUS and haplotype profile.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4537-9) contains supplementary material, which is available to authorized users.

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

Genetic testing of BRCA1 and BRCA2 identifies a large number of variants of uncertain clinical significance whose functional and clinical interpretations pose a challenge for genetic counseling. Interestingly, a relevant fraction of DNA variants can disrupt the splicing process in cancer susceptibility genes. We have tested more than 200 variants throughout 19 BRCA2 exons mostly by minigene assays, 54% of which displayed aberrant splicing, thus confirming the utility of this assay to check genetic variants in the absence of patient RNA. Our goal was to investigate BRCA2 exon 16 with a view to characterizing spliceogenic variants recorded at the mutational databases. Seventy-two different BIC and UMD variants were analyzed with NNSplice and Human Splicing Finder, 12 of which were selected because they were predicted to disrupt essential splice motifs: canonical splice sites (ss; eight variants) and exonic/intronic splicing enhancers (four variants). These 12 candidate variants were introduced into the BRCA2 minigene with seven exons (14–20) by site-directed mutagenesis and then transfected into MCF-7 cells. Seven variants (six intronic and one missense) induced complete abnormal splicing patterns: c.7618-2A>T, c.7618-2A>G, c.7618-1G>C, c.7618-1G>A, c.7805G>C, c.7805+1G>A, and c.7805+3A>C, as well as a partial anomalous outcome by c.7802A>G. They generated at least 10 different transcripts: Δ16p₄₄ (alternative 3’ss 44-nt downstream; acceptor variants), Δ16 (exon 16-skipping; donor variants), Δ16p₅₅ (alternative 3’ss 55-nt downstream), Δ16q₄ (alternative 5’ss 4-nt upstream), Δ16q₁₀₀ (alternative 5’ss 4-nt upstream), ▾16q²⁰ (alternative 5’ss 20-nt downstream), as well as minor (Δ16p₉₃ and Δ16,17p₆₉) and uncharacterized transcripts of 893 and 954 nucleotides. Isoforms Δ16p₄₄, Δ16, Δ16p₅₅, Δ16q₄, Δ16q₁₀₀, and ▾16q²⁰ introduced premature termination codons which presumably inactivate BRCA2. According to the guidelines the American College of Medical Genetics and Genomics these eight variants could be classified as pathogenic or likely pathogenic whereas the Evidence-based Network for the Interpretation of Germline Mutant Alleles rules suggested seven class 4 and one class 3 variants. In conclusion, our study highlights the relevance of splicing functional assays by hybrid minigenes for the clinical classification of genetic variations. Hence, we provide new data about spliceogenic variants of BRCA2 exon 16 that are directly correlated with breast cancer susceptibility.

Investigation of Experimental Factors That Underlie BRCA1/2 mRNA Isoform Expression Variation: Recommendations for Utilizing Targeted RNA Sequencing to Evaluate Potential Spliceogenic Variants

PCR-based RNA splicing assays are commonly used in diagnostic and research settings to assess the potential effects of variants of uncertain clinical significance in BRCA1 and BRCA2. The Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium completed a multicentre investigation to evaluate differences in assay design and the integrity of published data, raising a number of methodological questions associated with cell culture conditions and PCR-based protocols. We utilized targeted RNA-seq to re-assess BRCA1 and BRCA2 mRNA isoform expression patterns in lymphoblastoid cell lines (LCLs) previously used in the multicentre ENIGMA study. Capture of the targeted cDNA sequences was carried out using 34 BRCA1 and 28 BRCA2 oligonucleotides from the Illumina Truseq Targeted RNA Expression platform. Our results show that targeted RNA-seq analysis of LCLs overcomes many of the methodology limitations associated with PCR-based assays leading us to make the following observations and recommendations: (1) technical replicates (n > 2) of variant carriers to capture methodology induced variability associated with RNA-seq assays, (2) LCLs can undergo multiple freeze/thaw cycles and can be cultured up to 2 weeks without noticeably influencing isoform expression levels, (3) nonsense-mediated decay inhibitors are essential prior to splicing assays for comprehensive mRNA isoform detection, (4) quantitative assessment of exon:exon junction levels across BRCA1 and BRCA2 can help distinguish between normal and aberrant isoform expression patterns. Experimentally derived recommendations from this study will facilitate the application of targeted RNA-seq platforms for the quantitation of BRCA1 and BRCA2 mRNA aberrations associated with sequence variants of uncertain clinical significance.

Splice variants of the endonucleases XPF and XPG contain residual DNA repair capabilities and could be a valuable tool for personalized medicine

The two endonucleases XPF and XPG are essentially involved in nucleotide excision repair (NER) and interstrand crosslink (ICL) repair. Defects in these two proteins result in severe diseases like xeroderma pigmentosum (XP). We applied our newly CRISPR/Cas9 generated human XPF knockout cell line with complete loss of XPF and primary fibroblasts from an XP-G patient (XP20BE) to analyze until now uncharacterized spontaneous mRNA splice variants of these two endonucleases. Functional analyses of these variants were performed using luciferase-based reporter gene assays. Two XPF and XPG splice variants with residual repair capabilities in NER, as well as ICL repair could be identified. Almost all variants are severely C-terminally truncated and lack important protein-protein interaction domains. Interestingly, XPF-202, differing to XPF-003 in the first 12 amino acids only, had no repair capability at all, suggesting an important role of this region during DNA repair, potentially concerning protein-protein interaction. We also identified splice variants of XPF and XPG exerting inhibitory effects on NER. Moreover, we showed that the XPF and XPG splice variants presented with different inter-individual expression patterns in healthy donors, as well as in various tissues. With regard to their residual repair capability and dominant-negative effects, functionally relevant spontaneous XPF and XPG splice variants present promising prognostic marker candidates for individual cancer risk, disease outcome, or therapeutic success. This merits further investigations, large association studies, and translational research within clinical trials in the future.

Adding In Silico Assessment of Potential Splice Aberration to the Integrated Evaluation of BRCA Gene Unclassified Variants

Clinical mutation screening of the cancer susceptibility genes BRCA1 and BRCA2 generates many unclassified variants (UVs). Most of these UVs are either rare missense substitutions or nucleotide substitutions near the splice junctions of the protein coding exons. Previously, we developed a quantitative method for evaluation of BRCA gene UVs—the “integrated evaluation”—that combines a sequence analysis‐based prior probability of pathogenicity with patient and/or tumor observational data to arrive at a posterior probability of pathogenicity. One limitation of the sequence analysis‐based prior has been that it evaluates UVs from the perspective of missense substitution severity but not probability to disrupt normal mRNA splicing. Here, we calibrated output from the splice‐site fitness program MaxEntScan to generate spliceogenicity‐based prior probabilities of pathogenicity for BRCA gene variants; these range from 0.97 for variants with high probability to damage a donor or acceptor to 0.02 for exonic variants that do not impact a splice junction and are unlikely to create a de novo donor. We created a database http://priors.hci.utah.edu/PRIORS/ that provides the combined missense substitution severity and spliceogenicity‐based probability of pathogenicity for BRCA gene single‐nucleotide substitutions. We also updated the BRCA gene Ex‐UV LOVD, available at http://hci‐exlovd.hci.utah.edu, with 77 re‐evaluable variants.

Naturally occurring BRCA2 alternative mRNA splicing events in clinically relevant samples

BACKGROUND

BRCA1 and BRCA2 are the two principal tumour suppressor genes associated with inherited high risk of breast and ovarian cancer. Genetic testing of BRCA1/2 will often reveal one or more sequence variants of uncertain clinical significance, some of which may affect normal splicing patterns and thereby disrupt gene function. mRNA analyses are therefore among the tests used to interpret the clinical significance of some genetic variants. However, these could be confounded by the appearance of naturally occurring alternative transcripts unrelated to germline sequence variation or defects in gene function. To understand which novel splicing events are associated with splicing mutations and which are part of the normal BRCA2 splicing repertoire, a study was undertaken by members of the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium to characterise the spectrum of naturally occurring BRCA2 mRNA alternate-splicing events.

METHODS

mRNA was prepared from several blood and breast tissue-derived cells and cell lines by contributing ENIGMA laboratories. cDNA representing BRCA2 alternate splice sites was amplified and visualised using capillary or agarose gel electrophoresis, followed by sequencing.

RESULTS

We demonstrate the existence of 24 different BRCA2 mRNA alternate-splicing events in lymphoblastoid cell lines and both breast cancer and non-cancerous breast cell lines.

CONCLUSIONS

These naturally occurring alternate-splicing events contribute to the array of cDNA fragments that may be seen in assays for mutation-associated splicing defects. Caution must be observed in assigning alternate-splicing events to potential splicing mutations.

Functional classification of BRCA2 DNA variants by splicing assays in a large minigene with 9 exons

Numerous pathogenic DNA variants impair the splicing mechanism in human genetic diseases. Minigenes are optimal approaches to test variants under the splicing viewpoint without the need of patient samples. We aimed to design a robust minigene construct of the breast cancer gene BRCA2 in order to investigate the impact of variants on splicing. BRCA2 exons 19-27 (MGBR2_ex19-27) were cloned in the new vector pSAD. It produced a large transcript of the expected size (2,174 nucleotides) and exon structure (V1-ex19-27-V2). Splicing assays showed that 18 (17 splice-site and 1 silencer variants) out of 40 candidate DNA variants induced aberrant patterns. Twenty-four anomalous transcripts were accurately detected by fluorescent-RT-PCR that were generated by exon-skipping, alternative site usage, and intron-retention events. Fourteen variants induced major anomalies and were predicted to disrupt protein function so they could be classified as pathogenic. Furthermore, minigene mimicked previously reported patient RNA outcomes of seven variants supporting the reproducibility of minigene assays. Therefore, a relevant fraction of variants are involved in breast cancer through splicing alterations. MGBR2_ex19-27 is the largest reported BRCA2 minigene and constitutes a valuable tool for the functional and clinical classification of sequence variations.

The alternative splicing side of cancer

Alternative splicing emerges as a potent and pervasive mechanism of gene expression regulation that expands the coding capacity of the genome and forms an intermediate layer of regulation between transcriptional and post-translational networks. Indeed, alternative splicing occupies a pivotal position in developmental programs and in the cell response to external and internal stimuli. Not surprisingly, therefore, its deregulation frequently leads to human disease. In this review we provide an updated overview of the impact of alternative splicing on tumorigenesis. Moreover, we discuss the intricacy of the reciprocal interactions between alternative splicing programs and signal transduction pathways, which appear to be crucially linked to cancer progression in response to the tumor microenvironment. Finally, we focus on the recently described interplay between splicing and chromatin organization which is expected to shed new lights into gene expression regulation in normal and cancer cells.

BRCA1/2 sequence variants of uncertain significance: a primer for providers to assist in discussions and in medical management

INTRODUCTION

DNA variants of uncertain significance (VUS) are common outcomes of clinical genetic testing for susceptibility to cancer. A statistically rigorous model that provides a pathogenicity score for each variant has been developed to aid in the clinical management of patients undergoing genetic testing.

METHODS

The information in this article is derived from multiple publications on VUS in BRCA genes, distilled for communicating with clinicians who may encounter VUS in their practice.

RESULTS

The posterior probability scores for BRCA1 or BRCA2 VUS, calculated from a multifactorial likelihood model, are explained, and links for looking up specific VUS are provided. The International Agency on Cancer Research (IARC) of the World Health Organization has proposed a simple five-tier system for clinical management that is not widely known to clinicians. Classes 1 and 2 in this system are managed as neutral variants, classes 4 and 5 are managed as pathogenic variants, and class 3 variants still have insufficient evidence to move to either end of this scale and, thus, cannot be used in medical management.

CONCLUSIONS

Development of models that integrate multiple independent lines of evidence has allowed classification of a growing number of VUS in the BRCA1 and BRCA2 genes. The pathogenicity score that is generated by this model maps to the IARC system for clinical management, which will assist clinicians in the medical management of those patients who obtain a VUS result upon testing.

Comparative in vitro and in silico analyses of variants in splicing regions of BRCA1 and BRCA2 genes and characterization of novel pathogenic mutations

Several unclassified variants (UVs) have been identified in splicing regions of disease-associated genes and their characterization as pathogenic mutations or benign polymorphisms is crucial for the understanding of their role in disease development. In this study, 24 UVs located at BRCA1 and BRCA2 splice sites were characterized by transcripts analysis. These results were used to evaluate the ability of nine bioinformatics programs in predicting genetic variants causing aberrant splicing (spliceogenic variants) and the nature of aberrant transcripts. Eleven variants in BRCA1 and 8 in BRCA2, including 8 not previously characterized at transcript level, were ascertained to affect mRNA splicing. Of these, 16 led to the synthesis of aberrant transcripts containing premature termination codons (PTCs), 2 to the up-regulation of naturally occurring alternative transcripts containing PTCs, and one to an in-frame deletion within the region coding for the DNA binding domain of BRCA2, causing the loss of the ability to bind the partner protein DSS1 and ssDNA. For each computational program, we evaluated the rate of non-informative analyses, i.e. those that did not recognize the natural splice sites in the wild-type sequence, and the rate of false positive predictions, i.e., variants incorrectly classified as spliceogenic, as a measure of their specificity, under conditions setting sensitivity of predictions to 100%. The programs that performed better were Human Splicing Finder and Automated Splice Site Analyses, both exhibiting 100% informativeness and specificity. For 10 mutations the activation of cryptic splice sites was observed, but we were unable to derive simple criteria to select, among the different cryptic sites predicted by the bioinformatics analyses, those actually used. Consistent with previous reports, our study provides evidences that in silico tools can be used for selecting splice site variants for in vitro analyses. However, the latter remain mandatory for the characterization of the nature of aberrant transcripts.

A multifactorial likelihood model for MMR gene variant classification incorporating probabilities based on sequence bioinformatics and tumor characteristics: a report from the Colon Cancer Family Registry

Mismatch repair (MMR) gene sequence variants of uncertain clinical significance are often identified in suspected Lynch syndrome families, and this constitutes a challenge for both researchers and clinicians. Multifactorial likelihood model approaches provide a quantitative measure of MMR variant pathogenicity, but first require input of likelihood ratios (LRs) for different MMR variation-associated characteristics from appropriate, well-characterized reference datasets. Microsatellite instability (MSI) and somatic BRAF tumor data for unselected colorectal cancer probands of known pathogenic variant status were used to derive LRs for tumor characteristics using the Colon Cancer Family Registry (CFR) resource. These tumor LRs were combined with variant segregation within families, and estimates of prior probability of pathogenicity based on sequence conservation and position, to analyze 44 unclassified variants identified initially in Australasian Colon CFR families. In addition, in vitro splicing analyses were conducted on the subset of variants based on bioinformatic splicing predictions. The LR in favor of pathogenicity was estimated to be ~12-fold for a colorectal tumor with a BRAF mutation-negative MSI-H phenotype. For 31 of the 44 variants, the posterior probabilities of pathogenicity were such that altered clinical management would be indicated. Our findings provide a working multifactorial likelihood model for classification that carefully considers mode of ascertainment for gene testing.

BRCA1 R1699Q variant displaying ambiguous functional abrogation confers intermediate breast and ovarian cancer risk

BACKGROUND

Clinical classification of rare sequence changes identified in the breast cancer susceptibility genes BRCA1 and BRCA2 is essential for appropriate genetic counselling of individuals carrying these variants. We previously showed that variant BRCA1 c.5096G>A p.Arg1699Gln in the BRCA1 transcriptional transactivation domain demonstrated equivocal results from a series of functional assays, and proposed that this variant may confer low to moderate risk of cancer.

METHODS

Measures of genetic risk (report of family history, segregation) were assessed for 68 BRCA1 c.5096G>A p.Arg1699Gln (R1699Q) families recruited through family cancer clinics, comparing results with 34 families carrying the previously classified pathogenic BRCA1 c.5095C>T p.Arg1699Trp (R1699W) mutation at the same residue, and to 243 breast cancer families with no BRCA1 pathogenic mutation (BRCA-X).

RESULTS

Comparison of BRCA1 carrier prediction scores of probands using the BOADICEA risk prediction tool revealed that BRCA1 c.5096G>A p.Arg1699Gln variant carriers had family histories that were less 'BRCA1-like' than BRCA1 c.5095C>T p.Arg1699Trp mutation carriers (p<0.00001), but more 'BRCA1-like' than BRCA-X families (p=0.0004). Further, modified segregation analysis of the subset of 30 families with additional genotyping showed that BRCA1 c.5096G >A p.Arg1699Gln had reduced penetrance compared with the average truncating BRCA1 mutation penetrance (p=0.0002), with estimated cumulative risks to age 70 of breast or ovarian cancer of 24%.

CONCLUSIONS

Our results provide substantial evidence that the BRCA1 c.5096G>A p.Arg1699Gln (R1699Q) variant, demonstrating ambiguous functional deficiency across multiple assays, is associated with intermediate risk of breast and ovarian cancer, highlighting challenges for risk modelling and clinical management of patients of this and other potential moderate-risk variants.

Analysis of 30 putative BRCA1 splicing mutations in hereditary breast and ovarian cancer families identifies exonic splice site mutations that escape in silico prediction

Screening for pathogenic mutations in breast and ovarian cancer genes such as BRCA1/2, CHEK2 and RAD51C is common practice for individuals from high-risk families. However, test results may be ambiguous due to the presence of unclassified variants (UCV) in the concurrent absence of clearly cancer-predisposing mutations. Especially the presence of intronic or exonic variants within these genes that possibly affect proper pre-mRNA processing poses a challenge as their functional implications are not immediately apparent. Therefore, it appears necessary to characterize potential splicing UCV and to develop appropriate classification tools. We investigated 30 distinct BRCA1 variants, both intronic and exonic, regarding their spliceogenic potential by commonly used in silico prediction algorithms (HSF, MaxEntScan) along with in vitro transcript analyses. A total of 25 variants were identified spliceogenic, either causing/enhancing exon skipping or activation of cryptic splice sites, or both. Except from a single intronic variant causing minor effects on BRCA1 pre-mRNA processing in our analyses, 23 out of 24 intronic variants were correctly predicted by MaxEntScan, while HSF was less accurate in this cohort. Among the 6 exonic variants analyzed, 4 severely impair correct pre-mRNA processing, while the remaining two have partial effects. In contrast to the intronic alterations investigated, only half of the spliceogenic exonic variants were correctly predicted by HSF and/or MaxEntScan. These data support the idea that exonic splicing mutations are commonly disease-causing and concurrently prone to escape in silico prediction, hence necessitating experimental in vitro splicing analysis.

A guide for functional analysis of BRCA1 variants of uncertain significance

Germline mutations in the tumor suppressor gene BRCA1 confer an estimated lifetime risk of 56-80% for breast cancer and 15-60% for ovarian cancer. Since the mid 1990s when BRCA1 was identified, genetic testing has revealed over 1,500 unique germline variants. However, for a significant number of these variants, the effect on protein function is unknown making it difficult to infer the consequences on risks of breast and ovarian cancers. Thus, many individuals undergoing genetic testing for BRCA1 mutations receive test results reporting a variant of uncertain clinical significance (VUS), leading to issues in risk assessment, counseling, and preventive care. Here, we describe functional assays for BRCA1 to directly or indirectly assess the impact of a variant on protein conformation or function and how these results can be used to complement genetic data to classify a VUS as to its clinical significance. Importantly, these methods may provide a framework for genome-wide pathogenicity assignment.

Comprehensive splicing functional analysis of DNA variants of the BRCA2 gene by hybrid minigenes

Introduction

The underlying pathogenic mechanism of a large fraction of DNA variants of disease-causing genes is the disruption of the splicing process. We aimed to investigate the effect on splicing of the BRCA2 variants c.8488-1G > A (exon 20) and c.9026_9030del (exon 23), as well as 41 BRCA2 variants reported in the Breast Cancer Information Core (BIC) mutation database.

Methods

DNA variants were analyzed with the splicing prediction programs NNSPLICE and Human Splicing Finder. Functional analyses of candidate variants were performed by lymphocyte RT-PCR and/or hybrid minigene assays. Forty-one BIC variants of exons 19, 20, 23 and 24 were bioinformatically selected and generated by PCR-mutagenesis of the wild type minigenes.

Results

Lymphocyte RT-PCR of c.8488-1G > A showed intron 19 retention and a 12-nucleotide deletion in exon 20, whereas c.9026_9030del did not show any splicing anomaly. Minigene analysis of c.8488-1G > A displayed the aforementioned aberrant isoforms but also exon 20 skipping. We further evaluated the splicing outcomes of 41 variants of four BRCA2 exons by minigene analysis. Eighteen variants presented splicing aberrations. Most variants (78.9%) disrupted the natural splice sites, whereas four altered putative enhancers/silencers and had a weak effect. Fluorescent RT-PCR of minigenes accurately detected 14 RNA isoforms generated by cryptic site usage, exon skipping and intron retention events. Fourteen variants showed total splicing disruptions and were predicted to truncate or eliminate essential domains of BRCA2.

Conclusions

A relevant proportion of BRCA2 variants are correlated with splicing disruptions, indicating that RNA analysis is a valuable tool to assess the pathogenicity of a particular DNA change. The minigene system is a straightforward and robust approach to detect variants with an impact on splicing and contributes to a better knowledge of this gene expression step.

Three novel BRCA1/BRCA2 mutations in breast/ovarian cancer families in Croatia

BRCA1 and BRCA2 genes from 167 candidates (145 families) were scanned for mutations. We identified 14 pathogenic point mutations in 17 candidates, 9 in BRCA1 and 5 in BRCA2. Of those, 11 have been previously described and 3 were novel (c.5335C>T in BRCA1 and c.4139_4140dupTT and c.8175G>A in BRCA2). No large deletions or duplications involving BRCA1 and BRCA2 genes were identified. No founder mutations were detected for the Croatian population. Croatia shares most of the mutations with neighboring Slovenia and also with Germany, Austria and Poland. Two common sequence variants in BRCA1, c.2077G>A and c.4956G>A, were found more frequently in mutation carriers compared to healthy controls. No difference in BRCA2 variants was detected between the groups. Haplotype inference showed no difference in haplotype distributions between deleterious mutation carriers and non-carriers in neither BRCA1 nor BRCA2. In silico analyses identified one BRCA1 sequence variant (c.4039A>G) and two BRCA2 variants (c.5986G>A and c.6884G>C) as harmful with high probability, and inconclusive results were obtained for our novel BRCA2 variant c.3864_3866delTAA. Combination of QMPSF and HRMA methods provides high detection rate and complete coverage of BRCA1/2 genes. Benefit of BRCA1/2 mutation testing is clear, since we detected mutations in young unaffected women, who will be closely monitored for breast and ovarian cancer.

Detection of exon skipping events in BRCA1 RNA using MLPA kit P002

A rapid and easy method to screen for aberrant cDNA would be a very useful diagnostic tool in genetics since a fraction of the DNA variants found affect RNA splicing. The currently used RT-PCR methods require new primer combinations to study each variant that might affect splicing. Since MLPA is routinely used to detect large genomic deletions and successfully detected exon skipping events in Duchenne muscular dystrophy in cDNA, we performed a pilot study to evaluate its value for BRCA1 cDNA. The effect of puromycin, DNase I and two different DNA cleaning protocols were tested in the RNA analysis of lymphocyte cultures. We used two samples from unrelated families with two different BRCA1 exon deletion events, two healthy unrelated controls and six samples from hereditary breast/ovarian cancer syndrome (HBOC) patients without BRCA1/2 mutations. Using RNA treated with DNase I and cleaned in a column system from puromycin-treated fractions, we were able to identify the two BRCA1 deletions. Additional HBOC patients did not show additional splice events. However, we were not able to get reproducible results. Therefore, the cDNA-MLPA technique using kit BRCA1 P002 is in our hands currently not reliable enough for routine RNA analysis and needs further optimization.

A review of a multifactorial probability-based model for classification of BRCA1 and BRCA2 variants of uncertain significance (VUS)

Clinical mutation screening of the BRCA1 and BRCA2 genes for the presence of germline inactivating mutations is used to identify individuals at elevated risk of breast and ovarian cancer. Variants identified during screening are usually classified as pathogenic (increased risk of cancer) or not pathogenic (no increased risk of cancer). However, a significant proportion of genetic tests yields variants of uncertain significance (VUS) that have undefined risk of cancer. Individuals carrying these VUS cannot benefit from individualized cancer risk assessment. Recently, a quantitative "posterior probability model" for assessing the clinical relevance of VUS in BRCA1 or BRCA2, which integrates multiple forms of genetic evidence has been developed. Here, we provide a detailed review of this model. We describe the components of the model and explain how these can be combined to calculate a posterior probability of pathogenicity for each VUS. We explain how the model can be applied to public data and provide tables that list the VUS that have been classified as not pathogenic or pathogenic using this method. While we use BRCA1 and BRCA2 VUS as examples, the method can be used as a framework for classification of the pathogenicity of VUS in other cancer genes.

ENIGMA--evidence-based network for the interpretation of germline mutant alleles: an international initiative to evaluate risk and clinical significance associated with sequence variation in BRCA1 and BRCA2 genes

As genetic testing for predisposition to human diseases has become an increasingly common practice in medicine, the need for clear interpretation of the test results is apparent. However, for many disease genes, including the breast cancer susceptibility genes BRCA1 and BRCA2, a significant fraction of tests results in the detection of a genetic variant for which disease association is not known. The finding of an "unclassified" variant (UV)/variant of uncertain significance (VUS) complicates genetic test reporting and counseling. As these variants are individually rare, a large collaboration of researchers and clinicians will facilitate studies to assess their association with cancer predisposition. It was with this in mind that the ENIGMA consortium (www.enigmaconsortium.org) was initiated in 2009. The membership is both international and interdisciplinary, and currently includes more than 100 research scientists and clinicians from 19 countries. Within ENIGMA, there are presently six working groups focused on the following topics: analysis, clinical, database, functional, tumor histopathology, and mRNA splicing. ENIGMA provides a mechanism to pool resources, exchange methods and data, and coordinately develop and apply algorithms for classification of variants in BRCA1 and BRCA2. It is envisaged that the research and clinical application of models developed by ENIGMA will be relevant to the interpretation of sequence variants in other disease genes.

Contribution of bioinformatics predictions and functional splicing assays to the interpretation of unclassified variants of the BRCA genes

A large fraction of sequence variants of unknown significance (VUS) of the breast and ovarian cancer susceptibility genes BRCA1 and BRCA2 may induce splicing defects. We analyzed 53 VUSs of BRCA1 or BRCA2, detected in consecutive molecular screenings, by using five splicing prediction programs, and we classified them into two groups according to the strength of the predictions. In parallel, we tested them by using functional splicing assays. A total of 10 VUSs were predicted by two or more programs to induce a significant reduction of splice site strength or activation of cryptic splice sites or generation of new splice sites. Minigene-based splicing assays confirmed four of these predictions. Five additional VUSs, all at internal exon positions, were not predicted to induce alterations of splice sites, but revealed variable levels of exon skipping, most likely induced by the modification of exonic splicing regulatory elements. We provide new data in favor of the pathogenic nature of the variants BRCA1 c.212+3A>G and BRCA1 c.5194-12G>A, which induced aberrant out-of-frame mRNA forms. Moreover, the novel variant BRCA2 c.7977-7C>G induced in frame inclusion of 6 nt from the 3' end of intron 17. The novel variants BRCA2 c.520C>T and BRCA2 c.7992T>A induced incomplete skipping of exons 7 and 18, respectively. This work highlights the contribution of splicing minigene assays to the assessment of pathogenicity, not only when patient RNA is not available, but also as a tool to improve the accuracy of bioinformatics predictions.