Multiple myeloma (MM) is the second most common hematologic malignancy. The disease is defined by an uninhibited, clonal growth of plasma cells in the bone marrow1. It is preceded by monoclonal gammopathy of unknown significance (MGUS)2, a common condition defined as a clonal growth of plasma cells that does not yet satisfy the criteria for MM, but progresses to MM at a rate of ~1% per year3.
Increasing evidence supports that the biology of MM is influenced by inborn genetic variation. MM and MGUS show familial clustering, and genome-wide association studies have identified DNA sequence variants that influence MM risk4,5,6,7,8. Additionally, two recent studies indicate that genetic variation could also influence MM survival9,10.
In the first of these, Johnson et al.9 describe an association between overall survival in multiple myeloma (MM-OS) and rs12374648, located between the MTHFD1L and AKAP12 genes at chromosome 6q25.19. The protein encoded by MTHFD1L is involved in folate metabolism11, and AKAP12 is related to cell growth12. The association with MM-OS was detected in a meta-analysis of 3256 cases from four clinical trials: two from the UK, one from the USA, and one from Germany (combined P-value = 4.69 × 10–9, hazard ratio, HR = 1.34, 95% CI 1.22–1.48). The association was present in the sample sets from the UK and USA (P = 1.69 × 10–6 to 0.009; HR 1.06 to 1.75), but not in the one from Germany (P = 0.55; HR = 1.09, 95% CI 0.82–1.44). No replication in independent material was carried out after the discovery meta-analysis.
In a second study, Ziv et al.10 describe an association between MM-OS and rs72773978 near FOPNL at 16p13. The protein encoded by FOPNL has been implicated in centrosome function13. The association was detected by meta-analysis of 545 cases from two clinical trials in USA (P = 6 × 10–10; HR = 2.65, 95% CI 1.94–3.58). The association was also seen in a meta-analysis of seven other data sets (six from the IMMEnSE study totaling n = 772 and one from Utah, n = 315) (combined P = 0.044; HR = 1.34, 95% CI 1.01–1.78). Yet, the positive replication result was driven by a P-value of 0.004 with unrealistically large effect size (HR = 9.73) in a subset of 109 patients from Spain in IMMEnSE, whereas the other six subsets (Italy, Poland, Portugal, Denmark, Edmonton in IMMEnSE and the one from Utah) did not show any evidence of association (Supplementary Table 7 in ref. 10).
Given that the MTHFD1L association was not replicated after discovery analysis and that the FOPLN association was based on small sample sizes, it remains a possibility that these two associations are false discoveries due to a winners curse effect. We therefore looked for further support of the MTHFDL1-AKAP12 and FOPLN loci in a Swedish study population. For this, we retrieved clinical data for 871 patients diagnosed with MM between 2005 and 2015 from the Swedish Multiple Myeloma Registry (Sahlgrenska Hospital, Gothenburg) (Table 1), which records clinical data on MM patients in Sweden and has about 90% inclusion rate compared to the Swedish Cancer Registry. The patients had been previously genotyped in genome-wide association studies using population-based samples from the Swedish National Myeloma Biobank (Skåne University Hospital, Lund)6,7. The clinical data and samples were obtained subject to informed consent and ethical approval (Lund University, dnr 2013/540), and in accordance with the principles of the Declaration of Helsinki. The samples were genotyped using Illumina microarrays and imputed with phased reference haplotypes from 1000 Genomes6,14. To test for association between genotypes and MM-OS, we used log rank test implemented in R (v.2.8) with adjustment for age, sex, and International Staging System (ISS) score. Survival was calculated from the date treatment started until the date of death, or until 5 April 2016 (median follow-up time 39.5 months).
In our analysis, we did not see any evidence of association with MM-OS for either rs12374648 (P = 0.84; HR = 0.97, 95% CI = 0.81–1.2) or rs72773978 (P = 0.93; HR = 0.98, 95% CI = 0.7–1.4) (Fig. 1). For completeness, we also tested for associations between MM-OS and all variants with minor allele frequency (MAF) >5% located within 1 Mb of MTHFD1L-AKAP12 (6,515 variants) or FOPNL (3,892 variants) but could not identify any significant association with any of these variants. Thus, we could not replicate the associations between MM-OS and MTHFD1L-AKAP12 and FOPNL in a population-based series, nor identify any other alleles associations with MM-OS at these loci.
Our results, in conjunction with the small sample sizes, lack of robust replication in the original studies, and the fact that the original studies do not replicate each other, indicate that the reported associations are false positives. As for alternative explanations, a first possibility could be limited power of our data set. Yet, our sample is comparable in size (n = 871) to the largest of the reported individual sample sets, including UK-My9 (n = 1,163) and UK-My11 (n = 871) where rs12374648 at MTHFD1L was detected, and substantially larger than the data sets where rs72773978 at FOPNL was detected. Power calculations15,16 indicate that our sample set has about 71% chance to detect an effect with HR = 1.34 (the effect size of rs12374648 and the replication effect size of r72773978), and about 99% chance to detect an effect with HR = 2.65 (the discovery effect size of rs72773978), in our sample set. A second possibility could be differences in geographic origin. However, this also seems unlikely given that the two reported variants are common, both in our data (MAF 21.5 and 4.7%) and in the different populations of 1000 Genomes14. Finally, a third possibility could be differences in clinical characteristics between the study populations. One difference is that our material is population-based, whereas the studies by Johnsson et al.9 and Ziv et al.10 are based on patients recruited into clinical trials. As a result, our population is older (average 68 years vs 54–66 years), and has not been selected for patients without comorbidity, as is common in clinical trials. A higher incidence of comorbidity could dilute effects of DNA sequence variation on survival, and differences in age and comorbidity will carry differences in treatment. For example, some of the reported populations contain a high proportion of patients who received autologous stem cell transplantation (ASCT; 100% in the German and US sample sets in Johnson et al.9), whereas our study population contains 32.5% transplanted patients.
In summary, our results together with the limitations of the original studies indicate that the reported associations between the MTHFD1L and FOPNL loci and MM survival are false positives due to a winner’s curse effect. While there could be alternative explanations, these seem unlikely in comparison. Our results motivate the collection of larger data sets to understand the impact of genetic variation on clinical outcome in MM.
References
Rajkumar, S. V. et al. International Myeloma Working Group updated criteria for the diagnosis of multiple myeloma. Lancet Oncol. 15, e538–e548 (2014).
Weiss, B. M., Abadie, J., Verma, P., Howard, R. S. & Kuehl, W. M. A monoclonal gammopathy precedes multiple myeloma in most patients. Blood 113, 5418–5422 (2009).
Kyle, R. A. et al. Prevalence of monoclonal gammopathy of undetermined significance. N. Engl. J. Med. 354, 1362–1369 (2006).
Broderick, P. et al. Common variation at 3p22.1 and 7p15.3 influences multiple myeloma risk. Nat. Genet 44, 58–61 (2012).
Chubb, D. et al. Common variation at 3q26.2, 6p21.33, 17p11.2 and 22q13.1 influences multiple myeloma risk. Nat. Genet 45, 1221–1225 (2013).
Swaminathan, B. et al. Variants in ELL2 influencing immunoglobulin levels associate with multiple myeloma. Nat. Commun. 6, 7213 (2015).
Mitchell, J. S. et al. Genome-wide association study identifies multiple susceptibility loci for multiple myeloma. Nat. Commun. 7, 12050 (2016).
Weinhold, N. et al. The CCND1 c.870G>A polymorphism is a risk factor for t(11;14)(q13;q32) multiple myeloma. Nat. Genet 45, 522–525 (2013).
Johnson, D. C. et al. Genome-wide association study identifies variation at 6q25.1 associated with survival in multiple myeloma. Nat. Commun. 7, 10290 (2016).
Ziv, E. et al. Genome-wide association study identifies variants at 16p13 associated with survival in multiple myeloma patients. Nat. Commun. 6, 7539 (2015).
Christensen, K. E. & Mackenzie, R. E. Mitochondrial methylenetetrahydrofolate dehydrogenase, methenyltetrahydrofolate cyclohydrolase, and formyltetrahydrofolate synthetases. Vitam Horm. 79, 393–410 (2008).
Su, B., Bu, Y., Engelberg, D. & Gelman, I. H. SSeCKS/Gravin/AKAP12 inhibits cancer cell invasiveness and chemotaxis by suppressing a protein kinase C-Raf/MEK/ERK pathway. J. Biol. Chem. 285, 4578–4586 (2010).
Sedjai, F. et al. Control of ciliogenesis by FOR20, a novel centrosome and pericentriolar satellite protein. J. Cell Sci. 123, 2391–2401 (2010).
The Genomes Project, C. A global reference for human genetic variation. Nature 526, 68–74 (2015).
Latouche, A., Porcher, R. & Chevret, S. Sample size formula for proportional hazards modelling of competing risks. Stat. Med. 23, 3263–3274 (2004).
Schoenfeld, D. A. Sample-size formula for the proportional-hazards regression model. Biometrics 39, 499–503 (1983).
Acknowledgements
We thank Anna Collin and Cecilie Blimark for their assistance, and the patients who participated in the study.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ali, M., Lemonakis, K., Wihlborg, AK. et al. Sequence variation at the MTHFD1L-AKAP12 and FOPNL loci does not influence multiple myeloma survival in Sweden. Blood Cancer J. 9, 57 (2019). https://doi.org/10.1038/s41408-019-0222-8
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41408-019-0222-8
This article is cited by
-
Functional dissection of inherited non-coding variation influencing multiple myeloma risk
Nature Communications (2022)