Yanpeng Wang
Yinzhen Zhang
Changwei Zhao2- 1Department of Traditional Chinese Medicine, Changchun University of Chinese Medicine, Changchun, China
- 2Department of Orthopedics, Affiliated Hospital of Changchun University of Chinese Medicine, Changchun, China
- 3Department of Spinal Orthopedics, Weifang Hospital of Traditional Chinese Medicine, Shandong Second Medical University, Weifang, China
Background: The association between depression and musculoskeletal diseases has long been a subject of contentious debate. However, the causal relationship between the two remains uncertain. This study employs a two-sample Mendelian randomization (MR) analysis to investigate the causality between depression and six musculoskeletal diseases.
Methods: In this study, we performed MR analysis to systematically explore the causal relationship between depression and six musculoskeletal disorders. Single nucleotide polymorphisms (SNPs) that are linked to depression were employed as instrumental variables. To ensure robust and reliable conclusions, multiple analytical approaches were utilized, including inverse variance weighting(IVW), weighted median, and MR-Egger regression. Additionally, sensitivity analysis methods such as the MR-Egger intercept test, Cochran’s Q test, leave-one-out analysis, and funnel plot were employed.
Results: Our MR analysis revealed a significant association between depression and cervical spondylosis (depression: OR 1.003, 95% CI 1.002–1.005, P = 8.32E-05; major depressive disorder: OR 1.003, 95% CI 1.001–1.005, P = 0.0052). Furthermore, a strong correlation was noted between major depressive disorder (MDD) and knee osteoarthritis (KOA) (OR 1.299, 95% CI 1.154–1.463, P = 1.50E−5). Sensitivity analysis confirmed the robustness of these findings. Our independent validation study also corroborated these results.
Conclusion: The MR analysis conducted in this study provides evidence supporting a genetic link between depression and cervical spondylosis, as well as KOA. Targeted interventions to manage depression in susceptible populations may contribute to lowering the risk of cervical spondylosis and KOA in these cohorts.
1 Introduction
Depression is a prevalent, debilitating, and costly mental disorder worldwide, affecting more than 300 million individuals. MDD represents a severe form of depression, presenting significant challenges in terms of treatment and profoundly impacting patients’ emotional, physical well-being, daily functioning, learning, and work (1, 2). Previous research has not only demonstrated its adverse effects on cardiovascular and metabolic diseases (3), but also indicated its significant role in specific postural changes in the human body (4), highlighting the close relationship between depression and the occurrence of musculoskeletal diseases (5, 6).
Musculoskeletal diseases (MSD) primarily include chronic pain-related conditions such as cervical spondylosis (CS), KOA, hip osteoarthritis (HOA), rheumatoid arthritis (RA), lumbar disc herniation (LDH), and osteoporosis (OP). These conditions typically manifest as chronic pain and restricted mobility, imposing substantial economic and societal burdens (7), while profoundly impacting individuals’ quality of life (8).Previous observational studies have indicated that CS patients face an elevated risk of depression (9), with similar trends observed in KOA (10), RA (11) and OP patients. However, limited by experimental design and sample size, few studies have explored whether depression increases the risk of MSD. Observational studies are inherently limited by confounding factors and reverse causality, hindering the determination of associations at a genetic level. Furthermore, comprehensive research on the relationship between depression and MSD remains lacking. MR offers a new perspective to address these challenges.
MR is a precise medical analytical method that utilizes genetic variations to establish associations between risk factors and outcomes, thereby evaluating causal effects (12). By capitalizing on the random assignment of genetic variations during conception, MR analysis introduces randomness and mitigates potential confounding factors and reverse causality issues commonly encountered in observational studies. This method effectively integrates genetic information into traditional epidemiological investigations, offering a fresh approach to elucidate the associations between different factors. Therefore, we posit that utilizing the MR approach aids in comprehending the causal relationship between depression and MSD. This study aims to employ MR analysis to investigate the potential association between depression and MSD, thereby offering valuable insights for both preventive and therapeutic interventions.
2 Method and design
2.1 Research design
Figure 1 illustrates the research design and the hypothesis regarding MR utilized in this study. The study employed a two-sample MR analysis and aggregated data from genome-wide association studies (GWAS) to investigate the causal effects of depression on chronic sinusitis CS, KOA, HOA, RA, LDH, and OP. Additionally, independent verification tests were conducted to ensure the reliability of the results. The MR design is guided by three key assumptions: (1) instrumental variables demonstrate a strong correlation with depression; (2) the instrumental variable is not linked to confounding factors; and (3) the association with MSD risk is solely mediated through depression (13).
Figure 1. Basic hypothesis of mendelian randomization.
2.2 Data sources
The depression phenotypes were categorized as depression and MDD. Data on depression included 23,424 patients and 192,220 controls from the Fingen, while data on MDD comprised 170,756 patients and 329,443 controls from the IEU Open GWAS. Furthermore, self-reported datasets on CS with 3,292 patients and 459,641 controls, and LDH with 2,214 patients and 460,796 controls, were obtained through a meta-analysis of UK Biobank sources in GWAS. The datasets on KOA with 24,955 patients and 378,169 controls, HOA with 15,704 patients and 378,169 controls, RA with 14,361 patients and 43,923 controls, and OP with 3,203 patients and 209,575 controls were all sourced from the IEU Open GWAS. In the validation cohort, CS data were extracted from the IEU Open GWAS, encompassing 3,292 cases and 459,641 controls, while KOA data were obtained from the Genetics of Osteoarthritis (GO) (14) Consortium, comprising 62,497 cases and 333,557 controls. All genetic data included in this study were sourced from patients of European descent and were accessible from the open GWAS dataset.1 Since the data used in this study are derived from existing publications and public databases, no additional ethical approval or consent is required for participation.
2.3 Selection of instrumental variables
First, we select SNP at the genome-wide significant level (P < 5 × 10–8) (since only one SNP is associated with CS at the genome-wide significant level, so we set the threshold to the significant level p < 5 × 10–6); then, to remove the linkage disequilibrium, set the genetic distance to 10000kb, and set the linkage disequilibrium parameter r2 threshold = 0.001; third, remove ambiguity and palindromes SNP through the coordination process. Fourth, F statistics are used to evaluate whether the included SNPs is affected by weak tool variables, and the corresponding SNP of F < 10 is eliminated. Finally, we removed SNPs associated with BMI, obesity, lipids, and outcomes through online websites Phenocsnner.
2.4 Statistical analysis
In our study, we applied the IVW method with random effects as the primary analytical approach. Moreover, we employed the Weighted Median and MR Egger methods as supplementary techniques to the IVW method. The IVW method yields precise estimates under the assumption that all SNPs serve as valid instrumental variables (15). To assess estimates under the presence of directional pleiotropy, we utilized the MR Egger and Weighted Median methods, which offer increased robustness (16, 17). These methods were used for comparative analysis in conjunction with the IVW estimates. The MR-PRESSO method was employed in this study to detect any outliers, and any outliers identified were removed for reanalysis. The results were presented as odds ratios (ORs) with their corresponding 95% confidence intervals (CIs), providing an estimate of the associated risk following a one-standard deviation (SD) increase in the risk factors.
We have observed a considerable degree of sample overlap between MDD and the outcome variable. Previous MR analyses (18, 19) have indicated that, in two-sample MR studies, aside from the MR-Egger method, other methodologies can confidently employ the same biological sample repository for analysis purposes. Given our utilization of the IVW method as the primary analytical approach, our findings can to some degree discount the bias stemming from sample overlap.
2.5 Sensitivity testing
The quality control of the results calculated by MR was carried out by sensitivity test. It mainly includes heterogeneity test and pleiotropic test. In addition, in order to prevent bias in the analysis results, this study also conducted a ‘leave-one-out’ and funnel plot sensitivity test.
The analysis in this study was conducted using R4.2.1 software and the Two sample MR package. Following Bonferroni correction (0.05/6), statistical significance was determined at P < 0.008. Values between 0.05 and 0.08 were deemed nominally significant.
3 Results
3.1 Instrumental variables
Following screening, we initially incorporated 30 SNP instrumental variables for depression and 50 SNP instrumental variables for MDD. Subsequently removing repetitive and palindromic sequences (rs12491634, rs2876520, rs4730387, rs35469634), we ultimately opted for 29 and 47 SNPs as instrumental variables to estimate depression and MDD, respectively (Supplementary Table 1).
The instrumental variables utilized in this study accounted for approximately 4.1% and 1.7% of the genetic variance of depression and MDD, respectively. Importantly, the F values associated with the selected SNPs exceeded 10, thereby ruling out the presence of weak instrumental variables (Supplementary Table 2).
3.2 MR analysis of depression and MDD on MSD
In the analysis, it was found that depression is associated with an increased likelihood of CS (OR = 1.003, 95%CI: 1.002–1.005, P-value = 8.32E−05) (refer to Figure 2). Similarly, MDD also elevates the risk of CS, with an OR of 1.003 (95% CI: 1.001–1.005, P-value = 0.0052). Furthermore, the research indicated that MDD is linked to a higher risk of KOA, with an OR of 1.299 (95% CI: 1.154–1.463, P-value = 1.50E−5). Notably, no causal relationship was found between depression, MDD, and HOA,LDH, RA, OP. For detailed results, please refer to Figure 3.
Figure 2. MR analysis and sensitivity analysis results of depression and six musculoskeletal diseases.
Figure 3. MR analysis and sensitivity analysis results of major depressive disorder and six musculoskeletal diseases.
3.3 Sensitivity analysis results
The presence of pleiotropy in our study was assessed through the examination of the MR-Egger regression intercept. The MR Egger intercept values for depression, MDD, and CS were 0.0003 and −0.0001, respectively, indicating values close to 0. Both corresponding p-values, 0.2075 and 0.6167, exceeded 0.05. These results imply the absence of horizontal pleiotropy and multiple effects interference in the MR analysis. Conversely, the MR Egger intercept values for MDD and KOA were 0.015 and 0.183, both exceeding 0.05, suggesting no presence of horizontal pleiotropy (Figures 2, 3).
The heterogeneity tests were conducted to evaluate the variability within the dataset. For depression, MDD, and CS, the Cochran’s Q values and QP values for the IVW method were 6.4495 and 42.1315, respectively, both surpassing the significance threshold of 0.05. These findings indicate the absence of heterogeneity and confirm the reliability of the results. Similarly, in patients with MDD and KOA, the Cochran’s Q and QP values for IVW were 53.658 and 0.107, respectively, exceeding 0.05, suggesting no heterogeneity in the data (refer to Figures 2, 3).
Sensitivity analyses were conducted using the “Leave-one-out” method and funnel plots. These analyses did not reveal any significant abnormalities, as illustrated in Supplementary Figures 1–3.
3.4 Validation cohort
In our re-analysis of the validation cohorts for CS and KOA, consistent trends persist. The findings indicated a continued significant association between depression and CS at the genetic level [depression: (OR: 1.002, 95% CI: 1.001–1.004, p = 0.00003)], MDD [OR: 1.003, 95% CI: 1.001–1.004, p = 0.0047)], along with analogous results observed between MDP and KOA (OR: 1.270, 95% CI: 1.156–1.395, p = 5.7E−07). Furthermore, sensitivity analysis confirmed the robustness of the findings (Figure 4).
Figure 4. MR analysis and sensitivity analysis results of validation cohort.
4 Discussion
In this study, the GWAS dataset was utilized to thoroughly investigate the causal relationship between depression and six musculoskeletal diseases. The findings revealed a genetic-level association between depression and CS, albeit with a small effect size. Furthermore, our research indicated a link between depression and the susceptibility to developing KOA. Conversely, no causal relationship was found between depression, MDD, and other conditions such as HOA, LDH,RA, or OP.
Previous observational studies have reported an elevated risk of depression in patients with CS (20, 21). However, limited research exists on whether depression contributes to an increased risk of CS. A cohort study involving 68,000 individuals in Taiwan, utilizing data from the Taiwan Health Insurance Research Database, demonstrated a 76% higher risk of CS in individuals with depression compared to those without depression; notably, the use of specific antidepressant medications did not mitigate this risk (22). These findings align with our study results. Regarding the relationship between depression and knee pain, well-conducted studies (23) have indicated that depression heightens the risk of knee pain in adults, and antidepressant use is linked to relieving knee pain. Nevertheless, limited observational data exist on whether depression increases the likelihood of KOA. Hence, further validation through extensive cohort studies is essential.
The increased risk of CS and KOA among individuals with depression may be closely linked to poor posture. Research indicates that individuals with depression often exhibit poor posture, which can lead to neck muscle tension and alterations in spinal structure, potentially contributing to the development of CS (4–6). Moreover, individuals with depression often engage in decreased physical activity, leading to muscle atrophy, reduced joint stability, and potentially adverse effects on knee function and outcomes (24).
In exploring the underlying mechanism of the causal relationship between depression and MSD, pro-inflammatory factors are significantly implicated. Prior studies have established a positive connection between inflammatory markers and MDD (25–29), including IL-6, IL-1β, CRP, and TNF. The heightened levels of these inflammatory markers are also prominent drivers of CS by contributing to degeneration and disc deterioration in the cervical region (30). Yin et al. (31) have underscored that these inflammatory factors may incite intervertebral disc degeneration and neuropathic pain, ultimately culminating in the development of CS (32). Additionally, research has delineated a strong relationship between the inflammatory process and the pathophysiology of depression, indicating that inflammation may serve as a critical disease-modifying factor that increases susceptibility to depression (33). Furthermore, evidence suggests that medications prescribed for depression could potentially foster the development of CS. Lin et al. (22) conducted an investigation involving 34,166 patients with depression, segregating them into two groups: one treated with antidepressants and the other without. The outcomes indicated that the incidence of CS in the medication group was 1.8 times higher in comparison to the non-medication group.
The primary strength of our study lies in being the first MR investigation to comprehensively examine the association between depression and MSD. By leveraging the inherent strengths of MR methodology, this study minimizes biases arising from confounding variables and reverse causation. The findings also highlight the impact of depression on the genetic susceptibility to CS and KOA. Nonetheless, several limitations warrant acknowledgment. Firstly, the majority of participants in the MR analysis are of European descent, potentially introducing racial bias concerns. Secondly, due to constraints in the available data, a more thorough subgroup analysis assessing the influence of MDD on MSD outcomes was not feasible. Additionally, it is pertinent to investigate the possible causal relationships between inflammatory mediators and both CS and depression.
5 Conclusion
The MR analysis conducted in our study reveals a genetic link between depression and both CS and KOA. Improving the care and support for populations vulnerable to depression could potentially mitigate the incidence of MSD in these specific cohorts.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
YW: Conceptualization, Writing – original draft, Writing – review and editing. YZ: Conceptualization, Writing – original draft. CZ: Data curation, Investigation, Writing – review and editing. HY: Data curation, Software, Writing – review and editing. CA: Data curation, Software, Writing – review and editing. WZ: Funding acquisition, Supervision, Validation, Visualization, Writing – review and editing. JX: Funding acquisition, Supervision, Validation, Visualization, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Natural Science Foundation of Jilin Province (YDZJ202201ZYTS216), the Shandong Traditional Chinese Medicine Science and Technology Project (Q-2023034), and the 2021 project of Jilin Science and Technology Development Plan (20210303004SF).
Acknowledgments
We thank all those involved 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2024.1398203/full#supplementary-material
Footnotes
References
1. Berton O, Nestler EJ. New approaches to antidepressant drug discovery: Beyond monoamines. Nat Rev Neurosci. (2006) 7:137–51. doi: 10.1038/nrn1846
2. Daly EJ, Trivedi MH, Janik A, Li H, Zhang Y, Li X, et al. Efficacy of esketamine nasal spray plus oral antidepressant treatment for relapse prevention in patients with treatment-resistant depression: A randomized clinical trial. JAMA Psychiatry. (2019) 76:893–903. doi: 10.1001/jamapsychiatry.2019.1189
3. Zhang Y, Chen Y, Ma L. Depression and cardiovascular disease in elderly: Current understanding. J Clin Neurosci. (2018) 47:1–5. doi: 10.1016/j.jocn.2017.09.022
4. Rosario JL, Bezerra Diogenes MS, Mattei R, Leite JR. Differences and similarities in postural alterations caused by sadness and depression. J Bodyw Mov Ther. (2014) 18:540–4. doi: 10.1016/j.jbmt.2013.12.010
5. Canales JZ, Cordas TA, Fiquer JT, Cavalcante AF, Moreno RA. Posture and body image in individuals with major depressive disorder: A controlled study. Braz J Psychiatry. (2010) 32:375–80. doi: 10.1590/s1516-44462010000400010
6. Kumaresan S, Yoganandan N, Pintar FA, Maiman DJ, Goel VK. Contribution of disc degeneration to osteophyte formation in the cervical spine: A biomechanical investigation. J Orthop Res. (2001) 19:977–84. doi: 10.1016/S0736-0266(01)00010-9
7. Beaudart C, Biver E, Bruyere O, Cooper C, Al-Daghri N, Reginster JY, et al. Quality of life assessment in musculo-skeletal health. Aging Clin Exp Res. (2018) 30:413–8. doi: 10.1007/s40520-017-0794-8
8. Carstairs A, Genever P. Stem cell treatment for musculoskeletal disease. Curr Opin Pharmacol. (2014) 16:1–6. doi: 10.1016/j.coph.2014.01.005
9. He J, Xiong W, Li F, Luo W, Gao SC. Depression influences pain and function after cervical disc arthroplasty. J Neurosurg Sci. (2017) 61:39–45. doi: 10.23736/S0390-5616.16.03032-0
10. Li M, Nie Y, Zeng Y, Wu Y, Liu Y, Wu L, et al. The trajectories of depression symptoms and comorbidity in knee osteoarthritis subjects. Clin Rheumatol. (2022) 41:235–43. doi: 10.1007/s10067-021-05847-9
11. Matcham F, Rayner L, Steer S, Hotopf M. The prevalence of depression in rheumatoid arthritis: A systematic review and meta-analysis. Rheumatology. (2013) 52:2136–48. doi: 10.1093/rheumatology/ket169
12. Wu F, Huang Y, Hu J, Shao Z. Mendelian randomization study of inflammatory bowel disease and bone mineral density. BMC Med. (2020) 18:312. doi: 10.1186/s12916-020-01778-5
13. Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: A guide, glossary, and checklist for clinicians. BMJ. (2018) 362:k601. doi: 10.1136/bmj.k601
14. Boer CG, Hatzikotoulas K, Southam L, Stefansdottir L, Zhang Y, Coutinho de Almeida R, et al. Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell. (2021) 184:6003–5. doi: 10.1016/j.cell.2021.11.003
15. Hartwig FP, Davies NM, Hemani G, Davey Smith G. Two-sample Mendelian randomization: Avoiding the downsides of a powerful, widely applicable but potentially fallible technique. Int J Epidemiol. (2016) 45:1717–26. doi: 10.1093/ije/dyx028
16. 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:512–25. doi: 10.1093/ije/dyv080
17. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan NA, Thompson JR. Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: The role of the I2 statistic. Int J Epidemiol. (2016) 45:1961–74. doi: 10.1093/ije/dyw220
18. Minelli C, Del Greco MF, van der Plaat DA, Bowden J, Sheehan NA, Thompson J. The use of two-sample methods for Mendelian randomization analyses on single large datasets. Int J Epidemiol. (2021) 50:1651–9. doi: 10.1093/ije/dyab084
19. Zheng BK, Sun XY, Xian J, Niu PP. Maternal testosterone and offspring birth weight: A Mendelian randomization study. J Clin Endocrinol Metab. (2022) 107:2530–8. doi: 10.1210/clinem/dgac389
20. Stoffman MR, Roberts MS, King JT Jr. Cervical spondylotic myelopathy, depression, and anxiety: A cohort analysis of 89 patients. Neurosurgery. (2005) 57:307–313; discussion 307–313. doi: 10.1227/01.neu.0000166664.19662.43
21. Kelly JC, Groarke PJ, Butler JS, Poynton AR, O’Byrne JM. The natural history and clinical syndromes of degenerative cervical spondylosis. Adv Orthop. (2012) 2012:393642. doi: 10.1155/2012/393642
22. Lin SY, Sung FC, Lin CL, Chou LW, Hsu CY, Kao CH. Association of depression and cervical spondylosis: A nationwide retrospective propensity score-matched cohort study. J Clin Med. (2018) 7:387. doi: 10.3390/jcm7110387
23. Phyomaung PP, Dubowitz J, Cicuttini FM, Fernando S, Wluka AE, Raaijmaakers P, et al. Are depression, anxiety and poor mental health risk factors for knee pain? A systematic review. BMC Musculoskelet Disord. (2014) 15:10. doi: 10.1186/1471-2474-15-10
24. Hurley MV. The role of muscle weakness in the pathogenesis of osteoarthritis. Rheum Dis Clin North Am. (1999) 25:283–298,vi. doi: 10.1016/s0889-857x(05)70068-5
25. Dowlati Y, Herrmann N, Swardfager W, Liu H, Sham L, Reim EK, et al. A meta-analysis of cytokines in major depression. Biol Psychiatry. (2010) 67:446–57. doi: 10.1016/j.biopsych.2009.09.033
26. Figueroa-Hall LK, Paulus MP, Savitz J. Toll-like receptor signaling in depression. Psychoneuroendocrinology. (2020) 121:104843. doi: 10.1016/j.psyneuen.2020.104843
27. Goldsmith DR, Rapaport MH, Miller BJ. A meta-analysis of blood cytokine network alterations in psychiatric patients: Comparisons between schizophrenia, bipolar disorder and depression. Mol Psychiatry. (2016) 21:1696–709. doi: 10.1038/mp.2016.3
28. Haapakoski R, Mathieu J, Ebmeier KP, Alenius H, Kivimaki M. Cumulative meta-analysis of interleukins 6 and 1beta, tumour necrosis factor alpha and C-reactive protein in patients with major depressive disorder. Brain Behav Immun. (2015) 49:206–15. doi: 10.1016/j.bbi.2015.06.001
29. Lamers F, Milanescazhi Y, Smit JH, Schoevers RA, Wittenberg G, Penninx B. Longitudinal association between depression and inflammatory markers: Results from the Netherlands study of depression and anxiety. Biol Psychiatry. (2019) 85:829–37. doi: 10.1016/j.biopsych.2018.12.020
30. De Geer CM. Cytokine involvement in biological inflammation related to degenerative disorders of the intervertebral disk: A narrative review. J Chiropr Med. (2018) 17:54–62. doi: 10.1016/j.jcm.2017.09.003
31. Yin J, Huang Y, Gao G, Nong L, Xu N, Zhou D. Changes and significance of inflammatory cytokines in a rat model of cervical spondylosis. Exp Ther Med. (2018) 15:400–6. doi: 10.3892/etm.2017.5418
32. Wuertz K, Haglund L. Inflammatory mediators in intervertebral disk degeneration and discogenic pain. Global Spine J. (2013) 3:175–84. doi: 10.1055/s-0033-1347299
Keywords: depression, cervical spondylosis, musculoskeletal diseases, Mendelian randomization, GWAS
Citation: Wang Y, Zhang Y, Zhao C, Yang H, Ai C, Zhao W and Xu J (2024) Genetic link between depression and musculoskeletal disorders: insights from Mendelian randomization analysis. Front. Med. 11:1398203. doi: 10.3389/fmed.2024.1398203
Received: 09 March 2024; Accepted: 15 May 2024;
Published: 31 May 2024.
Edited by:
Rubén Queiro, Foundation for Biosanitary Research and Innovation of the Principality of Asturias (FINBA), SpainReviewed by:
Felipe Julio Ramírez, Hospital Clinic of Barcelona, SpainCarlos Montilla, University Hospital of Salamanca, Spain
Copyright © 2024 Wang, Zhang, Zhao, Yang, Ai, Zhao and Xu. 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: Ji Xu, eHVqaTEwNTEzODgxNzlAMTYzLmNvbQ==
†These authors have contributed equally to this work and share first authorship