Comprehensive Functional Annotation of 77 Prostate Cancer Risk Loci

Genome-wide association studies (GWAS) have revolutionized the field of cancer genetics, but the causal links between increased genetic risk and onset/progression of disease processes remain to be identified. Here we report the first step in such an endeavor for prostate cancer. We provide a comprehensive annotation of the 77 known risk loci, based upon highly correlated variants in biologically relevant chromatin annotations— we identified 727 such potentially functional SNPs. We also provide a detailed account of possible protein disruption, microRNA target sequence disruption and regulatory response element disruption of all correlated SNPs at . 88% of the 727 SNPs fall within putative enhancers, and many alter critical residues in the response elements of transcription factors known to be involved in prostate biology. We define as risk enhancers those regions with enhancer chromatin biofeatures in prostate-derived cell lines with prostate-cancer correlated SNPs. To aid the identification of these enhancers, we performed genomewide ChIP-seq for H3K27-acetylation, a mark of actively engaged enhancers, as well as the transcription factor TCF7L2. We analyzed in depth three variants in risk enhancers, two of which show significantly altered androgen sensitivity in LNCaP cells. This includes rs4907792, that is in linkage disequilibrium () with an eQTL for NUDT11 (on the X chromosome) in prostate tissue, and rs10486567, the index SNP in intron 3 of the JAZF1 gene on chromosome 7. Rs4907792 is within a critical residue of a strong consensus androgen response element that is interrupted in the protective allele, resulting in a 56% decrease in its androgen sensitivity, whereas rs10486567 affects both NKX3-1 and FOXA-AR motifs where the risk allele results in a 39% increase in basal activity and a 28% fold-increase in androgen stimulated enhancer activity. Identification of such enhancer variants and their potential target genes represents a preliminary step in connecting risk to disease process.

[1]  J. Carpten,et al.  Two genome-wide association studies of aggressive prostate cancer implicate putative prostate tumor suppressor gene DAB2IP. , 2007, Journal of the National Cancer Institute.

[2]  J. Trapman Molecular mechanisms of prostate cancer. , 2001, European journal of cancer.

[3]  P. Bork,et al.  A method and server for predicting damaging missense mutations , 2010, Nature Methods.

[4]  Irene K. Moore,et al.  A genomic code for nucleosome positioning , 2006, Nature.

[5]  S. Salzberg,et al.  Computational identification of developmental enhancers: conservation and function of transcription factor binding-site clusters in Drosophila melanogaster and Drosophila pseudoobscura , 2004, Genome Biology.

[6]  J. Carpten,et al.  Association of prostate cancer risk with snps in regions containing androgen receptor binding sites captured by ChIP‐On‐chip analyses , 2012, The Prostate.

[7]  Thomas Whitington,et al.  Transcription Factor Binding in Human Cells Occurs in Dense Clusters Formed around Cohesin Anchor Sites , 2013, Cell.

[8]  J. Carpten,et al.  Two Independent Prostate Cancer Risk–Associated Loci at 11q13 , 2009, Cancer Epidemiology Biomarkers & Prevention.

[9]  Bertram Ludäscher,et al.  Sole-Search: an integrated analysis program for peak detection and functional annotation using ChIP-seq data , 2009, Nucleic acids research.

[10]  M. Sadar,et al.  Cell lines used in prostate cancer research: a compendium of old and new lines--part 1. , 2005, The Journal of urology.

[11]  C. Jung,et al.  HOXB13 Induces Growth Suppression of Prostate Cancer Cells as a Repressor of Hormone-Activated Androgen Receptor Signaling , 2004, Cancer Research.

[12]  Kevin S. Xu,et al.  A Common Genetic Variant (97906C>A) of DAB2IP/AIP1 Is Associated with an Increased Risk and Early Onset of Lung Cancer in Chinese Males , 2011, PloS one.

[13]  Ali Amin Al Olama,et al.  Identification of seven new prostate cancer susceptibility loci through a genome-wide association study , 2009, Nature Genetics.

[14]  北村 聖 "The New England Journal of Medicine". , 1962, British medical journal.

[15]  Simon G. Coetzee,et al.  A rare variant, which destroys a FoxA1 site at 8q24, is associated with prostate cancer risk , 2013, Cell cycle.

[16]  P. Fearnhead,et al.  Genome-wide association study of prostate cancer identifies a second risk locus at 8q24 , 2007, Nature Genetics.

[17]  Peter A. Jones,et al.  Dynamic Nucleosome-Depleted Regions at Androgen Receptor Enhancers in the Absence of Ligand in Prostate Cancer Cells , 2011, Molecular and Cellular Biology.

[18]  M. Sadar,et al.  Cell lines used in prostate cancer research: a compendium of old and new lines--part 2. , 2005, The Journal of urology.

[19]  Gary K. Chen,et al.  Genome-wide association study of prostate cancer in men of African ancestry identifies a susceptibility locus at 17q21 , 2011, Nature Genetics.

[20]  C. Burge,et al.  Most mammalian mRNAs are conserved targets of microRNAs. , 2008, Genome research.

[21]  J. Carpten,et al.  A novel prostate cancer susceptibility locus at 19q13. , 2009, Cancer research.

[22]  A. Whittemore,et al.  HOXB13 is a susceptibility gene for prostate cancer: results from the International Consortium for Prostate Cancer Genetics (ICPCG) , 2012, Human Genetics.

[23]  J. Miller,et al.  Predicting the Functional Effect of Amino Acid Substitutions and Indels , 2012, PloS one.

[24]  Carol Moreno,et al.  Identifying multiple causative genes at a single GWAS locus , 2013, Genome research.

[25]  Yusuke Nakamura,et al.  Genome-wide association study identifies five new susceptibility loci for prostate cancer in the Japanese population , 2010, Nature Genetics.

[26]  Ali Amin Al Olama,et al.  Multiple newly identified loci associated with prostate cancer susceptibility , 2008, Nature Genetics.

[27]  Arcadi Navarro,et al.  High Trans-ethnic Replicability of GWAS Results Implies Common Causal Variants , 2013, PLoS genetics.

[28]  Noula Shembade,et al.  Essential role for TAX1BP1 in the termination of TNF‐α‐, IL‐1‐ and LPS‐mediated NF‐κB and JNK signaling , 2007 .

[29]  D. Gudbjartsson,et al.  Two variants on chromosome 17 confer prostate cancer risk, and the one in TCF2 protects against type 2 diabetes , 2007, Nature Genetics.

[30]  Simon G. Coetzee,et al.  Comprehensive Functional Annotation of Seventy-One Breast Cancer Risk Loci , 2013, PloS one.

[31]  C. Glass,et al.  Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. , 2010, Molecular cell.

[32]  M. Thun,et al.  Genome-wide association study identifies new prostate cancer susceptibility loci. , 2011, Human molecular genetics.

[33]  Benjamin P. Berman,et al.  FunciSNP: an R/bioconductor tool integrating functional non-coding data sets with genetic association studies to identify candidate regulatory SNPs , 2012, Nucleic acids research.

[34]  Bo Dai,et al.  Genome-wide association study in Chinese men identifies two new prostate cancer risk loci at 9q31.2 and 19q13.4 , 2012, Nature Genetics.

[35]  D M Crothers,et al.  Artificial nucleosome positioning sequences. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[36]  G. Murphy,et al.  LNCaP model of human prostatic carcinoma. , 1983, Cancer research.

[37]  K. Cooney,et al.  The Prostate 71 : 225 ^ 231 ( 2011 ) Evidence for anAssociationBetweenProstate Cancer andChromosome 8 q 24 and 10 q 11 GeneticVariants inAfricanAmericanMen : The FlintMen ’ sHealth Study , 2011 .

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

[39]  D. Reich,et al.  Functional Enhancers at the Gene-Poor 8q24 Cancer-Linked Locus , 2009, PLoS genetics.

[40]  M. Kattan,et al.  New variants at 10q26 and 15q21 are associated with aggressive prostate cancer in a genome-wide association study from a prostate biopsy screening cohort , 2011, Cancer biology & therapy.

[41]  B. Rost,et al.  SNAP: predict effect of non-synonymous polymorphisms on function , 2007, Nucleic acids research.

[42]  W. Hahn,et al.  Genetic and functional analyses implicate the NUDT11, HNF1B, and SLC22A3 genes in prostate cancer pathogenesis , 2012, Proceedings of the National Academy of Sciences.

[43]  I. Mills,et al.  The androgen receptor induces a distinct transcriptional program in castration-resistant prostate cancer in man. , 2013, Cancer cell.

[44]  Jin Woo Kim,et al.  Sequence variants at 22q13 are associated with prostate cancer risk. , 2009, Cancer research.

[45]  Beatrice S Knudsen,et al.  Mechanisms of prostate cancer initiation and progression. , 2010, Advances in cancer research.

[46]  Meng Li,et al.  Association of a Genetic Variation in a miR-191 Binding Site in MDM4 with Risk of Esophageal Squamous Cell Carcinoma , 2013, PloS one.

[47]  E. Dermitzakis,et al.  Passive and active DNA methylation and the interplay with genetic variation in gene regulation , 2013, eLife.

[48]  David B. Goldstein,et al.  Rare Variants Create Synthetic Genome-Wide Associations , 2010, PLoS biology.

[49]  Teri A Manolio,et al.  Genomewide association studies and assessment of the risk of disease. , 2010, The New England journal of medicine.

[50]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[51]  Daniel Levy,et al.  A genome-wide association study of breast and prostate cancer in the NHLBI's Framingham Heart Study , 2007, BMC Medical Genetics.

[52]  M. Taplin Drug Insight: role of the androgen receptor in the development and progression of prostate cancer , 2007, Nature Clinical Practice Oncology.

[53]  Kari Stefansson,et al.  Genome-wide association and replication studies identify four variants associated with prostate cancer susceptibility , 2009, Nature Genetics.

[54]  Peter Kraft,et al.  Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping array , 2013, Nature Genetics.

[55]  M. Vingron,et al.  Sequence-dependent nucleosome positioning. , 2009, Journal of molecular biology.

[56]  Data production leads,et al.  An integrated encyclopedia of DNA elements in the human genome , 2012 .

[57]  G. Rubin,et al.  Exploiting transcription factor binding site clustering to identify cis-regulatory modules involved in pattern formation in the Drosophila genome , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[58]  B. Bjerkehagen,et al.  Consistent rearrangement of chromosomal band 6p21 with generation of fusion genes JAZF1/PHF1 and EPC1/PHF1 in endometrial stromal sarcoma. , 2006, Cancer research.

[59]  Clifford A. Meyer,et al.  Androgen Receptor Regulates a Distinct Transcription Program in Androgen-Independent Prostate Cancer , 2009, Cell.

[60]  M. Thun,et al.  Refining the Prostate Cancer Genetic Association within the JAZF1 Gene on Chromosome 7p15.2 , 2010, Cancer Epidemiology, Biomarkers & Prevention.

[61]  Shane J. Neph,et al.  Systematic Localization of Common Disease-Associated Variation in Regulatory DNA , 2012, Science.

[62]  Kenny Q. Ye,et al.  An integrated map of genetic variation from 1,092 human genomes , 2012, Nature.

[63]  J. Carpten,et al.  Evidence for two independent prostate cancer risk–associated loci in the HNF1B gene at 17q12 , 2008, Nature Genetics.

[64]  ENCODEConsortium,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[65]  J. Carpten,et al.  Germline mutations in HOXB13 and prostate-cancer risk. , 2012, The New England journal of medicine.

[66]  Nathan C. Sheffield,et al.  Chromatin accessibility reveals insights into androgen receptor activation and transcriptional specificity , 2012, Genome Biology.

[67]  Liqing Zhou,et al.  Functional MDM4 rs4245739 genetic variant, alone and in combination with P53 Arg72Pro polymorphism, contributes to breast cancer susceptibility , 2013, Breast Cancer Research and Treatment.

[68]  Kevin M. Bradley,et al.  Common sequence variants on 2p15 and Xp11.22 confer susceptibility to prostate cancer , 2008, Nature Genetics.

[69]  R. Kittles,et al.  Prostate Cancer Susceptibility Loci Identified on Chromosome 12 in African Americans , 2011, PloS one.

[70]  Kari Stefansson,et al.  A study based on whole-genome sequencing yields a rare variant at 8q24 associated with prostate cancer , 2012, Nature Genetics.

[71]  C. Jung,et al.  HOXB13 homeodomain protein suppresses the growth of prostate cancer cells by the negative regulation of T-cell factor 4. , 2004, Cancer research.

[72]  Ilir Agalliu,et al.  The Prostate 68 : 740 ^ 747 ( 2008 ) Evaluationof aVariant intheTranscription Factor 7-Like 2 ( TCF 7 L 2 ) Gene andProstate CancerRisk in a Population-Based Study , 2008 .

[73]  M. Capecchi,et al.  Hoxb13 is required for normal differentiation and secretory function of the ventral prostate , 2003, Development.

[74]  Gene W. Yeo,et al.  Determination of tag density required for digital transcriptome analysis: Application to an androgen-sensitive prostate cancer model , 2008, Proceedings of the National Academy of Sciences.

[75]  Suzanne Chambers,et al.  Seven prostate cancer susceptibility loci identified by a multi-stage genome-wide association study , 2011, Nature Genetics.

[76]  D. Gudbjartsson,et al.  Genome-wide association study identifies a second prostate cancer susceptibility variant at 8q24 , 2007, Nature Genetics.

[77]  P. Farnham,et al.  Using ChIP-seq technology to generate high-resolution profiles of histone modifications. , 2011, Methods in molecular biology.

[78]  W. Willett,et al.  Multiple loci identified in a genome-wide association study of prostate cancer , 2008, Nature Genetics.

[79]  Eurie L. Hong,et al.  Annotation of functional variation in personal genomes using RegulomeDB , 2012, Genome research.

[80]  Andrew Gelman,et al.  The No-U-turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo , 2011, J. Mach. Learn. Res..

[81]  S. Henikoff,et al.  Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm , 2009, Nature Protocols.

[82]  J. Sklar,et al.  Frequent fusion of the JAZF1 and JJAZ1 genes in endometrial stromal tumors , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[83]  William Stafford Noble,et al.  Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors , 2012, Genome research.

[84]  Daniel L. Lakeland,et al.  Affinity Density: a novel genomic approach to the identification of transcription factor regulatory targets , 2009, Bioinform..

[85]  L. Jia,et al.  Persistent androgen receptor-mediated transcription in castration-resistant prostate cancer under androgen-deprived conditions , 2012, Nucleic acids research.

[86]  Wyeth W. Wasserman,et al.  Identification of functional clusters of transcription factor binding motifs in genome sequences: the MSCAN algorithm , 2003, ISMB.

[87]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

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

[89]  Samuel P. Dickson,et al.  Interpretation of association signals and identification of causal variants from genome-wide association studies. , 2010, American journal of human genetics.