Genetic correlation between multiple myeloma and chronic lymphocytic leukaemia provides evidence for shared aetiology

<jats:title>Abstract</jats:title><jats:p>The clustering of different types of B-cell malignancies in families raises the possibility of shared aetiology. To examine this, we performed cross-trait linkage disequilibrium (LD)-score regression of multiple myeloma (MM) and chronic lymphocytic leukaemia (CLL) genome-wide association study (GWAS) data sets, totalling 11,734 cases and 29,468 controls. A significant genetic correlation between these two B-cell malignancies was shown (<jats:italic>R</jats:italic><jats:sub>g</jats:sub> = 0.4, <jats:italic>P</jats:italic> = 0.0046). Furthermore, four of the 45 known CLL risk loci were shown to associate with MM risk and five of the 23 known MM risk loci associate with CLL risk. By integrating eQTL, Hi-C and ChIP-seq data, we show that these pleiotropic risk loci are enriched for B-cell regulatory elements and implicate B-cell developmental genes. These data identify shared biological pathways influencing the development of CLL and, MM and further our understanding of the aetiological basis of these B-cell malignancies.</jats:p>

section_content: Chronic lymphocytic leukaemia (CLL) and multiple myeloma (MM) are both B-cell malignancies, which arise from the clonal expansion of progenitor cells at different stages of B-cell maturity [1] [2] [3]. Evidence for inherited predisposition to CLL and MM comes from the six-and twofold increased risk of the respective diseases seen in relatives of patients 4. 

Recent genome-wide association studies (GWAS) have transformed our understanding of genetic susceptibility to the B-cell malignancies, identifying 45 CLL [5] [6] [7] [8] and 23 MM risk loci [9] [10] [11] [12] [13]. Furthermore, statistical modelling of GWAS data indicates that common genetic variation is likely to account for 34% of CLL and 15% of MM heritability 6, 14. 

Epidemiological observations on familial cancer risks across the different B-cell malignancies suggest an element of shared inherited susceptibility, especially between CLL and MM 4. 

Linkage disequilibrium (LD) score regression is a method which exploits the feature of a test statistic for a given single nucleotide polymorphism (SNP), whereby that test statistic will incorporate the effects of correlated SNPs 15. Conventional LD score regression regresses trait χ 2 statistics against the LD score for a given SNP, with the coefficient of the regression line providing an estimate of trait heritability. This method can be modified by instead regressing the product of SNP Z-scores from two traits against the SNP LD score, with the slope providing an estimate of genetic covariance between the two traits 16. This method can be applied to summary statistics, is not biased by sample overlap and does not require multiple traits to be measured for each individual. 

By analysis of GWAS data for MM and CLL and applying cross-trait LD score regression, we have been able to demonstrate a positive genetic correlation between CLL and MM. We find evidence of shared genetic susceptibility at 10 known risk loci and by integrating promoter capture Hi-C (PCHi-C) data, ChIP-seq and expression data we provide insight into the shared biological basis of CLL and MM.
section_name: Introduction

section_name: Methods

section_content: The data from six previously reported MM GWAS [9] [10] [11] [12] [13] are summarized in Supplementary Table 1. All these studies were based on individuals of European ancestry and comprised: Oncoarray-GWAS (878 cases 7054 controls) UK-GWAS (2282 cases, 5197 controls), Swedish-GWAS (1714 cases, 10,391 controls), German-GWAS (1508 cases, 2107 controls), Netherlands-GWAS (555 cases, 2669 controls) and US-GWAS (780 cases, 1857 controls). 

The data from three previously reported CLL GWAS 8- 13 are summarized in Supplementary The diagnosis of MM (ICD-10 C90. 0) in all cases was established in accordance with World Health Organization guidelines. All samples from patients for genotyping were obtained before treatment or at presentation. The diagnosis of CLL (ICD-10-CM C91. 10, ICD-O M9823/3 and 9670/3) was established in accordance with the International Workshop on Chronic Lymphocytic Leukaemia guidelines.
section_name: GWAS data sets

section_content: Standard quality-control measures were applied to the GWAS 17. Specifically, individuals with low SNP call rate (<95%) as well as individuals evaluated to be of non-European ancestry (using the HapMap version 2 CEU, JPT/CHB and YRI populations as a reference) were excluded. For apparent first-degree relative pairs, we excluded the control from a case-control pair; otherwise, we excluded the individual with the lower call rate. SNPs with a call rate <95% were excluded as were those with a MAF <0. 01 or displaying significant deviation from Hardy-Weinberg equilibrium (P < 10 -5 ). GWAS data were imputed to >10 million SNPs using IMPUTE2 v4 (for CLL) and IMPUTE2 v2. 3 (for MM) software in conjunction with a merged reference panel consisting of data from 1000 Genomes Project 18 (phase 1 integrated release 3 March 2012) and UK10K 19. Genotypes were aligned to the positive strand in both imputation and genotyping. We imposed predefined thresholds for imputation quality to retain potential risk variants with MAF >0. 01 for validation. Poorly imputed SNPs with an information measure <0. 80 were excluded. Tests of association between imputed SNPs and MM were performed under an additive model in SNPTESTv2. 5 20. The adequacy of the case-control matching and possibility of differential genotyping of cases and controls was evaluated using a Q-Q plot of test statistics. The inflation λ was based on the 90% least-significant SNPs and assessment of λ 1000. Details of SNP QC are provided in Supplementary Table 3 and 4. Four principal components, generated using common SNPs, were included to limit the effects of cryptic population stratification in the US-CLL data set. Eigenvectors for the GWAS data sets were inferred using smartpca (part of EIGENSOFT) by merging cases and controls with phase II HapMap samples.
section_name: Quality control

section_content: Meta-analyses were performed using the fixed-effects inverse-variance method using META v1. 6 21. Cochran's Q-statistic to test for heterogeneity and the I 2 statistic to quantify the proportion of the total variation due to heterogeneity was calculated.
section_name: Meta-analysis

section_content: To investigate genetic correlation between MM and CLL, we implemented cross-trait LD score regression by Bulik-Sullivan et al. 16. Using summary statistics from the GWAS meta-analysis we implemented filters as recommended by the authors 16. Specifically, filtering SNPs to INFO >0. 9, MAF >0. 01, and harmonizing to Hap Map3 SNPs with 1000 Genomes EUR MAF >0. 05, removing indels and structural variants, removing strandambiguous SNPs and removing SNPs where alleles did not match those in 1000 Genomes. This was performed by running the munge-sumstats. pr script included with ldsc. We ran ldsc. py, part of the ldsc package, excluding the HLA region. We report heritability estimates on the observed scale. There is no distinction between observed and liability scale genetic correlation for case/control traits 16.
section_name: LD score regression

section_content: To identify pleiotropic risk loci, that is genetic loci that influence two traits, we identified SNPs previously reported to be associated with each disease at genomewide significance (P < 5 × 10 -8 ), as well as highly correlated variants (r 2 > 0. 8) at the 45 and 23 known risk loci for CLL and MM, respectively. Within these correlated variant sets at each locus, we determined how many of the CLL susceptibility loci were associated with MM at region-wide significance after Bonferroni correction for multiple testing (i. e. P adj < 0. 05/45). We then repeated the process, examining MM susceptibility SNPs in CLL, applying a significance level of P adj < 0. 05/23. A full list of results is summarized in Supplementary Data File 1 and 2.
section_name: Shared risk loci

section_content: A variation of LD score regression, namely stratified LD score regression, can be used to partition heritability according to different genomic categories. For both MM and CLL we applied stratified LD score regression across the baseline model used in Finucane et al. 22. We plotted the enrichment of functional categories for each diseasethis is defined as proportion heritability divided by the total heritability. We excluded from our plot additional flanking regions around each functional category, which authors designed to allow observation of enrichment of SNP heritability in intermediary regions. A plot of the results is found in Supplementary Figure 1.
section_name: Partitioned heritability

section_content: To examine enrichment in specific histone mark binding across shared risk loci, we adapted the method of Cowper-Sal lari et al. 23. Briefly, for each risk locus, a region of strong LD (defined as r 2 > 0. 8 and D′ > 0. 8) was determined, and these SNPs were considered the associated variant set (AVS). Publically available ChIP-seq data for 6 histone marks from naive B cells was downloaded from Blueprint Epigenome Project 24. For each mark, the overlap of the SNPs in the AVS and the binding sites was assessed to generate a mapping tally. A null distribution was produced by randomly selecting SNPs with the same characteristics as the risk-associated SNPs, and the null mapping tally calculated. This process was repeated 10,000 times, and P-values calculated as the proportion of permutations where null mapping tally was greater or equal to the AVS mapping tally. An enrichment score was calculated by normalizing the tallies to the median of the null distribution. Thus, the enrichment score is the number of standard deviations of the AVS mapping tally from the median of the null distribution tallies. An enrichment plot for naive B cells is shown in Supplementary Figure 2.
section_name: Variant set enrichment

section_content: We considered chromatin mark overlap enrichment for genome-wide significant loci in different cell types using the methodology of Trynka et al. 25. This approach scores GWAS SNPs based on proximity to chromatin mark and fold-enrichment of respective chromatin mark, assessing significance using a tissue-specific permutation method. We obtained chip-seq data for H3K4me3 from primary blood cells and CLL samples downloaded from Blueprint Epigenome project 24. In addition, we included in our analysis 4 MM cell lines-KMS11, JJN3, MM1-S and L363 processed as previously described 26. A heat map of results is shown in Supplementary Figure 3. eQTL eQTL analyses were performed using publicly available whole-blood data downloaded from GTeX 27. The relationship between SNP genotype and gene expression we carried out using Summary-data-based Mendelian Randomization (SMR) analysis as per Zhu et al. 28. Briefly, if b xy is the effect size of x (gene expression) on y (slope of y regressed on the genetic value of x), b zx is the effect of z on x, and b zy be the effect of z on y, b xy (b zy /b zx ) is the effect of x on y. To distinguish pleiotropy from linkage where the top associated cis-eQTL is in LD with two causal variants, one affecting gene expression the other affecting trait we tested for heterogeneity in dependent instruments (HEIDI), using multiple SNPs in each cis-eQTL region. Under the hypothesis of pleiotropy b xy values for SNPs in LD with the causal variant should be identical. For each probe that passed significance threshold for the SMR test, we tested the heterogeneity in the b xy values estimated for multiple SNPs in the cis-eQTL region using HEIDI. 

GWAS summary statistics files were generated from the meta-analysis. For the disease discovery GWAS, we set a threshold for the SMR test of P SMR < 2. 5 × 10 -5 corresponding to a Bonferroni correction for the number of probes which demonstrated an association in the SMR test. For all genes passing this threshold we generated plots of the eQTL and GWAS associations at the locus, as well as plots of GWAS and eQTL effect sizes (i. e. input for the HEIDI heterogeneity test). HEIDI test P-values <0. 05 were considered as reflective of heterogeneity. This threshold is, however, conservative for gene discovery because it retains fewer genes than when correcting for multiple testing. SMR plots for significant eQTLs are shown in Supplementary Figures 4, 5 and a summary of results are shown in Supplementary Table 5.
section_name: Cell-type-specific analyses

section_name: Results

section_content: We performed cross-trait LD-score regression using summary statistics from two recent GWAS meta-analyses based on 7717 MM cases and 21,587 controls and 4017 CLL cases and 7881 controls (Fig. 1, Supplementary Table 1 2 3 4 ). While these data sets have been previously subject to quality control (QC) [5] [6] [7] [9] [10] [11] [12] for the current analysis we implemented additional filtering steps as per Bulik-Sullivan et al. 16, resulting in 1,055,728 harmonized SNPs between the two data sets. Heri ability estimates from cross-trait LD score regression of 9. 2 (±1. 8%) and 22 (±5. 9%) were comparable with previous estimates for MM 14 and CLL 6. LD-score regression revealed a significant-positive genetic correlation between MM and CLL with an R g value of 0. 44 (P = 4. 6 × 10 -3 ).
section_name: Genetic correlation and heritability

section_content: We identified SNPs previously reported to be associated with each disease at genome-wide significance (P < 5 × 10 -8 ), as well as highly correlated variants (r 2 > 0. 8) at the 45 and 23 known risk loci for CLL and MM, respectively. To identify pleiotropic risk loci, that is genetic loci that influence two traits, we determined how many of the CLL susceptibility loci were associated with MM at regionwide significance after Bonferroni correction for multiple testing (i. e. P adj < 0. 05/45). We then repeated the process, examining MM susceptibility SNPs in CLL, applying a significance level of P adj < 0. 05/23. Of the 45 CLL risk loci, four were associated with MM (P adj < 0. 0011) while, of 23 MM risk loci, five were significantly associated in CLL (P adj < 0. 0022) (Table 1, Fig. 2 ). Correlated SNPs (r 2 > 0. 8) Fig. 1 Schematic outlining the processing of data sets used in the genetic correlation at 3q26. 2 are associated with both CLL and MM at genome-wide significance (Fig. 2 ), bringing the total number of pleiotropic loci to 10.
section_name: Identification of pleiotropic risk loci

section_content: Trynka et al. have recently shown that chromatin marks highlighting active regulatory regions overlap with phenotype-associated variants in a cell-type-specific manner 25. As H3K4me3 was shown to be the most phenotypically cell-type-specific chromatin mark, we examined cell-type specificity of the 10 pleiotropic risk loci by analysing H3K4me3 chromatin marks in normal haematopoietic cells and CLL patient samples from Blueprint, and de novo data on KMS11, MM1S, JJN3 and L363 MM cell lines. Cell types showing the strongest enrichment of risk SNPs at H3K4me3 marks included naive B cells and CD38-B cells. Notably, variants at 2q31. 1, 6p25. 3, 8q24. 21, 16q23. 1 and 22q13. 33 were enriched for H3K4me3 in naive B cells (Supplementary Figure 3 ). 

Most GWAS signals map to non-coding regions of the genome 29, 30 and influence gene expression through chromatin looping interactions 31, 32. Application of partitioned heritability analysis, stratifying across 53 genomic categories demonstrated enrichment of CLL and MM heritability in functional elements of the genome, in particular FANTOM5 enhancers (CLL and MM) transcription start sites (CLL) and 5′ untranslated region and coding regions (MM) (Supplementary Figure 1 ). Furthermore, we found significant enrichment of SNPs in the shared loci within regions of active chromatin, as indicated by the presence of H3K27ac and H3K4Me3 marks in naive B cells, supporting the principle that SNPs in shared loci influence risk through regulatory effects (Supplementary Figure 2 ). To identify target genes we analysed PCHi-C data on naive B cells from Blueprint 24. We also sought to gain insight into the possible biological mechanisms for associations by performing an expression quantitative trait locus (eQTL) analysis using mRNA expression data on blood from GTEx. Applying Summary data-based Mendelian Randomization (SMR) methodology, we tested for pleiotropy between GWAS signal and cis-eQTL for genes to identify a causal relationship. Broadly, our analysis of the shared loci groups them into those which act on a B-cell regulation and differentiation and those which underpin the distinctive biology of cancer; specifically, loci relating to genome instability, angiogenesis and dysregulated apoptosis (Supplementary Table 6 ). 

Of the shared loci, three were related to B-cell regulation. This included composite evidence at 10q23. 31, from looping interactions in naive B cells and correlation in GWAS effect size and expression, which provide evidence for two candidate genes ACTA2, encoding smooth muscle (α)-2 actin, a protein involved in cell movement and contraction of muscles 33 and FAS, a member of the TNFreceptor superfamily. FAS, has a central role in regulating the immune response through apoptosis of B cells 34, 35. At 2q31. 1, looping interactions implicated transcription factor SP3, which has been shown to influence expression of germinal centre genes, 36, 37. Variants at 6p25. 3 reside in the 3′-UTR of IRF4, which has an established role in Bcell regulation 38, 39 and MM oncogenesis 40, 41. 

Three of the 10 loci contain genes with roles in maintenance of genomic stability. Specifically, evidence from expression and PCHi-C data implicated RFWD3 at 16q23. 1. This gene encodes an E3 ubiquitin-protein ligase, which has been shown to promote progression to late stage homologous recombination through ubiquitination and timely removal of RAD51 and RPA at sites of DNA damage 42 and is necessary for replication fork restart 43. Variants in this locus demonstrated enrichment of H3K4me3 marks in two samples of naive B cells, which represents a plausible cell of disease origin. rs58618031 (7q31. 33) maps 5′ of POT1, the protection of telomeres 1 gene, which is part of the shelterin complex and functions to maintain chromosomal stability 44, 45. Variant rs1317082 at 3q26. 2 is located proximal to TERC, a gene which has been shown to influence telomere length 46. Additionally, we observed looping interactions to a number of genes at 3q26. 2 including SEC62, which has been proposed as a cancer biomarker [47] [48] [49] [50]. Intriguingly, variants at 3q26. 2 this locus have been implicated in colorectal 51, thyroid 52 and bladder 53 cancer. 

Several genes were implicated at 22q13. 33 by looping interactions for SCO2, LMF2, ODF3B, TYMP/ECGF1, NCAPH2, SYCE3 and ARSA, with TYMP/ECGF1 and SCO2 demonstrating evidence of correlation in GWAS and eQTL effect size, albeit not significant after multiple testing (P SMR = 2. 38 × 10 -4 and 3. 19 × 10 -4 ). Variants within this locus were enriched in H3K4me3 chromatin marks in both CD38-B cells and inflammatory macrophages. TYMP (alias ECGF1) encodes thymidine phosphorylase, which is often overexpressed in tumours and has been linked to angiogenesis 54, 55. A detailed study on this gene has implicated TYMP in the development of lytic bone lesions in MM, via a mechanism involving activation of PI3K/Akt signalling and increased DNMT3A expression resulting in hypermethylation of RUNX2, osterix, and IRF8 56. Furthermore, SCO2 (synthesis of cytochrome c oxidase), also mapping to this locus, has been implicated in the development of breast 57, 58, gastric 59 and leukaemia 60, through glucose metabolism reprogramming 61, a hallmark of cancer 62. Tumour suppressor, p53, regulates metabolic pathways, p53transactivated TP53-induced glycolysis (TIGAR), and regulation of apoptosis in part through SCO2 58, 59, 61. 

Finally, whereas these data were indifferent to decipher 8q24. 21, this locus has also been shown to harbour risk SNPs for other cancers, which localize within distinct LD blocks and likely reflect tissue-specific effects on cancer risk through regulation of MYC 30.
section_name: Biological inference

section_content: Our analysis provides evidence of a genetic correlation between MM and CLL. Furthermore, we have identified shared genetic susceptibility at 10 known risk loci. While requiring biological validation, integration of data from PCHi-C, chromatin mark enrichment and eQTL at shared loci has provided insight into how these loci may confer susceptibility to both CLL and MM. Applying a working hypothesis that the loci may act in pleiotropic fashion, we selected relevant cells representing a common tissue of disease origin; namely naive B cells. 

A significant genetic correlation between MM and CLL, as well as the discovery of risk loci shared between them, Fig. 2 Overlap of loci in multiple myeloma and chronic lymphocytic leukaemia. *correlated variants at 3q26. 2 had been previously published as genome wide significant in each data set prior to this analysis supports epidemiological data demonstrating elevated familial risks between these B-cell malignancies 4. Furthermore, the shared loci we identified could be broadly grouped into those containing genes related to B-cell regulation and differentiation and those containing genes involved in angiogenesis, genome stability and apoptosis, supporting the tenet that these alleles can influence aetiology of either disease. With the expansion of GWAS of the B-cell malignancies, more detailed characterisation of common underlying risk alleles and affected pathways can inform the biology of B-cell oncogenesis. 

and other individuals who participated in the project. The Genome-Wide Association Study (GWAS) of Non-Hodgkin Lymphoma (NHL) project, from which US-CLL samples were obtained, was supported by the intramural program of the Division of Cancer Epidemiology and Genetics (DCEG), National Cancer Institute (NCI), National Institutes of Health (NIH). The data sets have been accessed through the NIH database for Genotypes and Phenotypes (dbGaP) under accession # phs000801. A full list of acknowledgements can be found in supplementary note (Berndt SI et al., Nature Genet., 2013, PMID: 23770605).
section_name: Discussion

section_name: Acknowledgements

section_content: In the United Kingdom, Myeloma UK and Bloodwise (# 05001, # 06002 and # 13044 ) provided principal funding. Additional funding was provided by Cancer Research UK ( C1298/A8362 supported by the Bobby Moore Fund) and The Rosetrees Trust. M. W. is supported by funding from Mr Ralph Stockwell. A. S. is supported by a clinical fellowship from Cancer Research UK and the Royal Marsden Haematology Research Fund. This study made use of genotyping data on the 1958 Birth Cohort generated by the Wellcome Trust Sanger Institute (http://www. wtccc. org. uk). We are grateful to all investigators who contributed to NSCCG and GELCAPS, from which controls in the replication were drawn. We also thank the staff of the CTRU University of Leeds and the NCRI Haematology Clinical Studies Group. We thank the High-throughput Genomics Group at the Wellcome Trust Centre for Human Genetics (funded by Wellcome Trust grant reference 090532/Z/09/Z ) for the generation of UK myeloma Oncoarray data. The BCAC study would not have been possible without the contributions of the following: Manjeet K. Bolla, Qin Wang, Kyriaki Michailidou and Joe Dennis. BCAC is funded by Cancer Research UK ( C1287/ A10118, C1287/A16563 ). For the BBCS study, we thank Eileen Williams, Elaine Ryder-Mills and Kara Sargus. The BBCS is funded by Cancer Research UK and Breast Cancer Now and acknowledges NHS funding to the National Institute of Health Research (NIHR) Biomedical Research Centre (BRC), and the National Cancer Research Network (NCRN). We thank the participants and the investigators of EPIC (European Prospective Investigation into Cancer and Nutrition). The coordination of EPIC is financially supported by the European Commission (DG-SANCO) and the International Agency for Research on Cancer. The national cohorts are supported by: Ligue Contre le Cancer, Institut Gustave Roussy, Mutuelle Générale de l'Education Nationale, Institut National de la Santé et de la Recherche Médicale (INSERM) (France) ; German Cancer Aid, German Cancer Research Center (DKFZ), Federal Ministry of Education and Research (BMBF) (Germany) ; the Hellenic Health Foundation, the Stavros Niarchos Foundation (Greece) ; Associazione Italiana per la Ricerca sul Cancro-AIRC-Italy and National Research Council (Italy) ; Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands (The Netherlands); Health Research Fund (FIS), PI13/00061 to Granada, PI13/01162 to EPIC-Murcia, Regional Governments of Andalucía, Asturias, Basque Country, Murcia and Navarra, ISCIII RETIC ( RD06/0020 ) (Spain); Cancer Research UK ( 14136 to EPIC-Norfolk ; C570/A16491 and C8221/A19170 to EPIC-Oxford ), Medical Research Council ( 1000143 to EPIC-Norfolk, MR/M012190/1 to EPIC-Oxford) (UK). We thank the SEARCH and EPIC teams which were funded by a programme grant from Cancer Research UK [A16561] and supported by the UK NIHR BRC at the University of Cambridge. We thank Breast Cancer Now and The Institute of Cancer Research ( ICR ) for support and funding of the UKBGS, and the study participants, study staff, and the doctors, nurses and other health care providers and health information sources who have contributed to the study. We acknowledge NHS funding to the Royal Marsden/ICR NIHR BRC. UKGPCS would like to thank The Institute of Cancer Research and The Everyman Campaign for funding support. The UKGPCS acknowledges The Prostate Cancer Research Foundation, Prostate Action, The Orchid Cancer Appeal, The National Cancer Research Network UK, The National Cancer Research Institute (NCRI), the NIHR funding to the NIHR Biomedical Research data managers and consultants for their work in the UKGPCS study and urologists and other persons involved in the planning, and data collection of the CAPS study. Genotyping of the OncoArray was funded by the US National Institutes of Health (NIH) [ U19 CA 148537 for ELucidating Loci Involved in Prostate cancer SuscEptibility (ELLIPSE) project and X01HG007492 to the Center for Inherited Disease Research (CIDR) under contract number HHSN268201200008I ]. Additional analytic support was provided by NIH NCI U01 CA188392 (PI: Schumacher). The PRACTICAL consortium was supported by Cancer Research UK Grants C5047/A7357, C1287/A10118, C1287/A16563, C5047/A3354, C5047/ A10692, C16913/A6135, European Commission's Seventh Framework Programme grant agreement n°223175 ( HEALTH-F2-2009-223175 ), and The National Institute of Health (NIH) Cancer Post-Cancer GWAS initiative grant: No. 1 U19 CA 148537-01 (the GAME-ON initiative). We would also like to thank the following for funding support: The Institute of Cancer Research and The Everyman Campaign, The Prostate Cancer Research Foundation, Prostate Research Campaign UK (now Prostate Action ), The Orchid Cancer Appeal, The National Cancer Research Network UK, The National Cancer Research Institute (NCRI) UK. We are grateful for support of NIHR funding to the NIHR Biomedical Research Centre at The Institute of Cancer Research and The Royal Marsden NHS Foundation Trust. 

The APBC BioResource, which forms part of the PRACTICAL consortium, consists of the following members: Wayne Tilley, Gail Risbridger, Renea Taylor, Judith A Clements, Lisa Horvath, Vanessa Hayes, Lisa Butler, Trina Yeadon, Allison Eckert, Pamela Saunders, Anne-Maree Haynes, Melissa Papargiris. The US GWAS was supported by a grant from the National Institutes of Health ( P01CA055819 ). The German study was supported by the Dietmar-Hopp-Stiftung, Germany, the German Cancer Aid ( 110, 131 ), the Federal Ministry of Education and Research ( CLIOMMICS 01ZX1309 ), The German Research Council (DFG; Project SI 236/81, SI 236/)-1, ER 155/6-1 and the DFG CRI 216 ), the Harald Huppert Foundation and the Multiple Myeloma Research Foundation. The patients were collected by the GMMG and DSMM studies. The German GWAS made use of genotyping data from the population-based HNR study, which is supported by the Heinz Nixdorf Foundation (Germany). The genotyping of the Illumina HumanOmni-1 Quad BeadChips of the HNR subjects was financed by the German Center for Neurodegenerative Disorders (DZNE), Bonn. We are grateful to all investigators who contributed to the generation of this data set. The German replication controls were collected by Peter Bugert, Institute of Transfusion Medicine and Immunology, Medical Faculty Mannheim, Heidelberg University, German Red Cross Blood Service of Baden-Württemberg-Hessen, Mannheim, Germany. This work was supported by research grants from the Swedish Foundation for Strategic Research ( KF10-0009 ), the Marianne and Marcus Wallenberg Foundation ( 2010. 0112 ), the Knut and Alice Wallenberg Foundation ( 2012. 0193 ), the Swedish Research Council ( 2012-1753 ), the Royal Swedish Academy of Science, ALF grants to the University and Regional Laboratories (Labmedicin Skåne), the Siv-Inger and Per-Erik Andersson Foundation, the Medical Faculty at Lund University, the Borås foundation for cancer research, and the Swedish Society of Medicine. We thank Jörgen Adolfsson, Tomas Axelsson, Anna Collin, Ildikó Frigyesi, Patrik Magnusson, Bertil Johansson, Jan Westin and Helga Ögmundsdóttir for their assistance. This work was supported by Center for Translational Molecular Medicine (BioCHIP), a clinical research grant from the European Hematology Association, an EMCR Translational Research Grant, a BMBF grant from CLIOMMICS ( 01ZX1309A ) and FP7 grant MSCNET (LSHC-Ct-2006-037602 ). N. W. was supported by the National Institute of General Medical Sciences of the National Institutes of Health under Award Number P20GM125503. We thank the staff of the HOVON, as well as patients and physicians at participating sites. In addition, we also thank Jasper Koenders, Michael Vermeulen, André Uitterlinden and Nathalie van der Velde for their assistance. We are indebted to the clinicians who contributed samples to Swedish, Norwegian and Danish biobanks. We are indebted to the patients
section_name: Acknowledgements

section_name: Data availability

section_content: SNP genotyping data that support the findings of this study have been deposited in Gene Expression Omnibus with accession codes GSE21349, GSE19784, GSE24080, GSE2658 and GSE15695; in the European Genomephenome Archive (EGA) with accession code EGAS00000000001; in the European Bioinformatics Institute (Part of the European Molecular Biology Laboratory) (EMBL-EBI) with accession code E-MTAB-362 and E-TABM-1138; and in the database of Genotypes and Phenotypes (dbGaP) with accession code phs000207. v1. p1. The remaining data are contained within the paper and Supplementary Files or available from the author upon request. Naive B-cell HiC data used in this work is publicly available from Blueprint Epigenome Project [https://osf. io/u8tzp/]. ChIP-seq data for H3K27ac, H3K4Me1, H3K27Me3, H3K9Me3, H3K36Me3 and H3K27Me3 from naive B cells are publicly available and was obtained from Blueprint Epigenome Project [http:// www. blueprint-epigenome. eu/].
section_name: Data availability

section_name: Conflict of interest

section_content: The authors declare that they have no conflict of interest. 

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Information accompanies this paper at (https://doi. org/ 10. 1038/s41408-018-0162-8).
section_name: Conflict of interest