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.

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.

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 review of a multifactorial probability-based model for classification of BRCA1 and BRCA2 variants of uncertain significance (VUS).

Clinical mutation screening of the BRCA1 and BRCA2 genes for the presence of germline inactivating mutations is used to identify individuals at elevated risk of breast and ovarian cancer. Variants identified during screening are usually classified as pathogenic (increased risk of cancer) or not pathogenic (no increased risk of cancer). However, a significant proportion of genetic tests yields variants of uncertain significance (VUS) that have undefined risk of cancer. Individuals carrying these VUS cannot benefit from individualized cancer risk assessment. Recently, a quantitative "posterior probability model" for assessing the clinical relevance of VUS in BRCA1 or BRCA2, which integrates multiple forms of genetic evidence has been developed. Here, we provide a detailed review of this model. We describe the components of the model and explain how these can be combined to calculate a posterior probability of pathogenicity for each VUS. We explain how the model can be applied to public data and provide tables that list the VUS that have been classified as not pathogenic or pathogenic using this method. While we use BRCA1 and BRCA2 VUS as examples, the method can be used as a framework for classification of the pathogenicity of VUS in other cancer genes.

Rare germline mutations in PALB2 and breast cancer risk: a population-based study.

Germline mutations in the PALB2 gene are associated with an increased risk of developing breast cancer but little is known about the frequencies of rare variants in PALB2 and the nature of the variants that influence risk. We selected participants recruited to the Women's Environment, Cancer, and Radiation Epidemiology (WECARE) Study and screened lymphocyte DNA from cases with contralateral breast cancer (n = 559) and matched controls with unilateral breast cancer (n = 565) for PALB2 mutations. Five pathogenic PALB2 mutations were identified among the cases (0.9%) versus none among the controls (P = 0.04). The first-degree female relatives of these five carriers demonstrated significantly higher incidence of breast cancer than relatives of noncarrier cases, indicating that pathogenic PALB2 mutations confer an estimated 5.3-fold increase in risk (95% CI: 1.8-13.2). The frequency of rare (<1% MAF) missense mutations was similar in both groups (23 vs. 21). Our findings confirm in a population-based study setting of women with breast cancer the strong risk associated with truncating mutations in PALB2 that has been reported in family studies. Conversely, there is no evidence from this study that rare PALB2 missense mutations strongly influence breast cancer risk.

Classification of missense substitutions in the BRCA genes: a database dedicated to Ex-UVs.

Unclassified sequence variants (UVs) arising from clinical mutation screening of cancer susceptibility genes present a frustrating issue to clinical genetics services and the patients that they serve. We created an open-access database holding missense substitutions from the breast and ovarian cancer susceptibility genes BRCA1 and BRCA2. The main inclusion criterion is that each variant should have been assessed in a published work that used the Bayesian integrated evaluation of unclassified BRCA gene variants. Transfer of data on these substitutions from the original publications to our database afforded an opportunity to analyze the missense substitutions under a single model and to remove inconsistencies that arose during the evolution of the integrated evaluation over the last decade. This analysis also afforded the opportunity to reclassify these missense substitutions according to the recently published IARC 5-Class system. From an initial set of 248 missense substitutions, 31 were set aside due to nonnegligible probability to interfere with splicing. Of the remaining substitutions, 28 fell into one of the two pathogenic classes (IARC Class 4 or 5), 174 fell into one of the two nonpathogenic classes (IARC Class 1 or 2), and 15 remain in IARC Class 3, "Uncertain." The database is available at http://brca.iarc.fr/LOVD.

Rare, evolutionarily unlikely missense substitutions in ATM confer increased risk of breast cancer.

The susceptibility gene for ataxia telangiectasia, ATM, is also an intermediate-risk breast-cancer-susceptibility gene. However, the spectrum and frequency distribution of ATM mutations that confer increased risk of breast cancer have been controversial. To assess the contribution of rare variants in this gene to risk of breast cancer, we pooled data from seven published ATM case-control mutation-screening studies, including a total of 1544 breast cancer cases and 1224 controls, with data from our own mutation screening of an additional 987 breast cancer cases and 1021 controls. Using an in silico missense-substitution analysis that provides a ranking of missense substitutions from evolutionarily most likely to least likely, we carried out analyses of protein-truncating variants, splice-junction variants, and rare missense variants. We found marginal evidence that the combination of ATM protein-truncating and splice-junction variants contribute to breast cancer risk. There was stronger evidence that a subset of rare, evolutionarily unlikely missense substitutions confer increased risk. On the basis of subset analyses, we hypothesize that rare missense substitutions falling in and around the FAT, kinase, and FATC domains of the protein may be disproportionately responsible for that risk and that a subset of these may confer higher risk than do protein-truncating variants. We conclude that a comparison between the graded distributions of missense substitutions in cases versus controls can complement analyses of truncating variants and help identify susceptibility genes and that this approach will aid interpretation of the data emerging from new sequencing technologies.

KRAS mutations in blood circulating cell-free DNA: a pancreatic cancer case-control.

The utility of KRAS mutations in plasma circulating cell-free DNA (cfDNA) samples as non-invasive biomarkers for the detection of pancreatic cancer has never been evaluated in a large case-control series. We applied a KRAS amplicon-based deep sequencing strategy combined with analytical pipeline specifically designed for the detection of low-abundance mutations to screen plasma samples of 437 pancreatic cancer cases, 141 chronic pancreatitis subjects, and 394 healthy controls. We detected mutations in 21.1% (N=92) of cases, of whom 82 (89.1%) carried at least one mutation at hotspot codons 12, 13 or 61, with mutant allelic fractions from 0.08% to 79%. Advanced stages were associated with an increased proportion of detection, with KRAS cfDNA mutations detected in 10.3%, 17,5% and 33.3% of cases with local, regional and systemic stages, respectively. We also detected KRAS cfDNA mutations in 3.7% (N=14) of healthy controls and in 4.3% (N=6) of subjects with chronic pancreatitis, but at significantly lower allelic fractions than in cases. Combining cfDNA KRAS mutations and CA19-9 plasma levels on a limited set of case-control samples did not improve the overall performance of the biomarkers as compared to CA19-9 alone. Whether the limited sensitivity and specificity observed in our series of KRAS mutations in plasma cfDNA as biomarkers for pancreatic cancer detection are attributable to methodological limitations or to the biology of cfDNA should be further assessed in large case-control series.

The 12p13.33/RAD52 locus and genetic susceptibility to squamous cell cancers of upper aerodigestive tract.

Genetic variants located within the 12p13.33/RAD52 locus have been associated with lung squamous cell carcinoma (LUSC). Here, within 5,947 UADT cancers and 7,789 controls from 9 different studies, we found rs10849605, a common intronic variant in RAD52, to be also associated with upper aerodigestive tract (UADT) squamous cell carcinoma cases (OR = 1.09, 95% CI: 1.04-1.15, p = 6x10(-4)). We additionally identified rs10849605 as a RAD52 cis-eQTL inUADT(p = 1x10(-3)) and LUSC (p = 9x10(-4)) tumours, with the UADT/LUSC risk allele correlated with increased RAD52 expression levels. The 12p13.33 locus, encompassing rs10849605/RAD52, was identified as a significant somatic focal copy number amplification in UADT(n = 374, q-value = 0.075) and LUSC (n = 464, q-value = 0.007) tumors and correlated with higher RAD52 tumor expression levels (p = 6x10(-48) and p = 3x10(-29) in UADT and LUSC, respectively). In combination, these results implicate increased RAD52 expression in both genetic susceptibility and tumorigenesis of UADT and LUSC tumors.

A rare truncating BRCA2 variant and genetic susceptibility to upper aerodigestive tract cancer.

Deleterious BRCA2 genetic variants markedly increase risk of developing breast cancer. A rare truncating BRCA2 genetic variant, rs11571833 (K3326X), has been associated with a 2.5-fold risk of lung squamous cell carcinoma but only a modest 26% increase in breast cancer risk. We analyzed the association between BRCA2 SNP rs11571833 and upper aerodigestive tract (UADT) cancer risk with multivariable unconditional logistic regression adjusted by sex and combinations of study and country for 5942 UADT squamous cell carcinoma case patients and 8086 control patients from nine different studies. All statistical tests were two-sided. rs11571833 was associated with UADT cancers (odds ratio = 2.53, 95% confidence interval = 1.89 to 3.38, P = 3x10(-10)) and was present in European, Latin American, and Indian populations but extremely rare in Japanese populations. The association appeared more apparent in smokers (current or former) compared with never smokers (P het = .026). A robust association between a truncating BRCA2 variant and UADT cancer risk suggests that treatment strategies orientated towards BRCA2 mutations may warrant further investigation in UADT tumors.

Modelling mutational landscapes of human cancers in vitro.

Experimental models that recapitulate mutational landscapes of human cancers are needed to decipher the rapidly expanding data on human somatic mutations. We demonstrate that mutation patterns in immortalised cell lines derived from primary murine embryonic fibroblasts (MEFs) exposed in vitro to carcinogens recapitulate key features of mutational signatures observed in human cancers. In experiments with several cancer-causing agents we obtained high genome-wide concordance between human tumour mutation data and in vitro data with respect to predominant substitution types, strand bias and sequence context. Moreover, we found signature mutations in well-studied human cancer driver genes. To explore endogenous mutagenesis, we used MEFs ectopically expressing activation-induced cytidine deaminase (AID) and observed an excess of AID signature mutations in immortalised cell lines compared to their non-transgenic counterparts. MEF immortalisation is thus a simple and powerful strategy for modelling cancer mutation landscapes that facilitates the interpretation of human tumour genome-wide sequencing data.

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.