scieee Open visual document viewer

Germline variation at 8q24 and prostate cancer risk in men of European ancestry

Matejcic, Marco,Saunders, Edward J.,Tammela, Teuvo

Full text

ARTICLE Ge mline a ia ion a 8q24 and p os a e cance isk in men o Eu opean ances y Ma co Ma ejcic1, Edwa d J. Saunde s e al. # Ch omosome 8q24 is a suscep ibili y locus o mul iple cance s, including p os a e cance . He e we combine gene ic da a ac oss he 8q24 suscep ibili y egion om 71,535 p os a e cance cases and 52,935 con ols o Eu opean ances y o define he o e all con ibu ion o ge mline a ia ion a 8q24 o p os a e cance isk. We iden i y 12 independen isk signals o p os a e cance (p< 4.28 × 10−15), including h ee isk a ian s ha ha e ye o be epo ed. F om a polygenic isk sco e (PRS) model, de i ed o assess he cumula i e e ec o isk a ian s a 8q24, men in he op 1% o he PRS ha e a 4- old (95%CI =3.62–4.40) g ea e isk compa ed o he popula ion a e age. These 12 a ian s accoun o ~25% o wha can be cu en ly explained o he amilial isk o p os a e cance by known gene ic isk ac o s. These findings highligh he o e whelming con ibu ion o ge mline a ia ion a 8q24 on p os a e cance isk which has implica ions o popula ion isk s a ifica ion. Ma co Ma ejcic, Edwa d J. Saunde s e al. # DOI: 10.1038/s41467-018-06863-1 OPEN Co espondence and eques s o ma e ials should be add essed o C.A.H. (email: Ch is op[email p o ec ed]). # A ull lis o au ho s and hei a filia ions appea s a he end o he pape . NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions 1 1234567890():,; P os a e cance (PCa) is he mos common cance among men in he US, wi h 161,360 new cases and 26,730 ela ed dea hs es ima ed in 20171. Familial and epidemiological s udies ha e p o ided e idence o subs an ial he i abili y o PCa2, and ~170 common isk loci ha e been iden ified h ough genome-wide associa ion s udies (GWAS)3. The suscep ibili y egion on ch omosome 8q24 has been shown o be a majo con ibu o o PCa isk, wi h mul iple a ian s clus e ed in fi e linkage disequilib ium (LD) blocks spanning ~600 Mb ha a e independen ly associa ed wi h isk4. Many o hese associa ion signals epo ed a 8q24 ha e been eplica ed ac oss acial/e hnic popula ions5,6, poin ing o common sha ed unc ional a ian s wi hin 8q24. Howe e , a e ances y-specific a ian s ha e also been de ec ed, which con e la ge ela i e isks o PCa (odds a ios [ORs] >2.0) han common isk a ian s in he egion and signi y allelic he e ogenei y in he con ibu ion o ge mline a - ia ion a 8q24 o PCa isk ac oss popula ions7. In he cu en s udy, we pe o m a comp ehensi e in es iga- ion o gene ic a ia ion ac oss he 1.4 Mb cance suscep ibili y egion a 8q24 (127.6–129.0 Mb) in ela ion o PCa isk. We combine geno yped and impu ed da a om wo la ge GWAS conso ia (PRACTICAL/ELLIPSE OncoA ay and iCOGS) including >124,000 indi iduals o Eu opean ances y o sea ch o no el isk a ian s, as well as o de e mine he o e all con ibu- ion o gene ic a ia ion a 8q24 o PCa he i abili y. Ou findings unde sco e he sizable impac o gene ic a ia ion in he 8q24 egion in explaining in e -indi idual di e ences in PCa isk, wi h po en ial clinical u ili y o gene ic isk p edic ion. Resul s Ma ginal and condi ional associa ion analysis. Geno ype da a om he Illumina OncoA ay and iCOGS a ay and impu a ion o 1000 Genomes P ojec (1KGP) we e gene a ed among 71,535 PCa cases and 52,935 con ols o Eu opean ances y om 86 case- con ol s udies (see Me hods). O he 5600 geno yped and impu ed a ian s a 8q24 (127.6–129.0 Mb) wi h mino allele equency (MAF) > 0.1% e ained o analysis (see Me hods), 1268 (23%) we e associa ed wi h PCa isk a p< 5×10−8while 2772 (49%) we e ma ginally associa ed a p< 0.05. These 5600 ma ke s cap u e, a 2> 0.8, 90% and 97% o all a ian s a 8q24 (127.6–129.0 Mb) wi h MAF ≥1% and ≥5%, espec i ely (based on 1KGP Phase 3 EUR panel). In a o wa d and backwa d s epwise selec ion model on a ian s ma ginally associa ed wi h PCa isk (p< 0.05, n=2772; see Me hods), we iden ified 12 a ian s wi h condi ional p- alues om he Wald es be ween 2.93 × 10−137 and 4.28 × 10−15 (Table 1). None o he o he a ian s we e s a is ically significan a p<5×10 −8a e adjus - men o he 12 independen hi s (Fig. 1). The 8q24 egion is shown in Supplemen a y Fig. 1. O hese 12 s epwise signals, h ee had alleles wi h ex eme isk allele equencies (RAFs) ha con- eyed la ge e ec s ( s77541621, RAF =2%, OR =1.85, 95%CI = 1.76–1.94; s183373024, RAF =1%, OR =2.67, 95%CI = 2.43–2.93; s190257175, RAF =99%, OR =1.60, 95%CI = 1.42–1.80). The emaining a ian s had RAFs be ween 0.11 and 0.92 and condi ional ORs ha we e mo e modes and anged om 1.10 o 1.37 (Table 1). Fo 8 o he 12 a ian s, he allele ound o be posi i ely associa ed wi h PCa isk was he p e- dominan allele (i.e., >50% in equency). Fo wo a ian s, s78511380 and s190257175, he ma ginal associa ions we e no genome-wide significan and subs an ially weake han hose in he condi ional model. Fo s78511380, he ma ginal OR was sligh ly p o ec i e (OR =0.97; p=0.027), bu e e sed di ec ion and was highly s a is ically significan when condi ioning on he o he 11 a ian s (OR =1.19; p=3.5 × 10−18; Table 1). Haplo ype analysis. The haplo ype analysis showed an addi i e e ec o he 12 independen isk a ian s consis en wi h ha p edic ed in he single a ian es ; co-occu ence o he 8q24 isk alleles on he same haplo ype does no u he inc ease he isk o PCa (Supplemen a y Table 1). The unique haplo ype ca ying he e e ence allele o s190257175 (GCTTAT, 0.5% equency) is also he sole haplo ype associa ed wi h a educed isk o PCa, sugges ing ha ha ing he C allele con e s a p o ec i e e ec . The e e ence allele o s78511380 (A, 8% equency) occu s on a haplo ype in block 2 oge he wi h he isk alleles o s190257175, s72725879 and s5013678 (haplo ype GTTTAA, 8%) which obscu es he posi i e associa ion wi h he T allele o s78511380. Thus, he ma ginal p o ec i e e ec associa ed wi h he isk allele o s78511380 eflec s an inc eased isk associa ed wi h he occu ence on a isk haplo ype wi h o he isk alleles (Supplemen a y Table 1). Table 1 Ma ginal and condi ional es ima es o gene ic ma ke s a 8q24 independen ly associa ed wi h p os a e cance isk Va ian IDaPosi ionbAllelecRAFdLD clus e eCondi ional associa ion Ma ginal associa ion OR (95%CI)gp- alue OR (95%CI)hp- alue s1914295 127910317 T/C 0.68 block 1 1.10 (1.08–1.12) 7.30 × 10−25 1.09 (1.07–1.11) 3.07 × 10−21 s1487240 128021752 A/G 0.74 block 1 1.20 (1.17–1.22) 2.77 × 10−66 1.16 (1.14–1.18) 2.97 × 10−54 s77541621 128077146 A/G 0.02 block 2 1.85 (1.76–1.94) 2.93 × 10−137 1.83 (1.74–1.92) 4.33 × 10−137 s190257175 128103466 T/C 0.99 block 2 1.60 (1.42–1.80) 4.28 × 10−15 1.36 (1.22–1.53) 6.90 × 10−8 s72725879 128103969 T/C 0.18 block 2 1.31 (1.28–1.35) 1.26 × 10−83 1.17 (1.14–1.19) 3.96 × 10−48 s5013678 128103979 T/C 0.78 block 2 1.10 (1.08–1.13) 1.58 × 10−19 1.20 (1.17–1.22) 4.44 × 10−68 s183373024 128104117 G/A 0.01 block 2 2.67 (2.43–2.93) 4.89 × 10−95 3.20 (2.92–3.50) 6.60 × 10−138 s78511380 128114146 T/A 0.92 block 2 1.19 (1.14–1.23) 3.48 × 10−18 0.97 (0.94–1.00) 0.027 s17464492 128342866 A/G 0.72 block 3 1.16 (1.14–1.18) 3.01 × 10−52 1.17 (1.15–1.19) 9.05 × 10−61 s6983267 128413305 G/T 0.51 block 4 1.18 (1.16–1.20) 5.68 × 10−84 1.23 (1.21–1.25) 3.15 × 10−135 s7812894 128520479 A/T 0.11 block 5 1.37 (1.33–1.40) 1.55 × 10−122 1.45 (1.41–1.49) 1.20 × 10−181 s12549761 128540776 C/G 0.87 block 5 1.21 (1.18–1.24) 1.61 × 10−45 1.28 (1.25–1.31) 1.38 × 10−78 aVa ian s ha emained genome-wide significan ly associa ed wi h PCa isk (p<10 −8) in he final s epwise model bCh omosome posi ion based on human genome build 37 cRisk allele/ e e ence allele dRisk allele equency eLD clus e s we e in e ed based on ecombina ion ho spo s using Haplo iew 4.229 and defined as p e iously epo ed by Al Olama e al.4 Each a ian was inco po a ed in he s epwise model based on he s eng h o ma ginal associa ion om he me a-analysis o OncoA ay and iCOGS da a gPe -allele odds a io and 95% confidence in e al adjus ed o coun y, 7(OncoA ay)/8(iCOGS) p incipal componen s and all o he a ian s in he able hPe -allele odds a io and 95% confidence in e al adjus ed o coun y and 7(OncoA ay)/8(iCOGS) p incipal componen s ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 2NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions Co ela ion wi h known isk loci. The 12 isk a ian s spanned ac oss he fi e LD blocks p e iously epo ed o ha bo isk a ian s o PCa a 8q244, wi h block 2 ha bo ing six signals, blocks 1 and 5 wo signals each, and blocks 3 and 4 only one (Supplemen a y Fig. 2). Excep o a weak co ela ion be ween s72725879 and s78511380 in block 2 ( 2=0.28), he isk a - ian s we e unco ela ed wi h each o he ( 2≤0.09; Supplemen a y Da a 1), which co obo a es hei independen associa ion wi h PCa isk. Eigh o he a ian s ( s1487240, s77541621, s72725879, s5013678, s183373024, s17464492, s6983267, s1914295 s1487240 s77541621 s190257175 s72725879 s5013678 s183373024 s78511380 s17464492 s6983267 s7812894 s12549761 0 20 40 60 80 100 120 140 160 180 200 −Log10(p) −Log10(p) Ma ginal s1914295 s1487240 s77541621 s190257175 s72725879 s5013678 s183373024 s78511380 s17464492 s6983267 s7812894 s12549761 0 20 40 60 80 100 120 140 Condi ional s1914295 s1487240 s77541621 s190257175 s72725879 s5013678 s183373024 s78511380 s17464492 s6983267 s7812894 s12549761 LD 2 Hi s His one 128,000,000 128,200,000 128,400,000 DNaseI Conse ed Ch omHMM eQTL Bio ea u es PCAT1 POU5F1B CASC8 128,000,000 128,200,000 128,400,000 Genes FAM84B POU5F1B NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 ARTICLE NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions 3 s7812894) ha e been p e iously epo ed ei he di ec ly (Sup- plemen a y Table 2) o a e co ela ed ( 2≥0.42) wi h known ma ke s o PCa isk om s udies in popula ions o Eu opean, A ican o Asian ances y (Supplemen a y Da a 1)4,7–10. The ma ginal es ima es o p e iously published PCa isk a ian s a 8q24 in he cu en s udy a e shown in Supplemen a y Table 2. The a ian s1914295 in block 1 is only weakly co ela ed wi h he p e iously epo ed isk a ian s s12543663 and s10086908 ( 2=0.17 and 0.14, espec i ely), while s7851380 is modes ly co ela ed wi h he p e iously epo ed isk a ian s1016343 ( 2=0.28). The emaining wo a ian s, s190257175 and s12549761, a e no co ela ed ( 2< 0.027) wi h any known PCa isk ma ke . Polygenic isk sco e and amilial ela i e isk. To es ima e he cumula i e e ec o ge mline a ia ion a 8q24 on PCa isk, a polygenic isk sco e (PRS) was calcula ed o he 12 independen isk alleles om he final model based on allele dosages weigh ed by he pe -allele condi ionally adjus ed ORs (see Me hods). Compa ed o he men a ‘a e age isk’(i.e., he 25 h–75 h PRS ange among con ols), men in he op 10% o he PRS dis ibu ion had a 1.93- old ela i e isk (95%CI =1.86–2.01) (Table 2), wi h he isk being 3.99- old highe (95%CI = 3.62–4.40) o men in he op 1%. Risk es ima es by PRS ca ego y we e no modified by amily his o y (FamHis -yes: OR =4.24, 95%CI =2.85–6.31; FamHis -no: OR =3.38, 95%CI = 2.88–3.97). To quan i y he impac o ge mline a ia ion a 8q24, we also es ima ed he p opo ion o amilial ela i e isk (FRR) and he i abili y o PCa con ibu ed by 8q24 and compa ed his o he p opo ions explained by all known PCa isk a ian s including 8q24 (see Me hods). The 175 es ablished PCa sus- cep ibili y loci iden ified o da e3,11 a e es ima ed o explain 37.08% (95%CI =32.89–42.49) o he FRR o PCa, while he 12 independen signals a 8q24 alone cap u e 9.42% (95%CI = 8.22–10.88), which is 25.4% o he o al FRR explained by known gene ic isk ac o s o PCa (Table 3). This is simila o he p opo ion o he i abili y explained by 8q24 a ian s (22.2%) compa ed o he o al explained he i abili y by he known isk a ian s (0.118). In compa ison, he nex highes con ibu ion o an indi idual suscep ibili y egion o he FRR o PCa is he TERT egion a ch omosome 5p15, whe e 5 independen signals con- ibu ed 2.63% (95%CI =2.34–3.00). No o he indi idual GWAS Fig. 1 LocusExplo e plo s o he 12 a ian s a 8q24 significan ly associa ed wi h PCa isk. ‘Ma ginal’and ‘Condi ional’Manha an plo panels show ma ginal and condi ional associa ion esul s, espec i ely. Va ian posi ions (x-axis) and −log 10 p- alues om he Wald es (y-axis) a e shown, wi h he ed line indica ing he h eshold o genome-wide significan associa ion wi h PCa isk (p≤5×10 −8) and blue peaks local es ima es o ecombina ion a es. The posi ion o he 12 independen a ian s is labeled in each plo . Clus e s o co ela ed a ian s o each independen signal a e dis inguished using di e en colo s and also depic ed on he ‘LD 2Hi s’ ack. S onge shading indica es g ea e co ela ion wi h he lead a ian , wi h a ian s no co ela ed a 2≥0.2 wi h any lead a ian uncolo ed. Pai wise co ela ions a e based on he Eu opean ances y (EUR) panel om he 1000 Genomes P ojec (1KGP) Phase 3. The ela i e posi ion o Re Seq genes and biological anno a ions a e shown in he ‘Genes’and ‘Bio ea u es’panels, espec i ely. Genes on he posi i e s and a e deno ed in g een and hose on he nega i e s and in pu ple. Anno a ions displayed a e: his one modifica ions in ENCODE ie 1 cell lines (His one ack), he posi ions o any a ian s ha we e eQTLs wi h p os a e umo exp ession in TCGA p os a e adenoca cinoma samples and he espec i e genes o which exp ession is al e ed (eQTL ack), ch oma in s a e ca ego iza ions in he P EC cell-line by Ch omHMM (Ch omHMM ack), he posi ion o conse ed elemen peaks (Conse ed ack) and he posi ion o DNaseI hype sensi i i y si e peaks in ENCODE p os a e cell-lines (DNaseI ack). The da a displayed in his plo may be explo ed in e ac i ely h ough he LocusExplo e applica ion (h p://www.oncogene ics.ic .ac.uk/8q24/) Table 2 Rela i e isk o PCa o polygenic isk sco e (PRS) g oups Risk ca ego y pe cen ileaNo. o indi iduals Risk es ima es o PRS g oups Con ols Cases OR (95% CI)bp- alue ≤1% 530 339 0.52 (0.45–0.59) 2.11 × 10−20 1%–10% 4771 3636 0.62 (0.59–0.65) 6.26 × 10−90 10%–25% 7936 7359 0.75 (0.72–0.78) 3.62 × 10−54 25%–75% 26,464 32,743 1.00 (Re ) 75%–90% 7940 13,431 1.37 (1.32–1.41) 6.55 × 10−77 90%–99% 4766 11,451 1.93 (1.86–2.01) 4.13 × 10−249 >99% 528 2576 3.99 (3.62–4.40) 5.64 × 10−172 No e: PRS we e calcula ed o a ian s om he final s epwise model wi h allele dosage om OncoA ay and iCOGS weigh ed by he pe -allele condi ionally adjus ed odds a ios om he me a-analysis aRisk ca ego y g oups we e based on he pe cen ile dis ibu ion o isk alleles in o e all con ols bEs ima ed e ec o each PRS g oup ela i e o he in e qua ile ange (25–75%) in OncoA ay and iCOGS da ase s sepa a ely, and hen me a-analyzed ac oss he wo s udies; odds a ios we e adjus ed o coun y and 7(OncoA ay)/8(iCOGS) p incipal componen s Table 3 P opo ion o amilial ela i e isk (FRR) and he i abili y (h g 2) o PCa explained by known isk a ian s Sou ce No. o a ian s P opo ion o FRR (95%CI) % o o al FRR h g 2(SE) % o o al h g 2 8q24a12 9.42 (8.22–10.88) 25.4 0.027 (0.011) 22.2 HOXB13b1 1.91 (1.20–2.85) 5.2 0.004 (0.005) 3.0 All o he a ian sb,c 162 25.77 (22.94–29.36) 69.5 0.092 (0.010) 74.9 To al 175 37.08 (32.89–42.49) 100 0.118 (0.012) 100 aCondi ional es ima es we e de i ed by fi ing a single model wi h all a ian s om OncoA ay da a bRisk es ima es and allele equencies o egions wi h a single a ian a e om a me a-analysis o OncoA ay, iCOGS and 6 addi ional GWAS3 cRisk a ian s included om fine-mapping o PCa suscep ibili y loci in Eu opean ances y popula ions11 ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 4NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions locus has been es ablished as explaining >2% o he FRR, including he low equency, non-synonymous, mode a e pene- ance HOXB13 a ian ( s138213197) a ch omosome 17q21 ha is es ima ed o explain only 1.91% (95%CI =1.20–2.85) o he FRR11. JAM analysis. We explo ed ou da a wi h a second fine-mapping app oach, JAM (Join Analysis o Ma ginal summa y s a is ics)12, which uses GWAS summa y s a is ics o iden i y c edible se s o a ian s ha define he independen associa ion signals in sus- cep ibili y egions (see Me hods). The 95% c edible se o he JAM analysis confi med all o he independen signals om s epwise analysis excep s190257175, o which e idence o an associa ion was weak ( a ian -specific Bayes ac o (BF) =1.17). The e we e 50 o al a ian s included in he 95% c edible se , and 174 a e including a ian s in high LD ( 2> 0.9) wi h hose in he c edible se (Supplemen a y Da a 2). Discussion In his la ge s udy o ge mline gene ic a ia ion ac oss he 8q24 egion, we iden ified 12 independen associa ion signals among men o Eu opean ances y, wi h h ee o he isk a ian s ( s1914295, s190257175, and s12549761) being weakly co e- la ed ( 2≤0.17) wi h known PCa isk ma ke s. The combina ion o hese 12 independen signals a 8q24 cap u e app oxima ely one qua e o he o al PCa FRR explained by known gene ic isk ac o s, which is subs an ially g ea e han any o he known PCa isk locus. The 8q24 egion is he majo suscep ibili y egion o PCa; howe e , he unde lying biological mechanism(s) h ough which ge mline a ia ion in his egion influences PCa isk emains unce ain. Fo each o he 12 isk a ian s a 8q24, he 95% c edible se defined no ewo hy (i.e., pu a i e unc ional) a ian s based on summa y s a is ics while accoun ing o LD. To in o m biological unc ionali y, we o e laid epigene ic unc ional anno- a ion using publicly a ailable da ase s (see Me hods) wi h he loca ion o he 12 independen signals (and co esponding 174 a ian s wi hin hei 95% c edible se s; Supplemen a y Da a 3). O he 12 independen lead a ian s, 6 a e si ua ed wi hin pu a i e ansc ip ional enhance s in p os a e cell-lines; ei he h ough in e sec ion wi h H3K27AC ( s72725879, s5013678, s78511380, s6983267 and s7812894) o h ough a Ch omHMM enhance anno a ion ( s17464492, s6983267, s7812894). Eigh o he 12 s epwise hi s ( s77541621, s190257175, s5013678, s183373024, s78511380, s17464492, s6983267, s7812894) also in e sec ansc ip ion ac o binding si e peaks om mul i- ple ChIP-seq da ase s ep esen ing he AR, ERG, FOXA1, GABPA, GATA2, HOXB13, and NKX3.1 ansc ip ion ac o s, wi h all 8 in e sec ing a FOXA1 ma k and hal an AR binding si e. These a ian s may he e o e exe hei e ec h ough egula ion o enhance ac i i y and long- ange exp ession o genes impo - an o cance umo igenesis and/o p og ession13. The a ian s6983267 has also been shown o ac in an allele-specific manne o egula e p os a e enhance ac i i y and exp ession o he p o o- oncogene MYC in i o and in i o14,15. Howe e , despi e he close p oximi y o he MYC locus, no di ec associa ion has been de ec ed be ween 8q24 isk alleles and MYC exp ession in no mal and umo human p os a e issues16. The a e a ian wi h he la ges e ec on isk, s183373024, shows high e idence o unc ionali y based on o e lap wi h mul iple DNaseI and ansc ip ion ac o binding si e peaks ( o AR, FOXA1, HOXB13, and NKX3.1), which suppo s p e ious findings o an allele-dependen e ec o his a ian on he dis up ion o a FOXA1 binding mo i 17. Se en independen signals ( s1914295, s1487240, s77541621, s72725879, s5013678, s183373024, s78511380) and a ian s co ela ed a 2> 0.9 wi h hese signals (Supplemen a y Da a 2) a e loca ed wi hin o nea a numbe o p os a e cance –associa ed long noncoding RNAs (lncRNAs), including PRNCR1,PCAT1, and CCAT2, p e iously epo ed o be up egula ed in human PCa cells18 and issues19,20. Based on eQTL anno a ions in p os a e adenoca cinoma cells, he inde- penden signal s1914295 and h ee co ela ed a ian s ( 2> 0.9; Supplemen a y Da a 2) a e associa ed wi h o e exp ession o FAM84B, a gene p e iously associa ed wi h p og ession and poo p ognosis o PCa in animal s udies21. Va ian s co ela ed a 2> 0.9 wi h s7812894 (n=9; Supplemen al Table 4) a e eQTLs o POU5F1B, a gene o e exp essed in cance cell lines and cance issues22,23, al hough i s ole in PCa de elopmen is unknown. Whils we ha e success ully efined he 8q24 egion and iden ified a subse o a ian s wi h pu a i e biological unc ion wi hin ou c edible se , mul i-e hnic compa isons may help efine he asso- cia ion signals e en u he and p ecisely iden i y he unc ional alleles and biological mechanisms ha modi y PCa isk. Whe eas he indi idual associa ions o he 8q24 a ian s wi h PCa isk a e ela i ely modes (ORs < 2.0, excep o s183373024), hei cumula i e e ec s a e subs an ial, wi h isk being 4- old highe o men in he op 1% o he 8q24-only PRS. The con ibu ion o he o e all FRR o PCa is subs an ially g ea e o he 8q24 egion (9.42%) han o any o he known GWAS locus, including he mode a e pene ance non-synonymous a - ian in HOXB13 (1.91%). The abili y o hese ma ke s o explain ~25.4% o wha can be cu en ly explained by all known PCa isk a ian s is a clea indica ion o he impo an con ibu ion o ge mline a ia ion a 8q24 on PCa isk. Ou s udy was p e- dominan ly powe ed o analyze a ian s wi h MAF > 1% as he impu ed a ian s wi h MAF =0.1-1% we e mos likely o ail quali y con ol (QC); howe e , he high densi y o geno yped ma ke s and haplo ypes a 8q24 in he OncoA ay and iCOGS s udies p o ided a obus backbone o impu a ion and inc eased he chances o impu e lowe MAF a ian s wi h high impu a ion quali y sco e. Unde s anding o he biology o hese a ian s and he unde lying gene ic basis o PCa could p o ide new insigh s in o he iden ifica ion o eliable isk-p edic ion bioma ke s o PCa, as well as enable he de elopmen o e ec i e s a egies o a ge ed sc eening and p e en ion. Me hods S udy subjec s, geno yping, and quali y con ol. We combined geno ype da a om he PRACTICAL/ELLIPSE OncoA ay and iCOGS conso ia3,24, which included 143,699 men o Eu opean ances y om 86 case-con ol s udies la gely based in ei he he US o Eu ope. In each s udy, cases p ima ily included men wi h inciden PCa while con ols we e men wi hou a p io diagnosis o he disease. Bo h o he OncoA ay and iCOGS cus om a ays we e designed o p o ide high co e age o common alleles (mino allele equency [MAF] > 5%) ac oss 8q24 (127.6–129.0 Mb) based on he 1000 Genomes P ojec (1KGP) Phase 3 o OncoA ay, and he Eu opean ances y (EUR) panel om HapMap Phase 2 o iCOGS. A o al o 57,580 PCa cases and 37,927 con ols o Eu opean ances y we e geno yped wi h he Illumina OncoA ay, and 24,198 PCa cases and 23,994 con ols o Eu opean ances y we e geno yped wi h he Illumina iCOGS a ay. Fo bo h s udies, sample exclusion c i e ia included duplica e samples, fi s -deg ee ela i es, samples wi h a call a e <95% o wi h ex eme he e ozygosi y (p<10 −6), and samples wi h an es ima ed p opo ion o Eu opean ances y <0.83,24. In o al, geno ype da a o 53,449 PCa cases and 36,224 con ols om OncoA ay and 18,086 PCa cases and 16,711 con ols om iCOGS we e included in he analysis. Gene ic a ian s wi h call a es <0.95, de ia ion om Ha dy-Weinbe g equilib ium (p<10 −7in con ols), and geno ype disc epancy in >2% o duplica e samples we e excluded. O he final 498,417 geno yped a ian s on he OncoA ay and 201,598 on he iCOGS a ay ha passed QC, 1581 and 1737 wi hin he 8q24 egion, espec i ely, we e e ained o impu a ion. All s udies complied wi h all ele an e hical egula ions and we e app o ed by he ins i u ional e iew boa ds a each o he pa icipa ing ins i u ions. In o med consen was ob ained om all s udy pa icipan s. Addi ional de ails o each s udy a e p o ided in he Supplemen a y No e 1. Impu a ion analysis. Impu a ion o bo h OncoA ay and iCOGS geno ype da a was pe o med using SHAPEIT25 and IMPUTE 226 o he Oc obe 2014 (Phase 3) NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 ARTICLE NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions 5 elease o he 1KGP e e ence panel. A o al o 10,136 a ian s om OncoA ay and 10,360 a ian s om iCOGS wi h MAF > 0.1% we e impu ed ac oss he isk egion a 8q24 (127.6-129.0 Mb). Va ian s wi h an impu a ion quali y sco e >0.8 we e e ained o a o al o 5600 o e lapping a ian s be ween he wo da ase s. S a is ical analysis. Uncondi ional logis ic eg ession was used o es ima e pe -allele odds a ios (ORs) and 95% confidence in e als (CIs) o he associa ion be ween gene ic a ian s (single nucleo ide polymo phisms and inse ion/ dele ion polymo phisms) and PCa isk adjus ing o coun y and p incipal componen s (7 o OncoA ay and 8 o iCOGS). Allele dosage e ec s we e es ed h ough a 1-deg ee o eedom wo- ailed Wald end es . The ma ginal isk es ima es o he 5600 a ian s a 8q24 ha passed QC we e combined by a fixed e ec me a-analysis wi h in e se a iance weigh ing using METAL27.A modified o wa d and backwa d s epwise model selec ion wi h inclusion and exclusion c i e ia o p≤5×10 −8was pe o med on a ian s ma ginally associa ed wi h PCa isk om he me a esul s (p< 0.05, n=2772). A each s ep, he e ec es ima es o he candida e a ian s om bo h s udies (OncoA ay and iCOGS) we e me a-analyzed and each a ian was inco po a ed in o he model based on he s eng h o associa ion. All emaining a ian s we e included one-a -a- ime in o he logis ic eg ession model condi ioning on hose al eady inco po a ed in he model. We applied a conse a i e h eshold o independen associa ions, wi h a ian s kep in he model i hei me a p- alue om he Wald es was genome- wide significan a p≤5×10 −8a e adjus men o he o he a ian s in he model. Co ela ions be ween a ian s in he final model and p e iously published PCa isk a ian s a 8q24 we e es ima ed using he 1KGP Phase 3 EUR panel (Supplemen a y Da a 1). Haplo ype analysis. Haplo ypes we e es ima ed in he Oncoa ay da a only using a ian s om he final s epwise model selec ion (n=12) and he EM algo i hm28 wi hin LD block egions in e ed based on ecombina ion ho spo s using Haplo iew 4.2 (B oad Ins i u e, Camb idge, MA, USA)29. Only haplo ypes wi h an es ima ed equency ≥0.5% we e es ed. Polygenic isk sco e and amilial ela i e isk. An 8q24-only polygenic isk sco e (PRS) was calcula ed o a ian s om he final model (n=12) wi h allele dosage om OncoA ay and iCOGS weigh ed by he pe -allele condi ionally adjus ed ORs om he me a-analysis. Ca ego iza ion o he PRS was based on he pe cen ile dis ibu ion in con ols, and he isk o each ca ego y was es ima ed ela i e o he in e qua ile ange (25–75%) in OncoA ay and iCOGS sepa a ely, and hen me a- analyzed ac oss he wo s udies. We es ima ed he con ibu ion o 8q24 a ian s o he amilial (fi s -deg ee) ela i e isk (FRR) o PCa (FRR =2.5)30 unde a mul- iplica i e model, and compa ed his o he FRR explained by all known PCa isk a ian s including 8q24 (Supplemen a y Da a 4). We also es ima ed he i abili y o PCa using he LMM app oach as implemen ed in GCTA31. Fo egions which ha e been fine-mapped using he OncoA ay me a-analysis da a, we used he upda ed ep esen a i e lead a ian s, o he wise he o iginally epo ed a ian was included p o ided ha i had eplica ed a genome-wide significance in he me a-analysis; his iden ified a o al o 175 independen ly associa ed PCa a ian s o he FRR and he i abili y calcula ions3,11. Fo hese analyses, we used condi ional es ima es om fi ing a single model wi h all a ian s in he OncoA ay da ase o egions wi h mul iple a ian s and he o e all ma ginal me a-analysis esul s om Schumache e al.3 o egions wi h a single a ian . To co ec o po en ial bias in e ec es ima ion o newly disco e ed a ian s, we implemen ed a Bayesian e sion o he weigh ed co ec ion32, which inco po a es he unce ain y in he e ec es ima e in o he final es ima es o he bias-co ec ed ORs, 95%CIs and he co esponding calcula ions o pe cen FRR explained. JAM analysis. To confi m he s epwise esul s and iden i y candida e a ian s o po en ial unc ional ollow-up, we used a second fine-mapping app oach, JAM (Join Analysis o Ma ginal summa y s a is ics)12. JAM is a mul i a ia e Bayesian a iable selec ion amewo k ha uses GWAS summa y s a is ics o iden i y he mos likely numbe o independen associa ions wi hin a locus and define c edible se s o a ian s d i ing hose associa ions. JAM was applied o summa y s a is ics om he me a-analysis esul s using LD es ima ed om impu ed indi idual le el da a om 20,000 cases and 20,000 con ols andomly selec ed om he OncoA ay sub-s udy. LD p uning was pe o med using P io i y P une (h p://p io i yp une . sou ce o ge.ne /) on he 2772 ma ginally associa ed a ian s a 2=0.9, esul ing in 825 ag a ian s analyzed in ou independen JAM uns wi h a ying s a ing seeds. C edible se s we e de e mined as he ag a ian s ha we e selec ed in he op models ha summed o a specific cumula i e pos e io p obabili y in all ou o he independen JAM uns, plus hei designa ed high LD p oxy a ian s om he p uning s ep. Func ional anno a ion. Va ian s in he 95% c edible se (n=50) plus a ian s co ela ed a 2> 0.9 wi h hose in he c edible se (n=174) we e anno a ed o pu a i e e idence o biological unc ionali y using publicly a ailable da ase s as desc ibed by Dadae e al.11. B iefly, a ian s we e anno a ed o p oximi y o gene (GENCODE 19), miRNA ansc ip s (miRBase elease 20), e olu iona y cons ain (acco ding o GERP++, SiPhy and Phas Cons algo i hms), likelihood o pa hogenici y (CADD 1.3) and o e lap wi h p ospec i e egula o y elemen s in p os a e-specific da ase s (DNaseI hype sensi i i y si es, H3K27Ac, H3K27me3 and H3K4me3 his one modifica ions, and o AR, CTCF, ERG, FOXA1, GABPA, GATA2, HOXB13, and NKX3.1 ansc ip ion ac o binding si es) in a mix u e o LNCaP, PC-3, P EC, RWPE1, and VCaP cell lines and human p os a e umo issues downloaded om he Cis ome Da a B owse (h p://cis ome.o g/db/). The ch oma in s a e in which each a ian esides was assessed using Ch omHMM anno a ions om wo p os a e cell lines (P EC and PC3). Cis-gene egula ion was e alua ed using 359 p os a e adenoma cases om The Cance Genome A las (TCGA PRAD; h ps://gdc-po al.nci.nih.go ) ha passed QC11. The eQTL analysis was pe o med using Fas QTL wi h 1000 pe mu a ions o each gene wi hin a 1Mb window. We hen used he me hod by Nica e al.33 ha in eg a es eQTLs and GWAS esul s in o de o e eal he subse o associa ion signals ha a e due o cis eQTLs. Fo each significan eQTL, we added he candida e a ian o he linea eg ession model o assess i he inclusion be e explains he change in exp ession o he gene. We e ie ed he p- alue o he model, assigning p- alue o 1 i he eQTL and a ian a e he same. Then we anked he p- alues in descending o de o each eQTL, and finally calcula ed he colocaliza ion sco e o each pai o eQTL and a ian s. In gene al, i an eQTL and candida e a ian ep esen he same signal, his will be eflec ed by he a ian ha ing a high p- alue, a low ank and consequen ly a high colocaliza ion sco e. Da a a ailabili y The au ho s decla e ha da a suppo ing he findings o his s udy a e a ailable wi hin he pape [and in he supplemen a y in o ma ion files]. Howe e , some o he da a used o gene a e he esul s o his s udy a e a ailable om he fi s au ho and he PRAC- TICAL Conso ium upon eques . Recei ed: 7 Feb ua y 2018 Accep ed: 1 Oc obe 2018 Re e ences 1. Siegel, R. L., Mille , K. D. & Jemal, A. Cance s a is ics, 2017. CA Cance J. Clin. 67,7–30 (2017). 2. Hjelmbo g, J. B. e al. The he i abili y o p os a e cance in he No dic Twin S udy o Cance . Cance Epidemiol. P e . Bioma k. 23, 2303–2310 (2014). 3. Schumache , F. R. e al. Associa ion analyses o mo e han 140,000 men iden i y 63 new p os a e cance suscep ibili y loci. Na . Gene .50, 928–936 (2018). 4. Al Olama, A. A. e al. Mul iple loci on 8q24 associa ed wi h p os a e cance suscep ibili y. Na . Gene . 41, 1058–1060 (2009). 5. Haiman, C. A. e al. Mul iple egions wi hin 8q24 independen ly a ec isk o p os a e cance . Na . Gene . 39, 638–644 (2007). 6. Han, Y. e al. Gene alizabili y o es ablished p os a e cance isk a ian s in men o A ican ances y. In . J. Cance 136, 1210–1217 (2015). 7. Gudmundsson, J. e al. A s udy based on whole-genome sequencing yields a a e a ian a 8q24 associa ed wi h p os a e cance . Na . Gene . 44, 1326–1329 (2012). 8. Han, Y. e al. P os a e cance suscep ibili y in men o A ican ances y a 8q24. J. Na l Cance Ins .108, dj 431 (2016). 9. Ho mann, T. J. e al. A la ge mul ie hnic genome-wide associa ion s udy o p os a e cance iden ifies no el isk a ian s and subs an ial e hnic di e ences. Cance Disco . 5, 878–891 (2015). 10. Con i, D. V. e al. Two no el suscep ibili y loci o p os a e cance in men o A ican ances y. J. Na l Cance . Ins .109, djx084 (2017). 11. Dadae , T. e al. Fine-mapping o p os a e cance suscep ibili y loci in a la ge me a-analysis iden ifies candida e causal a ian s. Na . Commun. 9, 2256 (2018). 12. Newcombe, P. J., Con i, D. V. & Richa dson, S. JAM: a scalable Bayesian amewo k o join analysis o ma ginal SNP e ec s. Gene . Epidemiol. 40, 188–201 (2016). 13. Jia, L. e al. Func ional enhance s a he gene-poo 8q24 cance -linked locus. PLoS Gene . 5, e1000597 (2009). 14. Pome an z, M. M. e al. The 8q24 cance isk a ian s6983267 shows long- ange in e ac ion wi h MYC in colo ec al cance . Na . Gene . 41, 882–884 (2009). 15. Wasse man, N. F., Aneas, I. & Nob ega, M. A. An 8q24 gene dese a ian associa ed wi h p os a e cance isk con e s di e en ial in i o ac i i y o a MYC enhance . Genome Res. 20, 1191–1197 (2010). 16. Pome an z, M. M. e al. E alua ion o he 8q24 p os a e cance isk locus and MYC exp ession. Cance Res. 69, 5568–5574 (2009). 17. Hazele , D. J., Coe zee, S. G. & Coe zee, G. A. A a e a ian , which des oys a FoxA1 si e a 8q24, is associa ed wi h p os a e cance isk. Cell Cycle Geo ge . Tex. 12, 379–380 (2013). ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 6NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions 18. Chung, S. e al. Associa ion o a no el long non-coding RNA in 8q24 wi h p os a e cance suscep ibili y. Cance Sci. 102, 245–252 (2011). 19. P ensne , J. R. e al. T ansc ip ome sequencing ac oss a p os a e cance coho iden ifies PCAT-1, an unanno a ed lincRNA implica ed in disease p og ession. Na . Bio echnol. 29, 742–749 (2011). 20. Zheng, J. e al. The up- egula ion o long non-coding RNA CCAT2 indica es a poo p ognosis o p os a e cance and p omo es me as asis by a ec ing epi helial-mesenchymal ansi ion. Biochem. Biophys. Res. Commun. 480, 508–514 (2016). 21. Wong, N. e al. Up egula ion o FAM84B du ing p os a e cance p og ession. Onco a ge 8, 19218–19235 (2017). 22. Suo, G. e al. Oc 4 pseudogenes a e ansc ibed in cance s. Biochem. Biophys. Res. Commun. 337, 1047–1051 (2005). 23. Hayashi, H. e al. The OCT4 pseudogene POU5F1B is amplified and p omo es an agg essi e pheno ype in gas ic cance . Oncogene 34, 199–208 (2015). 24. Eeles, R. A. e al. Iden ifica ion o 23 new p os a e cance suscep ibili y loci using he iCOGS cus om geno yping a ay. Na . Gene . 45, 385–391 (2013). 391e1–2. 25. Delaneau, O., Ma chini, J. & Zagu y, J.-F. A linea complexi y phasing me hod o housands o genomes. Na . Me hods 9, 179–181 (2011). 26. Howie, B. N., Donnelly, P. & Ma chini, J. A flexible and accu a e geno ype impu a ion me hod o he nex gene a ion o genome-wide associa ion s udies. PLoS Gene . 5, e1000529 (2009). 27. Wille , C. J., Li, Y. & Abecasis, G. R. METAL: as and e ficien me a-analysis o genomewide associa ion scans. Bioin o ma. Ox . Engl. 26, 2190–2191 (2010). 28. Exco fie , L. & Sla kin, M. Maximum-likelihood es ima ion o molecula haplo ype equencies in a diploid popula ion. Mol. Biol. E ol. 12, 921–927 (1995). 29. Ba e , J. C., F y, B., Malle , J. & Daly, M. J. Haplo iew: analysis and isualiza ion o LD and haplo ype maps. Bioin o ma. Ox . Engl. 21, 263–265 (2005). 30. Johns, L. E. & Houls on, R. S. A sys ema ic e iew and me a-analysis o amilial p os a e cance isk. BJU In . 91, 789–794 (2003). 31. Yang, J., Lee, S. H., Godda d, M. E. & Vissche , P. M. GCTA: a ool o genome-wide complex ai analysis. Am. J. Hum. Gene . 88,76–82 (2011). 32. Zhong, H. & P en ice, R. L. Bias- educed es ima o s and confidence in e als o odds a ios in genome-wide associa ion s udies. Bios a . Ox . Engl. 9, 621–634 (2008). 33. Nica, A. C. e al. Candida e causal egula o y e ec s by in eg a ion o exp ession QTLs wi h complex ai gene ic associa ions. PLoS Gene . 6, e1000895 (2010). Acknowledgemen s Geno yping o he OncoA ay was unded by he US Na ional Ins i u es o Heal h (NIH) [U19 CA 148537 o ELucida ing Loci In ol ed in P os a e Cance SuscEp ibili y (ELLIPSE) p ojec and X01HG007492 o he Cen e o Inhe i ed Disease Resea ch (CIDR) unde con ac numbe HHSN268201200008I]. Addi ional analy ic suppo was p o ided by NIH NCI U01 CA188392 (PI: Schumache ). The PRACTICAL conso ium (h p://p ac ical.ic .ac.uk/) was suppo ed by Cance Resea ch UK G an s C5047/A7357, C1287/A10118, C1287/A16563, C5047/A3354, C5047/A10692, C16913/A6135, Eu - opean Commission's Se en h F amewo k P og amme g an ag eemen n° 223175 (HEALTH-F2-2009-223175), and The Na ional Ins i u e o Heal h (NIH) Cance Pos - Cance GWAS ini ia i e g an : No. 1 U19 CA 148537-01 ( he GAME-ON ini ia i e). We wish o hank all GWAS s udy g oups con ibu ing o he da a se om which his s udy was conduc ed: OncoA ay; iCOGS; The PRACTICAL (P os a e Cance Associa ion G oup o In es iga e Cance -Associa ed Al e a ions in he Genome) Conso ium; and The GAME-ON/ELLIPSE Conso ium. De ailed acknowledgemen s and unding in o - ma ion o all GWAS s udy g oups and om all he indi idual s udies in ol ed in he PRACTICAL Conso ium a e included in Supplemen a y No e 1. We would also like o hank he ollowing o unding suppo : The Ins i u e o Cance Resea ch and The E e yman Campaign, The P os a e Cance Resea ch Founda ion, P os a e Resea ch Campaign UK (now P os a e Ac ion), The O chid Cance Appeal, The Na ional Cance Resea ch Ne wo k UK, The Na ional Cance Resea ch Ins i u e (NCRI) UK. We a e g a e ul o suppo o NIHR unding o he NIHR Biomedical Resea ch Cen e a The Ins i u e o Cance Resea ch and The Royal Ma sden NHS Founda ion T us . Au ho con ibu ions M.M. and E.J.S. con ibu ed equally o his wo k. R.A.E., Z.K.-J., D.V.C., and C.A.H. join ly supe ised his wo k. T.D. con ibu ed wi h JAM analysis. M.N.B. con ibu ed wi h FRR analysis. K.W. con ibu ed wi h o wa d and backwa d s epwise selec ion. X.S. con ibu ed wi h co e age analysis. A.A.A.O., F.R.S., S.A.I., K.G., S.B., S.I.B., D.A., S.K., K.M., V.L.S., S.M.G., C.M.T., J.B., J.C., H.G., N.P., J.S., A.W., C.W., L.Mu., P.K., G.C.-T., K.D.S., L.Ma., E.M.G., S.S.S., D.E.N., F.C.H., J.L.D., R.C.T., R.J.H., B.R., Y.-J.L., G.G.G., A.S.K., A.V., J.T.B., M.K., K.L.P., J.Y.P., J.L.S., C.C., B.G.N., H.B., C.M., J.K., M.R.T., S.L.N., K.D.R., A.R., L.F.N., D.L., R.K., N.U., F.C., P.A.T., M.G.D., M.J.R., F.M., K.-T.K., L.A.C.-A., H.P., S.N.T., D.J.S., The PRACTICAL Conso ium, F.W., S.J.C., and D.F.E. we e in ol ed in sample and da a collec ion. Addi ional in o ma ion Supplemen a y In o ma ion accompanies his pape a h ps://doi.o g/10.1038/s41467- 018-06863-1. Compe ing in e es s: The au ho s decla e no compe ing in e es s. Rep in s and pe mission in o ma ion is a ailable online a h p://npg.na u e.com/ ep in sandpe missions/ Publishe 's no e: Sp inge Na u e emains neu al wi h ega d o ju isdic ional claims in published maps and ins i u ional a filia ions. Open Access This a icle is licensed unde a C ea i e Commons A ibu ion 4.0 In e na ional License, which pe mi s use, sha ing, adap a ion, dis ibu ion and ep oduc ion in any medium o o ma , as long as you gi e app op ia e c edi o he o iginal au ho (s) and he sou ce, p o ide a link o he C ea i e Commons license, and indica e i changes we e made. The images o o he hi d pa y ma e ial in his a icle a e included in he a icle’s C ea i e Commons license, unless indica ed o he wise in a c edi line o he ma e ial. I ma e ial is no included in he a icle’s C ea i e Commons license and you in ended use is no pe mi ed by s a u o y egula ion o exceeds he pe mi ed use, you will need o ob ain pe mission di ec ly om he copy igh holde . To iew a copy o his license, isi h p://c ea i ecommons.o g/ licenses/by/4.0/. © The Au ho (s) 2018 Ma co Ma ejcic 1 , Edwa d J. Saunde s 2 , Tokhi Dadae 2 , Ma k N. B ook 2 , Kan Wang 1 , Xin Sheng 1 , Ali Amin Al Olama 3,4 , F ed ick R. Schumache 5,6 , Sue A. Ingles 1 , Ko eela Go indasami 2 , Sa a Benlloch 2,3 , Sonja I. Be nd 7 , Deme ius Albanes 7 , S ella Kou os 7 , Kenne h Mui 8,9 , Vic o ia L. S e ens 10 , Susan M. Gaps u 10 , Ca he ine M. Tangen 11 , Jyo sna Ba a 12,13 , Judi h Clemen s 12,13 , Hen ik G onbe g 14 , No a Pashayan 15,16 , Johanna Schleu ke 17,18,19 , Alicja Wolk 20 , Ca ha ine Wes 21 , Lo elei Mucci 22 , Pe e K a 23 , Gé aldine Cancel-Tassin 24,25 , Ka ina D. So ensen 26,27 , Lo ise Maehle 28 , Eli M. G indedal 28 , Sa a S. S om 29 , Da id E. Neal 30,31 , F eddie C. Hamdy 32 , Jenny L. Dono an 33 , Ru h C. T a is 34 , Robe J. Hamil on 35 , Ba y Rosens ein 36,37 , Yong-Jie Lu 38 , G aham G. Giles 39,40 , Adam S. Kibel 41 , Ana Vega 42 , Jeane e T. Bensen 43 , Manolis Koge inas 44,45,46,47 , Ka h yn L. Penney 48 , Jong Y. Pa k 49 , Jane L. S an o d 50,51 , Ceza y Cybulski 52 , Bø ge G. No des gaa d 53,54 , He mann B enne 55,56,57 , Ch is iane Maie 58 , Je i Kim 59 , NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 ARTICLE NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions 7 Manuel R. Teixei a 60,61 , Susan L. Neuhausen 62 , Kim De Ruyck 63 , Azad Razack 64 , Lisa F. Newcomb 50,65 , Da o Lessel 66 , Radka Kane a 67 , Nawaid Usmani 68,69 , F ank Claessens 70 , Paul A. Townsend 71 , Manuela G. Dominguez 72,73 , Monique J. Roobol 74 , Flo ence Menegaux 75 , Kay-Tee Khaw 76 , Lisa A. Cannon-Alb igh 77,78 , Ha de Pandha 79 , S ephen N. Thibodeau 80 , Daniel J. Schaid 81 , The PRACTICAL Conso ium, F ed ik Wiklund 14 , S ephen J. Chanock 7 , Douglas F. Eas on 3,15 , Rosalind A. Eeles 2,82 , Zsofia Ko e-Ja ai 2 , Da id V. Con i 1 & Ch is ophe A. Haiman 1 1 Depa men o P e en i e Medicine, Keck School o Medicine, Uni e si y o Sou he n Cali o nia/No is Comp ehensi e Cance Cen e , Los Angeles, CA 90033, USA. 2 The Ins i u e o Cance Resea ch, London SW7 3RP, UK. 3 Cen e o Cance Gene ic Epidemiology, Depa men o Public Heal h and P ima y Ca e, S angeways Resea ch Labo a o y, Uni e si y o Camb idge, Camb idge CB1 8RN, UK. 4 Depa men o Clinical Neu osciences, Uni e si y o Camb idge, Camb idge CB2 0QQ, UK. 5 Depa men o Popula ion and Quan i a i e Heal h Sciences, Case Wes e n Rese e Uni e si y, Cle eland, OH 44106-7219, USA. 6 Seidman Cance Cen e , Uni e si y Hospi als, Cle eland, OH 44106, USA. 7 Di ision o Cance Epidemiology and Gene ics, Na ional Cance Ins i u e, NIH, Be hesda, MD 20892, USA. 8 Ins i u e o Popula ion Heal h, Uni e si y o Manches e , Manches e M13 9PL, UK. 9 Wa wick Medical School, Uni e si y o Wa wick, Co en y CV4 7AL, UK. 10 Epidemiology Resea ch P og am, Ame ican Cance Socie y, 250 Williams S ee , A lan a, GA 30303, USA. 11 SWOG S a is ical Cen e , F ed Hu chinson Cance Resea ch Cen e , Sea le, WA 98109, USA. 12 Aus alian P os a e Cance Resea ch Cen e-Qld, Ins i u e o Heal h and Biomedical Inno a ion and School o Biomedical Science, Queensland Uni e si y o Technology, B isbane, QLD 4059, Aus alia. 13 T ansla ional Resea ch Ins i u e, B isbane, QLD 4102, Aus alia. 14 Depa men o Medical Epidemiology and Bios a is ics, Ka olinska Ins i u e, SE-171 77 S ockholm, Sweden. 15 Cen e o Cance Gene ic Epidemiology, Depa men o Oncology, S angeways Resea ch Labo a o y, Uni e si y o Camb idge, Camb idge CB1 8RN, UK. 16 Depa men o Applied Heal h Resea ch, Uni e si y College London, London WC1E 7HB, UK. 17 Depa men o Medical Biochemis y and Gene ics, Ins i u e o Biomedicine, Uni e si y o Tu ku, FI-20014 Tu ku, Finland. 18 Tyks Mic obiology and Gene ics, Depa men o Medical Gene ics, Tu ku Uni e si y Hospi al, 20521 Tu ku, Finland. 19 BioMediTech, Uni e si y o Tampe e, 33520 Tampe e, Finland. 20 Di ision o Nu i ional Epidemiology, Ins i u e o En i onmen al Medicine, Ka olinska Ins i u e , SE-171 77 S ockholm, Sweden. 21 Di ision o Cance Sciences, Manches e Academic Heal h Science Cen e, Radio he apy Rela ed Resea ch, Manches e NIHR Biomedical Resea ch Cen e, The Ch is ie Hospi al NHS Founda ion T us , Uni e si y o Manches e , Manches e M13 9PL, UK. 22 Depa men o Epidemiology, Ha a d School o Public Heal h, Bos on, MA 02115, USA. 23 P og am in Gene ic Epidemiology and S a is ical Gene ics, Depa men o Epidemiology, Ha a d T.H. Chan School o Public Heal h, Bos on, MA 02115, USA. 24 GRC N°5 ONCOTYPE-URO, UPMC Uni Pa is 06, Tenon Hospi al, F-75020 Pa is, F ance. 25 CeRePP, Tenon Hospi al, F-75020 Pa is, F ance. 26 Depa men o Molecula Medicine, Aa hus Uni e si y Hospi al, 8200 Aa hus N, Denma k. 27 Depa men o Clinical Medicine, Aa hus Uni e si y, 8200 Aa hus N, Denma k. 28 Depa men o Medical Gene ics, Oslo Uni e si y Hospi al, 0424 Oslo, No way. 29 Depa men o Epidemiology, The Uni e si y o Texas MD Ande son Cance Cen e , Hous on, TX 77030, USA. 30 Depa men o Oncology, Addenb ooke’s Hospi al, Uni e si y o Camb idge, Camb idge CB2 0QQ, UK. 31 Cance Resea ch UK Camb idge Resea ch Ins i u e, Li Ka Shing Cen e, Camb idge CB2 0RE, UK. 32 Nu field Depa men o Su gical Sciences, Uni e si y o Ox o d, Ox o d OX1 2JD, UK. 33 School o Social and Communi y Medicine, Uni e si y o B is ol, Canynge Hall, 39 Wha ley Road, B is ol BS8 2PS, UK. 34 Cance Epidemiology, Nu field Depa men o Popula ion Heal h, Uni e si y o Ox o d, Ox o d OX3 7LF, UK. 35 Depa men o Su gical Oncology, P incess Ma ga e Cance Cen e, To on o, ON M5G 2M9, Canada. 36 Depa men o Radia ion Oncology, Icahn School o Medicine a Moun Sinai, New Yo k, NY 10029, USA. 37 Depa men o Gene ics and Genomic Sciences, Icahn School o Medicine a Moun Sinai, New Yo k, NY 10029-5674, USA. 38 Cen e o Molecula Oncology, John Vane Science Cen e, Ba s Cance Ins i u e, Queen Ma y Uni e si y o London, London EC1M 6BQ, UK. 39 Cance Epidemiology & In elligence Di ision, Cance Council Vic o ia, Melbou ne, VIC 3004, Aus alia. 40 Cen e o Epidemiology and Bios a is ics, Melbou ne School o Popula ion and Global Heal h, The Uni e si y o Melbou ne, Melbou ne, VIC 3010, Aus alia. 41 Di ision o U ologic Su ge y, B igham and Womens Hospi al, Bos on, MA 02115, USA. 42 Fundación Pública Galega de Medicina Xenómica-SERGAS, G upo de Medicina Xenómica, CIBERER, IDIS, 15706 San iago de Compos ela, Spain. 43 Depa men o Epidemiology, Gillings School o Global Public Heal h, Uni e si y o No h Ca olina, Columbia, SC 29208, USA. 44 Cen e o Resea ch in En i onmen al Epidemiology (CREAL), Ba celona Ins i u e o Global Heal h (ISGlobal), 08003 Ba celona, Spain. 45 CIBER Epidemiología y Salud Pública (CIBERESP), 28029 Mad id, Spain. 46 IMIM (Hospi al del Ma Resea ch Ins i u e), 08003 Ba celona, Spain. 47 Uni e si a Pompeu Fab a (UPF), 08002 Ba celona, Spain. 48 Channing Di ision o Ne wo k Medicine, Depa men o Medicine, B igham and Women’s Hospi al/ Ha a d Medical School, Bos on, MA 02184, USA. 49 Depa men o Cance Epidemiology, Mo fi Cance Cen e , Tampa, FL 33612, USA. 50 Di ision o Public Heal h Sciences, F ed Hu chinson Cance Resea ch Cen e , Sea le, WA 98109-1024, USA. 51 Depa men o Epidemiology, School o Public Heal h, Uni e si y o Washing on, Sea le, WA 98195, USA. 52 In e na ional He edi a y Cance Cen e , Depa men o Gene ics and Pa hology, Pome anian Medical Uni e si y, 70-115 Szczecin, Poland. 53 Facul y o Heal h and Medical Sciences, Uni e si y o Copenhagen, 2200 Copenhagen, Denma k. 54 Depa men o Clinical Biochemis y, He le and Gen o e Hospi al, Copenhagen Uni e si y Hospi al, He le , 2200 Copenhagen, Denma k. 55 Di ision o Clinical Epidemiology and Aging Resea ch, Ge man Cance Resea ch Cen e (DKFZ), D-69120 Heidelbe g, Ge many. 56 Ge man Cance Conso ium (DKTK), Ge man Cance Resea ch Cen e (DKFZ), D-69120 Heidelbe g, Ge many. 57 Di ision o P e en i e Oncology, Ge man Cance Resea ch Cen e (DKFZ) and Na ional Cen e o Tumo Diseases (NCT), 69120 Heidelbe g, Ge many. 58 Ins i u e o Human Gene ics, Uni e si y Hospi al Ulm, 89075 Ulm, Ge many. 59 Depa men o Geni ou ina y Medical Oncology, The Uni e si y o Texas MD Ande son Cance Cen e , Hous on, TX 77030, USA. 60 Depa men o Gene ics, Po uguese Oncology Ins i u e o Po o, 4200-072 Po o, Po ugal. 61 Biomedical Sciences Ins i u e (ICBAS), Uni e si y o Po o, 4050-313 Po o, Po ugal. 62 Depa men o Popula ion Sciences, Beckman Resea ch Ins i u e o he Ci y o Hope, Dua e, CA 91010, USA. 63 Ghen Uni e si y, Facul y o Medicine and Heal h Sciences, Basic Medical Sciences, B-9000 Gen , Belgium. 64 Depa men o Su ge y, Facul y o Medicine, Uni e si y o Malaya, 50603 Kuala Lumpu , Malaysia. 65 Depa men o U ology, Uni e si y o Washing on, Sea le, WA 98195, USA. 66 Ins i u e o Human Gene ics, Uni e si y Medical Cen e Hambu g- Eppendo , D-20246 Hambu g, Ge many. 67 Molecula Medicine Cen e , Depa men o Medical Chemis y and Biochemis y, Medical Uni e si y o Sofia, 1431 Sofia, Bulga ia. 68 Depa men o Oncology, C oss Cance Ins i u e, Uni e si y o Albe a, Edmon on, AB T6G 1Z2, Canada. 69 Di ision o Radia ion Oncology, C oss Cance Ins i u e, Uni e si y o Albe a, Edmon on, AB T6G 1Z2, Canada. 70 Molecula Endoc inology Labo a o y, Depa men o Cellula and Molecula Medicine, KU Leu en, BE-3000 Leu en, Belgium. 71 Manches e Cance Resea ch Cen e, Facul y o Biology Medicine and Heal h, Manches e Academic Heal h Science Cen e, NIHR Manches e Biomedical Resea ch Cen e, Heal h Inno a ion Manches e , Uni e si y o Manches e , Manches e M13 9WL, UK. 72 Genomic Medicine G oup, Galician Founda ion o Genomic Medicine, Ins i u o de In es igacion Sani a ia de San iago de Compos ela (IDIS), Complejo Hospi ala io Uni e si a io de San iago, Se icio Galego de Saúde, SERGAS, ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 8NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions 15706 San iago de Compos ela, Spain. 73 Moo es Cance Cen e , Uni e si y o Cali o nia San Diego, La Jolla, CA 92037, USA. 74 Depa men o U ology, E asmus Uni e si y Medical Cen e , 3015 CE Ro e dam, The Ne he lands. 75 Cance and En i onmen G oup, Cen e o Resea ch in Epidemiology and Popula ion Heal h (CESP), INSERM, Uni e si y Pa is-Sud, Uni e si y Pa is-Saclay, 94807 Villejui Cédex, F ance. 76 Clinical Ge on ology Uni , Uni e si y o Camb idge, Camb idge CB2 2QQ, UK. 77 Di ision o Gene ic Epidemiology, Depa men o Medicine, Uni e si y o U ah School o Medicine, Sal Lake Ci y, UT 84112, USA. 78 Geo ge E. Wahlen Depa men o Ve e ans A ai s Medical Cen e , Sal Lake Ci y, UT 84148, USA. 79 The Uni e si y o Su ey, Guild o d, Su ey GU2 7XH, UK. 80 Depa men o Labo a o y Medicine and Pa hology, Mayo Clinic, Roches e , MN 55905, USA. 81 Di ision o Biomedical S a is ics and In o ma ics, Mayo Clinic, Roches e , MN 55905, USA. 82 Royal Ma sden NHS Founda ion T us , London SW3 6JJ, UK. These au ho s con ibu ed equally: Ma co Ma ejcic, Edwa d J. Saunde s. These au ho s join ly supe ised his wo k: Rosalind A. Eeles, Zsofia Ko e‑Ja ai, Da id V. Con i, Ch is ophe A. Haiman. A ull lis o conso ium membe s appea s a he end o he pape . The PRACTICAL (P os a e Cance Associa ion G oup o In es iga e Cance -Associa ed Al e a ions in he Genome) Conso ium B ian E. Hende son 1 , Ma iana C. S e n 1 , Alison Thwai es 2 , Michelle Guy 2 , Ian Whi mo e 2 , Angela Mo gan 2 , Cy il Fishe 2 , S e e Hazel 2 , Naomi Li ni 2 , Ma ga e Cook 3 , Lau a Fachal 3,42 , S ephanie Weins ein 7 , Lau a E. Beane F eeman 7 , Robe N. Hoo e 7 , Mi chell J. Machiela 7 , A i aya Lopha ananon 8,9 , B ian D. Ca e 10 , Phyllis Goodman 11 , Lei e Moya 12,13 , S ilakshmi S ini asan 12,13 , Ma y-Anne Kedda 12,13 , T ina Yeadon 12,13 , Allison Ecke 12,13 , Ma in Eklund 14 , Ca in Ca alli-Bjoe kman 14 , Alison M. Dunning 15 , Csilla Sipeky 17 , Niclas Hakansson 20 , Rebecca Ellio 21 , Ha deep Ranu 22 , Edwa d Gio annucci 22 , Cons ance Tu man 23 , Da id J. Hun e 23 , Oli ie Cusseno 24,25 , To ben Falck O n o 26,27 , A hene Lane 33 , Sa ah J. Lewis 33 , Michael Da is 33 , Tim J. Key 34 , Paul B own 35 , Gi ish S. Kulka ni 35 , Alexand e R. Zlo a 35 , Neil E. Fleshne 35 , An onio Finelli 35 , Xueying Mao 38 , Jacek Ma zec 38 , Robe J. MacInnis 39,40 , Roge Milne 39,40 , John L. Hoppe 40 , Miguel Aguado 42 , Ma iona Bus aman e 44 , Gemma Cas año-Vinyals 44,45,46,47 , Es he G acia-La edan 44,45,46,47 , Lluís Cecchini 46 , Mei S amp e 48 , Jing Ma 48 , Thomas A. Selle s 49 , Milan S. Geybels 49 , Hyun Pa k 49 , Babu Zacha iah 49 , Suzanne Kolb 50 , Dominika Wokolo czyk 52 , Jan Lubinski 52 , Wojciech Kluzniak 52 , Sune F. Nielsen 53,54 , Ma en Weishe 54 , Ka a ina Cuk 55 , Wal he Vogel 58 , Manuel Luedeke 58 , Ch is ophe J. Logo he is 59 , Paula Paulo 60 , Ma a Ca doso 60 ,Sofia Maia 60 , Ma ia P. Sil a 60 , Linda S eele 62 , Yuan Chun Ding 62 , Ge De Mee lee 63 ,Sofie De Langhe 63 , Hube Thie ens 63 , Jasmine Lim 64 , Meng H. Tan 64 , Aik T. Ong 64 , Daniel W. Lin 50,65 , Da ina Kachako a 67 , A anaska Mi ko a 67 , Vanio Mi e 67 , Ma hew Pa liamen 68,69 , Guido Jens e 74 , Ch is ophe Bangma 74 , F.H. Sch ode 74 , Thé èse T uong 75 , Y es Akoli Koudou 75 , Agnieszka Michael 79 , And zej Kie zek 79 , Ami Ka lsson 79 , Michael B oms 79 , Huihai Wu 79 , Clai e Aukim-Has ie 79 , Lo i Tillmans 80 , Shaun Riska 80 , Shannon K. McDonnell 81 , Da id Dea naley 2,82 , Amanda Spu dle 83 , Robe Ga dine 84,85 , Vanessa Hayes 86 , Lisa Bu le 87 , Renea Taylo 88 , Melissa Papa gi is 88 , Pamela Saunde s 89 , Paula Kujala 90 , Ki si Talala 91 , Kimmo Taa i 92 , Sø en Ben zen 93 , Belynda Hicks 94 , Au elie Vog 94 , Amy Hu chinson 95 , Angela Cox 96 , Anne Geo ge 97 , An s Toi 98 , And ew E ans 99 , Theodo us H. an de Kwas 99 , Takashi Imai 100 , Shi o Sai o 101 , Shan-Chao Zhao 102 , Guoping Ren 103 , Yangling Zhang 103 , Yongwei Yu 104 , Yudong Wu 105 ,JiWu 106 , Bo Zhou 107 , John Pede sen 108 , Ramón Loba o-Bus o 109 , José Manuel Ruiz-Dominguez 110 , Lou des Mengual 111,112 , An onio Alca az 113 , Julio Pow-Sang 114 , Ka hleen He komme 115 , Aleksand ina Vlaho a 116 , Tihomi Diko 116 , S e lana Ch is o a 116 , Angel Ca acedo 42,117,118 , B igi e T e a e 119 , Xa ie Rebilla d 120 , Clai e Mulo 121 , Jan Adol sson 122,123 , Pa S a in 124,125 , Jan-E ik Johansson 126 , Richa d M. Ma in 33,127,128 , Ian M. Thompson J . 129 , Suzanne Chambe s 130,131 , Joanne Ai ken 130,131 , Lisa Ho a h 132,133 , Anne-Ma ee Haynes 86,133 , Wayne Tilley 134 , Gail Risb idge 135,136 , Ma kus Aly 14,137 , Tobias No ds öm 14,138 , Paul Pha oah 3,139 , Teu o L.J. Tammela 140 , Teemu Mu ola 140,141 , Anssi Au inen 142 , Neil Bu ne 143 , Gill Ba ne 143 , Ge ald And iole 144 , Aleksand a Klim 144 , Be ina F. D ake 144 , Michael Bo e 27,145 , Sa ah Ke ns 146 , Ha y Os e 147 , Hong-Wei Zhang 148 , Guangwen Cao 148 , Ji Lin 148 , Jin Ling 148 , Meiling Li 148 , Ninghan Feng 149 , Jie Li 150 , Weiyang He 150 , Xin Guo 150,151 , Zan Sun 151 , NATURE COMMUNICATIONS | DOI: 10.1038/s41467-018-06863-1 ARTICLE NATURE COMMUNICATIONS | (2018) 9:4616 | DOI: 10.1038/s41467-018-06863-1 | www.na u e.com/na u ecommunica ions 9