Comprehensive splicing functional analysis of DNA variants of the BRCA2 gene by hybrid minigenes
Breast Cancer Research volume 14, Article number: R87 (2012)
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.
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.
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.
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.
Germline mutations in the BRCA1 (MIM# 113705) and BRCA2 (MIM# 600185) genes confer a high lifetime risk of developing breast/ovarian cancer  and account for about 16% of the breast cancer familial risk . Nearly 3,500 different DNA variants of BRCA1 and BRCA2 have been reported at the Breast Cancer Information Core Database (BIC) . Only truncating mutations (1,457 nonsense and frameshift variants, 41.6%), variants in natural splice sites (141, 4%) and a minor fraction of the 1,346 BIC missense variants have been catalogued as deleterious [4, 5]. Little is known about the pathogenicity of most of the remaining variants (approximately 50%), usually referred to as unclassified variants (UV), which complicates genetic counseling in UV carrier families . A more reliable interpretation and classification of these sequence changes will improve the clinical management of cancer patients and their families. Several methods have been developed to classify them, including functional assays  and integrated statistical analysis of several parameters [8, 9].
Genetic variants in disease-responsible genes that disrupt the splicing code have a key role in human hereditary disorders [10, 11]. In fact, it has been estimated that more than 60% of sequence variations may affect pre-mRNA processing . Therefore, splicing disruptions should be considered a core mechanism of gene inactivation to be investigated in UVs. Splicing mutations are traditionally considered those that affect essential nucleotides of the donor (GT) and the acceptor (AG) sites , although other intronic and exonic nucleotides are critical for splice site selection . Furthermore, exon recognition is supported by other essential cis-regulatory motifs, the so-called splicing enhancers (ESE, ISE) or silencers (ISE, ISS) that are usually bound by SR and hnRNP proteins, respectively. They are short (6 to 8 nucleotides) and degenerate sequences that can promote (enhancers) or repress (silencers) exon inclusion in the mature mRNA. Sequence variations in any of these elements can result in alterations of the pre-mRNA processing step that can also be associated with human disorders, including hereditary cancer [14–18].
Identification of the splicing regulatory elements (SRE) by computational predictions is not accurate because of the degeneracy of these motifs [11, 14, 19, 20], supporting the importance of RNA analysis to demonstrate the deleteriousness of a particular DNA change. However, it is often complex to obtain extra blood samples for RNA extraction from patients carrying a suspicious DNA change. Splicing functional analysis with hybrid minigenes is a simple and robust approach to study potential variants with impact on splicing without the need of collecting additional blood samples from patients. The genomic region of interest (exons and their flanking introns) from control and affected individuals can be cloned into a splicing reporter plasmid such as pSPL3. The construct is then transiently transfected into eukaryotic cells and the splicing pattern is analyzed by RT-PCR with specific primers of plasmid exons. The minigene assay constitutes a useful approach for identifying splicing anomalies and the study of their underlying functional mechanisms [21, 22]. We previously studied 57 different pre-selected DNA variants from BRCA1 and BRCA2 using a combined approach: bioinformatics analysis and splicing assays by lymphocyte RNA RT-PCR and/or hybrid minigenes. Twenty-eight out of the 57 pre-selected variants displayed abnormal splicing patterns, suggesting that disruption of this process is an important disease mechanism .
In this study, we have performed a splicing analysis of BRCA2 variants c.8488-1G > A (acceptor site of exon 20) and c.9026_9030del (exon23) by lymphocyte RT-PCR and splicing reporter minigenes of exons 19-20 and 23-24. Moreover, we have extended this analysis to another 41 candidate splicing DNA variants of these exons reported in the BIC mutation database. A total of 19 variants induced aberrant splicing and we identified and quantified up to 12 different aberrant RNA isoforms from minigenes of exons 19-20 and 23-24.
Materials and methods
Patients, nucleic acid isolation and mutation detection
Breast and/or ovarian cancer patients harboring variants c.8488-1G > A and c.9026_9030del (nine families) were selected by the Genetic Counseling Unit (Complejo Hospitalario de Burgos) . Written, informed consent was obtained from all patients prior to blood extraction. This study was approved by the ethics committee of the Faculty of Medicine, University of Valladolid (Spain). DNA and RNA were purified from peripheral blood lymphocytes of patients VA1220 (c.8488-1G > A) and VA1612 (c.9026_9030del) by using the QIAamp DNA and RNA blood mini kit (Qiagen, Hilden, Germany), respectively. RNA purification included a DNAse I treatment step. Mutation detection was performed by capillary heteroduplex analysis [24, 25]. Nomenclature of sequence variations follows the guidelines of the Human Genome Variation Society (HGVS) , and was based on the BRCA2 cDNA sequence NM_000059.
Splicing prediction programs
To identify potential variants with impact on splicing, mutant and reference sequences were analyzed with the following programs: a) NNSPLICE, which evaluates the strength of splice sites , and Human Splicing Finder (HSF) , which includes several matrices to analyze splice sites and splicing silencers and enhancers (for example, MaxEnt, ESEfinder, PESX and so on). In order to improve the specificity of ESEfinder predictions  we examined the evolutionary conservation of ESE motifs by aligning the BRCA2 sequences of different organisms with CLUSTALW2 .
Lymphocyte RNA from patients carrying variants c.8488-1G > A (exon 20) and c.9026_9030del (exon 23) was retrotranscribed with the High Capacity kit (Applied Biosystems, Carlsbad, CA, USA). Amplification of cDNA was conducted with GoTaq Hot Start DNA polymerase (Promega, Madison, WI, USA) and flanking exonic primers forward RTBR2_EX19-FW, 5'AAACTTGGATTCTTTCCTGAC 3' (exon 19), and reverse RTBR2_EX23-RV, 5'ATTCTGTATCTCTTTCCTTCTGTT 3' (exon 23), for variant c.8488-1G > A (amplicon length: 653 bp), and RTBR2_EX22-FW, 5' GTGAAGAGCAGTTAAGAGCCT 3' (exon 22), and RTBR2_EX26-RV, 5' CTTGTTTTCTGCTTCATTGC 3' (exon 26) for variant c.9026_9030del (amplicon length: 773 bp). Splicing profiles of patients were compared with those obtained from RT-PCR experiments of control lymphocytes and non-tumorigenic breast epithelial cells (MCF10A; ATCC, LGC Standards) . RT-PCR products were run on 1.5% agarose gels and stained with ethidium bromide. Isoforms induced by variant c.8488-1G > A were sequenced with the Big Dye 3.1 Sequencing Kit (Applied Biosystems) and primers RTBR2_EX19-FW and RTBR2_EX23-RV.
Construction of minigenes
Wild type (wt) exons 19-20 (intron 18-exon 19-intron 19-exon 20-intron 20) and 23-24 (intron 22-exon23-intron 23-exon 24-intron 24) of BRCA2 were amplified with Pfx50 high fidelity polymerase (Invitrogen, Carlsbad, CA, USA) and primers with a restriction site (underlined) for either XhoI or BamHI: MGBR2_ex19-20_XhoI-FW, 5' CACACACTCGAGATAGCATTAAGAACTTGTAGCA 3', and MGBR2_ex19-20-BamHI-RV, 5' CACACAGGATCCATTACAAATGGCTTAGACCTGA 3' (Exons 19-20; size: 1164 bp); and MGBR2_ex23-24_XhoI-FW, 5' CACACACTCGAGAATGCTTTGTTTTTATCAGTTTT 3', and MGBR2_ex23-24_BamHI-RV: 5' CACACAGGATCCAATTTGCCAACTGGTAGCTC 3' (exons 23-24; size: 717 bp).
After digestion with XhoI and BamHI (Fermentas, Vilnius, Lithuania), fragments were ligated to the exon trapping vector pSPL3 (Invitrogen, discontinued) to transform DH5α competent cells (Invitrogen). The new constructs were sequenced with the Big Dye 3.1 Sequencing Kit (Applied Biosystems) and primers PSPL3-SEQ-FW (5' CCTTGGGATGTTGATGAT 3') and PSPL3-SEQ-RV (5' TTGCTTCCTTCCACACAG 3').
Site directed mutagenesis
Mutagenesis was carried out according to the PCR mutagenesis protocol with Pfu Turbo DNA polymerase (Agilent, Santa Clara, CA, USA) . Wild type (wt) minigenes of BRCA2 exons 19-20 and 23-24 were used as templates to generate 41 candidate variants from the BIC database [see Additional file 1, Table S1]. In addition, several artificial variants were designed to target putative regulatory motifs: c.8484A > T and c.8486A > C, c.8512T > A, c.8518del and two combined BIC variants c.[8609A > G;8611G > T] and c.[8972G > A; 9006A > T] within the same construction. Finally, two conserved motifs of intron 19 were deleted separately and together: c.8487+31_8487+42del, c.8488-49_8488-44del, and c.[8487+31_8487+42del; c.8488-49_8488-44del].
Transfection of HeLa cells
Approximately 105 HeLa cells (human cervical carcinoma) were grown to 90% to 95% confluency in 0.5 mL of growth medium ((D)MEM, 10% fetal bovine serum, 1% glucose and 1% penicillin/streptomycin) in four-well plates (Nunc, Roskilde, Denmark). Cells were transfected with 1 μg of plasmid minigene and 2 μL of lipofectamine 2000 (Invitrogen) following the manufacturer's protocol. After 48 hours RNA was purified with Nucleospin-RNA-II (Macherey-Nagel, Düren, Germany) that includes on-column rDNAse treatment and quantified in a Nanodrop 1000 spectrophotometer (Fisher Scientific, Spain).
RT-PCR of minigenes and quantification of mRNA isoforms
Retrotranscription was carried out with 200 ng of RNA and the transcriptor first strand cDNA synthesis kit (Roche, Sant Cugat del Vallés, Barcelona, Spain). Semiquantitative fluorescent RT-PCRs were performed in triplicate in a final volume of 20 μL that contained 2 μL of cDNA and the flanking primers of constitutive exons of pSPL3, SD6-PSPL3_RTFW (FAM-5' TCACCTGGACAACCTCAAAG 3') and SA2-PSPL3_RTREV (5' TGAGGAGTGAATTGGTCGAA 3'). Samples were denatured at 94°C for five minutes, followed by 26 cycles consisting of 94°C for 20 seconds, 58°C for 20 seconds, and 72°C for 30 seconds, and a final extension step at 72°C for two minutes. Sizes of the RT-PCR products were 486 and 488 bp for minigenes 19-20 and 23-24, respectively. The RT-PCR products (1 μL of a 1/10 dilution) were mixed with 18 μL of Hi-Di Formamide (Applied Biosystems) and 0.2 μL of Genescan 500 Rox Size Standard (Applied Biosystems). Samples were run on an ABI3130 sequencer using POP7 polymer and analyzed with the Peak Scanner software (Applied Biosystems). Mean peak areas were used to calculate ratios of the different splicing isoforms generated by DNA variants from minigenes 19-20 and 23-24.
Sequencing of the minigene RT-PCR products was carried out as described above with primers SD6-PSPL3_RTFW and SA2-PSPL3_RTREV except for intron 19 retention isoform that was sequenced with a specific primer in the boundary intron 19-exon 20 (RTBR2_ivs19-ex20-RV: 5' TGTCTTCTCCATCCACTGTAAT 3').
Results and discussion
Variant c.8488-1G > A was detected in a patient who developed bilateral breast cancer at age 66 and 74 and ovarian cancer at age 77. This variant affects the essential nucleotide -1 (G to A) of the acceptor splice site of BRCA2 exon 20. This variant had not been reported at the BIC mutation database although it had formerly been found in a consanguineous Fanconi anemia patient . Analysis with HSF at the intron/exon boundary identified a weak acceptor splice site that was even weakened by the change -1 G > A. Lymphocyte RT-PCR of the carrier patient revealed an upper band with intron 19 retention (1053 bp) and a transcript with a 12-nucleotide deletion of exon 20 generated by the use of an alternative cryptic acceptor [see Additional file 1, Figure S1]. This variant was also evaluated by hybrid minigenes of exons 19-20. Splicing assay and fluorescent RT-PCR (Figure 1B and Additional file 1, Figure S1) revealed three main isoforms: 12-nucleotide deletion of exon 20 (9.7%), intron 19 retention (16.0%), exon 20 skipping (72.1%), which was not previously identified in lymphocytes, and other minor isoforms. The different splicing profiles between both assays may be due to the different genomic context where adjacent exons may be involved in the splicing efficiency of exons 19 and 20 . The most common events of alternative splicing are exon skipping and alternative splice site selection whereas intron retention is the less frequent phenomenon in physiological alternative splicing , indicating that its contribution to the complexity of the human proteome is low since this event is usually associated with the introduction of premature stop codons . Moreover, intron retention is increased by two-fold in cancerous cells, suggesting the disruption of essential repressor genes of tumor progression . Nevertheless, intron retention rate can be particularly high in some specific genes such as the Kallikrein gene family of serine proteases, in which six isoforms showed intron III exonization and were predicted to truncate all the resultant proteins .
Translation of the intron 19 retention and exon 20 skipping isoforms was predicted to result in truncated proteins (p.Gln2829fsX2 and p.Trp2830LysfsX13, respectively), whereas the 12-nucleotide deletion produced an in frame deletion of four amino acids (p.Trp2830_Lys2833del, W-M-E-K) of the BRCA2 DNA binding domain. Tryptophan 2830 and glutamic acid 2832 are strictly conserved in vertebrates, (IARC BRCA2 alignment), suggesting their functional relevance . Moreover, one Fanconi anemia patient of subtype D1 (associated with BRCA2 mutations) with a family history of consanguinity was previously reported to be homozygous for this variant . Altogether these data strongly support the pathogenicity of c.8488-1G > A.
BRCA2 c.9026_9030del of exon 23 is a deleterious frameshift mutation that would lead to a premature stop codon (p.Tyr3009SerfsX9). It is one of the most prevalent mutations in the Spanish population, accounting for 7.6% of BRCA2 families in Castilla y León (Spain) . We had already shown that other deletions, such as BRCA2 c.470_474del, can induce splicing defects . HSF analysis predicted changes in several presumed regulatory elements (Exon and Intron Identity Elements) . We performed lymphocyte RT-PCR of one carrier patient and we found the wt allele and the corresponding 5-nucleotide deletion transcript derived from the genomic variant (c.9026_9030del) without any splicing anomaly (data not shown). This result was confirmed in the minigene of exons 23 and 24 (Figure 1F).
Minigene analysis of BIC variants
We had previously shown that first, more than 20% of BRCA variants reported at the BIC database can induce splicing disruptions and, second, the combination of PCR mutagenesis and splicing reporter plasmids is a straightforward and strong approach that allows the analysis of the splicing effect of any sequence change . We, therefore, proceeded to analyze all the reported variants of exons 19, 20, 23 and 24 from the BIC database (155 different DNA changes) with NNSplice and HSF. We chose variants that met one of the following criteria: disruption of the natural splice site, creation of novel alternative donor or acceptor sites, alteration of evolutionarily conserved ESEs, and creation of splicing silencers. A total of 41 variants (26.5%) were preselected and generated by PCR mutagenesis in the wt minigenes 19-20 and 23-24. Splicing functional assays were carried out by semiquantitative fluorescent RT-PCR and the peak areas of the different isoforms were quantified. Eighteen variants (43.9%) produced aberrant splicing patterns (Table 1, Figure 1 and Additional file 1, Figure S2) that affected the natural splice sites (12 variants: c.8487+1G > A, c.8487+3A > G, c.8488-2A > G, c.8954-3C > G, c.8954-1_8955delinsAA, c.9117G > A, c.9117+1G > T, c.9117+1G > A, c.9118-2A > G, c.9248_9256+7del, c.9256G > T and c.9256+1G > A), splicing enhancers or silencers (four variants: c.8378G > A, c.8969G > A, c.9006A > T and c.9076C > T) or both of them simultaneously, splice site and enhancer (c.8486A > T and c.8487G > A). Interestingly, variant c.8488-2A > G (Figure 1C) replicated the splicing pattern of c.8488-1G > A (Figure 1B) although slight differences in the relative proportions of each isoform were observed [see Additional file 1, Figure S3]. Variants c.8486A > T and c.8487G > A of the penultimate and last nucleotides of exon 19, respectively, decreased the canonical donor splice site score of intron 19 (MaxEnt -480% and -502%, respectively) and NNSPLICE (0.95→0.87 and 0.95→0.4, respectively) but also disrupted one putative SF2/ASF enhancer. To ascertain the nature of these splicing anomalies two artificial variants were designed and generated: c.8484A > T (intact splice site and SF2 disruption) and c.8486A > C (weak alteration of splice site: 0.95 →0.92 and SF2 disruption; Additional file 1, Table S1). The first one, c.8484A > T, only revealed the wt isoform whereas c.8486A > C produced total exon 19 skipping [Additional file 1, Figure S3]. Consequently, only those variants that were predicted to affect the donor site, even slightly (c.8486A > T, c.8486A > C and c.8487G > A), altered the splicing process, suggesting that this is the causative mechanism. Actually, the last and penultimate exonic nucleotides of exons as well as the intronic positions +1 to +5 and -3 to -1 are also highly conserved  and should be considered potential targets of variants with impact on splicing . Finally, we also generated another five artificial variants, c.8512T > A and c.8518del (disruption of conserved SC35 and SRp55 ESEs), and three single and combined deletions of two conserved intron 19 motifs (c.8487+31_8487+42del, c.8488-49_8488-44del and c.[8487+31_8487+42del; c.8488-49_8488-44del]; Additional file 1, Table S1) that were genetically engineered in the respective minigenes but neither of them showed patent splicing anomalies.
Aberrant RNA isoforms and predicted effect on protein
A total of 14 different RNA isoforms, including the canonical transcripts, were detected in the minigene experiments of exons 19-20 (five isoforms) and 23-24 (nine isoforms; Additional file 1, Figure S2). The 12 aberrant isoforms (Figures 2 and 3) were caused by four different events including exon skipping (four isoforms: exons 19, 20, 23 and 24 skipping), double exon skipping (exons 23+24), intron 19 retention and use of alternative cryptic acceptor or donor splice sites (exon 20del12, exon 23ins2, exon 23del51, exon 24del7, exon 24del24 and exon 24del43). Moreover, most variants (12/19; Table 1) produced two or more distinct RNA isoforms. Likewise, the high resolution of the fluorescent RT-PCR technique in splicing reporter assays should be highlighted as it allowed accurate detection of isoforms that differed in size by as little as two (exon 23ins2) or seven nucleotides (exon 24del7) (Figure 3), which otherwise could be masked by the wt isoform in agarose gel electrophoresis . Furthermore, its high resolution is capable of detecting minor transcripts (less than 5% of total mRNA isoforms; Additional file 1, Figure S3) that could not be visualized in agarose gels. Finally, the minigene assay is a single-allele assay that allows a precise quantification of the different RNA isoforms without the interference of the normal allele as in lymphocyte RT-PCR. However, the Nonsense-Mediated mRNA Decay (NMD) mechanism selectively degrades mRNAs harboring premature termination codons (PTC) that can impair the relative proportions of each isoform. In fact, 13 out of 14 transcripts, including the wt ones (Figures 2 and 3), do not keep the ORF of exon 1 of pSPL3 and are, therefore, susceptible to undergo NMD unless this process is inhibited in cell culture , so special care should be taken in interpreting these results.
With regard to the putative pathogenicity of variants with an effect on splicing, two basic criteria are considered biological indicators [16–18]: first, total or almost total splicing disruption (absence of the canonical transcript); and second, the predicted effects of aberrant transcripts on protein translation, such as introduction of frameshifts and premature stop codons or loss of essential domains of BRCA2 (Table 1). Thirteen BIC variants of the splice sites (c.8486A > T, c.8487G > A, c.8487+1G > A, c.8487+3A > G, c.8488-2A > G, c.8954-3C > G, c.8954-1_8955delinsAA, c.9117G > A, c.9117+1G > T, c.9117+1G > A, c.9118-2A > G, c.9248_9256+7del and c.9256+1G > A) induced major splicing defects that consisted of exon skipping, double exon skipping, insertion of nucleotides and partial deletions of exons due to the use of cryptic acceptor or donor sites. Translation of the anomalous transcripts was predicted to result in protein truncation (exon 20 skipping, intron 19 retention, exon 23 skipping, exon 24 skipping, exon 23ins2, exon 24del7 and exon 24del43) or in-frame deletions (ex19 skipping, exon 20del12, exon 23del51, exons 23+24 skipping, exon 24del24) of conserved amino acids of the essential BRCA2 DNA binding domain (Table 1) which spans amino acids 2500 to 3098, suggesting that these variants might have a role in the disease. On the other hand, variant c.9256G > T of exon 24 produced a significant fraction of the canonical transcript (42.3%; Additional file 1, Figure S3). However, the canonical transcript contains the nonsense change r.9256g > u, thus supporting a double deleterious mechanism: protein truncation (p.G3086X) and splicing disruption (three aberrant isoforms; Table 1). Other variants, such as BRCA1 c.5123C > A and BRCA2 c.8168A > G (both missense changes), had already been shown to alter protein function and pre-mRNA processing [4, 5, 14, 44].
In contrast, the four SRE variants (ESE disruption and/or ESS creation; Table 1) had low to moderate effects on splicing since aberrant isoforms accounted for 8% (c.8969G > A) to 22% (c.9006A > T) of total transcripts. Up to now we have detected 16 SRE variants (two of them artificial variants)  of which five are placed in BRCA2 exon 5 suggesting that regulation of some exons strongly depends on supplementary cis-regulatory elements such as ESE or ESS. Thus, it is expected that each exon is regulated by different mechanisms that can only be determined experimentally. Another two putative silencer variants (c.9148C > T and c.9227G > T) showed minor splicing defects (< 2% of transcripts; data not shown) that were not considered relevant. Interestingly, the combination of two putative ESE variants, c.8972G > A (SRp55) and c.9006A > T (SRp40), in the same exon 23-24 minigene induced a greater splicing defect than each variant alone. In addition, variant c.9006A > T also strengthened the cryptic acceptor site 51 nucleotides downstream of the canonical one that is used in one aberrant isoform (exon 23del51; Figure 3). In single assays, c.8972G > A had no effect on splicing as previously reported  whereas c.9006A > T produced the canonical transcript (68.8%), exon 23del51 (19.7%), exon 23 skipping (9.2%) and exons 23+24 skipping (2.3%). Surprisingly, we found that the combination of both variants increased the fraction of the 51-nucleotide deletion isoform from 19.7% to 39.9% whereas the other anomalous isoforms decreased [See Additional file 1, Figures S3 and S4]. Furthermore, acceptor and donor sites of exon 23 are weak, suggesting that supplementary control elements are required for correct exon 23 recognition . Altogether these data suggest that exon 23 recognition depends, at least in part, on the splicing factors SRp55 and SRp40 that might act cooperatively although this issue should be confirmed by RNA binding assays. Conversely, other SRE variants were previously proven to trigger complete splicing defects, such as BRCA1 c.5080G > T and BRCA2 c.93G > T, c.145G > T or c.470_474del [14, 47].
Pathogenicity of variants with incomplete splicing outcomes is uncertain but they might constitute low-moderate penetrance alleles of breast/ovarian cancer as it occurs in CFTR-related disorders . Thus, the penetrance and severity of this group of diseases are correlated with an abnormal number of UG and U repeats located in the acceptor site of exon 9 of the CFTR gene, which increase the exon 9 skipping rate. Consequently, genetic variation can affect splicing efficiency that can modify the severity of the disease phenotype or be linked with disease susceptibility . In any case, additional epidemiological studies should be conducted to estimate accurately the breast/ovarian cancer risks of DNA variants with total and partial defective splicing outcomes.
Taking all these results together, 12.2% of all the variants (19/155) of BRCA2 exons 19, 20, 23 and 24 are associated with splicing defects. They comprise three nonsense, three missense, two synonymous and 11 splice site variants (including two deletions), indicating that any DNA change can disrupt pre-mRNA processing. Furthermore, in this and previous reports , we have analyzed in depth 14 exons of BRCA1 and BRCA2 (541 different BIC DNA variants) and detected 45 natural variants with an impact on splicing of which 21 altered splice sites, 14 modified SREs (enhancer disruption and/or silencer creation), eight created alternative sites and two affected the polypyrimidine tract. However, taking into account all the critical positions of the splice sites (intronic -3 to -1 and +1 to +5, and the conserved exonic nucleotides) of those 14 exons, the number of variants with impact on splicing would therefore increase to 85 (15.7% of all reported BIC variants) that would account for 33.9% of all presumed deleterious mutations (113 frameshift, 53 nonsense and 85 splicing variants), more than two-fold higher than the classical estimate .
Computational predictions of variants at the splice sites were precise by the three algorithms used in this work, but only HSF was able to identify all the cryptic alternative sites of exons 20, 23 and 24 (Figures 2 and 3). In contrast, bioinformatics evaluations of enhancer inactivation or silencer creation were not so accurate and had a high false-positive rate. A complete knowledge of the splicing code will enhance sensitivity of these bioinformatics predictions . However, the construction of one single prediction model seems to be an arduous task that will require comprehensive experimental validations .
Direct analysis of splicing anomalies in patient RNA should be the method of choice to identify variants with an effect on splicing but this approach presents two principal limitations. First, the patient sample might not be available and, second, leukocytes are the main source of RNA, so caution should be taken when interpreting these results since differential tissue-specific alternative splicing events could mask the real splicing outcome of a DNA variant [21, 34]. In this context, splicing reporter minigenes are valuable tools to corroborate bioinformatics data. This and previous reports have also shown its importance for the identification of variants with an effect on splicing without the need of patient RNA samples that are usually difficult to obtain [14, 18].
As a general rule, splicing profiles of patient and minigene RNA are very similar [14, 19, 49]. However, we have observed discrepancies between splicing anomalies of patient RNA and minigene assays, such as variants c.8488-1G > A (this work) and c.212+1G > A of BRCA1 exon 5 . The clinical relevance of c.8488-1G > A has been discussed above and that of c.212+1G > A is strongly supported by family, bioinformatics (donor site disruption) and functional data as well as by the fact that other variants of the donor site of exon5/intron 5 (c.211A > G and c.212+3A > G) were definitely classified as deleterious mutations by several approaches [19, 50, 51]. Therefore, it seems that the observed differences may probably be due to a technical issue of minigene constructs with particular exons where the lack of the natural genomic environment without its neighboring exons can be responsible for the different results obtained. Moreover, it was previously reported that a minigene containing only exon 37 of the NF1 gene with the pathogenic variant 6792C > G principally induced exon 37 skipping. Conversely, a minigene with exons 31 to 38 with the same variant replicated almost exactly the splicing pattern of patient lymphoblasts (canonical transcript, exon 37 skipping and exons 36+37 skipping) . Consequently, larger minigene constructions with more exons should be carried out in order to mimic the natural genomic background. The minigene system is a straightforward and robust assay that helps to classify DNA variants of unknown clinical significance under the splicing viewpoint, although these tests require further validation. As more data are collected, it will provide a more accurate risk estimation of breast and ovarian cancer associated with splicing alterations. Finally, analysis of minigene RT-PCR products in a DNA sequencer provides higher resolution than agarose electrophoresis [42, 52], since we have shown that a precise identification of minor or rare transcripts and quantification of all the isoforms generated by a specific DNA change is possible (Figures 2 and 3).
In conclusion, an important fraction of DNA variants are associated with splicing aberrations that should be considered as a primary mechanism of gene inactivation to be investigated in unclassified DNA variants. These studies provide insights into the basic regulatory mechanisms of this step of eukaryotic gene expression contributing to a better knowledge of the rules for exon definition. Hence, splicing functional assays supply essential information to distinguish between neutral variants and variants with an impact on splicing and should be incorporated in genotype screenings of human hereditary diseases.
The Breast Cancer Information Core
exonic splicing enhancer
exonic splicing silencer
Human Genome Variation Society
heterogeneous nuclear ribonucleoproteins
intronic splicing enhancer
Intronic splicing silencer
nonsense-mediated mRNA decay
- SR proteins:
splicing regulatory elements
Ford D, Easton DF, Stratton M, Narod S, Goldgar D, Devilee P, Bishop DT, Weber B, Lenoir G, Chang-Claude J, Sobol H, Teare MD, Struewing J, Arason A, Scherneck S, Peto J, Rebbeck TR, Tonin P, Neuhausen S, Barkardottir R, Eyfjord J, Lynch H, Ponder BA, Gayther SA, Zelada-Hedman M, the Breast Cancer Linkage Consortium: Genetic heterogeneity and penetrance analysis of the BRCA1 and BRCA2 genes in breast cancer families. The Breast Cancer Linkage Consortium. Am J Hum Genet. 1998, 62: 676-689. 10.1086/301749.
Stratton MR, Rahman N: The emerging landscape of breast cancer susceptibility. Nat Genet. 2008, 40: 17-22. 10.1038/ng.2007.53.
The Breast Cancer Information Core Database. (Accessed June 2011), [http://research.nhgri.nih.gov/projects/bic]
Vallon-Christersson J, Cayanan C, Haraldsson K, Loman N, Bergthorsson JT, Brondum-Nielsen K, Gerdes AM, Moller P, Kristoffersson U, Olsson H, Borg A, Monteiro AN: Functional analysis of BRCA1 C-terminal missense mutations identified in breast and ovarian cancer families. Hum Mol Genet. 2001, 10: 353-360. 10.1093/hmg/10.4.353.
Farrugia DJ, Agarwal MK, Pankratz VS, Deffenbaugh AM, Pruss D, Frye C, Wadum L, Johnson K, Mentlick J, Tavtigian SV, Goldgar DE, Couch FJ: Functional assays for classification of BRCA2 variants of uncertain significance. Cancer Res. 2008, 68: 3523-3531. 10.1158/0008-5472.CAN-07-1587.
Radice P, De Summa S, Caleca L, Tommasi S: Unclassified variants in BRCA genes: guidelines for interpretation. Ann Oncol. 2011, 22 (Suppl 1): i18-i23.
Kuznetsov SG, Liu P, Sharan SK: Mouse embryonic stem cell-based functional assay to evaluate mutations in BRCA2. Nat Med. 2008, 14: 875-881. 10.1038/nm.1719.
Goldgar DE, Easton DF, Deffenbaugh AM, Monteiro AN, Tavtigian SV, Couch FJ: Integrated evaluation of DNA sequence variants of unknown clinical significance: application to BRCA1 and BRCA2. Am J Hum Genet. 2004, 75: 535-544. 10.1086/424388.
Gomez Garcia EB, Oosterwijk JC, Timmermans M, van Asperen CJ, Hogervorst FB, Hoogerbrugge N, Oldenburg R, Verhoef S, Dommering CJ, Ausems MG, van Os TA, van der Hout AH, Ligtenberg M, van den OA, van der Luijt RB, Wijnen JT, Gille JJ, Lindsey PJ, Devilee P, Blok MJ, Vreeswijk MP: A method to assess the clinical significance of unclassified variants in the BRCA1 and BRCA2 genes based on cancer family history. Breast Cancer Res. 2009, 11: R8-10.1186/bcr2223.
Wang GS, Cooper TA: Splicing in disease: disruption of the splicing code and the decoding machinery. Nat Rev Genet. 2007, 8: 749-761. 10.1038/nrg2164.
Cartegni L, Chew SL, Krainer AR: Listening to silence and understanding nonsense: exonic mutations that affect splicing. Nat Rev Genet. 2002, 3: 285-298. 10.1038/nrg775.
Lopez-Bigas N, Audit B, Ouzounis C, Parra G, Guigo R: Are splicing mutations the most frequent cause of hereditary disease?. FEBS Lett. 2005, 579: 1900-1903. 10.1016/j.febslet.2005.02.047.
Krawczak M, Reiss J, Cooper DN: The mutational spectrum of single base-pair substitutions in mRNA splice junctions of human genes: causes and consequences. Hum Genet. 1992, 90: 41-54.
Sanz DJ, Acedo A, Infante M, Duran M, Perez-Cabornero L, Esteban-Cardenosa E, Lastra E, Pagani F, Miner C, Velasco EA: A high proportion of DNA variants of BRCA1 and BRCA2 is associated with aberrant splicing in breast/ovarian cancer patients. Clin Cancer Res. 2010, 16: 1957-1967. 10.1158/1078-0432.CCR-09-2564.
Pagani F, Stuani C, Tzetis M, Kanavakis E, Efthymiadou A, Doudounakis S, Casals T, Baralle FE: New type of disease causing mutations: the example of the composite exonic regulatory elements of splicing in CFTR exon 12. Hum Mol Genet. 2003, 12: 1111-1120. 10.1093/hmg/ddg131.
Brandao RD, van Roozendaal K, Tserpelis D, Garcia EG, Blok MJ: Characterisation of unclassified variants in the BRCA1/2 genes with a putative effect on splicing. Breast Cancer Res Treat. 2011, 129: 971-982. 10.1007/s10549-011-1599-7.
Menendez M, Castellsague J, Mirete M, Pros E, Feliubadalo L, Osorio A, Calaf M, Tornero E, Del Valle J, Fernandez-Rodriguez J, Quiles F, Salinas M, Velasco A, Teule A, Brunet J, Blanco I, Capella G, Lazaro C: Assessing the RNA effect of 26 DNA variants in the BRCA1 and BRCA2 genes. Breast Cancer Res Treat. 2012, 132: 979-992. 10.1007/s10549-011-1661-5.
Bonnet C, Krieger S, Vezain M, Rousselin A, Tournier I, Martins A, Berthet P, Chevrier A, Dugast C, Layet V, Rossi A, Lidereau R, Frebourg T, Hardouin A, Tosi M: Screening BRCA1 and BRCA2 unclassified variants for splicing mutations using reverse transcription PCR on patient RNA and an ex vivo assay based on a splicing reporter minigene. J Med Genet. 2008, 45: 438-446. 10.1136/jmg.2007.056895.
Thery JC, Krieger S, Gaildrat P, Revillion F, Buisine MP, Killian A, Duponchel C, Rousselin A, Vaur D, Peyrat JP, Berthet P, Frebourg T, Martins A, Hardouin A, Tosi M: Contribution of bioinformatics predictions and functional splicing assays to the interpretation of unclassified variants of the BRCA genes. Eur J Hum Genet. 2011, 19: 1052-1058. 10.1038/ejhg.2011.100.
Thomassen M, Blanco A, Montagna M, Hansen TV, Pedersen IS, Gutierrez-Enriquez S, Menendez M, Fachal L, Santamarina M, Steffensen AY, Jonson L, Agata S, Whiley P, Tognazzo S, Tornero E, Jensen UB, Balmana J, Kruse TA, Goldgar DE, Lazaro C, Diez O, Spurdle AB, Vega A: Characterization of BRCA1 and BRCA2 splicing variants: a collaborative report by ENIGMA consortium members. Breast Cancer Res Treat. 2012, 132: 1009-1023. 10.1007/s10549-011-1674-0.
Baralle D, Lucassen A, Buratti E: Missed threads. The impact of pre-mRNA splicing defects on clinical practice. EMBO Rep. 2009, 10: 810-816. 10.1038/embor.2009.170.
Cooper TA: Use of minigene systems to dissect alternative splicing elements. Methods. 2005, 37: 331-340. 10.1016/j.ymeth.2005.07.015.
Infante M, Duran M, Esteban-Cardenosa E, Miner C, Velasco E: High proportion of novel mutations of BRCA1 and BRCA2 in breast/ovarian cancer patients from Castilla-Leon (central Spain). J Hum Genet. 2006, 51: 611-617. 10.1007/s10038-006-0404-7.
Esteban-Cardenosa E, Duran M, Infante M, Velasco E, Miner C: High-throughput mutation detection method to scan BRCA1 and BRCA2 based on heteroduplex analysis by capillary array electrophoresis. Clin Chem. 2004, 50: 313-320. 10.1373/clinchem.2003.023614.
Velasco E, Infante M, Duran M, Perez-Cabornero L, Sanz DJ, Esteban-Cardenosa E, Miner C: Heteroduplex analysis by capillary array electrophoresis for rapid mutation detection in large multiexon genes. Nat Protoc. 2007, 2: 237-246.
Human Genome Variation Society: Guidelines for mutation nomenclature. (Accessed September 2011), [http://www.hgvs.org/rec.html]
Reese MG, Eeckman FH, Kulp D, Haussler D: Improved splice site detection in Genie. J Comput Biol. 1997, 4: 311-323. 10.1089/cmb.1997.4.311.
Desmet FO, Hamroun D, Lalande M, Collod-Beroud G, Claustres M, Beroud C: Human Splicing Finder: an online bioinformatics tool to predict splicing signals. Nucleic Acids Res. 2009, 37: e67-10.1093/nar/gkp215.
Pettigrew C, Wayte N, Lovelock PK, Tavtigian SV, Chenevix-Trench G, Spurdle AB, Brown MA: Evolutionary conservation analysis increases the colocalization of predicted exonic splicing enhancers in the BRCA1 gene with missense sequence changes and in-frame deletions, but not polymorphisms. Breast Cancer Res. 2005, 7: R929-R939. 10.1186/bcr1324.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.
Soule HD, Maloney TM, Wolman SR, Peterson WD, Brenz R, McGrath CM, Russo J, Pauley RJ, Jones RF, Brooks SC: Isolation and characterization of a spontaneously immortalized human breast epithelial cell line, MCF-10. Cancer Res. 1990, 50: 6075-6086.
METHODBOOK: PCR mutagenesis. (Accessed January 2010), [http://www.methodbook.net/pcr/pcrmut.html]
Howlett NG, Taniguchi T, Olson S, Cox B, Waisfisz Q, De Die-Smulders C, Persky N, Grompe M, Joenje H, Pals G, Ikeda H, Fox EA, D'Andrea AD: Biallelic inactivation of BRCA2 in Fanconi anemia. Science. 2002, 297: 606-9. 10.1126/science.1073834.
Buratti E, Baralle M, Baralle FE: Defective splicing, disease and therapy: searching for master checkpoints in exon definition. Nucleic Acids Res. 2006, 34: 3494-3510. 10.1093/nar/gkl498.
Ast G: How did alternative splicing evolve?. Nat Rev Genet. 2004, 5: 773-782.
Kim E, Goren A, Ast G: Insights into the connection between cancer and alternative splicing. Trends Genet. 2008, 24: 7-10. 10.1016/j.tig.2007.10.001.
Michael IP, Kurlender L, Memari N, Yousef GM, Du D, Grass L, Stephan C, Jung K, Diamandis EP: Intron retention: a common splicing event within the human kallikrein gene family. Clin Chem. 2005, 51: 506-515. 10.1373/clinchem.2004.042341.
The International Agency for Research on Cancer (IARC): Alignment of the BRCA2 protein (12 species). (Accessed September 2011), [http://agvgd.iarc.fr/BRCA2_Align.htm]
Infante M, Duran M, Acedo A, Perez-Cabornero L, Sanz DJ, Garcia-Gonzalez M, Beristain E, Esteban-Cardenosa E, de la Hoya M, Teule A, Vega A, Tejada MI, Lastra E, Miner C, Velasco EA: BRCA1 5272-1G > A and BRCA2 5374delTATG are founder mutations of high relevance for genetic counselling in breast/ovarian cancer families of Spanish origin. Clin Genet. 2010, 77: 60-69. 10.1111/j.1399-0004.2009.01272.x.
Zhang C, Li WH, Krainer AR, Zhang MQ: RNA landscape of evolution for optimal exon and intron discrimination. Proc Natl Acad Sci USA. 2008, 105: 5797-5802. 10.1073/pnas.0801692105.
Whiley PJ, Guidugli L, Walker LC, Healey S, Thompson BA, Lakhani SR, Da Silva LM, Tavtigian SV, Goldgar DE, Brown MA, Couch FJ, Spurdle AB: 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. Hum Mutat. 2011, 32: 678-687. 10.1002/humu.21495.
Hinzpeter A, Aissat A, Sondo E, Costa C, Arous N, Gameiro C, Martin N, Tarze A, Weiss L, de Becdelievre A, Costes B, Goossens M, Galietta LJ, Girodon E, Fanen P: Alternative splicing at a NAGNAG acceptor site as a novel phenotype modifier. PLoS Genet. 2010, 6: e100153-
Perrin-Vidoz L, Sinilnikova OM, Stoppa-Lyonnet D, Lenoir GM, Mazoyer S: The nonsense-mediated mRNA decay pathway triggers degradation of most BRCA1 mRNAs bearing premature termination codons. Hum Mol Genet. 2002, 11: 2805-2814. 10.1093/hmg/11.23.2805.
Millevoi S, Bernat S, Telly D, Fouque F, Gladieff L, Favre G, Vagner S, Toulas C: The c.5242C > A BRCA1 missense variant induces exon skipping by increasing splicing repressors binding. Breast Cancer Res Treat. 2010, 120: 391-399. 10.1007/s10549-009-0392-3.
Whiley PJ, Pettigrew CA, Brewster BL, Walker LC, Spurdle AB, Brown MA: Effect of BRCA2 sequence variants predicted to disrupt exonic splice enhancers on BRCA2 transcripts. BMC Med Genet. 2010, 11: 80-
Hertel KJ: Combinatorial control of exon recognition. J Biol Chem. 2008, 283: 1211-1215.
Goina E, Skoko N, Pagani F: Binding of DAZAP1 and hnRNPA1/A2 to an exonic splicing silencer in a natural BRCA1 exon 18 mutant. Mol Cell Biol. 2008, 28: 3850-3860. 10.1128/MCB.02253-07.
Steiner B, Truninger K, Sanz J, Schaller A, Gallati S: The role of common single-nucleotide polymorphisms on exon 9 and exon 12 skipping in nonmutated CFTR alleles. Hum Mutat. 2004, 24: 120-129. 10.1002/humu.20064.
Gaildrat P, Krieger S, Thery JC, Killian A, Rousselin A, Berthet P, Frebourg T, Hardouin A, Martins A, Tosi M: The BRCA1 c.5434C- > G (p.Pro1812Ala) variant induces a deleterious exon 23 skipping by affecting exonic splicing regulatory elements. J Med Genet. 2010, 47: 398-403. 10.1136/jmg.2009.074047.
Milne RL, Osorio A, Cajal TR, Vega A, Llort G, de la Hoya M, Diez O, Alonso MC, Lazaro C, Blanco I, Sanchez-de-Abajo A, Caldes T, Blanco A, Grana B, Duran M, Velasco E, Chirivella I, Cardenosa EE, Tejada MI, Beristain E, Miramar MD, Calvo MT, Martinez E, Guillen C, Salazar R, San RC, Antoniou AC, Urioste M, Benitez J: The average cumulative risks of breast and ovarian cancer for carriers of mutations in BRCA1 and BRCA2 attending genetic counseling units in Spain. Clin Cancer Res. 2008, 14: 2861-2869. 10.1158/1078-0432.CCR-07-4436.
Easton DF, Deffenbaugh AM, Pruss D, Frye C, Wenstrup RJ, Allen-Brady K, Tavtigian SV, Monteiro AN, Iversen ES, Couch FJ, Goldgar DE: A systematic genetic assessment of 1,433 sequence variants of unknown clinical significance in the BRCA1 and BRCA2 breast cancer-predisposition genes. Am J Hum Genet. 2007, 81: 873-883. 10.1086/521032.
Dosil V, Tosar A, Canadas C, Perez-Segura P, Diaz-Rubio E, Caldes T, de la Hoya M: Alternative splicing and molecular characterization of splice site variants: BRCA1 c.591C > T as a case study. Clin Chem. 2010, 56: 53-61. 10.1373/clinchem.2009.132274.
We are grateful to the breast/ovarian cancer patients and clinicians who collaborated in this study, particularly Dr. Enrique Lastra (Complejo Hospitalario de Burgos). This work was supported by grants PI10/2910 of the Instituto de Salud Carlos III (Spanish Ministry of Science and Innovation), and CSI004A10-2 (Consejería de Educación) and BIO39/VA27/10 (Consejería de Sanidad) of the Regional Government of Castilla y León. Alberto Acedo is supported by the European Social Fund and Consejeria de Educación de la Junta de Castilla y León under the P.O. Castilla y León 2007-2013.We are also grateful to Cristina Hernández Moro and Beatriz Díez Gómez for critical reading of the manuscript. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
The authors declare that they have no competing interests.
AA and DJS participated in the experimental design and bioinformatics analysis, and performed the splicing functional assays. MI and LP-C participated in the functional assays and bioinformatics analyses. MD and CM provided the breast/ovarian cancer patients included in this study as well as mutation data of patients from Castilla y León (Spain). EAV designed the study, conducted the bioinformatics analysis, supervised minigene construction and data analysis, and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Bioinformatics and RNA analyses of candidate DNA variants. Bioinformatics analysis of putative splicing variants of exons 19, 20, 23 and 24 of BRCA2 (Table S1). Splicing outcomes of DNA variant c.8488-1G > A (Figure S1). Sequence chromatograms of the main splicing outcomes of the minigenes 19-20 and 23-24 of BRCA2 (Figure S2). Quantification of RNA isoforms induced by variants of minigenes of exons 19-20 and 23-24 (Figure S3). Cumulative effect of ESE mutations from BRCA2 exon 23 on splicing (Figure S4). (PDF 854 KB)
About this article
Cite this article
Acedo, A., Sanz, D.J., Durán, M. et al. Comprehensive splicing functional analysis of DNA variants of the BRCA2 gene by hybrid minigenes. Breast Cancer Res 14, R87 (2012). https://doi.org/10.1186/bcr3202