- 1Department of Interventional Neuroradiology, Sanbo Brain Hospital, Capital Medical University, Beijing, China
- 2Department of Interventional Neuroradiology, Beijing Tiantan Hospital, Beijing Neurosurgical Institute, Capital Medical University, Beijing, China
- 3Department of Neurology, Beijing Pinggu Hospital, Beijing, China
- 4State Key Laboratory of Brain and Cognitive Science, Institute of Biophysics, Chinese Academy of Sciences, Beijing, China
- 5Department of Neurosurgery, Beijing Hospital, Beijing, China
Background: Glioblastoma multiforme (GBM) with synchronous metastasis(SM) is a rare occurrence. We extracted the data of GBM patients from the SEER database to look into the incidence of SM in GBM, determine the prognostic significance of SM in GBM, and assess therapeutic options for patients presenting with SM.
Methods: From 2004 to 2015, information on GBM patients was obtained from the Surveillance, Epidemiology, and End Results (SEER) database. The propensity score matching (PSM) method was employed to mitigate confounding factors between SM and non-SM groups, subsequently investigating the prognostic significance of SM in patients with GBM. Multivariate Cox proportional hazards regression analyses were employed to identify independent prognostic variables for GBM patients with SM. A forest plot was used to visualize the results.
Results: A cohort of 19,708 patients was obtained from the database, among which 272 (1.4%) had SM at the time of diagnosis. Following PSM at a 3:1 ratio, in both univariate and multivariate cox regression analysis, SM (HR = 1.27, 95% CI: 1.09–1.46) was found to be an independent predictive predictor for GBM patients. Furthermore, the Cox proportional hazard forest plot demonstrated that independent risk variables for GBM patients with SM included age (Old vs. Young, HR = 1.44, 95% CI: 1.11–1.88), surgery (biopsy vs. no surgery, HR = 0.67, 95% CI: 0.46–0.96;Subtotal resection vs. no surgery, HR = 0.47, 95% CI: 0.32–0.68;Gross total resection vs. no surgery, HR = 0.44, 95% CI: 0.31–0.62), radiotherapy (HR = 0.58, 95% CI: 0.41–0.83), and chemotherapy (HR = 0.51, 95% CI: 0.36–0.72).
Conclusion: The predictive value of SM in GBM was determined by this propensity-matched analysis using data from the SEER database. Radiotherapy, chemotherapy, and surgery constitute an effective treatment regimen for patients with SM. A more positive approach toward the use of aggressive treatment for GBM patients with SM may be warranted.
Introduction
Glioblastoma multiforme (GBM) is the most common type of primary malignant brain tumor in adults, distinguished by a dismal prognosis and poor quality of life. Intracranial and extracranial metastases have been documented in 1–2% (1) and 0.4–0.5% (2, 3) of cases, occurring late in the course of the disease, respectively. The prevalence rates might be underestimated since metastatic screening is not a standard practice. Synchronous metastasis(SM) of GBM is even more uncommon, with only a few case reports in the literature (4). The precise incidence of SM in GBM remains elusive.
Metastasis holds critical importance for the prognosis of GBM, as evidenced by a comprehensive literature review (5). Several studies have shown the intensive treatment can prolong survival time for certain individuals with metastases (6–8). However, these researches are confined to tiny sample size, and susceptible to selection bias. Therefore, it is essential to uncover additional information regarding prognostic factors and treatment strategies for these patients.
The present study aimed to use large, population-based cancer registry data to investigate both the incidence and risk factors of SM in patients with GBM. Additionally, it seeks to identify prognostic factors and formulate treatment strategies for patients with SM.
Patients and methods
Patients
The patient cases were enrolled in the study based on the latest version of the publicly available SEER17 database (published in November 2022) by using SEER*Stat 8.4.2.1
Patients diagnosed with primary GBM, both with and without SM, between the years 2004 and 2015, were retrospectively identified in this study. The inclusion criteria were as follows: (1) histology diagnosis with GBM (SEER Brain and CNS Recode 1.1.2 GB); (2) year of diagnosis (range: 2004–2015). While the exclusion criteria included: (1) information on SM (CS mets dx) unknown; (2) not only primary tumor; (3) patient with unknown survival time; (4) incomplete or unconfirmed diagnoses, and (5) tumor size unknown. The patient selection flowchart from the SEER database is delineated in Figure 1.
The following variables were incorporated in the current study: demographic characteristics (year, age, sex, race, marital status, household income, rural–urban area); clinicopathological information (primary site, laterality, tumor size, SM, surgery, radiotherapy and chemotherapy), and follow-up data (cause-specific death and survival time).
Age at diagnosis as a continuous variable was separated into young group (<60 years) and old group (≥60 years). Median household income inflation-adjusted to 2021 was used to classify individuals into the low-income ($60,000 or less) and high-income ($60,001 or more) groups. Tumor size was characterized by the dimensions of the primary tumor. The surgical interventions, which pertain to the resection of the primary tumor, were systematically documented using a specific coding schema.: gross total resection (GTR) was designated with codes 30 or 55; subtotal resection (STR) was represented by codes 22, 21, or 40; biopsy was assigned codes 10, 20, or 90; and instances where no surgery was performed were denoted with code 0. SM: Defined as intracranial and/or extracranial metastatic lesions resulting from the spread of a primary GBM at the initial diagnosis. In the SEER database, SM, which includes both extracranial and intracranial metastasis, is denoted by the term “CS mets dx.” Extracranial metastasis is denoted by codes 30 or 50, whereas intracranial metastasis by codes 10 or 20. Cancer-specific survival (CSS) is defined as the duration between the primary diagnosis and the date of death associated with GBM.
Statistical analysis
A Propensity Score Matching (PSM) analysis was utilized to adjust for all other variables between patients with and without SM. The “Matchit” package in R was employed to match the propensity scores between the two groups, employing the nearest neighbor algorithm with a matching ratio of 1:3. Then, Pearson’s Chi-squared test (“gtsummary” package) was used to compare baseline characteristics between groups. Categorical variables are presented as proportions and percentages of the total. Cox proportional hazards regression analysis was performed to evaluate the prognostic effect of SM in GBM (“survival” package). The effect of these confounders was quantified by calculating the changes in the effect size of SM for GBM when each variable is added to the model sequentially in a step-wise fashion (“chest” package) (9). Based on the observed effects, a multivariable Cox proportional hazards regression analysis was conducted. Furthermore, the Kaplan–Meier method and log-rank analysis were utilized to depict the CSS curves before and after PSM in both groups (“survminer” package).
Cox proportional hazards regression models, both univariable and multivariable, were used to identify the prognostic factors for GBM patients with SM. Moreover, a forest plot (“forestmodel” package) was constructed to visualize the findings.
In our study, statistical significance was defined as a two-tailed p value<0.05. The statistical analysis was performed with R program (version 4.3.2).
Results
Baseline characteristics before and after PSM
Based on inclusion and exclusion criteria, a total of 19,708 GBM patients were selected from the SEER database. Demographical and clinical characteristics for GBM with or without SM were shown in Table 1. Patients diagnosed with SM tended to be non-white, infratentorial, and not to have a paired site (p < 0.05); additionally, they were less likely to have undergone radiotherapy, chemotherapy, or surgery (p < 0.001).
Table 1. Characteristics of glioblastoma multiforme patients before and after propensity score matching.
After 1:3 PSM, 1088 patients were matched, which included 272 patients with SM and 816 patients without SM. All demographic and clinical characteristics were well-matched, suggesting that the PSM effectively minimized potential selection bias (p > 0.05).
The prognostic effect of SM in patients after PSM
The univariate Cox proportional hazards model demonstrated a significant increment in mortality risk for GBM patients with SM after PSM (HR: 1.27, 95% CI:1.10–1.46, p < 0.00; Table 2). Moreover, Kaplan–Meier analysis curves comparing CSS for SM and non-SM groups before and after PSM are depicted in Figure 2.
Table 2. Univariate and multivariate analyses of cancer-specific survival (CSS) in the cohort after propensity score matching.
Figure 2. Kaplan–Meier curves of cause-specific survival (CSS) by synchronous metastasis (SM) before (A) and after (B) propensity score matching (PSM).
The confounding effect of each remaining variable on the association of SM on the CSS of GBM patients was quantified and is presented in Figure 3. The results remain stable with each variable being added to multivariable Cox proportional hazards model sequentially in a step-wise fashion. Based on the effects and clinical experience, multivariable Cox proportional hazards model was established. After adjusting for potential confounding variables including age, primary site, tumor size, surgery, chemotherapy and radiotherapy, SM was still found to be an independent risk factor (HR:1.28, 95% CI:1.11–1.48, p < 0.001) in GBM patients (Table 2).
Figure 3. Graphical representation of the change in the estimate of the effect of synchronous metastasis (SM) on cause-specific survival (CSS) with each variable added to the multivariable Cox regression model.
Prognostic factors for GBM patients with SM
In the present research, a cohort of 272 patients diagnosed with GBM exhibiting SM were analyzed to identify and evaluate potential prognostic factors. As delineated in Table 1, the majority were male (56%), white (83%), and diagnosed between 2012 and 2015 (44%). And a total of 191 (70%) patients underwent surgical intervention, 179 (66%) received radiotherapy, and 152 (56%) were administered chemotherapy. Prognostic factors were identified using multivariate Cox regression analyses, revealing that for GBM patients with SM, age (p = 0.007), surgery (p < 0.001), chemotherapy (p < 0.001), and radiotherapy (p = 0.003) were independent prognostic factors (Figure 4). The Kaplan–Meier curves of the age, surgery, radiotherapy and chemotherapy subgroups are displayed in Figure 5.
Figure 4. Forest plot of the multivariate Cox analysis of CSS in GBM patients with SM. NS, No surgery; STR, Subtotal resection; GTR, Gross total resection; SM, Synchronous metastasis.
Figure 5. Kaplan–Meier survival curves comparing the cancer-specific survival (CSS) of glioblastoma multiforme (GBM) patients with synchronous metastasis (SM). (A), Age; (B), Surgery; (C), Radiotherapy; (D), Chemotherapy; NS, No surgery; STR, Subtotal resection; GTR, Gross total resection.
The prognostic effect of different types of SM in patients
Of these, 272 (1.4%) had SM, with 29 (0.1%) presenting extracranial metastasis and 243 (1.3%) presenting intracranial metastasis. Compared to those with extracranial metastasis, patients with intracranial metastasis were more likely to undergo GTR (30% vs. 14%) (Supplementary Table S1). The Kaplan–Meier curve and forest plot generated using the chest package (Supplementary Figure S1) revealed no statistically significant differences between the two groups.
Discussion
GBM represents the most prevalent form of malignant primary brain tumors in adults (10), with a 5-year overall relative survival of 6.9% (11). GBM patient with SM is uncommon (12), and the prognostic factors and optimal therapeutic approaches for these patients have not yet been fully elucidated. In this population-based study, we employed the SEER database to conduct a comprehensive analysis of patients diagnosed with GBM, both with and without SM, between the years 2004 and 2015. SM is extremely rare, affecting 1.4% of all patients with GBM in our cohort. We found that SM was an independent prognostic factor for CSS in patients with GBM before and after PSM. Additionally, the analysis also revealed age, surgery, chemotherapy, and radiation therapy as significant prognostic indicators for patients with SM. Notably, the statistical analysis did not discern a significant difference between intracranial and extracranial metastasis, potentially due to the constrained sample size of patients with SM.
Previous studies (12–14) have demonstrated SM is a significant prognostic factor for GBM, yet the exact value of this prognosis remains unknown. Amelot et al. (15) conducted a retrospective analysis of GBM patients with spinal cord metastasis in a French database between January 2004 and 2016, accompanied by a review of the pertinent literature. They found that spinal cord metastasis is associated with a poor prognosis. In our study, to balance baseline confounding factors, we employed PSM analysis to ascertain the exact value of SM for patients with GBM. After adjusting for various covariates, SM consistently emerged as an independent risk factor in GBM patients. This finding underscores the importance of conducting a comprehensive assessment of synchronous metastasis status at the time of initial diagnosis.
The underlying mechanism accounting for SM remains elusive (16). The majority of GBM metastases occur within the central nervous system, which is believed to be attributed to white matter tract infiltration and cerebrospinal fluid seeding (17–19). Extracranial metastases are documented in approximately 0.4–0.5% of all patients (2, 3, 20), suggesting that such interventions may facilitate the dissemination of tumor cells beyond the central nervous system. The presence of extracranial metastasis may be indicative of the existence of circulating tumor cells (CTCs) within the bloodstream. These cells exhibit advanced characteristics, such as epithelial-to-mesenchymal transition and dormancy (21, 22), which are crucial to survive in the bloodstream. And the detection of glioblastoma CTCs holds significant clinical promise for the early diagnosis and prognosis (23, 24).
The standard treatment of primary GBM is maximal safe resection followed by concomitant radiotherapy and temozolomide chemotherapy and then adjuvant temozolomide (25). However, due to the scarcity of case reports, the optimal treatment strategy for patients with SM has yet to be fully established. A meta-analysis of individual patient data, encompassing 115 younger patients from 1928 to 2013, indicated that while a survival benefit could not be statistically validated, an aggressive treatment approach may be ethically justified whenever feasible (6). Furthermore, another meta-analysis revealed that surgical intervention, chemotherapy, and radiotherapy significantly prolongs survival in GBM patients with metastasis (7). In our study, we discovered that surgical intervention exerted a significant impact on the management of patients with SM. Gross total resection, subtotal resection, and biopsy can significantly enhance the cancer-specific survival of patients compared to those no surgical intervention. This indicated that surgical resection can effectively decrease the tumor mass, consequently restoring neurological functions. Additionally, we identified both radiation therapy and chemotherapy as significant prognostic factors. Consistent with a previous study, patients who underwent an aggressive treatment, including surgery, radiation, chemotherapy, and cerebrospinal fluid shunting, demonstrated favorable prognosis (26). Overall, our results suggested that a combination of surgical intervention, chemotherapy, and radiation constitutes an effective therapeutic strategy for patients with SM. Therefore, GBM patients with SM may experience enhanced benefits from more aggressive therapeutic regimens. Nonetheless, additional studies are required to validate these preliminary findings.
Several limitations of our study need to be recognized. Firstly, in the retrospective study, the restricted sample size of patients with SM (N = 272) might have a risk of selection bias. Therefore, prospective registry studies with larger sample sizes are necessary to validate our findings. Secondly, as our study included patients from 2004 to 2015 and utilized diagnostic criteria established prior to 2016, these were solely based on pathological diagnosis without integrating molecular diagnosis. Thirdly, information was gathered when the diagnosis was initially made from the SEER database, excluding any metastasis that developed subsequently. Fourthly, the SEER database lacks comprehensive details regarding specific chemotherapeutic agents used, including dosages, treatment duration, and patient responses. Details pertaining to radiation therapy and surgical procedures for metastasis are also missing. These limitations may hinder a thorough evaluation of the efficacy of various treatment modalities.
Conclusion
In conclusion, the propensity-matched analysis has identified the prognostic value of SM in patients with GBM based on the SEER database. Our findings suggest that surgery, chemotherapy, and radiotherapy constitute an effective therapeutic strategy for GBM patients with SM. However, further research is necessary to confirm these findings through prospective registry studies with larger sample size.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
HS: Conceptualization, Data curation, Formal analysis, Software, Visualization, Writing – original draft, Writing – review & editing. QM: Conceptualization, Data curation, Formal analysis, Visualization, Writing – original draft, Writing – review & editing. XC: Conceptualization, Visualization, Writing – original draft, Writing – review & editing. YJ: Data curation, Funding acquisition, Writing – original draft, Writing – review & editing. AL: Conceptualization, Project administration, Supervision, Writing – original draft, Writing – review & editing. JL: Conceptualization, Project administration, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Natural Science Foundation of Beijing Municipality (grant no. 7222053), National Natural Science Foundation of China (grant no. 82003134) and Beijing Science and Technology Commission (Project no. Z171100000417046).
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.
The handling editor BH declared a shared parent affiliation with the authors JL, HS, XC, YJ, and AL at the time of review.
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/fneur.2024.1429826/full#supplementary-material
Footnotes
References
1. Choi, PP, and Shapera, S. What’s your call? Drop metastases. CMAJ Can Med Assoc J. (2006) 175:475, 477. doi: 10.1503/cmaj.060308
2. Smith, DR, Hardman, JM, and Earle, KM. Metastasizing neuroectodermal tumors of the central nervous system. J Neurosurg. (1969) 31:50–8. doi: 10.3171/jns.1969.31.1.0050
3. Pasquier, B, Pasquier, D, N’Golet, A, Panh, MH, and Couderc, P. Extraneural metastases of astrocytomas and glioblastomas clinicopathological study of two cases and review of literature. Cancer. (1980) 45:112–25. doi: 10.1002/1097-0142(19800101)45:1<112::AID-CNCR2820450121>3.0.CO;2-9
4. Witoonpanich, P, Bamrungrak, K, Jinawath, A, Wongwaisayawan, S, Phudhichareonrat, S, and Witoonpanich, R. Glioblastoma multiforme at the corpus callosum with spinal leptomeningeal metastasis. Clin Neurol Neurosurg. (2011) 113:407–10. doi: 10.1016/j.clineuro.2010.12.001
5. Yang, R, Isaacs, AM, Cadieux, M, Hirmer, TJ, Creve Coeur, TS, Lapointe, AP, et al. Primary and metastatic glioblastoma of the spine in the pediatric population: a systematic review. Childs Nerv Syst. (2021) 37:1849–58. doi: 10.1007/s00381-021-05098-8
6. Pietschmann, S, Von Bueren, AO, Kerber, MJ, Baumert, BG, Kortmann, RD, and Müller, K. An individual patient data meta-analysis on characteristics, treatments and outcomes of glioblastoma/gliosarcoma patients with metastases outside of the central nervous system. PLoS One. (2015) 10:e0121592. doi: 10.1371/journal.pone.0121592
7. da Cunha, MLV, and Maldaun, MVC. Metastasis from glioblastoma multiforme: a meta-analysis. Rev Assoc Med Bras. (2019) 65:424–33. doi: 10.1590/1806-9282.65.3.424
8. Ray, A, Manjila, S, Hdeib, AM, Radhakrishnan, A, Nock, CJ, Cohen, ML, et al. Extracranial metastasis of gliobastoma: three illustrative cases and current review of the molecular pathology and management strategies. Mol Clin Oncol. (2015) 3:479–86. doi: 10.3892/mco.2015.494
9. Wang, Z. Two postestimation commands for assessing confounding effects in epidemiological studies. Stata J. (2007) 7:183–96. doi: 10.1177/1536867X0700700203
10. Wen, PY, Weller, M, Lee, EQ, Alexander, BM, Barnholtz-Sloan, JS, Barthel, FP, et al. Glioblastoma in adults: a Society for Neuro-Oncology (SNO) and European Society of Neuro-Oncology (EANO) consensus review on current management and future directions. Neuro-Oncology. (2020) 22:1073–113. doi: 10.1093/neuonc/noaa106
11. Ostrom, QT, Price, M, Neff, C, Cioffi, G, Waite, KA, Kruchko, C, et al. CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2016-2020. Neuro-Oncology. (2023) 25:iv1–iv99. doi: 10.1093/neuonc/noad149
12. Sibanda, Z, Farahani, N, Ogbonnaya, E, and Albanese, E. Glioblastoma Multiforme: a rare case of spinal drop metastasis. World Neurosurg. (2020) 144:24–7. doi: 10.1016/j.wneu.2020.08.086
13. Shibahara, I, Saito, R, Osada, Y, Kanamori, M, Sonoda, Y, Kumabe, T, et al. Incidence of initial spinal metastasis in glioblastoma patients and the importance of spinal screening using MRI. J Neuro-Oncol. (2019) 141:337–45. doi: 10.1007/s11060-018-03036-4
14. Dardis, C, Milton, K, Ashby, L, and Shapiro, W. Leptomeningeal metastases in high-grade adult glioma: development, diagnosis, management and outcomes in a series of 34 patients. Front Neurol. (2014) 5 JUN:) 5:220. doi: 10.3389/fneur.2014.00220
15. Amelot, A, Terrier, LM, Cognacq, G, Jecko, V, Marlier, B, Seizeur, R, et al. Natural history of spinal cord metastasis from brain glioblastomas. J Neuro-Oncol. (2023) 162:373–82. doi: 10.1007/s11060-023-04298-3
16. Georgescu, MM, and Olar, A. Genetic and histologic spatiotemporal evolution of recurrent, multifocal, multicentric and metastatic glioblastoma. Acta Neuropathol Commun. (2020) 8:10. doi: 10.1186/s40478-020-0889-x
17. Hu, Y, Yin, W, Ma, J, Liu, J, Huang, S, and Chen, H. Glioblastoma multiforme of the conus medullaris with leptomeningeal dissemination that presented as intracranial hypertension: a case report and literature review. Glioma. (2021) 4:57–61. doi: 10.4103/glioma.glioma_15_21
18. Zhang, K, Yang, Y, Zhuang, J, Guo, G, Chao, X, and Zhang, Z. Intracranial dissemination of glioblastoma multiforme: a case report and literature review. J Int Med Res. (2022) 50:030006052211120. doi: 10.1177/03000605221112047
19. Wright, CH, Wright, J, Onyewadume, L, Raghavan, A, Lapite, I, Casco-Zuleta, A, et al. Diagnosis, treatment, and survival in spinal dissemination of primary intracranial glioblastoma: systematic literature review. J Neurosurg Spine. (2019) 31:723–32. doi: 10.3171/2019.5.SPINE19164
20. Hamilton, JD, Rapp, M, Schneiderhan, TM, Sabel, M, Hayman, A, Scherer, A, et al. Glioblastoma multiforme metastasis outside the CNS: three case reports and possible mechanisms of escape. J Clin Oncol. (2014) 32:e80–4. doi: 10.1200/JCO.2013.48.7546
21. Chistiakov, DA, and Chekhonin, VP. Circulating tumor cells and their advances to promote cancer metastasis and relapse, with focus on glioblastoma multiforme. Exp Mol Pathol. (2018) 105:166–74. doi: 10.1016/j.yexmp.2018.07.007
22. Lombard, A, Goffart, N, and Rogister, B. Glioblastoma circulating cells: reality, trap or illusion? Stem Cells Int. (2015) 2015:1–11. doi: 10.1155/2015/182985
23. Senhaji, N, Houssaini, AS, Lamrabet, S, Louati, S, and Bennis, S. Molecular and circulating biomarkers in patients with glioblastoma. Int J Mol Sci. (2022) 23:474. doi: 10.3390/ijms23137474
24. Müller Bark, J, Kulasinghe, A, Chua, B, Day, BW, and Punyadeera, C. Circulating biomarkers in patients with glioblastoma. Br J Cancer. (2020) 122:295–305. doi: 10.1038/s41416-019-0603-6
25. McKinnon, C, Nandhabalan, M, Murray, SA, and Plaha, P. Glioblastoma: Clinical presentation, diagnosis, and management. The BMJ. (2021) 374:15060. doi: 10.1136/bmj.n1560
Keywords: glioblastoma multiforme, synchronous metastasis, propensity score matching, cancer-specific survival, SEER
Citation: Shen H, Mei Q, Chai X, Jiang Y, Liu A and Liu J (2024) The prognostic significance of synchronous metastasis in glioblastoma multiforme patients: a propensity score-matched analysis using SEER data. Front. Neurol. 15:1429826. doi: 10.3389/fneur.2024.1429826
Edited by:
Bo Yuan Huang, Electric Power Teaching Hospital, Capital Medical University, ChinaReviewed by:
Peng-Bo Zhou, Shengjing Hospital of China Medical University, ChinaMin Bao, China Medical University, China
Liang Zhang, Chinese Academy of Medical Sciences and Peking Union Medical College, China
Theodora Katsila, National Hellenic Research Foundation, Greece
Jinn-Rung Kuo, Chi Mei Medical Center, Taiwan
Gang Lin, Peking University, China
Copyright © 2024 Shen, Mei, Chai, Jiang, Liu and Liu. 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: Jiachun Liu, jiachunliu2005@163.com; Aihua Liu, 15901398688@163.com
†These authors have contributed equally to this work