BRAF and NRAS locus-specific variants have different outcomes on survival to colorectal cancer

Purpose: Somatic mutation status at KRAS, BRAF and NRAS is associated with prognosis in patients with advanced colorectal cancer (aCRC); however, it remains unclear whether there are intra-locus, variant-specific differences in survival and other clinicopathological parameters. Experimental design: We profiled 2,157 aCRCs for somatic mutations in KRAS, BRAF and NRAS and determined microsatellite instability status. We sought inter- and intra-locus correlations between mutations, and variant-specific associations with survival and clinicopathology. Results: KRAS mutations were rarely found together and those in codons 12 and 13 conferred poor prognosis (HR 1.44, 95% CI 1.28-1.61, p=6.4e-10 and HR 1.53, 95% CI 1.26-1.86, p=1.5e-05, respectively). For BRAF, more c.1781A>G (p.D594G) CRCs carried RAS mutations (14% [3/21]) compared to c.1799T>A (p.V600E) CRCs (1% [2/178], p=9.0e-03). c.1799T>A (p.V600E) was associated with poor prognosis (HR 2.60, 95% CI 2.06-3.28, p=1.0e-15), whereas c.1781A>G (p.D594G) was not (HR 1.30, 95% CI 0.73-2.31, p=0.37); this intra-locus difference was significant (p=0.04). More c.1799T>A (p.V600E) CRCs were found in the right colon (47% [47/100]), compared to c.1781A>G (p.D594G) CRCs (7% [1/15], p=3.7e-03). For NRAS, 5% (3/60) of codon 61 mutant CRCs had KRAS mutations compared to 44% (10/23) of codons 12 and 13 mutant CRCs (p=7.9e-05). Codon 61 mutations conferred poor prognosis (HR 1.47, 95% CI 1.09-1.99, p=0.01), whereas codons 12 and 13 mutations did not (HR 1.29, 95% CI 0.64-2.58, p=0.48). Conclusions: Our data show considerable intra-locus variation in the outcomes of mutations in BRAF and NRAS. These data need to be considered in patient management and personalised cancer therapy

Functional characterization of two human MutY homolog (hMYH) missense mutations (R227W and V232F) that lie within the putative hMSH6 binding domain and are associated with hMYH polyposis

The base excision repair DNA glycosylase MutY homolog (MYH) is responsible for removing adenines misincorporated into DNA opposite guanine or 7,8-dihydro-8-oxo-guanine (8-oxoG), thereby preventing G:C to T:A mutations. Biallelic germline mutations in the human MYH gene predispose individuals to multiple colorectal adenomas and carcinoma. We have recently demonstrated that hMYH interacts with the mismatch repair protein hMSH6, and that the hMSH2/hMSH6 (hMutSα) heterodimer stimulates hMYH activity. Here, we characterize the functional effect of two missense mutations (R227W and V232F) associated with hMYH polyposis that lie within, or adjacent to, the putative hMSH6 binding domain. Neither missense mutation affects the physical interaction between hMYH and hMSH6. However, hMYH(R227W) has a severe defect in A/8-oxoG binding and glycosylase activities, while hMYH(V232F) has reduced A/8-oxoG binding and glycosylase activities. The glycosylase activity of the V232F mutant can be partially stimulated by hMutSα but cannot be restored to the wild-type level. Both mutants also fail to complement mutY-deficiency in Escherichia coli. These data define the pathogenic mechanisms underlying two further hMYH polyposis-associated mutations

Increased frequency of the k-ras G12C mutation in MYH polyposis colorectal adenomas

Colorectal tumours from MYH polyposis patients display an excess of somatic G : C right arrow T : A transversions in the adenomatous polyposis coli gene. Here, we identify k-ras mutations in nine out of 54 (16.7%) MYH polyposis tumours. Their presence was associated with increased dysplasia and tubulovillous morphology (P=0.005). G : C right arrow T : A transversions in k-ras were significantly more frequent in MYH polyposis adenomas than in sporadic or familial adenomatous polyposis-associated tumours (Pless than or equal to0.002), and all resulted in a glycine-to-cysteine substitution at codon 12

Pan-cancer interrogation of MUTYH variants reveals biallelic inactivation and defective base excision repair across a spectrum of solid tumors

Purpose Biallelic germline pathogenic variants of the base excision repair (BER) pathway gene MUTYH predispose to colorectal cancer (CRC) and other cancers. The possible association of heterozygous variants with broader cancer susceptibility remains uncertain. This study investigated the prevalence and consequences of pathogenic MUTYH variants and MUTYH loss of heterozygosity (LOH) in a large pan-cancer analysis. Materials and Methods Data from 354,366 solid tumor biopsies that were sequenced as part of routine clinical care were analyzed using a validated algorithm to distinguish germline from somatic MUTYH variants. Results Biallelic germline pathogenic MUTYH variants were identified in 119 tissue biopsies. Most were CRCs and showed increased tumor mutational burden (TMB) and a mutational signature consistent with defective BER (COSMIC Signature SBS18). Germline heterozygous pathogenic variants were identified in 5,991 biopsies and their prevalence was modestly elevated in some cancer types. About 12% of these cancers (738 samples: including adrenal gland cancers, pancreatic islet cell tumors, nonglioma CNS tumors, GI stromal tumors, and thyroid cancers) showed somatic LOH for MUTYH, higher rates of chromosome 1p loss (where MUTYH is located), elevated genomic LOH, and higher COSMIC SBS18 signature scores, consistent with BER deficiency. Conclusion This analysis of MUTYH alterations in a large set of solid cancers suggests that in addition to the established role of biallelic pathogenic MUTYH variants in cancer predisposition, a broader range of cancers may possibly arise in MUTYH heterozygotes via a mechanism involving somatic LOH at the MUTYH locus and defective BER. However, the effect is modest and requires confirmation in additional studies before being clinically actionable

Comprehensive analysis of colorectal cancer-risk loci and survival outcome: a prognostic role for CDH1 variants.

PURPOSE: Genome-wide association studies have identified common single nucleotide polymorphisms (SNPs) at 83 loci associated with colorectal cancer (CRC) risk in European populations. Because germline variation can also influence patient outcome, we studied the relationship between these SNPs and CRC survivorship. EXPERIMENTAL DESIGN: For the 83 risk loci, 10 lead SNPs were directly genotyped, 72 were imputed and 1 was not genotyped nor imputed, in 1948 unrelated patients with advanced CRC from the clinical trials COIN and COIN-B (oxaliplatin and fluoropyrimidine chemotherapy ± cetuximab). A Cox survival model was used for each variant, and variants classified by pathway, adjusting for known prognostic factors. We imposed a Bonferroni threshold of P = 6.6 × 10-4 for multiple testing. We carried out meta-analyses of published risk SNPs associated with survival. RESULTS: Univariate analysis identified six SNPs associated with overall survival (OS) (P < 0.05); however, only rs9939049 in CDH1 remained significant beyond the Bonferroni threshold (Hazard Ratio [HR] 1.44, 95% Confidence Intervals [CI]: 1.21-1.71, P = 5.0 × 10-5). Fine mapping showed that rs12597188 was the most significant SNP at this locus and remained significant after adjustment for known prognostic factors beyond multiple testing thresholds (HR 1.23, 95% CI: 1.13-1.34, P = 1.9 × 10-6). rs12597188 was also associated with poor response to therapy (OR 0.61, 95% CI: 0.42-0.87, P = 6.6 × 10-3). No combinations of SNPs within pathways were more significantly associated with survival compared with single variants alone, and no other risk SNPs were associated with survival in meta-analyses. CONCLUSIONS: The CRC susceptibility SNP rs9939049 in CDH1 influences patient survival and warrants further evaluation as a prognostic biomarker

Comprehensive pharmacogenetic profiling of the epidermal growth factor receptor pathway for biomarkers of response to, and toxicity from, cetuximab

Background Somatic mutations in the epidermal growth factor receptor (EGFR) intracellular signalling pathways predict non-response to cetuximab in the treatment of advanced colorectal cancer (aCRC). We hypothesized that common germline variants within these pathways may also play similar roles. Methods We analysed 54 potentially functional, common, inherited EGFR pathway variants in 815 aCRC patients treated with oxaliplatin-fluoropyrimidine chemotherapy +cetuximab. Primary endpoints were response and skin rash (SR). We had >85% power to detect ORs=1.6 for variants with minor allele frequencies >20%. Results We identified five potential biomarkers for response and four for SR, although none remained significant after correction for multiple testing. Our initial data supported a role for Ser313Pro in PIK3R2 in modulating response to cetuximab - in patients with KRAS wild type CRCs, 36.4% of patients with one allele encoding proline responded, as compared to 71.2% of patients homozygous for alleles encoding serine (OR 0.23, 95% CI 0.09-0.56, P=0.0014) and this association was predictive for cetuximab (Pinteraction=0.017); however, independent replication failed to validate this association. No previously proposed predictive biomarkers were validated. Conclusions Our study highlights the need to validate potential pharmacogenetic biomarkers. We did not find strong evidence for common germline biomarkers of cetuximab response and toxicity

Relationship between 233 colorectal cancer risk loci and survival in 1926 patients with advanced disease

Background: Genome, transcriptome and methylome-wide association studies have identified single-nucleotide polymorphisms (SNPs) or genes at 258 loci associated with colorectal cancer (CRC) risk. We studied the relationship between these and patient outcome. Methods: We studied 1926 unrelated patients with advanced CRC from COIN and COIN-B. Of 205 CRC-risk SNPs, 19 were directly genotyped and 162 were imputed, and of 53 risk genes, 52 were tested. An additive Cox model for overall survival was adjusted for known prognostic factors. For nominally significant SNPs or genes, we considered a recessive model with a Bonferroni corrected threshold of P = 2.1 × 10−4. We examined SNPs as expression quantitative trait loci (eQTL) and the relationship between gene expression in colorectal tumours and survival in 597 unrelated patients. Results: Eleven SNPs or genes were nominally associated with survival under an additive model. Only rs117079142 mapping to UTP23 and EIF3H (Hazard Ratio [HR] = 2.79, 95% Confidence Intervals [CI] = 1.70–4.58, P = 4.7 × 10−5) and rs9924886 mapping to CDH1 and CDH3 (HR = 1.24, 95% CI = 1.12–1.38, P = 5.2 × 10−5) passed the multiple testing threshold under a recessive model. rs117079142 was an eQTL for UTP23 and rs9924886 for CDH1, CDH3 and ZFP90. Decreased CDH1 expression in CRCs was associated with worse survival (HR = 2.18, 95% CI = 1.3–3.5, P = 1.8 × 10−3). Conclusion: rs117079142 and rs9924886 may represent potential prognostic biomarkers for CRC

Lack of an association between gallstone disease and bilirubin levels with risk of colorectal cancer : a Mendelian randomisation analysis

BACKGROUND: Epidemiological studies of the relationship between gallstone disease and circulating levels of bilirubin with risk of developing colorectal cancer (CRC) have been inconsistent. To address possible confounding and reverse causation, we examine the relationship between these potential risk factors and CRC using Mendelian randomisation (MR). METHODS: We used two-sample MR to examine the relationship between genetic liability to gallstone disease and circulating levels of bilirubin with CRC in 26,397 patients and 41,481 controls. We calculated the odds ratio per genetically predicted SD unit increase in log bilirubin levels (ORSD) for CRC and tested for a non-zero causal effect of gallstones on CRC. Sensitivity analysis was applied to identify violations of estimator assumptions. RESULTS: No association between either gallstone disease (P value = 0.60) or circulating levels of bilirubin (ORSD = 1.00, 95% confidence interval (CI) = 0.96-1.03, P value = 0.90) with CRC was shown. CONCLUSIONS: Despite the large scale of this study, we found no evidence for a causal relationship between either circulating levels of bilirubin or gallstone disease with risk of developing CRC. While the magnitude of effect suggested by some observational studies can confidently be excluded, we cannot exclude the possibility of smaller effect sizes and non-linear relationships.Peer reviewe

Modifiable pathways for colorectal cancer : a mendelian randomisation analysis

Background Epidemiological studies have linked lifestyle, cardiometabolic, reproductive, developmental, and inflammatory factors to the risk of colorectal cancer. However, which specific factors affect risk and the strength of these effects are unknown. We aimed to examine the relationship between potentially modifiable risk factors and colorectal cancer. Methods We used a random-effects model to examine the relationship between 39 potentially modifiable risk factors and colorectal cancer in 26 397 patients with colorectal cancer and 41 481 controls (ie, people without colorectal cancer). These population data came from a genome-wide association study of people of European ancestry, which was amended to exclude UK BioBank data. In the model, we used genetic variants as instruments via two-sample mendelian randomisation to limit bias from confounding and reverse causation. We calculated odds ratios per genetically predicted SD unit increase in each putative risk factor (OR SD) for colorectal cancer risk. We did mendelian randomisation Egger regressions to identify evidence of potential violations of mendelian randomisation assumptions. A Bonferroni-corrected threshold of p=1.3 x 10(-3) was considered significant, and p values less than 0.05 were considered to be suggestive of an association. Findings No putative risk factors were significantly associated with colorectal cancer risk after correction for multiple testing. However, suggestive associations with increased risk were noted for genetically predicted body fat percentage (OR SD 1.14 [95% CI 1.03-1.25]; p=0.0086), body-mass index (1.09 [1.01-1.17]; p=0.023), waist circumference (1.13 [1.02-1.26]; p=0.018), basal metabolic rate (1.10 [1.03-1.18]; p=0.0079), and concentrations of LDL cholesterol (1.14 [1.04-1.25]; p=0.0056), total cholesterol (1.09 [1.01-1.18]; p=0.025), circulating serum iron (1.17 [1.00-1.36]; p=0.049), and serum vitamin B12 (1.21 [1.04-1.42]; p=0.016), although potential pleiotropy among genetic variants used as instruments for vitamin B12 constrains the finding. A suggestive association was also noted between adult height and increased risk of colorectal cancer (OR SD 1.04 [95% CI 1.00-1.08]; p=0.032). Low blood selenium concentration had a suggestive association with decreased risk of colorectal cancer (OR SD 0.85 [95% CI 0.75-0.96]; p=0.0078) based on a single variant, as did plasma concentrations of interleukin-6 receptor subunit a (also based on a single variant; 0.98 [0.96-1.00]; p=0.035). Risk of colorectal cancer was not associated with any sex hormone or reproductive factor, serum calcium, or circulating 25-hydroxyvitamin D concentrations. Interpretation This analysis identified several modifiable targets for primary prevention of colorectal cancer, including lifestyle, obesity, and cardiometabolic factors, that should inform public health policy. Copyright (C) 2019 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license.Peer reviewe

Genetic variation in ST6GAL1 is a determinant of capecitabine and oxaliplatin induced hand-foot syndrome

Cancer patients treated with capecitabine and oxaliplatin (XELOX) often develop hand-foot syndrome (HFS) or palmar-plantar erythrodysesthesia. Genetic variation in ST6GAL1 is a risk factor for type-2 diabetes (T2D), a disease also associated with HFS. We analysed genome-wide association data for ten toxicities in advanced colorectal cancer (CRC) patients from the COIN and COIN-B trials. 1,055 patients were treated with XELOX ±cetuximab and 745 with folinic acid, fluorouracil and oxaliplatin ±cetuximab. We also analysed rs6783836 in ST6GAL1 with HFS in CRC patients from QUASAR2. Using UK Biobank data, we sought to confirm an association between ST6GAL1 and T2D (17,384 cases, 317,887 controls) and analysed rs6783836 against markers of diabetes, inflammation and psoriasis. We found that 68% of patients from COIN and COIN-B with grade 2-3 HFS responded to treatment as compared to 58% with grade 0-1 HFS (Odds Ratio [OR]=1.1, 95% Confidence Interval [CI]=1.02-1.2, P=2.0x10-4). HFS was also associated with improved overall survival (Hazard Ratio=0.92, 95%CI=0.84-0.99, P=4.6x10-2). rs6783836 at ST6GAL1 was associated with HFS in patients treated with XELOX (OR=3.1, 95%CI=2.1-4.6, P=4.3x10-8) and was borderline significant in patients receiving capecitabine from QUASAR2, but with an opposite allele effect (OR=0.66, 95% CI=0.42-1.03, P=0.05). ST6GAL1 was associated with T2D (lead SNP rs3887925, OR=0.94, 95%CI=0.92-0.96, P=1.2x10-8) and the rs6783836-T allele was associated with lowered HbA1c levels (P=5.9x10-3) and lymphocyte count (P=2.7x10-3), and psoriasis (P=7.5x10-3) beyond thresholds for multiple testing. In conclusion, HFS is a biomarker of treatment outcome and rs6783836 in ST6GAL1 is a potential biomarker for HFS with links to T2D and inflammation