- Department of Neurosurgery, Affiliated Hangzhou Chest Hospital, Zhejiang University School of Medicine, Hangzhou, China
Purpose: This study was launched to ascertain the independent prognostic factors influencing the overall survival (OS) prognosis of intracranial subependymoma and construct a prognostic model to predict OS time.
Materials and methods: We collected data from patients with intracranial subependymoma, including treatment data, follow-up data, and clinical and pathological characteristics from the SEER database within 2004 to 2016, and patients were randomly classified into training and validation cohorts. Univariate and multivariate analyses were applied to the training group through building a Cox proportional hazards model. According to the results of multivariate analysis, we established a nomogram to forecast the OS rate of the per-case patient graphically, then calculated the accuracy of verification in both training and validation cohorts by concordance index (C-index). Univariate and multivariate analyses were used for different subgroups of unoperated versus operated, gross total resection (GTR), subtotal resection (STR), and biopsy after using the propensity score matching (PSM) analyses.
Results: A total of 667 patients were enrolled, and we randomly assigned 535 patients (80.21%) into the training cohort and 132 patients (19.79%) into the validation cohort. Age [hazard ratio (HR) = 6.355; 95% confidence interval (CI), 2.240–18.029; p = 0.001] and sex (HR = 0.475; 95% CI, 0.232–0.974; p = 0.042) were the independent prognostic factors in the training cohort. On the basis of age and sex, the nomogram was established to predict the OS for every patient (C-index = 0.733 ± 0.065 in the training cohort and 0.850 ± 0.065 in the validation cohort), and calibration plots reflected the reliability of the nomogram. Age, gender, or laterality was the independent prognostic factor for OS in the different matched subgroups of unoperated versus operated, GTR, STR, and biopsy. Surgical treatment, race, year of diagnosis, insurance, tumor location, tumor size, pathology, tumor grade, and radiation were not statistically significantly different in OS for subependymoma in our research.
Conclusion: Age and sex were the independent prognostic variables for OS in intracranial subependymoma. According to our research, we should not be more inclined to choose conservative or surgical treatment. Nonetheless, the information that we present might be useful to suggest potential hypotheses to be tested in the clinical research setting.
Introduction
Subependymoma is a neoplasm with a low incidence and low degree of malignancy (1–3).
Middle-aged and elderly men were the most affected age group by this type of cancer (4). The location of neoplasms was more likely to occur in the ventricle system than in the brain parenchyma or spinal cord (5). Generally speaking, surgical intervention has been recommended once symptoms occur, such as hydrocephalus (6), and conservative treatment has been used for incidental asymptomatic subependymoma. However, there is no detailed analysis of the different prognosis between conservative and surgical treatment including biopsy, STR, and GTR. Consequently, our study was launched to ascertain the independent factors influencing the OS prognosis of intracranial subependymoma and construct a prognostic model to predict OS time through exploring the SEER database.
Methods
Data
The data of 667 patients with intracranial subependymoma were investigated, including treatment and follow-up data, and clinical and pathological characteristics between 2004 and 2016 from the SEER (1975–2016 varying) database, by the SEER*Stat software (version 8.3.9.2).
This study’s inclusion criteria included the following: (1) patient’s ICD-O-3 histology codes in accordance with 9383/0 (subependymoma, benign), 9383/1 (subependymoma), or 9383/2(subependymoma, malignant); and (2) patients with definite information on the vital status and OS.
This study’s exclusion criteria included the following: (1) patients with no specific survival time or with an OS time of less than 1 month; (2) tumor location involving pineal gland (C75.3), spinal cord (C72.0), or optic nerve (C72.3); and (3) the patient had no other specific information or unknown treatment, only a death certificate or an autopsy. The data on age, sex, race, year of diagnosis, insurance, marital status, primary site, tumor size, pathology, grade, laterality, primary site surgery, radiation, vital status, and OS were obtained. The method of retrieving data from the database is shown in Figure 1.
Endpoints
We used OS defined from the time of diagnosis to death or last investigation as the primary endpoint.
Statistical analysis
The whole sample was divided into a training and a validation cohort. Age, a continuous variable, was changed to an ordered classification variable. Disordered classification variables were analyzed by using the χ2 test or Fisher’s exact test, including sex, race, year of diagnosis, marital status, primary site, tumor size, grade, laterality, primary site surgery, and radiation. Ordered classification variables were analyzed by using the Mann-Whitney U test, including age, insurance, and pathology. Different survival rates of variables were graphically evaluated by using the Kaplan-Meier method.
The Cox proportional hazards model was used to perform univariate and multivariate analyses on the training group. According to clinical independent prognostic factors, a nomogram predicting survival probabilities at 3, 5, and 10 years for subependymoma patients was constructed through using the rms package in R (version 4.1.2) in the training group. The model’s C-index, and 3 - and 5-year calibration curves in the training cohort were calculated. The nomogram was further validated by calculating C-index in the validation cohort.
The Cox proportional hazards model method and PSM were used in different subgroups of 667 patients including unoperated versus operated, GTR, STR, and biopsy. A logistic regression model was constructed with operation status as the dependent variable for calculating the propensity scores. One-to-one matching without replacement was performed using the nearest-neighbor match on the logit of the propensity score for confounding factors (derived from age, sex, race, marital status, primary site, tumor size, pathology, grade, laterality, and radiation). The χ2 test, Fisher’s exact test, and Mann-Whitney U test were used to inspect the statistical differences of subgroups before and after matching. Cox proportional hazards model was used to perform univariate and multivariate analyses of various subgroups’ data after PSM. Equilibrium of covariables between subgroups was indicated by p > 0.05.
Various statistical methods were finished in this paper by SPSS (SPSS 26.0, IBM, Inc., Armonk, NY, United States) and R software (R 4.1.2, Vienna, Austria). The p-value < 0.10 of the factor in the univariate analysis was included in the multivariate analysis. Two-tailed p-value < 0.05 was indicated statistically significant (7).
Results
Patient characteristics
Our study included 667 cases of intracranial subependymoma, which randomly assigned 535 patients (80.21%) into the training cohort and 132 patients (19.79%) into the validation cohort (Figure 1). The median OS for all of patients, training cohort, and validation cohort was 56 months [interquartile range (IQR), 24–93], 56 months (IQR, 22–93), and 57 months (IQR,25–90), respectively (Table 1).
The survival factors of the training cohort
The survival curves of age (p < 0.0001; Figure 2A) and sex (p = 0.038; Figure 2B) were compared using a log-rank test. The Cox proportional hazards model was used to perform univariate and multivariate analysis for the training cohort. As exhibited in Table 2, age (p < 0.001) and sex were (p = 0.042) independent prognostic predictors. Patients who were male or more than 60 years of age had less OS time compared with patients who were female or less than 60 years of age. Race, year of diagnosis, insurance, marital status, primary site, tumor size, pathology, grade, laterality, primary site surgery, and radiation had no statistically significant differences in OS for subependymoma in our research (Table 2).
Figure 2 Prognosis of intracranial subependymoma for OS in the training cohort. (A) OS between the different age groups. (B) OS between the different sex groups. OS, overall survival.
Construction and validation of the nomogram
The Cox proportional hazards model uncovered two significant factors that were used to build the nomogram in the training cohort at last (Figure 3A). As exhibited in Figures 3B, C, calibration diagrams complement the internal validation of training queue. The C-index of the training cohort was 0.733 ± 0.065. The C-index and calibration plots confirmed the dependability of the nomograms. Then, the C-index of the validation cohort was 0.850 ± 0.065. Therefore, the 3-year, 5-year, and 10-year predictions of OS by the nomograms were reliable.
Figure 3 Nomogram analyses for patients with intracranial subependymoma. (A) A nomogram for predicting 3-, 5- and 10-year OS of patients. (B) Calibration curve of the nomogram predicting 3-year OS in training cohort. (C) Calibration curve of the nomogram predicting 5-year OS in training cohort. OS, overall survival.
Different subgroups after COX regression analysis and PSM
We performed different subgroup analyses to determine whether surgery was an independent predictor for OS. After PSM of 667 patients in unoperated versus operated, 210 non-operative patients were matched with 210 surgical patients (Tables 3, S1). In the matched cohort, there was no significant difference in OS between the non-surgical and surgical groups (HR = 0.788; 95% CI, 0.457–1.359; p = 0.391; Table 3). In the multivariable regression analysis, age (HR = 8.870; 95% CI, 2.106–22.410; p = 0.001; Table 3) and sex (HR = 0.380; 95% CI, 0.170–0.846; p = 0.018; Table 3) were independent risk prognostic factors for OS.
Table 3 The characteristics of 420 patients from the 667 patients grouped according to no surgery and surgery after PSM.
After PSM of 470 patients in unoperated versus GTR, 164 non-operative patients were matched with 164 GTR patients (Tables 4, S2). In the matched cohort, there was no significant difference in OS between the non-surgical and GTR groups (HR = 0.562; 95% CI, 0.299–1.054; p = 0.072; Table 4). In the multivariable regression analysis, sex (HR = 0.211; 95% CI, 0.065–0.684; p = 0.010; Table 4) was an independent risk prognostic factor for OS.
Table 4 The characteristics of 328 patients from the 470 patients grouped according to no surgery and GTR after PSM.
After PSM of 339 patients in unoperated versus STR, 92 non-operative patients were matched with 92 STR patients (Tables 5, S3). In the matched cohort, there was no significant difference in OS between the non-surgical and STR groups (HR = 0.765; 95% CI, 0.330–1.772; p = 0.532; Table 5). In the multivariable regression analysis, laterality (HR = 0.300; 95% CI, 0.106–0.847; p = 0.023; Table 5) was an independent risk prognostic factor for OS.
Table 5 The characteristics of 184 patients from the 339 patients grouped according to no surgery and STR after PSM.
After PSM of 344 patients in unoperated versus surgery NOS or excisional biopsy, 87 non-operative patients were matched with 87 patients with surgery NOS or excisional biopsy (Tables 6, S4). In the matched cohort, there was no significant difference in OS between the non-surgical and surgery NOS or excisional biopsy groups (HR = 0.596; 95% CI, 0.258–1.377; p = 0.225; Table 6). In the multivariable regression analysis, age (HR = 10.758; 95% CI,2.377–48.693; p = 0.002; Table 6) was an independent risk prognostic factor for OS.
Table 6 The characteristics of 174 patients from the 344 patients grouped according to no surgery and surgery NOS or excisional biopsy after PSM.
Discussion
Scheinker reported a case of a newly recognized tumor derived from the fourth subependymal zone in a 56-year-old man and firstly named subependymoma (8). To date, subependymoma was sporadically reported on case reports (4, 9–11) and accounted for 0.07%-0.7% of all brain tumors (9, 12). Subependymomas were brain neoplasms that tended to be benign, to be less aggressive, to grow slowly, and to be histologically classified as World Health Organization (WHO) grade 1 (13). D’Amico et al. reported a case that was diagnosed with subependymoma by pathological biopsy; CT and MRI confirmed no significant tumor progression after a 36-year follow-up, highlighting the extremely indolent nature of subependymoma (14). The pathogenesis of subependymoma may be related to potential precursor cells (13, 15). Zhiyong et al. reported that 43 patients with subependymoma were found in 60,000 cases of surgically intracranial tumors and the incidence of intracranial subependymoma was about 0.07%. The lesions were mostly located in lateral ventricles accounting for 65% of cases, followed by the fourth ventricle and third ventricle accounting for 19% and 7% of cases, respectively. Tumors were less common in the brain parenchyma and stem (2, 13). The occurrence of symptoms, such as initial clinical manifestations of increased intracranial pressure, was related to the disturbance of cerebrospinal fluid circulation caused by the tumor. Uncommon clinical symptoms including epilepsy, memory loss, ataxia, tremor, blurred vision, and subarachnoid hemorrhage have been reported in some cases (2, 7).
However, there is no prediction model for the OS of subependymoma and no large sample study about the impact of different surgical methods on patient prognosis.
We conducted a study for subependymoma based on the SEER database. The National Cancer Institute’s SEER database collected cancer diagnosis, treatment, and survival data for approximately 30% of the U.S. population. SEER database is an important population-based resource that has become a unique research resource for oncology practice in the United States. The SEER database had the following advantages: representative and universal responses to disease in the United States population, long data collection time, large number of cases, and collection of specific cancer outcomes.
However, the SEER database had the following limitations: individual-level data on specific cancer risks and treatments are incomplete. The accuracy and completeness of raw data collected from the registry needed to be improved. SEER database could not evaluate the progression-free survival (PFS) of tumors.
In this study, we combined the treatment data, follow-up data, and clinical and pathological data from 535 patients in the training group to construct a nomogram for the prediction of the OS of each patient.
Independent prognostic predictors and nomogram
Nguyen et al. drew a conclusion that age < 40 years, female sex, and location within ventricles or near brain stem were positive factors with OS by analyzing 466 cases of intracranial subependymomas from 2004 to 2013 in the SEER database (16). The authors suggest that surgery remains a mainstay treatment. Like prior studies, our study supported that age and sex were significant independent predictors of OS. After our statistical analysis, the prognostic model constructed by age and sex was in good agreement with the reality.
However, whether in the training or subgroup cohort, we revealed that surgery, tumor size, and location were not independent prognostic factors for OS. This seems to challenge the choice of surgical treatment.
D’Amico et al. found that the presence of early malignant lesions in subependymoma cannot be confirmed by early imaging examination and drew a conclusion that early resection was preferred by immunohistochemical analysis of 31 patients with pathologically proven subependymomas (9).
However, some scholars also proposed that conservative treatment was the main treatment for subependymoma. Kammer et al. reviewed 33 cases and showed that subependymomas were usually symptomless; 29 patients were discovered by chance. Subependymoma with no obvious growth tendency seldom led to decompensation of cerebrospinal fluid circulation by blocking the interventricular foramen or Magendie foramen.
In other words, hydrocephalus was relatively rare in subependymoma, which recommended expectant treatment or longer imaging follow-up than other lesions at the same location (10). With a retrospective analysis of 13 patients with intracranial WHO grade 1 subependymoma from 1990 to 2015, Varma revealed that occasional intraventricular subependymoma could be treated conservatively with MRI monitoring. Because there was no significant change in disease during a mean follow-up of 46 months, long-term follow-up was not necessary (14). The author further expounded that hydrocephalus was the main complication of surgical treatment of hydrocephalic subependymoma (2, 14). The appeal suggests that conservative treatment was also an appropriate approach. This may seem counterintuitive. Due to the low degree of malignancy, there were few reports of death caused by subependymoma in a short period. Patients are more likely to die from accidents or other factors.
We recognized that the established nomogram had some value in evaluating patient prognosis and were inclined to use models to predict the prognosis of the subependymoma.
Although nomograms had certain predictive accuracy in the training and validation groups in our study, the treatment strategy still needed to be further improved through subsequent studies, considering the inherent limitations of the SEER databases.
This may not mean that surgery was meaningless for subependymoma. We conducted a further subgroup analysis of the benefits of surgical treatment.
Subgroups analyses of different surgical methods
No large sample data analysis has reported the prognostic impact of different surgical methods for subependymoma.
Reviewing 466 patients with intracranial subependymoma, Nguyen et al. concluded that surgery was a significantly positive prognostic factor. However, the author further elaborated that GTR was not a significant prognostic factor and locations within “ventricles, NOS” or near “brain stem” were low-risk predictors factors for OS (16). Although there was the most extensive data analysis of intracranial subependymoma before 2017, the conclusion seemed to be counterintuitive.
Considering the importance of data quality, we screened the data of higher quality from the large sample and tried to make it equally comparable, through the application of statistical methods, such as eliminating incomplete data, univariate and multivariate analyses, PSM, and subgroup analysis.
To study the influence of different surgical methods on OS, we completed subgroup analyses through Cox regression analysis and PSM. Exhibited in Tables 3 and S1, patients matched after PSM had no significant difference in OS between the non-surgical and surgical groups. In the multivariable regression analysis, age and sex were significant prognostic variables for OS. The subgroup analysis of non-surgical versus surgical groups confirmed this finding in the training cohort.
Surgical treatment, race, year of diagnosis, insurance, tumor location, tumor size, pathology, tumor grade, and radiation had no statistically significant differences in OS for subependymoma in our research.
However, it was diacritical that age, sex, and laterality were the significant prognostic variables for OS in the different matched subgroups of unoperated versus GTR, STR, and biopsy. However, it did not mean that conservative treatment had a better prognosis.
The study of Nguyen et al. included fewer patients, lacking subgroup analyses, not clear prognostic factors of different treatment modalities, and the confounding factors were not matched. Improving deficiencies of previous studies, we should not be more inclined to choose conservative or surgical treatment.
Considering the rarity of the disease, our study is a retrospective analysis of the largest sample size of subependymoma to date, taking advantage of the SEER database’s wide population coverage.
We predicted the 3-, 5-, and 10-year survival rates for subependymoma using a nomogram model based on age and sex as prognostic factors. Although the prognostic model performed well in the experimental and validation groups, the two prognostic factors might not be sufficient to clinical use.
Nonetheless, the information that we present might be useful to suggest potential hypotheses to be tested in the clinical research setting. Doctors needed to evaluate the indications, contraindications, and risks of surgery comprehensively, and then made recommendations based on the wishes of patients’ families.
We suggested that the following measures needed to be adopted before the clinical implementation. Due to the extremely indolent nature of subependymoma, longer follow-up time was required to assess the outcome of the operative treatment. We need to expand the sample further and include more prognostic variables, such as immunohistochemical information. Prospective multicenter randomized controlled studies of subependymoma were needed to develop models with greater sensitivity and specificity.
Limitations
In our study, there was a particular patient selection bias based on the SEER database. Therefore, data quality is also a limitation of this study. Considering the rigor of the data, our study excluded patients with a survival time of less than 1 month, which might have skewed the results by excluding acute deaths from severe hydrocephalus without surgery. In addition, the median follow-up was only 56 months and the sample size was small after PSM. Longer follow-up and further multicenter studies with more sample sizes are needed.
Conclusion
Age and sex were the independent prognostic variables for OS in intracranial subependymoma. According to our research, we should not be more inclined to choose conservative or surgical treatment. Nonetheless, the information that we present might be useful to suggest potential hypotheses to be tested in the clinical research setting.
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 author.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent from the participants’ legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.
Author contributions
HC and ZZ: conception of the project. XP and YW: data reduction and screening. QL and HJ: statistical analysis and processing of data. ZZ, HC, and YW: writing and revising article. All authors contributed to the article and approved the submitted version.
Funding
This research project was funded by the Key Science Research Foundation of Traditional Chinese Medicine of Zhejiang Province of China (2019ZA089).
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/fonc.2022.939816/full#supplementary-material
Supplementary Table 1 | The χ2, Fisher’s exact or Mann-Whitney U test was used to inspect the differences of the non-surgical and surgical groups before and after matching. Age, insurance and pathology were analyzed by using the Mann-Whitney U test. The grade was analyzed by using Fisher’s exact test.
Supplementary Table 2 | The χ2, Fisher’s exact or Mann-Whitney U test was used to inspect the differences of the no surgery and GTR before and after matching. Age, insurance and pathology were analyzed by using the Mann-Whitney Utest. The grade was analyzed by using Fisher’s exact test.
Supplementary Table 3 | The χ2, Fisher’s exact or Mann-Whitney U test was used to inspect the differences of the no surgery and STR in before and after matching. Age, insurance and pathology were analyzed by using the Mann-Whitney U test. The grade was analyzed by using Fisher’s exact test.
Supplementary Table 4 | The χ2, Fisher’s exact or Mann-Whitney U test was used to inspect the differences of the no surgery and Surgery NOS or excisional biopsy before and after matching. Age, insurance and pathology were analyzed by using the Mann-Whitney U test. The grade was analyzed by using Fisher’s exact test.
References
1. Wu L, Yang T, Deng X, Yang C, Zhao L, Fang J, et al. Surgical outcomes in spinal cord subependymomas: an institutional experience. J J Neuro-Oncol (2014) 116(1):99–106. doi: 10.1007/s11060-013-1256-6
2. Bi Z, Ren X, Zhang J, Jia W, Jon J. Clinical, radiological, and pathological features in 43 cases of intracranial subependymoma. J Neurosurg (2015) 122(1):49–60. doi: 10.3171/2014.9.JNS14155
3. Tiwari N, Powell SZ, Takei H. Recurrent subependymoma of fourth ventricle with unusual atypical histological features: A case report. Pathol Int (2015) 65(8):438–42. doi: 10.1111/pin.12316
4. Vitanovics D, Áfra D, Nagy G, Hanzely Z, Banczerowski PJIS. Symptomatic subependymomas of the ventricles. Review of twenty consecutive cases. Ideggyogy Sz (2014) 67(11-12):415–9.
5. Hou Z, Wu Z, Zhang J, Zhang L, Tian R, Liu B, et al. Clinical features and management of intracranial subependymomas in children. J Clin Neurosci (2013) 20(1):84–8. doi: 10.1016/j.jocn.2012.05.026
6. Kandenwein JA, Bostroem A, Feuss M, Pietsch T, Simon MJAN. Surgical management of intracranial subependymomas. Acta Neurochir (2011) 153(7):1469–75. doi: 10.1007/s00701-011-1007-4
7. Rushing EJ, Cooper PB, Quezado M, Begnami M, Crespo A, Smirniotopoulos JG, et al. Subependymoma revisited: clinicopathological evaluation of 83 cases. J Neuro-oncol (2007) 85(3):297–305. doi: 10.1007/s11060-007-9411-6
8. Scheinker IM. Subependymoma: A newly recognized tumor of subependymal derivation. J Neurosurg (1945) 2(3):232–40. doi: 10.3171/jns.1945.2.3.0232
9. D'Amico RS, Praver M, Zanazzi GJ, Englander Z, Sims JS, Samanamud JL, et al. Subependymomas are low-grade heterogeneous glial neoplasms defined by subventricular zone lineage markers. World Neurosurg (2017) 107:451–63. doi: 10.1016/j.wneu.2017.08.009
10. Kammerer S, Mueller-Eschner M, Lauer A, Luger A-L, Quick-Weller J, Franz K, et al. Subependymomas–characteristics of a “Leave me alone” lesion. RöFo-Fortschritte Auf Dem Gebiet Röntgenstrahlen Bildgebenden Verfahren (2018) 190(10):955–66. doi: 10.1055/a-0576-1028
11. Zhou S, Xiong J, Pan J, Zhang T, Geng D, Zhang J. Neuroradiological features of cervical and cervicothoracic intraspinal subependymomas: a study of five cases. Clin Radiol (2016) 71(5):499.e9–.e15. doi: 10.1016/j.crad.2016.01.017
12. Jain A, Amin AG, Jain P, Burger P, Jallo GI, Lim M, et al. Subependymoma: clinical features and surgical outcomes. Neurol Res (2012) 34(7):677–84. doi: 10.1179/1743132812Y.0000000064
13. Louis DN, Arie P, Pieter W, Brat DJ, Cree IA, Dominique FB, et al. The 2021 WHO classification of tumors of the central nervous system: a summary. Neuro-Oncology (2021) 23(8):1231–51. doi: 10.1093/neuonc/noab106
14. Varma A, Giraldi D, Mills S, Brodbelt AR, Jenkinson MDJAN. Surgical management and long-term outcome of intracranial subependymoma. 2. Acta Neurochir (018) 160(9):1793–9. doi: 10.1007/s00701-018-3570-4
15. Kurian KM, Jones DTW, Marsden F, Openshaw S, Pearson DM, Ichimura K, et al. Genome-wide analysis of subependymomas shows underlying chromosomal copy number changes involving chromosomes 6, 7, 8 and 14 in a proportion of cases. Brain Pathol (2010) 18(4):469–73. doi: 10.1111/j.1750-3639.2008.00148.x
Keywords: intracranial subependymoma, SEER, nomogram, surgery, prognosis
Citation: Zhang Z, Pang X, Wei Y, Lv Q, Jin X and Chen H (2022) Clinical independent prognostic factors and overall survival prognostic nomogram for intracranial subependymoma: A SEER population-based analysis 2004–2016. Front. Oncol. 12:939816. doi: 10.3389/fonc.2022.939816
Received: 09 May 2022; Accepted: 25 July 2022;
Published: 22 August 2022.
Edited by:
Erik P Sulman, New York University, United StatesReviewed by:
Andrea Carai, Bambino Gesù Children’s Hospital (IRCCS), ItalySilvia Chiesa, Agostino Gemelli University Polyclinic (IRCCS), Italy
Copyright © 2022 Zhang, Pang, Wei, Lv, Jin and Chen. 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: Huai Chen, Y2g5OHpqQDE2My5jb20=