A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk

Pleiotropy, 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.

Does a Multiple Myeloma Polygenic Risk Score Predict Overall Survival of Myeloma Patients?

BACKGROUND

Genome-wide association studies (GWAS) of multiple myeloma (MM) 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 MM overall survival (OS) in multiple populations of EA (IMMEnSE consortium, InterLymph consortium, CoMMpass and the German GWAS) for a total of 3748 MM 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 MM risk SNPs were associated with OS (p<0.05), specifically POT1-AS1-rs2170352 (HR=1.37, 95% C.I.=1.09-1.73, p=0.007) and TNFRSF13B-rs4273077 (HR=1.19, 95% C.I.=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 MM risk SNPs and OS.

IMPACT

This is the first study to investigate the association between MM PRS and OS in MM.

Genetically determined telomere length and multiple myeloma risk and outcome

Telomeres are involved in processes like cellular growth, chromosomal stability, and proper segregation to daughter cells. Telomere length measured in leukocytes (LTL) has been investigated in different cancer types, including multiple myeloma (MM). However, LTL measurement is prone to heterogeneity due to sample handling and study design (retrospective vs. prospective). LTL is genetically determined; genome-wide association studies identified 11 SNPs that, combined in a score, can be used as a genetic instrument to measure LTL and evaluate its association with MM risk. This approach has been already successfully attempted in various cancer types but never in MM. We tested the "teloscore" in 2407 MM patients and 1741 controls from the International Multiple Myeloma rESEarch (IMMeNSE) consortium. We observed an increased risk for longer genetically determined telomere length (gdTL) (OR = 1.69; 95% CI 1.36-2.11; P = 2.97 × 10^-6 for highest vs. lowest quintile of the score). Furthermore, in a subset of 1376 MM patients we tested the relationship between the teloscore and MM patients survival, observing a better prognosis for longer gdTL compared with shorter gdTL (HR = 0.93; 95% CI 0.86-0.99; P = 0.049). In conclusion, we report convincing evidence that longer gdTL is a risk marker for MM risk, and that it is potentially involved in increasing MM survival.

Common gene variants within 3'-untranslated regions as modulators of multiple myeloma risk and survival

We evaluated the association between germline genetic variants located within the 3'-untranlsated region (polymorphic 3'UTR, ie, p3UTR) of candidate genes involved in multiple myeloma (MM). We performed a case-control study within the International Multiple Myeloma rESEarch (IMMEnSE) consortium, consisting of 3056 MM patients and 1960 controls recruited from eight countries. We selected p3UTR of six genes known to act in different pathways relevant in MM pathogenesis, namely KRAS (rs12587 and rs7973623), VEGFA (rs10434), SPP1 (rs1126772), IRF4 (rs12211228) and IL10 (rs3024496). We found that IL10-rs3024496 was associated with increased risk of developing MM and with a worse overall survival of MM patients. The variant allele was assayed in a vector expressing eGFP chimerized with the IL10 3'-UTR and it was found functionally active following transfection in human myeloma cells. In this experiment, the A-allele caused a lower expression of the reporter gene and this was also in agreement with the in vivo expression of mRNA measured in whole blood as reported in the GTEx portal. Overall, these data are suggestive of an effect of the IL10-rs3024496 SNP on the regulation of IL10 mRNA expression and it could have clinical implications for better characterization of MM patients in terms of prognosis.

Inherited variation in the xenobiotic transporter pathway and survival of multiple myeloma patients

Over the past four decades, remarkable progress has been made in the treatment and prognosis of multiple myeloma (MM), although it remains an incurable disease. Chemotherapy resistance is a major hurdle for treatment efficacy. Drug resistance can be innate and so driven by genes involved in the drug metabolism pathways. We performed an association study of 71 germline variants within the major genes in those pathways (ABCB1, ABCC2, ABCG2, and their regulators NR1I2/PXR and NR1I3/CAR) in the International Multiple Myeloma rESEarch (IMMEnSE) consortium, consisting of 1365 MM cases with survival information recruited in 5 European countries. Two of the SNPs showed a significant association with the survival of MM patients, namely rs2235013, located in ABCB1 [Hazard ratio (HR) = 1·52, 95% confidence interval (CI) = 1·18-1·95, P = 0·00087], and rs4148388, located in ABCC2 (HR = 2·15, 95% CI = 1·44-3·22, P = 0·0001). ABCC2 plays an essential role in transporting various anticancer drugs, including several used against MM, out of the cell. In silico analyses predict that the variant alleles of four SNPs in linkage disequilibrium with ABCC2-rs4148388 are associated with increased gene expression. Overexpression of ABCC2 increases drug clearance and therefore may induce drug resistance mechanisms. In conclusion, we found a promising association between ABCC2-rs4148388 and MM outcome that is supported by a plausible biological explanation.

Abstract 5078: Genome wide association study identifies variants at 16p13 associated with survival in multiple myeloma patients

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