Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization

Autophagy; Genetic variants; Multiple myelomaAutofagia; Variantes genéticas; Mieloma múltipleAutofàgia; Variants genètiques; Mieloma múltipleMultiple 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.This work was supported by the European Union’s Horizon 2020 research and innovation program, N° 856620 and by grants from the Instituto de Salud Carlos III and FEDER (Madrid, Spain; PI17/02256 and PI20/01845), Consejería de Transformación Económica, Industria, Conocimiento y Universidades and FEDER (PY20/01282), from the CRIS foundation against cancer, from the Cancer Network of Excellence (RD12/10 Red de Cáncer), from the Dietmar Hopp Foundation and the German Ministry of Education and Science (BMBF: CLIOMMICS [01ZX1309]), and from National Cancer Institute of the National Institutes of Health under award numbers: R01CA186646, U01CA249955 (EEB). This work was also funded d by Portuguese National funds, through the Foundation for Science and Technology (FCT)—project UIDB/50026/2020 and UIDP/50026/2020 and by the project NORTE-01-0145-FEDER-000055, supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF)

GWAS-Identified Variants for Obesity Do Not Influence the Risk of Developing Multiple Myeloma: A Population-Based Study and Meta-Analysis

Genetic variants; Multiple myeloma; ObesityVariantes genéticas; Mieloma múltiple; ObesidadVariants genètiques; Mieloma múltiple; ObesitatMultiple myeloma (MM) is an incurable disease characterized by the presence of malignant plasma cells in the bone marrow that secrete specific monoclonal immunoglobulins into the blood. Obesity has been associated with the risk of developing solid and hematological cancers, but its role as a risk factor for MM needs to be further explored. Here, we evaluated whether 32 genome-wide association study (GWAS)-identified variants for obesity were associated with the risk of MM in 4189 German subjects from the German Multiple Myeloma Group (GMMG) cohort (2121 MM cases and 2068 controls) and 1293 Spanish subjects (206 MM cases and 1087 controls). Results were then validated through meta-analysis with data from the UKBiobank (554 MM cases and 402,714 controls) and FinnGen cohorts (914 MM cases and 248,695 controls). Finally, we evaluated the correlation of these single nucleotide polymorphisms (SNPs) with cQTL data, serum inflammatory proteins, steroid hormones, and absolute numbers of blood-derived cell populations (n = 520). The meta-analysis of the four European cohorts showed no effect of obesity-related variants on the risk of developing MM. We only found a very modest association of the POC5rs2112347G and ADCY3rs11676272G alleles with MM risk that did not remain significant after correction for multiple testing (per-allele OR = 1.08, p = 0.0083 and per-allele OR = 1.06, p = 0.046). No correlation between these SNPs and functional data was found, which confirms that obesity-related variants do not influence MM risk.This work was supported by grants from the Instituto de Salud Carlos III (Madrid, Spain; PI17/02256 and PI20/01845), from the Consejería de Salud y Familia de la Junta de Andalucía (PY20/01282) and from the Dietmar Hopp Foundation and the German Ministry of Education and Science (BMBF: CLIOMMICS (01ZX1309))

Identification of Recessively Inherited Genetic Variants Potentially Linked to Pancreatic Cancer Risk

Although 21 pancreatic cancer susceptibility loci have been identified in individuals of European ancestry through genome-wide association studies (GWASs), much of the heritability of pancreatic cancer risk remains unidentified. A recessive genetic model could be a powerful tool for identifying additional risk variants. To discover recessively inherited pancreatic cancer risk loci, we performed a re-analysis of the largest pancreatic cancer GWAS, the Pancreatic Cancer Cohort Consortium (PanScan) and the Pancreatic Cancer Case-Control Consortium (PanC4), including 8,769 cases and 7,055 controls of European ancestry. Six single nucleotide polymorphisms (SNPs) showed associations with pancreatic cancer risk according to a recessive model of inheritance. We replicated these variants in 3,212 cases and 3,470 controls collected from the PANcreatic Disease ReseArch (PANDoRA) consortium. The results of the meta-analyses confirmed that rs4626538 (7q32.2), rs7008921 (8p23.2) and rs147904962 (17q21.31) showed specific recessive effects (p10-3), although none of the six SNPs reached the conventional threshold for genome-wide significance (p < 5×10-8). Additional bioinformatic analysis explored the functional annotations of the SNPs and indicated a possible relationship between rs36018702 and expression of the BCL2L11 and BUB1 genes, which are known to be involved in pancreatic biology. Our findings, while not conclusive, indicate the importance of considering non-additive genetic models when performing GWAS analysis. The SNPs associated with pancreatic cancer in this study could be used for further meta-analysis for recessive association of SNPs and pancreatic cancer risk and might be a useful addiction to improve the performance of polygenic risk scores

A polygenic risk score for multiple myeloma risk prediction

This work was partially supported by intramural funds of the University of Pisa, DKFZ, and University Hospital of Southern Jutland, Denmark, and by a grant of the French National Cancer Institute (INCA). The authors wish to thank Dr. Dominic Edelmann (Division of Biostatistics, DKFZ) for helpful advice about data analysis.There is overwhelming epidemiologic evidence that the risk of multiple myeloma (MM) has a solid genetic background. Genome-wide association studies (GWAS) have identified 23 risk loci that contribute to the genetic susceptibility of MM, but have low individual penetrance. Combining the SNPs in a polygenic risk score (PRS) is a possible approach to improve their usefulness. Using 2361 MM cases and 1415 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium, we computed a weighted and an unweighted PRS. We observed associations with MM risk with OR = 3.44, 95% CI 2.53-4.69, p = 3.55 x 10(-15) for the highest vs. lowest quintile of the weighted score, and OR = 3.18, 95% CI 2.1 = 34-4.33, p = 1.62 x 10(-13) for the highest vs. lowest quintile of the unweighted score. We found a convincing association of a PRS generated with 23 SNPs and risk of MM. Our work provides additional validation of previously discovered MM risk variants and of their combination into a PRS, which is a first step towards the use of genetics for risk stratification in the general population.University of Pisa, DKFZUniversity Hospital of Southern Jutland, DenmarkInstitut National du Cancer (INCA) Franc

Identification of miRSNPs associated with the risk of multiple myeloma

Accepted articleMultiple myeloma (MM) is a malignancy of plasma cells usually infiltrating the bone marrow, associated with the production of a monoclonal immunoglobulin (M protein) which can be detected in the blood and/or urine. Multiple lines of evidence suggest that genetic factors are involved in MM pathogenesis, and several studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to the disease. SNPs within miRNA-binding sites in target genes (miRSNPs) may alter the strength of miRNA-mRNA interactions, thus deregulating protein expression. MiRSNPs are known to be associated with risk of various types of cancer, but they have never been investigated in MM. We performed an in silico genome-wide search for miRSNPs predicted to alter binding of miRNAs to their target sequences. We selected 12 miRSNPs and tested their association with MM risk. Our study population consisted of 1,832 controls and 2,894 MM cases recruited from seven European countries and Israel in the context of the IMMEnSE (International Multiple Myeloma rESEarch) consortium. In this population two SNPs showed an association with p<0.05: rs286595 (located in gene MRLP22) and rs14191881 (located in gene TCF19). Results from IMMEnSE were meta-analyzed with data from a previously published genome-wide association study (GWAS). The SNPs rs13409 (located in the 3UTR of the POU5F1 gene), rs1419881 (TCF19), rs1049633, rs1049623 (both in DDR1) showed significant associations with MM risk. In conclusion, we sought to identify genetic polymorphisms associated with MM risk starting from genome-wide prediction of miRSNPs. For some mirSNPs, we have shown promising associations with MM risk. What's new? Even though deregulation of miRNA expression has been associated with human cancers little information is available regarding their relation with MM susceptibility. We performed an in silico genome-wide search for miRSNPs and selected the most promising ones for an association study. The SNPs with the strongest associations with MM risk are localized in genes which have never been related with MM.This work was partially funded by: intramural funds of German Cancer Research Center (DKFZ), Grant ref. HUS412A1271 from the “Gerencia Regional de Salud de la Junta de Castilla y Léon”. This work was supported by grants from the Instituto de Salud Carlos III (Madrid, Spain; PI12/02688). Catalan Government DURSI grant 2014SGR647 and Instituto de Salud Carlos III, co7funded by FEDER funds –a way to build Europe– grants PI11701439 and PIE13/00022info:eu-repo/semantics/publishedVersio

Genetically determined telomere length and multiple myeloma risk and outcome

This work was partially supported by intramural funds of Univerity of Pisa and DKFZ; by Fondo de Investigaciones Sanitarias (Madrid, Spain) [PI12/02688 to J. S., PI17/02276 to J.S.]; by Instituto de Salud Carlos III, co-funded by FEDER funds —a way to build Europe—[PI14-00613 to V.M.] and by Agency for Management of University and Research Grants (AGAUR) of the Catalan Government (Barcelona, Spain) [2017SGR723 to V.M.]. Open Access funding enabled and organized by Projekt DEAL.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 x 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.Univerity of PisaHelmholtz AssociationInstituto de Salud Carlos III PI12/02688 PI17/02276Instituto de Salud Carlos IIIEuropean CommissionFEDER funds-a way to build Europe PI14-00613Agency for Management of University and Research Grants (AGAUR) of the Catalan Government (Barcelona, Spain) 2017SGR723Projekt DEA

Polymorphisms within autophagy-related genes as susceptibility biomarkers for multiple myeloma: a meta-analysis of three large cohorts and functional characterization

Functional data used in this project have been meticulously catalogued and archived in the BBMRI-NL data infrastructure (https://hfgp.bbmri.nl/, accessed on 12 February 2020) using the MOLGENIS open-source platform for scientific data.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 hemoattractant Protein (MCP)-2 (p = 3.6 × 10−4). We also found that the CD46rs1142469 SNP corre lated 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.This work was supported by the European Union’s Horizon 2020 research and innovation program, N° 856620 and by grants from the Instituto de Salud Carlos III and FEDER (Madrid, Spain; PI17/02256 and PI20/01845), Consejería de Transformación Económica, Industria, Conocimiento y Universidades and FEDER (PY20/01282), from the CRIS foundation against cancer, from the Cancer Network of Excellence (RD12/10 Red de Cáncer), from the Dietmar Hopp Foundation and the German Ministry of Education and Science (BMBF: CLIOMMICS [01ZX1309]), and from National Cancer Institute of the National Institutes of Health under award numbers: R01CA186646, U01CA249955 (EEB).This work was also funded d by Portuguese National funds, through the Foundation for Science and Technology (FCT)—project UIDB/50026/2020 and UIDP/50026/2020 and by the project NORTE-01-0145-FEDER-000055, supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF)

Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization

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 \u3c 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

Identification of germline variants in risk and survival of multiple myeloma

Multiple myeloma (MM) is a plasma cell malignancy, originating from the bone marrow. Over the past decades, remarkable progress has been made in the understanding of the biology and pathogenesis of MM which, in turn, led to significant improvements in the disease treatment and patients survival, but most of the patients eventually relapse, making MM still an incurable disease. A small number of environmental and lifestyle-related risk factors for MM have been identified. Familial aggregation of MM and its precursor monoclonal gammopathy of undetermined significance (MGUS) suggest that genetic factors play a role in risk of MM as well. Several genetic loci affecting MM risk have been identified so far with both a candidate gene approach and genome-wide analysis, even if they are still considered few comparing with other better studied and more common cancer types. The genetics behind the differences in prognosis of MM patients is still poorly known and only two genome-wide association studies (GWAS) have been attempted so far. The main goal of this project was to discover new variants and new key genes that affect risk and outcome of MM. To achieve it, we performed several association studies in a case-control population, carried out with a candidate gene approach. The association studies were performed in the context of the IMMEnSE (International Multiple Myeloma rESEarch) consortium and for some studies a replication population was used. In particular, we used data from the InterLymph consortium, the German myeloma group and controls from the ESTHER consortium, for a total of over 3000 cases and 2000 controls. In particular, we investigated: SNPs in genes involved in the DNA repair mechanisms and their relation with both risk and survival of MM patients, expression quantitative traits loci (eQTLs) of genes whose expression is known to affect MM prognosis, SNPs in genes involved in the xenobiotic metabolisms, SNPs known to affect MM risk in association with MM survival. Additionally, we considered genetic variants not only individually but also combined in scores, and explored their performance in predicting risk and outcome of MM. Regarding MM risk we found a new association with SNPs within genes involved in DNA repair system. We also successfully replicated the 23 risk loci emerged from GWASs that were also significantly associated with MM risk when combined in a polygenic score. In particular carriers of more than 20 risk alleles showed an increased risk of MM of almost 3 fold (OR=2.99, 95% C.I.=2.14-4.18, p=1.58×10-10). Regarding MM survival with our approaches we found 14 new SNP associations, which were obtained using overall survival as endpoint. Combining those SNPs in a score showed promising results, indeed, we found that carriers of more than 15 “survival alleles” had shorter survival compared with carriers of less than 10 alleles (HR=1.80, 95% C.I=1.27-2.54, p=0.001). In conclusion, our results contribute to expand the knowledge of the genetic architecture of MM. This will lead to the development of useful tools to facilitate early diagnosis by screening the general population based on their genetic risk, and to stratify patients based on their response to therapy and outcome

Role of mirSNP in Multiple Myeloma risk

Multiple myeloma (MM) is a malignancy of plasma cells usually infiltrating the bone marrow, associated with the production of a monoclonal immunoglobulin (M protein) which can be detected in the blood and/or urine. About 63 000 subjects are reported to die from the disease each year, accounting for 0.9% of all cancer deaths and nearly 10% of all haematological neoplastic diseases. Prognosis is usually unfavourable, with a mean overall survival ranging between 20 and 60 months. It has been shown that MM usually evolves from an asymptomatic premalignant condition termed monoclonal gammopathy of undetermined significance (MGUS). In some patients, an intermediate asymptomatic, but more advanced premalignant stage, defined as smouldering multiple myeloma (SMM) could be clinically recognized. Converging evidence of MM in monozygotic twins and familial aggregation of MM strongly suggest that MM aetiology has a robust genetic component. Single Nucleotide Polymorphisms (SNPs) are the major source of genetic variation in humans and thought to be responsible, at least in part, for the individual differences in genetic susceptibility to complex diseases as tumors. In the last ten years several studies have identified SNPs associated with the disease susceptibility. Genetic polymorphisms in miRNA-binding sites in target genes may alter the strength of miRNA– mRNA interactions, thus deregulating protein levels. This category of SNPs are called miR-SNPs. I analyzed 12 SNPs located in the 3' UTR region of miRNAs target genes, with the aim of testing whether they are associated with MM risk. Those MiR-SNP were previously selected from an in silico genome-wide search for their potential ability to alter binding of miRNAs to their target sequences. My study population consisted of 1935 controls and 2457 cases recruited from 7 European countries and from Israel and Japan in the context of the IMMEnSE (International Multiple Myeloma rESEarch) consortium. I performed the genotyping with TaqMan technology. Association between SNPs and multiple myeloma risk was assessed with unconditional logistic regression using allelic, codominant, dominant and recessive inheritance models, adjusting by age, gender, and region of origin. Afterwords I performed a meta-analysis between my data and the data from a previously published genome-wide association study (GWAS). The SNPs rs13409 (located in the 3’UTR of the PUOF5 gene), rs1419881 (CCHCR1), rs1049633, rs1049623 (both in DDR1) have shown significant associations with MM risk, with no heterogeneity between IMMEnSE and the GWAS in most of inheritance models tested