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.

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,

INTRODUCTION
Since the discovery of the breast cancer genes BRCA1 (OMIM #113705) and BRCA2 (OMIM #600185) (Miki et al., 1994;Wooster et al., 1995), nearly 17,000 different variants of both genes have been recorded at the ClinVar database 1 (date last accessed; November 2018). Germline inactivating variants in BRCA1 and BRCA2 confer high lifetime risks of breast and ovarian cancers (Mavaddat et al., 2013). Also, other cancer types, such as prostate, pancreatic and melanoma, are associated with pathogenic variants in these genes (Petrucelli et al., 2013). Despite the high penetrance of BRCA pathogenic variants, they are responsible for only ∼15-20% of hereditary breast and ovarian cancer (HBOC) (Stratton and Rahman, 2008). In fact, HBOC is a highly genetically heterogeneous disease with about 25 known or proposed susceptibility genes (Nielsen et al., 2016). Apart from the BRCA genes, PALB2 (OMIM #610355), ATM (OMIM #607585), and CHEK2 (OMIM #604373) have a prominent contribution since, in a recent study, more than 30% of pathogenic variants were found in these genes (Buys et al., 2017).
Commonly, the variants are classified attending to their predicted effect on the protein so that truncating variants (frameshift and nonsense) are directly classified as pathogenic, while intronic, missense and synonymous variants are usually considered to be variants of uncertain clinical significance (VUS). In fact, VUS are identified by a relevant proportion of BRCA genetic tests (∼20%), which hamper genetic counseling and subsequent preventive or therapeutic actions, since risk assessment is solely based on family history (Radice et al., 2011;Eccles et al., 2015;Ricks et al., 2015).
Furthermore, other upstream gene-expression processes, such as transcription or splicing, can be impaired if regulatory motifs are targeted by nucleotide variations (Wang and Cooper, 2007). Splicing is the process by which introns are removed from a pre-mRNA and exons are consecutively joined. This mechanism is performed in the nucleus by the spliceosome, a macrocomplex constituted by 5 small nuclear ribonucleoproteins (snRNPs) and many other associated proteins (De Conti et al., 2012). The spliceosome recognizes in the pre-mRNA specific sequences which define the exons/introns boundaries and other elements needed to carry out the process. These sequences are: the acceptor or 3 splice site (3 ss), the donor or 5 splice site (5 ss), the branch point, the polypyrimidine tract and the auxiliary cis sequences known as splicing regulatory elements (SREs) where enhancer or silencer trans factors can bind. Therefore, any change in the sequence may disrupt splicing (Cartegni et al., 2002). Splicing variants usually break the 3 ss or 5 ss leading to abnormal splicing events such as exon skipping, alternative site usage or intron retention. However, they may also create new splicing sites or strengthen cryptic ones that would then be recognized. Other mechanism that may alter splicing is the disruption of exonic/intronic splicing enhancers (ESEs/ISEs) or the creation of exonic/intronic splicing silencers (ESSs/ISSs) (Abramowicz and Gos, 2018). Nevertheless, it is extremely difficult to identify active 1 https://www.ncbi.nlm.nih.gov/clinvar SREs and predict the impact of the DNA variants on splicing given the low accuracy of SRE-detection softwares. Therefore, splicing variants can induce abnormal transcripts that either introduce premature termination codons (PTC), in-frame loss of essential protein domains or even inclusion of new translated sequences. Consequently, variants with impact on splicing (or spliceogenic variants) may be associated with an increased risk of a given disease. This ethiopathogenic mechanism has been so far underestimated, even though some authors have suggested that spliceogenic variants may represent more than 60% of diseasecausing mutations (López-Bigas et al., 2005).
Previous studies have shown that a significant number of splicing variants have been detected in BRCA2 (Spurdle et al., 2008;Rebbeck et al., 2018). In fact, previous results from our group showed that more than 50% of tested variants of BRCA2 exons 16-27 impaired splicing (Acedo et al., 2012(Acedo et al., , 2015Fraile-Bethencourt et al., 2017. Likewise, at least 24 different BRCA2 alternative transcripts have been identified. They are helpful to interpret the splicing outcomes of genetic variations (Fackenthal et al., 2016) and suggest a fine regulation of BRCA2 exon processing. This feature is supported by the fact that several ESE-rich regions have been functionally mapped by exonic deletions throughout most BRCA2 exons. Thus, these motifs would be involved in precise exon recognition and alternative splicing events (Acedo et al., 2015;Fraile-Bethencourt et al., 2017). Moreover, we showed that functional mapping is an optimal approach that improves ESE-software predictions and facilitates the identification of spliceogenic mutations of this sort of cis-elements.
In this work, we have extended our analysis to BRCA2 exons 14 and 15 by carrying out an in-depth study of candidate spliceogenic variants. We have explored the presence of splicing enhancers in exons 14 and 15 and have undertaken RNA assays of 52 selected variants from both exons.

MATERIALS AND METHODS
Ethical approval for this study was obtained from the Ethics Review Committee of the Hospital Universitario Río Hortega de Valladolid (6/11/2014).

Bioinformatics: Databases and in silico Analysis
We collected BRCA2 variants from the main databases: ClinVar 2 , the BRCA Share Database (UMD 3 ) (Beroud et al., 2016) and the Breast Cancer Information Core (BIC 4 ) ( Supplementary  Table S1). Variants and transcripts were annotated according to the Human Genome Variation Society (HGVS) guidelines on basis of the BRCA2 GenBank sequence NM000059.1. In order to simplify, we identified transcripts with a shortened code that combines the following symbols (Lopez-Perolio et al., 2019): (skipping of reference exonic sequences), (inclusion of reference intronic sequences), E (exon), p (acceptor shift), q (donor shift). When necessary, the exact number of skipped or retained nucleotides is indicated. For example, transcript (E14p10) indicates the use of an alternative acceptor site 10-nt downstream that causes a 10-nt deletion.

Minigene Construction and Mutations
The minigene MGBR2_14-20 was built as previously reported (Fraile-Bethencourt et al., 2017). A total of 52 variants and 8 microdeletions were introduced into the wild type (wt) minigene by site-directed mutagenesis with the QuikChange Lightning Kit (Agilent, Santa Clara, CA, United States), following the manufacturer's instructions (Supplementary Table S2). All mutant clones were confirmed by sequencing (Macrogen, Madrid, Spain).

RT-PCR and Transcripts Amplification
Retrotranscription was carried out with 400 ng of RNA and the RevertAid First Strand cDNA Synthesis Kit (Life Technologies), using the specific minigene primer RTPSPL3-RV (5 -TGAGGAGTGAATTGGTCGAA-3 ). Samples were incubated at 42 • C for 1 h, followed by 5 min at 70 • C. Transcripts were amplified with Platinum Taq DNA polymerase (Life Technologies) using 40 ng of cDNA and the primers pMAD_607FW (Patent P201231427, CSIC) and RTBR2_ex17RV2 (5 -GGCTTAGGCATCTATTAGCA-3 ). PCR consisted of: denaturation step at 94 • C for 2 min, followed by 35 cycles 94 • C-30 s, 60 • C-30 s and 72 • C-1 min/kb, and a final extension step at 72 • C for 5 min. Transcripts were sequenced at the Macrogen Spain facility.
In order to relatively quantify all transcripts, semi-quantitative fluorescent RT-PCRs were undertaken in triplicate with the primers pMAD_607FW (FAM-labeled) and RTBR2_ex17RV2 and Platinum Taq DNA polymerase (Life Technologies) under standard conditions except that 26 cycles were herein applied (Acedo et al., 2015). FAM-labeled products were run with LIZ-1200 Size Standard at the Macrogen facility and analyzed with the Peak Scanner software V1.0 (Life Technologies). Only peak heights ≥ 50 RFU (Relative Fluorescence Units) were considered.

Minigene Construction and ESE Mapping
The minigene MGBR2_14-20 had previously been used to study variants of exons 16, 17, and 18, proving that it is a reliable and robust tool to functionally assay splicing variants (Fraile-Bethencourt et al., 2017Montalban et al., 2018). The MGBR2_14-20 is a 10.7 Kb construct which, after transfection in MCF-7 cells, produces a transcript with the following structure: V1-BRCA2_exons from 14 to 20-V2 (1,806 nt) (Figure 1). To study exons 14 and 15, cDNA was amplified with a forward primer located in V1 (pMAD_607FW) and a reverse primer located in exon 17 (RTBR2_Ex17RV2), with an expected transcript size of 1028 nt (Figures 1B,C).
To map ESEs, 30-nt overlapping microdeletions were performed along the first and the last 55-nt of exons 14 and 15 (Acedo et al., 2015), always preserving the splice site conserved positions (the first 2 nt and the last 3 nt of the exon). None of the deletions but one altered splicing, suggesting the absence of cis-regulatory motifs within these segments. Only exon 15 deletion c.7463_7492

Variant Collection and Bioinformatics Analysis
A total of 294 different variants, spread throughout exons 14 and 15 and flanking introns, were collected from the main databases (ClinVar, UMD and BIC) (Supplementary Table S1). They were analyzed in silico with MES and NNSplice for splice site prediction, and with ESE/ESS estimation algorithms integrated in Human Splicing Finder (HSF). Potential splicing variants were selected following these criteria: creation or disruption of splice sites (according to MES or NNSplice); disruption of the branch point; disruption of the polypyrimidine tract; elimination of enhancers or creation of silencers. Some of the selected variants had a combined effect, for example, they were predicted to simultaneously create an ESS and removed an ESE. A total of 53 candidate variants (∼18%) that included 19 intronic, 18 missense, 5 nonsense, 8 synonymous, and 3 frameshift variants were selected ( Table 1). According to their previous clinical classification, the selection contained: 8 benign or likely benign variants, 30 VUS and 15 pathogenic or likely pathogenic variants.
Bioinformatics indicated that 13 variants disrupted the natural splice sites, three decreased their scores (one disrupted the polypyrimidine tract), 11 created new splice sites, one decreased the branch point score (HSF: 79.34→55.5), three altered intronic splicing elements (ISEs and ISSs) and 22 altered exonic splicing elements (ESEs and ESSs) ( Table 1). Exceptionally, variants c.7008-5T > C (ivs13-5T > C) and c.7435+5T > C (ivs14+5T > C) were also selected because, even though the bioinformatics did not show significant score changes (MES: 10.37→8.48 and 5.64→5.56, respectively), these positions are relevant for exon recognition. Thus, c.7008-5T > C is the closest position of the polypyrimidine tract to the canonical acceptor site, and c.7435+5T > C is part of the consensus 5 ss sequence and +5 changes were previously associated with disease (Sharma et al., 2014;Montalban et al., 2018).

Functional Splicing Assays Using the Minigene MGBR2_14-20
A total of 52 variants were genetically engineered in the wt MGBR2_14-20 by using specific primers (Supplementary Table  S2). Despite 53 variants were initially selected, mutagenesis experiments did not work for c.7598C > G variant. The 52 mutant minigenes were checked by Sanger sequencing and assayed in MCF-7 cells. Results showed that 12 of them (23%) altered splicing (Table 2 and Figure 3), seven of which had previously been classified as pathogenic and five as VUS ( Table 3). Among these 12 variants, there were 9 intronic, 2 missense and 1 nonsense changes. Functionally, the 9 intronic variants (c.7008 7617+1G > T, and c.7617+2T > G) disrupted the natural splice sites, the two missense changes (c.7177A > G and c.7447A > G) triggered the use of de novo splice sites and originated other transcripts, and the nonsense one (c.7501C > T) probably altered SREs despite it was primarily selected because of the creation of a new 5 ss ( Table 1). According to Vallée et al. (2016), nine variants induced severe splicing disruptions as they produced more than 60% of aberrant transcripts, ranging from 92.8 to 100% (Table 3).

Splicing Regulatory Element-Variants
A total of 26 SRE variants were selected according to the criteria above described and assayed in MCF-7 cells ( Table 1). Only

Analysis of Transcripts
The so called full-length or canonical transcript (expected size: 1028-nt) was amplified with primers placed on vector exon V1 and BRCA2 exon 17. Apart from the canonical transcript, we have detected other ten different ones (Figure 3). Aberrant exon 14 splicing produced six different isoforms, but only five of them were characterized: (E14), (E14p10), (E14p246),

DISCUSSION
Due to the implementation of Next Generation Sequencing in the clinical setting (Slavin et al., 2015), a large number of variants have been detected in disease-responsible genes. HBOC and the breast cancer susceptibility genes are not exceptions, where thousands of different variants have been reported although many of them are considered as VUS (Spurdle et al., 2012;Slavin et al., 2017). In this context, the functional and clinical classifications pose a challenge for Medical Genetics. We have herein functionally assayed 52 BRCA2 variants using the minigene MGBR2_14-20, whose reliability had been previously proven (Fraile-Bethencourt et al., 2017Montalban et al., 2018). We found 12 variants that altered splicing, nine of which would severely alter the protein. This study forms part of a comprehensive study of our group concerning potential splicing BRCA2 (Vreeswijk et al., 2009;Colombo et al., 2013;de Garibay et al., 2014). Moreover, another 31 variants of this and other constructs replicated previously reported patient splicing outcomes (Acedo et al., 2015;Fraile-Bethencourt et al., 2017, confirming that splicing reporters are robust and valuable tools to test the impact of variants on splicing, especially when patient RNA is not available. Here, we have shown that nine variants drastically disrupted the splicing pattern. We found five 3 ss disrupting variants, one of which (c.7008-3C > G) provoked exon 14 skipping, whereas the rest of them induced the use of cryptic sites (c.7008-2A > T, c.7008-1G > A, c.7436-2A > T, and c.7436-1G > A). Moreover, variant c.7435+1G > C, which disrupt the exon 14 donor site, provoked the use of a cryptic donor. Curiously, neither c.7435+3A > G, c.7435+5T > C nor c.7435+6G > A, which are also part of the consensus 5 ss, affected splicing. This may be due to the low frequency of +5T and +6G at these positions, so that any nucleotide change only equals or improves splice site strength. Conversely, +3A is the main nucleotide at this position (71%) but +3G is also relatively frequent (24%) (Zhang, 1998). Thus, a substitution A to G might have a reduced or no splicing impact, as it is the case of the c.7435+3A > G variant. As a matter of fact, variant c.9501+3A > T produced 87% of the canonical transcript (Acedo et al., 2015). On the other hand, the three variants of exon 15 donor site c.7617+1G > A, c.7617+1G > T, and c.7617+2T > G resulted in exon skipping. In this regard, it was recently recommended the use of a combination of the computational tools HSF plus Splice Site Finder-like to select candidate splice site variants with high sensitivity and specificity (Moles-Fernández et al., 2018). According to HSF, variants c.7435+3A > G, c.7435+5T > C, and c.7435+6G > A showed only minimal changes of the splice site scores (±1%) so that they should have been excluded from subsequent functional tests.
Concerning other splicing motifs, minigenes also allow the identification of regulatory sequences (splicing enhancers and silencers) and splicing factors involved in the specific regulation (Baralle and Buratti, 2017). Indeed, the SRE mapping constitutes an interesting experimental approach since it identifies critical regions for exon recognition. In this context, our group previously found exonic variants that disrupted splicing through elimination of ESEs or creation of de novo silencers, such as BRCA2 variants c.7985C > G (predicted missense p.Thr2662Arg) or c.8009C > A (predicted nonsense p.Ser2670 * ) (Fraile-Bethencourt et al., 2017). The splicing assays showed how both variants elicited complete splicing aberrations (mainly exon 18 skipping). However, they were a priori classified as missense and nonsense variants, respectively, due to their predicted protein effect. After microdeletion mapping, we have identified a putative ESE-region in exon 15 (c.7463_7492). These ESEs might be involved in exon 15 3 ss recognition since their loss produced the use of a cryptic acceptor 83-nt downstream (Figure 2). Curiously, we did not find any ESE-variants that affected pre-mRNA processing. Only, c.7501C > T, which lays near to this presumed ESE interval, provoked a similar outcome to that of r.7463_7492 deletion ( Table 2). In summary, we have tested a total of 117 different variants in minigene MGBR2_14-20, from exons 14-18 ( Table 4), 51 of which (43.6%) induced abnormal splicing patterns: 31 disrupted the natural splice sites (16 3 ss and 15 5 ss), 11 affected SREs, six created de novo splice sites and three altered the polypyrimidine tract.

Transcript Interpretation
The BRCA2 exons 14 (c.7008_7435) and 15 (c.7436_7617) encode for amino acids 2336 to p.2539, which are part of the DNA Binding Domain (DBD; p.2459_p.3190). The DBD is the largest conserved region of BRCA2 and is composed of a helical domain (HD), three oligonucleotide binding sites (OB) and a tower domain (TD) (Guidugli et al., 2014). Specifically, exons 14 and 15 are part of the HD (p.2481_2667) that binds to the protein DSS1 (deleted in split-hand/split-foot syndrome) in the region comprised by the residues 2472-2957 (Marston et al., 1999). Among these residues, a total of 125 are strictly conserved from human to sea urchin. DSS1 is important for BRCA2 stability, since its loss leads to a   (Li et al., 2006). Moreover, exons 14 and 15 coding region is also recognized by FANCD2 (Fanconi anemia group D2) protein, which binds to the BRCA2 protein between codons 2350 and 2545 (Hussain et al., 2004). FANCD2, like BRCA2, is one of the 16 proteins that form the Fanconi Anemia complex, aimed to repair DNA interstrand crosslinks. However, it was shown that BRCA2-FANCD2 association has an independent function in the Fanconi Anemia pathway. The BRCA2-FANCD2 complex is involved in the restart of the replication fork, by protecting the nascent DNA strands from degradation (Raghunandan et al., 2015). Taken together, exons 14 and 15 contain crucial sequences of BRCA2, owing to its function in homologous recombination and other relevant biological processes. Moreover, the biological relevance of exons 14 and 15 is supported by the presence of numerous pathological nonsense and frameshift variants at the mutation databases. Hence, the exon 14-15 spliceogenic variants that induce PTC-transcripts may be associated with an increased risk of HBOC.

Clinical Interpretation of Spliceogenic Variants
Twelve variants altered splicing with different patterns. While some variants caused the total absence of canonical transcript, others originated just ∼5% of aberrant transcripts (  Table 3). Variant c.7008-3C > G produced (E14) as the major transcript, but the full-length (∼7%) and the in-frame transcript (E14p246) (∼13%) were also identified. The (E14p246) isoform contains a deletion of 82 non-conserved amino acids (p.Thr2337_Arg2418del) that form part of the FANCD2 binding site (p.Thr2350_Val2545). At the UMD database, c.7008-3C > G is classified as a "causal" variant because of the skipping of exon 14, but ClinVar shows it as VUS. According to the ACMG criteria, this variant should be cataloged as likely pathogenic (Table 3). On the other hand, according to the guidelines of the ENIGMA consortium 7 eight variants should be classified as Class 4 (Likely pathogenic), three as VUS and one as pathogenic, though this is due to its predicted nonsense effect (Table 3).
On the other hand, missense variants c.7177A > G (p.Met2393Val) and c.7447A > G (p.Ser2483Gly) produced the canonical transcript as the main outcome and only 5 and 10% of aberrant isoforms, respectively ( Table 2). The canonical transcript generated by these two variants carried missense changes, but in silico predictions with the PolyPhen tool 8 suggested that p.Met2393Val and p.Ser2483Gly were both benign changes. So, following the ACMG criteria, c.7177A > G and c.7447A > G remain as VUS (Table 3). Finally, variant c.7501C > T generated mainly the canonical transcript (84%) that includes a nonsense change (pathogenic according to ClinVar) so this change should be classified as pathogenic under a combined splicing-protein viewpoint. Altogether, we have re-classified three variants (c.7008-3C > G, c.7008-1G > A, and c.7436-1G > A) from VUS to pathogenic or likely pathogenic, and we have provided further support for the classification of six spliceogenic variants (c.7008-2A > T, c.7435+1G > C, c.7436-2A > T, c.7617+1G > A, c.7617+1G > T, and c.7617+2T > G). Interestingly, c.7617G > A is classified as "causal" in the UMD database 9 and indicated that causes exon 15 skipping but no functional proofs were provided. However, the functional assay of MGBR2_14-20-c.7617G > A only showed the canonical transcript (Supplementary Figure  S2). In fact, NNSplice, HSF and MES estimated just a small decrease (−19, −11.6, and −24%, respectively) of the donor site score. Therefore, c.7617G > A behaves as a neutral variant from the splicing perspective. Finally, the minigene MGBR2_14-20 results of 5 exons suggest that a total of 39 spliceogenic variants should be classified as pathogenic or likely pathogenic (Table 4), lending further support to this strategy for the clinical interpretation of variants.
In summary, splicing is a finely regulated mechanism which can be altered by any change in the sequence. The disruption of this process might cause serious effects on the protein, from the loss of important domains to the generation of a PTC. Thus, splicing alteration is a common mechanism of gene inactivation, which is often involved in human disease. Here, we have revealed 12 spliceogenic variants of BRCA2 exons 14 and 15. The minigene based assays offer a relevant information about effects of splicing variants, since they allow to functionally assay almost any DNA change, to quantify all generated transcripts, including very rare 7 https://enigmaconsortium.org/library/general-documents/enigmaclassification-criteria/ 8 http://genetics.bwh.harvard.edu/pph2/ 9 https://goo.gl/bG2V4n ones, and to initially study the splicing regulation. Moreover, we have detected an ESE region that seems to be regulating exon 15 splicing, and therefore constitutes a hypothetical hotspot for putative ESE-mutations. Indeed, pSAD-based minigenes have constituted an invaluable technology to functionally test variants of other disease genes such as GRN (Frontotemporal Dementia), SERPINA1 (Severe alpha-1 antitrypsin deficiency) and CHD7 (Charge Syndrome) genes (Lara et al., 2014;Villate et al., 2018). Altogether, these results highlight, once more, the importance of RNA assays to know the splicing effects of DNA variants to give support to their clinical interpretation and consequently to activate preventive and/or therapeutic interventions.

AUTHOR CONTRIBUTIONS
EF-B contributed to the bioinformatics analysis, minigene construction, manuscript writing, and performed most of the splicing functional assays. BD-G and AV-P participated in minigene construction, mutagenesis experiments, and functional assays. SG-B and MC carried out the data collection of variants and their computer analysis, as well as manuscript editing. EV conceived the study and the experimental design, supervised all the experiments, and wrote the manuscript. All authors contributed to data interpretation, revisions of the manuscript, and approved the final version of the manuscript. FUNDING EV's lab was supported by grants from the Spanish Ministry of Science, Innovation and Universities, Plan Nacional de I+D+I 2013-2016, Instituto de Salud Carlos III (PI13/01749 and PI17/00227) co-funded by FEDER from Regional Development European Funds (European Union), and grants CSI090U14 (ORDEN EDU/122/2014) and CSI242P18 (actuación cofinanciada P.O. FEDER 2014-2020 de Castilla y León) from the Consejería de Educación, Junta de Castilla y León. EF-B was supported by a predoctoral fellowship from the University of Valladolid and Banco de Santander (2015-2019). AV-P was supported by a predoctoral fellowship from Consejería de Educación, Junta de Castilla y León (2018)(2019)(2020)(2021)(2022).