Genome-wide association studies (GWAS) have identified 19 risk variants associated with colorectal cancer. As most of these risk variants reside outside the coding regions of genes, we conducted cis-expression quantitative trait loci (cis-eQTL) analyses to investigate possible regulatory functions on the expression of neighboring genes. Forty microsatellite stable and CpG island methylator phenotype-negative colorectal tumors and paired adjacent normal colon tissues were used for genome-wide SNP and gene expression profiling. We found that three risk variants (rs10795668, rs4444235 and rs9929218, using near perfect proxies rs706771, rs11623717 and rs2059252, respectively) were significantly associated (FDR q-value ≤0.05) with expression levels of nearby genes (<2 Mb up- or down-stream). We observed an association between the low colorectal cancer risk allele (A) for rs10795668 at 10p14 and increased expression of ATP5C1 (q = 0.024) and between the colorectal cancer high risk allele (C) for rs4444235 at 14q22.2 and increased expression of DLGAP5 (q = 0.041), both in tumor samples. The colorectal cancer low risk allele (A) for rs9929218 at 16q22.1 was associated with a significant decrease in expression of both NOL3 (q = 0.017) and DDX28 (q = 0.046) in the adjacent normal colon tissue samples. Of the four genes, DLGAP5 and NOL3 have been previously reported to play a role in colon carcinogenesis and ATP5C1 and DDX28 are mitochondrial proteins involved in cellular metabolism and division, respectively. The combination of GWAS findings, prior functional studies, and the cis-eQTL analyses described here suggest putative functional activities for three of the colorectal cancer GWAS identified risk loci as regulating the expression of neighboring genes.
References
[1]
Broderick P, Carvajal-Carmona L, Pittman AM, Webb E, Howarth K, et al. (2007) A genome-wide association study shows that common alleles of SMAD7 influence colorectal cancer risk. Nat Genet 39: 1315–1317.
[2]
Houlston RS, Cheadle J, Dobbins SE, Tenesa A, Jones AM, et al. (2010) Meta-analysis of three genome-wide association studies identifies susceptibility loci for colorectal cancer at 1q41, 3q26.2, 12q13.13 and 20q13.33. Nat Genet 42: 973–977.
[3]
Houlston RS, Webb E, Broderick P, Pittman AM, Di Bernardo MC, et al. (2008) Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer. Nat Genet 40: 1426–1435.
[4]
Tenesa A, Farrington SM, Prendergast JG, Porteous ME, Walker M, et al. (2008) Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21. Nat Genet 40: 631–637.
[5]
Tomlinson I, Webb E, Carvajal-Carmona L, Broderick P, Kemp Z, et al. (2007) A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21. Nat Genet 39: 984–988.
[6]
Tomlinson IP, Webb E, Carvajal-Carmona L, Broderick P, Howarth K, et al. (2008) A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3. Nat Genet 40: 623–630.
[7]
Zanke BW, Greenwood CM, Rangrej J, Kustra R, Tenesa A, et al. (2007) Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24. Nat Genet 39: 989–994.
[8]
Schadt EE, Lamb J, Yang X, Zhu J, Edwards S, et al. (2005) An integrative genomics approach to infer causal associations between gene expression and disease. Nat Genet 37: 710–717.
[9]
Cheung VG, Jen KY, Weber T, Morley M, Devlin JL, et al. (2003) Genetics of quantitative variation in human gene expression. Cold Spring Harb Symp Quant Biol 68: 403–407.
[10]
Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, et al. (2004) Genetic analysis of genome-wide variation in human gene expression. Nature 430: 743–747.
[11]
Schadt EE, Monks SA, Drake TA, Lusis AJ, Che N, et al. (2003) Genetics of gene expression surveyed in maize, mouse and man. Nature 422: 297–302.
[12]
Duan S, Huang RS, Zhang W, Bleibel WK, Roe CA, et al. (2008) Genetic architecture of transcript-level variation in humans. Am J Hum Genet 82: 1101–1113.
[13]
Stranger BE, Forrest MS, Clark AG, Minichiello MJ, Deutsch S, et al. (2005) Genome-wide associations of gene expression variation in humans. PLoS Genet 1: e78.
[14]
Nicolae DL, Gamazon E, Zhang W, Duan S, Dolan ME, et al. (2010) Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS. PLoS Genet 6: e1000888.
[15]
Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, et al. (2007) A genome-wide association study of global gene expression. Nat Genet 39: 1202–1207.
[16]
Veyrieras JB, Kudaravalli S, Kim SY, Dermitzakis ET, Gilad Y, et al. (2008) High-resolution mapping of expression-QTLs yields insight into human gene regulation. PLoS Genet 4: e1000214.
[17]
Musunuru K, Strong A, Frank-Kamenetsky M, Lee NE, Ahfeldt T, et al. (2010) From noncoding variant to phenotype via SORT1 at the 1p13 cholesterol locus. Nature 466: 714–719.
[18]
Pomerantz MM, Shrestha Y, Flavin RJ, Regan MM, Penney KL, et al. (2010) Analysis of the 10q11 cancer risk locus implicates MSMB and NCOA4 in human prostate tumorigenesis. PLoS Genet 6: e1001204.
[19]
Schadt EE, Molony C, Chudin E, Hao K, Yang X, et al. (2008) Mapping the genetic architecture of gene expression in human liver. PLoS Biol 6: e107.
[20]
Freedman ML, Monteiro AN, Gayther SA, Coetzee GA, Risch A, et al. (2011) Principles for the post-GWAS functional characterization of cancer risk loci. Nat Genet 43: 513–518.
[21]
Newcomb PA, Baron J, Cotterchio M, Gallinger S, Grove J, et al. (2007) Colon Cancer Family Registry: an international resource for studies of the genetic epidemiology of colon cancer. Cancer Epidemiol Biomarkers Prev 16: 2331–2343.
[22]
Lindor NM, Burgart LJ, Leontovich O, Goldberg RM, Cunningham JM, et al. (2002) Immunohistochemistry versus microsatellite instability testing in phenotyping colorectal tumors. J Clin Oncol 20: 1043–1048.
[23]
Weisenberger DJ, Siegmund KD, Campan M, Young J, Long TI, et al. (2006) CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer. Nat Genet 38: 787–793.
[24]
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, et al. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4: 249–264.
[25]
Benjamini Y, Hochber Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological) 57: 289–300.
[26]
Cuezva JM, Krajewska M, de Heredia ML, Krajewski S, Santamaria G, et al. (2002) The bioenergetic signature of cancer: a marker of tumor progression. Cancer Res 62: 6674–6681.
[27]
Gudmundsson KO, Thorsteinsson L, Sigurjonsson OE, Keller JR, Olafsson K, et al. (2007) Gene expression analysis of hematopoietic progenitor cells identifies Dlg7 as a potential stem cell gene. Stem Cells 25: 1498–1506.
[28]
Tsou AP, Yang CW, Huang CY, Yu RC, Lee YC, et al. (2003) Identification of a novel cell cycle regulated gene, HURP, overexpressed in human hepatocellular carcinoma. Oncogene 22: 298–307.
[29]
Willers IM, Isidoro A, Ortega AD, Fernandez PL, Cuezva JM (2010) Selective inhibition of beta-F1-ATPase mRNA translation in human tumours. Biochem J 426: 319–326.
[30]
Fearon ER, Cho KR, Nigro JM, Kern SE, Simons JW, et al. (1990) Identification of a chromosome 18q gene that is altered in colorectal cancers. Science 247: 49–56.
[31]
Funakoshi S, Ezaki T, Kong J, Guo RJ, Lynch JP (2008) Repression of the desmocollin 2 gene expression in human colon cancer cells is relieved by the homeodomain transcription factors Cdx1 and Cdx2. Mol Cancer Res 6: 1478–1490.
[32]
Powell SM, Zilz N, Beazer-Barclay Y, Bryan TM, Hamilton SR, et al. (1992) APC mutations occur early during colorectal tumorigenesis. Nature 359: 235–237.
[33]
Stein U, Walther W, Arlt F, Schwabe H, Smith J, et al. (2009) MACC1, a newly identified key regulator of HGF-MET signaling, predicts colon cancer metastasis. Nat Med 15: 59–67.
[34]
Warburg O (1930) The Metabolism of Tumors. London: Arnold Constable.
[35]
Levine AJ, Puzio-Kuter AM (2010) The control of the metabolic switch in cancers by oncogenes and tumor suppressor genes. Science 330: 1340–1344.
[36]
Koffa MD, Casanova CM, Santarella R, Kocher T, Wilm M, et al. (2006) HURP is part of a Ran-dependent complex involved in spindle formation. Curr Biol 16: 743–754.
[37]
Wong J, Fang G (2006) HURP controls spindle dynamics to promote proper interkinetochore tension and efficient kinetochore capture. J Cell Biol 173: 879–891.
[38]
Wong J, Lerrigo R, Jang CY, Fang G (2008) Aurora A regulates the activity of HURP by controlling the accessibility of its microtubule-binding domain. Mol Biol Cell 19: 2083–2091.
[39]
Yu CT, Hsu JM, Lee YC, Tsou AP, Chou CK, et al. (2005) Phosphorylation and stabilization of HURP by Aurora-A: implication of HURP as a transforming target of Aurora-A. Mol Cell Biol 25: 5789–5800.
[40]
Bassal S, Nomura N, Venter D, Brand K, McKay MJ, et al. (2001) Characterization of a novel human cell-cycle-regulated homologue of Drosophila dlg1. Genomics 77: 5–7.
[41]
Lin F, Wang R, Shen JJ, Wang X, Gao P, et al. (2008) Knockdown of RCK/p54 expression by RNAi inhibits proliferation of human colorectal cancer cells in vitro and in vivo. Cancer Biol Ther 7: 1669–1676.
[42]
Nakagawa Y, Morikawa H, Hirata I, Shiozaki M, Matsumoto A, et al. (1999) Overexpression of rck/p54, a DEAD box protein, in human colorectal tumours. Br J Cancer 80: 914–917.
[43]
Shin S, Rossow KL, Grande JP, Janknecht R (2007) Involvement of RNA helicases p68 and p72 in colon cancer. Cancer Res 67: 7572–7578.
[44]
Foo RS, Nam YJ, Ostreicher MJ, Metzl MD, Whelan RS, et al. (2007) Regulation of p53 tetramerization and nuclear export by ARC. Proc Natl Acad Sci U S A 104: 20826–20831.
[45]
Heikaus S, Kempf T, Mahotka C, Gabbert HE, Ramp U (2008) Caspase-8 and its inhibitors in RCCs in vivo: the prominent role of ARC. Apoptosis 13: 938–949.
[46]
Mercier I, Vuolo M, Jasmin JF, Medina CM, Williams M, et al. (2008) ARC (apoptosis repressor with caspase recruitment domain) is a novel marker of human colon cancer. Cell Cycle 7: 1640–1647.
[47]
Wu L, Nam YJ, Kung G, Crow MT, Kitsis RN (2010) Induction of the apoptosis inhibitor ARC by Ras in human cancers. J Biol Chem 285: 19235–19245.
[48]
Carvajal-Carmona LG, Cazier JB, Jones AM, Howarth K, Broderick P, et al. (2011) Fine-mapping of colorectal cancer susceptibility loci at 8q23.3, 16q22.1 and 19q13.11: refinement of association signals and use of in silico analysis to suggest functional variation and unexpected candidate target genes. Hum Mol Genet 20: 2879–2888.
[49]
Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, et al. (2007) Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science 315: 848–853.
[50]
Prokunina-Olsson L, Hall JL (2009) No effect of cancer-associated SNP rs6983267 in the 8q24 region on co-expression of MYC and TCF7L2 in normal colon tissue. Mol Cancer 8: 96.
