Comparative in vitro and in silico analyses of variants in splicing regions of BRCA1 and BRCA2 genes and characterization of novel pathogenic mutations

Determination of genetic predisposition to early breast cancer in women of Kazakh ethnicity

Breast cancer (BC) is the most common type of cancer among women in Kazakhstan. To date, little data are available on the spectrum of genetic variation in Kazakh women with BC. We aimed to identify population-specific genetic markers associated with the risk of developing early-onset BC and test their association with clinical and prognostic factors. The study included 224 Kazakh women diagnosed with BC (≤40 age). Entire coding regions (>1700 exons) and the flanking noncoding regions of 94 cancer-associated genes were sequenced from blood DNA using MiSeq platform. We identified 38 unique pathogenic variants (PVs) in 13 different cancer-predisposing genes among 57 patients (25.4%), of which 6 variants were novel. In total, 12 of the 38 distinct PVs were detected recurrently, including BRCA1 c.5266dup, c.5278-2del, and c.2T>C, and BRCA2 c.9409dup and c.9253del that may be founder in this population. BRCA1 carriers were significantly more likely to develop triple-negative BC (OR = 6.61, 95% CI 2.44-17.91, p = 0.0002) and have family history of BC (OR = 3.17, 95% CI 1.14-8.76, p = 0.03) compared to non-carriers. This study allowed the identification of PVs specific to early-onset BC, which may be used as a foundation to develop regional expertise and diagnostic tools for early detection of BC in young Kazakh women.

Unexpected finding of a rare pathogenic germline BRCA1 variant in an intrahepatic cholangiocarcinoma using the Oncomine Focus DNA assay: clinical and diagnostic implications

BACKGROUND

Cholangiocarcinoma (CCA) is a malignant tumor arising from the epithelial cells of the bile ducts and is the second most common liver cancer after hepatocellular carcinoma. Recently, our Institution launched a Comprehensive Genomic Profiling (CGP) program (named FPG500 program), set up to provide a complete molecular characterization through the TruSight Oncology 500 High Throughput (TSO500HT) solution and samples that do not reach pre-set sample quantity and/or quality thresholds required for TSO500HT, are addressed to Oncomine Focus DNA Assay (OFA) and the Archer's FusionPlex Lung Panel (AFL).

METHODS AND RESULTS

Here we report the case of a patient with iCCA enrolled in the FPG500 program and screened by the orthogonal workflow (OFA/AFL). Although BRCA1 is not among the genes declared in the OFA panel, we unexpectedly detected a pathogenic variant in this gene (c.5278-2del, rs878853285).

CONCLUSIONS

This case highlights the diagnostic capabilities of CGP, now widely used in both clinical practice and academic setting. The incidental involvement of BRCA1 focuses attention on the role of BRCA genes in biliary tract cancers. Finally, as an orthogonal test confirmed the germline origin of BRCA1 c.5278-2del variant, the germline implications of CGP need to be considered.

BRCA1 and BRCA2 mutations and their clinical relevance in selected women diagnosed with triple-negative breast cancer in Kenya: a descriptive cross-sectional study

Introduction

triple-negative breast cancer (TNBC) is a heterogeneous breast cancer type with a poor prognosis. About 25% of TNBC patients carry breast cancer susceptibility genes 1 and 2 (BRCA1 and BRCA2) mutations. Screening for BRCA mutations would facilitate early detection and initiation of personalized therapy, thus improving prognosis. However, this has not been explored in our population. We aimed at identifying BRCA1 and BRCA2 gene mutations and their clinical relevance among selected women with TNBC in Kenya.

Methods

six participants enrolled in a larger descriptive cross-sectional study who met the inclusion criteria were selected. Structured questionnaires were used to obtain qualitative data. Deoxyribonucleic acid (DNA) was extracted from saliva. Whole exome sequencing of BRCA1 and BRCA2 genes using a next-generation sequencer was done.

Results

overall, 83.3% of BRCA1 and BRCA2 gene mutations with clinical relevance were detected. Most of the variants (63%) were found in BRCA1 whereas 37% were found in BRCA2. Pathogenic mutations in BRCA1 gene included c.5513T>A, c.5291T>C, c.5297T>G, c.110C>A, c.5212G>C, c.122A>C, c.5117G>A, c.5095C>T, c.5054C>T, c.5053A>G, c.115T>A, c.5143A>G, and c.130T>G. Those in BRCA2 gene were c.7878G>A, c.9154C>T, c.8243G>A, c.7976G>A, c.8165C>G, c.8167G>C, and c.8168A>T. One variant (c.5352delG: p. Leu1785Terfs) not matching any in the BRCA Exchange and ClinVar databases was detected.

Conclusion

our study revealed BRCA mutations that could be common among our population. Further, it has shown that BRCA1 and BRCA2 genetic mutations identified are of clinical relevance and there is a need to screen for these mutations in breast cancer patients to understand their implication in patient management outcomes.

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.

Analyzing the effects of BRCA1/2 variants on mRNA splicing by minigene assay

As BRCA1/2 gene sequencing become more extensive, a large number VUS (variants of uncertain significance) emerge rapidly. Verifying the splicing effect is an effective means for VUS reclassification. The Minigene Assay platform was established and its reliability was verified in this article. 47 BRCA1 or BRCA2 variants were selected and performed to validate their effect on mRNA splicing. The results showed that, a total of 16 variants were experimentally proved to have effects on mRNA splicing, among which 14 variants were shown to cause truncated proteins by Sanger sequencing. While the other two variants, BRCA2 c.7976 + 3 A > G and BRCA1 c.5152 + 3_5152 + 4insT was analyzed to cause 57 bp and 26 bp base in-frame deletion, respectively. The remaining 31 variants were not shown to cause mRNA splicing abnormity, including several sites at the edge of exons, which were predicted to affect splicing of mRNA by multiple bioinformatic software. Based on our experimental results, 37 variants were reclassified by ACMG rules. Our study showed that experimental splicing analysis was effectual for variants classification, and multiple functional assay or clinical data were also necessary for comprehensive judgment of variants.

Empirical prediction of variant-activated cryptic splice donors using population-based RNA-Seq data

Predicting which cryptic-donors may be activated by a splicing variant in patient DNA is notoriously difficult. Through analysis of 5145 cryptic-donors (versus 86,963 decoy-donors not used; any GT or GC), we define an empirical method predicting cryptic-donor activation with 87% sensitivity and 95% specificity. Strength (according to four algorithms) and proximity to the annotated-donor appear important determinants of cryptic-donor activation. However, other factors such as splicing regulatory elements, which are difficult to identify, play an important role and are likely responsible for current prediction inaccuracies. We find that the most frequently recurring natural mis-splicing events at each exon-intron junction, summarised over 40,233 RNA-sequencing samples (40K-RNA), predict with accuracy which cryptic-donor will be activated in rare disease. 40K-RNA provides an accurate, evidence-based method to predict variant-activated cryptic-donors in genetic disorders, assisting pathology consideration of possible consequences of a variant for the encoded protein and RNA diagnostic testing strategies.

Genetic variants affecting the consensus splicing motifs can alter binding of spliceosomal components and induce mis-splicing. Here, the authors develop a method, showing that ranking the most common recurring mis-splicing events in public RNA-Seq data can predict the activation of cryptic-donors.

Consistency of variant interpretations among bioinformaticians and clinical geneticists in hereditary cancer panels

Next-generation sequencing (NGS) is used increasingly in hereditary cancer patients' (HCP) management. While enabling evaluation of multiple genes simultaneously, the technology brings to light the dilemma of variant interpretation. Here, we aimed to reveal the underlying reasons for the discrepancy in the evidence titles used during variant classification according to ACMG guidelines by two different bioinformatic specialists (BIs) and two different clinical geneticists (CGs). We evaluated final reports of 1920 cancer patients and 189 different variants from 285 HCP were enrolled to the study. A total of 173 of these variants were classified as pathogenic (n = 132) and likely pathogenic (n = 41) by the BI and an additional 16 variants, that were classified as VUS by at least one interpreter and their classification would change the clinical management, were compared for their evidence titles between different specialists. The attributed evidence titles and the final classification of the variants among BIs and CGs were compared. The discrepancy between P/LP final reports was 22.5%. The discordance between CGs was 30% whereas the discordance between two BIs was almost 75%. The use of PVS1, PS3, PP3, PP5, PM1, PM2, BP1, BP4 criteria markedly varied from one expert to another. This difference was particularly noticeable in PP3, PP5, and PM1 evidence and mostly in the variants affecting splice sites like BRCA1(NM_007294.4) c.4096 + 1 G > A and CHEK2(NM_007194.4) c.592 + 3 A > T. With recent advancements in precision medicine, the importance of variant interpretations is emerging. Our study shows that variant interpretation is subjective process that is in need of concrete definitions for accurate and standard interpretation.

Background splicing as a predictor of aberrant splicing in genetic disease

Mutations of splice sites, auxiliary splicing elements and the splicing machinery cause a wide range of genetic disease. Here we report that many of the complex effects of splicing mutations can be predicted from background splicing information, with emphasis on BRCA1, BRCA2 and DMD. Background splicing arises from very low level splicing between rarely used background splice sites and from low-level exon skipping between intron splice sites. We show how this information can be downloaded from the Snaptron database of spliced RNA, which we then compared with databases of human splice site mutations. We report that inactivating mutations of intron splice sites typically caused the non-mutated partner splice site to splice to a known background splice site in over 90% of cases and to the strongest background splice site in the large majority of cases. Consequently, background splicing information can usefully predict the effects of splice site mutations, which include cryptic splice activation and single or multiple exon skipping. In addition, de novo splice sites and splice sites involved in pseudoexon formation, recursive splicing and aberrant splicing in cancer show a 90% match to background splice sites, so establishing that the enhancement of background splicing causes a wide range of splicing aberrations. We also discuss how background splicing information can identify cryptic splice sites that might be usefully targeted by antisense oligonucleotides (ASOs) and how it might indicate possible multiple exon skipping side effects of ASOs designed to induce single exon skipping.

Exploration of Alternative Splicing (AS) Events in MDV-Infected Chicken Spleens

Marek’s disease (MD) was an immunosuppression disease induced by Marek’s disease virus (MDV). MD caused huge economic loss to the global poultry industry, but it also provided an ideal model for studying diseases induced by the oncogenic virus. Alternative splicing (AS) simultaneously produced different isoform transcripts, which are involved in various diseases and individual development. To investigate AS events in MD, RNA-Seq was performed in tumorous spleens (TS), spleens from the survivors (SS) without any lesion after MDV infection, and non-infected chicken spleens (NS). In this study, 32,703 and 25,217 AS events were identified in TS and SS groups with NS group as the control group, and 1198, 1204, and 348 differently expressed (DE) AS events (p-value < 0.05 and FDR < 0.05) were identified in TS vs. NS, TS vs. SS, SS vs. NS, respectively. Additionally, Function enrichment analysis showed that ubiquitin-mediated proteolysis, p53 signaling pathway, and phosphatidylinositol signaling system were significantly enriched (p-value < 0.05). Small structural variations including SNP and indel were analyzed based on RNA-Seq data, and it showed that the TS group possessed more variants on the splice site region than those in SS and NS groups, which might cause more AS events in the TS group. Combined with previous circRNA data, we found that 287 genes could produce both circular and linear RNAs, which suggested these genes were more active in MD lymphoma transformation. This study has expanded the understanding of the MDV infection process and provided new insights for further analysis of resistance/susceptibility mechanisms.

Solid Pseudopapillary Neoplasm of the Pancreas and Abdominal Desmoid Tumor in a Patient Carrying Two Different BRCA2 Germline Mutations: New Horizons from Tumor Molecular Profiling

This case report describes the history of a 41 year-old woman with a solid pseudopapillary neoplasm (SPN) of the pancreas and a metachronous abdominal desmoid tumor (DT) that occurred two years after the SPN surgical resection. At next-generation sequencing of 174 cancer-related genes, both neoplasms harbored a CTNNB1 somatic mutation which was different in each tumor. Moreover, two BRCA2 pathogenic mutations were found in both tumors, confirmed as germline by the sequencing of normal tissue. The BRCA2 mutations were c.631G>A, resulting in the amino-acid change p.V211I, and c.7008-2A>T, causing a splice acceptor site loss. However, as the two neoplasms showed neither loss of heterozygosity nor somatic mutation in the second BRCA2 allele, they cannot be considered as BRCA-dependent tumors. Nevertheless, this study highlights the important opportunities opened by extensive tumor molecular profiling. In this particular case, it permitted the detection of BRCA2-germline mutations, essential for addressing the necessary BRCA-related genetic counseling, surveillance, and screening for the patient and her family.

Variation benchmark datasets: update, criteria, quality and applications

Development of new computational methods and testing their performance has to be carried out using experimental data. Only in comparison to existing knowledge can method performance be assessed. For that purpose, benchmark datasets with known and verified outcome are needed. High-quality benchmark datasets are valuable and may be difficult, laborious and time consuming to generate. VariBench and VariSNP are the two existing databases for sharing variation benchmark datasets used mainly for variation interpretation. They have been used for training and benchmarking predictors for various types of variations and their effects. VariBench was updated with 419 new datasets from 109 papers containing altogether 329 014 152 variants; however, there is plenty of redundancy between the datasets. VariBench is freely available at http://structure.bmc.lu.se/VariBench/. The contents of the datasets vary depending on information in the original source. The available datasets have been categorized into 20 groups and subgroups. There are datasets for insertions and deletions, substitutions in coding and non-coding region, structure mapped, synonymous and benign variants. Effect-specific datasets include DNA regulatory elements, RNA splicing, and protein property for aggregation, binding free energy, disorder and stability. Then there are several datasets for molecule-specific and disease-specific applications, as well as one dataset for variation phenotype effects. Variants are often described at three molecular levels (DNA, RNA and protein) and sometimes also at the protein structural level including relevant cross references and variant descriptions. The updated VariBench facilitates development and testing of new methods and comparison of obtained performances to previously published methods. We compared the performance of the pathogenicity/tolerance predictor PON-P2 to several benchmark studies, and show that such comparisons are feasible and useful, however, there may be limitations due to lack of provided details and shared data.

Database URL: http://structure.bmc.lu.se/VariBench

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.

Functional evaluation of five BRCA2 unclassified variants identified in a Sri Lankan cohort with inherited cancer syndromes using a mouse embryonic stem cell-based assay

Next-generation sequencing of Sri Lankan families with inherited cancer syndromes resulted in the identification of five BRCA2 variants of unknown clinical significance. Interpreting such variants poses significant challenges for both clinicians and patients. Using a mouse embryonic stem cell-based functional assay, we found I785V, N830D, and K2077N to be functionally indistinguishable from wild-type BRCA2. Specific but mild sensitivity to olaparib and reduction in homologous recombination (HR) efficiency suggest partial loss of function of the A262T variant. This variant is located in the N-terminal DNA binding domain of BRCA2 that can facilitate HR by binding to dsDNA/ssDNA junctions. P3039P is clearly pathogenic because of premature protein truncation caused by exon 23 skipping. These findings highlight the value of mouse embryonic stem cell-based assays for determining the functional significance of variants of unknown clinical significance and provide valuable information regarding risk estimation and genetic counseling of families carrying these BRCA2 variants.

Rare germline genetic variants and risk of aggressive prostate cancer

Few genetic risk factors have been demonstrated to be specifically associated with aggressive prostate cancer (PrCa). Here, we report a case-case study of PrCa comparing the prevalence of germline pathogenic/likely pathogenic (P/LP) genetic variants in 787 men with aggressive disease and 769 with nonaggressive disease. Overall, we observed P/LP variants in 11.4% of men with aggressive PrCa and 9.8% of men with nonaggressive PrCa (two-tailed Fisher's exact tests, P = .28). The proportion of BRCA2 and ATM P/LP variant carriers in men with aggressive PrCa exceeded that observed in men with nonaggressive PrCa; 18/787 carriers (2.3%) and 4/769 carriers (0.5%), P = .004, and 14/787 carriers (0.02%) and 5/769 carriers (0.01%), P = .06, respectively. Our findings contribute to the extensive international effort to interpret the genetic variation identified in genes included on gene-panel tests, for which there is currently an insufficient evidence-base for clinical translation in the context of PrCa risk.

Expression Changes Confirm Genomic Variants Predicted to Result in Allele-Specific, Alternative mRNA Splicing

Splice isoform structure and abundance can be affected by either noncoding or masquerading coding variants that alter the structure or abundance of transcripts. When these variants are common in the population, these nonconstitutive transcripts are sufficiently frequent so as to resemble naturally occurring, alternative mRNA splicing. Prediction of the effects of such variants has been shown to be accurate using information theory-based methods. Single nucleotide polymorphisms (SNPs) predicted to significantly alter natural and/or cryptic splice site strength were shown to affect gene expression. Splicing changes for known SNP genotypes were confirmed in HapMap lymphoblastoid cell lines with gene expression microarrays and custom designed q-RT-PCR or TaqMan assays. The majority of these SNPs (15 of 22) as well as an independent set of 24 variants were then subjected to RNAseq analysis using the ValidSpliceMut web beacon (http://validsplicemut.cytognomix.com), which is based on data from the Cancer Genome Atlas and International Cancer Genome Consortium. SNPs from different genes analyzed with gene expression microarray and q-RT-PCR exhibited significant changes in affected splice site use. Thirteen SNPs directly affected exon inclusion and 10 altered cryptic site use. Homozygous SNP genotypes resulting in stronger splice sites exhibited higher levels of processed mRNA than alleles associated with weaker sites. Four SNPs exhibited variable expression among individuals with the same genotypes, masking statistically significant expression differences between alleles. Genome-wide information theory and expression analyses (RNAseq) in tumor exomes and genomes confirmed splicing effects for 7 of the HapMap SNP and 14 SNPs identified from tumor genomes. q-RT-PCR resolved rare splice isoforms with read abundance too low for statistical significance in ValidSpliceMut. Nevertheless, the web-beacon provides evidence of unanticipated splicing outcomes, for example, intron retention due to compromised recognition of constitutive splice sites. Thus, ValidSpliceMut and q-RT-PCR represent complementary resources for identification of allele-specific, alternative splicing.

BRCA and PALB2 mutations in a cohort of male breast cancer with one bilateral case

INTRODUCTION

Male Breast Cancer (MBC) is a rare disease, about 1% of all breast cancers worldwide and less than 1% of cancers occurring in men. The bilateral male breast cancer (bMBC) is extremely rare. Germline mutations of BRCA1/BRCA2 genes are associated with a significantly increased risk of cancer in MBC; the role of PALB2 remains to be clarified. Our main goal was to provide contribution on characterization of BRCA1/BRCA2 and PALB2 mutations in MBC patients.

METHODS

We observed 28 MBC cases; one of them was a bMBC. Screening for BRCA1, BRCA2 and PALB2 genes was performed on all 28 MBC patients. Mutational analysis was extended to family members of mutated patients.

RESULTS

In our study, the MBC incidence was 5.2% and for bMBC was 3.6%. Mutation analysis showed pathogenic mutations in 11/28 (39.3%) patients; 2/28 (7.1%) displayed a mutation in BRCA1, 8/28 (28.6%) in BRCA2 and 1/28 (3.6%) in PALB2. Out of 11 mutated patients, one (9.1%) reported a double mutation in BRCA2. Personal history of other cancers was reported in 2/28 (7.1%) patients affected by bladder cancer. A first/second degree family history of breast/ovarian and other cancers occurred in 23/28 (82.1%) patients.

CONCLUSION

Our findings indicate BRCA2 as the main MBC susceptibility gene and describe an increased risk of bMBC and bladder cancer in mutated patients. The identification of mutations in MBC susceptibility genes supports the usage of oncology prevention programs in affected patients and their relatives carrying the mutation.

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.

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.

Extensive Alternative Splicing of KIR Transcripts

The killer-cell Ig-like receptors (KIR) form a multigene entity involved in modulating immune responses through interactions with MHC class I molecules. The complexity of the KIR cluster is reflected by, for instance, abundant levels of allelic polymorphism, gene copy number variation, and stochastic expression profiles. The current transcriptome study involving human and macaque families demonstrates that KIR family members are also subjected to differential levels of alternative splicing, and this seems to be gene dependent. Alternative splicing may result in the partial or complete skipping of exons, or the partial inclusion of introns, as documented at the transcription level. This post-transcriptional process can generate multiple isoforms from a single KIR gene, which diversifies the characteristics of the encoded proteins. For example, alternative splicing could modify ligand interactions, cellular localization, signaling properties, and the number of extracellular domains of the receptor. In humans, we observed abundant splicing for KIR2DL4, and to a lesser extent in the lineage III KIR genes. All experimentally documented splice events are substantiated by in silico splicing strength predictions. To a similar extent, alternative splicing is observed in rhesus macaques, a species that shares a close evolutionary relationship with humans. Splicing profiles of Mamu-KIR1D and Mamu-KIR2DL04 displayed a great diversity, whereas Mamu-KIR3DL20 (lineage V) is consistently spliced to generate a homolog of human KIR2DL5 (lineage I). The latter case represents an example of convergent evolution. Although just a single KIR splice event is shared between humans and macaques, the splicing mechanisms are similar, and the predicted consequences are comparable. In conclusion, alternative splicing adds an additional layer of complexity to the KIR gene system in primates, and results in a wide structural and functional variety of KIR receptors and its isoforms, which may play a role in health and disease.

Additional molecular and clinical evidence open the way to definitive IARC classification of the BRCA1 c.5332G > A (p.Asp1778Asn) variant

OBJECTIVES

In silico splicing analysis, a mini-gene assay and splicing data, obtained using RNA from blood samples, have shown that the BRCA1 c.5332G > A variant induces exon 21 skipping. However, despite these evidences, up to date, this variant is unclassified. The aim of this study is to provide further molecular and clinical evidence for the BRCA1 c.5332G > A variant in a patient with high grade serous ovarian carcinoma (HGSOC) to allow a definitive classification of this variant.

DESIGN AND METHOD

The effect of the BRCA1 c.5332G > A variant on RNA splicing was evaluated by amplifying regions of BRCA1 from the cDNA of the patient. Loss of heterozygosity (LOH) in tumor tissue was also investigated.

RESULTS

The c.5332G > A allele causes significantly aberrant splicing of the BRCA1 exon 21. Evaluation of the c.5332A allele in tumor tissue highlights a possible loss of heterozygosity, supporting her pathogenic effect.

CONCLUSIONS

Our results regarding the c.5332G > A variant confirm that it contributed to predisposition and onset of ovarian carcinoma in the patient. We propose to classify this variant as 'likely-pathogenic' (class IV).

Germline pathogenic variants of 11 breast cancer genes in 7,051 Japanese patients and 11,241 controls

Pathogenic variants in highly penetrant genes are useful for the diagnosis, therapy, and surveillance for hereditary breast cancer. Large-scale studies are needed to inform future testing and variant classification processes in Japanese. We performed a case-control association study for variants in coding regions of 11 hereditary breast cancer genes in 7051 unselected breast cancer patients and 11,241 female controls of Japanese ancestry. Here, we identify 244 germline pathogenic variants. Pathogenic variants are found in 5.7% of patients, ranging from 15% in women diagnosed <40 years to 3.2% in patients ≥80 years, with BRCA1/2, explaining two-thirds of pathogenic variants identified at all ages. BRCA1/2, PALB2, and TP53 are significant causative genes. Patients with pathogenic variants in BRCA1/2 or PTEN have significantly younger age at diagnosis. In conclusion, BRCA1/2, PALB2, and TP53 are the major hereditary breast cancer genes, irrespective of age at diagnosis, in Japanese women.

An intronic mutation c.6430-3C>G in the F8 gene causes splicing efficiency and premature termination in hemophilia A

: Hemophilia A is a bleeding disorder caused by coagulation factor VIII protein deficiency or dysfunction, which is classified into severe, moderate, and mild according to factor clotting activity. An overwhelming majority of missense and nonsense mutations occur in exons of F8 gene, whereas mutations in introns can also be pathogenic. This study aimed to investigate the effect of an intronic mutation, c.6430-3C>G (IVS22-3C>G), on pre-mRNA splicing of the F8 gene. We applied DNA and cDNA sequencing in a Chinese boy with hemophilia A to search if any pathogenic mutation in the F8 gene. Functional analysis was performed to investigate the effect of an intronic mutation at the transcriptional level. Human Splicing Finder and PyMol were also used to predict its effect. We found the mutation c.6430-3C>G (IVS22-3C>G) in the F8 gene in the affected boy, with his mother being a carrier. cDNA from the mother and pSPL3 splicing assay showed that the mutation IVS22-3C>G results in a two-nucleotide AG inclusion at the 3' end of intron 22 and leads to a truncated coagulation factor VIII protein, with partial loss of the C1 domain and complete loss of the C2 domain. The in-silico tool predicted that the mutation induces altered pre-mRNA splicing by using a cryptic acceptor site in intron 22. The IVS22-3C>G mutation was confirmed to affect pre-mRNA splicing and produce a truncated protein, which reduces the stability of binding between the F8 protein and von Willebrand factor carrier protein due to the loss of an interaction domain.

Computational Tools for Splicing Defect Prediction in Breast/Ovarian Cancer Genes: How Efficient Are They at Predicting RNA Alterations?

In silico tools for splicing defect prediction have a key role to assess the impact of variants of uncertain significance. Our aim was to evaluate the performance of a set of commonly used splicing in silico tools comparing the predictions against RNA in vitro results. This was done for natural splice sites of clinically relevant genes in hereditary breast/ovarian cancer (HBOC) and Lynch syndrome. A study divided into two stages was used to evaluate SSF-like, MaxEntScan, NNSplice, HSF, SPANR, and dbscSNV tools. A discovery dataset of 99 variants with unequivocal results of RNA in vitro studies, located in the 10 exonic and 20 intronic nucleotides adjacent to exon-intron boundaries of BRCA1, BRCA2, MLH1, MSH2, MSH6, PMS2, ATM, BRIP1, CDH1, PALB2, PTEN, RAD51D, STK11, and TP53, was collected from four Spanish cancer genetic laboratories. The best stand-alone predictors or combinations were validated with a set of 346 variants in the same genes with clear splicing outcomes reported in the literature. Sensitivity, specificity, accuracy, negative predictive value (NPV) and Mathews Coefficient Correlation (MCC) scores were used to measure the performance. The discovery stage showed that HSF and SSF-like were the most accurate for variants at the donor and acceptor region, respectively. The further combination analysis revealed that HSF, HSF+SSF-like or HSF+SSF-like+MES achieved a high performance for predicting the disruption of donor sites, and SSF-like or a sequential combination of MES and SSF-like for predicting disruption of acceptor sites. The performance confirmation of these last results with the validation dataset, indicated that the highest sensitivity, accuracy, and NPV (99.44%, 99.44%, and 96.88, respectively) were attained with HSF+SSF-like or HSF+SSF-like+MES for donor sites and SSF-like (92.63%, 92.65%, and 84.44, respectively) for acceptor sites. We provide recommendations for combining algorithms to conduct in silico splicing analysis that achieved a high performance. The high NPV obtained allows to select the variants in which the study by in vitro RNA analysis is mandatory against those with a negligible probability of being spliceogenic. Our study also shows that the performance of each specific predictor varies depending on whether the natural splicing sites are donors or acceptors.

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.

RNA analysis of cancer predisposing genes in formalin-fixed paraffin-embedded tissue determines aberrant splicing

High-throughput sequencing efforts in molecular tumour diagnostics detect increasing numbers of novel variants, including variants predicted to affect splicing. In silico prediction tools can reliably predict the effect of variant disrupting canonical splice sites; however, experimental validation is required to confirm aberrant splicing. Here, we present RNA analysis performed for 13 canonical splice site variants predicted or known to result in splicing in the cancer predisposition genes MLH1, MSH2, MSH6, APC and BRCA1. Total nucleic acid was successfully isolated for 10 variants from eight formalin-fixed paraffin-embedded (FFPE) tumour tissues and two B-cell lines. Aberrant splicing was confirmed in all six variants known to result in splicing. Of one known variant in the B-cell line, aberrant splicing could only be detected after formalin fixation, which indicated that formalin fixation could possibly inhibit RNA degradation. Aberrant splicing was concluded in three of four predicted splice variants of uncertain significance, supporting their pathogenic effect. With this assay, somatic splice variants can be easily and rapidly analysed, enabling retrospective analysis to support the pathogenicity of variants predicted to result in splicing when only FFPE material is available.

Characterization of a novel germline BRCA1 splice variant, c.5332+4delA

PURPOSE

Germline mutations in BRCA1 and BRCA2 confer a significant increase in risk for cancer, and determining pathogenicity of a BRCA variant can guide the clinical management of the disease. About 1/3 of BRCA1 variants reported in the public databases have uncertain clinical significance due to lack of conclusive evidence. This study aims to characterize a novel BRCA1 deletion affecting the + 4 splice donor site identified in an individual with early-onset breast cancer.

METHODS

The effect of BRCA1 c.5332+4delA variant on RNA splicing was evaluated by amplifying regions of BRCA1 from cDNA derived from the patient. The proportion of abnormal transcript in the total transcripts was quantified. Loss of heterozygosity (LOH) in tumor tissue was investigated using Sanger sequencing and fragment analysis.

RESULTS

BRCA1 c.5332+4delA caused skipping of exon 21 in patient-derived samples. Semi-quantitative analysis indicated that this aberrant RT-PCR product accounts for about 40% of the total transcript levels. Loss of heterozygosity (LOH) was observed in patient's tumor tissue.

CONCLUSIONS

Our results indicate that the BRCA1 c.5332+4delA variant contributes to cancer predisposition through disruption of normal mRNA splicing. We classify this variant as likely pathogenic.

Leveraging splice-affecting variant predictors and a minigene validation system to identify Mendelian disease-causing variants among exon-captured variants of uncertain significance

The genetic heterogeneity of Mendelian disorders results in a significant proportion of patients that are unable to be assigned a confident molecular diagnosis after conventional exon sequencing and variant interpretation. Here, we evaluated how many patients with an inherited retinal disease (IRD) have variants of uncertain significance (VUS) that are disrupting splicing in a known IRD gene by means other than affecting the canonical dinucleotide splice site. Three in silico splice-affecting variant predictors were leveraged to annotate and prioritize variants for splicing functional validation. An in vitro minigene system was used to assay each variant's effect on splicing. Starting with 745 IRD patients lacking a confident molecular diagnosis, we validated 23 VUS as splicing variants that likely explain disease in 26 patients. Using our results, we optimized in silico score cutoffs to guide future variant interpretation. Variants that alter base pairs other than the canonical GT-AG dinucleotide are often not considered for their potential effect on RNA splicing but in silico tools and a minigene system can be utilized for the prioritization and validation of such splice-disrupting variants. These variants can be overlooked causes of human disease but can be identified using conventional exon sequencing with proper interpretation guidelines.

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.

Classification of the spliceogenic BRCA1 c.4096+3A>G variant as likely benign based on cosegregation data and identification of a healthy homozygous carrier

BRCA1, c.4096+3A>G was identified in a consanguineous Danish family with several cases of breast/ovarian cancer. In silico analysis and splicing assays indicated that the variant caused aberrant splicing. However, based on segregation data and the finding of a healthy homozygous carrier, we classify the BRCA1 c.4096+3A>G variant as likely benign.

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, BRCA2 and PALB2 mutations and CHEK2 c.1100delC in different South African ethnic groups diagnosed with premenopausal and/or triple negative breast cancer

Background

Current knowledge of the aetiology of hereditary breast cancer in the four main South African population groups (black, coloured, Indian and white) is limited. Risk assessments in the black, coloured and Indian population groups are challenging because of restricted information regarding the underlying genetic contributions to inherited breast cancer in these populations. We focused this study on premenopausal patients (diagnosed with breast cancer before the age of 50; n = 78) and triple negative breast cancer (TNBC) patients (n = 30) from the four South African ethnic groups. The aim of this study was to determine the frequency and spectrum of germline mutations in BRCA1, BRCA2 and PALB2 and to evaluate the presence of the CHEK2 c.1100delC allele in these patients.

Methods

In total, 108 South African breast cancer patients underwent mutation screening using a Next-Generation Sequencing (NGS) approach in combination with Multiplex Ligation-dependent Probe Amplification (MLPA) to detect large rearrangements in BRCA1 and BRCA2.

Results

In 13 (12 %) patients a deleterious mutation in BRCA1/2 was detected, three of which were novel mutations in black patients. None of the study participants was found to have an unequivocal pathogenic mutation in PALB2. Two (white) patients tested positive for the CHEK2 c.1100delC mutation, however, one of these also carried a deleterious BRCA2 mutation. Additionally, six variants of unknown clinical significance were identified (4 in BRCA2, 2 in PALB2), all in black patients. Within the group of TNBC patients, a higher mutation frequency was obtained (23.3 %; 7/30) than in the group of patients diagnosed before the age of 50 (7.7 %; 6/78).

Conclusion

This study highlights the importance of evaluating germline mutations in major breast cancer genes in all of the South African population groups. This NGS study shows that mutation analysis is warranted in South African patients with triple negative and/or in premenopausal breast cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1913-6) contains supplementary material, which is available to authorized users.

FANCM c.5791C>T nonsense mutation (rs144567652) induces exon skipping, affects DNA repair activity and is a familial breast cancer risk factor

Numerous genetic factors that influence breast cancer risk are known. However, approximately two-thirds of the overall familial risk remain unexplained. To determine whether some of the missing heritability is due to rare variants conferring high to moderate risk, we tested for an association between the c.5791C>T nonsense mutation (p.Arg1931*; rs144567652) in exon 22 of FANCM gene and breast cancer. An analysis of genotyping data from 8635 familial breast cancer cases and 6625 controls from different countries yielded an association between the c.5791C>T mutation and breast cancer risk [odds ratio (OR) = 3.93 (95% confidence interval (CI) = 1.28-12.11; P = 0.017)]. Moreover, we performed two meta-analyses of studies from countries with carriers in both cases and controls and of all available data. These analyses showed breast cancer associations with OR = 3.67 (95% CI = 1.04-12.87; P = 0.043) and OR = 3.33 (95% CI = 1.09-13.62; P = 0.032), respectively. Based on information theory-based prediction, we established that the mutation caused an out-of-frame deletion of exon 22, due to the creation of a binding site for the pre-mRNA processing protein hnRNP A1. Furthermore, genetic complementation analyses showed that the mutation influenced the DNA repair activity of the FANCM protein. In summary, we provide evidence for the first time showing that the common p.Arg1931* loss-of-function variant in FANCM is a risk factor for familial breast cancer.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Splicing analysis of 14 BRCA1 missense variants classifies nine variants as pathogenic

Pathogenic germline mutations in the BRCA1 gene predispose carriers to early onset breast and ovarian cancer. Clinical genetic screening of BRCA1 often reveals variants with uncertain clinical significance, complicating patient and family management. Therefore, functional examinations are urgently needed to classify whether these uncertain variants are pathogenic or benign. In this study, we investigated 14 BRCA1 variants by in silico splicing analysis and mini-gene splicing assay. All 14 alterations were missense variants located within the BRCT domain of BRCA1 and had previously been examined by functional analysis at the protein level. Results from a validated mini-gene splicing assay indicated that nine BRCA1 variants resulted in splicing aberrations leading to truncated transcripts and thus can be considered pathogenic (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5072C>T/p.Thr1691Ile, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala), whereas five BRCA1 variants had no effect on splicing (c.4985T>C/p.Phe1662Ser, c.5072C>A/p.Thr1691Lys, c.5153G>C/p.Trp1718Ser, c.5154G>T/p.Trp1718Cys, and c.5333A>G/p.Asp1778Gly). Eight of the variants having an effect on splicing (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala) were previously determined to have no or an uncertain effect on the protein level, whereas one variant (c.5072C>T/p.Thr1691Ile) were shown to have a strong effect on the protein level as well. In conclusion, our study emphasizes that in silico splicing prediction and mini-gene splicing analysis are important for the classification of BRCA1 missense variants located close to exon/intron boundaries.

Diversity and impact of rare variants in genes encoding the platelet G protein-coupled receptors

Platelet responses to activating agonists are influenced by common population variants within or near G protein-coupled receptor (GPCR) genes that affect receptor activity. However, the impact of rare GPCR gene variants is unknown. We describe the rare single nucleotide variants (SNVs) in the coding and splice regions of 18 GPCR genes in 7,595 exomes from the 1,000-genomes and Exome Sequencing Project databases and in 31 cases with inherited platelet function disorders (IPFDs). In the population databases, the GPCR gene target regions contained 740 SNVs (318 synonymous, 410 missense, 7 stop gain and 6 splice region) of which 70 % had global minor allele frequency (MAF) < 0.05 %. Functional annotation using six computational algorithms, experimental evidence and structural data identified 156/740 (21 %) SNVs as potentially damaging to GPCR function, most commonly in regions encoding the transmembrane and C-terminal intracellular receptor domains. In 31 index cases with IPFDs (Gi-pathway defect n=15; secretion defect n=11; thromboxane pathway defect n=3 and complex defect n=2) there were 256 SNVs in the target regions of 15 stimulatory platelet GPCRs (34 unique; 12 with MAF< 1 % and 22 with MAF≥ 1 %). These included rare variants predicting R122H, P258T and V207A substitutions in the P2Y12 receptor that were annotated as potentially damaging, but only partially explained the platelet function defects in each case. Our data highlight that potentially damaging variants in platelet GPCR genes have low individual frequencies, but are collectively abundant in the population. Potentially damaging variants are also present in pedigrees with IPFDs and may contribute to complex laboratory phenotypes.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

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.

PALB2 sequencing in Italian familial breast cancer cases reveals a high-risk mutation recurrent in the province of Bergamo

PURPOSE

Monoallelic germ-line deleterious mutations of PALB2 (partner and localizer of BRCA2) are associated with breast cancer risk and have been found in several populations, with carrier frequencies of ~1-2%. Initially, these mutations were considered to have moderate penetrance, but accumulating evidence now indicates that they are associated with much higher risk.

METHODS

In this study, we sequenced the PALB2 coding regions unlinked to BRCA (breast cancer) genes in 575 probands from Italian breast cancer families recruited in Milan.

RESULTS

We found 12 carriers (2.1%) of deleterious mutations, and none of the mutations was found in 784 controls collected in Milan. One of these mutations, the c.1027C>T (p.Gln343X), was found to be recurrent in the province of Bergamo in northern Italy, being detected in 6/113 (5.3%) familial breast cancer cases and 2/477 (0.4%) controls recruited in this area (Fisher's exact test: P < 0.01).

CONCLUSIONS

Our data provide confirmatory findings that, in the Italian population also, deleterious mutations of PALB2 are relatively frequent predisposing factors for breast cancer and may be associated with high risk of the disease.

Characterization of an Italian founder mutation in the RING-finger domain of BRCA1

The identification of founder mutations in cancer predisposing genes is important to improve risk assessment in geographically defined populations, since it may provide specific targets resulting in cost-effective genetic testing. Here, we report the characterization of the BRCA1 c.190T>C (p.Cys64Arg) mutation, mapped to the RING-finger domain coding region, that we detected in 43 hereditary breast/ovarian cancer (HBOC) families, for the large part originating from the province of Bergamo (Northern Italy). Haplotype analysis was performed in 21 families, and led to the identification of a shared haplotype extending over three BRCA1-associated marker loci (0.4 cM). Using the DMLE+2.2 software program and regional population demographic data, we were able to estimate the age of the mutation to vary between 3,100 and 3,350 years old. Functional characterization of the mutation was carried out at both transcript and protein level. Reverse transcriptase-PCR analysis on lymphoblastoid cells revealed expression of full length mRNA from the mutant allele. A green fluorescent protein (GFP)-fragment reassembly assay showed that the p.Cys64Arg substitution prevents the binding of the BRCA1 protein to the interacting protein BARD1, in a similar way as proven deleterious mutations in the RING-domain. Overall, 55 of 83 (66%) female mutation carriers had a diagnosis of breast and/or ovarian cancer. Our observations indicate that the BRCA1 c.190T>C is a pathogenic founder mutation present in the Italian population. Further analyses will evaluate whether screening for this mutation can be suggested as an effective strategy for the rapid identification of at-risk individuals in the Bergamo area.

In silico tools for splicing defect prediction: a survey from the viewpoint of end users

RNA splicing is the process during which introns are excised and exons are spliced. The precise recognition of splicing signals is critical to this process and mutations affecting splicing comprise a considerable proportion of genetic disease etiology. Analysis of RNA samples from the patient is the most straightforward and reliable method to detect splicing defects. However, currently the technical limitation prohibits its use in routine clinical practice. In silico tools that predict potential consequences of splicing mutations may be useful in daily diagnostic activities. In this review, we provide medical geneticists with some basic insights into some of the most popular in silico tools for splicing defect prediction, from the viewpoint of end-users. Bioinformaticians in relevant areas who are working on huge datasets may also benefit from this review. Specifically, we focus on those tools whose primary goal is to predict the impact of mutations within the 5′ and 3′ splicing consensus regions: the algorithms used by different tools as well as their major advantages and disadvantages are briefly introduced; the formats of their input and output are summarized; and the interpretation, evaluation, and prospection are also discussed.

Comparison of mRNA splicing assay protocols across multiple laboratories: recommendations for best practice in standardized clinical testing

BACKGROUND

Accurate evaluation of unclassified sequence variants in cancer predisposition genes is essential for clinical management and depends on a multifactorial analysis of clinical, genetic, pathologic, and bioinformatic variables and assays of transcript length and abundance. The integrity of assay data in turn relies on appropriate assay design, interpretation, and reporting.

METHODS

We conducted a multicenter investigation to compare mRNA splicing assay protocols used by members of the ENIGMA (Evidence-Based Network for the Interpretation of Germline Mutant Alleles) consortium. We compared similarities and differences in results derived from analysis of a panel of breast cancer 1, early onset (BRCA1) and breast cancer 2, early onset (BRCA2) gene variants known to alter splicing (BRCA1: c.135-1G>T, c.591C>T, c.594-2A>C, c.671-2A>G, and c.5467+5G>C and BRCA2: c.426-12_8delGTTTT, c.7988A>T, c.8632+1G>A, and c.9501+3A>T). Differences in protocols were then assessed to determine which elements were critical in reliable assay design.

RESULTS

PCR primer design strategies, PCR conditions, and product detection methods, combined with a prior knowledge of expected alternative transcripts, were the key factors for accurate splicing assay results. For example, because of the position of primers and PCR extension times, several isoforms associated with BRCA1, c.594-2A>C and c.671-2A>G, were not detected by many sites. Variation was most evident for the detection of low-abundance transcripts (e.g., BRCA2 c.8632+1G>A Δ19,20 and BRCA1 c.135-1G>T Δ5q and Δ3). Detection of low-abundance transcripts was sometimes addressed by using more analytically sensitive detection methods (e.g., BRCA2 c.426-12_8delGTTTT ins18bp).

CONCLUSIONS

We provide recommendations for best practice and raise key issues to consider when designing mRNA assays for evaluation of unclassified sequence variants.

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