Pathway analysis of genome-wide association studies (GWAS) offer a unique opportunity to collectively evaluate genetic variants with effects that are too small to be detected individually. We applied a pathway analysis to a bladder cancer GWAS containing data from 3,532 cases and 5,120 controls of European background (n = 5 studies). Thirteen hundred and ninety-nine pathways were drawn from five publicly available resources (Biocarta, Kegg, NCI-PID, HumanCyc, and Reactome), and we constructed 22 additional candidate pathways previously hypothesized to be related to bladder cancer. In total, 1421 pathways, 5647 genes and ~90,000 SNPs were included in our study. Logistic regression model adjusting for age, sex, study, DNA source, and smoking status was used to assess the marginal trend effect of SNPs on bladder cancer risk. Two complementary pathway-based methods (gene-set enrichment analysis [GSEA], and adapted rank-truncated product [ARTP]) were used to assess the enrichment of association signals within each pathway. Eighteen pathways were detected by either GSEA or ARTP at P≤0.01. To minimize false positives, we used the I2 statistic to identify SNPs displaying heterogeneous effects across the five studies. After removing these SNPs, seven pathways (‘Aromatic amine metabolism’ [PGSEA = 0.0100, PARTP = 0.0020], ‘NAD biosynthesis’ [PGSEA = 0.0018, PARTP = 0.0086], ‘NAD salvage’ [PARTP = 0.0068], ‘Clathrin derived vesicle budding’ [PARTP = 0.0018], ‘Lysosome vesicle biogenesis’ [PGSEA = 0.0023, PARTP<0.00012], ’Retrograde neurotrophin signaling’ [PGSEA = 0.00840], and ‘Mitotic metaphase/anaphase transition’ [PGSEA = 0.0040]) remained. These pathways seem to belong to three fundamental cellular processes (metabolic detoxification, mitosis, and clathrin-mediated vesicles). Identification of the aromatic amine metabolism pathway provides support for the ability of this approach to identify pathways with established relevance to bladder carcinogenesis.
References
[1]
Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, et al. (2009) Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A 106: 9362–9367.
[2]
Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, et al. (2009) Finding the missing heritability of complex diseases. Nature 461: 747–753.
[3]
Eichler EE, Flint J, Gibson G, Kong A, Leal SM, et al. (2010) Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet 11: 446–450.
[4]
Park JH, Wacholder S, Gail MH, Peters U, Jacobs KB, et al. (2010) Estimation of effect size distribution from genome-wide association studies and implications for future discoveries. Nat Genet 42: 570–575.
[5]
Elbers CC, van Eijk KR, Franke L, Mulder F, van der Schouw YT, et al. (2009) Using genome-wide pathway analysis to unravel the etiology of complex diseases. Genet Epidemiol.
[6]
Wang K, Li M, Hakonarson H (2010) Analysing biological pathways in genome-wide association studies. Nat Rev Genet 11: 843–854.
[7]
Baranzini SE, Galwey NW, Wang J, Khankhanian P, Lindberg R, et al. (2009) Pathway and network-based analysis of genome-wide association studies in multiple sclerosis. Hum Mol Genet.
[8]
Perry JR, McCarthy MI, Hattersley AT, Zeggini E, Weedon MN, et al. (2009) Interrogating Type 2 Diabetes Genome-Wide Association Data Using a Biological Pathway-Based Approach. Diabetes.
[9]
Zhong H, Yang X, Kaplan LM, Molony C, Schadt EE (2010) Integrating pathway analysis and genetics of gene expression for genome-wide association studies. Am J Hum Genet 86: 581–591.
[10]
Wang K, Zhang H, Kugathasan S, Annese V, Bradfield JP, et al. (2009) Diverse genome-wide association studies associate the IL12/IL23 pathway with Crohn Disease. Am J Hum Genet 84: 399–405.
[11]
Chen X, Wang L, Hu B, Guo M, Barnard J, et al. (2010) Pathway-based analysis for genome-wide association studies using supervised principal components. Genet Epidemiol 34: 716–724.
[12]
Wang K, Li M, Bucan M (2007) Pathway-Based Approaches for Analysis of Genomewide Association Studies. Am J Hum Genet 81: 1278–1283.
[13]
Lesnick TG, Papapetropoulos S, Mash DC, Ffrench-Mullen J, Shehadeh L, et al. (2007) A genomic pathway approach to a complex disease: axon guidance and Parkinson disease. PLoS Genet 3: e98.
[14]
Chen LS, Hutter CM, Potter JD, Liu Y, Prentice RL, et al. (2010) Insights into colon cancer etiology via a regularized approach to gene set analysis of GWAS data. Am J Hum Genet 86: 860–871.
[15]
Menashe I, Maeder D, Garcia-Closas M, Figueroa JD, Bhattacharjee S, et al. (2010) Pathway analysis of breast cancer genome-wide association study highlights three pathways and one canonical signaling cascade. Cancer Res 70: 4453–4459.
[16]
Silverman DT, Devesa SS, Moore LE, Rothman N (2006) Bladder cancer. In: Schottenfeld D, Fraumeni JF Jr, editors. Cancer Epidemiology and Prevention. 3 ed. New York: Oxford University Press. pp. 1101–1127.
[17]
Silverman DT, Hartge P, Morrison AS, Devesa SS (1992) Epidemiology of bladder cancer. Hematol Oncol Clin North Am 6: 1–30.
[18]
Vineis P, Pirastu R (1997) Aromatic amines and cancer. Cancer Causes Control 8: 346–355.
[19]
Talaska G (2003) Aromatic amines and human urinary bladder cancer: exposure sources and epidemiology. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev 21: 29–43.
[20]
Garcia-Closas M, Malats N, Silverman D, Dosemeci M, Kogevinas M, et al. (2005) NAT2 slow acetylation, GSTM1 null genotype, and risk of bladder cancer: results from the Spanish Bladder Cancer Study and meta-analyses. Lancet 366: 649–659.
[21]
Moore LE, Baris DR, Figueroa JD, Garcia-Closas M, Karagas MR, et al. (2011) GSTM1 null and NAT2 slow acetylation genotypes, smoking intensity and bladder cancer risk: results from the New England bladder cancer study and NAT2 meta-analysis. Carcinogenesis 32: 182–189.
[22]
Rothman N, Garcia-Closas M, Chatterjee N, Malats N, Wu X, et al. (2010) A multi-stage genome-wide association study of bladder cancer identifies multiple susceptibility loci. Nat Genet 42: 978–984.
[23]
Kiemeney LA, Sulem P, Besenbacher S, Vermeulen SH, Sigurdsson A, et al. (2010) A sequence variant at 4p16.3 confers susceptibility to urinary bladder cancer. Nat Genet 42: 415–419.
[24]
Kiemeney LA, Thorlacius S, Sulem P, Geller F, Aben KK, et al. (2008) Sequence variant on 8q24 confers susceptibility to urinary bladder cancer. Nat Genet 40: 1307–1312.
[25]
Wu X, Ye Y, Kiemeney LA, Sulem P, Rafnar T, et al. (2009) Genetic variation in the prostate stem cell antigen gene PSCA confers susceptibility to urinary bladder cancer. Nat Genet 41: 991–995.
[26]
BioCarta website. Available: http://www.biocarta.com/genes/allpathway？s.asp. Accessed 2011 Mar 10.
[27]
KEGG website. Available: http://www.genome.jp/kegg/pathway.html. Accessed 2011 Mar 10.
[28]
Pathway Interaction Database. Available: http://pid.nci.nih.gov. Accessed 2011 Mar 10.
[29]
Reactome website. Available: http://www.reactome.org. Accessed 2011 Mar 10.
[30]
HumanCyc website. Available: http://humancyc.org. Accessed 2011 Mar 10.
[31]
Figueroa JD, Malats N, Rothman N, Real FX, Silverman D, et al. (2007) Evaluation of genetic variation in the double-strand break repair pathway and bladder cancer risk. Carcinogenesis 28: 1788–1793.
[32]
Figueroa JD, Malats N, Real FX, Silverman D, Kogevinas M, et al. (2007) Genetic variation in the base excision repair pathway and bladder cancer risk. Hum Genet 121: 233–242.
[33]
Figueroa JD, Garcia-Closas M, Rothman N (2010) Case studies: Cumulative assessment of the role of human genome variation in specific diseases - bladder cancer. In: Khoury M, Bedrosian S, Gwinn M, Higgins J, Ioannidis J, et al., editors. Human Genome Epidemiology, 2nd Edition Building the evidence for using genetic information to improve health and prevent disease. Oxford University Press.
[34]
Skipper PL, Tannenbaum SR (1994) Molecular dosimetry of aromatic amines in human populations. Environ Health Perspect 102: Suppl 617–21.
[35]
Skipper PL, Kim MY, Sun HL, Wogan GN, Tannenbaum SR (2010) Monocyclic aromatic amines as potential human carcinogens: old is new again. Carcinogenesis 31: 50–58.
[36]
Ingenuity. website. Available: http://www.ingenuity.com/. Accessed 2010 Feb 24.
[37]
The Gene Ontology website. Available: http://www.geneontology.org/. Accessed 2010 Feb 24.
[38]
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.
[39]
Yu K, Li Q, Bergen AW, Pfeiffer RM, Rosenberg PS, et al. (2009) Pathway analysis by adaptive combination of P-values. Genet Epidemiol 33: 700–709.
[40]
Dudbridge F, Koeleman BP (2003) Rank truncated product of P-values, with application to genomewide association scans. Genet Epidemiol 25: 360–366.
[41]
Hoh J, Wille A, Ott J (2001) Trimming, weighting, and grouping SNPs in human case-control association studies. Genome Res 11: 2115–2119.
[42]
Higgins JP, Thompson SG (2002) Quantifying heterogeneity in a meta-analysis. Stat Med 21: 1539–1558.
[43]
Magni G, Orsomando G, Raffelli N, Ruggieri S (2008) Enzymology of mammalian NAD metabolism in health and disease. Front Biosci 13: 6135–6154.
[44]
Burkle A (2005) Poly(ADP-ribose). The most elaborate metabolite of NAD+. FEBS J 272: 4576–4589.
[45]
Hein DW (2002) Molecular genetics and function of NAT1 and NAT2: role in aromatic amine metabolism and carcinogenesis. Mutat Res 506–507: 65–77.
[46]
Royle SJ (2006) The cellular functions of clathrin. Cell Mol Life Sci 63: 1823–1832.
[47]
Kirisits A, Pils D, Krainer M (2007) Epidermal growth factor receptor degradation: an alternative view of oncogenic pathways. Int J Biochem Cell Biol 39: 2173–2182.
[48]
Royle SJ, Bright NA, Lagnado L (2005) Clathrin is required for the function of the mitotic spindle. Nature 434: 1152–1157.
[49]
Li X, Chen J, Hu X, Huang Y, Li Z, et al. (2011) Comparative mRNA and microRNA expression profiling of three genitourinary cancers reveals common hallmarks and cancer-specific molecular events. PLoS One 6: e22570.
