Causal relationship between particulate matter 2.5 and diabetes: two sample Mendelian randomization
- 1Department of Dermatovenereology, Chengdu Second People's Hospital, Chengdu, Sichuan, China
- 2Health Science Center, Yangtze University, Jingzhou, Hubei, China
- 3Department of Dermatology, Yunnan Provincial Hospital of Traditional Chinese Medicine, The First Affiliated Hospital of Yunnan University of Chinese Medicine, Kunming, Yunnan, China
A Commentary on
Causal relationship between particulate matter 2.5 and diabetes: two sample Mendelian randomization
by Kim, J. M., Kim, E., Song, D. K., Kim, Y.-J., Lee, J. H., and Ha, E. (2023). Front. Public Health 11:1164647. doi: 10.3389/fpubh.2023.1164647
Diabetes mellitus (DM) is a significant metabolic disorder characterized by chronic hyperglycemia and is showing an increasing trend globally (1). Recent research has indicated that particulate matter 2.5 (PM2.5) can induce a pro-inflammatory state, increasing oxidative stress, thereby elevating the risk of DM, a condition detectable through the measurement of homocysteine levels (2–4). However, due to the inherent limitations of traditional research methodologies, the causal relationship remains unconfirmed. Mendelian randomization (MR) uses single nucleotide polymorphisms (SNPs) from genome-wide association studies (GWASs) as instrumental variables (IVs) to assess the impact of exposure on outcomes. This method can avoid the influence of confounding factors and reverse results (5). Consequently, we were intrigued by the recent study conducted by Kim et al., which employed the two-sample MR (TSMR) approach in investigating the causal relationship between PM2.5 and DM (6).
Kim et al.'s TSMR analysis primarily indicated that genetic susceptibility to PM2.5 is associated with a higher risk of DM. However, the dataset used for the analysis had significant shortcomings, leading to potential false-positive results and considerable bias in their findings. Specifically, the DM dataset originated from the UK Biobank (UKB), and notably, the exposure dataset was also sourced from UKB. This contravenes the principles of TSMR studies based on summary-level GWAS data, with an astonishingly high sample overlap rate of 91.81% calculated. Additionally, since Type 1 DM (T1DM) is an autoimmune disease and existing observational studies primarily focus on PM2.5 and Type 2 DM (T2DM) (4), the MR study's use of generalized DM data is problematic. The UKB describes DM as “Has a doctor ever told you that you have diabetes?” which lacks specificity in diabetes typing. Furthermore, according to the publication date of the article, the most up-to-date and comprehensive GWAS data should be used to ensure scientific advancement, a detail the authors overlooked. We reanalyzed the data using Kim et al.'s method and other approaches, differentiating DM types. We obtained the most current and comprehensive data on European ancestry T2DM (80,154 cases/853,816 controls) from the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) consortium (7), and T1DM data (4,196 cases/308,252 controls) from the FinnGen consortium (8). Sensitivity analyses were conducted which included additional calculations of I2 to assess the suitability of fixed effects in the inverse-variance weighted (IVW) method (Table 1). Following the application of the Bonferroni correction (P < 0.05/2), our MR analysis validated a significant causal association between PM2.5 and T2DM (Figure 1). Discrepancies in the outcomes can be attributed to the inadequate consideration of the DM phenotype and misunderstanding of the TSMR principle in the initial MR study conducted by Joyce Mary Kim et al., which subsequently resulted in false-positive findings. We acknowledge the challenges often faced when analyzing risk factors such as PM2.5 and appreciate the authors' endeavors in advancing our comprehension of the intricate connection between PM2.5 and DM. However, in this intricate field, precise identification of datasets and adherence to the TSMR principle are imperative in order to substantiate conclusions regarding potential causality between clinical characteristics and/or diseases.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1353982/fpubh-12-1353982-HTML/image_m/fpubh-12-1353982-g001.jpg)
Figure 1. Summary of estimates of genetically predicted PM2.5 for causality of T1DM or T2DM as assessed by IVW, the primary method of Mendelian randomization. MR, Mendelian randomization; IVW, Inverse-Variance-Weighted; CML, Conditional Maximum likelihood; OR, odds ratio; CI, confidence interval; T2DM, Type 2 Diabetes Mellitus; T1DM, Type 1 Diabetes Mellitus; SNP, Single Nucleotide Polymorphism; PM2.5, particulate matter 2.5.
In the analysis of MR, the F-statistic plays a crucial role in evaluating the instrumental variables (IVs) and their ability to explain the exposure variable. Additionally, the variables R2 and POWER are important indicators used to measure the proportion of explained variance and the statistical test's power, respectively. In MR study, the presence of weak IVs can potentially skew the estimates and results. To mitigate this issue, researchers rely on the F-statistic, which should ideally exceed a value of 10 to ensure the IVs possess sufficient strength for accurate estimation of causal associations between exposure and outcome (9). Additionally, ensuring the validity and dependability of outcomes requires adequate statistical power. Insufficient power frequently leads to an incapability of precisely identifying the extent of causal consequences, amplifying the likelihood of erroneous positive findings. Hence, power analysis is of utmost significance. Considering that the majority of genetic variations account for a minute portion of phenotypic variability, statistical power is considered a prominent obstacle in MR investigations (10). However, the authors, Kim et al., overlooked these crucial factors, and the absence of such critical statistical information could potentially lead to misinterpretation of the study findings. This, in turn, may have a significant impact on the overall quality and credibility of the research conducted. To address this concern, we conducted additional calculations, including evaluating the F-statistic [F = R2 × (N - 2)/(1 - R2)] (9), determining the R2 value [R2 = 2 × MAF × (1 - MAF) × beta2, where MAF represents the minor allele frequency for each single nucleotide polymorphism] (11), and assessing statistical power (https://shiny.cnsgenomics.com/mRnd/). The results revealed that in the analysis of PM2.5 with T1DM and T2DM, six (Outlier:rs77205736) and seven IVs were identified (Figure 1), respectively. The average F-statistic for both exceeded 37, indicating robust instrument strength. The genetic variance explained by PM2.5 was 0.06% for T1DM and 0.05% for T2DM. Using odds ratios derived from the IVW method, the power values for detecting associations in PM2.5 with T1DM and T2DM were 96 and 100%, respectively. This confirms the reliability of the causality.
This commentary thoroughly investigates the intricate connections between PM2.5 and DM. The original study conducted by Kim et al. established a crucial groundwork for comprehending these associations through MR analysis. They posit that there exists a genetic-level causal relationship between PM2.5 and diabetes in the European population. However, our reevaluation, utilizing refined definitions of phenotypes and expanded analytical techniques, which encompass meticulous assessment of instrumental variables and power analysis, presents a more nuanced standpoint. Our findings suggest that PM2.5 lacks a noteworthy causal association with type 1 diabetes, but does elevate the risk of type 2 diabetes. In order to uphold scientific transparency, we express our gratitude to the editors of “Frontiers in Public Health” and the authors of the paper for recognizing and addressing our concerns.
Author contributions
YN: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Supervision, Writing – original draft, Writing – review & editing. YZ: Conceptualization, Formal analysis, Investigation, Methodology, Resources, Supervision, Writing – original draft, Writing – review & editing. JY: Conceptualization, Formal analysis, Methodology, Resources, Supervision, Writing – review & editing. XY: Conceptualization, Formal analysis, Methodology, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
Special thanks to the related investigators for sharing the GWAS summary statistics included in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Zhang Y, Tang Z, Tong L, Wang Y, Li L. Serum uric acid and risk of diabetic neuropathy: a genetic correlation and mendelian randomization study. Front Endocrinol. (2023) 14:1277984. doi: 10.3389/fendo.2023.1277984
2. Liu L, Zheng B, Jiang Z, Wu S, Jin Q, Lin P, et al. Association of elevated cord blood oxidative stress biomarkers with neonatal outcomes in mothers with pre-eclampsia: a case-control study. Gynecol Obstet Invest. (2021) 86:361–9. doi: 10.1159/000517996
3. Li S, Guo B, Jiang Y, Wang X, Chen L, Wang X, et al. Long-term exposure to ambient PM25 and its components associated with diabetes: evidence from a large population-based cohort from China. Diabet Care. (2023) 46:111–9. doi: 10.2337/dc22-1585
4. Rajagopalan S, Brook RD. Air pollution and type 2 diabetes: mechanistic insights. Diabetes. (2012) 61:3037–45. doi: 10.2337/db12-0190
5. Tin A, Kottgen A. Mendelian randomization analysis as a tool to gain insights into causes of diseases: a primer. J Am Soc Nephrol. (2021) 32:2400–7. doi: 10.1681/ASN.2020121760
6. Kim JM, Kim E, Song DK, Kim Y-J, Lee JH, Ha E. Causal relationship between particulate matter 25 and diabetes: two sample Mendelian randomization. Front Public Health. (2023) 11:1164647. doi: 10.3389/fpubh.2023.1164647
7. Mahajan A, Spracklen CN, Zhang W, Ng MCY, Petty LE, Kitajima H, et al. Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation. Nat Genet. (2022) 54:560–72. doi: 10.1038/s41588-022-01058-3
8. Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. (2023) 613:508–18. doi: 10.1038/s41586-022-05473-8
9. Teumer A. Common methods for performing Mendelian randomization. Front Cardiovasc Med. (2018) 5:51. doi: 10.3389/fcvm.2018.00051
10. Brion M-JA, Shakhbazov K, Visscher PM. Calculating statistical power in Mendelian randomization studies. Int J Epidemiol. (2013) 42:1497–501. doi: 10.1093/ije/dyt179
Keywords: Mendelian randomization, particulate matter 2.5, diabetes, environmental epidemiology, genome-wide association study (GWAS)
Citation: Ni Y, Zhang Y, Ye J and Yang X (2024) Commentary: Causal relationship between particulate matter 2.5 and diabetes: two sample Mendelian randomization. Front. Public Health 12:1353982. doi: 10.3389/fpubh.2024.1353982
Received: 11 December 2023; Accepted: 13 February 2024;
Published: 26 February 2024.
Edited by:
Bo Zhu, Boston Children's Hospital and Harvard Medical School, United StatesReviewed by:
Cosmin Mihai Vesa, University of Oradea, RomaniaCopyright © 2024 Ni, Zhang, Ye and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianzhou Ye, kmyjz63@sina.com; Xuesong Yang, yangxuesong@ynutcm.edu.cn
†These authors share first authorship