Naturally occurring BRCA2 alternative mRNA splicing events in clinically relevant samples

Systematic Minigene-Based Splicing Analysis and Tentative Clinical Classification of 52 CHEK2 Splice-Site Variants

BACKGROUND

Disrupted pre-mRNA splicing is a frequent deleterious mechanism in hereditary cancer. We aimed to functionally analyze candidate spliceogenic variants of the breast cancer susceptibility gene CHEK2 by splicing reporter minigenes.

METHODS

A total of 128 CHEK2 splice-site variants identified in the Breast Cancer After Diagnostic Gene Sequencing (BRIDGES) project (https://cordis.europa.eu/project/id/634935) were analyzed with MaxEntScan and subsetted to 52 variants predicted to impact splicing. Three CHEK2 minigenes, which span all 15 exons, were constructed and validated. The 52 selected variants were then genetically engineered into the minigenes and assayed in MCF-7 (human breast adenocarcinoma) cells.

RESULTS

Of 52 variants, 46 (88.5%) impaired splicing. Some of them led to complex splicing patterns with up to 11 different transcripts. Thirty-four variants induced splicing anomalies without any trace or negligible amounts of the full-length transcript. A total of 89 different transcripts were annotated, which derived from different events: single- or multi-exon skipping, alternative site-usage, mutually exclusive exon inclusion, intron retention or combinations of the abovementioned events. Fifty-nine transcripts were predicted to introduce premature termination codons, 7 kept the original open-reading frame, 5 removed the translation start codon, 6 affected the 5'UTR (Untranslated Region), and 2 included missense variations. Analysis of variant c.684-2A > G revealed the activation of a non-canonical TG-acceptor site and exon 6 sequences critical for its recognition.

CONCLUSIONS

Incorporation of minigene read-outs into an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme allowed us to classify 32 CHEK2 variants (27 pathogenic/likely pathogenic and 5 likely benign). However, 20 variants (38%) remained of uncertain significance, reflecting in part the complex splicing patterns of this gene.

Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

Analysis of BRCA1 and BRCA2 alternative splicing in predisposition to ovarian cancer

BACKGROUND

The mRNA splicing is regulated on multiple levels, resulting in the proper distribution of genes' transcripts in each cell and maintaining cell homeostasis. At the same time, the expression of alternative transcripts can change in response to underlying genetic variants, often missed during routine diagnostics.

AIM

The main aim of this study was to define the frequency of aberrant splicing in BRCA1 and BRCA2 genes in blood RNA extracted from ovarian cancer patients who were previously found negative for the presence of pathogenic alterations in the 25 most commonly analysed ovarian cancer genes, including BRCA1 and BRCA2.

MATERIAL AND METHODS

Frequency and spectrum of splicing alterations in BRCA1 and BRCA2 genes were analysed in blood RNA from 101 ovarian cancer patients and healthy controls (80 healthy women) using PCR followed by gel electrophoresis and Sanger sequencing. The expression of splicing events was examined using RT-qPCR.

RESULTS

We did not identify any novel, potentially pathogenic splicing alterations. Nevertheless, we detected six naturally occurring transcripts, named BRCA1ΔE9-10, BRCA1ΔE11, BRCA1ΔE11q, and BRCA2ΔE3, BRCA2ΔE12 and BRCA2ΔE17-18 of which three (BRCA1ΔE11q, BRCA1ΔE11 and BRCA2ΔE3) were significantly higher expressed in the ovarian cancer cohort than in healthy controls (p ≤ 0.0001).

CONCLUSIONS

This observation indicates that the upregulation of selected naturally occurring transcripts can be stimulated by non-genetic mechanisms and be a potential systemic response to disease progression and/or treatment. However, this hypothesis requires further examination.

APPLICATION OF THE ACMG/AMP FRAMEWORK TO CAPTURE EVIDENCE RELEVANT TO PREDICTED AND OBSERVED IMPACT ON SPLICING: RECOMMENDATIONS FROM THE CLINGEN SVI SPLICING SUBGROUP

The American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1 (null variant in a gene where loss-of-function is the mechanism of disease), PS3 (functional assays show damaging effect on splicing), PP3 (computational evidence supports a splicing effect), BS3 (functional assays show no damaging effect on splicing), BP4 (computational evidence suggests no splicing impact), and BP7 (silent change with no predicted impact on splicing). However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation (SVI) Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. Our study utilised empirically derived splicing evidence to: 1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, 2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and 3) exemplify methodology to calibrate bioinformatic splice prediction tools. We propose repurposing of the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely BP7 may be used to capture RNA results demonstrating no impact on splicing for both intronic and synonymous variants, and for missense variants if protein functional impact has been excluded. Furthermore, we propose that the PS3 and BS3 codes are applied only for well-established assays that measure functional impact that is not directly captured by RNA splicing assays. We recommend the application of PS1 based on similarity of predicted RNA splicing effects for a variant under assessment in comparison to a known Pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA assay evidence described aim to help standardise variant pathogenicity classification processes and result in greater consistency when interpreting splicing-based evidence.

Identification of a large intra-exonic deletion in BRCA2 exon 18 in a pancreatic ductal adenocarcinoma

By 2030, pancreatic cancer will become the second leading cause of cancer-related deaths in the United States and in Europe. The management of patients with advanced pancreatic cancer relies on chemotherapy and poly (ADP-ribose) polymerase inhibitors for patients who carry BRCA1/2 inactivating alterations. Some variants, such as large insertion/deletions (Indels), inactivating BRCA1/2 and therefore of clinical relevance can be hard to detect by next-generation sequencing techniques. Here we report a 47-year-old patient presenting with pancreatic cancer whose tumour harbours a large somatic intra-exonic deletion of BRCA2 of 141 bp. This BRCA2 deletion, located in the C-terminal domain, can be considered as pathogenic and consequently affect tumorigenesis because it is involved in the interaction between the DSS1 protein and DNA. Thanks to the optimized bioinformatics algorithm, this intermediate size deletion in BRCA2 was identified, enabling personalized patient management via the inclusion of the patients in a clinical trial.

Clinical, splicing, and functional analysis to classify BRCA2 exon 3 variants: Application of a points-based ACMG/AMP approach

Skipping of BRCA2 exon 3 (∆E3) is a naturally occurring splicing event, complicating clinical classification of variants that may alter ∆E3 expression. This study used multiple evidence types to assess pathogenicity of 85 variants in/near BRCA2 exon 3. Bioinformatically predicted spliceogenic variants underwent mRNA splicing analysis using minigenes and/or patient samples. ∆E3 was measured using quantitative analysis. A mouse embryonic stem cell (mESC) based assay was used to determine the impact of 18 variants on mRNA splicing and protein function. For each variant, population frequency, bioinformatic predictions, clinical data, and existing mRNA splicing and functional results were collated. Variant class was assigned using a gene-specific adaptation of ACMG/AMP guidelines, following a recently proposed points-based system. mRNA and mESC analysis combined identified six variants with transcript and/or functional profiles interpreted as loss of function. Cryptic splice site use for acceptor site variants generated a transcript encoding a shorter protein that retains activity. Overall, 69/85 (81%) variants were classified using the points-based approach. Our analysis shows the value of applying gene-specific ACMG/AMP guidelines using a points-based approach and highlights the consideration of cryptic splice site usage to appropriately assign PVS1 code strength.

Fanconi Anaemia, Childhood Cancer and the BRCA Genes

Fanconi anaemia (FA) is an inherited chromosomal instability disorder characterised by congenital and developmental abnormalities and a strong cancer predisposition. In less than 5% of cases FA can be caused by bi-allelic pathogenic variants (PGVs) in BRCA2/FANCD1 and in very rare cases by bi-allelic PGVs in BRCA1/FANCS. The rarity of FA-like presentation due to PGVs in BRCA2 and even more due to PGVs in BRCA1 supports a fundamental role of the encoded proteins for normal development and prevention of malignant transformation. While FA caused by BRCA1/2 PGVs is strongly associated with distinct spectra of embryonal childhood cancers and AML with BRCA2-PGVs, and also early epithelial cancers with BRCA1 PGVs, germline variants in the BRCA1/2 genes have also been identified in non-FA childhood malignancies, and thereby implying the possibility of a role of BRCA PGVs also for non-syndromic cancer predisposition in children. We provide a concise review of aspects of the clinical and genetic features of BRCA1/2-associated FA with a focus on associated malignancies, and review novel aspects of the role of germline BRCA2 and BRCA1 PGVs occurring in non-FA childhood cancer and discuss aspects of clinical and biological implications.

BRCA1 and BRCA2 whole cDNA analysis in unsolved hereditary breast/ovarian cancer patients

Germline pathogenic variants in BRCA1 and BRCA2 genes (BRCA1/2) explain an important fraction of hereditary breast/ovarian cancer (HBOC) cases. Genetic testing generally involves examining coding regions and exon/intron boundaries, thus the frequency of deleterious variants in non-coding regions is unknown. Here we analysed BRCA1/2 whole cDNA in a large cohort of 320 unsolved high-risk HBOC cases in order to identify potential splicing alterations explained by variants in BRCA1/2 deep intronic regions. Whole RNA splicing profiles were analysed by RT-PCR using Sanger sequencing or high-resolution electrophoresis in a QIAxcel instrument. Known predominant BRCA1/2 alternative splicing events were detected, together with two novel events BRCA1 ▼21 and BRCA2 Δ18q_27p. BRCA2 exon 3 skipping was detected in one patient (male) affected with breast cancer, caused by a known Portuguese founder mutation (c.156_157insAluYa5). An altered BRCA2 splicing pattern was detected in three patients, consisting in the up-regulation of ▼20A, Δ22 and ▼20A+Δ22 transcripts. In silico analysis and semi-quantitative data identified the polymorphism BRCA2 c.8755-66T>C as a potential modifier of Δ22 levels. Our findings suggest that mRNA alterations in BRCA1/2 caused by deep intronic variants are rare in Spanish population. However, RNA analysis complements DNA-based strategies allowing the identification of alterations that could go undetected by conventional testing.

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.

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

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

Implementing NGS-based BRCA tumour tissue testing in FFPE ovarian carcinoma specimens: hints from a real-life experience within the framework of expert recommendations

AIMS

Next Generation Sequencing (NGS)-based BRCA tumour tissue testing poses several challenges. As a first step of its implementation within a regional health service network, an in-house validation study was compared with published recommendations.

METHODS

Epithelial ovarian cancer (EOC) formalin-fixed paraffin-embedded specimens stored in the archives of the eight regional pathology units were selected from a consecutive series of patients with known BRCA germline status. Two expert pathologists evaluated tumour cell content for manual macrodissection. DNA extraction, library preparation and NGS analyses were performed blinded to the germinal status. Parameters used in the study were confronted with guidelines for the validation of NGS-based oncology panels and for BRCA tumour tissue testing.

RESULTS

NGS analyses were successful in 66 of 67 EOC specimens, with good quality metrics and high reproducibility among different runs. In all, 19 BRCA pathogenic variants were identified: 12 were germline and 7 were somatic. A 100% concordance with blood tests was detected for germline variants. A BRCA1 variant showed a controversial classification. In different areas of two early stage EOCs showing somatic variants, intratumour heterogeneity not relevant for test results (variant allele frequency >5%) was observed. Compared with expert recommendations, main limitations of the study were absence of controls with known somatic BRCA status and exclusion from the validation of BRCA copy number variations (CNV).

CONCLUSIONS

A close collaboration between pathology and genetics units provides advantages in the implementation of BRCA tumour tissue testing. The development of tools for designing and interpreting complex testing in-house validation could improve process quality.

Lessons learned from understanding chemotherapy resistance in epithelial tubo-ovarian carcinoma from BRCA1and BRCA2mutation carriers

BRCA1 and BRCA2 are multi-functional proteins and key factors for maintaining genomic stability through their roles in DNA double strand break repair by homologous recombination, rescuing stalled or damaged DNA replication forks, and regulation of cell cycle DNA damage checkpoints. Impairment of any of these critical roles results in genomic instability, a phenotypic hallmark of many cancers including breast and epithelial ovarian carcinomas (EOC). Damaging, usually loss of function germline and somatic variants in BRCA1 and BRCA2, are important drivers of the development, progression, and management of high-grade serous tubo-ovarian carcinoma (HGSOC). However, mutations in these genes render patients particularly sensitive to platinum-based chemotherapy, and to the more innovative targeted therapies with poly-(ADP-ribose) polymerase inhibitors (PARPis) that are targeted to BRCA1/BRCA2 mutation carriers. Here, we reviewed the literature on the responsiveness of BRCA1/2-associated HGSOC to platinum-based chemotherapy and PARPis, and propose mechanisms underlying the frequent development of resistance to these therapeutic agents.

Calibration of Pathogenicity Due to Variant-Induced Leaky Splicing Defects by Using BRCA2 Exon 3 as a Model System

BRCA2 is a clinically actionable gene implicated in breast and ovarian cancer predisposition that has become a high priority target for improving the classification of variants of unknown significance (VUS). Among all BRCA2 VUS, those causing partial/leaky splicing defects are the most challenging to classify because the minimal level of full-length (FL) transcripts required for normal function remains to be established. Here, we explored BRCA2 exon 3 (BRCA2e3) as a model for calibrating variant-induced spliceogenicity and estimating thresholds for BRCA2 haploinsufficiency. In silico predictions, minigene splicing assays, patients' RNA analyses, a mouse embryonic stem cell (mESC) complementation assay and retrieval of patient-related information were combined to determine the minimal requirement of FL BRCA2 transcripts. Of 100 BRCA2e3 variants tested in the minigene assay, 64 were found to be spliceogenic, causing mild to severe RNA defects. Splicing defects were also confirmed in patients' RNA when available. Analysis of a neutral leaky variant (c.231T>G) showed that a reduction of approximately 60% of FL BRCA2 transcripts from a mutant allele does not cause any increase in cancer risk. Moreover, data obtained from mESCs suggest that variants causing a decline in FL BRCA2 with approximately 30% of wild-type are not pathogenic, given that mESCs are fully viable and resistant to DNA-damaging agents in those conditions. In contrast, mESCs producing lower relative amounts of FL BRCA2 exhibited either null or hypomorphic phenotypes. Overall, our findings are likely to have broader implications on the interpretation of BRCA2 variants affecting the splicing pattern of other essential exons. SIGNIFICANCE: These findings demonstrate that BRCA2 tumor suppressor function tolerates substantial reduction in full-length transcripts, helping to determine the pathogenicity of BRCA2 leaky splicing variants, some of which may not increase cancer risk.

HSF2BP negatively regulates homologous recombination in DNA interstrand crosslink repair

The tumor suppressor BRCA2 is essential for homologous recombination (HR), replication fork stability and DNA interstrand crosslink (ICL) repair in vertebrates. We show that ectopic production of HSF2BP, a BRCA2-interacting protein required for meiotic HR during mouse spermatogenesis, in non-germline human cells acutely sensitize them to ICL-inducing agents (mitomycin C and cisplatin) and PARP inhibitors, resulting in a phenotype characteristic of cells from Fanconi anemia (FA) patients. We biochemically recapitulate the suppression of ICL repair and establish that excess HSF2BP compromises HR by triggering the removal of BRCA2 from the ICL site and thereby preventing the loading of RAD51. This establishes ectopic expression of a wild-type meiotic protein in the absence of any other protein-coding mutations as a new mechanism that can lead to an FA-like cellular phenotype. Naturally occurring elevated production of HSF2BP in tumors may be a source of cancer-promoting genomic instability and also a targetable vulnerability.

Criteria of the German Consortium for Hereditary Breast and Ovarian Cancer for the Classification of Germline Sequence Variants in Risk Genes for Hereditary Breast and Ovarian Cancer

More than ten years ago, the German Consortium for Hereditary Breast and Ovarian Cancer (GC-HBOC) set up a panel of experts (VUS Task Force) which was tasked with reviewing the classifications of genetic variants reported by individual centres of the GC-HBOC to the central database in Leipzig and reclassifying them, where necessary, based on the most recent data. When it evaluates variants, the VUS Task Force must arrive at a consensus. The resulting classifications are recorded in a central database where they serve as a basis for ensuring the consistent evaluation of previously known and newly identified variants in the different centres of the GC-HBOC. The standardised VUS evaluation by the VUS Task Force is a key element of the recall system which has also been set up by the GC-HBOC. The system will be used to pass on information to families monitored and managed by GC-HBOC centres in the event that previously classified variants are reclassified based on new information. The evaluation algorithm of the VUS Task Force was compiled using internationally established assessment methods (IARC, ACMG, ENIGMA) and is presented here together with the underlying evaluation criteria used to arrive at the classification decision using a flow chart. In addition, the characteristics and special features of specific individual risk genes associated with breast and/or ovarian cancer are discussed in separate subsections. The URLs of relevant databases have also been included together with extensive literature references to provide additional information and cover the scope and dynamism of the current state of knowledge on the evaluation of genetic variants. In future, if criteria are updated based on new information, the update will be published on the website of the GC-HBOC ( https://www.konsortium-familiaerer-brustkrebs.de/ ).

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.

Germline investigation in male breast cancer of DNA repair genes by next-generation sequencing

PURPOSE

In order to better define the breast cancer (BC) genetic risk factors in men, a germline investigation was carried out on 81 Male BC cases by screening the 24 genes involved in BC predisposition, genome stability maintenance and DNA repair mechanisms by next-generation sequencing.

METHODS

Germline DNAs were tested in a custom multi-gene panel focused on all coding exons and exon-intron boundaries of 24 selected genes using two amplicon-based assays on PGM-Ion Torrent (ThermoFisher Scientific) and MiSeq (Illumina) platforms. All variants were recorded and classified by using a custom pipeline.

RESULTS

Clinical pathological data and the family history of 81 Male BC cases were gathered and analysed, revealing the average age of onset to be 61.3 years old and that in 35 cases there was a family history of BC. Our genetic screening allowed us to identify a germline mutation in 22 patients (23%) in 4 genes: BRCA2, BRIP1, MUTYH and PMS2. Moreover, 12 variants of unknown clinical significance (VUS) in 9 genes (BARD1, BRCA1, BRIP1, CHEK2, ERCC1, NBN, PALB2, PMS1, RAD50) were predicted as potentially pathogenic by in silico analysis bringing the mutation detection rate up to 40%.

CONCLUSION

As expected, a positive family history is a strong predictor of germline BRCA2 mutations in male BC. Understanding the potential pathogenicity of VUS represents an extremely urgent need for the management of BC risk in Male BC cases and their own families.

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.

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.

Diagnostic mRNA splicing assay for variants in BRCA1 and BRCA2 identified two novel pathogenic splicing aberrations

Background

Pathogenic variants in BRCA1 and BRCA2 cause hereditary breast and ovarian cancer. Screening of these genes has become easily accessible in diagnostic laboratories. Sequencing and copy number analyses are used to detect pathogenic variants, but also lead to identification of variants of unknown clinical significance (VUS). If the effect of a VUS can be clarified, it has direct consequence for the clinical management of the patient and family members. A splicing assay is one of several tools that might help in the classification of VUS. We therefore established mRNA analyses for BRCA1 and BRCA2 in the diagnostic laboratory in 2015. We hereby report the results of mRNA analysis variants in BRCA1 and BRCA2 after three years.

Methods

Variants predicted to alter splicing and variants within the canonical splice sites were selected for splicing analyses. Splicing assays were performed by reverse transcription-PCR of patient RNA. A biallalic expression analysis was carried out whenever possible.

Results

Twenty-five variants in BRCA1 and BRCA2 were analyzed by splicing assays; nine showed altered transcripts and 16 showed normal splicing patterns. The two novel pathogenic variants in BRCA1 c.4484 + 3 A > C and c.5407–10G > A were characterized.

Conclusions

We conclude that mRNA analyses are useful in characterization of variants that may affect splicing. The results can guide classification of variants from unknown clinical significance to pathogenic or benign in a diagnostic laboratory, and thus be of direct clinical importance.

BRCA1-No Matter How You Splice It

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

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

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

[Hereditary breast and ovarian cancer]

Hereditary breast and ovarian carcinomas are frequently caused by germline mutations of the BRCA1 and BRCA2 genes (BRCA1/2 syndromes) and are often less associated with other hereditary syndromes such as Li-Fraumeni and Peutz-Jeghers. The BRCA1/2 proteins have a special role in DNA repair. Therefore, loss of function due to mutation causes an accumulation of mutations in other genes and subsequent tumorigenesis at an early age. BRCA1/2 mutations are irregularly distributed over the length of the genes without hot spots, although special mutations are known. Breast and ovarian cancer occur far more frequently in women with BRCA1/2 germline mutations compared with the general population. Breast cancer occurs increasingly from the age of 30, ovarian cancer in BRCA1 syndrome from the age of 40 and BRCA2 from the age of 50. Suspicion of a BRCA syndrome should be prompted in the case of clustering of breast cancer in 1st degree relatives, in particular at a young age, if breast and ovarian cancer have occurred, and if cases of male breast cancer are known. Breast carcinomas with medullary differentiation seem to predominate in BRCA syndromes, but other carcinoma types may also occur. BRCA germline mutations seem to occur frequently in triple-negative breast carcinomas, whereas an association with ductal carcinoma in situ (DCIS) is rare. Ovarian carcinomas in BRCA syndromes are usually high-grade serous, mucinous carcinomas and borderline tumors are unusual. Pathology plays a special role within the multidisciplinary team in the recognition of patients with hereditary cancer syndromes.

Alternative splicing and ACMG-AMP-2015-based classification of PALB2 genetic variants: an ENIGMA report

BACKGROUND

PALB2 monoallelic loss-of-function germ-line variants confer a breast cancer risk comparable to the average BRCA2 pathogenic variant. Recommendations for risk reduction strategies in carriers are similar. Elaborating robust criteria to identify loss-of-function variants in PALB2-without incurring overprediction-is thus of paramount clinical relevance. Towards this aim, we have performed a comprehensive characterisation of alternative splicing in PALB2, analysing its relevance for the classification of truncating and splice site variants according to the 2015 American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines.

METHODS

Alternative splicing was characterised in RNAs extracted from blood, breast and fimbriae/ovary-related human specimens (n=112). RNAseq, RT-PCR/CE and CloneSeq experiments were performed by five contributing laboratories. Centralised revision/curation was performed to assure high-quality annotations. Additional splicing analyses were performed in PALB2 c.212-1G>A, c.1684+1G>A, c.2748+2T>G, c.3113+5G>A, c.3350+1G>A, c.3350+4A>C and c.3350+5G>A carriers. The impact of the findings on PVS1 status was evaluated for truncating and splice site variant.

RESULTS

We identified 88 naturally occurring alternative splicing events (81 newly described), including 4 in-frame events predicted relevant to evaluate PVS1 status of splice site variants. We did not identify tissue-specific alternate gene transcripts in breast or ovarian-related samples, supporting the clinical relevance of blood-based splicing studies.

CONCLUSIONS

PVS1 is not necessarily warranted for splice site variants targeting four PALB2 acceptor sites (exons 2, 5, 7 and 10). As a result, rare variants at these splice sites cannot be assumed pathogenic/likely pathogenic without further evidences. Our study puts a warning in up to five PALB2 genetic variants that are currently reported as pathogenic/likely pathogenic in ClinVar.

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

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

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

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

Quantifying BRCA1 and BRCA2 mRNA Isoform Expression Levels in Single Cells

BRCA1 and BRCA2 spliceogenic variants are often associated with an elevated risk of breast and ovarian cancers. Analyses of BRCA1 and BRCA2 splicing patterns have traditionally used technologies that sample a population of cells but do not account for the variation that may be present between individual cells. This novel proof of concept study utilises RNA in situ hybridisation to measure the absolute expression of BRCA1 and BRCA2 mRNA splicing events in single lymphoblastoid cells containing known spliceogenic variants (BRCA1c.671-2 A>G or BRCA2c.7988 A>T). We observed a large proportion of cells (>42%) in each sample that did not express mRNA for the targeted gene. Increased levels (average mRNA molecules per cell) of BRCA2 ∆17_18 were observed in the cells containing the known spliceogenic variant BRCA2c.7988 A>T, but cells containing BRCA1c.671-2 A>G were not found to express significantly increased levels of BRCA1 ∆11, as had been shown previously. Instead, we show for each variant carrier sample that a higher proportion of cells expressed the targeted splicing event compared to control cells. These results indicate that BRCA1/2 mRNA is expressed stochastically, suggesting that previously reported results using RT-PCR may have been influenced by the number of cells with BRCA1/2 mRNA expression and may not represent an elevation of constitutive mRNA expression. Detection of mRNA expression in single cells allows for a more comprehensive understanding of how spliceogenic variants influence the expression of mRNA isoforms. However, further research is required to assess the utility of this technology to measure the expression of predicted spliceogenic BRCA1 and BRCA2 variants in a diagnostic setting.

Quantitative Analysis of BRCA1 and BRCA2 Germline Splicing Variants Using a Novel RNA-Massively Parallel Sequencing Assay

Clinical genetic testing for hereditary breast and ovarian cancer (HBOC) is becoming widespread. However, the interpretation of variants of unknown significance (VUS) in HBOC genes, such as the clinically actionable genes BRCA1 and BRCA2, remain a challenge. Among the variants that are frequently classified as VUS are those with unclear effects on splicing. In order to address this issue we developed a high-throughput RNA-massively parallel sequencing assay—CloneSeq—capable to perform quantitative and qualitative analysis of transcripts in cell lines and HBOC patients. This assay is based on cloning of RT-PCR products followed by massive parallel sequencing of the cloned transcripts. To validate this assay we compared it to the RNA splicing assays recommended by members of the ENIGMA (Evidence-based Network for the Interpretation of Germline Mutant Alleles) consortium. This comparison was performed using well-characterized lymphoblastoid cell lines (LCLs) generated from carriers of the BRCA1 or BRCA2 germline variants that have been previously described to be associated with splicing defects. CloneSeq was able to replicate the ENIGMA results, in addition to providing quantitative characterization of BRCA1 and BRCA2 germline splicing alterations in a high-throughput fashion. Furthermore, CloneSeq was used to analyze blood samples obtained from carriers of BRCA1 or BRCA2 germline sequence variants, including the novel uncharacterized alteration BRCA1 c.5152+5G>T, which was identified in a HBOC family. CloneSeq provided a high-resolution picture of all the transcripts induced by BRCA1 c.5152+5G>T, indicating it results in significant levels of exon skipping. This analysis proved to be important for the classification of BRCA1 c.5152+5G>T as a clinically actionable likely pathogenic variant. Reclassifications such as these are fundamental in order to offer preventive measures, targeted treatment, and pre-symptomatic screening to the correct individuals.

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.

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

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

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

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

The BRCA2 c.68-7T > A variant is not pathogenic: A model for clinical calibration of spliceogenicity

Although the spliceogenic nature of the BRCA2 c.68-7T > A variant has been demonstrated, its association with cancer risk remains controversial. In this study, we accurately quantified by real-time PCR and digital PCR (dPCR), the BRCA2 isoforms retaining or missing exon 3. In addition, the combined odds ratio for causality of the variant was estimated using genetic and clinical data, and its associated cancer risk was estimated by case-control analysis in 83,636 individuals. Co-occurrence in trans with pathogenic BRCA2 variants was assessed in 5,382 families. Exon 3 exclusion rate was 4.5-fold higher in variant carriers (13%) than controls (3%), indicating an exclusion rate for the c.68-7T > A allele of approximately 20%. The posterior probability of pathogenicity was 7.44 × 10-115 . There was neither evidence for increased risk of breast cancer (OR 1.03; 95% CI 0.86-1.24) nor for a deleterious effect of the variant when co-occurring with pathogenic variants. Our data provide for the first time robust evidence of the nonpathogenicity of the BRCA2 c.68-7T > A. Genetic and quantitative transcript analyses together inform the threshold for the ratio between functional and altered BRCA2 isoforms compatible with normal cell function. These findings might be exploited to assess the relevance for cancer risk of other BRCA2 spliceogenic variants.

Thorough in silico and in vitro cDNA analysis of 21 putative BRCA1 and BRCA2 splice variants and a complex tandem duplication in BRCA2 allowing the identification of activated cryptic splice donor sites in BRCA2 exon 11

For 21 putative BRCA1 and BRCA2 splice site variants, the concordance between mRNA analysis and predictions by in silico programs was evaluated. Aberrant splicing was confirmed for 12 alterations. In silico prediction tools were helpful to determine for which variants cDNA analysis is warranted, however, predictions for variants in the Cartegni consensus region but outside the canonical sites, were less reliable. Learning algorithms like Adaboost and Random Forest outperformed the classical tools. Further validations are warranted prior to implementation of these novel tools in clinical settings. Additionally, we report here for the first time activated cryptic donor sites in the large exon 11 of BRCA2 by evaluating the effect at the cDNA level of a novel tandem duplication (5' breakpoint in intron 4; 3' breakpoint in exon 11) and of a variant disrupting the splice donor site of exon 11 (c.6841+1G > C). Additional sites were predicted, but not activated. These sites warrant further research to increase our knowledge on cis and trans acting factors involved in the conservation of correct transcription of this large exon. This may contribute to adequate design of ASOs (antisense oligonucleotides), an emerging therapy to render cancer cells sensitive to PARP inhibitor and platinum therapies.

Characterization of BRCA1 and BRCA2 variants found in a Norwegian breast or ovarian cancer cohort

Germline mutations in BRCA1 and BRCA2 cause hereditary breast and ovarian cancer. Molecular screening of these two genes in patients with a family history of breast or ovarian cancer has revealed pathogenic variants as well as genetic variants of unknown significance (VUS). These VUS may cause a challenge in the genetic counseling process regarding clinical management of the patient and the family. We investigated 32 variants previously detected in 33 samples from patients with a family history of breast or ovarian cancer. cDNA was analyzed for alternative transcripts and selected missense variants located in the BRCT domains of BRCA1 were assessed for their trans-activation ability. Although an extensive cDNA analysis was done, only three of the 32 variants appeared to affect the splice-process (BRCA1 c.213-5T>A, BRCA1 c.5434C>G and BRCA2 c.68-7T>A). In addition, two variants located in the BRCT domains of BRCA1 (c.5075A>C p.Asp1692Ala and c.5513T>G p.Val1838Gly) were shown to abolish the BRCT domain trans-activation ability, whereas BRCA1 c.5125G>A p.Gly1709Arg exhibited equal trans-activation capability as the WT domain. These functional studies may offer further insights into the pathogenicity of certain identified variants; however, this assay is only applicable for a subset of missense variants.

Detecting splicing patterns in genes involved in hereditary breast and ovarian cancer

Interpretation of variants of unknown significance (VUS) is a major challenge for laboratories performing molecular diagnosis of hereditary breast and ovarian cancer (HBOC), especially considering that many genes are now known to be involved in this syndrome. One important way these VUS can have a functional impact is through their effects on RNA splicing. Here we present a custom RNA-Seq assay plus bioinformatics and biostatistics pipeline to analyse specifically alternative and abnormal splicing junctions in 11 targeted HBOC genes. Our pipeline identified 14 new alternative splices in BRCA1 and BRCA2 in addition to detecting the majority of known alternative spliced transcripts therein. We provide here the first global splicing pattern analysis for the other nine genes, which will enable a comprehensive interpretation of splicing defects caused by VUS in HBOC. Previously known splicing alterations were consistently detected, occasionally with a more complex splicing pattern than expected. We also found that splicing in the 11 genes is similar in blood and breast tissue, supporting the utility and simplicity of blood splicing assays. Our pipeline is ready to be integrated into standard molecular diagnosis for HBOC, but it could equally be adapted for an integrative analysis of any multigene disorder.

Acquired cross-linker resistance associated with a novel spliced BRCA2 protein variant for molecular phenotyping of BRCA2 disruption

BRCA2 encodes a protein with a fundamental role in homologous recombination that is essential for normal development. Carrier status of mutations in BRCA2 is associated with familial breast and ovarian cancer, while bi-allelic BRCA2 mutations can cause Fanconi anemia (FA), a cancer predisposition syndrome with cellular cross-linker hypersensitivity. Cancers associated with BRCA2 mutations can acquire chemo-resistance on relapse. We modeled acquired cross-linker resistance with an FA-derived BRCA2-mutated acute myeloid leukemia (AML) platform. Associated with acquired cross-linker resistance was the expression of a functional BRCA2 protein variant lacking exon 5 and exon 7 (BRCA2^ΔE5+7), implying a role for BRCA2 splicing for acquired chemo-resistance. Integrated network analysis of transcriptomic and proteomic differences for phenotyping of BRCA2 disruption infers impact on transcription and chromatin remodeling in addition to the DNA damage response. The striking overlap with transcriptional profiles of FA patient hematopoiesis and BRCA mutation associated ovarian cancer helps define and explicate the ‘BRCAness’ profile.

Functional classification of DNA variants by hybrid minigenes: Identification of 30 spliceogenic variants of BRCA2 exons 17 and 18

Mutation screening of the breast cancer genes BRCA1 and BRCA2 identifies a large fraction of variants of uncertain clinical significance (VUS) whose functional and clinical interpretations pose a challenge for genomic medicine. Likewise, an increasing amount of evidence indicates that genetic variants can have deleterious effects on pre-mRNA splicing. Our goal was to investigate the impact on splicing of a set of reported variants of BRCA2 exons 17 and 18 to assess their role in hereditary breast cancer and to identify critical regulatory elements that may constitute hotspots for spliceogenic variants. A splicing reporter minigene with BRCA2 exons 14 to-20 (MGBR2_ex14-20) was constructed in the pSAD vector. Fifty-two candidate variants were selected with splicing prediction programs, introduced in MGBR2_ex14-20 by site-directed mutagenesis and assayed in triplicate in MCF-7 cells. Wild type MGBR2_ex14-20 produced a stable transcript of the expected size (1,806 nucleotides) and structure (V1-[BRCA2_exons_14–20]–V2). Functional mapping by microdeletions revealed essential sequences for exon recognition on the 3’ end of exon 17 (c.7944-7973) and the 5’ end of exon 18 (c.7979-7988, c.7999-8013). Thirty out of the 52 selected variants induced anomalous splicing in minigene assays with >16 different aberrant transcripts, where exon skipping was the most common event. A wide range of splicing motifs were affected including the canonical splice sites (15 variants), novel alternative sites (3 variants), the polypyrimidine tract (3 variants) and enhancers/silencers (9 variants). According to the guidelines of the American College of Medical Genetics and Genomics (ACMG), 20 variants could be classified as pathogenic (c.7806-2A>G, c.7806-1G>A, c.7806-1G>T, c.7806-1_7806-2dup, c.7976+1G>A, c.7977-3_7978del, c.7977-2A>T, c.7977-1G>T, c.7977-1G>C, c.8009C>A, c.8331+1G>T and c.8331+2T>C) or likely pathogenic (c.7806-9T>G, c.7976G>C, c.7976G>A, c.7977-7C>G, c.7985C>G, c.8023A>G, c.8035G>T and c.8331G>A), accounting for 30.8% of all pathogenic/likely pathogenic variants of exons 17–18 at the BRCA Share database. The remaining 8 variants (c.7975A>G, c.7977-6T>G, c.7988A>T, c.7992T>A, c.8007A>G, c.8009C>T, c.8009C>G, and c.8072C>T) induced partial splicing anomalies with important ratios of the full-length transcript (≥70%), so that they remained classified as VUS. Aberrant splicing is therefore especially prevalent in BRCA2 exons 17 and 18 due to the presence of active ESEs involved in exon recognition. Splicing functional assays with minigenes are a valuable strategy for the initial characterization of the splicing outcomes and the subsequent clinical interpretation of variants of any disease-gene, although these results should be checked, whenever possible, against patient RNA.

A significant proportion of disease-causing mutations of inherited disorders impair splicing. Massive sequencing projects of genetic diseases generate thousands of sequence variations that require functional and clinical interpretations. We have shown that splicing reporter minigenes of the breast cancer genes BRCA1 and BRCA2 are useful tools to functionally test DNA variants. In this work, we have constructed a 7-exon BRCA2 minigene (exons 14 to 20) where we mapped critical splicing regulatory sequences and tested 52 selected variants of exons 17 and 18 detected in breast cancer patients. We finely located three DNA segments on both exons that presumably contain splicing enhancer sequences. We observed that a total of 30 variants of any type disrupted the splicing patterns and, given the severity of their outcomes, we classified 20 of them as pathogenic or likely pathogenic. We also showed that a wide range of splicing elements were affected including canonical and novel 5’ and 3’ splice sites, the polypyrimidine tract and enhancer and silencer sequences. We concluded that splicing aberrations are frequent in Hereditary Breast and Ovarian Cancer and that minigenes are valuable tools to functionally classify DNA variants of any human disease gene under the splicing viewpoint.

Hereditary breast and ovarian cancer: new genes in confined pathways

Genetic abnormalities in the DNA repair genes BRCA1 and BRCA2 predispose to hereditary breast and ovarian cancer (HBOC). However, only approximately 25% of cases of HBOC can be ascribed to BRCA1 and BRCA2 mutations. Recently, exome sequencing has uncovered substantial locus heterogeneity among affected families without BRCA1 or BRCA2 mutations. The new pathogenic variants are rare, posing challenges to estimation of risk attribution through patient cohorts. In this Review article, we examine HBOC genes, focusing on their role in genome maintenance, the possibilities for functional testing of putative causal variants and the clinical application of new HBOC genes in cancer risk management and treatment decision-making.