Multifactorial likelihood assessment of BRCA1 and BRCA2 missense variants confirms that BRCA1:c.122A>G(p.His41Arg) is a pathogenic mutation

Neo-Darwinian Principles Exemplified in Cancer Genomics

Within the last two decades, the advent of next-generation sequencing accompanied by single-cell technologies has enabled cancer researchers to study in detail mutations and other genetic aberrations that transpire during transformation of cells to a neoplastic state. This article covers the insights gained through these extensive studies where neo-Darwinian principles can be inferred to play roles throughout neoplastic transformation. The cells promoted during cancer development exhibit cancer hallmarks combined with the related enabling characteristics as outlined by Hanahan and Weinberg, analogous to natural selection and survival of the fittest. Selection of driver mutations that inactivate proteins encoded by tumor suppressor genes differs in profound ways from mutations that activate tumor promoter proteins. In most cases, the later stages of cancer development are characterized by sudden, extensive damage to chromosomes in a process that is not Darwinian in nature. Nevertheless, cells that survive these cataclysmic events remain subject to Darwinian selection promoting clones exhibiting the greatest rates of progression. Duplications of chromosomal segments containing oncogenes, deletions of segments harboring tumor suppressor genes, or distinctive chromosomal rearrangements are often found in cells progressing into later stages of cancer. In summary, the technological developments in genome sequencing since the start of the century have given us clear insights into genomic alterations promoting tumor progression where neo-Darwinian mechanisms of clonal selection can be inferred to play a primary role.

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.

Contribution of mRNA Splicing to Mismatch Repair Gene Sequence Variant Interpretation

Functional assays that assess mRNA splicing can be used in interpretation of the clinical significance of sequence variants, including the Lynch syndrome-associated mismatch repair (MMR) genes. The purpose of this study was to investigate the contribution of splicing assay data to the classification of MMR gene sequence variants. We assayed mRNA splicing for 24 sequence variants in MLH1, MSH2, and MSH6, including 12 missense variants that were also assessed using a cell-free in vitro MMR activity (CIMRA) assay. Multifactorial likelihood analysis was conducted for each variant, combining CIMRA outputs and clinical data where available. We collated these results with existing public data to provide a dataset of splicing assay results for a total of 671 MMR gene sequence variants (328 missense/in-frame indel), and published and unpublished repair activity measurements for 154 of these variants. There were 241 variants for which a splicing aberration was detected: 92 complete impact, 33 incomplete impact, and 116 where it was not possible to determine complete versus incomplete splicing impact. Splicing results mostly aided in the interpretation of intronic (72%) and silent (92%) variants and were the least useful for missense substitutions/in-frame indels (10%). MMR protein functional activity assays were more useful in the analysis of these exonic variants but by design they were not able to detect clinically important splicing aberrations identified by parallel mRNA assays. The development of high throughput assays that can quantitatively assess impact on mRNA transcript expression and protein function in parallel will streamline classification of MMR gene sequence variants.

Systematic misclassification of missense variants in BRCA1 and BRCA2 "coldspots"

Purpose

Guidelines for variant interpretation incorporate variant hotspots in critical functional domains as evidence for pathogenicity (e.g., PM1 and PP2), but do not use “coldspots,” that is, regions without essential functions that tolerate variation, as evidence a variant is benign. To improve variant classification we evaluated BRCA1 and BRCA2 missense variants reported in ClinVar to identify regions where pathogenic missenses are extremely infrequent, defined as coldspots.

Methods

We used Bayesian approaches to model variant classification in these regions.

Results

BRCA1 exon 11 (~60% of the coding sequence), and BRCA2 exons 10 and 11 (~65% of the coding sequence), are coldspots. Of 89 pathogenic (P) or likely pathogenic (LP) missense variants in BRCA1, none are in exon 11 (odds <0.01, 95% confidence interval [CI] 0.0–0.01). Of 34 P or LP missense variants in BRCA2, none are in exons 10–11 (odds <0.01, 95% CI 0.0–0.01). More than half of reported missense variants of uncertain significance (VUS) in BRCA1 and BRCA2 are in coldspots (3115/5301 = 58.8%). Reclassifying these 3115 VUS as likely benign would substantially improve variant classification.

Conclusion

In BRCA1 and BRCA2 coldspots, missense variants are very unlikely to be pathogenic. Classification schemes that incorporate coldspots can reduce the number of VUS and mitigate risks from reporting benign variation as VUS.

Whole-genome sequencing reveals clinically relevant insights into the aetiology of familial breast cancers

BACKGROUND

Whole-genome sequencing (WGS) is a powerful method for revealing the diversity and complexity of the somatic mutation burden of tumours. Here, we investigated the utility of tumour and matched germline WGS for understanding aetiology and treatment opportunities for high-risk individuals with familial breast cancer.

PATIENTS AND METHODS

We carried out WGS on 78 paired germline and tumour DNA samples from individuals carrying pathogenic variants in BRCA1 (n = 26) or BRCA2 (n = 22) or from non-carriers (non-BRCA1/2; n = 30).

RESULTS

Matched germline/tumour WGS and somatic mutational signature analysis revealed patients with unreported, dual pathogenic germline variants in cancer risk genes (BRCA1/BRCA2; BRCA1/MUTYH). The strategy identified that 100% of tumours from BRCA1 carriers and 91% of tumours from BRCA2 carriers exhibited biallelic inactivation of the respective gene, together with somatic mutational signatures suggestive of a functional deficiency in homologous recombination. A set of non-BRCA1/2 tumours also had somatic signatures indicative of BRCA-deficiency, including tumours with BRCA1 promoter methylation, and tumours from carriers of a PALB2 pathogenic germline variant and a BRCA2 variant of uncertain significance. A subset of 13 non-BRCA1/2 tumours from early onset cases were BRCA-proficient, yet displayed complex clustered structural rearrangements associated with the amplification of oncogenes and pathogenic germline variants in TP53, ATM and CHEK2.

CONCLUSIONS

Our study highlights the role that WGS of matched germline/tumour DNA and the somatic mutational signatures can play in the discovery of pathogenic germline variants and for providing supporting evidence for variant pathogenicity. WGS-derived signatures were more robust than germline status and other genomic predictors of homologous recombination deficiency, thus impacting the selection of platinum-based or PARP inhibitor therapy. In this first examination of non-BRCA1/2 tumours by WGS, we illustrate the considerable heterogeneity of these tumour genomes and highlight that complex genomic rearrangements may drive tumourigenesis in a subset of cases.

GFP-Fragment Reassembly Screens for the Functional Characterization of Variants of Uncertain Significance in Protein Interaction Domains of the BRCA1 and BRCA2 Genes

Genetic testing for BRCA1 and BRCA2 genes has led to the identification of many unique variants of uncertain significance (VUS). Multifactorial likelihood models that predict the odds ratio for VUS in favor or against cancer causality, have been developed, but their use is conditioned by the amount of necessary data, which are difficult to obtain if a variant is rare. As an alternative, variants mapping to the coding regions can be examined using in vitro functional assays. BRCA1 and BRCA2 proteins promote genome protection by interacting with different proteins. In this study, we assessed the functional effect of two sets of variants in BRCA genes by exploiting the green fluorescent protein (GFP)-reassembly in vitro assay, which was set-up to test the BRCA1/BARD1, BRCA1/UbcH5a, and BRCA2/DSS1 interactions. Based on the findings observed for the validation panels of previously classified variants, BRCA1/UbcH5a and BRCA2/DSS1 binding assays showed 100% sensitivity and specificity in identifying pathogenic and non-pathogenic variants. While the actual efficiency of these assays in assessing the clinical significance of BRCA VUS has to be verified using larger validation panels, our results suggest that the GFP-reassembly assay is a robust method to identify variants affecting normal protein functioning and contributes to the classification of VUS.

Reclassification of BRCA1 and BRCA2 variants of uncertain significance: a multifactorial analysis of multicentre prospective cohort

BACKGROUND

BRCA1 and BRCA2 (BRCA1/2) variants classified ambiguously as variants of uncertain significance (VUS) are a major challenge for clinical genetic testing in breast cancer; their relevance to the cancer risk is unclear and the association with the response to specific BRCA1/2-targeted agents is uncertain. To minimise the proportion of VUS in BRCA1/2, we performed the multifactorial likelihood analysis and validated this method using an independent cohort of patients with breast cancer.

METHODS

We used a data set of 2115 patients with breast cancer from the nationwide multicentre prospective Korean Hereditary Breast Cancer study. In total, 83 BRCA1/2 VUSs (BRCA1, n=26; BRCA2, n=57) were analysed. The multifactorial probability was estimated by combining the prior probability with the overall likelihood ratio derived from co-occurrence of each VUS with pathogenic variants, personal and family history, and tumour characteristics. The classification was compared with the interpretation according to the American College of Medical Genetics and Genomics-Association for Molecular Pathology (ACMG/AMP) guidelines. An external validation was conducted using independent data set of 810 patients.

RESULTS

We were able to redefine 38 VUSs (BRCA1, n=10; BRCA2, n=28). The revised classification was highly correlated with the ACMG/AMP guideline-based interpretation (BRCA1, p for trend=0.015; BRCA2, p=0.001). Our approach reduced the proportion of VUS from 19% (154/810) to 8.9% (72/810) in the retrospective validation data set.

CONCLUSION

The classification in this study would minimise the 'uncertainty' in clinical interpretation, and this validated multifactorial model can be used for the reliable annotation of BRCA1/2 VUSs.

Performance of in silico prediction tools for the classification of rare BRCA1/2 missense variants in clinical diagnostics

Background

The use of next-generation sequencing approaches in clinical diagnostics has led to a tremendous increase in data and a vast number of variants of uncertain significance that require interpretation. Therefore, prediction of the effects of missense mutations using in silico tools has become a frequently used approach. Aim of this study was to assess the reliability of in silico prediction as a basis for clinical decision making in the context of hereditary breast and/or ovarian cancer.

Methods

We tested the performance of four prediction tools (Align-GVGD, SIFT, PolyPhen-2, MutationTaster2) using a set of 236 BRCA1/2 missense variants that had previously been classified by expert committees. However, a major pitfall in the creation of a reliable evaluation set for our purpose is the generally accepted classification of BRCA1/2 missense variants using the multifactorial likelihood model, which is partially based on Align-GVGD results. To overcome this drawback we identified 161 variants whose classification is independent of any previous in silico prediction. In addition to the performance as stand-alone tools we examined the sensitivity, specificity, accuracy and Matthews correlation coefficient (MCC) of combined approaches.

Results

PolyPhen-2 achieved the lowest sensitivity (0.67), specificity (0.67), accuracy (0.67) and MCC (0.39). Align-GVGD achieved the highest values of specificity (0.92), accuracy (0.92) and MCC (0.73), but was outperformed regarding its sensitivity (0.90) by SIFT (1.00) and MutationTaster2 (1.00). All tools suffered from poor specificities, resulting in an unacceptable proportion of false positive results in a clinical setting. This shortcoming could not be bypassed by combination of these tools. In the best case scenario, 138 families would be affected by the misclassification of neutral variants within the cohort of patients of the German Consortium for Hereditary Breast and Ovarian Cancer.

Conclusion

We show that due to low specificities state-of-the-art in silico prediction tools are not suitable to predict pathogenicity of variants of uncertain significance in BRCA1/2. Thus, clinical consequences should never be based solely on in silico forecasts. However, our data suggests that SIFT and MutationTaster2 could be suitable to predict benignity, as both tools did not result in false negative predictions in our analysis.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0353-y) contains supplementary material, which is available to authorized users.

Clinical Variant Classification: A Comparison of Public Databases and a Commercial Testing Laboratory

BACKGROUND

There is a growing move to consult public databases following receipt of a genetic test result from a clinical laboratory; however, the well-documented limitations of these databases call into question how often clinicians will encounter discordant variant classifications that may introduce uncertainty into patient management. Here, we evaluate discordance in BRCA1 and BRCA2 variant classifications between a single commercial testing laboratory and a public database commonly consulted in clinical practice.

MATERIALS AND METHODS

BRCA1 and BRCA2 variant classifications were obtained from ClinVar and compared with the classifications from a reference laboratory. Full concordance and discordance were determined for variants whose ClinVar entries were of the same pathogenicity (pathogenic, benign, or uncertain). Variants with conflicting ClinVar classifications were considered partially concordant if ≥1 of the listed classifications agreed with the reference laboratory classification.

RESULTS

Four thousand two hundred and fifty unique BRCA1 and BRCA2 variants were available for analysis. Overall, 73.2% of classifications were fully concordant and 12.3% were partially concordant. The remaining 14.5% of variants had discordant classifications, most of which had a definitive classification (pathogenic or benign) from the reference laboratory compared with an uncertain classification in ClinVar (14.0%).

CONCLUSION

Here, we show that discrepant classifications between a public database and single reference laboratory potentially account for 26.7% of variants in BRCA1 and BRCA2. The time and expertise required of clinicians to research these discordant classifications call into question the practicality of checking all test results against a database and suggest that discordant classifications should be interpreted with these limitations in mind.

IMPLICATIONS FOR PRACTICE

With the increasing use of clinical genetic testing for hereditary cancer risk, accurate variant classification is vital to ensuring appropriate medical management. There is a growing move to consult public databases following receipt of a genetic test result from a clinical laboratory; however, we show that up to 26.7% of variants in BRCA1 and BRCA2 have discordant classifications between ClinVar and a reference laboratory. The findings presented in this paper serve as a note of caution regarding the utility of database consultation.

DNA methylation profiling to assess pathogenicity of BRCA1 unclassified variants in breast cancer

Germline pathogenic mutations in BRCA1 increase risk of developing breast cancer. Screening for mutations in BRCA1 frequently identifies sequence variants of unknown pathogenicity and recent work has aimed to develop methods for determining pathogenicity. We previously observed that tumor DNA methylation can differentiate BRCA1-mutated from BRCA1-wild type tumors. We hypothesized that we could predict pathogenicity of variants based on DNA methylation profiles of tumors that had arisen in carriers of unclassified variants. We selected 150 FFPE breast tumor DNA samples [47 BRCA1 pathogenic mutation carriers, 65 BRCAx (BRCA1-wild type), 38 BRCA1 test variants] and analyzed a subset (n=54) using the Illumina 450K methylation platform, using the remaining samples for bisulphite pyrosequencing validation. Three validated markers (BACH2, C8orf31, and LOC654342) were combined with sequence bioinformatics in a model to predict pathogenicity of 27 variants (independent test set).  Predictions were compared with standard multifactorial likelihood analysis. Prediction was consistent for c.5194-12G>A (IVS 19-12 G>A) (P>0.99); 13 variants were considered not pathogenic or likely not pathogenic using both approaches. We conclude that tumor DNA methylation data alone has potential to be used in prediction of BRCA1 variant pathogenicity but is not independent of estrogen receptor status and grade, which are used in current multifactorial models to predict pathogenicity.

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.

Comparative analysis of BRCA1 and BRCA2 variants of uncertain significance in patients with breast cancer: a multifactorial probability-based model versus ACMG standards and guidelines for interpreting sequence variants

PURPOSE

To investigate variants of uncertain significance (VUS) in BRCA1 and BRCA2, we assessed the multifactorial posterior probability of VUS in BRCA1 and BRCA2 and compared these analyses with interpretations according to the recently released American College of Medical Genetics and Genomics (ACMG) standards and guidelines.

METHODS

The analysis involved 715 Korean patients with breast cancer. The multifactorial probability of a VUS was analyzed using the prior probability and combined likelihoods of personal and family history, the pathologic profile of the breast cancer, and co-occurrence with pathogenic variants. Results were compared with those obtained according to the ACMG standards/guidelines.

RESULTS

Sixteen VUS from 51 BRCA1 VUS carriers and 28 VUS from 62 BRCA2 VUS carriers were analyzed. There was a slight agreement between the two analyses, with a kappa value of 0.14 (95% confidence interval (CI) = -0.34 to 0.62) for the BRCA1 VUS and a kappa value of 0.17 (95% CI = -0.10 to 0.49) for the BRCA2 VUS.

CONCLUSION

We propose that genetic counseling should be based on the concordant results between these two analyses. When discrepancies are found, those variants are still considered VUS and careful counseling should be provided.Genet Med 18 12, 1250-1257.

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.