- 1Zhejiang Chinese Medical University, Hangzhou, Zhejiang, China
- 2Department of Burns and Plastic Surgery, Affiliated Hospital of Jiangsu University, Zhenjiang, Jiangsu, China
- 3Department of Burns and Plastic Surgery, Jiaxing Hospital of Traditional Chinese Medicine, Zhenjiang, Jiaxing, China
The goal was to explore the effect of interleukin-6 (IL6) and C reactive protein (CRP) on malignant melanoma (MM) using two-sample Mendelian randomization.
Methods: Data for this study were obtained from the IEU Open GWAS project website for genome-wide association study data (GWAS) on interleukin-6, C reactive protein levels and malignant melanoma. Inverse variance weighted (IVW) method was mainly used and supplemented with MR-Egger regression and weighted median. Finally, horizontal multivariate validity and heterogeneity tests were performed to assess the stability and reliability of the results.
Results: The results of univariate two-sample MR analyses showed no significant effect of CRP on MM: inverse variance weighting method (OR=0.999, 95% CI: 0.998–1.001, P=0.343), MR-Egger regression (OR= 1.000, 95% CI: 0.998–1.001, P= 0.180), and weighted median method (OR= 0.999, 95% CI: 0.997 to 1.000, P= 0.583), and weighted model (OR= 0.999, 95% CI: 0.998 to 1.001, P= 0.328). Also,IL-6 had no significant effect on MM: inverse variance weighting method (OR= 1.001, 95% CI: 0.999 to 1.002, P=0.461), MR-Egger regression (OR= 1.000, 95% CI: 0.997 to 1.004, P= 0.910), weighted median method (OR= 1.000, 95% CI: 0.998 to 1.002, P= 0.749), and weighted mode (OR= 1.000, 95% CI: 0.998 to 1.002, P= 0.820).
Conclusion: There was no causal relationship between C-reactive protein and IL-6 on the risk of malignant melanoma.
1 Introduction
Malignant melanoma is an aggressive tumor differentiated from melanocytes (1). Epidemiology studies indicated that the incidences of MM are increasing significantly worldwide. The lethality rate reaches high up to 90% (2, 3). It is estimated that 1 in 63 Americans have melanoma (4, 5). MM generally implicates skin tissue, eyes and mucous membranes. Surgery is the main treatment method (2, 6), the margins of resection were determined based on the depth of the tumor. Other treatments such as targeted therapies and immunotherapy have changed the paradigm for advanced melanoma significantly with the development of therapeutic options. However, the prognosis of MM remains less favorable (7). Early diagnosis and timely treatment of MM increases the 5-year overall survival up to 95%. Previously only 5% of MM patients are long-term survivors after metastasis (7, 8). However, in stage IV metastatic melanoma, immunotherapy had excellent results, reaching 50% survival 5 years after treatment with the nivo+ipi combination, but these patients may be at high risk for new exacerbations or drug toxicity while continuing co-administration (9). Therefore, how to effectively diagnose and prevent MM at early stage is crucial to the prognosis of patients.
IL-6 is a pleiotropic cytokine that is almost universally expressed in stromal and immune cells (10). It not only involved in immune response, but also in basic processes such as inflammation, hematopoiesis, bone metabolism and embryonic development (11). The producing process of IL-6 is regulated by a variety of signals in normal cells such as tumor necrosis factor-a (TNF-a), interleukin-1, interferon-b, DNA viruses and RNA viruses (12). IL-6 plays an important role in the pathogenesis and progression of malignant tumors by promoting tumor growth through inhibiting apoptosis and induces tumor angiogenesis (13, 14). It can induce hepatocytes to generate CRP (15). The latter is a member of the pentapeptide protein family secreted by liver in response to a variety of inflammatory cytokines. CRP binds to lysophosphatidylcholine on the surfaces of dead or dying eukaryotic cells and bacteria, and then activates complement via C1q (16). Levels of CRP in serum are markedly increasing in acute inflammatory responses and have little difference in chronic inflammatory responses (17). CRP greater than 10 mg/L is usually associated with infections, certain inflammatory diseases or malignancies (18). A study has shown that high levels of CRP induce an immunosuppressive milieu in melanoma and supports blocking CRP as a therapeutic strategy to augment cancer immune checkpoint therapies (19). Several studies have suggested that IL-6 and CRP may be associated with diagnosis and prognosis of melanoma (12, 17, 20–22). However, it is difficult for traditional observational studies to avoid the interference of potential confounders and reverse causation, which may produce some bias in the results. Therefore, the effects of IL-6 or CRP on MM remain somewhat controversial.
2 Materials and methods
2.1 Research and design
Traditional observational epidemiologic studies face many challenges in discovering disease etiology and inferring causality, such as reverse causation, potential confounders and secondary exposure factors. It is difficult to explore the causal relationship between C-reactive protein and MM in observational and retrospective cohort studies due to the presence of treatment or other confounding factors beyond our control. Mendelian randomization (MR) is a statistical method which uses genetic variants as proxy markers for risk factors for instrumental variable analysis and is based on whole genome sequencing data used to expose causal relationships (23). Thereby, bias in traditional observational studies is circumvented as well as reduced exposure to confounding and reverse causation. Two-sample MR must satisfy the following three core assumptions: (1) IVs is required to be strongly associated with exposure; (2) IVs are not associated with confounders affecting “exposure-endpoints”; (3) IVs affect endpoints through exposure only but not through other pathways as shown in Figure 1.
2.2 Data sources
The data for this study were acquired from the GWAS of IL-6, CRP and MM from the website, gwas.mrcieu.ac.uk, and standardized to the European population database. The sample size of IL-6 database is 21,758 cases and the number of SNPs is 11,782,139; the sample size of CRP database is 353,466 cases and the number of SNPs is 19,057,467; the sample size of MM database is 462,933 cases and the number of SNPs is 9,851,867, as shown in Table 1.
2.3 Selection of SNPs
To satisfy assumption 1-IVs are required closely associated with exposure-P<5 × 10-8 is generally used as the screening condition, SNPs which are significantly associated with exposure were extracted from the GWAS database as the IVs. The linkage disequilibrium is removed (r2 ≤ 0.001; genetic distance<10,000 kilobases) to ensure independence of instrumental variable. Major and minor allele data is extracted for each SNP as well as allele frequencies, β coefficients, standard error (SE) of β coefficients and p-values of relevant associations (24). The F-value of IVs is calculated in the end and the corresponding IVs are included in the study to avoid bias caused by weak IVs if F>10.
In the study of CRP and MM, screening was performed according to the conditions indicated above and a total of 135 SNPs were included for MR analysis of CRP and MM at last. In the study of IL-6 and MM, P<5 × 10-6 as the screening condition was merged with the MM data set due to the small number of SNPs and incorporated eight SNPs for MR analysis of IL-6 and MM, as shown in Table 2.
2.4 MR analysis
In this study, four analytical methods is used to assess the causal relationship between exposure factors and outcomes, including four methods: IVW, MR-Egger regression, weighted median method and weighted mode method (25). IVW was the predominant analytical method. The results of it are the most reliable of all when the heterogeneity and horizontal polytomous effects are absent (26). MR-Egger regression was performed under the InSIDE assumption. Thus, assessing the presence of multinomial using the intercept term is allowed. When the intercept term was 0, it indicated the absence of horizontal multinomial and led to the inference that the results of the MR-Egger regression were consistent with IVW (27). Weighted median method provided a more accurate estimation of causality when the instrumental variable was more than 50% of invalidation. Weighted mode estimation method can detect causality better and with less deviation than MR-Egger regression method When the InSIDE assumption was not satisfied (28).
3 Results
3.1 Results of MR analysis
In the association study between CRP and MM, the results of IVW method indicated that CRP had no significant effect on MM (OR= 0.999, 95% CI: 0.998–1.001, P= 0.343), and MR-Egger regression (OR= 1.000, 95% CI: 0.998–1.001, P= 0.180). The results of the other 2 methods (weighted median method and weighted mode) were on the same trends as the IVW results as shown in Figure 2A.
Figure 2 The Scatter plot of the association between C reactive protein and malignant melanoma. (A) The Scatter plot of the association between interleukin-6 and malignant melanoma (B).
In the association study between IL-6 and MM, the results of IVW method showed that IL-6 had no significant effect on MM (OR= 1.001, 95% CI: 0.999–1.002, P=0.461), and the MR-Egger regression (OR= 1.000, 95% CI: 0.997–1.004, P=0.910). The results of the other 2 methods (weighted median method and weighted mode) were in the same trends as the IVW results as shown in Figure 2B.
3.2 Sensitivity analysis results
The P-value of Cochran Q was used to infer the heterogeneity of MR results in this study (27). If P < 0.05, it indicates that there were no heterogeneity in the included instrumental variables. In the association study between CRP and MM, Cochran Q resulted in the absence of heterogeneity in both IVW analysis (Cochran Q=169.1, P=0.113) and MR-Egger regression analysis (Cochran Q=169.1, P=0.103). The results of the MR regression analysis showed that horizontal multiplicity had no significant influence (intercept=9×10–7, P=0.971).
In the association study between IL-6 and MM, the results of Cochran Q for IVW analysis (Cochran Q=2.866, P=0.897) and MR-Egger regression analysis (Cochran Q=2.822, P=0.831) were not heterogeneous. The results of the MR regression analysis showed that horizontal pleiotropy had no significant effect on the results of MR analysis (intercept=3.6×10–5, P=0.842) as well. The results of “leave-one-out” method suggested that excluding SNPs separately did not significantly affect the causal associations which indicated the reliability of the MR results as shown in Figure 3. The symmetry between the two sides of the distribution of the funnel plot also indicated that the causal associations were less affected by the potential bias as shown in Figure 4.
Figure 3 “Leave-one-out” sensitivity analysis of the causal effect of CRP and melanoma (left). “Leave-one-out” sensitivity analysis of the causal effect of IL-6 and melanoma (right).
Figure 4 A funnel plot of the association between C reactive protein and malignant melanoma (left). A funnel plot of the association between interleukin-6 and malignant melanoma (right).
4 Discussion
Exploring the causal relationship between IL-6 or CRP and MM is vital for the early diagnosis and prognosis of the disease. No current study has demonstrated the underlying mechanisms clearly between IL-6 or CRP and MM.
The relationship between IL-6 and MM has been discussed for decades. However, the underlying mechanisms remain unclear. A link between inflammation and cancer has long been suspected and it has been hypothesized that at least 15% of cancers are developed from chronic inflammation or infection (29). In clinical practice, physicians always seek a balance between protecting their own immunity and ensuring a good tumor response, prolonged use of immunosuppressive agents or use at ICI initiation may impede antitumor immune response, but Inhibitors of IL-6 can synergize with ICIs on tumor response and also prevent severe immune-related adverse events (irAEs) (30). In an animal study, IL-6 was found to upregulate the expression of CCR5 and arginase 1 in MDSC through a STAT3-dependent mechanism. MDSC differentiated in the presence of IL-6 strongly suppressed the function of CD8+ T cells and IL-6 overexpressing tumors grew at a significantly slower rate in CD8+ T cell-activated mice. Meanwhile, a correlation between IL-6 levels, phosphorylated STAT3 and CCR5 expression in tumor-infiltrating MDSC was confirmed in a RET transgenic melanoma mouse model (31). It has been shown that interleukin-6 (IL-6) inhibits the growth of melanocytes and early-stage melanoma cells but has no significant inhibitory effect on late-stage melanoma cells. However, the results of this MR analysis do not support a causal relationship between IL-6 and MM.
CRP has been used as an independent prognostic indicator in studies of a variety of diseases, including multiple myeloma, lymphoma, melanoma and ovarian, renal, pancreatic and gastrointestinal tumors (32). One study noted the importance of CRP as a predictor of survival in melanoma (33). In a study of 1144 melanoma patients (587 initial and 557 confirmed), CRP was found to be an independent prognostic marker in melanoma patients. The measurements of CRP should be considered for inclusion in prospective studies of prognosis in melanoma patients and clinical trials of systemic therapy for melanoma patients (21). A prospective study has shown that CRP can be used as a biomarker to assist immune checkpoint inhibitors in treatment of immune-related adverse events in melanoma patients (34). Several studies have shown that subsequent CRP levels are usually decreased in cancer patients with solid tumors and treated infections (35). Patients with advanced cancer are susceptible to a variety of infections that can negatively impact prognosis and increase mortality.
IL-6 and CRP are the two most commonly used biomarkers. Large-scale epidemiologic studies and meta-analyses have demonstrated a strong association between IL-6 or CRP and mortality outcomes from a variety of causes (including diseases such as cancer, cardiovascular disease and depression) (36–38). However, the potential mechanisms of IL-6 or CRP on MM remain to be elucidated. In past large randomized trials, higher IL-6 and CRP were associated with shorter overall survival in patients with metastatic melanoma who received ICI or chemotherapy (21). Low levels of IL6 and CRP may indicate a less inflammatory tumor microenvironment, which can enhance the effectiveness of treatment with ICIs. Monitoring these molecules in cancer patients undergoing immunotherapy may help identify those who are more likely to benefit from the treatment.
The aim of the this study was to investigate the causal relationship between IL-6 or CRP and MM using two-sample Mendelian randomization. The research showed that there was no significant causal relationship between IL-6 or CRP and MM. Besides, the four different methods of MR analysis came to consistent results.
The advantage of this study is that the selected GWAS sample size was large and two-sample MR analysis was performed, which can effectively avoid bias and reverse causality interference from observational studies. The reliability of the results was assessed by four methods: IVW, MR-Egger regression, weighted median method and weighted mode method. There are some limitations in this study without doubt. Firstly, all GWAS data were derived from European populations and the applicability to other populations still requires further research. Secondly, this study used the statistical results as a conclusion which does not apply for further in-depth research.
In conclusion, this study investigated the causal relationship between IL-6 or CRP and MM. The results showed that there was no significant causal relationship between IL-6 or CRP and MM.
Author contributions
QW: Writing – original draft, Writing – review & editing. WZ: Data curation, Writing – original draft, Writing – review & editing. SP: Funding acquisition, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Zhejiang Province Medicine and Health Science and Technology Program (2021KY1124) and Jiaxing Public Welfare Research Program Projects (2021AD30019).
Conflict of interest
The authors declare 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. Lodde G, Zimmer L, Livingstone E, SChadendorf D, Ugurel S. Malignant melanoma. Hautarzt. (2020) 71:63–77. doi: 10.1007/s00105–019-04514–0
2. Lens MB, Dawes M. Global perspectives of contemporary epidemiological trends of cutaneous Malignant melanoma. Br J Dermatol. (2004) 150:179–85. doi: 10.1111/j.1365-2133.2004.05708.x
3. Aubuchon MM, Bolt LJ, Janssen-Heijnen ML, Verleisdonk-Bolhaar ST, van Marion A, van Berlo CL. Epidemiology, management and survival outcomes of primary cutaneous melanoma: a ten-year overview. Acta Chir Belg. (2016) 117:29–35. doi: 10.1080/00015458.2016.1242214
4. Markovic SN, Erickson LA, Rao RD, Weenig RH, Pockaj BA, Bardia A, et al. Malignant melanoma in the 21st century, part 1: epidemiology, risk factors, screening, prevention, and diagnosis. Mayo Clin Proc. (2007) 82:364–80. doi: 10.4065/82.3.364
5. Dzwierzynski WW. Melanoma risk factors and prevention. Clin Plast Surg. (2021) 48:543–50. doi: 10.1016/j.cps.2021.05.001
6. Kyrgidis A, Tzellos T, Mocellin S, Apalla Z, Lallas A, Pilati P, et al. Sentinel lymph node biopsy followed by lymph node dissection for localised primary cutaneous melanoma. Cochrane Database Syst Rev. (2015) 5):CD010307. doi: 10.1002/14651858.CD010307.pub2
7. Hartman RI, Lin JY. Cutaneous melanoma-A review in detection, staging, and management. Hematol Oncol Clin North Am. (2019) 33:25–38. doi: 10.1016/j.hoc.2018.09.005
8. Cummins DL, Cummins JM, Pantle H, Silverman MA, Leonard AL, Chanmugam A. Cutaneous Malignant melanoma. Mayo Clin Proc. (2006) 81:500–7. doi: 10.4065/81.4.500
9. Postow MA, Goldman DA, Shoushtari AN, Betof Warner A, Callahan MK, Momtaz P, et al. Adaptive dosing of nivolumab + Ipilimumab immunotherapy based upon early, interim radiographic assessment in advanced melanoma (The ADAPT-IT study). J Clin Oncol. (2021) 40:1059–67. doi: 10.1200/JCO.21.01570
10. Montesarchio V, Parrela R, Iommelli C, Bianco A, Manzillo E, Fraganza F, et al. Outcomes and biomarker analyses among patients with COVID-19 treated with interleukin 6 (IL-6) receptor antagonist sarilumab at a single institution in Italy. J Immunother Cancer. (2020) 8:2. doi: 10.1136/jitc-2020–001089
11. Hirano T. IL-6 in inflammation, autoimmunity and cancer. Int Immunol. (2021) 33:127–48. doi: 10.1093/intimm/dxaa078
12. Hoejberg L, Bastholt L, Schmidt H. Interleukin-6 and melanoma. Melanoma Res. (2012) 22:327–33. doi: 10.1097/CMR.0b013e3283543d72
13. Naugler WE, Karin M. The wolf in sheep’s clothing: the role of interleukin-6 in immunity, inflammation and cancer. Trends Mol Med. (2008) 14:109–19. doi: 10.1016/j.molmed.2007.12.007
14. Aggarwal BB, Sethi G, Ahn KS, Sandur SK, Pandey MK, Kunnumakkara AB, et al. Targeting signal-transducer-and-activator-of-transcription-3 for prevention and therapy of cancer: modern target but ancient solution. Ann N Y Acad Sci. (2006) 1091:151–69. doi: 10.1196/annals.1378.063
15. Weinhold B, Rüther U. Interleukin-6-dependent and -independent regulation of the human C-reactive protein gene. Biochem J. (1997) 327:425–9. doi: 10.1042/bj3270425
16. Du Clos TW. Function of C-reactive protein. Ann Med. (2000) 32:274–8. doi: 10.3109/07853890009011772
17. Schildbach VA, Horn S, Hidalgo-Gadea G, Johannis W, Mauch C, Franklin C. C-reactive protein and lymphocyte-to-monocyte ratio predict recurrence in stage III melanoma patients with microscopic sentinel lymph node metastasis. Cancers (Basel). (2023) 15:2. doi: 10.3390/cancers15030702
18. Li Z, Feng J, Sun X. Is C-reactive protein a specific marker in melanoma. J Clin Oncol. (2015) 33:3673–4. doi: 10.1200/JCO.2015.62.2696
19. Yoshida T, Ichikawa J, Giuroiu I, Laino AS, Hao Y, Krogsgaard M, et al. C reactive protein impairs adaptive immunity in immune cells of patients with melanoma. J Immunother Cancer. (2020) 8:12. doi: 10.1136/jitc-2019–000234
20. Findeisen P, Zapatka M, Peccerella T, Matzk H, Neumaier M, Schadendorf D, et al. Serum amyloid A as a prognostic marker in melanoma identified by proteomic profiling. J Clin Oncol. (2009) 27:2199–208. doi: 10.1200/JCO.2008.18.0554
21. Fang S, Wang Y, Sui D, Liu H, Ross MI, Gershenwald JE, et al. C-reactive protein as a marker of melanoma progression. J Clin Oncol. (2015) 33:1389–96. doi: 10.1200/JCO.2014.58.0209
22. Laino AS, Woods D, Vassallo M, Qian X, Tang H, Wind-Rotolo M, et al. Serum interleukin-6 and C-reactive protein are associated with survival in melanoma patients receiving immune checkpoint inhibition. J Immunother Cancer. (2020) 8:8. doi: 10.1136/jitc-2020–000842
23. Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. (2008) 27:1133–63. doi: 10.1002/sim.3034
24. Dusingize JC, Olsen CM, An J, Pandeya N, Law MH, Thompson BS, et al. Body mass index and height and risk of cutaneous melanoma: Mendelian randomization analyses. Int J Epidemiol. (2020) 49:1236–45. doi: 10.1093/ije/dyaa009
25. Zheng J, Baird D, Borges M, Bowden J, Hemani G, Haycock P, et al. Recent developments in mendelian randomization studies. Curr Epidemiol Rep. (2017) 4:330–45. doi: 10.1007/s40471–017-0128–6
26. Burgess S, Dudbridge F, Thompson SG. Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat Med. (2015) 35:1880–906. doi: 10.1002/sim.6835
27. 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
28. 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:1985–98. doi: 10.1093/ije/dyx102
29. Kuper H, Adami HO, Trichopoulos D. Infections as a major preventable cause of human cancer. J Intern Med. (2000) 248:171–83. doi: 10.1046/j.1365-2796.2000.00742.x
30. Tison A, Garaud S, Chiche L, Cornec D, Kostine M. Immune-checkpoint inhibitor use in patients with cancer and pre-existing autoimmune diseases. Nat Rev Rheumatol. (2022) 18:641–56. doi: 10.1038/s41584–022-00841–0
31. Weber R, Riester Z, Hüser L, Sticht C, Siebenmorgen A, Groth C, et al. IL-6 regulates CCR5 expression and immunosuppressive capacity of MDSC in murine melanoma. J Immunother Cancer. (2020) 8:11. doi: 10.1136/jitc-2020–000949
32. Mahmoud FA, Rivera NI. The role of C-reactive protein as a prognostic indicator in advanced cancer. Curr Oncol Rep. (2002) 4:250–5. doi: 10.1007/s11912–002-0023–1
33. Tartour E, Dorval T, Mosseri V, Deneux L, Mathiot C, Brailly H, et al. Serum interleukin 6 and C-reactive protein levels correlate with resistance to IL-2 therapy and poor survival in melanoma patients. Br J Cancer. (1994) 69:911–3. doi: 10.1038/bjc.1994.175
34. Kallio R, Bloigu A, Surcel HM, Syrjälä H. C-reactive protein and erythrocyte sedimentation rate in differential diagnosis between infections and neoplastic fever in patients with solid tumours and lymphomas. Support Care Cancer. (2001) 9:124–8. doi: 10.1007/s005200000181
35. Lauwyck J, Beckwée A, Santens A, Schwarze JK, Awada G, Vandersleyen V, et al. C-reactive protein as a biomarker for immune-related adverse events in melanoma patients treated with immune checkpoint inhibitors in the adjuvant setting. Melanoma Res. (2021) 31:371–7. doi: 10.1097/CMR.0000000000000748
36. Schnabel RB, Yin X, Larson MG, Yamamoto JF, Fontes JD, Kathiresan S, et al. Multiple inflammatory biomarkers in relation to cardiovascular events and mortality in the community. Arterioscler Thromb Vasc Biol. (2013) 33:1728–33. doi: 10.1161/ATVBAHA.112.301174
37. Singh-Manoux A, Shipley MJ, Bell JA, Canonico M, Elbaz A, Kivimäki M. Association between inflammatory biomarkers and all-cause, cardiovascular and cancer-related mortality. CMAJ. (2016) 189:E384–90. doi: 10.1503/cmaj.160313
Keywords: Mendelian randomization, malignant melanoma, C-reactive protein, interleukin-6, causal association
Citation: Wang QJ, Zheng W and Pan SF (2024) Exploring the causal relationship between interleukin-6 or C reactive protein and malignant melanoma using a two-sample Mendelian randomization approach. Front. Oncol. 14:1375362. doi: 10.3389/fonc.2024.1375362
Received: 23 January 2024; Accepted: 06 June 2024;
Published: 17 June 2024.
Edited by:
Ari VanderWalde, Caris Life Sciences Inc., United StatesReviewed by:
Domenico Mallardo, G. Pascale National Cancer Institute Foundation (IRCCS), ItalyHui Lu, Zhejiang University, China
Copyright © 2024 Wang, Zheng and Pan. 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: Sun Feng Pan, pansunfeng@163.com
†These authors have contributed equally to this work