BRCA1‐ and BRCA2‐specific in silico tools for variant interpretation in the CAGI 5 ENIGMA challenge

BRCA1 and BRCA2 (BRCA1/2) germline variants disrupting the DNA protective role of these genes increase the risk of hereditary breast and ovarian cancers. Correct identification of these variants then becomes clinically relevant, because it may increase the survival rates of the carriers. Unfortunately, we are still unable to systematically predict the impact of BRCA1/2 variants. In this article, we present a family of in silico predictors that address this problem, using a gene‐specific approach. For each protein, we have developed two tools, aimed at predicting the impact of a variant at two different levels: Functional and clinical. Testing their performance in different datasets shows that specific information compensates the small number of predictive features and the reduced training sets employed to develop our models. When applied to the variants of the BRCA1/2 (ENIGMA) challenge in the fifth Critical Assessment of Genome Interpretation (CAGI 5) we find that these methods, particularly those predicting the functional impact of variants, have a good performance, identifying the large compositional bias towards neutral variants in the CAGI sample. This performance is further improved when incorporating to our prediction protocol estimates of the impact on splicing of the target variant.     

Variant effect on splicing regulatory elements,branchpoint usage,and pseudoexonization: Strategies to enhance bioinformatic prediction using hereditary cancer genes as exemplars

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre‐existing large‐scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.     

Effect of <Emphasis Type="Italic">BRCA2</Emphasis> sequence variants predicted to disrupt exonic splice enhancers on <Emphasis Type="Italic">BRCA2</Emphasis>transcripts

Background  

Genetic screening of breast cancer patients and their families have identified a number of variants of unknown clinical significance in the breast cancer susceptibility genes, BRCA1 and BRCA2. Evaluation of such unclassified variants may be assisted by web-based bioinformatic prediction tools, although accurate prediction of aberrant splicing by unclassified variants affecting exonic splice enhancers (ESEs) remains a challenge.     

Functional redundancy of exon 12 of BRCA2 revealed by a comprehensive analysis of the c.6853A>G (p.I2285V) variant

Variants of unknown significance (VUS) in BRCA1 and BRCA2 are common, and present significant challenges for genetic counseling. We observed that BRCA2: c.6853A>G (p.I2285V) (Breast Cancer Information Core [BIC] name: 7081A>G; http://research.nhgri.nih.gov/bic/ ) co-occurs in trans with the founder mutation c.5946delT (p.S1982RfsX22) (BIC name: 6174delT), supporting the published classification of p.I2285V as a neutral variant. However, we also noted that when compared with wild-type BRCA2, p.I2285V resulted in increased exclusion of exon 12. Functional assay using allelic complementation in Brca2-null mouse embryonic stem cells revealed that p.I2285V, an allele with exon 12 deleted and wild-type BRCA2 were all phenotypically indistinguishable, as measured by sensitivity to DNA-damaging agents, effect on irradiation-induced Rad51 foci formation, homologous recombination, and overall genomic integrity. An allele frequency study showed the p.I2285V variant was identified in 15 out of 722 (2.1%) Ashkenazi Jewish cases and 10 out of 475 (2.1%) ethnically-matched controls (odds ratio, 0.99; 95% confidence interval: 0.44–2.21; P=0.97). Thus the p.I2285V variant is not associated with an increased risk for breast cancer. Taken together, our clinical and functional studies strongly suggest that exon 12 is functionally redundant and therefore missense variants in this exon are likely to be neutral. Such comprehensive functional studies will be important adjuncts to genetic studies of variants. Hum Mutat 30:1–8, 2009. Published 2009 Wiley-Liss, Inc.     

Genotype-first approach to the detection of hereditary breast and ovarian cancer risk,and effects of risk disclosure to biobank participants

Genotype-first approach allows to systematically identify carriers of pathogenic variants in BRCA1/2 genes conferring a high risk of familial breast and ovarian cancer. Participants of the Estonian biobank have expressed support for the disclosure of clinically significant findings. With an Estonian biobank cohort, we applied a genotype-first approach, contacted carriers, and offered return of results with genetic counseling. We evaluated participants’ responses to and the clinical utility of the reporting of actionable genetic findings. Twenty-two of 40 contacted carriers of 17 pathogenic BRCA1/2 variants responded and chose to receive results. Eight of these 22 participants qualified for high-risk assessment based on National Comprehensive Cancer Network criteria. Twenty of 21 counseled participants appreciated being contacted. Relatives of 10 participants underwent cascade screening. Five of 16 eligible female BRCA1/2 variant carriers chose to undergo risk-reducing surgery, and 10 adhered to surveillance recommendations over the 30-month follow-up period. We recommend the return of results to population-based biobank participants; this approach could be viewed as a model for population-wide genetic testing. The genotype-first approach permits the identification of individuals at high risk who would not be identified by application of an approach based on personal and family histories only.Subject terms: Genetic counselling, Population screening, Cancer genomics     

Loss‐of‐function variants of the filaggrin gene are associated with clinical reactivity to foods

The aim of this study was to assess the genetic association of Filaggrin loss‐of‐function (FLG LOF) genetic variants with food allergy, and to investigate the added value of this test in diagnosing food allergy. Clinical reactivity to foods was diagnosed by the gold standard, the double‐blind, placebo‐controlled food challenge. Of 155 children, 33 (21.3%) children had at least one FLG LOF variant, and of these, 29 (87.9%) were clinically reactive to at least one food, compared to 73 of 122 children (59.8%) carrying wild‐type alleles. The odds ratio for having at least one FLG LOF variant and clinical reactivity to at least one food was 4.9 (CI = 1.6–14.7, P = 0.005), corresponding to a relative risk of 1.5, compared to carriers of wild‐type alleles. Prediction of food allergy improved when FLG LOF variants were included in the model. Therefore, genetic markers may be useful as an addition to clinical assessment in the diagnosis of food allergy.     

Functional Analysis of a Large set of BRCA2 exon 7 Variants Highlights the Predictive Value of Hexamer Scores in Detecting Alterations of Exonic Splicing Regulatory Elements

Exonic variants can alter pre‐mRNA splicing either by changing splice sites or by modifying splicing regulatory elements. Often these effects are difficult to predict and are only detected by performing RNA analyses. Here, we analyzed, in a minigene assay, 26 variants identified in the exon 7 of BRCA2, a cancer predisposition gene. Our results revealed eight new exon skipping mutations in this exon: one directly altering the 5′ splice site and seven affecting potential regulatory elements. This brings the number of splicing regulatory mutations detected in BRCA2 exon 7 to a total of 11, a remarkably high number considering the total number of variants reported in this exon (n = 36), all tested in our minigene assay. We then exploited this large set of splicing data to test the predictive value of splicing regulator hexamers’ scores recently established by Ke et al. ( 2011 ). Comparisons of hexamer‐based predictions with our experimental data revealed high sensitivity in detecting variants that increased exon skipping, an important feature for prescreening variants before RNA analysis. In conclusion, hexamer scores represent a promising tool for predicting the biological consequences of exonic variants and may have important applications for the interpretation of variants detected by high‐throughput sequencing.     

Functional characterization of BRCA1 gene variants by mini-gene splicing assay

Mutational screening of the breast cancer susceptibility gene BRCA1 leads to the identification of numerous pathogenic variants such as frameshift and nonsense variants, as well as large genomic rearrangements. The screening moreover identifies a large number of variants, for example, missense, silent, and intron variants, which are classified as variants of unknown clinical significance owing to the lack of causal evidence. Variants of unknown clinical significance can potentially have an impact on splicing and therefore functional examinations are warranted to classify whether these variants are pathogenic or benign. Here we validate a mini-gene splicing assay by comparing the results of 24 variants with previously published data from RT-PCR analysis on RNA from blood samples/lymphoblastoid cell lines. The analysis showed an overall concordance of 100%. In addition, we investigated 13 BRCA1 variants of unknown clinical significance or putative variants affecting splicing by in silico analysis and mini-gene splicing assay. Both the in silico analysis and mini-gene splicing assay classified six BRCA1 variants as pathogenic (c.80+1G>A, c.132C>T (p.=), c.213−1G>A, c.670+1delG, c.4185+1G>A, and c.5075−1G>C), whereas six BRCA1 variants were classified as neutral (c.-19-22_-19-21dupAT, c.302−15C>G, c.547+14delG, c.4676−20A>G, c.4987−21G>T, and c.5278−14C>G) and one BRCA1 variant remained unclassified (c.670+16G>A). In conclusion, our study emphasizes that in silico analysis and mini-gene splicing assays are important for the classification of variants, especially if no RNA is available from the patient. This knowledge is crucial for proper genetic counseling of patients and their family members.     

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. Hum Mutat 33:1526–1537, 2012. © 2012 Wiley Periodicals, Inc.     

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

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

Somatic genetic alterations in BRCA2‐associated and sporadic male breast cancer

The genetic changes underlying the development and progression of male breast cancer are poorly understood. Germline BRCA2 mutations account for a significant part of male breast cancer, but the majority of patients lack a known inherited predisposition. We recently demonstrated that the progression of breast cancer in female carriers of a germline BRCA1 or BRCA2 mutation follows specific genetic pathways, distinct from each other and from sporadic breast cancer. In the present study, we performed a genome‐wide survey by comparative genomic hybridization (CGH) of somatic genetic aberrations in 26 male breast cancers, including five tumors from BRCA2 mutation carriers. BRCA2 tumors exhibited a significantly higher number of chromosomal aberrations than sporadic tumors. The most common alterations in sporadic male breast cancer were +1q (38%), +8q (33%), +17q (33%), –13q (29%), and –8p (24%). In tumors from BRCA2 mutation carriers, the five most common genetic changes were +8q (100%), +20q (100%), +17q (80%), –13q (80%), and –6q (60%). The CGH results in these two groups of male breast cancers are almost identical to those identified in the corresponding sporadic and BRCA2‐associated female breast cancers. The results suggest that despite substantial hormonal differences between females and males, similar genetic changes are selected for during tumor progression. Furthermore, the presence of a highly penetrant germline BRCA2 mutation apparently leads to a characteristic somatic tumor progression pathway, again shared between affected male and female mutation carriers. Genes Chromosomes Cancer 24:56–61, 1999. © 1999 Wiley‐Liss, Inc.