Sáinz, Juan

Summary Evidence for an association between insulin‐like growth factors (IGF) and multiple myeloma (MM) is inconsistent. We examined total IGF‐I concentrations and risk of MM by combining baseline serological data among UK Biobank participants ( n = 444 187; 732 incident MM) with a two‐sample Mendelian randomisation (MR) analysis using identified genetic variants associated with circulating total IGF‐I and IGF‐binding protein 3 (IGFBP‐3) in the InterLymph consortium (2434 MM and their 2567 controls). Finally, additional lymphoid neoplasm (LN) subtypes were included for comparison with the main hypothesis. Circulating IGF‐I level was positively associated with MM risk Hazard ratio‐HR‐per one standard deviation‐SD‐increase (HR 1‐SD = 1.11, 95% confidence interval [CI]: 1.01–1.22; p ‐value = 0.03), especially closer to diagnosis. Genetically inferred IGF‐I levels were associated with increased MM risk (odds ratio [OR] = 1.27, 95% CI: 1.05–1.54) but not with any other LNs. Genetically inferred IGFBP‐3 levels showed no associations with any LN evaluated. Corroborating previous findings, in a secondary analysis, IGF‐I levels were associated with the risk of chronic lymphocytic leukaemia/small lymphocytic lymphoma (CLL/SLL) in males with higher body mass index (HR 1‐SD in obese male = 1.36, 95% CI: 1.14–1.61). Our serological and MR analyses suggest a contributing role of IGF‐I in the susceptibility of MM; the CLL/SLL findings warrant further investigation considering sex‐specific adiposity.

Abstract Monoclonal gammopathy of undetermined significance (MGUS) is a necessary precursor condition to multiple myeloma (MM). Given the role of autophagy in modulating MM risk, we investigated whether genetic variation in autophagy‐related genes influences susceptibility to MGUS. We analyzed the association of 34,042 common autophagy‐related single nucleotide polymorphisms (SNPs) with MGUS across six independent cohorts, five from Europe and one from North America, comprising 2317 MGUS cases and 282,358 controls. We also assessed their impact on immune parameters, including absolute counts of 91 blood‐derived immune cell subsets and 103 circulating immunological proteins. Meta‐analysis revealed a genome‐wide significant association between the ULK4 rs6599175C allele and increased MGUS risk ( p = 3.35 × 10 −8 ). Carriers of this allele showed reduced counts of memory B cell subsets (IgM+CD38+CD27+ and IgD+IgM+CD27+; p = .0038 and p = .0056, respectively) and natural effector B cells (CD24+CD38+IgD+IgM+ cells; p = .0060). Although these associations were not statistically significant after multiple testing correction, they suggest a role of ULK4 in early B‐cell differentiation. Additionally, the CDKN2A rs2811710 variant showed a suggestive association with MGUS risk ( p = 2.17 × 10 −4 ), affecting transcription factor binding involved in B‐cell proliferation and differentiation, although it lacked association with immune markers. In conclusion, we confirm a genome‐wide significant association of the ULK4 locus and MGUS risk, supporting its role in early B‐cell differentiation, and identify CDKN2A as a candidate susceptibility locus warranting further investigation.

Abstract We investigated the influence of 55 583 autophagy-related single-nucleotide polymorphisms (SNPs) on chronic lymphocytic leukemia (CLL) risk across 4 independent populations comprising 5472 CLL cases and 726 465 controls. We also examined their impact on overall survival (OS), time to first treatment (TTFT), autophagy flux, and immune responses. A meta-analysis of the 4 populations identified, to our knowledge, for the first time, significant associations between CDKN2A (rs3731204) and BCL2 (rs4940571, rs12457371, and rs1026825) SNPs and CLL risk, with CDKN2A showing the strongest association (P = 1.57 × 10−12). We also validated previously reported associations for FAS, BCL2, and BAK1 SNPs with CLL risk (P = 4.73 × 10−21 to 3.39 × 10−9). The CDKN2Ars3731204 and FASrs1926194 SNPs associated with increased CDKN2A and ACTA2 messenger RNA expression levels in the whole blood and/or lymphocytes (P = 5.1 × 10−7, P = 1.58 × 10−21, and P = 7.8 × 10−41), although no significant effect on autophagy flux was observed. However, associations were found between CDKN2A, BCL2, and FAS SNPs and various T-cell subsets, cytokine production, and circulating concentrations of interferon gamma, tumor necrosis factor–related apoptosis-inducing ligand, CD40, chemokine ligand 20, and interleukin-2 receptor subunit β proteins (P ≤ .005). No significant association was detected between autophagy variants and OS or TTFT, suggesting that these variants drive disease initiation rather than progression. In conclusion, this study identified 4 novel associations for CLL and provided insights into the biological pathways that influence CLL development.

Multiple myeloma (MM) arises following malignant proliferation of plasma cells in the bone marrow, that secrete high amounts of specific monoclonal immunoglobulins or light chains, resulting in the massive production of unfolded or misfolded proteins. Autophagy can have a dual role in tumorigenesis, by eliminating these abnormal proteins to avoid cancer development, but also ensuring MM cell survival and promoting resistance to treatments. To date no studies have determined the impact of genetic variation in autophagy-related genes on MM risk. We performed meta-analysis of germline genetic data on 234 autophagy-related genes from three independent study populations including 13,387 subjects of European ancestry (6863 MM patients and 6524 controls) and examined correlations of statistically significant single nucleotide polymorphisms (SNPs; p < 1 × 10−9) with immune responses in whole blood, peripheral blood mononuclear cells (PBMCs), and monocyte-derived macrophages (MDM) from a large population of healthy donors from the Human Functional Genomic Project (HFGP). We identified SNPs in six loci, CD46, IKBKE, PARK2, ULK4, ATG5, and CDKN2A associated with MM risk (p = 4.47 × 10−4−5.79 × 10−14). Mechanistically, we found that the ULK4rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10−4), whereas the IKBKErs17433804 SNP correlated with the number of transitional CD24+CD38+ B cells (p = 4.8 × 10−4) and circulating serum concentrations of Monocyte Chemoattractant Protein (MCP)-2 (p = 3.6 × 10−4). We also found that the CD46rs1142469 SNP correlated with numbers of CD19+ B cells, CD19+CD3− B cells, CD5+IgD− cells, IgM− cells, IgD−IgM− cells, and CD4−CD8− PBMCs (p = 4.9 × 10−4−8.6 × 10−4) and circulating concentrations of interleukin (IL)-20 (p = 0.00082). Finally, we observed that the CDKN2Ars2811710 SNP correlated with levels of CD4+EMCD45RO+CD27− cells (p = 9.3 × 10−4). These results suggest that genetic variants within these six loci influence MM risk through the modulation of specific subsets of immune cells, as well as vitamin D3−, MCP-2−, and IL20-dependent pathways.

Chronic lymphocytic leukemia (CLL) is the most common leukemia among adults worldwide. Although genome-wide association studies (GWAS) have uncovered the germline genetic component underlying CLL susceptibility, the potential use of GWAS-identified risk variants to predict disease progression and patient survival remains unexplored. Here, we evaluated whether 41 GWAS-identified risk variants for CLL could influence overall survival (OS) and disease progression, defined as time to first treatment (TTFT) in a cohort of 1039 CLL cases ascertained through the CRuCIAL consortium. Although this is the largest study assessing the effect of GWAS-identified susceptibility variants for CLL on OS, we only found a weak association of ten single nucleotide polymorphisms (SNPs) with OS (p < 0.05) that did not remain significant after correction for multiple testing. In line with these results, polygenic risk scores (PRSs) built with these SNPs in the CRuCIAL cohort showed a modest association with OS and a low capacity to predict patient survival, with an area under the receiver operating characteristic curve (AUROC) of 0.57. Similarly, seven SNPs were associated with TTFT (p < 0.05); however, these did not reach the multiple testing significance threshold, and the meta-analysis with previous published data did not confirm any of the associations. As expected, PRSs built with these SNPs showed reduced accuracy in prediction of disease progression (AUROC = 0.62). These results suggest that susceptibility variants for CLL do not impact overall survival and disease progression in CLL patients.

AbstractPleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome‐wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10−8) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4‐rs34517439‐A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13‐1.32, P = 4.81 × 10−7). rs34517439‐A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA‐binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4‐rs34517439 as a potentially novel MM risk locus.

Abstract Background: Genome-wide association studies (GWAS) of multiple myeloma in populations of European ancestry (EA) identified and confirmed 24 susceptibility loci. For other cancers (e.g., colorectum and melanoma), risk loci have also been associated with patient survival. Methods: We explored the possible association of all the known risk variants and their polygenic risk score (PRS) with multiple myeloma overall survival (OS) in multiple populations of EA [the International Multiple Myeloma rESEarch (IMMEnSE) consortium, the International Lymphoma Epidemiology consortium, CoMMpass, and the German GWAS] for a total of 3,748 multiple myeloma cases. Cox proportional hazards regression was used to assess the association between each risk SNP with OS under the allelic and codominant models of inheritance. All analyses were adjusted for age, sex, country of origin (for IMMEnSE) or principal components (for the others) and disease stage (ISS). SNP associations were meta-analyzed. Results: SNP associations were meta-analyzed. From the meta-analysis, two multiple myeloma risk SNPs were associated with OS (P &lt; 0.05), specifically POT1-AS1-rs2170352 [HR = 1.37; 95% confidence interval (CI) = 1.09–1.73; P = 0.007] and TNFRSF13B-rs4273077 (HR = 1.19; 95% CI = 1.01–1.41; P = 0.04). The association between the combined 24 SNP MM-PRS and OS, however, was not significant. Conclusions: Overall, our results did not support an association between the majority of multiple myeloma risk SNPs and OS. Impact: This is the first study to investigate the association between multiple myeloma PRS and OS in multiple myeloma.

AbstractGene expression profiling can be used for predicting survival in multiple myeloma (MM) and identifying patients who will benefit from particular types of therapy. Some germline single nucleotide polymorphisms (SNPs) act as expression quantitative trait loci (eQTLs) showing strong associations with gene expression levels. We performed an association study to test whether eQTLs of genes reported to be associated with prognosis of MM patients are directly associated with measures of adverse outcome. Using the genotype‐tissue expression portal, we identified a total of 16 candidate genes with at least one eQTL SNP associated with their expression with P < 10−7 either in EBV‐transformed B‐lymphocytes or whole blood. We genotyped the resulting 22 SNPs in 1327 MM cases from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and examined their association with overall survival (OS) and progression‐free survival (PFS), adjusting for age, sex, country of origin and disease stage. Three polymorphisms in two genes (TBRG4‐rs1992292, TBRG4‐rs2287535 and ENTPD1‐rs2153913) showed associations with OS at P < .05, with the former two also associated with PFS. The associations of two polymorphisms in TBRG4 with OS were replicated in 1277 MM cases from the International Lymphoma Epidemiology (InterLymph) Consortium. A meta‐analysis of the data from IMMEnSE and InterLymph (2579 cases) showed that TBRG4‐rs1992292 is associated with OS (hazard ratio = 1.14, 95% confidence interval 1.04‐1.26, P = .007). In conclusion, we found biologically a plausible association between a SNP in TBRG4 and OS of MM patients.

Abstract Background: Genome wide association studies (GWAS) have recently identified 8 loci identified with multiple myeloma (MM) susceptibility. We hypothesized that germ line genetic variants may also affect MM progression. We performed the first GWAS of MM survival, by conducting a meta-analysis of two existing studies at University of California San Francisco (UCSF) and the Mayo Clinic, and replicating findings in the International Multiple Myeloma rESEarch (IMMEnSE) consortium. Methods: The UCSF study included 353 MM patients genotyped on an Illumina Omni5 and Illumina 660 array. The Mayo clinic study included 239 MM patients genotyped on an Affymetrix 6.0 array. We imputed missing genotypes using the 1000 Genomes dataset and performed a GWAS for survival using proportional hazards models adjusting for age and gender and genetic ancestry using principal components analysis. We conducted a meta-analysis of the results from the two GWAS. We replicated the top SNPs in the IMMEnSE cohort which included 772 patients with survival data from 7 European countries and one North American site. We used the dataset from Grundberg et al (Nat Genetics 2012) to test for an association between SNPs associated with survival and gene expression in immortalized lymphocytes. We used gene expression data from Zhan et al (Blood 2006) to analyze the association between expression and survival. Results: In meta-analysis, we found a genome wide significant association between SNPs on a region at 16p13 and survival (rs72773978; p=2.9x10e-10). Patients with the minor allele were at increased risk for mortality (HR 2.73; 95% CI: 1.99 - 3.73). The top associated SNPs were all in strong linkage disequilibrium and were centered over the gene FOPNL. We replicated the association in the IMMEnSE cohort and found a significant association between the top SNP from the discovery dataset (rs72773978 p=0.037). We used previously published data to determine the association between gene expression and the top associated SNPs from the GWAS. The top SNP (rs72773978) was not in previously published eQTL databases, so we used rs7201759 which was also associated with survival in our data (p=3.8x10e-10) as a proxy. The minor allele for this SNP was associated with increased levels of expression of FOPNL (p=9x10e-5), but not with expression of any other genes within a 1 megabase region. Finally, we also observed an association between increased expression of FOPNL and shorter survival (HR: 2.62 per standard deviation.; 95% CI: 1.08 - 6.32). Conclusion: Germ line variants at 16p13 near the FOPNL gene are associated with survival among MM patients. FOPNL is thought to be involved in centrosomal function, and centrosomal activity has previously been associated with survival in MM patients. Our results support this previous observation. Furthermore, they suggest that germ line variants could also be useful in assessing prognosis for patients with MM. Citation Format: Elad Ziv, Eric Dean, Donglei Hu, Alessandro Martino, Daniel Serie, Daniele Campa, Blake Aftab, Paige Bracci, Gabriele Buda, Jennifer Caswell, Charles Dumontet, Marek Dudziński, Laura Fejerman, Alexandra Greenberg, Scott Hunstman, Artur Jurczyszyn, Krzysztof Jamroziak, Shaji Kumar, Herlander Marques, Thomas Martin, Joaquin Martinez-Lopez, Vincent Rajkumar, Juan Sainz, Annette Juul Vangsted, Marzena Watek, Jeffrey Wolf, Susan Slager, Federico Canzian, Celine Vachon. Genome wide association study identifies variants at 16p13 associated with survival in multiple myeloma patients. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5078. doi:10.1158/1538-7445.AM2014-5078