Challenges and approaches to calibrating patient phenotype as evidence for cancer gene variant classification under ACMG/AMP guidelines

Since first publication of the American College of Medical Genetics and Genomics/Association for Medical Pathology (ACMG/AMP) variant classification guidelines, additional recommendations for application of certain criteria have been released (https://clinicalgenome.org/docs/), to improve their application in the diagnostic setting. However, none have addressed use of the PS4 and PP4 criteria, capturing patient presentation as evidence towards pathogenicity. Application of PS4 can be done through traditional case-control studies, or "proband counting" within or across clinical testing cohorts. Review of the existing PS4 and PP4 specifications for Hereditary Cancer Gene Variant Curation Expert Panels revealed substantial differences in the approach to defining specifications. Using BRCA1, BRCA2 and TP53 as exemplar genes, we calibrated different methods proposed for applying the "PS4 proband counting" criterion. For each approach, we considered limitations, non-independence with other ACMG/AMP criteria, broader applicability, and variability in results for different datasets. Our findings highlight inherent overlap of proband-counting methods with ACMG/AMP frequency codes, and the importance of calibration to derive dataset-specific code weights that can account for potential between-dataset differences in ascertainment and other factors. Our work emphasizes the advantages and generalizability of logistic regression analysis over simple proband-counting approaches to empirically determine the relative predictive capacity and weight of various personal clinical features in the context of multigene panel testing, for improved variant interpretation. We also provide a general protocol, including instructions for data formatting and a web-server for analysis of personal history parameters, to facilitate dataset-specific calibration analyses required to use such data for germline variant classification.

CAGI, the Critical Assessment of Genome Interpretation, establishes progress and prospects for computational genetic variant interpretation methods

BACKGROUND

The Critical Assessment of Genome Interpretation (CAGI) aims to advance the state-of-the-art for computational prediction of genetic variant impact, particularly where relevant to disease. The five complete editions of the CAGI community experiment comprised 50 challenges, in which participants made blind predictions of phenotypes from genetic data, and these were evaluated by independent assessors.

RESULTS

Performance was particularly strong for clinical pathogenic variants, including some difficult-to-diagnose cases, and extends to interpretation of cancer-related variants. Missense variant interpretation methods were able to estimate biochemical effects with increasing accuracy. Assessment of methods for regulatory variants and complex trait disease risk was less definitive and indicates performance potentially suitable for auxiliary use in the clinic.

CONCLUSIONS

Results show that while current methods are imperfect, they have major utility for research and clinical applications. Emerging methods and increasingly large, robust datasets for training and assessment promise further progress ahead.

Gene-specific ACMG/AMP classification criteria for germline APC variants: Recommendations from the ClinGen InSiGHT Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel

PURPOSE

The Hereditary Colorectal Cancer/Polyposis Variant Curation Expert Panel (VCEP) was established by the International Society for Gastrointestinal Hereditary Tumours and the Clinical Genome Resource, who set out to develop recommendations for the interpretation of germline APC variants underlying Familial Adenomatous Polyposis, the most frequent hereditary polyposis syndrome.

METHODS

Through a rigorous process of database analysis, literature review, and expert elicitation, the APC VCEP derived gene-specific modifications to the ACMG/AMP (American College of Medical Genetics and Genomics and Association for Molecular Pathology) variant classification guidelines and validated such criteria through the pilot classification of 58 variants.

RESULTS

The APC-specific criteria represented gene- and disease-informed specifications, including a quantitative approach to allele frequency thresholds, a stepwise decision tool for truncating variants, and semiquantitative evaluations of experimental and clinical data. Using the APC-specific criteria, 47% (27/58) of pilot variants were reclassified including 14 previous variants of uncertain significance (VUS).

CONCLUSION

The APC-specific ACMG/AMP criteria preserved the classification of well-characterized variants on ClinVar while substantially reducing the number of VUS by 56% (14/25). Moving forward, the APC VCEP will continue to interpret prioritized lists of VUS, the results of which will represent the most authoritative variant classification for widespread clinical use.

ClinGen guidance for use of the PP1/BS4 co-segregation and PP4 phenotype specificity criteria for sequence variant pathogenicity classification

The 2015 American College of Medical Genetics and Genomics and the Association for Molecular Pathology variant classification publication established a standard employed internationally to guide laboratories in variant assessment. Those recommendations included both pathogenic (PP1) and benign (BS4) criteria for evaluating the inheritance patterns of variants, but details of how to apply those criteria at appropriate evidence levels were sparse. Several publications have since attempted to provide additional guidance, but anecdotally, this issue is still challenging. Additionally, it is not clear that those prior efforts fully distinguished disease-gene identification considerations from variant pathogenicity considerations nor did they address autosomal-recessive and X-linked inheritance. Here, we have taken a mixed inductive and deductive approach to this problem using real diseases as examples. We have developed a practical heuristic for genetic co-segregation evidence and have also determined that the specific phenotype criterion (PP4) is inseparably coupled to the co-segregation criterion. We have also determined that negative evidence at one locus constitutes positive evidence for other loci for disorders with locus heterogeneity. Finally, we provide a points-based system for evaluating phenotype and co-segregation as evidence types to support or refute a locus and show how that can be integrated into the Bayesian framework now used for variant classification and consistent with the 2015 guidelines.

Clinical Multigene Panel Testing Identifies Racial and Ethnic Differences in Germline Pathogenic Variants Among Patients With Early-Onset Colorectal Cancer

PURPOSE

The early-onset colorectal cancer (EOCRC) burden differs across racial/ethnic groups, yet the role of germline genetic predisposition in EOCRC disparities remains uncharacterized. We defined the prevalence and spectrum of inherited colorectal cancer (CRC) susceptibility gene variations among patients with EOCRC by race and ethnicity.

PATIENTS AND METHODS

We included individuals diagnosed with a first primary CRC between age 15 and 49 years who identified as Ashkenazi Jewish, Asian, Black, Hispanic, or White and underwent germline genetic testing of 14 CRC susceptibility genes performed by a clinical testing laboratory. Variant comparisons by racial and ethnic groups were evaluated using chi-square tests and multivariable logistic regression adjusted for sex, age, CRC site, and number of primary colorectal tumors.

RESULTS

Among 3,980 patients with EOCRC, a total of 530 germline pathogenic or likely pathogenic variants were identified in 485 individuals (12.2%). By race/ethnicity, 12.7% of Ashkenazim patients, 9.5% of Asian patients, 10.3% of Black patients, 14.0% of Hispanic patients, and 12.4% of White patients carried a germline variant. The prevalence of Lynch syndrome (P = .037), as well as APC, CHEK2, MLH1, monoallelic MUTYH, and PTEN variants, varied by race/ethnicity among patients with EOCRC (all P < .026). Ashkenazim and Hispanic patients had significantly higher odds of presenting with a pathogenic APC variant, which included p.I1307K (odds ratio [OR], 2.67; 95% CI, 1.30 to 5.49; P = .007) and MLH1 variant (OR, 8.69; 95% CI, 2.68 to 28.20; P = .0003), respectively, versus White patients in adjusted models.

CONCLUSION

Germline genetic features differed by race/ethnicity in young patients with CRC, suggesting that current multigene panel tests may not be representative of EOCRC risk in diverse populations. Further study is needed to optimize genes selected for genetic testing in EOCRC via ancestry-specific gene and variant discovery to yield equitable clinical benefits for all patients and to mitigate inequities in disease burden.

Racial/ethnic differences of germline pathogenic variants in cancer susceptibility genes among patients with early-onset colorectal cancer

19

Background: The burden of early-onset colorectal cancer (CRC diagnosed before age 50 years) differs across racial/ethnic groups. Yet the distinct prevalence and spectrum of germline pathogenic variants associated with early-onset CRC by race/ethnicity is uncharacterized. Methods: We identified 3,980 early-onset CRC patients unselected for family CRC history who underwent clinical multi-gene panel testing of 14 CRC susceptibility genes ( MLH1, MSH2, MSH6, PMS2, SMAD4, APC, BMPR1, CHEK2, EPCAM, MUTYH, PTEN, STK11, CDH1, and TP53) performed by a nationwide clinical testing laboratory from 2012 to 2016 and identified as Ashkenazi Jewish, Asian, Black, Hispanic, or White. Multi-gene panel testing was performed by targeted custom capture and sequencing and chromosome rearrangement analysis. Clinical and individual/family histories were ascertained from clinician-completed test requisition forms. Germline variant prevalence and spectrum in early-onset CRC was determined by race/ethnicity. Pathogenic variant comparisons by race/ethnicity were evaluated using chi-square tests and multivariable logistic regression adjusted for sex, age, CRC site and sequence. Results: Among 3,980 early-onset CRC patients, a total of 503 germline genetic variants were identified in 481 individuals (12.1%). By race/ethnicity, 12.7% of Ashkenazi Jewish, 8.6% of Asian, 10.3% of Black, 14.0% of Hispanic, and 12.3% of White individuals presented with a pathogenic variant in a CRC susceptibility gene. Nearly 7% of early-onset CRC patients (n = 268) presented with a germline pathogenic mutation in a mismatch repair gene—with the prevalence of Lynch Syndrome in early-onset CRC differing across racial/ethnic groups ( P= 0.037). The prevalence of germline APC, CHEK2, MLH1, monoallelic MUTYH, and PTEN variants also varied by race/ethnicity (all P< 0.02). Overall, Asian individuals with early-onset CRC had 49% lower odds of carrying a pathogenic germline variant in a CRC susceptibility gene versus White individuals in adjusted models (Odds Ratio [OR] = 0.51, 95% Confidence Interval [CI]: 0.26-0.99, P= 0.046). Ashkenazim and Hispanic individuals had higher odds of presenting with pathogenic APC and MLH1 variants, respectively, compared with White individuals in adjusted models ( APC: OR = 8.69, 95% CI: 2.68-28.21, P= 0.0003; MLH1: OR = 2.67, 95% CI: 1.30-5.49, P= 0.007). Conclusions: Of 3,980 patients with early-onset CRC unselected for family history, the prevalence and spectrum of pathogenic germline variants differed across racial/ethnic populations. Beyond multigene panel testing for all early-onset CRC patients, further study is needed to optimize genes selected for germline genetic testing in early-onset CRC across diverse groups—which would have clinical implications for disease surveillance and intervention.

Calibration of computational tools for missense variant pathogenicity classification and ClinGen recommendations for PP3/BP4 criteria

Recommendations from the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) for interpreting sequence variants specify the use of computational predictors as "supporting" level of evidence for pathogenicity or benignity using criteria PP3 and BP4, respectively. However, score intervals defined by tool developers, and ACMG/AMP recommendations that require the consensus of multiple predictors, lack quantitative support. Previously, we described a probabilistic framework that quantified the strengths of evidence (supporting, moderate, strong, very strong) within ACMG/AMP recommendations. We have extended this framework to computational predictors and introduce a new standard that converts a tool's scores to PP3 and BP4 evidence strengths. Our approach is based on estimating the local positive predictive value and can calibrate any computational tool or other continuous-scale evidence on any variant type. We estimate thresholds (score intervals) corresponding to each strength of evidence for pathogenicity and benignity for thirteen missense variant interpretation tools, using carefully assembled independent data sets. Most tools achieved supporting evidence level for both pathogenic and benign classification using newly established thresholds. Multiple tools reached score thresholds justifying moderate and several reached strong evidence levels. One tool reached very strong evidence level for benign classification on some variants. Based on these findings, we provide recommendations for evidence-based revisions of the PP3 and BP4 ACMG/AMP criteria using individual tools and future assessment of computational methods for clinical interpretation.

A calibrated cell-based functional assay to aid classification of MLH1 DNA mismatch repair gene variants

Functional assays provide important evidence for classifying the disease significance of germline variants in DNA mismatch repair genes. Numerous laboratories, including our own, have developed functional assays to study mismatch repair gene variants. However, previous assays are limited due to the model system employed, the manner of gene expression, or the environment in which function is assessed. Here, we developed a human cell-based approach for testing the function of variants of uncertain significance (VUS) in the MLH1 gene. Using clustered regularly interspaced short palindromic repeats gene editing, we knocked in MLH1 VUS into the endogenous MLH1 loci in human embryonic stem cells. We examined their impact on RNA and protein, including their ability to prevent microsatellite instability and instigate a DNA damage response. A statistical clustering analysis determined the range of functions associated with known pathogenic or benign variants, and linear regression was performed using existing odds in favor of pathogenicity scores for these control variants to calibrate our functional assay results. By converting the functional outputs into a single odds in favor of pathogenicity score, variant classification expert panels can use these results to readily reassess these VUS. Ultimately, this information will guide proper diagnosis and disease management for suspected Lynch syndrome patients.

Inherited Cancer Susceptibility Gene Sequence Variations Among Patients With Appendix Cancer

Importance

Germline sequence variations in APC, BMPR1A, CDH1, CHEK2, EPCAM, MLH1, MSH2, MSH6, MUTYH, PMS2, PTEN, SMAD4, STK11, and TP53 genes are associated with susceptibility to gastrointestinal cancers. As a rare cancer, the evaluation of appendiceal cancer (AC) predisposition has been limited.

Objective

To assess the prevalence and spectrum of inherited cancer susceptibility gene sequence variations in patients with AC and the utility of germline genetic testing for this population.

Design, Setting, and Participants

This cohort study included patients with AC who underwent germline genetic testing of 14 cancer susceptibility genes performed by a clinical testing laboratory between March 1, 2012, and December 31, 2016. Data were analyzed from March to August 2022. Clinical, individual, and family histories were obtained from clinician-completed test requisition forms. Multigene panel testing was performed by targeted custom capture and sequencing and chromosome rearrangement analysis.

Main Outcomes and Measures

The main outcomes were germline sequence variation prevalence and spectrum in patients with AC.

Results

Among the 131 patients with AC in the cohort (90 female [68.7%]), a total of 16 deleterious sequence variations were identified in 15 patients (11.5%). Similarly, when limited to the 74 patients with AC as the first and only primary tumor, a total of 8 patients (10.8%) had at least 1 deleterious sequence variation in a cancer susceptibility gene. Overall, 6 patients (4.6%) had a deleterious sequence variation observed in MUTYH (5 with monoallelic MUTYH and 1 with biallelic MUTYH). All 4 patients with Lynch syndrome (3.1%) had a sequence variation in the MLH1 gene, of whom 3 were aged 50 years or older at AC diagnosis. Five patients (3.8%) had deleterious sequence variations in other cancer predisposition genes (1 with APC [c.3920T>A, p.I1307K], 2 with CHEK2 [c.470T>C, p.I157T], 1 with SMAD4 [c.263 287dup, p.L98IFS*14], and 1 with TP53 [c.524G>A, p.R175H]).

Conclusions and Relevance

In this cohort study, 1 in every 10 patients with AC who underwent testing for hereditary cancer predisposition carried an inherited gene sequence variation associated with cancer susceptibility. Given the high frequency and broad spectrum of germline gene sequence variations, these data suggest that genetic evaluation might be warranted for all patients diagnosed with this rare malignant tumor. A systemic sequencing effort for all patients with AC may also identify cancer vulnerabilities to exploit for therapeutic development in a cancer type for which clinical trials are limited.

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.

An updated quantitative model to classify missense variants in the TP53 gene: A novel multifactorial strategy

Multigene panel testing has led to an increase in the number of variants of uncertain significance identified in the TP53 gene, associated with Li-Fraumeni syndrome. We previously developed a quantitative model for predicting the pathogenicity of P53 missense variants based on the combination of calibrated bioinformatic information and somatic to germline ratio. Here, we extended this quantitative model for the classification of P53 predicted missense variants by adding new pieces of evidence (personal and family history parameters, loss-of-function results, population allele frequency, healthy individual status by age 60, and breast tumor pathology). We also annotated which missense variants might have an effect on splicing based on bioinformatic predictions. This updated model plus annotation led to the classification of 805 variants into a clinically relevant class, which correlated well with existing ClinVar classifications, and resolved a large number of conflicting and uncertain classifications. We propose this model as a reliable approach to TP53 germline variant classification and emphasize its use in contributing to optimize TP53-specific ACMG/AMP guidelines.

Mobile element insertions and associated structural variants in longitudinal breast cancer samples

While mobile elements are largely inactive in healthy somatic tissues, increased activity has been found in cancer tissues, with significant variation among different cancer types. In addition to insertion events, mobile elements have also been found to mediate many structural variation events in the genome. Here, to better understand the timing and impact of mobile element insertions and associated structural variants in cancer, we examined their activity in longitudinal samples of four metastatic breast cancer patients. We identified 11 mobile element insertions or associated structural variants and found that the majority of these occurred early in tumor progression. Most of the variants impact intergenic regions; however, we identified a translocation interrupting MAP2K4 involving Alu elements and a deletion in YTHDF2 involving mobile elements that likely inactivate reported tumor suppressor genes. The high variant allele fraction of the translocation, the loss of the other copy of MAP2K4, the recurrent loss-of-function mutations found in this gene in other cancers, and the important function of MAP2K4 indicate that this translocation is potentially a driver mutation. Overall, using a unique longitudinal dataset, we find that most variants are likely passenger mutations in the four patients we examined, but some variants impact tumor progression.

Pathogenic Germline Variants in Cancer Susceptibility Genes in Children and Young Adults With Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is the most common pediatric soft-tissue sarcoma and accounts for 3% of all pediatric cancer. In this study, we investigated germline sequence and structural variation in a broad set of genes in two large, independent RMS cohorts.

MATERIALS AND METHODS

Genome sequencing of the discovery cohort (n = 273) and exome sequencing of the secondary cohort (n = 121) were conducted on germline DNA. Analyses were performed on 130 cancer susceptibility genes (CSG). Pathogenic or likely pathogenic (P/LP) variants were predicted using the American College of Medical Genetics and Genomics (ACMG) criteria. Structural variation and survival analyses were performed on the discovery cohort.

RESULTS

We found that 6.6%-7.7% of patients with RMS harbored P/LP variants in dominant-acting CSG. An additional approximately 1% have structural variants (ATM, CDKN1C) in CSGs. CSG variants did not influence survival, although there was a significant correlation with an earlier age of tumor onset. There was a nonsignificant excess of P/LP variants in dominant inheritance genes in the patients with FOXO1 fusion-negative RMS patients versus the patients with FOXO1 fusion-positive RMS. We identified pathogenic germline variants in CSGs previously (TP53, NF1, DICER1, mismatch repair genes), rarely (BRCA2, CBL, CHEK2, SMARCA4), or never (FGFR4) reported in RMS. Numerous genes (TP53, BRCA2, mismatch repair) were on the ACMG Secondary Findings 2.0 list.

CONCLUSION

In two cohorts of patients with RMS, we identified pathogenic germline variants for which gene-specific therapies and surveillance guidelines may be beneficial. In families with a proband with an RMS-risk P/LP variant, genetic counseling and cascade testing should be considered, especially for ACMG Secondary Findings genes and/or with gene-specific surveillance guidelines.

FANCM c5791C>T stopgain mutation (rs144567652) is a familial colorectal cancer risk factor

Purpose

While familial aggregation of colorectal cancer (CRC) is recognized, the majority of the germline predisposition factors remain unidentified, and many high‐risk CRC pedigrees remain unexplained by known risk variants. Fanconi Anemia genes have been recognized to be associated with cancer risk. Notably, FANCM (OMIM 609644) variants have been reported to confer risk for CRC and breast cancer.

Methods

Exome sequencing of CRC‐affected cousins in a set of 47 independent extended high‐risk CRC pedigrees identified a candidate set of rare, shared variants. Variants were tested for association with risk in 744 Utah CRC cases and 1525 controls, and for segregation with CRC in affected relatives.

Results

A FANCM stopgain variant was observed in two CRC‐affected cousin pairs, each from an independent Utah high‐risk pedigree, and yielded a nonsignificant, but elevated OR = 2.05 in a set of Utah cases and controls. Segregation of the variant to other related CRC‐affected cases was observed in the two extended pedigrees.

Conclusion

A rare stopgain variant in FANCM (rs144567652) that is recognized as a breast cancer predisposition variant, and that has previously been proposed, but not confirmed, as a CRC predisposition variant, is validated here as a risk factor for familial CRC.

Fitting a naturally scaled point system to the ACMG/AMP variant classification guidelines

Recently, we demonstrated that the qualitative American College of Medical Genetics and Genomics/Association for Medical Pathology (ACMG/AMP) guidelines for evaluation of Mendelian disease gene variants are fundamentally compatible with a quantitative Bayesian formulation. Here, we show that the underlying ACMG/AMP "strength of evidence categories" can be abstracted into a point system. These points are proportional to Log(odds), are additive, and produce a system that recapitulates the Bayesian formulation of the ACMG/AMP guidelines. The strengths of this system are its simplicity and that the connection between point values and odds of pathogenicity allows empirical calibration of the strength of evidence for individual data types. Weaknesses include that a narrow range of prior probabilities is locked in and that the Bayesian nature of the system is inapparent. We conclude that a points-based system has the practical attribute of user-friendliness and can be useful so long as the underlying Bayesian principles are acknowledged.

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.

Recommendations for application of the functional evidence PS3/BS3 criterion using the ACMG/AMP sequence variant interpretation framework

BACKGROUND

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) clinical variant interpretation guidelines established criteria for different types of evidence. This includes the strong evidence codes PS3 and BS3 for "well-established" functional assays demonstrating a variant has abnormal or normal gene/protein function, respectively. However, they did not provide detailed guidance on how functional evidence should be evaluated, and differences in the application of the PS3/BS3 codes are a contributor to variant interpretation discordance between laboratories. This recommendation seeks to provide a more structured approach to the assessment of functional assays for variant interpretation and guidance on the use of various levels of strength based on assay validation.

METHODS

The Clinical Genome Resource (ClinGen) Sequence Variant Interpretation (SVI) Working Group used curated functional evidence from ClinGen Variant Curation Expert Panel-developed rule specifications and expert opinions to refine the PS3/BS3 criteria over multiple in-person and virtual meetings. We estimated the odds of pathogenicity for assays using various numbers of variant controls to determine the minimum controls required to reach moderate level evidence. Feedback from the ClinGen Steering Committee and outside experts were incorporated into the recommendations at multiple stages of development.

RESULTS

The SVI Working Group developed recommendations for evaluators regarding the assessment of the clinical validity of functional data and a four-step provisional framework to determine the appropriate strength of evidence that can be applied in clinical variant interpretation. These steps are as follows: (1) define the disease mechanism, (2) evaluate the applicability of general classes of assays used in the field, (3) evaluate the validity of specific instances of assays, and (4) apply evidence to individual variant interpretation. We found that a minimum of 11 total pathogenic and benign variant controls are required to reach moderate-level evidence in the absence of rigorous statistical analysis.

CONCLUSIONS

The recommendations and approach to functional evidence evaluation described here should help clarify the clinical variant interpretation process for functional assays. Further, we hope that these recommendations will help develop productive partnerships with basic scientists who have developed functional assays that are useful for interrogating the function of a variety of genes.

Association of Sulindac and Erlotinib vs Placebo With Colorectal Neoplasia in Familial Adenomatous Polyposis: Secondary Analysis of a Randomized Clinical Trial

Importance

Patients with familial adenomatous polyposis (FAP) are at markedly increased risk for colorectal polyps and cancer. A combination of sulindac and erlotinib led to a 71% reduction in duodenal polyp burden in a phase 2 trial.

Objective

To evaluate effect of sulindac and erlotinib on colorectal adenoma regression in patients with FAP.

Design, Setting, and Participants

Prespecified secondary analysis for colorectal adenoma regression was carried out using data from a double-blind, randomized, placebo-controlled trial, enrolling 92 patients with FAP, conducted from July 2010 to June 2014 in Salt Lake City, Utah.

Interventions

Patients were randomized to sulindac, 150 mg twice daily, and erlotinib, 75 mg daily (n = 46), vs placebo (n = 46) for 6 months.

Main Outcomes and Measurements

The total number of polyps in the intact colorectum, ileal pouch anal anastomosis, or ileo-rectum were recorded at baseline and 6 months. The primary outcomes were change in total colorectal polyp count and percentage change in colorectal polyps, following 6 months of treatment.

Results

Eighty-two randomized patients (mean [SD] age, 40 [13] years; 49 [60%] women) had colorectal polyp count data available for this secondary analysis: 22 with intact colon, 44 with ileal pouch anal anastomosis and 16 with ileo-rectal anastomosis; 41 patients received sulindac/erlotinib and 41 placebo. The total colorectal polyp count was significantly different between the placebo and sulindac-erlotinib group at 6 months in patients with net percentage change of 69.4% in those with an intact colorectum compared with placebo (95% CI, 28.8%-109.2%; P = .009).

Conclusion and Relevance

In this double-blind, placebo-controlled, randomized trial we showed that combination treatment with sulindac and erlotinib compared with placebo resulted in significantly lower colorectal polyp burden after 6 months of treatment. There was a reduction in polyp burden in both those with an entire colorectum and those with only a rectal pouch or rectum.

Trial Registration

clinicaltrials.gov Identifier: NCT01187901.

Assessing the performance of in silico methods for predicting the pathogenicity of variants in the gene CHEK2, among Hispanic females with breast cancer

The availability of disease-specific genomic data is critical for developing new computational methods that predict the pathogenicity of human variants and advance the field of precision medicine. However, the lack of gold standards to properly train and benchmark such methods is one of the greatest challenges in the field. In response to this challenge, the scientific community is invited to participate in the Critical Assessment for Genome Interpretation (CAGI), where unpublished disease variants are available for classification by in silico methods. As part of the CAGI-5 challenge, we evaluated the performance of 18 submissions and three additional methods in predicting the pathogenicity of single nucleotide variants (SNVs) in checkpoint kinase 2 (CHEK2) for cases of breast cancer in Hispanic females. As part of the assessment, the efficacy of the analysis method and the setup of the challenge were also considered. The results indicated that though the challenge could benefit from additional participant data, the combined generalized linear model analysis and odds of pathogenicity analysis provided a framework to evaluate the methods submitted for SNV pathogenicity identification and for comparison to other available methods. The outcome of this challenge and the approaches used can help guide further advancements in identifying SNV-disease relationships.

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.

A quantitative model to predict pathogenicity of missense variants in the TP53 gene

Germline pathogenic variants in the TP53 gene cause Li-Fraumeni syndrome, a condition that predisposes individuals to a wide range of cancer types. Identification of individuals carrying a TP53 pathogenic variant is linked to clinical management decisions, such as the avoidance of radiotherapy and use of high-intensity screening programs. The aim of this study was to develop an evidence-based quantitative model that integrates independent in silico data (Align-GVGD and BayesDel) and somatic to germline ratio (SGR), to assign pathogenicity to every possible missense variant in the TP53 gene. To do this, a likelihood ratio for pathogenicity (LR) was derived from each component calibrated using reference sets of assumed pathogenic and benign missense variants. A posterior probability of pathogenicity was generated by combining LRs, and algorithm outputs were validated using different approaches. A total of 730 TP53 missense variants could be assigned to a clinically interpretable class. The outputs of the model correlated well with existing clinical information, functional data, and ClinVar classifications. In conclusion, these quantitative outputs provide the basis for individualized assessment of cancer risk useful for clinical interpretation. In addition, we propose the value of the novel SGR approach for use within the ACMG/AMP guidelines for variant classification.

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.

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.

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.

Panel sequencing of 264 candidate susceptibility genes and segregation analysis in a cohort of non-BRCA1, non-BRCA2 breast cancer families

PURPOSE

The main aim of this study was to screen epigenetic modifier genes and known breast cancer driver genes for germline mutations in non-BRCA1/2 (BRCAx) breast cancer families in order to identify novel susceptibility genes of moderate-high penetrance.

METHODS

We screened 264 candidate susceptibility genes in 656 index cases from non-BRCA1/2 families. Potentially pathogenic candidate mutations were then genotyped in all available family members for the assessment of co-segregation of the variant with disease in the family in order to estimate the breast cancer risks associated with these mutations. For 11 of the candidate susceptibility genes, we screened an additional 800 non-BRCA1/2 breast cancer cases and 787 controls.

RESULTS

Only two genes, CHD8 and USH2A showed any evidence of an increased risk of breast cancer (RR = 2.40 (95% CI 1.0-7.32) and 2.48 (95% CI 1.11-6.67), respectively).

CONCLUSIONS

We found no convincing evidence that epigenetic modifier and known breast cancer driver genes carry germline mutations that increase breast cancer risk. USH2A is no longer regarded as a breast cancer driver gene and seems an implausible candidate given its association with Usher syndrome. However, somatic mutations in CHD8 have been recently reported, making it an even more promising candidate, but further analysis of CHD8 in very large cohorts of families or case-control studies would be required to determine if it is a moderate-risk breast cancer susceptibility gene.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

ABRAXAS (FAM175A) and Breast Cancer Susceptibility: No Evidence of Association in the Breast Cancer Family Registry

Approximately half of the familial aggregation of breast cancer remains unexplained. This proportion is less for early-onset disease where familial aggregation is greater, suggesting that other susceptibility genes remain to be discovered. The majority of known breast cancer susceptibility genes are involved in the DNA double-strand break repair pathway. ABRAXAS is involved in this pathway and mutations in this gene impair BRCA1 recruitment to DNA damage foci and increase cell sensitivity to ionizing radiation. Moreover, a recurrent germline mutation was reported in Finnish high-risk breast cancer families. To determine if ABRAXAS could be a breast cancer susceptibility gene in other populations, we conducted a population-based case-control mutation screening study of the coding exons and exon/intron boundaries of ABRAXAS in the Breast Cancer Family Registry. In addition to the common variant p.Asp373Asn, sixteen distinct rare variants were identified. Although no significant difference in allele frequencies between cases and controls was observed for the identified variants, two variants, p.Gly39Val and p.Thr141Ile, were shown to diminish phosphorylation of gamma-H2AX in MCF7 human breast adenocarcinoma cells, an important biomarker of DNA double-strand breaks. Overall, likely damaging or neutral variants were evenly represented among cases and controls suggesting that rare variants in ABRAXAS may explain only a small proportion of hereditary breast cancer.

Multigene testing of moderate-risk genes: be mindful of the missense

BACKGROUND

Moderate-risk genes have not been extensively studied, and missense substitutions in them are generally returned to patients as variants of uncertain significance lacking clearly defined risk estimates. The fraction of early-onset breast cancer cases carrying moderate-risk genotypes and quantitative methods for flagging variants for further analysis have not been established.

METHODS

We evaluated rare missense substitutions identified from a mutation screen of ATM, CHEK2, MRE11A, RAD50, NBN, RAD51, RINT1, XRCC2 and BARD1 in 1297 cases of early-onset breast cancer and 1121 controls via scores from Align-Grantham Variation Grantham Deviation (GVGD), combined annotation dependent depletion (CADD), multivariate analysis of protein polymorphism (MAPP) and PolyPhen-2. We also evaluated subjects by polygenotype from 18 breast cancer risk SNPs. From these analyses, we estimated the fraction of cases and controls that reach a breast cancer OR≥2.5 threshold.

RESULTS

Analysis of mutation screening data from the nine genes revealed that 7.5% of cases and 2.4% of controls were carriers of at least one rare variant with an average OR≥2.5. 2.1% of cases and 1.2% of controls had a polygenotype with an average OR≥2.5.

CONCLUSIONS

Among early-onset breast cancer cases, 9.6% had a genotype associated with an increased risk sufficient to affect clinical management recommendations. Over two-thirds of variants conferring this level of risk were rare missense substitutions in moderate-risk genes. Placement in the estimated OR≥2.5 group by at least two of these missense analysis programs should be used to prioritise variants for further study. Panel testing often creates more heat than light; quantitative approaches to variant prioritisation and classification may facilitate more efficient clinical classification of variants.

No evidence that protein truncating variants in BRIP1 are associated with breast cancer risk: implications for gene panel testing

BACKGROUND

BRCA1 interacting protein C-terminal helicase 1 (BRIP1) is one of the Fanconi Anaemia Complementation (FANC) group family of DNA repair proteins. Biallelic mutations in BRIP1 are responsible for FANC group J, and previous studies have also suggested that rare protein truncating variants in BRIP1 are associated with an increased risk of breast cancer. These studies have led to inclusion of BRIP1 on targeted sequencing panels for breast cancer risk prediction.

METHODS

We evaluated a truncating variant, p.Arg798Ter (rs137852986), and 10 missense variants of BRIP1, in 48 144 cases and 43 607 controls of European origin, drawn from 41 studies participating in the Breast Cancer Association Consortium (BCAC). Additionally, we sequenced the coding regions of BRIP1 in 13 213 cases and 5242 controls from the UK, 1313 cases and 1123 controls from three population-based studies as part of the Breast Cancer Family Registry, and 1853 familial cases and 2001 controls from Australia.

RESULTS

The rare truncating allele of rs137852986 was observed in 23 cases and 18 controls in Europeans in BCAC (OR 1.09, 95% CI 0.58 to 2.03, p=0.79). Truncating variants were found in the sequencing studies in 34 cases (0.21%) and 19 controls (0.23%) (combined OR 0.90, 95% CI 0.48 to 1.70, p=0.75).

CONCLUSIONS

These results suggest that truncating variants in BRIP1, and in particular p.Arg798Ter, are not associated with a substantial increase in breast cancer risk. Such observations have important implications for the reporting of results from breast cancer screening panels.

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

Clinical mutation screening of the cancer susceptibility genes BRCA1 and BRCA2 generates many unclassified variants (UVs). Most of these UVs are either rare missense substitutions or nucleotide substitutions near the splice junctions of the protein coding exons. Previously, we developed a quantitative method for evaluation of BRCA gene UVs—the “integrated evaluation”—that combines a sequence analysis‐based prior probability of pathogenicity with patient and/or tumor observational data to arrive at a posterior probability of pathogenicity. One limitation of the sequence analysis‐based prior has been that it evaluates UVs from the perspective of missense substitution severity but not probability to disrupt normal mRNA splicing. Here, we calibrated output from the splice‐site fitness program MaxEntScan to generate spliceogenicity‐based prior probabilities of pathogenicity for BRCA gene variants; these range from 0.97 for variants with high probability to damage a donor or acceptor to 0.02 for exonic variants that do not impact a splice junction and are unlikely to create a de novo donor. We created a database http://priors.hci.utah.edu/PRIORS/ that provides the combined missense substitution severity and spliceogenicity‐based probability of pathogenicity for BRCA gene single‐nucleotide substitutions. We also updated the BRCA gene Ex‐UV LOVD, available at http://hci‐exlovd.hci.utah.edu, with 77 re‐evaluable variants.

Effect of Sulindac and Erlotinib vs Placebo on Duodenal Neoplasia in Familial Adenomatous Polyposis: A Randomized Clinical Trial

IMPORTANCE

Patients with familial adenomatous polyposis (FAP) are at markedly increased risk for duodenal polyps and cancer. Surgical and endoscopic management of duodenal neoplasia is difficult and chemoprevention has not been successful.

OBJECTIVE

To evaluate the effect of a combination of sulindac and erlotinib on duodenal adenoma regression in patients with FAP.

DESIGN, SETTING, AND PARTICIPANTS

Double-blind, randomized, placebo-controlled trial, enrolling 92 participants with FAP, conducted from July 2010 through June 2014 at Huntsman Cancer Institute in Salt Lake City, Utah.

INTERVENTIONS

Participants with FAP were randomized to sulindac (150 mg) twice daily and erlotinib (75 mg) daily (n = 46) vs placebo (n = 46) for 6 months.

MAIN OUTCOMES AND MEASURES

The total number and diameter of polyps in the proximal duodenum were mapped at baseline and 6 months. The primary outcome was change in total polyp burden at 6 months. Polyp burden was calculated as the sum of the diameters of polyps. The secondary outcomes were change in total duodenal polyp count, change in duodenal polyp burden or count stratified by genotype and initial polyp burden, and percentage of change from baseline in duodenal polyp burden.

RESULTS

Ninety-two participants (mean age, 41 years [range, 24-55]; women, 56 [61%]) were randomized when the trial was stopped by the external data and safety monitoring board because the second preplanned interim analysis met the prespecified stopping rule for superiority. Grade 1 and 2 adverse events were more common in the sulindac-erlotinib group, with an acne-like rash observed in 87% of participants receiving treatment and 20% of participants receiving placebo (P < .001). Only 2 participants experienced grade 3 adverse events. [table: see text].

CONCLUSIONS AND RELEVANCE

Among participants with FAP, the use of sulindac and erlotinib compared with placebo resulted in a lower duodenal polyp burden after 6 months. Adverse events may limit the use of these medications at the doses used in this study. Further research is necessary to evaluate these preliminary findings in a larger study population with longer follow-up to determine whether the observed effects will result in improved clinical outcomes.

TRIAL REGISTRATION

clinicaltrials.gov Identifier: NCT 01187901.

BRCA1 Circos: a visualisation resource for functional analysis of missense variants

BACKGROUND

Inactivating germline mutations in the tumour suppressor gene BRCA1 are associated with a significantly increased risk of developing breast and ovarian cancer. A large number (>1500) of unique BRCA1 variants have been identified in the population and can be classified as pathogenic, non-pathogenic or as variants of unknown significance (VUS). Many VUS are rare missense variants leading to single amino acid changes. Their impact on protein function cannot be directly inferred from sequence information, precluding assessment of their pathogenicity. Thus, functional assays are critical to assess the impact of these VUS on protein activity. BRCA1 is a multifunctional protein and different assays have been used to assess the impact of variants on different biochemical activities and biological processes.

METHODS AND RESULTS

To facilitate VUS analysis, we have developed a visualisation resource that compiles and displays functional data on all documented BRCA1 missense variants. BRCA1 Circos is a web-based visualisation tool based on the freely available Circos software package. The BRCA1 Circos web tool (http://research.nhgri.nih.gov/bic/circos/) aggregates data from all published BRCA1 missense variants for functional studies, harmonises their results and presents various functionalities to search and interpret individual-level functional information for each BRCA1 missense variant.

CONCLUSIONS

This research visualisation tool will serve as a quick one-stop publically available reference for all the BRCA1 missense variants that have been functionally assessed. It will facilitate meta-analysis of functional data and improve assessment of pathogenicity of VUS.

Growing recognition of the role for rare missense substitutions in breast cancer susceptibility

Most cancer susceptibility genes function as tumor suppressors; accordingly, the focus of mutation screening in breast cancer families has been to identify protein-truncating mutations. However, it is now clear that, for some breast cancer susceptibility genes, a significant proportion of the burden of disease comes from rare missense substitutions. Among genes that have been extensively evaluated, BRCA1, BRCA2, PALB2 and BRIP1 stand as examples where the majority of mutations lead to protein truncation;TP53 provides a counter example, where the majority of pathogenic variants are missense substitutions. In ATM and CHEK2, missense substitutions are probably equally or more important in terms of their frequency and attributable risk. Therefore, ongoing efforts to identify new susceptibility genes should not ignore missense variation.

Rare mutations in RINT1 predispose carriers to breast and Lynch syndrome-spectrum cancers

Approximately half of the familial aggregation of breast cancer remains unexplained. A multiple-case breast cancer family exome-sequencing study identified three likely pathogenic mutations in RINT1 (NM_021930.4) not present in public sequencing databases: RINT1 c.343C>T (p.Q115X), c.1132_1134del (p.M378del), and c.1207G>T (p.D403Y). On the basis of this finding, a population-based case-control mutation-screening study was conducted that identified 29 carriers of rare (minor allele frequency < 0.5%), likely pathogenic variants: 23 in 1,313 early-onset breast cancer cases and six in 1,123 frequency-matched controls [OR, 3.24; 95% confidence interval (CI), 1.29-8.17; P = 0.013]. RINT1 mutation screening of probands from 798 multiple-case breast cancer families identified four additional carriers of rare genetic variants. Analysis of the incidence of first primary cancers in families of women carrying RINT1 mutations estimated that carriers were at increased risk of Lynch syndrome-spectrum cancers [standardized incidence ratio (SIR), 3.35; 95% CI, 1.7-6.0; P = 0.005], particularly for relatives diagnosed with cancer under the age of 60 years (SIR, 10.9; 95% CI, 4.7-21; P = 0.0003).

SIGNIFICANCE

The work described in this study adds RINT1 to the growing list of genes in which rare sequence variants are associated with intermediate levels of breast cancer risk. Given that RINT1 is also associated with a spectrum of cancers with mismatch repair defects, these findings have clinical applications and raise interesting biological questions.

Rare key functional domain missense substitutions in MRE11A, RAD50, and NBN contribute to breast cancer susceptibility: results from a Breast Cancer Family Registry case-control mutation-screening study

INTRODUCTION

The MRE11A-RAD50-Nibrin (MRN) complex plays several critical roles related to repair of DNA double-strand breaks. Inherited mutations in the three components predispose to genetic instability disorders and the MRN genes have been implicated in breast cancer susceptibility, but the underlying data are not entirely convincing. Here, we address two related questions: (1) are some rare MRN variants intermediate-risk breast cancer susceptibility alleles, and if so (2) do the MRN genes follow a BRCA1/BRCA2 pattern wherein most susceptibility alleles are protein-truncating variants, or do they follow an ATM/CHEK2 pattern wherein half or more of the susceptibility alleles are missense substitutions?

METHODS

Using high-resolution melt curve analysis followed by Sanger sequencing, we mutation screened the coding exons and proximal splice junction regions of the MRN genes in 1,313 early-onset breast cancer cases and 1,123 population controls. Rare variants in the three genes were pooled using bioinformatics methods similar to those previously applied to ATM, BRCA1, BRCA2, and CHEK2, and then assessed by logistic regression.

RESULTS

Re-analysis of our ATM, BRCA1, and BRCA2 mutation screening data revealed that these genes do not harbor pathogenic alleles (other than modest-risk SNPs) with minor allele frequencies>0.1% in Caucasian Americans, African Americans, or East Asians. Limiting our MRN analyses to variants with allele frequencies of <0.1% and combining protein-truncating variants, likely spliceogenic variants, and key functional domain rare missense substitutions, we found significant evidence that the MRN genes are indeed intermediate-risk breast cancer susceptibility genes (odds ratio (OR)=2.88, P=0.0090). Key domain missense substitutions were more frequent than the truncating variants (24 versus 12 observations) and conferred a slightly higher OR (3.07 versus 2.61) with a lower P value (0.029 versus 0.14).

CONCLUSIONS

These data establish that MRE11A, RAD50, and NBN are intermediate-risk breast cancer susceptibility genes. Like ATM and CHEK2, their spectrum of pathogenic variants includes a relatively high proportion of missense substitutions. However, the data neither establish whether variants in each of the three genes are best evaluated under the same analysis model nor achieve clinically actionable classification of individual variants observed in this study.

A unified test of linkage analysis and rare-variant association for analysis of pedigree sequence data

High-throughput sequencing of related individuals has become an important tool for studying human disease. However, owing to technical complexity and lack of available tools, most pedigree-based sequencing studies rely on an ad hoc combination of suboptimal analyses. Here we present pedigree-VAAST (pVAAST), a disease-gene identification tool designed for high-throughput sequence data in pedigrees. pVAAST uses a sequence-based model to perform variant and gene-based linkage analysis. Linkage information is then combined with functional prediction and rare variant case-control association information in a unified statistical framework. pVAAST outperformed linkage and rare-variant association tests in simulations and identified disease-causing genes from whole-genome sequence data in three human pedigrees with dominant, recessive and de novo inheritance patterns. The approach is robust to incomplete penetrance and locus heterogeneity and is applicable to a wide variety of genetic traits. pVAAST maintains high power across studies of monogenic, high-penetrance phenotypes in a single pedigree to highly polygenic, common phenotypes involving hundreds of pedigrees.

Application of a 5-tiered scheme for standardized classification of 2,360 unique mismatch repair gene variants in the InSiGHT locus-specific database

The clinical classification of hereditary sequence variants identified in disease-related genes directly affects clinical management of patients and their relatives. The International Society for Gastrointestinal Hereditary Tumours (InSiGHT) undertook a collaborative effort to develop, test and apply a standardized classification scheme to constitutional variants in the Lynch syndrome-associated genes MLH1, MSH2, MSH6 and PMS2. Unpublished data submission was encouraged to assist in variant classification and was recognized through microattribution. The scheme was refined by multidisciplinary expert committee review of the clinical and functional data available for variants, applied to 2,360 sequence alterations, and disseminated online. Assessment using validated criteria altered classifications for 66% of 12,006 database entries. Clinical recommendations based on transparent evaluation are now possible for 1,370 variants that were not obviously protein truncating from nomenclature. This large-scale endeavor will facilitate the consistent management of families suspected to have Lynch syndrome and demonstrates the value of multidisciplinary collaboration in the curation and classification of variants in public locus-specific databases.

BRCA1/2 sequence variants of uncertain significance: a primer for providers to assist in discussions and in medical management

INTRODUCTION

DNA variants of uncertain significance (VUS) are common outcomes of clinical genetic testing for susceptibility to cancer. A statistically rigorous model that provides a pathogenicity score for each variant has been developed to aid in the clinical management of patients undergoing genetic testing.

METHODS

The information in this article is derived from multiple publications on VUS in BRCA genes, distilled for communicating with clinicians who may encounter VUS in their practice.

RESULTS

The posterior probability scores for BRCA1 or BRCA2 VUS, calculated from a multifactorial likelihood model, are explained, and links for looking up specific VUS are provided. The International Agency on Cancer Research (IARC) of the World Health Organization has proposed a simple five-tier system for clinical management that is not widely known to clinicians. Classes 1 and 2 in this system are managed as neutral variants, classes 4 and 5 are managed as pathogenic variants, and class 3 variants still have insufficient evidence to move to either end of this scale and, thus, cannot be used in medical management.

CONCLUSIONS

Development of models that integrate multiple independent lines of evidence has allowed classification of a growing number of VUS in the BRCA1 and BRCA2 genes. The pathogenicity score that is generated by this model maps to the IARC system for clinical management, which will assist clinicians in the medical management of those patients who obtain a VUS result upon testing.

FAVR (Filtering and Annotation of Variants that are Rare): methods to facilitate the analysis of rare germline genetic variants from massively parallel sequencing datasets

Background

Characterising genetic diversity through the analysis of massively parallel sequencing (MPS) data offers enormous potential to significantly improve our understanding of the genetic basis for observed phenotypes, including predisposition to and progression of complex human disease. Great challenges remain in resolving genetic variants that are genuine from the millions of artefactual signals.

Results

FAVR is a suite of new methods designed to work with commonly used MPS analysis pipelines to assist in the resolution of some of the issues related to the analysis of the vast amount of resulting data, with a focus on relatively rare genetic variants. To the best of our knowledge, no equivalent method has previously been described. The most important and novel aspect of FAVR is the use of signatures in comparator sequence alignment files during variant filtering, and annotation of variants potentially shared between individuals. The FAVR methods use these signatures to facilitate filtering of (i) platform and/or mapping-specific artefacts, (ii) common genetic variants, and, where relevant, (iii) artefacts derived from imbalanced paired-end sequencing, as well as annotation of genetic variants based on evidence of co-occurrence in individuals. We applied conventional variant calling applied to whole-exome sequencing datasets, produced using both SOLiD and TruSeq chemistries, with or without downstream processing by FAVR methods. We demonstrate a 3-fold smaller rare single nucleotide variant shortlist with no detected reduction in sensitivity. This analysis included Sanger sequencing of rare variant signals not evident in dbSNP131, assessment of known variant signal preservation, and comparison of observed and expected rare variant numbers across a range of first cousin pairs. The principles described herein were applied in our recent publication identifying XRCC2 as a new breast cancer risk gene and have been made publically available as a suite of software tools.

Conclusions

FAVR is a platform-agnostic suite of methods that significantly enhances the analysis of large volumes of sequencing data for the study of rare genetic variants and their influence on phenotypes.

RAD51 and breast cancer susceptibility: no evidence for rare variant association in the Breast Cancer Family Registry study

Background

Although inherited breast cancer has been associated with germline mutations in genes that are functionally involved in the DNA homologous recombination repair (HRR) pathway, including BRCA1, BRCA2, TP53, ATM, BRIP1, CHEK2 and PALB2, about 70% of breast cancer heritability remains unexplained. Because of their critical functions in maintaining genome integrity and already well-established associations with breast cancer susceptibility, it is likely that additional genes involved in the HRR pathway harbor sequence variants associated with increased risk of breast cancer. RAD51 plays a central biological function in DNA repair and despite the fact that rare, likely dysfunctional variants in three of its five paralogs, RAD51C, RAD51D, and XRCC2, have been associated with breast and/or ovarian cancer risk, no population-based case-control mutation screening data are available for the RAD51 gene. We thus postulated that RAD51 could harbor rare germline mutations that confer increased risk of breast cancer.

Methodology/Principal Findings

We screened the coding exons and proximal splice junction regions of the gene for germline sequence variation in 1,330 early-onset breast cancer cases and 1,123 controls from the Breast Cancer Family Registry, using the same population-based sampling and analytical strategy that we developed for assessment of rare sequence variants in ATM and CHEK2. In total, 12 distinct very rare or private variants were characterized in RAD51, with 10 cases (0.75%) and 9 controls (0.80%) carrying such a variant. Variants were either likely neutral missense substitutions (3), silent substitutions (4) or non-coding substitutions (5) that were predicted to have little effect on efficiency of the splicing machinery.

Conclusion

Altogether, our data suggest that RAD51 tolerates so little dysfunctional sequence variation that rare variants in the gene contribute little, if anything, to breast cancer susceptibility.

Calibration of multiple in silico tools for predicting pathogenicity of mismatch repair gene missense substitutions

Classification of rare missense substitutions observed during genetic testing for patient management is a considerable problem in clinical genetics. The Bayesian integrated evaluation of unclassified variants is a solution originally developed for BRCA1/2. Here, we take a step toward an analogous system for the mismatch repair (MMR) genes (MLH1, MSH2, MSH6, and PMS2) that confer colon cancer susceptibility in Lynch syndrome by calibrating in silico tools to estimate prior probabilities of pathogenicity for MMR gene missense substitutions. A qualitative five-class classification system was developed and applied to 143 MMR missense variants. This identified 74 missense substitutions suitable for calibration. These substitutions were scored using six different in silico tools (Align-Grantham Variation Grantham Deviation, multivariate analysis of protein polymorphisms [MAPP], MutPred, PolyPhen-2.1, Sorting Intolerant From Tolerant, and Xvar), using curated MMR multiple sequence alignments where possible. The output from each tool was calibrated by regression against the classifications of the 74 missense substitutions; these calibrated outputs are interpretable as prior probabilities of pathogenicity. MAPP was the most accurate tool and MAPP + PolyPhen-2.1 provided the best-combined model (R(2)  = 0.62 and area under receiver operating characteristic = 0.93). The MAPP + PolyPhen-2.1 output is sufficiently predictive to feed as a continuous variable into the quantitative Bayesian integrated evaluation for clinical classification of MMR gene missense substitutions.

A multifactorial likelihood model for MMR gene variant classification incorporating probabilities based on sequence bioinformatics and tumor characteristics: a report from the Colon Cancer Family Registry

Mismatch repair (MMR) gene sequence variants of uncertain clinical significance are often identified in suspected Lynch syndrome families, and this constitutes a challenge for both researchers and clinicians. Multifactorial likelihood model approaches provide a quantitative measure of MMR variant pathogenicity, but first require input of likelihood ratios (LRs) for different MMR variation-associated characteristics from appropriate, well-characterized reference datasets. Microsatellite instability (MSI) and somatic BRAF tumor data for unselected colorectal cancer probands of known pathogenic variant status were used to derive LRs for tumor characteristics using the Colon Cancer Family Registry (CFR) resource. These tumor LRs were combined with variant segregation within families, and estimates of prior probability of pathogenicity based on sequence conservation and position, to analyze 44 unclassified variants identified initially in Australasian Colon CFR families. In addition, in vitro splicing analyses were conducted on the subset of variants based on bioinformatic splicing predictions. The LR in favor of pathogenicity was estimated to be ~12-fold for a colorectal tumor with a BRAF mutation-negative MSI-H phenotype. For 31 of the 44 variants, the posterior probabilities of pathogenicity were such that altered clinical management would be indicated. Our findings provide a working multifactorial likelihood model for classification that carefully considers mode of ascertainment for gene testing.

Abstract 4444: BRCA1 Circos: A visualization tool for BRCA1 missense variants

Germline mutations in BRCA1 have been associated with a significantly increased risk of developing breast and ovarian cancers. A large number of variants have been identified in BRCA1 and are reported as pathogenic, non-pathogenic, or as variants of uncertain significance (VUS). Test results reporting pathogenic or non-pathogenic variants can be used by physicians and genetic counselors to help patients in making important clinical decisions. Around 10-15% of the reported BRCA1 and BRCA2 variants are clearly pathogenic and 5% of these reported variants are VUS. Results reporting a VUS can cause anxiety and make it difficult to advise the patient on further clinical intervention to reduce their risk of developing the disease. Many of these VUSs are non-synonymous single amino acid changes leading to missense changes. A number of previous studies have focused on assessing the functional impact of these changes on the BRCA1 protein. BRCA1 is a gene product with pleiotropic functions and therefore, different functional assays are used in these studies to analyze the variant's impact on function. However, information on the functional impact of VUS is currently not consolidated and collecting available information on missense variants depends on multiple database searches and literature review. Moreover, methods to incorporate functional data into multifactorial models have recently been developed highlighting the need for a common resource to facilitate analysis. Realizing the need for a resource that facilitates the global visualization and analysis of these variants we have developed a visualization tool for multidimensional data using Circos (http://circos.ca) to aggregate the data from all published studies that have assessed these missense variants. Circos plots have been used in the literature to represent genomic data and we have used the same concept to illustrate variants in BRCA1. We have included ∼700 variants for which we have 30 different functional experiments from various studies assessing over 150 missense variants. Results are color coded to uniformly depict various degrees of functional impact on the BRCA1 protein function so that information for any variant can be obtained in a snapshot without having to perform extensive literature searches. Therefore, we believe clinical geneticists and genetic counselors can use this research tool as a quick reference to analyze all the missense variants that have been identified along with associated information to characterize them when faced with the challenge of advising a patient with these types of VUSs.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4444. doi:1538-7445.AM2012-4444