Iina Niittymäki

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 × 10−10), 16q22.1 (rs9929218, CDH1; P = 1.2 × 10−8), 19q13.1 (rs10411210, RHPN2; P = 4.6 × 10−9) and 20p12.3 (rs961253; P = 2.0 × 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.

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 ∼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.

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−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−10 and rs7136702, OR = 1.06, P = 4.02 × 10−8) and 20q13.33 (rs4925386, OR = 0.93, P = 1.89 × 10−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.

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.

Refinement of the basis and impact of common 11q23.1 variation to the risk of developing colorectal cancer

The common single-nucleotide polymorphism (SNP) rs3802842 at 11q23.1 has recently been reported to be associated with risk of colorectal cancer (CRC). To examine this association in detail we genotyped rs3802842 in eight independent case-control series comprising a total of 10 638 cases and 10 457 healthy individuals. A significant association between the C allele of rs3802842 and CRC risk was found (per allele OR = 1.17; 95% confidence interval [CI]: 1.12-1.22; P = 1.08 x 10(-12)) with the risk allele more frequent in rectal than colonic disease (P = 0.02). In combination with 8q21, 8q24, 10p14, 11q, 15q13.3 and 18q21 variants, the risk of CRC increases with an increasing numbers of variant alleles for the six loci (OR(per allele) = 1.19; 95% CI: 1.15-1.23; P(trend) = 7.4 x 10(-24)). Using the data from our genome-wide association study of CRC, LD mapping and imputation, we were able to refine the location of the causal locus to a 60 kb region and screened for coding changes. The absence of exonic mutations in any of the transcripts (FLJ45803, LOC120376, C11orf53 and POU2AF1) mapping to this region makes the association likely to be a consequence of non-coding effects on gene expression.

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.

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.

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 G2-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. Mol Cancer Res; 8(7); 952–60. ©2010 AACR.

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 ∼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. Cancer Epidemiol Biomarkers Prev; 19(6); 1478–83. ©2010 AACR.

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

BackgroundCommon 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.MethodsPossible 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.ResultsNo 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.ConclusionsThe 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.