Adding In Silico Assessment of Potential Splice Aberration to the Integrated Evaluation of BRCA Gene Unclassified Variants

Meta-EA: a gene-specific combination of available computational tools for predicting missense variant effects

Computational methods for estimating missense variant impact suffer from inconsistent performance across genes, which poses a major challenge for their reliable use in clinical practice. While ensemble scores leverage multiple prediction methods to enhance consistency, the overrepresentation of certain genes in the training data can bias their outcomes. To address this critical limitation, we propose a gene-specific ensemble framework trained on reference computational annotations rather than on clinical or experimental data. Accordingly, we generate Meta-EA ensemble scores that achieve comparable performance to the top individual predicting method for each gene set. Incorporating the effects of splicing and the allele frequency of human polymorphisms further enhances the performance of Meta-EA, achieving an area under the receiver operating characteristic curve of 0.97 for both gene-balanced and imbalanced clinical assessments. In conclusion, this work leverages the wealth of existing variant impact prediction approaches to generate improved estimations for clinical interpretation.

Sequencing-based functional assays for classification of BRCA2 variants in mouse ESCs

Sequencing of genes, such as BRCA1 and BRCA2, is recommended for individuals with a personal or family history of early onset and/or bilateral breast and/or ovarian cancer or a history of male breast cancer. Such sequencing efforts have resulted in the identification of more than 17,000 BRCA2 variants. The functional significance of most variants remains unknown; consequently, they are called variants of uncertain clinical significance (VUSs). We have previously developed mouse embryonic stem cell (mESC)-based assays for functional classification of BRCA2 variants. We now developed a next-generation sequencing (NGS)-based approach for functional evaluation of BRCA2 variants using pools of mESCs expressing 10-25 BRCA2 variants from a given exon. We use this approach for functional evaluation of 223 variants listed in ClinVar. Our functional classification of BRCA2 variants is concordant with the classification reported in ClinVar or those reported by other orthogonal assays.

PDIVAS: Pathogenicity predictor for Deep-Intronic Variants causing Aberrant Splicing

BACKGROUND

Deep-intronic variants that alter RNA splicing were ineffectively evaluated in the search for the cause of genetic diseases. Determination of such pathogenic variants from a vast number of deep-intronic variants (approximately 1,500,000 variants per individual) represents a technical challenge to researchers. Thus, we developed a Pathogenicity predictor for Deep-Intronic Variants causing Aberrant Splicing (PDIVAS) to easily detect pathogenic deep-intronic variants.

RESULTS

PDIVAS was trained on an ensemble machine-learning algorithm to classify pathogenic and benign variants in a curated dataset. The dataset consists of manually curated pathogenic splice-altering variants (SAVs) and commonly observed benign variants within deep introns. Splicing features and a splicing constraint metric were used to maximize the predictive sensitivity and specificity, respectively. PDIVAS showed an average precision of 0.92 and a maximum MCC of 0.88 in classifying these variants, which were the best of the previous predictors. When PDIVAS was applied to genome sequencing analysis on a threshold with 95% sensitivity for reported pathogenic SAVs, an average of 27 pathogenic candidates were extracted per individual. Furthermore, the causative variants in simulated patient genomes were more efficiently prioritized than the previous predictors.

CONCLUSION

Incorporating PDIVAS into variant interpretation pipelines will enable efficient detection of disease-causing deep-intronic SAVs and contribute to improving the diagnostic yield. PDIVAS is publicly available at https://github.com/shiro-kur/PDIVAS .

Functional Analyses of Rare Germline Missense BRCA1 Variants Located within and outside Protein Domains with Known Functions

The BRCA1 protein is implicated in numerous important cellular processes to prevent genomic instability and tumorigenesis, and pathogenic germline variants predispose carriers to hereditary breast and ovarian cancer (HBOC). Most functional studies of missense variants in BRCA1 focus on variants located within the Really Interesting New Gene (RING), coiled-coil and BRCA1 C-terminal (BRCT) domains, and several missense variants in these regions have been shown to be pathogenic. However, the majority of these studies focus on domain specific assays, and have been performed using isolated protein domains and not the full-length BRCA1 protein. Furthermore, it has been suggested that BRCA1 missense variants located outside domains with known function are of no functional importance, and could be classified as (likely) benign. However, very little is known about the role of the regions outside the well-established domains of BRCA1, and only a few functional studies of missense variants located within these regions have been published. In this study, we have, therefore, functionally evaluated the effect of 14 rare BRCA1 missense variants considered to be of uncertain clinical significance, of which 13 are located outside the well-established domains and one within the RING domain. In order to investigate the hypothesis stating that most BRCA1 variants located outside the known protein domains are benign and of no functional importance, multiple protein assays including protein expression and stability, subcellular localisation and protein interactions have been performed, utilising the full-length protein to better mimic the native state of the protein. Two variants located outside the known domains (p.Met297Val and p.Asp1152Asn) and one variant within the RING domain (p.Leu52Phe) were found to make the BRCA1 protein more prone to proteasome-mediated degradation. In addition, two variants (p.Leu1439Phe and p.Gly890Arg) also located outside known domains were found to have reduced protein stability compared to the wild type protein. These findings indicate that variants located outside the RING, BRCT and coiled-coiled domains could also affect the BRCA1 protein function. For the nine remaining variants, no significant effects on BRCA1 protein functions were observed. Based on this, a reclassification of seven variants from VUS to likely benign could be suggested.

Comprehensive evaluation and efficient classification of BRCA1 RING domain missense substitutions

BRCA1 is a high-risk susceptibility gene for breast and ovarian cancer. Pathogenic protein-truncating variants are scattered across the open reading frame, but all known missense substitutions that are pathogenic because of missense dysfunction are located in either the amino-terminal RING domain or the carboxy-terminal BRCT domain. Heterodimerization of the BRCA1 and BARD1 RING domains is a molecularly defined obligate activity. Hence, we tested every BRCA1 RING domain missense substitution that can be created by a single nucleotide change for heterodimerization with BARD1 in a mammalian two-hybrid assay. Downstream of the laboratory assay, we addressed three additional challenges: assay calibration, validation thereof, and integration of the calibrated results with other available data, such as computational evidence and patient/population observational data to achieve clinically applicable classification. Overall, we found that 15%-20% of BRCA1 RING domain missense substitutions are pathogenic. Using a Bayesian point system for data integration and variant classification, we achieved clinical classification of 89% of observed missense substitutions. Moreover, among missense substitutions not present in the human observational data used here, we find an additional 45 with concordant computational and functional assay evidence in favor of pathogenicity plus 223 with concordant evidence in favor of benignity; these are particularly likely to be classified as likely pathogenic and likely benign, respectively, once human observational data become available.

Risks of breast and ovarian cancer for women harboring pathogenic missense variants in BRCA1 and BRCA2 compared with those harboring protein truncating variants

PURPOSE

Germline genetic testing for BRCA1 and BRCA2 variants has been a part of clinical practice for >2 decades. However, no studies have compared the cancer risks associated with missense pathogenic variants (PVs) with those associated with protein truncating (PTC) variants.

METHODS

We collected 582 informative pedigrees segregating 1 of 28 missense PVs in BRCA1 and 153 pedigrees segregating 1 of 12 missense PVs in BRCA2. We analyzed 324 pedigrees with PTC variants in BRCA1 and 214 pedigrees with PTC variants in BRCA2. Cancer risks were estimated using modified segregation analysis.

RESULTS

Estimated breast cancer risks were markedly lower for women aged >50 years carrying BRCA1 missense PVs than for the women carrying BRCA1 PTC variants (hazard ratio [HR] = 3.9 [2.4-6.2] for PVs vs 12.8 [5.7-28.7] for PTC variants; P = .01), particularly for missense PVs in the BRCA1 C-terminal domain (HR = 2.8 [1.4-5.6]; P = .005). In case of BRCA2, for women aged >50 years, the HR was 3.9 (2.0-7.2) for those heterozygous for missense PVs compared with 7.0 (3.3-14.7) for those harboring PTC variants. BRCA1 p.[Cys64Arg] and BRCA2 p.[Trp2626Cys] were associated with particularly low risks of breast cancer compared with other PVs.

CONCLUSION

These results have important implications for the counseling of at-risk women who harbor missense PVs in the BRCA1/2 genes.

Classification of 101 BRCA1 and BRCA2 variants of uncertain significance by cosegregation study: A powerful approach

Up to 80% of BRCA1 and BRCA2 genetic variants remain of uncertain clinical significance (VUSs). Only variants classified as pathogenic or likely pathogenic can guide breast and ovarian cancer prevention measures and treatment by PARP inhibitors. We report the first results of the ongoing French national COVAR (cosegregation variant) study, the aim of which is to classify BRCA1/2 VUSs. The classification method was a multifactorial model combining different associations between VUSs and cancer, including cosegregation data. At this time, among the 653 variants selected, 101 (15%) distinct variants shared by 1,624 families were classified as pathogenic/likely pathogenic or benign/likely benign by the COVAR study. Sixty-six of the 101 (65%) variants classified by COVAR would have remained VUSs without cosegregation data. Of note, among the 34 variants classified as pathogenic by COVAR, 16 remained VUSs or likely pathogenic when following the ACMG/AMP variant classification guidelines. Although the initiation and organization of cosegregation analyses require a considerable effort, the growing number of available genetic tests results in an increasing number of families sharing a particular variant, and thereby increases the power of such analyses. Here we demonstrate that variant cosegregation analyses are a powerful tool for the classification of variants in the BRCA1/2 breast-ovarian cancer predisposition genes.

5' Region Large Genomic Rearrangements in the BRCA1 Gene in French Families: Identification of a Tandem Triplication and Nine Distinct Deletions with Five Recurrent Breakpoints

Simple Summary

Large genomic rearrangements in BRCA1 consisting of deletions/duplications of one or several exons are complex events, often occurring in the 5′ region. We characterized 10 events in 20 families: one large triplication classified as benign and nine large deletions classified as pathogenic. The breakpoint localization will certainly help to further understand the chromatin structure in regions sensitive to rearrangement.

Abstract

Background: Large genomic rearrangements (LGR) in BRCA1 consisting of deletions/duplications of one or several exons have been found throughout the gene with a large proportion occurring in the 5′ region from the promoter to exon 2. The aim of this study was to better characterize those LGR in French high-risk breast/ovarian cancer families. Methods: DNA from 20 families with one apparent duplication and nine deletions was analyzed with a dedicated comparative genomic hybridization (CGH) array, high-resolution BRCA1 Genomic Morse Codes analysis and Sanger sequencing. Results: The apparent duplication was in fact a tandem triplication of exons 1 and 2 and part of intron 2 of BRCA1, fully characterized here for the first time. We calculated a causality score with the multifactorial model from data obtained from six families, classifying this variant as benign. Among the nine deletions detected in this region, eight have never been identified. The breakpoints fell in six recurrent regions and could confirm some specific conformation of the chromatin. Conclusions: Taken together, our results firmly establish that the BRCA1 5′ region is a frequent site of different LGRs and highlight the importance of the segmental duplication and Alu sequences, particularly the very high homologous region, in the mechanism of a recombination event. This also confirmed that those events are not systematically deleterious.

A Population-Based Analysis of BRCA1/2 Genes and Associated Breast and Ovarian Cancer Risk in Korean Patients: A Multicenter Cohort Study

Simple Summary

Although it has been suggested that cancer risk and genetic variation vary by population, there is still a lack of research on non-European populations. In this study, we applied Korean patients as a model to find out the way to conduct BRCA1/2-related clinical studies in non-European populations who do not have as much clinical data as Europeans. The BRCA1/2 variants were classified following the 2015 ACMG standards/guidelines and using a multifactorial probability-based approach. To estimate the additional sample numbers needed to resolve BRCA1/2 unclassified status, we applied a simulation analysis considering population-specific clinical characteristics. In addition, we estimated the risks of breast or ovarian cancer for BRCA1/2 carriers by mutation regions. Data from this study reveal that BRCA1/2 variants in the non-European population are highly specific; therefore, population-specific study is essential for clinical application of treatment or prevention for breast or ovarian cancer.

Abstract

In this study, we performed a comprehensive analysis of BRCA1/2 variants and associated cancer risk in Korean patients considering two aspects: variants of uncertain significance (VUS) and pathogenic or likely pathogenic variants (PLPVs) in BRCA1 and BRCA2. This study included 5433 Korean participants who were tested for BRCA1/2 genes. The BRCA1/2 variants were classified following the standards/guidelines for interpretation of genetic variants and using a multifactorial probability-based approach. In Korea, 15.8% of participants had BRCA1 or BRCA2 PLPVs. To estimate the additional sample numbers needed to resolve unclassified status, we applied a simulation analysis. The simulation study for VUS showed that the smaller the number of samples, the more the posterior probability was affected by the prior probability; in addition, more samples for BRCA2 VUS than those of BRCA1 VUS were required to resolve the unclassified status, and the presence of clinical information associated with their VUS was an important factor. The cumulative lifetime breast cancer risk was 59.1% (95% CI: 44.1–73.6%) for BRCA1 and 58.3% (95% CI: 43.2–73.0%) for BRCA2 carriers. The cumulative lifetime ovarian cancer risk was estimated to be 36.9% (95% CI: 23.4–53.9%) for BRCA1 and 14.9% (95% CI: 7.4–28.5%) for BRCA2 carriers.

Epidemiological and ES cell-based functional evaluation of BRCA2 variants identified in families with breast cancer

The discovery of high-risk breast cancer susceptibility genes, such as Breast cancer associated gene 1 (BRCA1) and Breast cancer associated gene 2 (BRCA2) has led to accurate identification of individuals for risk management and targeted therapy. The rapid decline in sequencing costs has tremendously increased the number of individuals who are undergoing genetic testing world-wide. However, given the significant differences in population-specific variants, interpreting the results of these tests can be challenging especially for novel genetic variants in understudied populations. Here we report the characterization of novel variants in the Malaysian and Singaporean population that consist of different ethnic groups (Malays, Chinese, Indian, and other indigenous groups). We have evaluated the functional significance of 14 BRCA2 variants of uncertain clinical significance by using multiple in silico prediction tools and examined their frequency in a cohort of 7840 breast cancer cases and 7928 healthy controls. In addition, we have used a mouse embryonic stem cell (mESC)-based functional assay to assess the impact of these variants on BRCA2 function. We found these variants to be functionally indistinguishable from wild-type BRCA2. These variants could fully rescue the lethality of Brca2-null mESCs and exhibited no sensitivity to six different DNA damaging agents including a poly ADP ribose polymerase inhibitor. Our findings strongly suggest that all 14 evaluated variants are functionally neutral. Our findings should be valuable in risk assessment of individuals carrying these variants.

A computational model for classification of BRCA2 variants using mouse embryonic stem cell-based functional assays

Sequencing-based genetic tests to identify individuals at increased risk of hereditary breast and ovarian cancers have resulted in the identification of more than 40,000 sequence variants of BRCA1 and BRCA2. A majority of these variants are considered to be variants of uncertain significance (VUS) because their impact on disease risk remains unknown, largely due to lack of sufficient familial linkage and epidemiological data. Several assays have been developed to examine the effect of VUS on protein function, which can be used to assess their impact on cancer susceptibility. In this study, we report the functional characterization of 88 BRCA2 variants, including several previously uncharacterized variants, using a well-established mouse embryonic stem cell (mESC)-based assay. We have examined their ability to rescue the lethality of Brca2 null mESC as well as sensitivity to six DNA damaging agents including ionizing radiation and a PARP inhibitor. We have also examined the impact of BRCA2 variants on splicing. In addition, we have developed a computational model to determine the probability of impact on function of the variants that can be used for risk assessment. In contrast to the previous VarCall models that are based on a single functional assay, we have developed a new platform to analyze the data from multiple functional assays separately and in combination. We have validated our VarCall models using 12 known pathogenic and 10 neutral variants and demonstrated their usefulness in determining the pathogenicity of BRCA2 variants that are listed as VUS or as variants with conflicting functional interpretation.

Segregation analysis of the BRCA2 c.9227G>T variant in multiple families suggests a pathogenic role in breast and ovarian cancer predisposition

Classification of variants in the BRCA1 and BRCA2 genes has a major impact on the clinical management of subjects at high risk for breast and ovarian cancer. The identification of a pathogenic variant allows for early detection/prevention strategies in healthy carriers as well as targeted treatments in patients affected by BRCA-associated tumors. The BRCA2 c.9227G>T p.(Gly3076Val) variant recurs in families from Northeast Italy and is rarely reported in international databases. This variant substitutes the evolutionary invariant glycine 3076 with a valine in the DNA binding domain of the BRCA2 protein, thus suggesting a high probability of pathogenicity. We analysed clinical and genealogic data of carriers from 15 breast/ovarian cancer families in whom no other pathogenic variants were detected. The variant was shown to co-segregate with breast and ovarian cancer in the most informative families. Combined segregation data led to a likelihood ratio of 81,527:1 of pathogenicity vs. neutrality. We conclude that c.9227G>T is a BRCA2 pathogenic variant that recurs in Northeast Italy. It can now be safely used for the predictive testing of healthy family members to guide preventive surgery and/or early tumor detection strategies, as well as for PARP inhibitors treatments in patients with BRCA2-associated tumors.

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 Categorization of BRCA1 Variants of Uncertain Clinical Significance in Homologous Recombination Repair Complementation Assays

PURPOSE

Because BRCA1 is a high-risk breast/ovarian cancer susceptibility gene, BRCA1 sequence variants of uncertain clinical significance (VUS) complicate genetic counseling. As most VUS are rare, reliable classification based on clinical and genetic data is often impossible. However, all pathogenic BRCA1 variants analyzed result in defective homologous recombination DNA repair (HRR). Thus, BRCA1 VUS may be categorized based on their functional impact on this pathway.

EXPERIMENTAL DESIGN

Two hundred thirty-eight BRCA1 VUS-comprising most BRCA1 VUS known in the Netherlands and Belgium-were tested for their ability to complement Brca1-deficient mouse embryonic stem cells in HRR, using cisplatin and olaparib sensitivity assays and a direct repeat GFP (DR-GFP) HRR assay. Assays were validated using 25 known benign and 25 known pathogenic BRCA1 variants. For assessment of pathogenicity by a multifactorial likelihood analysis method, we collected clinical and genetic data for functionally deleterious VUS and VUS occurring in three or more families.

RESULTS

All three assays showed 100% sensitivity and specificity (95% confidence interval, 83%-100%). Out of 238 VUS, 45 showed functional defects, 26 of which were deleterious in all three assays. For 13 of these 26 variants, we could calculate the probability of pathogenicity using clinical and genetic data, resulting in the identification of 7 (likely) pathogenic variants.

CONCLUSIONS

We have functionally categorized 238 BRCA1 VUS using three different HRR-related assays. Classification based on clinical and genetic data alone for a subset of these variants confirmed the high sensitivity and specificity of our functional assays.

A plugin for the Ensembl Variant Effect Predictor that uses MaxEntScan to predict variant spliceogenicity

Summary

Assessing the pathogenicity of genetic variants can be a complex and challenging task. Spliceogenic variants, which alter mRNA splicing, may yield mature transcripts that encode non-functional protein products, an important predictor of Mendelian disease risk. However, most variant annotation tools do not adequately assess spliceogenicity outside the native splice site and thus the disease-causing potential of variants in other intronic and exonic regions is often overlooked. Here, we present a plugin for the Ensembl Variant Effect Predictor that packages MaxEntScan and extends its functionality to provide splice site predictions using a maximum entropy model. The plugin incorporates a sliding window algorithm to predict splice site loss or gain for any variant that overlaps a transcript feature. We also demonstrate the utility of the plugin by comparing our predictions to two mRNA splicing datasets containing several cancer-susceptibility genes.

Availability and implementation

Source code is freely available under the Apache License, Version 2.0: https://github.com/Ensembl/VEP_plugins.

Supplementary information

Supplementary data are available at Bioinformatics online.

HRness in Breast and Ovarian Cancers

Ovarian and breast cancers are currently defined by the main pathways involved in the tumorigenesis. The majority are carcinomas, originating from epithelial cells that are in constant division and subjected to cyclical variations of the estrogen stimulus during the female hormonal cycle, therefore being vulnerable to DNA damage. A portion of breast and ovarian carcinomas arises in the context of DNA repair defects, in which genetic instability is the backdrop for cancer initiation and progression. For these tumors, DNA repair deficiency is now increasingly recognized as a target for therapeutics. In hereditary breast/ovarian cancers (HBOC), tumors with BRCA1/2 mutations present an impairment of DNA repair by homologous recombination (HR). For many years, BRCA1/2 mutations were only screened on germline DNA, but now they are also searched at the tumor level to personalize treatment. The reason of the inactivation of this pathway remains uncertain for most cases, even in the presence of a HR-deficient signature. Evidence indicates that identifying the mechanism of HR inactivation should improve both genetic counseling and therapeutic response, since they can be useful as new biomarkers of response.

Mutational Screening of BRCA1/2 Genes as a Predictive Factor for Therapeutic Response in Epithelial Ovarian Cancer: A Consensus Guide from the Spanish Society of Pathology (SEAP-IAP) and the Spanish Society of Human Genetics (AEGH)

Germline/somatic BRCA-mutated ovarian carcinomas (OC) are associated to have better response with platinum-based chemotherapy and long-term prognosis than non-BRCA-associated OCs. In addition, these mutations are predictive factors to response to Poly(ADP-ribose) polymerase (PARP) inhibitors. Different positioning papers have addressed the clinical recommendations for BRCA testing in OC. This consensus guide represents a collection of technical recommendations to address the detection of BRCA1/2 mutations in the molecular diagnostic testing strategy for OC. Under the coordination of Spanish Society of Pathology (SEAP-IAP) and the Spanish Society of Human Genetics (AEGH), these recommendations have been developed by pathologists and geneticists taking into account previously published recommendations and their experience in the molecular characterization of these genes. Since the implementation of BRCA testing as a predictive factor can initiate the workflow by testing germline mutations in the blood or by testing both germline and somatic mutations in tumor tissue, distinctive features of both strategies are discussed. Additionally, the recommendations included in this paper provide some references, quality parameters, and genomic tools aimed to standardize and facilitate the clinical genomic diagnosis of OC.

Impact of proactive high-throughput functional assay data on BRCA1 variant interpretation in 3684 patients with breast or ovarian cancer

The clinical utility of BRCA1/2 genotyping was recently extended from the selection of subjects at high risk for hereditary breast and ovary cancer to the identification of candidates for poly (ADP-ribose) polymerase (PARP) inhibitor treatment. This underscores the importance of accurate interpretation of BRCA1/2 genetic variants and of reducing the number of variants of uncertain significance (VUSs). Two recent studies by Findlay et al. and Starita et al. introduced high-throughput functional assays, and proactively analyzed variants in specific regions regardless of whether they had been previously observed. We retrospectively reviewed all BRCA1 and BRCA2 germline genetic test reports from patients with breast or ovarian cancer examined at Asan Medical Center (Seoul, Korea) between September 2011 and December 2018. Variants were assigned pathogenic or benign strong evidence codes according to the functional classification and were reclassified according to the ACMG/AMP 2015 guidelines. Among 3684 patients with available BRCA1 and BRCA2 germline genetic test reports, 429 unique variants (181 from BRCA1) were identified. Of 34 BRCA1 variants intersecting with the data reported by Findlay et al., three missense single-nucleotide variants from four patients (0.11%, 4/3684) were reclassified from VUSs to likely pathogenic variants. Four variants scored as functional were reclassified into benign or likely benign variants. Three variants that overlapped with the data reported by Starita et al. could not be reclassified. In conclusion, proactive high-throughput functional study data are useful for the reclassification of clinically observed VUSs. Integrating additional evidence, including functional assay results, may help reduce the number of VUSs.

Classification of variants of uncertain significance in BRCA1 and BRCA2 using personal and family history of cancer from individuals in a large hereditary cancer multigene panel testing cohort

Purpose

Genetic testing of individuals often results in identification of genomic variants of unknown significance (VUS). Multiple lines of evidence are used to help determine the clinical significance of these variants.

Methods

We analyzed ~138,000 individuals tested by multigene panel testing (MGPT). We used logistic regression to predict carrier status based on personal and family history of cancer. This was applied to 4644 tested individuals carrying 2383 BRCA1/2 variants to calculate likelihood ratios informing pathogenicity for each. Heterogeneity tests were performed for specific classes of variants defined by in silico predictions.

Results

Twenty-two variants labeled as VUS had odds of >10:1 in favor of pathogenicity. The heterogeneity analysis found that among variants in functional domains that were predicted to be benign by in silico tools, a significantly higher proportion of variants were estimated to be pathogenic than previously indicated; that missense variants outside of functional domains should be considered benign; and that variants predicted to create de novo donor sites were also largely benign.

Conclusion

The evidence presented here supports the use of personal and family history from MGPT in the classification of VUS and will be integrated into ongoing efforts to provide large-scale multifactorial classification.

Characterization of splice-altering mutations in inherited predisposition to cancer

Mutations responsible for inherited disease may act by disrupting normal transcriptional splicing. Such mutations can be difficult to detect, and their effects difficult to characterize, because many lie deep within exons or introns where they may alter splice enhancers or silencers or introduce new splice acceptors or donors. Multiple mutation-specific and genome-wide approaches have been developed to evaluate these classes of mutations. We introduce a complementary experimental approach, cBROCA, which yields qualitative and quantitative assessments of the effects of genomic mutations on transcriptional splicing of tumor suppressor genes. cBROCA analysis is undertaken by deriving complementary DNA (cDNA) from puromycin-treated patient lymphoblasts, hybridizing the cDNA to the BROCA panel of tumor suppressor genes, and then multiplex sequencing to very high coverage. At each splice junction suggested by split sequencing reads, read depths of test and control samples are compared. Significant Z scores indicate altered transcripts, over and above naturally occurring minor transcripts, and comparisons of read depths indicate relative abundances of mutant and normal transcripts. BROCA analysis of genomic DNA suggested 120 rare mutations from 150 families with cancers of the breast, ovary, uterus, or colon, in >600 informative genotyped relatives. cBROCA analysis of their transcripts revealed a wide variety of consequences of abnormal splicing in tumor suppressor genes, including whole or partial exon skipping, exonification of intronic sequence, loss or gain of exonic and intronic splicing enhancers and silencers, complete intron retention, hypomorphic alleles, and combinations of these alterations. Combined with pedigree analysis, cBROCA sequencing contributes to understanding the clinical consequences of rare inherited mutations.

Large scale multifactorial likelihood quantitative analysis of BRCA1 and BRCA2 variants: An ENIGMA resource to support clinical variant classification

The multifactorial likelihood analysis method has demonstrated utility for quantitative assessment of variant pathogenicity for multiple cancer syndrome genes. Independent data types currently incorporated in the model for assessing BRCA1 and BRCA2 variants include clinically calibrated prior probability of pathogenicity based on variant location and bioinformatic prediction of variant effect, co-segregation, family cancer history profile, co-occurrence with a pathogenic variant in the same gene, breast tumor pathology, and case-control information. Research and clinical data for multifactorial likelihood analysis were collated for 1,395 BRCA1/2 predominantly intronic and missense variants, enabling classification based on posterior probability of pathogenicity for 734 variants: 447 variants were classified as (likely) benign, and 94 as (likely) pathogenic; and 248 classifications were new or considerably altered relative to ClinVar submissions. Classifications were compared with information not yet included in the likelihood model, and evidence strengths aligned to those recommended for ACMG/AMP classification codes. Altered mRNA splicing or function relative to known nonpathogenic variant controls were moderately to strongly predictive of variant pathogenicity. Variant absence in population datasets provided supporting evidence for variant pathogenicity. These findings have direct relevance for BRCA1 and BRCA2 variant evaluation, and justify the need for gene-specific calibration of evidence types used for variant classification.

Assessment of blind predictions of the clinical significance of BRCA1 and BRCA2 variants

Testing for variation in BRCA1 and BRCA2 (commonly referred to as BRCA1/2), has emerged as a standard clinical practice and is helping countless women better understand and manage their heritable risk of breast and ovarian cancer. Yet the increased rate of BRCA1/2 testing has led to an increasing number of Variants of Uncertain Significance (VUS), and the rate of VUS discovery currently outpaces the rate of clinical variant interpretation. Computational prediction is a key component of the variant interpretation pipeline. In the CAGI5 ENIGMA Challenge, six prediction teams submitted predictions on 326 newly-interpreted variants from the ENIGMA Consortium. By evaluating these predictions against the new interpretations, we have gained a number of insights on the state of the art of variant prediction and specific steps to further advance this state of the art.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5′ and 3′ consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Incorporation of semi-quantitative analysis of splicing alterations for the clinical interpretation of variants in BRCA1 and BRCA2 genes

BRCA1 and BRCA2 (BRCA1/2) genetic variants that disrupt messenger RNA splicing are commonly associated with increased risks of developing breast/ovarian cancer. The majority of splicing studies published to date rely on qualitative methodologies (i.e., Sanger sequencing), but it is necessary to incorporate semi-quantitative or quantitative approaches to accurately interpret the clinical significance of spliceogenic variants. Here, we characterize the splicing impact of 31 BRCA1/2 variants using semi-quantitative capillary electrophoresis of fluorescent amplicons (CE), Sanger sequencing and allele-specific assays. A total of 14 variants were found to disrupt splicing. Allelic-specific assays could be performed for BRCA1 c.302-1G>A and BRCA2 c.516+2T>A, c.1909+1G>A, c.8332-13T>G, c.8332-2A>G, c.8954-2A>T variants, showing a monoallelic contribution to full-length transcript expression that was concordant with semi-quantitative data. The splicing fraction of alternative and aberrant transcripts was also measured by CE, facilitating variant interpretation. Following Evidence-based Network for the Interpretation of Germline Mutant Alleles criteria, we successfully classified eight variants as pathogenic (Class 5), five variants as likely pathogenic (Class 4), and 14 variants as benign (Class 1). We also provide splicing data for four variants classified as uncertain (Class 3), which produced a "leaky" splicing effect or introduced a missense change in the protein sequence, that will require further assessment to determine their clinical significance.

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

Towards controlled terminology for reporting germline cancer susceptibility variants: an ENIGMA report

The vocabulary currently used to describe genetic variants and their consequences reflects many years of studying and discovering monogenic disease with high penetrance. With the recent rapid expansion of genetic testing brought about by wide availability of high-throughput massively parallel sequencing platforms, accurate variant interpretation has become a major issue. The vocabulary used to describe single genetic variants in silico, in vitro, in vivo and as a contributor to human disease uses terms in common, but the meaning is not necessarily shared across all these contexts. In the setting of cancer genetic tests, the added dimension of using data from genetic sequencing of tumour DNA to direct treatment is an additional source of confusion to those who are not experienced in cancer genetics. The language used to describe variants identified in cancer susceptibility genetic testing typically still reflects an outdated paradigm of Mendelian inheritance with dichotomous outcomes. Cancer is a common disease with complex genetic architecture; an improved lexicon is required to better communicate among scientists, clinicians and patients, the risks and implications of genetic variants detected. This review arises from a recognition of, and discussion about, inconsistencies in vocabulary usage by members of the ENIGMA international multidisciplinary consortium focused on variant classification in breast-ovarian cancer susceptibility genes. It sets out the vocabulary commonly used in genetic variant interpretation and reporting, and suggests a framework for a common vocabulary that may facilitate understanding and clarity in clinical reporting of germline genetic tests for cancer susceptibility.

Hereditary Cancer Syndromes-A Primer on Diagnosis and Management: Part 1: Breast-Ovarian Cancer Syndromes

Cancer is the second leading cause of death in both men and women in the United States, with colorectal cancer and breast cancer being two of the most frequent cancer types. Hereditary causes occurring due to pathogenic sequence variants and defects in certain genes makes up roughly 5% of all colorectal cancers and breast-ovarian cancers. High-risk hereditary predisposition syndromes have been associated with a substantially increased lifetime risk for the development of colorectal cancers and breast-ovarian cancers depending on the genetic syndrome, and many people also carry an increased risk of several other cancers compared with the general population. The aim of this review was to provide comprehensive literature on the most commonly encountered hereditary predisposition syndromes, including Lynch syndrome, familial adenomatous polyposis, MUTYH-associated polyposis, hamartomatous polyposis, and breast-ovarian cancer conditions. This will be presented as a 2-part series: the first part will cover the breast-ovarian cancer syndromes, and the second will focus on the inherited colorectal cancer and polyposis conditions.

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.

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

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

Genetic predisposition to colorectal cancer: syndromes, genes, classification of genetic variants and implications for precision medicine

This article reviews genes and syndromes associated with predisposition to colorectal cancer (CRC), with an overview of gene variant classification. We include updates on the application of preventive and therapeutic measures, focusing on the use of non-steroidal anti-inflammatory drugs (NSAIDs) and immunotherapy. Germline pathogenic variants in genes conferring high or moderate risk to cancer are detected in 6-10% of all CRCs and 20% of those diagnosed before age 50. CRC syndromes can be subdivided into nonpolyposis and polyposis entities, the most common of which are Lynch syndrome and familial adenomatous polyposis, respectively. In addition to known and novel genes associated with highly penetrant CRC risk, identification of pathogenic germline variants in genes associated with moderate-penetrance cancer risk and/or hereditary cancer syndromes not traditionally linked to CRC may have an impact on genetic testing, counseling, and surveillance. The use of multigene panels in genetic testing has exposed challenges in the classification of variants of uncertain significance. We provide an overview of the main classification systems and strategies for improving these. Finally, we highlight approaches for integrating chemoprevention in the care of individuals with genetic predisposition to CRC and use of targeted agents and immunotherapy for treatment of mismatch repair-deficient and hypermutant tumors. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

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

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

BRCA Challenge: BRCA Exchange as a global resource for variants in BRCA1 and BRCA2

The BRCA Challenge is a long-term data-sharing project initiated within the Global Alliance for Genomics and Health (GA4GH) to aggregate BRCA1 and BRCA2 data to support highly collaborative research activities. Its goal is to generate an informed and current understanding of the impact of genetic variation on cancer risk across the iconic cancer predisposition genes, BRCA1 and BRCA2. Initially, reported variants in BRCA1 and BRCA2 available from public databases were integrated into a single, newly created site, www.brcaexchange.org. The purpose of the BRCA Exchange is to provide the community with a reliable and easily accessible record of variants interpreted for a high-penetrance phenotype. More than 20,000 variants have been aggregated, three times the number found in the next-largest public database at the project’s outset, of which approximately 7,250 have expert classifications. The data set is based on shared information from existing clinical databases—Breast Cancer Information Core (BIC), ClinVar, and the Leiden Open Variation Database (LOVD)—as well as population databases, all linked to a single point of access. The BRCA Challenge has brought together the existing international Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium expert panel, along with expert clinicians, diagnosticians, researchers, and database providers, all with a common goal of advancing our understanding of BRCA1 and BRCA2 variation. Ongoing work includes direct contact with national centers with access to BRCA1 and BRCA2 diagnostic data to encourage data sharing, development of methods suitable for extraction of genetic variation at the level of individual laboratory reports, and engagement with participant communities to enable a more comprehensive understanding of the clinical significance of genetic variation in BRCA1 and BRCA2.

The goal of this study and paper has been to develop an international resource to generate an informed and current understanding of the impact of genetic variation on cancer risk across the cancer predisposition genes, BRCA1 and BRCA2. Reported variants in BRCA1 and BRCA2 available from public databases were integrated into a single, newly created site, www.brcaexchange.org, to provide a reliable and easily accessible record of variants interpreted for a high-penetrance phenotype.

BRCA1 and BRCA2 5' noncoding region variants identified in breast cancer patients alter promoter activity and protein binding

The widespread use of next generation sequencing for clinical testing is detecting an escalating number of variants in noncoding regions of the genome. The clinical significance of the majority of these variants is currently unknown, which presents a significant clinical challenge. We have screened over 6,000 early-onset and/or familial breast cancer (BC) cases collected by the ENIGMA consortium for sequence variants in the 5' noncoding regions of BC susceptibility genes BRCA1 and BRCA2, and identified 141 rare variants with global minor allele frequency < 0.01, 76 of which have not been reported previously. Bioinformatic analysis identified a set of 21 variants most likely to impact transcriptional regulation, and luciferase reporter assays detected altered promoter activity for four of these variants. Electrophoretic mobility shift assays demonstrated that three of these altered the binding of proteins to the respective BRCA1 or BRCA2 promoter regions, including NFYA binding to BRCA1:c.-287C>T and PAX5 binding to BRCA2:c.-296C>T. Clinical classification of variants affecting promoter activity, using existing prediction models, found no evidence to suggest that these variants confer a high risk of disease. Further studies are required to determine if such variation may be associated with a moderate or low risk of BC.

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.

Improved, ACMG-compliant, in silico prediction of pathogenicity for missense substitutions encoded by TP53 variants

Clinical interpretation of germline missense variants represents a major challenge, including those in the TP53 Li-Fraumeni syndrome gene. Bioinformatic prediction is a key part of variant classification strategies. We aimed to optimize the performance of the Align-GVGD tool used for p53 missense variant prediction, and compare its performance to other bioinformatic tools (SIFT, PolyPhen-2) and ensemble methods (REVEL, BayesDel). Reference sets of assumed pathogenic and assumed benign variants were defined using functional and/or clinical data. Area under the curve and Matthews correlation coefficient (MCC) values were used as objective functions to select an optimized protein multisequence alignment with best performance for Align-GVGD. MCC comparison of tools using binary categories showed optimized Align-GVGD (C15 cut-off) combined with BayesDel (0.16 cut-off), or with REVEL (0.5 cut-off), to have the best overall performance. Further, a semi-quantitative approach using multiple tiers of bioinformatic prediction, validated using an independent set of nonfunctional and functional variants, supported use of Align-GVGD and BayesDel prediction for different strength of evidence levels in ACMG/AMP rules. We provide rationale for bioinformatic tool selection for TP53 variant classification, and have also computed relevant bioinformatic predictions for every possible p53 missense variant to facilitate their use by the scientific and medical community.

Using Somatic Mutations from Tumors to Classify Variants in Mismatch Repair Genes

Present guidelines for classification of constitutional variants do not incorporate inferences from mutations seen in tumors, even when these are associated with a specific molecular phenotype. When somatic mutations and constitutional mutations lead to the same molecular phenotype, as for the mismatch repair genes, information from somatic mutations may enable interpretation of previously unclassified variants. To test this idea, we first estimated likelihoods that somatic variants in MLH1, MSH2, MSH6, and PMS2 drive microsatellite instability and characteristic IHC staining patterns by calculating likelihoods of high versus low normalized variant read fractions of 153 mutations known to be pathogenic versus those of 760 intronic passenger mutations from 174 paired tumor-normal samples. Mutations that explained the tumor mismatch repair phenotype had likelihood ratio for high variant read fraction of 1.56 (95% CI 1.42-1.71) at sites with no loss of heterozygosity and of 26.5 (95% CI 13.2-53.0) at sites with loss of heterozygosity. Next, we applied these ratios to 165 missense, synonymous, and splice variants observed in tumors, combining in a Bayesian analysis the likelihood ratio corresponding with the adjusted variant read fraction with pretest probabilities derived from published analyses and public databases. We suggest classifications for 86 of 165 variants: 7 benign, 31 likely benign, 22 likely pathogenic, and 26 pathogenic. These results illustrate that for mismatch repair genes, characterization of tumor mutations permits tumor mutation data to inform constitutional variant classification. We suggest modifications to incorporate molecular phenotype in future variant classification guidelines.

Pancreatic cancer as a sentinel for hereditary cancer predisposition

BACKGROUND

Genes associated with hereditary breast and ovarian cancer (HBOC) and colorectal cancer (CRC) predisposition have been shown to play a role in pancreatic cancer susceptibility. Growing evidence suggests that pancreatic cancer may be useful as a sentinel cancer to identify families that could benefit from HBOC or CRC surveillance, but to date pancreatic cancer is only considered an indication for genetic testing in the context of additional family history.

METHODS

Preliminary data generated at the Huntsman Cancer Hospital (HCH) included variants identified on a custom 34-gene panel or 59-gene panel including both known HBOC and CRC genes for respective sets of 66 and 147 pancreatic cancer cases, unselected for family history. Given the strength of preliminary data and corresponding literature, 61 sequential pancreatic cancer cases underwent a custom 14-gene clinical panel. Sequencing data from HCH pancreatic cancer cases, pancreatic cancer cases of the Cancer Genome Atlas (TCGA), and an unselected pancreatic cancer screen from the Mayo Clinic were combined in a meta-analysis to estimate the proportion of carriers with pathogenic and high probability of pathogenic variants of uncertain significance (HiP-VUS).

RESULTS

Approximately 8.6% of unselected pancreatic cancer cases at the HCH carried a variant with potential HBOC or CRC screening recommendations. A meta-analysis of unselected pancreatic cancer cases revealed that approximately 11.5% carry a pathogenic variant or HiP-VUS.

CONCLUSION

With the inclusion of both HBOC and CRC susceptibility genes in a panel test, unselected pancreatic cancer cases act as a useful sentinel cancer to identify asymptomatic at-risk relatives who could benefit from relevant HBOC and CRC surveillance measures.

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.

Full in-frame exon 3 skipping of BRCA2 confers high risk of breast and/or ovarian cancer

Germline pathogenic variants in the BRCA2 gene are associated with a cumulative high risk of breast/ovarian cancer. Several BRCA2 variants result in complete loss of the exon-3 at the transcript level. The pathogenicity of these variants and the functional impact of loss of exon 3 have yet to be established.

As a collaboration of the COVAR clinical trial group (France), and the ENIGMA consortium for investigating breast cancer gene variants, this study evaluated 8 BRCA2 variants resulting in complete deletion of exon 3. Clinical information for 39 families was gathered from Portugal, France, Denmark and Sweden. Multifactorial likelihood analyses were conducted using information from 293 patients, for 7 out of the 8 variants (including 6 intronic). For all variants combined the likelihood ratio in favor of causality was 4.39*10²⁵. These results provide convincing evidence for the pathogenicity of all examined variants that lead to a total exon 3 skipping, and suggest that other variants that result in complete loss of exon 3 at the molecular level could be associated with a high risk of cancer comparable to that associated with classical pathogenic variants in BRCA1 or BRCA2 gene. In addition, our functional study shows, for the first time, that deletion of exon 3 impairs the ability of cells to survive upon Mitomycin-C treatment, supporting lack of function for the altered BRCA2 protein in these cells.

Finally, this study demonstrates that any variant leading to expression of only BRCA2 delta-exon 3 will be associated with an increased risk of breast and ovarian cancer.

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

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

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.

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

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

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

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

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

Hereditary breast and ovarian cancer: new genes in confined pathways

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

Implementation of next-generation sequencing for molecular diagnosis of hereditary breast and ovarian cancer highlights its genetic heterogeneity

Molecular diagnosis of hereditary breast and ovarian cancer (HBOC) by standard methodologies has been limited to the BRCA1 and BRCA2 genes. With the recent development of new sequencing methodologies, the speed and efficiency of DNA testing have dramatically improved. The aim of this work was to validate the use of next-generation sequencing (NGS) for the detection of BRCA1/BRCA2 point mutations in a diagnostic setting and to study the role of other genes associated with HBOC in Portuguese families. A cohort of 94 high-risk families was included in the study, and they were initially screened for the two common founder mutations with variant-specific methods. Fourteen index patients were shown to carry the Portuguese founder mutation BRCA2 c.156_157insAlu, and the remaining 80 were analyzed in parallel by Sanger sequencing for the BRCA1/BRCA2 genes and by NGS for a panel of 17 genes that have been described as involved in predisposition to breast and/or ovarian cancer. A total of 506 variants in the BRCA1/BRCA2 genes were detected by both methodologies, with a 100 % concordance between them. This strategy allowed the detection of a total of 39 deleterious mutations in the 94 index patients, namely 10 in BRCA1 (25.6 %), 21 in BRCA2 (53.8 %), four in PALB2 (10.3 %), two in ATM (5.1 %), one in CHEK2 (2.6 %), and one in TP53 (2.6 %), with 20.5 % of the deleterious mutations being found in genes other than BRCA1/BRCA2. These results demonstrate the efficiency of NGS for the detection of BRCA1/BRCA2 point mutations and highlight the genetic heterogeneity of HBOC.