PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Nat Genet. Author manuscript; available in PMC Mar 15, 2013.
Published in final edited form as:
Published online Feb 10, 2008. doi:  10.1038/ng.89
PMCID: PMC3598012
NIHMSID: NIHMS371835
Common sequence variants on 2p15 and Xp11.22 confer susceptibility to prostate cancer
Julius Gudmundsson,1,26 Patrick Sulem,1,26 Thorunn Rafnar,1 Jon T Bergthorsson,1 Andrei Manolescu,1 Daniel Gudbjartsson,1 Bjarni A Agnarsson,2 Asgeir Sigurdsson,1 Kristrun R Benediktsdottir,2 Thorarinn Blondal,1 Margret Jakobsdottir,1 Simon N Stacey,1 Jelena Kostic,1 Kari T Kristinsson,1 Birgitta Birgisdottir,1 Shyamali Ghosh,1 Droplaug N Magnusdottir,1 Steinunn Thorlacius,1 Gudmar Thorleifsson,1 S Lilly Zheng,3 Jielin Sun,3 Bao-Li Chang,3 J Bradford Elmore,4 Joan P Breyer,4 Kate M McReynolds,4 Kevin M Bradley,4 Brian L Yaspan,4 Fredrik Wiklund,5 Par Stattin,6 Sara Lindström,5 Hans-Olov Adami,5,7 Shannon K McDonnell,8 Daniel J Schaid,8 Julie M Cunningham,9 Liang Wang,9 James R Cerhan,8 Jennifer L St Sauver,8 Sara D Isaacs,10 Kathleen E Wiley,10 Alan W Partin,10 Patrick C Walsh,10 Sonia Polo,11 Manuel Ruiz-Echarri,12 Sebastian Navarrete,13 Fernando Fuertes,13 Berta Saez,14,15 Javier Godino,14 Philip C Weijerman,16 Dorine W Swinkels,17 Katja K Aben,18,19 J Alfred Witjes,20 Brian K Suarez,21 Brian T Helfand,22 Michael L Frigge,1 Kristleifur Kristjansson,1 Carole Ober,23 Eirikur Jonsson,24 Gudmundur V Einarsson,24 Jianfeng Xu,3 Henrik Gronberg,5 Jeffrey R Smith,4,25 Stephen N Thibodeau,9 William B Isaacs,10 William J Catalona,22 Jose I Mayordomo,12,14,15 Lambertus A Kiemeney,18,19,20 Rosa B Barkardottir,2 Jeffrey R Gulcher,1 Unnur Thorsteinsdottir,1 Augustine Kong,1 and Kari Stefansson1
1deCODE genetics, 101 Reykjavik, Iceland
2Department of Pathology, Landspitali-University Hospital, 101 Reykjavik, Iceland
3Center for Human Genomics, Wake Forest University School of Medicine, Winston-Salem, North Carolina 27157, USA
4Department of Medicine, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, Tennessee 37232, USA
5Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, S-17177 Stockholm, Sweden
6Department of Urology, Umea University Hospital, S-90185 Umea, Sweden
7Department of Epidemiology, Harvard School of Public Health, Boston, Massachusetts 02115, USA
8Department of Health Sciences Research, Mayo Clinic, Rochester, Minnesota 55905, USA
9Department of Laboratory Medicine and Pathology, Mayo Clinic, Rochester, Minnesota 55905, USA
10Brady Urological Institute, Johns Hopkins Medical Institution, Baltimore, Maryland 21209, USA
11Division of Radiation Oncology, Quiron Hospital, 50009 Zaragoza, Spain
12Division of Medical Oncology, Lozano Blesa University Hospital, University of Zaragoza, 50009 Zaragoza, Spain
13Division of Radiation Oncology, Lozano Blesa University Hospital, University of Zaragoza, 50009 Zaragoza, Spain
14Health Sciences Institute, 50009 Zaragoza, Spain
15Nanotechnology Institute of Aragon, 50009 Zaragoza, Spain
16Department of Urology, Rijnstate Hospital, NL-6800 TA Arnhem, The Netherlands
17Department of Clinical Chemistry, Radboud University Nijmegen Medical Center, 6500 HB Nijmegen, The Netherlands
18Comprehensive Cancer Center East, 6501 BG Nijmegen, The Netherlands
19Department of Epidemiology and Biostatistics, Radboud University Nijmegen Medical Center, 6500 HB Nijmegen, The Netherlands
20Department of Urology, Radboud University Nijmegen Medical Center, 6500 HB Nijmegen, The Netherlands
21Department of Psychiatry, Washington University School of Medicine, St Louis, Missouri 63110, USA
22Department of Urology, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA
23Department of Human Genetics, University of Chicago, Chicago, Illinois 60637, USA
24Department of Urology, Landspitali-University Hospital, 101 Reykjavik, Iceland
25Medical Research Service, VA Tennessee Valley Healthcare System, Nashville, Tennessee 37212, USA
Correspondence should be addressed to K.S. (kstefans/at/decode.is) or J.G. (julius/at/decode.is)
26These authors contributed equally to this work.
Abstract
We conducted a genome-wide SNP association study on prostate cancer on over 23,000 Icelanders, followed by a replication study including over 15,500 individuals from Europe and the United States. Two newly identified variants were shown to be associated with prostate cancer: rs5945572 on Xp11.22 and rs721048 on 2p15 (odds ratios (OR) = 1.23 and 1.15; P = 3.9 × 10−13 and 7.7 × 10−9, respectively). The 2p15 variant shows a significantly stronger association with more aggressive, rather than less aggressive, forms of the disease.
Recently, we presented results from a genome-wide association (GWA) analysis of prostate cancer in which we used the Illumina Hap300 SNP chip to genotype Icelandic samples1,2. In subsequent efforts to identify additional variants, we used two criteria to guide our follow-up analysis of the GWA data. First, several studies on prostate cancer risk in relatives of affected men have reported that brothers are at greater risk compared to fathers3,4, thereby implicating an X-linked or recessive mode of inheritance. Second, identifying genetic factors that preferentially associate to aggressive prostate cancer is of utmost importance for improving treatment selection for early stage disease. Hence, we decided to fast track the variant from our GWA analysis (Supplementary Methods online) showing the strongest association among the Illumina Hap300 SNPs located on the X chromosome, as well as a sequence variant on 2p15 that was selected after evaluating variants showing a greater risk in more aggressive prostate cancer (Gleason ≥7 and/or T3 or higher and/or node positive and/or metastatic disease) than in less aggressive disease (Gleason <7 and T2 or lower).
Allele A of the SNP rs5945572 (rs5945572 A) showed the most significant P value of markers located on the X chromosome in a GWA analysis of 23,205 Icelandic samples, with an allelic specific OR of 1.21 (P = 3.36 × 10−4; Table 1). The SNP is located on Xp11.22, a region that has not been previously implicated in prostate cancer and that is correlated with neither the androgen receptor locus on Xq12 nor the HPCX locus on Xq27–28. On chromosome 2p15, allele A of rs2710646 (rs2710646 A) had an OR of 1.16 (P = 7.79 × 10−4; Table 2) for all Icelandic prostate cancer cases and an OR of 1.33 (P = 3.73 × 10−5) for individuals diagnosed with aggressive disease (Supplementary Table 1 online). By examining the Utah CEPH (CEU) HapMap data, we identified numerous SNPs that are substantially correlated with either of the two anchor SNPs, rs5945572 and rs2710646 (D′ ≥ 0.9 and r2 > 0.4). From this list we selected several SNPs not present on the Illumina Hap300 chip, representing different degrees of correlation with the anchor SNPs, for further genotyping in 1,500 and 800 Icelandic cases and controls, respectively. None of the additional SNPs were found to be more significantly associated with the disease than the anchor SNPs (Supplementary Table 2 online).
Table 1
Table 1
Association results for rs5945572 on Xp11.22 and prostate cancer in Iceland, The Netherlands, Spain, Sweden and the United States
Table 2
Table 2
Association results for rs721048 on 2p15 and prostate cancer in Iceland, The Netherlands, Spain, Sweden and the United States
We proceeded to genotype rs5945572 and rs2710646 in seven prostate cancer study groups of European descent, with populations from The Netherlands, Spain, Sweden and the United States. However, as the TaqMan assay for rs2710646 on 2p15 failed in design, we replaced it with an assay for a fully correlated SNP, rs721048 (linkage disequilibrium (LD) characteristics between rs2710646 and rs721048 in Icelanders and in the four HapMap populations (CEU, CHB, JPT and YRI; Supplementary Table 3 online); D′ = 1; r2 ≥ 0.99), which was used for genotyping in the replication study populations. When results from all seven case-control replication groups were combined, using all prostate cancer cases, they were highly significant for both SNPs with OR = 1.24 (P = 2.57 × 10−10) for rs5945572 A on Xp11, and OR = 1.15 (P = 2.23 ×10−6) for rs721048 A on 2p15 (Tables 1 and and2).2). When we combined these data with the Icelandic data, the signals at both loci achieved genome-wide significance: rs5945572 A on Xp11.22 had an OR = 1.23 (P = 3.95 ×10−13), and rs721048 A on 2p15 had an OR = 1.15 (P = 7.66 × 10−9) (Tables 1 and and2).2). Removing all females (n = 14,135) from the control group in the analysis of the combined groups yielded essentially identical results for both loci (OR = 1.23 and 1.15 for rs5945572 A and rs721048 A, respectively; Supplementary Methods). In all of the replication groups, the observed effect for the two loci was in the same direction as in the Icelandic samples, except for rs721048 A, which in the Baltimore group showed no effect (OR = 1). However, a test of heterogeneity in the OR of the eight different study groups showed no significant difference for the two loci (P = 0.89 and 0.19 for Xp11 and 2p15, respectively). We note that in the results released by the Cancer Genetics Markers of Susceptibility study group5 (http://cgems.cancer.gov/data/), the two original anchor SNPs, rs5945572 and rs2710646, show nominal, but not genome-wide, significant association to prostate cancer, further supporting the data presented here.
For rs5945572 A on Xp11, the OR seen for cases with younger age at onset (≤65) or aggressive phenotype was the same as for the whole group. The frequency of rs721048 A, on the other hand, was significantly higher among individuals diagnosed with aggressive prostate cancer than among those with less aggressive disease (OR = 1.11; P = 2.6 × 10−3). Comparing the group of individuals with aggressive tumor (n = 4,787) to controls gave an allelic OR of 1.22 (P = 8.7 ×10−10) when combined for all study groups (Supplementary Table 1). Specifically, the heterozygous and homozygous carriers of rs721048 A, which are close to 31% and 4% of the general population, have a genotypic OR of 1.22 and 1.49 of being diagnosed with aggressive prostate cancer compared to the noncarriers, respectively. However, rs721048 A, like rs5945572 A, did not show a stronger association in individuals with younger age at onset of prostate cancer.
Both of the variants, rs5945572 on Xp11 and rs2710646 and rs721048 on 2p15, are located within regions characterized by extended linkage disequilibrium, according to the CEU HapMap results. On Xp11.22, the LD block spans the genomic region from 50.832 Mb to 51.553 Mb, covering 721 kb (NCBI Build 35). The 2p15 LD region is even larger—about 1.9 Mb (62.621–64.534 Mb; NCBI Build 35)—and contains several genes. Both of the 2p15 SNPs, rs2710646 and rs721048, are located within one of the introns of the EHBP1 gene. This gene is thought to be involved in endocytic trafficking and has not previously been implicated in cancer6. Through an RT-PCR analysis of various cDNA libraries, we detected expression of EHBP1 in several different tissue libraries, including those derived from the prostate (data not shown).
On Xp11, several genes have been localized to the 721-kb region of interest, but none have been previously linked with prostate cancer. Possible cancer candidate genes, based on previously described functions, are GSPT2 and MAGED1. GSPT2 is related to GSPT1, a gene encoding a GTP-binding protein that has an essential role at the G1-to S-phase transition of the cell cycle. MAGED1 has been implicated in programmed cell death through a JNK- and JUN-dependent mitochondrial pathway7,8. The genes closest to the SNP rs5945572 are NUDT10 and NUDT11, along with a single exon transcript, LOC340602, with unknown function (University of California Santa Cruz Genome Browser, May 2004). NUDT10 and NUDT11 belong to a subgroup of phosphohydrolases that preferentially hydrolyze diphosphoinositol polyphosphates (DIPPS)9. It has been proposed that members of this protein superfamily may be involved in vesicle trafficking, stress responses, DNA repair and apoptosis10,11. Using RT-PCR analysis of various cDNA libraries, we detected expression of GSPT2, MAGED1, LOC340602, NUDT10 and NUDT11 in several different tissue libraries, including those derived from the prostate (data not shown). Which one of these genes, if any, confers the risk identified by the association reported here remains to be shown.
Recently, several sequence variants have been identified that together account for a substantial proportion of the population attributable risk (PAR) of prostate cancer1,2,5,12,13. With the identification of the two variants described here, yet another piece has been added to the puzzle of prostate cancer genetic susceptibility. Both the Xp11 and the 2p15 variants are common and confer a moderate risk, resulting in an estimated PAR of about 7% and 5% in individuals of European descent, respectively. However, the underlying causative biological perturbation linked to these variants remains to be elucidated.
Supplementary Material
Suppl Data
Acknowledgments
We thank the individuals who participated in the study and whose contribution made this work possible. This project was funded in part by contract number 018827 (Polygene) from the 6th Framework Program of the EU to deCODE genetics, T.R. and L.A.K. The study was supported in part by National Cancer Institute CA105055, CA106523 and CA95052 to J.X., CA112517 and CA58236 to W.B.I., CA86323 to A.W.P., and Department of Defense grant PC051264 to J.X., and in part by a V Foundation award and US Department of Veterans Affairs grants to J.R.S. The principal investigators and corresponding authors for the respective replication study populations are L.A.K. (The Netherlands), J.I.M. (Spain), H.G. (Sweden), W.B.I. (Baltimore), W.J.C. (Chicago), J.R.S. (Nashville) and S.N.T. (Rochester).
Footnotes
Note: Supplementary information is available on the Nature Genetics website.
COMPETING INTERESTS STATEMENT
The authors declare competing financial interests: details accompany the full-text HTML version of the paper at http://www.nature.com/naturegenetics/.
References
1. Gudmundsson J, et al. Nat. Genet. 2007;39:631–637. [PubMed]
2. Gudmundsson J, et al. Nat. Genet. 2007;39:977–983. [PubMed]
3. Narod SA, et al. Nat. Med. 1995;1:99–101. [PubMed]
4. Monroe KR, et al. Nat. Med. 1995;1:827–829. [PubMed]
5. Yeager M, et al. Nat. Genet. 2007;39:645–649. [PubMed]
6. Guilherme A, et al. J. Biol. Chem. 2004;279:10593–10605. [PubMed]
7. Salehi AH, Xanthoudakis S, Barker PA. J. Biol. Chem. 2002;277:48043–48050. [PubMed]
8. Salehi AH, et al. Neuron. 2000;27:279–288. [PubMed]
9. Hidaka K, et al. J. Biol. Chem. 2002;277:32730–32738. [PubMed]
10. Dubois E, et al. J. Biol. Chem. 2002;277:23755–23763. [PubMed]
11. Morrison BH, Bauer JA, Kalvakolanu DV, Lindner DJ. J. Biol. Chem. 2001;276:24965–24970. [PMC free article] [PubMed]
12. Amundadottir LT, et al. Nat. Genet. 2006;38:652–658. [PubMed]
13. Haiman CA, et al. Nat. Genet. 2007;39:638–644. [PMC free article] [PubMed]