Insights into Pancreatic Cancer Etiology from Pathway Analysis of Genome-Wide Association Study Data

Background Pancreatic cancer is the fourth leading cause of cancer death in the U.S. and the etiology of this highly lethal disease has not been well defined. To identify genetic susceptibility factors for pancreatic cancer, we conducted pathway analysis of genome-wide association study (GWAS) data in 3,141 pancreatic cancer patients and 3,367 controls with European ancestry. Methods Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways identified from the Kyoto Encyclopedia of Genes and Genomes database. We used the logistic kernel machine (LKM) test to identify major contributing genes to each pathway. We conducted functional enrichment analysis of the most significant genes (P<0.01) using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). Results Two pathways were significantly associated with risk of pancreatic cancer after adjusting for multiple comparisons (P<0.00025) and in replication testing: neuroactive ligand-receptor interaction, (Ps<0.00002), and the olfactory transduction pathway (P = 0.0001). LKM test identified four genes that were significantly associated with risk of pancreatic cancer after Bonferroni correction (P<1×10−5): ABO, HNF1A, OR13C4, and SHH. Functional enrichment analysis using DAVID consistently found the G protein-coupled receptor signaling pathway (including both neuroactive ligand-receptor interaction and olfactory transduction pathways) to be the most significant pathway for pancreatic cancer risk in this study population. Conclusion These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer.

[1]  Paul Brennan,et al.  Comparison of Pathway Analysis Approaches Using Lung Cancer GWAS Data Sets , 2012, PloS one.

[2]  M. Nei,et al.  Evolution of olfactory receptor genes in the human genome , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[3]  G. Parmigiani,et al.  Core Signaling Pathways in Human Pancreatic Cancers Revealed by Global Genomic Analyses , 2008, Science.

[4]  S. Dangi‐Garimella,et al.  Crosstalk between Mast Cells and Pancreatic Cancer Cells Contributes to Pancreatic Tumor Progression , 2010, Clinical Cancer Research.

[5]  L. Miller,et al.  A novel secretin receptor splice variant potentially useful for early diagnosis of pancreatic carcinoma. , 2007, Gastroenterology.

[6]  Lin S. Chen,et al.  Insights into colon cancer etiology via a regularized approach to gene set analysis of GWAS data. , 2010, American journal of human genetics.

[7]  Geoffrey S. Tobias,et al.  Genome-wide association study identifies variants in the ABO locus associated with susceptibility to pancreatic cancer , 2009, Nature Genetics.

[8]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.

[9]  Stephen W. Hartley,et al.  Fetal hemoglobin in sickle cell anemia: genome-wide association studies suggest a regulatory region in the 5' olfactory receptor gene cluster. , 2010, Blood.

[10]  George C Patton,et al.  Variation in the gene coding for the M5 Muscarinic receptor (CHRM5) influences cigarette dose but is not associated with dependence to drugs of addiction: evidence from a prospective population based cohort study of young adults , 2007, BMC Genetics.

[11]  Zhaohui S. Qin,et al.  A second generation human haplotype map of over 3.1 million SNPs , 2007, Nature.

[12]  Deanne M. Taylor,et al.  Powerful SNP-set analysis for case-control genome-wide association studies. , 2010, American journal of human genetics.

[13]  K. Kaestner,et al.  The diabetes gene Pdx1 regulates the transcriptional network of pancreatic endocrine progenitor cells in mice. , 2009, The Journal of clinical investigation.

[14]  Hanns Hatt,et al.  Activation of an Olfactory Receptor Inhibits Proliferation of Prostate Cancer Cells* , 2009, The Journal of Biological Chemistry.

[15]  R. Collins,et al.  Common variants at 30 loci contribute to polygenic dyslipidemia , 2009, Nature Genetics.

[16]  M. Maggiolini,et al.  G protein-coupled receptors: novel targets for drug discovery in cancer , 2010, Nature Reviews Drug Discovery.

[17]  Xihong Lin,et al.  Hypothesis testing in semiparametric additive mixed models. , 2003, Biostatistics.

[18]  Zhongming Zhao,et al.  A bias-reducing pathway enrichment analysis of genome-wide association data confirmed association of the MHC region with schizophrenia , 2011, Journal of Medical Genetics.

[19]  A. Whittemore,et al.  A genome-wide association study identifies susceptibility loci for ovarian cancer at 2q31 and 8q24 , 2010, Nature Genetics.

[20]  Kai Wang,et al.  Pathway-based approaches for analysis of genomewide association studies. , 2007, American journal of human genetics.

[21]  Karin E. Borgmann-Winter,et al.  Scents and nonsense: olfactory dysfunction in schizophrenia. , 2009, Schizophrenia bulletin.

[22]  R. F. Luco,et al.  Distinct Roles of HNF1 Β , HNF1 α , and HNF4 α in Regulating Pancreas Development, Β -Cell Function and Growth , 2007 .

[23]  D. Reich,et al.  Principal components analysis corrects for stratification in genome-wide association studies , 2006, Nature Genetics.

[24]  M. Hebrok,et al.  KRAS, Hedgehog, Wnt and the twisted developmental biology of pancreatic ductal adenocarcinoma , 2010, Nature Reviews Cancer.

[25]  M. Bohlooly-y,et al.  Phenotypic screening of hepatocyte nuclear factor (HNF) 4-gamma receptor knockout mice. , 2006, Biochemical and biophysical research communications.

[26]  Simon C. Potter,et al.  Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls , 2007, Nature.

[27]  P. Elliott,et al.  Meta-Analysis of Genome-Wide Association Studies in >80 000 Subjects Identifies Multiple Loci for C-Reactive Protein Levels , 2011, Circulation.

[28]  Tanya M. Teslovich,et al.  Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index , 2010 .

[29]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[30]  M. Xiong,et al.  Genome-wide gene and pathway analysis , 2010, European Journal of Human Genetics.

[31]  H. Hakonarson,et al.  Analysing biological pathways in genome-wide association studies , 2010, Nature Reviews Genetics.

[32]  P. Camacho.  Clinical Endocrinology and Metabolism , 2011 .

[33]  Wei Zheng,et al.  A genome-wide association study identifies pancreatic cancer susceptibility loci on chromosomes 13q22.1, 1q32.1 and 5p15.33 , 2010, Nature Genetics.

[34]  Susumu Goto,et al.  KEGG for representation and analysis of molecular networks involving diseases and drugs , 2009, Nucleic Acids Res..

[35]  Chronic Disease Division.  Cancer facts and figures , 2010 .

[36]  Manuel A. R. Ferreira,et al.  Practical aspects of imputation-driven meta-analysis of genome-wide association studies. , 2008, Human molecular genetics.

[37]  J. Barrett,et al.  Pathway-Based Analysis of a Melanoma Genome-Wide Association Study: Analysis of Genes Related to Tumour-Immunosuppression , 2011, PloS one.

[38]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[39]  D. Lancet,et al.  Widespread ectopic expression of olfactory receptor genes , 2006, BMC Genomics.

[40]  U. G. Dailey.  Cancer,Facts and Figures about. , 2022, Journal of the National Medical Association.

[41]  Manuel A. R. Ferreira,et al.  Gene ontology analysis of GWA study data sets provides insights into the biology of bipolar disorder. , 2009, American journal of human genetics.

[42]  J. Rehfeld.  Clinical endocrinology and metabolism. Cholecystokinin. , 2004, Best practice & research. Clinical endocrinology & metabolism.

[43]  R. Gibson,et al.  EDNRA variants associate with smooth muscle mRNA levels, cell proliferation rates, and cystic fibrosis pulmonary disease severity. , 2010, Physiological genomics.

[44]  Wei Wang,et al.  Genomics Meets Glycomics—The First GWAS Study of Human N-Glycome Identifies HNF1α as a Master Regulator of Plasma Protein Fucosylation , 2010, PLoS genetics.

[45]  M. Bohlooly-y,et al.  Phenotypic screening of hepatocyte nuclear factor (HNF) 4-γ receptor knockout mice , 2006 .

[46]  J. Habener,et al.  Neurogenin3: A master regulator of pancreatic islet differentiation and regeneration , 2009, Islets.

[47]  Terrence S. Furey,et al.  The UCSC Table Browser data retrieval tool , 2004, Nucleic Acids Res..

[48]  Ayellet V. Segrè,et al.  Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis , 2010, Nature Genetics.

[49]  E. Schadt.  Molecular networks as sensors and drivers of common human diseases , 2009, Nature.

[50]  Yun Li,et al.  METAL: fast and efficient meta-analysis of genomewide association scans , 2010, Bioinform..

[51]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[52]  Geoffrey S. Tobias,et al.  Pathway analysis of genome-wide association study data highlights pancreatic development genes as susceptibility factors for pancreatic cancer. , 2012, Carcinogenesis.

[53]  R. F. Luco,et al.  Distinct roles of HNF1beta, HNF1alpha, and HNF4alpha in regulating pancreas development, beta-cell function and growth. , 2007, Endocrine development.

[54]  N. Funel,et al.  Inflammatory cells contribute to the generation of an angiogenic phenotype in pancreatic ductal adenocarcinoma , 2004, Journal of Clinical Pathology.