Skip to main content

Mendelian randomisation highlights type 1 diabetes as a causal determinant of idiopathic pulmonary fibrosis

Abstract

Background

It is unclear whether type 1 diabetes (T1D) causes idiopathic pulmonary fibrosis (IPF), despite observational research linking the two conditions. Therefore, our study aimed to examine the causal link between T1D and the likelihood of IPF by employing the Mendelian randomization (MR) technique of two-sample Mendelian randomization.

Methods

Using data from two genome-wide association studies (GWAS) with European ancestry, we performed a two-sample MR analysis. These studies involved 18,856 individuals (6,683 cases and 12,173 controls) for T1D and 198,014 individuals (10,028 cases and 196,986 controls) for IPF. We utilized inverse-variance weighted (IVW) analysis as our main approach to determine the association between the risk of IPF and T1D. To evaluate multidirectionality, the MR-Egger regression test was utilized, whereas heterogeneity was assessed using Cochran’s Q test. Additionally, a leave-one-out analysis was performed to assess the reliability of the results.

Results

38 SNPs linked to T1D were employed as instrumental variables (IVs). Multiple MR methods yielded consistent results, and the MR analysis reveals a significant and positive causal impact of T1D on IPF (MR-IVW, odds ratio [OR] = 1.128, 95% confidence interval [CI] 1.034–1.230; P = 0.006). The limitations of the study include the lack of data from non-European groups and the inability to rule out the possibility of small links. Larger MR experiments are necessary to investigate minute impacts.

Conclusions

The results of this study provide evidence that T1D contributes to the onset and advancement of IPF. This finding may provide important insights into the cause of IPF and possible treatments in the future.

Introduction

IPF is an incurable and worsening respiratory condition of unknown cause, characterized by the progressive deterioration of lung function. The median life expectancy for individuals with IPF is a mere 3.8 years [1, 2]. The progression of disease involves damage to alveolar epithelial cells, a transition from epithelial to mesenchymal cells, and the activation of fibroblasts, resulting in an elevated accumulation of extracellular matrix proteins. Alveolar architecture is remarkably disrupted and the extracellular cell matrix is changed, replacing normal healthy lung tissues. These alterations eventually result in a significant disturbance of the carefully regulated gas exchange mechanism and a decrease in lung compliance, which finally induce respiratory failure and mortality [3]. There is no cure for IPF and treatment options are limited. Therefore, a comprehensive comprehension of coexisting conditions in the IPF population, specifically their correlation, has the potential to greatly enhance medical treatment, ultimately leading to better chances of survival and improved quality of life [4]. Common comorbidities in IPF include respiratory diseases (e.g., emphysema, COPD, obstructive sleep apnea, and lung cancer), cardiovascular diseases (e.g., ischemic heart disease, pulmonary hypertension, and coronary artery disease), gastrointestinal diseases (e.g., gastroesophageal reflux disease and esophageal hiatal hernia), and metabolic diseases (e.g., diabetes mellitus and hypothyroidism) [5]. The occurrence of T1D is significantly greater among individuals with IPF compared to the overall populace, and T1D is linked to an elevated susceptibility to IPF [6]. There are clinical studies showing significant clinical and histopathological correlations between IPF and T1D.The formation of fibrosis is thought to be associated with reactive oxygen species and advanced glycation end products that arise from high blood sugar levels [7]. T1D and IPF are characterised by disorders of the immune system that can trigger chronic inflammation in any organ system. Abnormal leukocyte telomere length has the potential to increase the risk of T1D and IPF [8]. Therefore, it can serve as a predictor and may provide new potential therapeutic targets for T1D and IPF. However, altered leukocyte telomere length may not be a direct cause of T1D and IPF [8]. Previous observational studies investigating the relationship between endocrine and metabolic factors and IPF have yielded inconsistent results. Therefore, there is a need for a more comprehensive study of the effects of endocrine metabolic factors on IPF using Mendelian randomization (MR) analysis [9].

Although randomized controlled trials are considered the benchmark for evaluating causality, their application is restricted due to practical constraints, expenses, and ethical factors [10]. IPF is relatively rare, making it unfeasible to collect a large sample in a longitudinal study to adequately examine the relationship between T1D and IPF. In addition, many known and unknown confounding factors pose a great challenge in studying the causal effect of T1D on IPF. The emergence of causal inferences between exposures and observational findings offers this possibility. The study of causality using genetic instruments (e.g., single nucleotide polymorphisms [SNPs]) can be effectively conducted with the help of MR framework [11]. It is similar to a “natural” randomized controlled trial, where the random assignment of genetic alleles affecting exposure largely eliminates the effects of unobserved confounders and avoids reverse causality and measurement error, which may affect other study designs [12]. The availability of publicly accessible GWAS data has resulted in the generation of comprehensive summary-level information. As a result, two-sample MR has emerged as a cost-effective and efficient approach to investigate causal connections between health risk factors and disease outcomes [13].

Using publicly accessible GWAS summary statistics, we investigated the causal association between genetically predicted T1D and IPF within an MR framework in this study. The results may offer additional evidence for the etiology of IPF.

Method

The methodology and the origin of the data

We performed a two-sample MR analysis using combined GWAS data to evaluate the possible causal association between T1D and IPF. All GWAS studies included unrelated individuals of European ancestry. A total of 18,856 participants (6,683 cases and 12,173 controls) from the European Bioinformatics Institute (EBI) were part of this study, and genome-wide genotypes were measured. The statistical data for IPF were acquired from the GWAS summarized by Dhindsa et al. [14] and colleagues. The GWAS analysis consisted of 1028 IPF patients and 196,986 controls. To estimate and test for causal effects on outcomes (i.e., IPF), the MR framework employed independent instrumental SNPs as instrumental variables (IVs) for exposures (e.g., T1D). To ensure the validity of instrumental variables, three modeling assumptions are necessary for standard MR analysis. (i) IV is associated with T1D at a genome-wide significant level; (ii) IV must be independent of any confounding factors; (iii.) IV affects IPF only through T1D (Fig. 1).

Fig. 1
figure 1

Study design of univariable MR to identify the causal association between T1D and IPF. IPF, idiopathic pulmonary fibrosis; MR, Mendelian randomisation

Genetic instrumental variables

In order to satisfy the initial hypothesis of the MR analysis, which states that instrumental variables (IVs) have a strong connection with T1D biomarkers, we carefully chose independent IVs that exhibited statistically significant correlations with T1D at the genome-wide level (P < 5 × 10 − 8, linkage disequilibrium < 0.001, genetic distance = 10, 000 KB) [15]. To avoid potential confounding effects of genetic variation, we employed RStudio 4.3.2 and the TwoSampleMR package to detect 38 SNPs that are significantly linked to T1D (Additional file 1), meeting the initial hypothesis. To evaluate if the IVs included were linked to any recognized confounders, we conducted a search in the PhenoScanner database (http://www.phenoscanner.medschl.cam.ac.uk). Risk factors for IPF have been identified, such as smoking, dust exposure, and reflux esophagitis [1]. Hence, in this research, we eliminated the SNP associated with smoking (rs3184504), and subsequently collected information from IPF GWAS for 37 out of the 38 mentioned SNPs. Finally, we used 37 SNPs as instrumental variables for T1D in our study. Furthermore, we computed the F-value (F = beta^2/se^2) [16] to verify that the included independent variables are not affected by feeble IVs [17].

Statistical analysis

Many robust statistical methods were utilized to ensure the accuracy and dependability of the findings, and sensitivity analyses were conducted to assess the potential impact of various sources of bias. To assess the impact of a 1 standard deviation (SD) rise in standardized logarithmically transformed T1D on IPF vulnerability, the primary analytical approach employed was the MR inverse variance weighting (MR-IVW) technique [18].We conducted four additional sensitivity analyses using the weighted-median approach [19], the MR-Egger method [20], the weighted mode [21], and the simple mode [22] to assess the strength of the results. Reliable estimates were obtained using the weighted median method when valid IV accounted for over 50% of the information. We evaluated various instrumental variables (IVs) to examine horizontal pleiotropy, employing the MR-Egger method [20]. Furthermore, we investigated the diversity among particular independent variables by utilizing Cochrane’s Q-value [23]. Additionally, we conducted a sensitivity analysis by leaving out one SNP at a time to investigate any potential unequal impacts of individual SNPs on the overall estimation. P-values ranging from the Bonferroni correction to 0.05 were deemed indicative of a potential association warranting additional scrutiny. We utilized the “MendelianRandomization” [24], “MRPRESSO“ [25], and “TwoSampleMR” R packages for the analysis, which was carried out in the R environment (version 4.3.2). A significance threshold of 0.05 was applied.

Results

To assess the genetic relationship between T1D and IPF, we obtained 37 SNPs as IVs and used them in an MR analysis. The IVW method yielded significant evidence of a causal link between T1D and the likelihood of IPF, with an OR of 1.1137 (MR-IVW) and a 95% CI of 1.019089–1.217114 (P = 0.017434) (Table 1). This indicates that a self-reported history of T1D could lead to an average 11.37% increased risk of IPF, as depicted in the scatter plot (Fig. 2).The F-statistics > 10 for each SNP suggested a low level of instrumental variable bias. Furthermore, our results were validated through a sensitivity analysis, which revealed no variation among independent variables according to the Cochran’s Q test results (PIVW = 0.282977, PMR Egger = 0.259668, Table 2). The symmetrical nature of the funnel plot provided additional evidence for the absence of heterogeneity(Fig. 3).Moreover, the MR-Egger regression findings indicated the absence of significant pleiotropy for any instrumental variable (P = 0.54392, Table 2). This implies that the influence of instrumental variables on IPF through mechanisms unrelated to T1D is unlikely. Consistently similar results were obtained from the leave-one-out sensitivity analysis, even when removing one SNP at a time (Fig. 4). Overall, several MR techniques consistently estimated the causal link between T1D and IPF.

Table 1 MR analysis results between T1D and IPF.
Fig. 2
figure 2

Scatter plot for the causal effect of T1D on IPF risk The extent of the cause-and-effect relationship is demonstrated by the incline of the linear graph…

Table 2 Heterogeneity and Pleiotropy tests of MR.
Fig. 3
figure 3

Funnel plot for the the general heterogeneity in the impact of T1D on the risk of IPF.

Fig. 4
figure 4

Forest plot for leave-one-out analysis, where each point on the left side represents the causal effect determined by IVW after removing the specific SNP. IVW, inverse variance weighted method

Discussion

Using a two-sample MR design, we demonstrate the causal connection between T1D and IPF in this research. The results are robust and convincing, supported by the confirmation of the findings through different MR approaches with diverse modeling assumptions and an independent validation set. Furthermore, no indications of reverse causation were found in the two-way investigation. This study benefited from the extensive GWAS sample size, including 1028 cases and 196,986 controls, and fully utilized MR analysis with genetic instruments. Our research confirms the results of previous observational studies that associated T1D with pulmonary fibrosis [6].

Mounting evidence suggests that diabetes can exert a substantial influence on pulmonary health. The cause of lung damage in diabetic individuals is still a topic of debate. It is suggested that the development of T1D may involve complex underlying processes IPF. Studies have shown that a significant increase in blood sugar can lead to a gradual decline in the function of pancreatic β-cells and an increase in insulin resistance. This process is referred to as “glucose toxicity“ [26]. This can trigger oxidative stress, activate the JNK pathway, and potentially lead to lung fibroblast activation, further promoting pulmonary fibrosis [27]. Another study suggests that diabetes-related interactions between platelets and endothelial cells may contribute to oxidative damage and persistent vascular inflammation, leading to pulmonary fibrosis [28]. TGF-β’s pro-fibrotic characteristics are mediated by connective tissue growth factor (CTGF), which is crucial for fibrosis and tissue remodeling. The kidney, heart, liver, skin, and lungs are among the fibrotic organs where CTGF expression has been shown to be upregulated [29]. A study on the detrimental consequences of diabetes on pulmonary fibrosis used a streptozotocin-induced diabetic rat (STZ rat) model, in which the lung tissues had a considerably higher expression level of CTGF. Furthermore, the levels of CTGF and transcription activator in the lung tissues of diabetic rats were reversed through the treatment of hyperglycemia [30]. A study by Fariña et al. [31] suggests that diabetes can disrupt the alveolar-capillary barrier’s basement membrane, leading to excessive collagen and extracellular matrix deposition, as well as inflammatory cell infiltration, all of which can exacerbate fibrotic changes.

There are multiple advantages to our research. This is the initial study to investigate the possible connection between T1D and the likelihood of IPF, establishing a cause-and-effect association between the two ailments. Additionally, the GWAS data from two large European population samples formed the basis of this MR analysis, providing sufficient power to establish causality. Furthermore, T1D has a lasting influence on the risk of IPF, as demonstrated through MR analysis, and this impact is improbable to be affected by any confounding variables.

Nevertheless, the study has its constraints. Our findings may not be applicable to populations of diverse ethnic backgrounds, as they primarily rely on participants of European descent. Additionally, while horizontal pleiotropy was not detected, the precise purpose of most of these SNPs remains unclear. This reduces statistical efficiency by increasing the residuals, even though it may not significantly affect overall causality. Ultimately, the frequency of IPF rises considerably with advancing age [32], and relying on combined information from GWAS rather than individual-level data complicates the establishment of a direct cause-and-effect connection between T1D and age-specific IPF. Furthermore, as this information was not relevant to each person in all datasets, the IPF GWAS in this research could not be modified for age and gender. In conclusion, the MR findings suggest that genetic predisposition to T1D has an influence on the risk of IPF throughout one’s lifetime. Consequently, additional investigation is required to ascertain the immediate consequences of T1D on the risk of IPF.

Conclusions

In conclusion, our study provides compelling evidence that the self-reported history of T1D is a causative factor that influences the risk of IPF. This discovery may offer new insights into the etiology of IPF. Importantly, T1D is a medically significant mechanism that requires careful monitoring in IPF patients. However, treating IPF based on these findings may not be easy, and further pathological and biochemical research is definitely needed to fully comprehend the complex connection between T1D and an increased risk of IPF.

Data availability

No datasets were generated or analysed during the current study.

References

  1. Richeldi L, Collard HR, Jones MG. Idiopathic pulmonary fibrosis. Lancet. 2017;389(10082):1941–52.

    Article  PubMed  Google Scholar 

  2. Neumark N, Cosme C Jr, Rose KA, et al. The idiopathic pulmonary fibrosis cell Atlas. Am J Physiol Lung Cell Mol Physiol. 2020;319(6):L887–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Lederer DJ, Martinez FJ. Idiopathic pulmonary fibrosis. N Engl J Med. 2018;378(19):1811–23.

    Article  CAS  PubMed  Google Scholar 

  4. Raghu G, Amatto VC, Behr J, Stowasser S. Comorbidities in idiopathic pulmonary fibrosis patients: a systematic literature review. Eur Respir J. 2015;46(4):1113–30.

    Article  CAS  PubMed  Google Scholar 

  5. Oldham JM, Collard HR. Comorbid conditions in Idiopathic Pulmonary Fibrosis: Recognition and Management. Front Med (Lausanne). 2017;4:123.

    Article  PubMed  Google Scholar 

  6. Enomoto T, Usuki J, Azuma A, Nakagawa T, Kudoh S. Diabetes mellitus may increase risk for idiopathic pulmonary fibrosis. Chest. 2003;123(6):2007–11.

    Article  PubMed  Google Scholar 

  7. Kim YJ, Park JW, Kyung SY, et al. Clinical characteristics of idiopathic pulmonary fibrosis patients with diabetes mellitus: the national survey in Korea from 2003 to 2007. J Korean Med Sci. 2012;27(7):756–60.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Liu M, Luo P, Liu L, Wei X, Bai X, Li J, Wu L, Luo M. Immune-mediated inflammatory diseases and leukocyte telomere length: a mendelian randomization study. Front Genet. 2023;14:1129247.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Jiang Y, Chen R, Xu S, Ding Y, Zhang M, Bao M, He B, Li S. Endocrine and metabolic factors and the risk of idiopathic pulmonary fibrosis: a mendelian randomization study. Front Endocrinol (Lausanne). 2024;14:1321576.

    Article  PubMed  Google Scholar 

  10. Bothwell LE, Podolsky SH. The emergence of the Randomized, Controlled Trial. N Engl J Med. 2016;375(6):501–4.

    Article  PubMed  Google Scholar 

  11. Sekula P, Del Greco MF, Pattaro C, Köttgen A. Mendelian randomization as an Approach to assess causality using Observational Data. J Am Soc Nephrol. 2016;27(11):3253–65.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Smith GD, Ebrahim S. Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.

    Article  PubMed  Google Scholar 

  13. Burgess S, Butterworth A, Malarstig A, Thompson SG. Use of mendelian randomisation to assess potential benefit of clinical intervention. BMJ. 2012;345:e7325.

    Article  PubMed  Google Scholar 

  14. Dhindsa RS, Mattsson J, Nag A, et al. Identification of a missense variant in SPDL1 associated with idiopathic pulmonary fibrosis. Commun Biol. 2021;4(1):392.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Chen Z, Boehnke M, Wen X, Mukherjee B. Revisiting the genome-wide significance threshold for common variant GWAS. G3 (Bethesda). 2021;11(2):jkaa056.

    Article  PubMed  Google Scholar 

  16. Burgess S, Thompson SG, CRP CHD Genetics Collaboration. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.

    Article  PubMed  Google Scholar 

  17. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, Davey Smith G, Sterne JA. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21(3):223–42.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some Invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Zhu G, Zhou S, Xu Y, Gao R, Li H, Zhai B, Liu X, He Y, Wang X, Han G, Su W, Wang R. Mendelian randomization study on the causal effects of omega-3 fatty acids on rheumatoid arthritis. Clin Rheumatol. 2022;41(5):1305–12.

    Article  PubMed  Google Scholar 

  23. Hemani G, Bowden J, Davey Smith G. Evaluating the potential role of pleiotropy in mendelian randomization studies. Hum Mol Genet. 2018;27(R2):R195–208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Yavorska OO, Burgess S. MendelianRandomization: an R package for performing mendelian randomization analyses using summarized data. Int J Epidemiol. 2017;46(6):1734–9.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Kaneto H, Katakami N, Kawamori D, Miyatsuka T, Sakamoto K, Matsuoka TA, Matsuhisa M, Yamasaki Y. Involvement of oxidative stress in the pathogenesis of diabetes. Antioxid Redox Signal. 2007;9(3):355–66.

    Article  CAS  PubMed  Google Scholar 

  27. Tang R, Zhou Y, Mei S, Xu Q, Feng J, Xing S, Gao Y, Qin S, He Z. Fibrotic extracellular vesicles contribute to mechanical ventilation-induced pulmonary fibrosis development by activating lung fibroblasts via JNK signalling pathway: an experimental study. BMJ Open Respir Res. 2023;10(1):e001753.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Jagadapillai R, Rane MJ, Lin X, Roberts AM, Hoyle GW, Cai L, Gozal E. Diabetic Microvascular Disease and Pulmonary Fibrosis: the contribution of platelets and systemic inflammation. Int J Mol Sci. 2016;17(11):1853.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Gressner OA, Gressner AM. Connective tissue growth factor: a fibrogenic master switch in fibrotic liver diseases. Liver Int. 2008;28(8):1065–79.

    Article  CAS  PubMed  Google Scholar 

  30. Wang CM, Hsu CT, Niu HS, Chang CH, Cheng JT, Shieh JM. Lung damage induced by hyperglycemia in diabetic rats: the role of signal transducer and activator of transcription 3 (STAT3). J Diabetes Complications. 2016 Nov-Dec;30(8):1426–33.

  31. Fariña J, Furió V, Fernandez-Aceñero MJ, Muzas MA. Nodular fibrosis of the lung in diabetes mellitus. Virchows Arch. 1995;427(1):61–3.

    Article  PubMed  Google Scholar 

  32. Wolters PJ, Collard HR, Jones KD. Pathogenesis of idiopathic pulmonary fibrosis. Annu Rev Pathol. 2014;9:157–79.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

The establishment of this research was financially supported by the National Project for the Construction of a renowned traditional Chinese medicine experts’ heritage studio (National Chinese Medicine Human Education Letter [2022] No.75), Shandong lung disease heat and poison school inheritance studio construction project (Lu Wei letter [2021] No.45) and the Shandong Province Taishan Scholar Project [grant number tsqn202306392].The study design was independent of the sponsors, as was the collection, analysis, and interpretation of data. The sponsors were not involved in the writing of the report or in the decision to submit the article for publication.

Author information

Authors and Affiliations

Authors

Contributions

XM and YZ designed the study and drafted the manuscript. Data collection and analysis were conducted by FW, XL. WZ provided assistance with methodology guidance and revising the manuscript. The final manuscript was read and approved by all authors.

Corresponding author

Correspondence to Wei Zhang.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethical approval

The summary statistics for the MR study were acquired from GWAS (https://www.ebi.ac.uk). All of this data is ethically approved and freely accessible.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, X., Zhang, Y., Wu, F. et al. Mendelian randomisation highlights type 1 diabetes as a causal determinant of idiopathic pulmonary fibrosis. Diabetol Metab Syndr 16, 90 (2024). https://doi.org/10.1186/s13098-024-01331-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13098-024-01331-x

Keywords