Splicing and multifactorial analysis of intronic BRCA1 and BRCA2 sequence variants identifies clinically significant splicing aberrations up to 12 nucleotides from the intron/exon boundary

Clinical management of breast cancer families is complicated by identification of BRCA1 and BRCA2 sequence alterations of unknown significance. Molecular assays evaluating the effect of intronic variants on native splicing can help determine their clinical relevance. Twenty-six intronic BRCA1/2 variants ranging from the consensus dinucleotides in the splice acceptor or donor to 53 nucleotides into the intron were identified in multiple-case families. The effect of the variants on splicing was assessed using HSF matrices, MaxEntScan and NNsplice, followed by analysis of mRNA from lymphoblastoid cell lines. A total of 12 variants were associated with splicing aberrations predicted to result in production of truncated proteins, including a variant located 12 nucleotides into the intron. The posterior probability of pathogenicity was estimated using a multifactorial likelihood approach, and provided a pathogenic or likely pathogenic classification for seven of the 12 spliceogenic variants. The apparent disparity between experimental evidence and the multifactorial predictions is likely due to several factors, including a paucity of likelihood information and a nonspecific prior probability applied for intronic variants outside the consensus dinucleotides. Development of prior probabilities of pathogenicity incorporating bioinformatic prediction of splicing aberrations should improve identification of functionally relevant variants and enhance multifactorial likelihood analysis of intronic variants.     

Molecular Outcome,Prediction, and Clinical Consequences of Splice Variants in COL1A1, Which Encodes the proα1(I) Chains of Type I Procollagen

Approximately 10%–20% of germline pathogenic variants alter mRNA splicing, with phenotypes often dependent on the stability of the mRNA produced by the mutant allele. To better understand the relationships between genotype, mRNA splicing, and phenotype, we examined clinical and molecular data from 243 probands with osteogenesis imperfecta (OI) representing 145 unique splicing variants within the type I procollagen gene, COL1A1. All individuals with IVSX‐1G>A mutations had OI type I because the substitution shifted the splice acceptor site 1 nt downstream and destabilized the mRNA. OI phenotypes were not consistent for any other splice variant identified. We sequenced all cDNA species from cultured dermal fibroblasts from 40 individuals to identify splice outcome and compared those results to splice predictions from Human Splice Finder (HSF), Spliceport (SP), and Automatic Splice Site and Exon Definition Analyses (ASSEDA). Software‐based splice predictions were correct in 42%, 55%, and 74% instances for HSF, SP, and ASSEDA, respectively. As molecular diagnostics move increasingly to DNA sequence analysis, the need to understand the effects of splice site variants will increase. These data demonstrate that caution must be exercised when using splice prediction software to predict splice outcome.     

Classifying MLH1 and MSH2 variants using bioinformatic prediction,splicing assays,segregation, and tumor characteristics

Reliable methods for predicting functional consequences of variants in disease genes would be beneficial in the clinical setting. This study was undertaken to predict, and confirm in vitro, splicing aberrations associated with mismatch repair (MMR) variants identified in familial colon cancer patients. Six programs were used to predict the effect of 13 MLH1 and 6 MSH2 gene variants on pre‐mRNA splicing. mRNA from cycloheximide‐treated lymphoblastoid cell lines of variant carriers was screened for splicing aberrations. Tumors of variant carriers were tested for microsatellite instability and MMR protein expression. Variant segregation in families was assessed using Bayes factor causality analysis. Amino acid alterations were examined for evolutionary conservation and physicochemical properties. Splicing aberrations were detected for 10 variants, including a frameshift as a minor cDNA product, and altered ratio of known alternate splice products. Loss of splice sites was well predicted by splice‐site prediction programs SpliceSiteFinder (90%) and NNSPLICE (90%), but consequence of splice site loss was less accurately predicted. No aberrations correlated with ESE predictions for the nine exonic variants studied. Seven of eight missense variants had normal splicing (88%), but only one was a substitution considered neutral from evolutionary/physicochemical analysis. Combined with information from tumor and segregation analysis, and literature review, 16 of 19 variants were considered clinically relevant. Bioinformatic tools for prediction of splicing aberrations need improvement before use without supporting studies to assess variant pathogenicity. Classification of mismatch repair gene variants is assisted by a comprehensive approach that includes in vitro, tumor pathology, clinical, and evolutionary conservation data. Hum Mutat 0, 1–14, 2009. © 2009 Wiley‐Liss, Inc.     

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.     

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.     

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.     

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.     

Exceptions to the rule: Case studies in the prediction of pathogenicity for genetic variants in hereditary cancer genes

Based on current consensus guidelines and standard practice, many genetic variants detected in clinical testing are classified as disease causing based on their predicted impact on the normal expression or function of the gene in the absence of additional data. However, our laboratory has identified a subset of such variants in hereditary cancer genes for which compelling contradictory evidence emerged after the initial evaluation following the first observation of the variant. Three representative examples of variants in BRCA1, BRCA2 and MSH2 that are predicted to disrupt splicing, prematurely truncate the protein, or remove the start codon were evaluated for pathogenicity by analyzing clinical data with multiple classification algorithms. Available clinical data for all three variants contradicts the expected pathogenic classification. These variants illustrate potential pitfalls associated with standard approaches to variant classification as well as the challenges associated with monitoring data, updating classifications, and reporting potentially contradictory interpretations to the clinicians responsible for translating test outcomes to appropriate clinical action. It is important to address these challenges now as the model for clinical testing moves toward the use of large multi‐gene panels and whole exome/genome analysis, which will dramatically increase the number of genetic variants identified.     

Prediction of Mutant mRNA Splice Isoforms by Information Theory‐Based Exon Definition

Mutations that affect mRNA splicing often produce multiple mRNA isoforms, resulting in complex molecular phenotypes. Definition of an exon and its inclusion in mature mRNA relies on joint recognition of both acceptor and donor splice sites. This study predicts cryptic and exon‐skipping isoforms in mRNA produced by splicing mutations from the combined information contents (R_i, which measures binding‐site strength, in bits) and distribution of the splice sites defining these exons. The total information content of an exon (R_i,total) is the sum of the R_i values of its acceptor and donor splice sites, adjusted for the self‐information of the distance separating these sites, that is, the gap surprisal. Differences between total information contents of an exon (ΔR_i,total) are predictive of the relative abundance of these exons in distinct processed mRNAs. Constraints on splice site and exon selection are used to eliminate nonconforming and poorly expressed isoforms. Molecular phenotypes are computed by the Automated Splice Site and Exon Definition Analysis ( http://splice.uwo.ca ) server. Predictions of splicing mutations were highly concordant (85.2%; n = 61) with published expression data. In silico exon definition analysis will contribute to streamlining assessment of abnormal and normal splice isoforms resulting from mutations.     

Prevalence of pathogenic germline cancer risk variants in high-risk urothelial carcinoma

PurposeTo date, there has not been a large, systematic evaluation of the prevalence of germline risk variants in urothelial carcinoma (UC).MethodsWe evaluated the frequency of germline pathogenic and likely pathogenic variants in 1038 patients with high-risk UC who underwent targeted clinical germline testing. Case–control enrichment analysis was performed to screen for pathogenic variant enrichment in 17 DNA repair genes in 1038 UC patients relative to cancer-free individuals.ResultsAmong 1038 patients with UC, the cumulative frequency of patients with pathogenic variants was 24%; 18.6% of patients harbored ≥1 actionable germline variant with preventive or therapeutic utility. MSH2 (34/969, 3.5%) and BRCA1/2 (38/867, 4.4%) germline variants had the highest frequency. Germline variants in DNA damage repair genes accounted for 78% of pathogenic germline variants. Compared to the cancer-free cohort, UC patients had significant variant enrichment in MSH2 (odds ratio [OR]: 15.4, 95% confidence interval [CI]: 7.1–32.7, p < 0.0001), MLH1 (OR: 15.9, 95% CI: 4.4–67.7, p < 0.0001), BRCA2 (OR: 5.7, 95% CI: 3.2–9.6, p < 0.0001), and ATM (OR: 3.8, 95% CI: 1.8–8.3, p = 0.02).ConclusionIn this study, 24% of UC patients harbored pathogenic germline variants and 18.6% had clinically actionable variants. MLH1 and MSH2 were validated as UC risk genes while ATM and BRCA2 were highlighted as potential UC predisposition genes. This work emphasizes the utility of germline testing in selected high-risk UC cohorts.     

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.     

An evaluation of the challenges to developing tumor BRCA1 and BRCA2 testing methodologies for clinical practice

Ovarian cancer patients with germline or somatic pathogenic variants benefit from treatment with poly ADP ribose polymerase (PARP) inhibitors. Tumor BRCA1/2 testing is more challenging than germline testing as the majority of samples are formalin‐fixed paraffin embedded (FFPE), the tumor genome is complex, and the allelic fraction of somatic variants can be low. We collaborated with 10 laboratories testing BRCA1/2 in tumors to compare different approaches to identify clinically important variants within FFPE tumor DNA samples. This was not a proficiency study but an inter‐laboratory comparison to identify common issues. Each laboratory received the same tumor DNA samples ranging in genotype, quantity, quality, and variant allele frequency (VAF). Each laboratory performed their preferred next‐generation sequencing method to report on the variants. No false positive results were reported in this small study and the majority of methods detected the low VAF variants. A number of variants were not detected due to the bioinformatics analysis, variant classification, or insufficient DNA. The use of hybridization capture or short amplicon methods are recommended based on a bioinformatic assessment of the data. The study highlights the importance of establishing standards and standardization for tBRCA testing particularly when the test results dictate clinical decisions regarding life extending therapies.     

Comprehensive profiling of BRCA1 and BRCA2 variants in breast and ovarian cancer in Chinese patients

The identification and interpretation of germline BRCA1/2 variants become increasingly important in breast and ovarian cancer (OC) treatment. However, there is no comprehensive analysis of the germline BRCA1/2 variants in a Chinese population. Here we performed a systematic review and meta‐analysis on such variants from 94 publications. A total of 2,128 BRCA1/2 variant records were extracted, including 601 from BRCA1 and 632 from BRCA2. In addition, 414, 734, 449, and 307 variants were also recorded in the BIC, ClinVar, ENIGMA, and UMD databases, respectively, and 579 variants were newly reported. Subsequent analysis showed that the overall germline BRCA1/2 pathogenic variant frequency was 5.7% and 21.8% in Chinese breast and OC, respectively. Populations with high‐risk factors exhibited a higher pathogenic variant percentage. Furthermore, the variant profile in Chinese is distinct from that in other ethnic groups with no distinct founder pathogenic variants. We also tested our in‐house American College of Medical Genetics‐guided pathogenicity interpretation procedure for Chinese BRCA1/2 variants. Our results achieved a consistency of 91.2–97.6% (5‐grade classification) or 98.4–100% (2‐grade classification) with public databases. In conclusion, this study represents the first comprehensive meta‐analysis of Chinese BRCA1/2 variants and validates our in‐house pathogenicity interpretation procedure, thereby providing guidance for further PARP inhibitor development and companion diagnostics in the Chinese population.