Citations

Abstract Background: Mantle cell lymphoma (MCL) is an aggressive and rare B-cell non-Hodgkin lymphoma with poor prognosis, even with treatment. Despite its clinical importance, the etiology of MCL remains largely unknown. Some studies have reported a positive association with family history, suggesting that genetic factors could contribute to risk. Methods: To identify germline genetic variants associated with risk, we conducted the first genome-wide association studies (GWAS) and GWAS meta-analysis for MCL, comprising 1,163 cases and 61,271 controls. The most significant loci were taken forward for replication in three independent studies, including 576 cases and 771,773 controls. To identify potential target genes and pathways of discovered loci, we conducted F-MAGMA, colocalization, and transcriptome-wide association study (TWAS) analyses. To gain insight into underlying regulatory mechanisms, we performed integrative epigenomic analyses using DNase I hotspots and histone mark ChIP-seq data across diverse cell types applying the FORGE2 framework. Results: In the joint analyses, we identified eight genome-wide significant loci (P<5x10-8) associated with MCL risk, several of which were located near genes with known DNA repair (ATM), telomere (TERT), or RNA-binding functions (RBM20). Further analyses using F-MAGMA, which uses DNase-seq data to prioritize genes from GWAS, identified genes that were significantly associated with risk, including SP140, a chromatin reader and immune regulator of pathogen response. This gene also displayed strong colocalization with our lead variant in whole blood (COLOC posterior probability PP4=0.97). TWAS demonstrated significant associations with SP140 and ACTA2, which has been linked to cell proliferation (P<3x10-6). Enrichment for DNase I hotspots was observed in B and T lymphocytes (q-value<0.05), providing evidence for a role in immune cell-specific regulatory regions. Enrichment for enhancer-associated histone mark H3K4me1 was observed in B cells (q-value<0.01), highlighting a role for immune cell-specific enhancers. Transcription factor (TF) motif enrichment analysis pointed to TFs critical for B-cell development and differentiation, notably TCF3 and EP300 (q<0.01 for both), suggesting an association between B-cell lineage regulation and MCL risk. Finally, genome-wide SNP-based heritability was estimated to be 17%, underscoring the importance of common variants in MCL risk. Conclusion: Our study provides novel insight into the inherited susceptibility of MCL by identifying eight significant loci and estimating substantial common variant heritability. These findings highlight a role for several distinct biological pathways -specifically DNA repair, telomerase function, and B-cell differentiation- in the etiology of MCL, offering new targets for mechanistic investigation. Citation Format: Charles E. Breeze, Alyssa Clay-Gilmour, Hanla A. Park, Ólafur B. Davíðsson, Angelica Macauda, Brenda M. Birmann, Elizabeth E. Brown, Murat Güler, Michelle A. T. Hildebrandt, Alexandra Nieters, Krystle Ong, Jojo Biel-Nielsen Dietz, Karin E. Smedby, Karl Smith-Byrne, Rosalie Griffin, Sophia S. Wang, Susan Slager, James R. Cerhan, Jonathan N. Hofmann, Qing Lan, Nathaniel Rothman, Ingrid Glimelius, Henrik Hjalgrim, James Mckay, Sonja I. Berndt. Genome-wide association study of mantle cell lymphoma identifies novel loci suggesting a critical role for B-cell chromatin readers and DNA repair mechanisms [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 3602.

Abstract Background. Genome-wide association studies (GWAS) have identified 35 genetic susceptibility single-nucleotide polymorphisms (SNPs) for multiple myeloma (MM) in individuals of European ancestry (EA) and shown strong genetic correlation between MM and its precursor, monoclonal gammopathy of undetermined significance (MGUS). We evaluate the contribution of the 35 MM variants to MGUS susceptibility overall and by prognostic subgroups. Methods. The study included 20,756 participants (14,486 controls, 1,883 with MGUS, and 2,163 with MM) from the Mayo Clinic. Logistic regression assuming an additive model estimated odds ratios (ORs) and 95% confidence intervals (CIs) for individual SNPs and for the MM-PRS, adjusted for age, sex, study, and principal components. The PRS was a weighted sum of 35 SNPs with effect estimates from the largest MM GWAS and was modeled continuously (per SD) and by quintiles (Q1-Q5). SNPnexus annotated variants that replicated in MGUS (P < 0.05 and OR > 1.01) versus those that did not (P > 0.05 and OR < 1.01). Results. The 35-SNP MM-PRS was strongly associated with MM risk and modestly with MGUS. Compared with the middle quintile (Q3), MM odds rose from OR=0.55 (CI=0.46-0.65, P=2.2e-11) in Q1 to OR = 1.82 (CI=1.59-2.09, P=6.2e-18) in Q5. MGUS showed a similar, attenuated pattern (Q1 OR = 0.74, CI=0.64-0.84, P=6.4e-6 / Q5 OR=1.36, CI=1.21-1.53, P=3.4e-7). Each SD increase in PRS corresponded to OR=1.52 (CI=1.45-1.59, P=1.3e-66) for MM and OR = 1.22 (CI=1.18-1.27, P=1.0e-23) for MGUS. Higher PRS values were linked to larger M-protein ((0.1-1.5 g/dL: OR=1.40, CI=1.18-1.65, P=5.8e-12) vs <0.1 g/dL: OR = 1.19, CI=1.13-1.24, P=5.8e-12)) and abnormal free light chain (FLC) ratio ((OR = 1.33, CI=1.21-1.46, P=1.0e-9) vs normal ratio (OR=1.19, CI=1.14-1.25, P=1.4e-13)). Ten risk loci replicated in MGUS, mapping to genes involved in plasma-cell function, immune regulation, and DNA repair and enriched for Rho GTPase signaling, NF-κB-mediated apoptosis, and RNA polymerase II transcription, implicating early plasma-cell activation and transcriptional control. The 12 non-replicating loci, including PHC3, ATG5, and NFIC, mapped to genes involved in chromatin remodeling, autophagy, and SUMOylation, suggesting roles in stress response and genomic maintenance. Conclusions. The MM-PRS captures shared heritability between MM and MGUS and correlates with MGUS subtype and severity. Replicating variants highlight immune and cell-cycle pathways relevant to MGUS onset, whereas non-replicating loci cluster in DNA-repair and stress-response processes, underscoring their potential role in progression. Citation Format: Alyssa Ione Clay-Gilmour, Angelica Macauda, Cristine Allmer, Danelle Moonen, Aaron D. Norman, Nicholas Boddicker, Janet E. Olson, Elizabeth E. Brown, Vincent S. Rajkumar, Esteban Braggio, David Murray, Susan Slager, Shaji Kunnathu Kumar, Celine Vachon. Implications of multiple myeloma polygenic risk scores (PRS) for MGUS [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 3603.

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 Monoclonal gammopathy of undetermined significance (MGUS) is a precursor to multiple myeloma (MM). Clonal hematopoiesis of indeterminate potential (CHIP) is a precursor to either myeloid or lymphoid malignancy, depending on whether the mutation is in a gene associated with lymphoid (L-CHIP) or myeloid (M-CHIP) outcome. Little is known about the association between CHIP and MGUS and their combined effect on risk of MM. We examined this question within Mayo Clinic Biobank cohort of participants aged 50 years or older with no prior history of hematological malignancies. Using whole blood, individuals were screened for MGUS using a matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry assay. M-CHIP and L-CHIP were defined based on mutations in 56 genes associated with myeloid malignancies or mutations in 235 genes associated with lymphoid malignancies, respectively, using whole-exome-sequencing. Incident MM was identified using International Classification of Diseases (ICD) codes and confirmed via medical record review. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals for the association of CHIP and MGUS. Cox proportional hazards regression was used to estimate hazard ratios (HR), with time defined as date between blood sample and the first incident MM, date of death, or 12/31/2023. All models were adjusted for age and sex. Of the 16, 479 individuals (median age 65 years, 45% male) screened, 12.3% were positive for MGUS, 9.8% were positive for M-CHIP, and 1.5% were positive for L-CHIP. There was no evidence of an association between M-CHIP and MGUS (OR=0.97, 95% CI:0.83-1.12) nor L-CHIP and MGUS (OR=1.12, 95% CI:0.78-1.57), including with any MGUS isotype (IgA, IgG, or IgM; P-values>0.22). Next, we investigated the association of these precursors on risk of MM. Median follow-up was 11.7 years (range 0.1-15.1), and 46 individuals developed incident MM. Both MGUS (HR=29.7, 95% CI: 14.1-62.4) and L-CHIP (HR=4.1, 95% CI: 1.3-13.2) were significantly associated with incident MM, but not M-CHIP (HR-1.0, 95% CI: 0.4-2.5). Investigating the combined effect of MGUS and L-CHIP, we found those with MGUS alone had a 27.5-fold (95% CI: 13.0-58.2) increased risk of incident MM compared to those with no precursors. No individuals with L-CHIP only developed MM, but those with both L-CHIP and MGUS had a 129.4-fold (95% CI: 34.4-486.4) increased risk of incident MM. In the largest cohort with screened MGUS and CHIP, we found no evidence of an association between either definition of CHIP with MGUS. However, individuals with both L-CHIP and MGUS were at the highest risk of developing MM. Citation Format: Nicholas J. Boddicker, Cristine Allmer, Danelle H. Moonen, Aaron D. Norman, Angelica Macauda, Alyssa Clay-Gilmour, Janet E. Olson, Mrinal S. Patnaik, Vincent Rajkumar, Esteban Braggio, David L. Murray, Shaji Kumar, Susan L. Slager, Celine M. Vachon. The association between MGUS and CHIP and their impact on multiple myeloma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 4917.

Abstract Introduction: With over 80, 000 new cases of non-Hodgkin's lymphoma (NHL) diagnosed annually, modifiable risk factors remain unclear. Circadian disruption, linked to cancers like breast and prostate, affects immune cells (e.g., natural killer cells and T-helper cells) by altering their trafficking and proliferation. Some studies (PMID: 30566672, 27611440) suggest long sleep may increase overall NHL risk, including a Mendelian randomization (MR) study (PMID: 32895918). However, no research has examined the causal relationship between sleep traits, including long sleep, and NHL risk, by subtypes. We hypothesize that sleep traits are causally associated with NHL subtypes in European populations, testing this using MR with data from the International Lymphoma Epidemiology Consortium (InterLymph). Methods: We conducted two-sample MR using genome-wide association studies (GWAS) summary data for sleep traits and NHL subtypes. Only independent, genome-wide significant (p < 5×10-8) Single Nucleotide Polymorphisms (SNPs) were selected as valid instruments, sourced from the UK Biobank. Sleep traits included chronotype (153 SNPs), insomnia (48 SNPs), sleep duration (78 SNPs), excessive daytime sleepiness (37 SNPs), short sleep (27 SNPs), and long sleep (8 SNPs). NHL subtype data were from InterLymph: follicular lymphoma (6, 508 cases / 64, 183 controls), mantle cell lymphoma (1, 169 cases / 61, 603 controls), Waldenstrom macroglobulinemia/lymphoplasmacytic lymphoma (WM/LPL: 1, 697 cases / 59, 333 controls), and chronic lymphocytic leukemia (8, 522 cases / 67, 653 controls). The primary analysis used inverse-variance weighted (IVW) random-effects, with MR-Egger for pleiotropy adjustment. Results: Using the IVW method, we did not identify any statistically significant causal association between sleep traits and the risk of NHL subtypes, but did yield several intriguing trends. Long sleep (sleep duration >= 9 h per night) showed a marginally increased NHL risk [Odds ratio (OR) range: 1.51 - 247.9] across each subtype. Short sleep (sleep duration < 7 h per night) showed a similar trend [OR range: 1.25 - 1.81], except for WM/LPL, where we noticed a negative trend. The confidence limit for long sleep was imprecise due to the small number of variants. We did not notice any significant directional horizontal pleiotropy or weak instrument bias. Conclusions: We found no conclusive evidence of a relationship between sleep traits and four NHL subtypes. However, we see suggestive trends of long and short sleep associated with higher risk of some NHL subtypes. These findings align with the previous studies, which showed long sleep was associated with higher risk of NHL overall. Future directions include expanding our analysis to additional NHL subtypes, generating polygenic risk scores, and conducting stratified analyses. Citation Format: Pankhil Shah, Brittany Crawford, Anwar Merchant, Brenda Birmann, Angelica Macauda, Michelle A. Hildebrandt, Aaron Norman, Neil E. Caporaso, Meredith Yeager, Michael Dean, Immaculata De Vivo, Lynn Goldin, Nicola J. Camp, Rosalie Griffin, Delphine Casabonne, Federico Canzian, Pelin Unal, Elad Ziv, Catherine R. Marinac, Alexandra Nieters, Stephen J. Chanock, Mitchell J. Machiela, Michael Conry, Charlie Zhong, Hanla A. Park, Simon Cheah, Jonathan N. Hofmann, Elizabeth E. Brown, Celine Vachon, Susan Slager, Sonja Berndt, Sophia S. Wang, Vijai Joseph, James McKay, Henrik Hjalgrim, Lara Sucheston-Campbell, Karl Smith-Byrne, Alyssa Clay-Gilmour. Association between sleep traits and risk of non-Hodgkin's lymphoma subtypes: a mendelian randomization study in the International Lymphoma Epidemiology (InterLymph) Consortium [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 7415.

Genome-wide association studies (GWAS) of multiple myeloma (MM) in individuals of genetic European ancestry (EA) have identified 35 susceptibility loci. Co-heritability analyses have revealed a strong genetic correlation between MM and its precursor condition, monoclonal gammopathy of undetermined significance (MGUS). Previous research validated the association between a polygenic risk score (PRS) based on 23 MM risk loci and the risk of MGUS. Went and colleagues recently reported 12 new MM risk loci. We constructed a MM-PRS using all 35 MM risk variants and assessed its association with MGUS. Additionally, among individuals with MGUS, we assessed associations of MM-PRS with established prognostic factors [M-protein (<0.2 g/dL; 0.2-1.5 g/dL; >1.5 g/dL); isotype (IgG, IgA, IgM); free light chain (FLC) ratio (Normal (0.26-1.65), Abnormal (<0.26 or >1.65))]. Our study population included 1,723 individuals with MGUS seen at Mayo Clinic, Rochester, MN and 14,175 controls from the Mayo Clinic Biobank who lived in the 27-county region around Rochester, MN, and screened negative for MGUS using MALDI-TOF MS. The majority of MGUS patients and controls were of EA. We used logistic or multinomial regression to calculate odds ratios (ORs) and 95% confidence intervals (CIs). All models were adjusted for age, sex, and the first three principal components. The MM-PRS was significantly associated with MGUS when assessed continuously (OR, 1.28 per standard deviation (SD); 95% CI, 1.22-1.35). When comparing to the middle quintile of the MM-PRS distribution, individuals in the highest quintile had 1.41-fold increased risk of MGUS (95% CI, 1.21-1.64), and individuals in the lowest quintile had 0.68 decreased risk of MGUS (95% CI, 0.57-0.81). Of the 35 SNPs evaluated, 27 demonstrated consistent directions of effect for risk of MGUS with published risk estimates for MM. The effect of MM-PRS differed by isotype (P = 2.2 x 105), with higher MM-PRS associated with IgA MGUS (OR, 1.40 per SD; 95% CI, 1.08-1.81) and a lower MM-PRS associated with IgM (OR, 0.66 per SD; 95% CI, 0.52-0.84) compared to IgG MGUS. The MM-PRS was higher in MGUS patients with M-protein > 1.5 g/dL compared to M-protein of size 0.2-1.5 g/dL, although not statistically significant (OR, 1.41 per SD; 95% CI, 0.95-2.08); there were no differences in the MM-PRS between MGUS patients with M-protein size <0.2 g/dL compared to 0.2-1.5 g/dL (OR, 1.01 per SD;95% CI, 0.81-1.26). There were also no differences in MM-PRS by FLC ratio, (OR, 0.95 per SD; 95% CI, 0.852-1.10 for abnormal vs. normal). Our findings of an association of the MM-PRS with risk of MGUS provides further evidence for shared heritability between MGUS and MM. Future studies are needed to examine whether the expanded MM-PRS adds to existing models for MGUS progression. Citation Format: Angelica Macauda, Alyssa Clay-Gilmour, Cristine Allmer, Danelle Moonen, Aaron D. Normann, Nicholas J. Boddicker, Joselle M. Cook, Linda Baughn, Elizabeth E. Brown, Vincent S. Rajkumar, Esteban Braggio, David L. Murray, Susan L. Slager, Shaji Kumar, Celine M. Vachon. Implications of multiple myeloma polygenic risk scores (PRS) for MGUS [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 2274.

Introduction: Over 80,000 new cases of Non-Hodgkin's Lymphoma (NHL) are diagnosed annually, but modifiable risk factors are largely unknown. Circadian disruption has been identified as a risk factor for several cancers, including breast and prostate cancer, and affects immune cells such as Natural Killer (NK) cells and T-helper cells. Disruptions in circadian rhythms can lead to aberrant immune cell trafficking and proliferation. Despite this biological connection, evidence linking circadian disruptions to NHL remains inconclusive. Only single Mendelian Randomization (MR) study (PMID: 32895918) has explored the causal relationship between sleep traits and overall NHL risk, without examining NHL subtypes. Studies have shown distinct risk factor profiles among NHL subtypes. Therefore, this study aims to investigate the causal relationship between sleep traits and NHL subtypes in European populations using MR. Methods: We conducted a two-sample MR analyses using genome-wide association studies (GWAS) summary data on sleep traits and NHL subtypes. The GWAS summary statistics on sleep traits were obtained from publicly available UKBiobank datasets. The sleep traits analyzed included chronotype (N=449,734), insomnia (N=237,627), sleep duration (N=446,118), daytime sleepiness (N=452,071), short sleep (N=411,934), and long sleep (N=339,926). Only independent, genome-wide significant SNPs (p < 5×10-8) were considered as valid instruments. The NHL subtype-specific summary statistics came from the InterLymph including follicular lymphoma (FL: N=6,508 cases / 64,183 controls), mantle cell lymphoma (MCL: N=1,169 cases / 61,603 controls), and Waldenström macroglobulinemia/lymphoplasmacytic lymphoma (WM/LPL: N=1,697 cases / 59,333 controls). Our primary analysis used the random-effects inverse-variance weighted (IVW) method. In addition, we employed MR-Egger to adjust for directional pleiotropy. Heterogeneity in the MR results was evaluated with Cochran's Q statistic. Results: The IVW method did not identify any statistically significant causal association between sleep traits and the risk of NHL subtypes (FL,MCL,WM/LPL), but did yield several intriguing associations. Chronotype (early rising) showed a tendency towards a marginally an increased NHL risk across subtypes [Odds ratio (OR) range: 1.07 - 2.56]. Conversely, a decreased risk was noted for long sleep (≥ 9 hours vs. 7-8 hours) [OR range: 0.19 - 0.62], insomnia [OR range: 0.72 - 0.88], and short sleep (excludes MCL) [OR range: 0.12 - 0.72]. Confidence limit for long sleep was wide due to small number of variants (N=8). The MR-Egger regression analysis did not provide evidence of significant directional horizontal pleiotropy. The F-statistics for all sleep traits were greater than 10 indicating absences of weak instrument bias. Conclusions: Our study found no statistically significant associations between sleep traits and three NHL subtypes (FL, MCL, WM/LPL) in European populations. However, we observed potential influences of chronotype, long and short sleep, and insomnia on NHL risk. These findings align with and extend previous MR study by providing subtype-specific insights. Future directions include expanding our analysis to additional NHL subtypes, generating polygenic risk scores, and conducting stratified analyses using the InterLymph population. These efforts aim to further elucidate the complex relationship between sleep characteristics and NHL risk across diverse subtypes and populations.

AbstractWe conducted the first large genome‐wide association study to identify novel genetic variants that predict better (or poorer) prognosis in colorectal cancer patients receiving standard first‐line oxaliplatin‐based chemotherapy vs chemotherapy without oxaliplatin. We used data from two phase III trials, NCCTG N0147 and NCCTG N9741 and a population‐based patient cohort, DACHS. Multivariable Cox proportional hazards models were employed, including an interaction term between each SNP and type of treatment for overall survival (OS) and progression‐free survival. The analysis was performed for studies individually, and the results were combined using fixed‐effect meta‐analyses separately for resected stage III colon cancer (3098 patients from NCCTG N0147 and 549 patients from DACHS) and mCRC (505 patients from NCCTG N9741 and 437 patients from DACHS). We further performed gene‐based analysis as well as in silico bioinformatics analysis for CRC‐relevant functional genomic annotation of identified loci. In stage III colon cancer patients, a locus on chr22 (rs11912167) was associated with significantly poorer OS after oxaliplatin‐based chemotherapy vs chemotherapy without oxaliplatin (Pinteraction < 5 × 10−8). For mCRC patients, three loci on chr1 (rs1234556), chr12 (rs11052270) and chr15 (rs11858406) were found to be associated with differential OS (P < 5 × 10−7). The locus on chr1 located in the intronic region of RCSD1 was replicated in an independent cohort of 586 mCRC patients from ALGB/SWOG 80405 (Pinteraction = .04). The GWA gene‐based analysis yielded for RCSD1 the most significant association with differential OS in mCRC (P = 6.6 × 10−6). With further investigation into its biological mechanisms, this finding could potentially be used to individualize first‐line treatment and improve clinical outcomes.

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.