Back to Authors

Machiela, Mitchell J

0000-0001-6538-9705

Mayo Clinic

Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health, Bethesda, Maryland, USA.

Citations 10

Circulating levels of insulin‐like growth factor I ( IGF ‐I) and risk of multiple myeloma: An observational and Mendelian randomisation study

British Journal of Haematology | March 24, 2026 | Volume: 208 | Issue: 6 | Pages: 1942-1953

10.1111/bjh.70444 | PMID 41873017

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.
Genetically determined body mass index is associated with diffuse large B‐cell lymphoma in polygenic and Mendelian randomization analyses

International Journal of Cancer | September 5, 2025 | Volume: 158 | Issue: 1 | Pages: 45-59

10.1002/ijc.70039 | PMID 40910475

Abstract Obesity has been associated with non‐Hodgkin lymphoma (NHL), but the evidence is inconclusive. We examined the association between genetically determined adiposity and four common NHL subtypes: diffuse large B‐cell lymphoma (DLBCL), follicular lymphoma, chronic lymphocytic leukemia, and marginal zone lymphoma, using eight genome‐wide association studies of European ancestry ( N = 10,629 cases, 9505 controls) and constructing polygenic scores for body mass index (BMI), waist‐to‐hip ratio (WHR), and waist‐to‐hip ratio adjusted for BMI (WHRadjBMI). Higher genetically determined BMI was associated with an increased risk of DLBCL [odds ratio (OR) per standard deviation (SD) = 1.18, 95% confidence interval (95% CI): 1.05–1.33, p = .005]. This finding was consistent with Mendelian randomization analyses, which demonstrated a similar increased risk of DLBCL with higher genetically determined BMI (OR per SD = 1.12, 95% CI: 1.02–1.23, p = .03). No significant associations were observed with other NHL subtypes. Our study demonstrates a positive link between a genetically determined BMI and an increased risk of DLBCL, providing additional support for increased adiposity as a risk factor for DLBCL.
Abstract 7415: Association between sleep traits and risk of non-Hodgkin's lymphoma subtypes: a mendelian randomization study in the International Lymphoma Epidemiology (InterLymph) Consortium

Cancer Research | April 21, 2025 | Volume: 85 | Issue: 8_Supplement_1 | Pages: 7415-7415

10.1158/1538-7445.am2025-7415

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.
Association between Sleep Traits and Risk of Non-Hodgkin's Lymphoma Subtypes: A Mendelian Randomization Study in the International Lymphoma (InterLymph) Epidemiology Consortium

Blood | December 5, 2024 | Volume: 144 | Issue: Supplement 1 | Pages: 6182-6182

10.1182/blood-2024-199373

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.
Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization

International Journal of Molecular Sciences | May 10, 2023 | Volume: 24 | Issue: 10 | Pages: 8500

10.3390/ijms24108500 | PMID 37239846

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.
A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk

International Journal of Cancer | September 9, 2022 | Volume: 152 | Issue: 2 | Pages: 239-248

10.1002/ijc.34278 | PMID 36082445

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.
Does a Multiple Myeloma Polygenic Risk Score Predict Overall Survival of Patients with Myeloma?

Cancer Epidemiology, Biomarkers & Prevention | June 14, 2022 | Volume: 31 | Issue: 9 | Pages: 1863-1866

10.1158/1055-9965.epi-22-0043 | PMID 35700034

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.
Coinherited genetics of multiple myeloma and its precursor, monoclonal gammopathy of undetermined significance

Blood Advances | June 22, 2020 | Volume: 4 | Issue: 12 | Pages: 2789-2797

10.1182/bloodadvances.2020001435 | PMID 32569378

Abstract So far, 23 germline susceptibility loci have been associated with multiple myeloma (MM) risk. It is unclear whether the genetic variation associated with MM susceptibility also predisposes to its precursor, monoclonal gammopathy of undetermined significance (MGUS). Leveraging 2434 MM cases, 754 MGUS cases, and 2 independent sets of controls (2567/879), we investigated potential shared genetic susceptibility of MM and MGUS by (1) performing MM and MGUS genome-wide association studies (GWAS); (2) validating the association of a polygenic risk score (PRS) based on 23 established MM loci (MM-PRS) with risk of MM, and for the first time with MGUS; and (3) examining genetic correlation of MM and MGUS. Heritability and genetic estimates yielded 17% (standard error [SE] ±0.04) and 15% (SE ±0.11) for MM and MGUS risk, respectively, and a 55% (SE ±0.30) genetic correlation. The MM-PRS was associated with risk of MM when assessed continuously (odds ratio [OR], 1.17 per SD; 95% confidence interval [CI], 1.13-1.21) or categorically (OR, 1.70; 95% CI, 1.38-2.09 for highest; OR, 0.71; 95% CI, 0.55-0.90 for lowest compared with middle quintile). The MM-PRS was similarly associated with MGUS (OR, 1.19 per SD; 95% CI, 1.14-1.26 as a continuous measure, OR, 1.77, 95%CI: 1.29-2.43 for highest and OR, 0.70, 95%CI: 0.50-0.98 for lowest compared with middle quintile). MM and MGUS associations did not differ by age, sex, or MM immunoglobulin isotype. We validated a 23-SNP MM-PRS in an independent series of MM cases and provide evidence for its association with MGUS. Our results suggest shared common genetic susceptibility to MM and MGUS.
Association between a Polygenic Risk Score for Multiple Myeloma Risk and Overall Survival

Blood | November 14, 2019 | Volume: 134 | Issue: Supplement_1 | Pages: 4366-4366

10.1182/blood-2019-126088

Background Genome-wide association studies (GWAS) conducted in populations of European ancestry (EA) have identified and confirmed 23 germline susceptibility loci for multiple myeloma (MM). The effect sizes of single nucleotide polymorphisms (SNPs) at these loci are small, therefore combining them into a single summary measure, known as a polygenic risk score (PRS), may provide a more meaningful risk factor. We have previously shown a PRS comprised of the 23 SNPs for MM contributes to increased risk of MM, with a 2.7-fold increase for highest vs. lowest PRS quintiles. Whether the MM-PRS is also associated with overall survival (OS) in MM cases has not been evaluated. We examined the association between MM-PRS and OS in two EA studies. Methods The first study consisted of 2,179 EA MM cases from ten studies included in the Multiple Myeloma Working Group within the International Lymphoma Consortium (InterLymph). Cases were diagnosed between 1970 and 2015 and genotyped using multiple platforms (Oncoarray, Affymetrix, Human660W-quad Beadchip, and Illumina arrays); 885 cases also had stage [based on International Staging System (ISS)] available. Each of the GWAS was subjected to rigorous standard quality control independently (prior to imputation via the Michigan imputation server based on the Haplotype Reference Consortium (HRC). The second study consisted of 515 newly diagnosed EA MM cases from CoMMpass (Relating Clinical Outcomes in Multiple Myeloma to Personal Assessment of Genetic Profile), diagnosed from 2011-2013, who had whole genome sequencing (WGS) performed on germline DNA. The WGS data was used to call common germline genetic variants through the Mayo Clinic bioinformatics pipeline. Briefly, genetic variants were detected with GenomeGPS, aligned to the hg19 reference genome, called using the GATK (V3.6) Haplotype Caller, and merged for multiple-sample joint calling. To reduce the false positive variants, variant quality score recalibration (VQSR) was applied for both SNPs and INDELs. After quality control, 458 EA samples remained. Follow-up was available for both studies and consisted of time from MM diagnosis date until death or date of last known follow-up. The PRS was constructed from the 23 MM SNPs using the published per allele odds ratio associated with MM risk. The published log odds ratios for each SNP were multiplied by the number of risk alleles (0, 1, 2) for the corresponding SNP, and summed, resulting in a unique score per person. Kaplan-Meier curves and Cox proportional hazard models were used to assess the association between PRS with MM OS considering two models: 1) adjusted for age, sex, study and 2) additional adjustment by stage (ISS). Hazard ratios (OR) and 95% confidence intervals (CI) were estimated. The PRS was evaluated both as a continuous variable, per standard deviation (SD), and as a categorical variable (quintiles). Results MM cases (N=2,179) in the InterLymph study were 59% male and 41% female and the median age was 61.0 years (26-90 years). Median follow-up time was 57.2 months (1.0-509.0 months) with 868 reported deaths. MM cases with stage information available consisted of 20% stage I (n=178), 53% stage II (n=466), and 27% stage III (n=241). No association was observed between PRS and OS in MM patients regardless of adjustment for stage (continuous PRS (HR: 1.03, 95% CI: 0.83-1.28, P=0.80) or by quintile PRS (p>0.05)) (Table). The CoMMpass EA MM cases (n=458) had similar distributions for sex (61% male and 39% females) but were slightly older 65 years (27-93 years) and had shorter follow-up time (median=39.75 months (0.13-77.2)) with 117 deaths. Stage was available for 96% of CoMMpass cases including 36% stage I (n=159), 33% stage II (n=146), and 31% stage III (n=134). We also observed no association of PRS and OS in the CoMMpass study (HR=1.02, 95% CI: 0.72 -1.46, P= 0.89), adjusted for age, sex, and stage (Table). Discussion A PRS score for MM risk is not associated with OS for MM cases in two EA populations. Given that prior studies have shown association of genetic variation with MM survival, efforts to identify additional loci associated with OS or MM specific survival are warranted. Future studies should also consider germline variants impact on molecular subtypes, specific therapies, and outcomes. Disclosures Kumar: Celgene: Consultancy, Research Funding; Janssen: Consultancy, Research Funding; Takeda: Research Funding.
Body mass index and breast cancer survival: a Mendelian randomization analysis

International Journal of Epidemiology | June 29, 2017 | Volume: 46 | Issue: 6 | Pages: 1814-1822

10.1093/ije/dyx131 | PMID 29232439