Characterisation of unclassified variants in the BRCA1/2 genes with a putative effect on splicing

Synonymous Variants: Necessary Nuance in our Understanding of Cancer Drivers and Treatment Outcomes

Once called "silent mutations" and assumed to have no effect on protein structure and function, synonymous variants are now recognized to be drivers for some cancers. There have been significant advances in our understanding of the numerous mechanisms by which synonymous single nucleotide variants (sSNVs) can affect protein structure and function by affecting pre-mRNA splicing, mRNA expression, stability, folding, miRNA binding, translation kinetics, and co-translational folding. This review highlights the need for considering sSNVs in cancer biology to gain a better understanding of the genetic determinants of human cancers and to improve their diagnosis and treatment. We surveyed the literature for reports of sSNVs in cancer and found numerous studies on the consequences of sSNVs on gene function with supporting in vitro evidence. We also found reports of sSNVs that have statistically significant associations with specific cancer types but for which in vitro studies are lacking to support the reported associations. Additionally, we found reports of germline and somatic sSNVs that were observed in numerous clinical studies and for which in silico analysis predicts possible effects on gene function. We provide a review of these investigations and discuss necessary future studies to elucidate the mechanisms by which sSNVs disrupt protein function and are play a role in tumorigeneses, cancer progression, and treatment efficacy. As splicing dysregulation is one of the most well recognized mechanisms by which sSNVs impact protein function, we also include our own in silico analysis for predicting which sSNVs may disrupt pre-mRNA splicing.

Pathogenic neurofibromatosis type 1 (NF1) RNA splicing resolved by targeted RNAseq

Neurofibromatosis type 1 (NF1) is caused by loss-of-function variants in the NF1 gene. Approximately 10% of these variants affect RNA splicing and are either missed by conventional DNA diagnostics or are misinterpreted by in silico splicing predictions. Therefore, a targeted RNAseq-based approach was designed to detect pathogenic RNA splicing and associated pathogenic DNA variants. For this method RNA was extracted from lymphocytes, followed by targeted RNAseq. Next, an in-house developed tool (QURNAs) was used to calculate the enrichment score (ERS) for each splicing event. This method was thoroughly tested using two different patient cohorts with known pathogenic splice-variants in NF1. In both cohorts all 56 normal reference transcript exon splice junctions, 24 previously described and 45 novel non-reference splicing events were detected. Additionally, all expected pathogenic splice-variants were detected. Eleven patients with NF1 symptoms were subsequently tested, three of which have a known NF1 DNA variant with a putative effect on RNA splicing. This effect could be confirmed for all 3. The other eight patients were previously without any molecular confirmation of their NF1-diagnosis. A deep-intronic pathogenic splice variant could now be identified for two of them (25%). These results suggest that targeted RNAseq can be successfully used to detect pathogenic RNA splicing variants in NF1.

Interpretation of BRCA2 Splicing Variants: A Case Series of Challenging Variant Interpretations and the Importance of Functional RNA Analysis

A substantial proportion of pathogenic variants associated with an increased risk of hereditary cancer are sequence variants affecting RNA splicing. The classification of these variants can be complex when both non-functional and functional transcripts are produced from the variant allele. We present four BRCA2 splice site variants with complex variant interpretations (BRCA2 c.68-3T>G, c.68-2A>G, c.425G>T, c.8331+2T>C). Evidence supporting a pathogenic classification is available for each variant, including in silico models, absence in population databases, and published functional data. However, comprehensive RNA analysis showed that some functional transcript may be produced by each variant. BRCA2 c.68-3T>G results in a partial splice defect. For BRCA2 c.68-2A>G and c.425G>T, aberrant splicing was shown to produce a potentially functional, in-frame transcript. BRCA2 c.8331+2T>C may utilize a functional GC donor in place of the wild-type GT donor. The severity of cancer history for carriers of these variants was also assessed using a history weighting algorithm and was not consistent with pathogenic controls (carriers of known pathogenic variants in BRCA2). Due to the conflicting evidence, our laboratory classifies these BRCA2 variants as variants of uncertain significance. This highlights the importance of evaluating new and existing evidence to ensure accurate variant classification and appropriate patient care.

A non-invasive diagnostic assay for rapid detection and characterization of aberrant mRNA-splicing by nonsense mediated decay inhibition

BACKGROUND

Interpretation of genetic variants detected by sequencing of genomic DNA, which may cause splicing defects, regularly requires mRNA analysis. Usually, only bioinformatic testing is provided, because simple and non-invasive assay protocols are lacking. Furthermore, the detection of mis-splicing is often hampered by nonsense mediated mRNA decay (NMD).

METHODS

Starting from a case of Pompe disease with two potential splicing variants an assay for the analysis of splice defects in general was developed. We analyzed the transcripts from the gene of interest by standard methods after short-term culture of the patient's lymphocytes in the presence and absence of a NMD inhibitor. Variant and wild type transcript expression were quantified by allele specific PCR in the patient and both parents and the expression ratio with/without NMD inhibition was calculated for each transcript.

RESULTS

NMD detection in lymphocytes was optimized and evaluated by analyzing a naturally occurring NMD transcript. Several compounds inhibited NMD successfully, including potential therapeutic agents. Sample storage for up to 4 days at room temperature prior to lymphocyte isolation did not affect results. In a proof of concept we identified two candidate variants as severe splicing variants in a patient with Pompe disease, but the strategy can also be used to screen for any mis-spliced transcripts prone to NMD.

CONCLUSIONS

We developed a simple, non-invasive assay for the detection and characterization of potential splicing variants. This is essential, because early and near-term diagnosis and disease classification is required to facilitate therapy in many genetic diseases.

Skipping Nonsense to Maintain Function: The Paradigm of BRCA2 Exon 12

Germline nonsense and canonical splice site variants identified in disease-causing genes are generally considered as loss-of-function (LoF) alleles and classified as pathogenic. However, a fraction of such variants could maintain function through their impact on RNA splicing. To test this hypothesis, we used the alternatively spliced BRCA2 exon 12 (E12) as a model system because its in-frame skipping leads to a potentially functional protein. All E12 variants corresponding to putative LoF variants or predicted to alter splicing (n = 40) were selected from human variation databases and characterized for their impact on splicing in minigene assays and, when available, in patient lymphoblastoid cell lines. Moreover, a selection of variants was analyzed in a mouse embryonic stem cell-based functional assay. Using these complementary approaches, we demonstrate that a subset of variants, including nonsense variants, induced in-frame E12 skipping through the modification of splice sites or regulatory elements and, consequently, led to an internally deleted but partially functional protein. These data provide evidence, for the first time in a cancer-predisposition gene, that certain presumed null variants can retain function due to their impact on splicing. Further studies are required to estimate cancer risk associated with these hypomorphic variants. More generally, our findings highlight the need to exercise caution in the interpretation of putative LoF variants susceptible to induce in-frame splicing modifications. SIGNIFICANCE: This study presents evidence that certain presumed loss-of-function variants in a cancer predisposition gene can retain function due to their direct impact on RNA splicing.

The Missing "lnc" between Genetics and Cardiac Disease

Cardiovascular disease (CVD) is one of the biggest threats to public health worldwide. Identifying key genetic contributors to CVD enables clinicians to assess the most effective treatment course and prognosis, as well as potentially inform family members. This often involves either whole exome sequencing (WES) or targeted panel analysis of known pathogenic genes. In the future, tailored or personalized therapeutic strategies may be implemented, such as gene therapy. With the recent revolution in deep sequencing technologies, we know that up to 90% of the human genome is transcribed, despite only 2% of the 6 billion DNA bases coding for proteins. The long non-coding RNA (lncRNA) “genes” make up an important and significant fraction of this “dark matter” of the genome. We highlight how, despite lncRNA genes exceeding that of classical protein-coding genes by number, the “non-coding” human genome is neglected when looking for genetic components of disease. WES platforms and pathogenic gene panels still do not cover even characterized lncRNA genes that are functionally involved in the pathophysiology of CVD. We suggest that the importance of lncRNAs in disease causation and progression be taken as seriously as that of pathogenic protein variants and mutations, and that this is maybe a new area of attention for clinical geneticists.

Systematic misclassification of missense variants in BRCA1 and BRCA2 "coldspots"

Purpose

Guidelines for variant interpretation incorporate variant hotspots in critical functional domains as evidence for pathogenicity (e.g., PM1 and PP2), but do not use “coldspots,” that is, regions without essential functions that tolerate variation, as evidence a variant is benign. To improve variant classification we evaluated BRCA1 and BRCA2 missense variants reported in ClinVar to identify regions where pathogenic missenses are extremely infrequent, defined as coldspots.

Methods

We used Bayesian approaches to model variant classification in these regions.

Results

BRCA1 exon 11 (~60% of the coding sequence), and BRCA2 exons 10 and 11 (~65% of the coding sequence), are coldspots. Of 89 pathogenic (P) or likely pathogenic (LP) missense variants in BRCA1, none are in exon 11 (odds <0.01, 95% confidence interval [CI] 0.0–0.01). Of 34 P or LP missense variants in BRCA2, none are in exons 10–11 (odds <0.01, 95% CI 0.0–0.01). More than half of reported missense variants of uncertain significance (VUS) in BRCA1 and BRCA2 are in coldspots (3115/5301 = 58.8%). Reclassifying these 3115 VUS as likely benign would substantially improve variant classification.

Conclusion

In BRCA1 and BRCA2 coldspots, missense variants are very unlikely to be pathogenic. Classification schemes that incorporate coldspots can reduce the number of VUS and mitigate risks from reporting benign variation as VUS.

Massive computational identification of somatic variants in exonic splicing enhancers using The Cancer Genome Atlas

Owing to the development of next-generation sequencing (NGS) technologies, a large number of somatic variants have been identified in various types of cancer. However, the functional significance of most somatic variants remains unknown. Somatic variants that occur in exonic splicing enhancer (ESE) regions are thought to prevent serine and arginine-rich (SR) proteins from binding to ESE sequence motifs, which leads to exon skipping. We computationally identified somatic variants in ESEs by compiling numerous open-access datasets from The Cancer Genome Atlas (TCGA). Using somatic variants and RNA-seq data from 9635 patients across 32 TCGA projects, we identified 646 ESE-disrupting variants. The false positive rate of our method, estimated using a permutation test, was approximately 1%. Of these ESE-disrupting variants, approximately 71% were located in the binding motifs of four classical SR proteins. ESE-disrupting variants occurred in proportion to the number of somatic variants, but not necessarily in the specific genes associated with the biological processes of cancer. Existing bioinformatics tools could not predict the pathogenicity of ESE-disrupting variants identified in this study, although these variants could cause exon skipping. We demonstrated that ESE-disrupting nonsense variants tended to escape nonsense-mediated decay surveillance. Using integrated analyses of open access data, we could specifically identify ESE-disrupting variants. We have generated a powerful tool, which can handle datasets without normal samples or raw data, and thus contribute to reducing variants of uncertain significance because our statistical approach only uses the exon-junction read counts from the tumor samples.

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.

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.

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.

Next-Generation Sequencing in Oncology: Genetic Diagnosis, Risk Prediction and Cancer Classification

Next-generation sequencing (NGS) technology has expanded in the last decades with significant improvements in the reliability, sequencing chemistry, pipeline analyses, data interpretation and costs. Such advances make the use of NGS feasible in clinical practice today. This review describes the recent technological developments in NGS applied to the field of oncology. A number of clinical applications are reviewed, i.e., mutation detection in inherited cancer syndromes based on DNA-sequencing, detection of spliceogenic variants based on RNA-sequencing, DNA-sequencing to identify risk modifiers and application for pre-implantation genetic diagnosis, cancer somatic mutation analysis, pharmacogenetics and liquid biopsy. Conclusive remarks, clinical limitations, implications and ethical considerations that relate to the different applications are provided.

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.

BRCA1 EXON 11, a CERES (composite regulatory element of splicing) element involved in splice regulation

Unclassified variants (UV) of BRCA1 can affect normal pre-mRNA splicing. Here, we investigate the UV c.693G>A, a "silent" change in BRCA1 exon 11, which we have found induces aberrant splicing in patient carriers and in vitro. Using a minigene assay, we show that the UV c.693G>A has a strong effect on the splicing isoform ratio of BRCA1. Systematic site-directed mutagenesis of the area surrounding the nucleotide position c.693G>A induced variable changes in the level of exon 11 inclusion/exclusion in the mRNA, pointing to the presence of a complex regulatory element with overlapping enhancer and silencer functions. Accordingly, protein binding analysis in the region detected several splicing regulatory factors involved, including SRSF1, SRSF6 and SRSF9, suggesting that this sequence represents a composite regulatory element of splicing (CERES).

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 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.

Evolutionary constraint helps unmask a splicing regulatory region in BRCA1 exon 11

Background

Alternative splicing across exon 11 produces several BRCA1 isoforms. Their proportion varies during the cell cycle, between tissues and in cancer suggesting functional importance of BRCA1 splicing regulation around this exon. Although the regulatory elements driving exon 11 splicing have never been identified, a selective constraint against synonymous substitutions (silent nucleotide variations that do not alter the amino acid residue sequence) in a critical region of BRCA1 exon 11 has been reported to be associated with the necessity to maintain regulatory sequences.

Methodology/Principal Findings

Here we have designed a specific minigene to investigate the possibility that this bias in synonymous codon usage reflects the need to preserve the BRCA1 alternative splicing program. We report that in-frame deletions and translationally silent nucleotide substitutions in the critical region affect splicing regulation of BRCA1 exon 11.

Conclusions/Significance

Using a hybrid minigene approach, we have experimentally validated the hypothesis that the need to maintain correct alternative splicing is a selective pressure against translationally silent sequence variations in the critical region of BRCA1 exon 11. Identification of the trans-acting factors involved in regulating exon 11 alternative splicing will be important in understanding BRCA1-associated tumorigenesis.

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.

Characterization of BRCA1 and BRCA2 splicing variants: a collaborative report by ENIGMA consortium members

Mutations in BRCA1 and BRCA2 predispose carriers to early onset breast and ovarian cancer. A common problem in clinical genetic testing is interpretation of variants with unknown clinical significance. The Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium was initiated to evaluate and implement strategies to characterize the clinical significance of BRCA1 and BRCA2 variants. As an initial project of the ENIGMA Splicing Working Group, we report splicing and multifactorial likelihood analysis of 25 BRCA1 and BRCA2 variants from seven different laboratories. Splicing analysis was performed by reverse transcriptase PCR or mini gene assay, and sequencing to identify aberrant transcripts. The findings were compared to bioinformatic predictions using four programs. The posterior probability of pathogenicity was estimated using multifactorial likelihood analysis, including co-occurrence with a deleterious mutation, segregation and/or report of family history. Abnormal splicing patterns expected to lead to a non-functional protein were observed for 7 variants (BRCA1 c.441+2T>A, c.4184_4185+2del, c.4357+1G>A, c.4987-2A>G, c.5074G>C, BRCA2 c.316+5G>A, and c.8754+3G>C). Combined interpretation of splicing and multifactorial analysis classified an initiation codon variant (BRCA2 c.3G>A) as likely pathogenic, uncertain clinical significance for 7 variants, and indicated low clinical significance or unlikely pathogenicity for another 10 variants. Bioinformatic tools predicted disruption of consensus donor or acceptor sites with high sensitivity, but cryptic site usage was predicted with low specificity, supporting the value of RNA-based assays. The findings also provide further evidence that clinical RNA-based assays should be extended from analysis of invariant dinucleotides to routinely include all variants located within the donor and acceptor consensus splicing sites. Importantly, this study demonstrates the added value of collaboration between laboratories, and across disciplines, to collate and interpret information from clinical testing laboratories to consolidate patient management.