Cumulative impact of common genetic variants and other risk factors on colorectal cancer risk in 42,103 individuals.

OBJECTIVE: Colorectal cancer (CRC) has a substantial heritable component. Common genetic variation has been shown to contribute to CRC risk. A study was conducted in a large multi-population study to assess the feasibility of CRC risk prediction using common genetic variant data combined with other risk factors. A risk prediction model was built and applied to the Scottish population using available data. DESIGN: Nine populations of European descent were studied to develop and validate CRC risk prediction models. Binary logistic regression was used to assess the combined effect of age, gender, family history (FH) and genotypes at 10 susceptibility loci that individually only modestly influence CRC risk. Risk models were generated from case-control data incorporating genotypes alone (n=39,266) and in combination with gender, age and FH (n=11,324). Model discriminatory performance was assessed using 10-fold internal cross-validation and externally using 4187 independent samples. The 10-year absolute risk was estimated by modelling genotype and FH with age- and gender-specific population risks. RESULTS: The median number of risk alleles was greater in cases than controls (10 vs 9, p<2.2 × 10(-16)), confirmed in external validation sets (Sweden p=1.2 × 10(-6), Finland p=2 × 10(-5)). The mean per-allele increase in risk was 9% (OR 1.09; 95% CI 1.05 to 1.13). Discriminative performance was poor across the risk spectrum (area under curve for genotypes alone 0.57; area under curve for genotype/age/gender/FH 0.59). However, modelling genotype data, FH, age and gender with Scottish population data shows the practicalities of identifying a subgroup with >5% predicted 10-year absolute risk. CONCLUSION: Genotype data provide additional information that complements age, gender and FH as risk factors, but individualised genetic risk prediction is not currently feasible. Nonetheless, the modelling exercise suggests public health potential since it is possible to stratify the population into CRC risk categories, thereby informing targeted prevention and surveillance.

Candidate driver genes in microsatellite-unstable colorectal cancer.

Defects in the mismatch repair system lead to microsatellite instability (MSI), a feature observed in ∼ 15% of all colorectal cancers (CRCs). Microsatellite mutations that drive tumourigenesis, typically inactivation of tumour suppressors, are selected for and are frequently detected in MSI cancers. Here, we evaluated somatic mutations in microsatellite repeats of 790 genes chosen based on reduced expression in MSI CRC and existence of a coding mononucleotide repeat of 6-10 bp in length. All the repeats were initially sequenced in 30 primary MSI CRC samples and whenever frameshift mutations were identified in >20%, additional 70 samples were sequenced. To distinguish driver mutations from passengers, we similarly analyzed the occurrence of frameshift mutations in 121 intronic control repeats and utilized a statistical regression model to determine cut-off mutation frequencies for repeats of all types (A/T and C/G, 6-10 bp). Along with several know target genes, including TGFBR2, ACVR2, and MSH3, six novel candidate driver genes emerged that harbored significantly more mutations than identical control repeats. The mutation frequencies in 100 MSI CRC samples were 51% in G8 of GLYR1, 47% in T9 of ABCC5, 43% in G8 of WDTC1, 33% in A8 of ROCK1, 30% in T8 of OR51E2, and 28% in A8 of TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in tumors carrying biallelic mutations, supporting a loss of function hypothesis. This is a large scale, unbiased effort to identify genes that when mutated are likely to contribute to MSI CRC development.

High frequency of TTK mutations in microsatellite-unstable colorectal cancer and evaluation of their effect on spindle assembly checkpoint.

Frameshift mutations frequently accumulate in microsatellite-unstable colorectal cancers (MSI CRCs) typically leading to downregulation of the target genes due to nonsense-mediated messenger RNA decay. However, frameshift mutations that occur in the 3' end of the coding regions can escape decay, which has largely been ignored in previous works. In this study, we characterized nonsense-mediated decay-escaping frameshift mutations in MSI CRC in an unbiased, genome wide manner. Combining bioinformatic search with expression profiling, we identified genes that were predicted to escape decay after a deletion in a microsatellite repeat. These repeats, located in 258 genes, were initially sequenced in 30 MSI CRC samples. The mitotic checkpoint kinase TTK was found to harbor decay-escaping heterozygous mutations in exon 22 in 59% (105/179) of MSI CRCs, which is notably more than previously reported. Additional novel deletions were found in exon 5, raising the mutation frequency to 66%. The exon 22 of TTK contains an A(9)-G(4)-A(7) locus, in which the most common mutation was a mononucleotide deletion in the A(9) (c.2560delA). When compared with identical non-coding repeats, TTK was found to be mutated significantly more often than expected without selective advantage. Since TTK inhibition is known to induce override of the mitotic spindle assembly checkpoint (SAC), we challenged mutated cancer cells with the microtubule-stabilizing drug paclitaxel. No evidence of checkpoint weakening was observed. As a conclusion, heterozygous TTK mutations occur at a high frequency in MSI CRCs. Unexpectedly, the plausible selective advantage in tumourigenesis does not appear to be related to SAC.

Systematic search for enhancer elements and somatic allelic imbalance at seven low-penetrance colorectal cancer predisposition loci.

BACKGROUND: Common single-nucleotide polymorphisms (SNPs) in ten chromosomal loci have been shown to predispose to colorectal cancer (CRC) in genome-wide association studies. A plausible biological mechanism of CRC susceptibility associated with genetic variation has so far only been proposed for three loci, each pointing to variants that affect gene expression through distant regulatory elements. In this study, we aimed to gain insight into the molecular basis of seven low-penetrance CRC loci tagged by rs4779584 at 15q13, rs10795668 at 10p14, rs3802842 at 11q23, rs4444235 at 14q22, rs9929218 at 16q22, rs10411210 at 19q13, and rs961253 at 20p12. METHODS: Possible somatic gain of the risk allele or loss of the protective allele was studied by analyzing allelic imbalance in tumour and corresponding normal tissue samples of heterozygous patients. Functional variants were searched from in silico predicted enhancer elements locating inside the CRC-associating linkage-disequilibrium regions. RESULTS: No allelic imbalance targeting the SNPs was observed at any of the seven loci. Altogether, 12 SNPs that were predicted to disrupt potential transcription factor binding sequences were genotyped in the same population-based case-control series as the seven tagging SNPs originally. None showed association with CRC. CONCLUSIONS: The results of the allelic imbalance analysis suggest that the seven CRC risk variants are not somatically selected for in the neoplastic progression. The bioinformatic approach was unable to pinpoint cancer-causing variants at any of the seven loci. While it is possible that many of the predisposition loci for CRC are involved in control of gene expression by targeting transcription factor binding sites, also other possibilities, such as regulatory RNAs, should be considered.

Low-penetrance susceptibility variants in familial colorectal cancer.

BACKGROUND: Genomewide association studies have identified 10 low-penetrance loci that confer modestly increased risk for colorectal cancer (CRC). Although they underlie a significant proportion of CRC in the general population, their impact on the familial risk for CRC has yet to be formally enumerated. The aim of this study was to examine the combined contribution of the 10 variants, rs6983267, rs4779584, rs4939827, rs16892766, rs10795668, rs3802842, rs4444235, rs9929218, rs10411210, and rs961253, on familial CRC. METHODS: The population-based series of CRC samples included in this study consisted of 97 familial cases and 691 sporadic cases. Genotypes in the 10 loci and clinical data, including family history of cancer verified from the Finnish Cancer Registry, were available. The overall number of risk alleles (0-20) was determined, and its association with familial CRC was analyzed. Excess familial risk was estimated using cancer incidence data from the first-degree relatives of the cases. RESULTS: A linear association between the number of risk alleles and familial CRC was observed (P = 0.006). With each risk-allele addition, the odds of having an affected first-degree relative increased by 1.16 (95% confidence interval, 1.04-1.30). The 10 low-penetrance loci collectively explain approximately 9% of the variance in familial risk for CRC. CONCLUSIONS: This study provides evidence to support the previous indirect estimations that these low-penetrance variants account for a relatively small proportion of the familial aggregation of CRC. IMPACT: Our results emphasize the need to characterize the remaining molecular basis of familial CRC, which should eventually yield in individualized targeting of preventive interventions.

Mutations in the circadian gene CLOCK in colorectal cancer.

The circadian clock regulates daily variations in physiologic processes. CLOCK acts as a regulator in the circadian apparatus controlling the expression of other clock genes, including PER1. Clock genes have been implicated in cancer-related functions; in this work, we investigated CLOCK as a possible target of somatic mutations in microsatellite unstable colorectal cancers. Combining microarray gene expression data and public gene sequence information, we identified CLOCK as 1 of 790 putative novel microsatellite instability (MSI) target genes. A total of 101 MSI colorectal carcinomas (CRC) were sequenced for a coding microsatellite in CLOCK. The effect of restoring CLOCK expression was studied in LS180 cells lacking wild-type CLOCK by stably expressing GST-CLOCK or glutathione S-transferase empty vector and testing the effects of UV-induced apoptosis and radiation by DNA content analysis using flow cytometry. Putative novel CLOCK target genes were searched by using ChIP-seq. CLOCK mutations occurred in 53% of MSI CRCs. Restoring CLOCK expression in cells with biallelic CLOCK inactivation resulted in protection against UV-induced apoptosis and decreased G(2)-M arrest in response to ionizing radiation. Using ChIP-Seq, novel CLOCK-binding elements were identified near DNA damage genes p21, NBR1, BRCA1, and RAD50. CLOCK is shown to be mutated in cancer, and altered response to DNA damage provides one plausible mechanism of tumorigenesis.

Meta-analysis of three genome-wide association studies identifies susceptibility loci for colorectal cancer at 1q41, 3q26.2, 12q13.13 and 20q13.33.

Genome-wide association studies (GWAS) have identified ten loci harboring common variants that influence risk of developing colorectal cancer (CRC). To enhance the power to identify additional CRC risk loci, we conducted a meta-analysis of three GWAS from the UK which included a total of 3,334 affected individuals (cases) and 4,628 controls followed by multiple validation analyses including a total of 18,095 cases and 20,197 controls. We identified associations at four new CRC risk loci: 1q41 (rs6691170, odds ratio (OR) = 1.06, P = 9.55 × 10⁻¹° and rs6687758, OR = 1.09, P = 2.27 × 10⁻9, 3q26.2 (rs10936599, OR = 0.93, P = 3.39 × 10⁻8), 12q13.13 (rs11169552, OR = 0.92, P = 1.89 × 10⁻¹° and rs7136702, OR = 1.06, P = 4.02 × 10⁻8) and 20q13.33 (rs4925386, OR = 0.93, P = 1.89 × 10⁻¹°). In addition to identifying new CRC risk loci, this analysis provides evidence that additional CRC-associated variants of similar effect size remain to be discovered.

The common colorectal cancer predisposition SNP rs6983267 at chromosome 8q24 confers potential to enhanced Wnt signaling.

Homozygosity for the G allele of rs6983267 at 8q24 increases colorectal cancer (CRC) risk approximately 1.5 fold. We report here that the risk allele G shows copy number increase during CRC development. Our computer algorithm, Enhancer Element Locator (EEL), identified an enhancer element that contains rs6983267. The element drove expression of a reporter gene in a pattern that is consistent with regulation by the key CRC pathway Wnt. rs6983267 affects a binding site for the Wnt-regulated transcription factor TCF4, with the risk allele G showing stronger binding in vitro and in vivo. Genome-wide ChIP assay revealed the element as the strongest TCF4 binding site within 1 Mb of MYC. An unambiguous correlation between rs6983267 genotype and MYC expression was not detected, and additional work is required to scrutinize all possible targets of the enhancer. Our work provides evidence that the common CRC predisposition associated with 8q24 arises from enhanced responsiveness to Wnt signaling.

Allelic imbalance at rs6983267 suggests selection of the risk allele in somatic colorectal tumor evolution.

A common single nucleotide polymorphism (SNP), rs6983267, at 8q24.21 has recently been shown to associate with colorectal cancer (CRC). Three independent SNP association studies showed that rs6983267 contributes to CRC with odds ratios (OR) of 1.17 to 1.22. Here, we genotyped a population-based series of 1,042 patients with CRC and 1,012 healthy controls for rs6983267 and determined the contribution of SNP to CRC in Finland, using germ line DNA, as well as the respective cancer DNA in heterozygous patients. The comprehensive clinical data available from the 1,042 patients and their first-degree relatives enabled us to thoroughly examine the possible association of this variant with different clinical features. As expected, a significant association between the G allele of rs6983267 and CRC [OR, 1.22; 95% confidence interval (CI), 1.08-1.38; P = 0.0018] was found, confirming the previous observations. A trend towards association of the G allele with microsatellite-stable cancer (OR, 1.37; 95% CI, 1.02-1.85; P = 0.04) and family history of cancers other than CRC was seen (OR, 1.20; 95% CI, 1-1.43; P = 0.05). Four hundred and sixty-six GT heterozygotes identified in this study were analyzed for allelic imbalance at rs6983267 in the respective cancer DNA. One hundred and one tumors showed allelic imbalance (22%). The risk allele G was favored in 67 versus 34 tumors (P = 0.0007). This finding implicates that the underlying germ line genetic defect in 8q24.21 is a target in the somatic evolution of CRC.

Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer.

Genome-wide association (GWA) studies have identified multiple loci at which common variants modestly influence the risk of developing colorectal cancer (CRC). To enhance power to identify additional loci with similar effect sizes, we conducted a meta-analysis of two GWA studies, comprising 13,315 individuals genotyped for 38,710 common tagging SNPs. We undertook replication testing in up to eight independent case-control series comprising 27,418 subjects. We identified four previously unreported CRC risk loci at 14q22.2 (rs4444235, BMP4; P = 8.1 x 10(-10)), 16q22.1 (rs9929218, CDH1; P = 1.2 x 10(-8)), 19q13.1 (rs10411210, RHPN2; P = 4.6 x 10(-9)) and 20p12.3 (rs961253; P = 2.0 x 10(-10)). These findings underscore the value of large sample series for discovery and follow-up of genetic variants contributing to the etiology of CRC.