- 1Thyroid Center, Shanghai Tenth People's Hospital, Tongji University School of Medicine, Shanghai, China
- 2Shanghai Center for Thyroid Diseases, Shanghai, China
- 3Department of Biochemistry and Molecular Biology, The Faculty of Basic Medical Science, Second Military Medical University, Shanghai, China
- 4Department of Gastroenterology, Gansu Provincial Hospital, Lanzhou, China
- 5Department of Nuclear Medicine, Shanghai Tenth People's Hospital, Tongji University School of Medicine, Shanghai, China
- 6Department of Thoracic Surgery, Xuzhou Central Hospital, The Affiliated Xuzhou Hospital of Medical College of Southeast University, Xuzhou, China
Purpose: Currently, of the two most common staging systems of pancreatic neuroendocrine tumors (pNETs) one is from the European Neuroendocrine Tumor Society (ENETS) and the other is from the American Joint Committee on Cancer (AJCC). However, there are imperfections in both these staging systems.
Patients and methods: Patients were selected retrospectively from the Surveillance Epidemiology and End Results (SEER) database (2004 to 2013). The effect of age on the hazard ratio (HR) was evaluated using restricted cubic splines. The discriminatory power of the staging systems was determined using the concordance index (C-index).
Results: A total of 3,034 patients with pNETs were included in the final analyses. The risk of death increased slowly along with age for patients under 60 years of age, but the risk of death rose sharply for those over 60 years of age, forming a mirrored L-shaped survival curve. In the current AJCC tumor-node-metastasis (TNM) staging system, no statistical significance was observed between stages IA and IB (p = 0.105). Patients with stage IIB even had longer OS than patients with IIA, although there was no statistical significance (p = 0.574). The proportion of stage III patients was small (2.7%). In the proposed aTNM staging system, significant survival differences could be observed among stage I, IIA, and IIB (p < 0.001) and the proportion of stage III rose from 2.7 to 25.7%.
Conclusion: Our findings suggest that age has a critical influence on the survival of patients with pNETs. Age should be considered as a factor in future staging systems of pNETs.
Introduction
Pancreatic neuroendocrine tumors (pNETs) are rare, with an incidence rate of >0.5/100,000 in 2004 in the United States (1), and represent ~1–2% of all pancreatic neoplasms (2). Most pNETs are sporadic and may also occur with some inherited genetic syndromes, the most two common of which were multiple endocrine neoplasia (MEN) types 1 and MEN2. The former are related to the mutated menin gene and the latter are associated with the mutated RET gene. pNETs are usually indolent and have a more favorable outcome (3). However, factors such as the diagnosis, treatment and prognosis of pNETs have not been clearly understood compared with other pancreatic exocrine tumors. Therefore, a well-defined, useful classification is needed to better manage patients with pNETs.
pNETs can be classified into different groups depending on different characteristics. According to the molecular mutations, pNETs are divided into sporadic pNETs and hereditary pNETs. According to the functionality, pNETs are distinguished in functional and non-functional. Among the functional pNETs, insulinomas account for 70% approximately and 15% are glucagonomas. Other rare functional pNETs include pancreatic polypeptidoma (PPoma), VIPomas secreting vasoactive intestinal polypeptide (VIP) and cholecystokininoma (CCKoma). According to the histological types, pNETs are divided into Grade 1, Grade 2 and Grade 3 based on ki-67 index, mitotic count, and proliferative activity of tumor cells (4). To guide clinical practice, of the two most common staging systems for pNETs, one is from the European Neuroendocrine Tumor Society (ENETS) and the other is from the American Joint Committee on Cancer (AJCC). However, through the ENETS staging system a similar outcome is observed for stage I and stage IIA (5, 6). Luo et al. (7) attempted to modify the ENETS (mENETS) staging system by adopting AJCC staging definitions and ENETS T, N, and M definitions. Unfortunately, in the mENETS staging system patients with stage IA still had a similar prognosis to those with stage IB. In the AJCC staging system, only a small percentage of patients, ranging from 4 to 5.3%, had a stage III prognosis (8, 9). Even in the mENETS staging system, stage III prognosis accounted for merely 8.9% of patients (7). These observations suggest that current staging systems (AJCC, ENETS and mENETS) still need improvement.
Multiple studies have shown that the main prognostic factors of pNETs include tumor stage, histological grade (6), age (10, 11), functional status (11), surgical margin (12), and metastatic pattern (e.g., diffuse liver metastases, extrahepatic metastases, bone metastases) (13–16). Halfdanarson et al. (11) divided patients with pNETs into 4 groups based on age: 18–50, 51–60, 61–70, and 71–95 years. Their median overall survival (OS) was 52, 44, 19, 9.5 months, respectively. Therefore, we hypothesized that the older patients got, the worse their prognosis was. Additionally, we noticed that age was included in the AJCC tumor-node-metastasis (TNM) staging system of papillary thyroid cancer (PTC) because of the observation of a dramatic increase in the risk of death beginning around the age interval of 50–60 years (17). Thus, this study aims to validate our hypothesis and explore whether age has an influence on the AJCC TNM classification of pNETs similar to the influence it had on PTC, using the Surveillance Epidemiology and End Results (SEER) database.
Materials and Methods
Ethics Statement
This study was approved by the institutional review board of Shanghai Tenth People's Hospital, Tongji University School of Medicine. Patients from the Surveillance, Epidemiology, and End Results (SEER) database had previously consented to participate in any scientific research worldwide.
Patients
The SEER database (2004 to 2013) was used to identify pNETs patients. We defined pNETs to include the following International Classification of Diseases for Oncology third edition (ICD-O3) codes: islet cell carcinoma (8150), insulinoma (8151), glucagonoma (8152), gastrinoma (8153), vipoma (8155), somatostatinoma (8156), enteroglucagonoma (8157), carcinoid (8240), enterochromaffin cell carcinoid (8241), enterochromaffin-like cell tumors (8242), goblet cell carcinoid (8243), composite carcinoid (8244), adenocarcinoid (8245), neuroendocrine carcinoma (8246), and atypical carcinoid (8249) (18). TNM classifications based on two codes: the derived AJCC stage group (6th edition; 2004+) and the derived AJCC stage group (7th edition; 2010+) were retrieved. Study enrollment criteria were patients who had a positive histological diagnosis of pNETs, definite TNM information and patients who had survival data available. The study cohort didn't distinguish radical or palliative patients. Patients were excluded if there were coexisting pancreatic adenocarcinoma. Clinical and pathological features included age, sex, race, functional status, location of the primary tumor, grade, and TNM stage.
Statistical Analysis
Multivariate analysis, including sex, age, tumor type, stage, grade, tumor location and race, was done using the Cox proportional hazard regression model and log-rank tests were used to evaluate prognostic factors. The stepwise backward procedure based on the likelihood ratio was used in the Cox model. The hazard ratio (HR) and a 95% confidence interval (CI) were calculated. The effect of age on the HRs was evaluated using restricted cubic splines in R i386 3.3.2 software (running rms package), in which age was regarded as continuous variable. Survival analysis was conducted using Kaplan-Meier methods in PASW Statistics 18 software to assess the prognosis of pNETs patients with different stages. The concordance index (C-index) was used to determine the discriminatory ability of the staging systems. Statistical significance was defined as a two-sided p < 0.05.
Results
Patient Characteristics
Clinical and pathological features of patients enrolled in the study are listed in Table 1. A total of 3,034 patients with pNETs were included and comprised of 1,655 males and 1,379 females. The median age of all patients was 59 years (ranging from 11 to 97 years). Primary tumors were most frequently located at the body-tail of the pancreas (n = 1,350, 44.5%), followed by the head (n = 907, 29.9%). The primary tumor site, mainly included the pancreatic duct and islets of Langerhans while other specific parts of the pancreas were defined as others (n = 777, 25.6%). Among the available data, 1,547 cases were described on pathologic reports as well-differentiated (Grade I) or moderately differentiated (Grade II) (I-II−52.3%), and 281 were poorly differentiated (Grade III) or undifferentiated (Grade IV) (III-IV−9.3%). Apart from patients with unknown grade information, Grade I-II accounted for a majority (85.0%, 1588/1869). According to the 7th AJCC staging system, the frequency, and proportion of stage I, II, III, IV was 879 (29.0%), 607 (20.0%), 83 (2.7%), and 1,465 (48.3%), respectively.
Prognostic Factors of Survival
Multivariate analysis demonstrated that among factors being considered age (≤60 vs. >60), grade (I-II vs. III-IV), functional status (functional type vs. non-functional type), race (white vs. black), location (head vs. body-tail), and stage (Table 2) were independently associated with survival.
Table 2. Prognostic significance for OS by multivariate analysis of variables for patients with pNETs using the Cox proportional hazard regression model.
Effect of Age on Staging
Mean age of patients was 58.7 years (interquartile range, 50 to 69 years, Figure 1), and approximately equaled median age. In a univariable Cox model for OS, the effect of age was shown by a mirrored L-shaped curve (Figure 2), the first half steadily increasing and the latter half showing a sharp increase. The turning point was approximately around the age of 60 years. As shown in Table 2, the hazard ratio (HR) of death for patients over 60 years was 1.753 in reference to patients under 60 (p < 0.001). Notably, the oldest group (>75 years) showed an increase in the risk of death as high as 235.6% over that of the youngest group (<40 years) (HR- 3.356; 95% CI−2.544–4.426, p < 0.001). Thus, we included age as a factor in the TNM staging system. The principle of modifying the TNM staging system was that patients over 60 years of age were upstaged as the next advanced stage (Figure 3), resulting in modified stages of I, IIA, IIB, III, and IV. To increase the discriminatory ability, stage IA and stage IB were classified as stage I. Supplementary Figure 1 provided another more detailed version of modified stages in which stage IA, stage IB, stage IIA, and stage IIB were regarded as an individual group, respectively.
Figure 2. Risk of death with age growing by R software (blue line stands for HR and gray area for 95% CI).
Figure 3. The principle diagram of modifying AJCC TNM staging system based on age: stage I patients under 60 years were classified as stage I in the modified staging system; stage II patients under 60 years and stage I patients older than 60 years were classified as stage IIA in the modified staging system; stage III patients under 60 years and stage II patients older than 60 years were classified as stage IIB in the modified staging system; stage IV patients under 60 years and stage III patients older than 60 years were classified as stage III in the modified staging system; stage IV patients older than 60 years were classified as stage IV in the modified staging system.
The difference in survival between each stage was assessed according to the AJCC TNM staging system. Between stage IA and IB, no statistical significance was observed (p = 0.105, Figure 4A). Patients with stage IIB even had longer OS than patients with IIA, although there was no statistical significance (p = 0.574, Figure 4A). Significant differences were observed among stage II, III, and IV patients (p < 0.001), but the survival gap between stage II and III was huge (Figure 4A). According to the modified TNM staging system that included age (aTNM), a significant difference could be observed among all stages (p < 0.001, Figure 4B). Survival gaps between neighboring stages were generally uniformly distributed. The C-index of the TNM staging system was 0.720 (95% CI−0.685–0.755) while the C-index of the aTNM staging system was 0.753 (95% CI−0.718–0.788). Additionally, the proportion of stage III patients rose from 2.7 to 25.7% (Table 3).
Figure 4. (A) Survival curve of patients with pNETs according to the AJCC TNM staging system. (B) Survival curve of patients with pNETs according to the aTNM staging system.
Table 3. Comparison of pNET patient distributions between TNM and aTNM staging system by each stage.
Discussion
pNETs are neoplasms that arise from hormone producing pancreatic islet cells (19). Until now, the influence of age on pNETs has remained to be further elucidated. In Halfdanarson's report (11), age was defined as categorical variable rather than continuous variable, which cannot fully reflect the true effect of age on the survival. This is the first study regarding age as continuous variable to investigate the effect of age on the prognosis of pNETs. The results of this study suggest that age is an independent prognostic factor of pNETs. The effect of age on the survival of patients with pNETs was specific showing a mirrored L-shaped curve. The risk of death increased slowly along with age for patients under 60 years of age. However, for patients over 60 years of age the risk of death rose sharply. Our findings indicate that patients over 60 years are a different population from those under 60 years of age. Thus, it was necessary to incorporate age into the staging system. The clinical guidelines recommended that patinets with non-functional pNETs <1 cm could be safely followed (20), which might be questioned since older patients (>60 years) showed a higher death risk. Additionally, older patients should undergo more frequent surveillance than younger patients.
Bilimoria et al. (21) first introduced the AJCC staging system for pancreatic adenocarcinoma into the staging for pNETs. However, it couldn't stratify patients with pNETs by death risk as properly as it did among patients with pancreatic adenocarcinoma. In the AJCC staging system for pNETs, the first two stages (stage IA and IB) had similar survival outcomes. When taking age into consideration, significant survival differences could be observed between the first two stages of the aTNM staging system (with stage I as the reference: HR for stage IIA−2.591, p < 0.001). In the AJCC staging system, there were a relatively low proportion of patients with stage III, but with the aTNM staging system, the proportion of stage III rose from 2.7 to 25.7%. Stage III patients with distant metastasis under 60 years that caused this increase were classified as stage IV in the TNM staging system.
Some studies (22–24) have argued that age has no significant impact on the survival of pNET patients, which has mainly resulted from studies with small sample sizes. In studies with large sample sizes or multicenter data, advanced age was closely related with shorter survival time (11, 25). The mechanism of how advanced age increases mortality is not yet clear. The mTOR pathway has been reported to play a critical role in senescence and senescence-related diseases (26, 27) and was frequently altered in pNETs (28). Everolimus, an mTOR inhibitor, has been proven effective in the treatment of pNETs in clinical practice (29). These findings suggest that the mTOR pathway may be the potential mechanism.
There were some limitations in this study. First, our results were based on retrospective data. Second, external data are needed to validate current conclusions. Third, a lack of genetic information hampered further exploration of the mechanisms underlying the influence of age on the prognosis of patients with pNETs. Fourth, progression-free survival data were needed to further evaluate the true effect of age on disease recurrence. Therefore, more prospective data and basic studies are needed to investigate in depth the effect of age on pNETs.
Data Availability
Publicly available datasets were analyzed in this study. This data can be found here: https://seer.cancer.gov/.
Ethics Statement
The studies involving human participants were reviewed and approved by the institutional review board of Shanghai Tenth People's Hospital, Tongji University School of Medicine. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
ZW and YZ made substantial contributions to the design of the study, carried out the analysis, interpreted the data. WJ and LZ contributed to the review of previous literature. JY and JD contributed substantially to the data discussion and critically commented on the manuscript for scientific content. All authors made substantial contributions to the conception and design of the study, data interpretation and drafting of the manuscript, were responsible for the quality of the overall manuscript. All authors approved the final version of the manuscript.
Funding
This research was supported by the National Science Foundation of China (41776140). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest Statement
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.00906/full#supplementary-material
Supplementary Figure S1. Survival curve of patients with pNETs according to the aTNM staging system, in which stage IA, stage IB, stage IIA, and stage IIB were regarded as an individual group respectively.
References
1. Yao JC, Hassan M, Phan A, Dagohoy C, Leary C, Mares JE, et al. One hundred years after “carcinoid”: epidemiology of and prognostic factors for neuroendocrine tumors in 35,825 cases in the United States. J Clin Oncol. (2008) 26:3063–7. doi: 10.1200/JCO.2007.15.4377
2. Milan SA, Yeo CJ. Neuroendocrine tumors of the pancreas. Curr Opin Oncol. (2012) 24:46–55. doi: 10.1097/CCO.0b013e32834c554d
3. Bilimoria KY, Talamonti MS, Tomlinson JS, Stewart AK, Winchester DP, Ko CY, et al. Prognostic score predicting survival after resection of pancreatic neuroendocrine tumors: analysis of 3851 patients. Ann Surg. (2008) 247:490–500. doi: 10.1097/SLA.0b013e31815b9cae
4. Bosman FT, Carneiro F, Hruban RH, Theise ND (eds). WHO Classification of Tumors of the Digestive System. Lyon: International Agency for Research on Cancer (IRAC) (2010). p. 13–14.
5. Ekeblad S, Skogseid B, Dunder K, Oberg K, Eriksson B. Prognostic factors and survival in 324 patients with pancreatic endocrine tumor treated at a single institution. Clin Cancer Res. (2008) 14:7798–803. doi: 10.1158/1078-0432.CCR-08-0734
6. Scarpa A, Mantovani W, Capelli P, Beghelli S, Boninsegna L, Bettini R, et al. Pancreatic endocrine tumors: improved TNM staging and histopathological grading permit a clinically efficient prognostic stratification of patients. Mod Pathol. (2010) 23:824–33. doi: 10.1038/modpathol.2010.58
7. Luo G, Javed A, Strosberg JR, Jin K, Zhang Y, Liu C, et al. Modified staging classification for pancreatic neuroendocrine tumors on the basis of the american joint committee on cancer and european neuroendocrine tumor society systems. J Clin Oncol. (2017) 35:274–80. doi: 10.1200/JCO.2016.67.8193
8. Rindi G, Falconi M, Klersy C, Albarello L, Boninsegna L, Buchler MW, et al. TNM staging of neoplasms of the endocrine pancreas: results from a large international cohort study. J Natl Cancer Inst. (2012) 104:764–77. doi: 10.1093/jnci/djs208
9. Strosberg JR, Cheema A, Weber J, Han G, Coppola D, Kvols LK, et al. Prognostic validity of a novel American Joint Committee on Cancer Staging Classification for pancreatic neuroendocrine tumors. J Clin Oncol. (2011) 29:3044–9. doi: 10.1200/JCO.2011.35.1817
10. Lepage C, Rachet B, Coleman MP. Survival from malignant digestive endocrine tumors in England and Wales: a population-based study. Gastroenterology. (2007) 132:899–904. doi: 10.1053/j.gastro.2007.01.006
11. Halfdanarson TR, Rabe KG, Rubin J, Petersen GM. Pancreatic neuroendocrine tumors (PNETs): incidence, prognosis and recent trend toward improved survival. Ann Oncol. (2008) 19:1727–33. doi: 10.1093/annonc/mdn351
12. Kleine M, Schrem H, Vondran FW, Krech T, Klempnauer J, Bektas H. Extended surgery for advanced pancreatic endocrine tumours. Br J Surg. (2012) 99:88–94. doi: 10.1002/bjs.7681
13. Panzuto F, Merola E, Rinzivillo M, Partelli S, Campana D, Iannicelli E, et al. Advanced digestive neuroendocrine tumors: metastatic pattern is an independent factor affecting clinical outcome. Pancreas. (2014) 43:212–8. doi: 10.1097/MPA.0000000000000032
14. Solorzano CC, Lee JE, Pisters PW, Vauthey JN, Ayers GD, Jean ME, et al. Nonfunctioning islet cell carcinoma of the pancreas: survival results in a contemporary series of 163 patients. Surgery. (2001) 130:1078–85. doi: 10.1067/msy.2001.118367
15. Yu F, Venzon DJ, Serrano J, Goebel SU, Doppman JL, Gibril F, et al. Prospective study of the clinical course, prognostic factors, causes of death, and survival in patients with long-standing Zollinger-Ellison syndrome. J Clin Oncol. (1999) 17:615–30. doi: 10.1200/JCO.1999.17.2.615
16. Gupta S, Johnson MM, Murthy R, Ahrar K, Wallace MJ, Madoff DC, et al. Hepatic arterial embolization and chemoembolization for the treatment of patients with metastatic neuroendocrine tumors: variables affecting response rates and survival. Cancer. (2005) 104:1590–602. doi: 10.1002/cncr.21389
17. Amin MB, Edge S, Greene F, Byrd DR, Brookland RK, Washington MK, et al. AJCC Cancer Staging Manual, 8th edn. New York, NY: Springer (2017). doi: 10.1007/978-3-319-40618-3_2
18. Conrad C, Kutlu OC, Dasari A, Chan JA, Vauthey JN, Adams DB, et al. Prognostic value of lymph node status and extent of lymphadenectomy in pancreatic neuroendocrine tumors confined to and extending beyond the pancreas. J Gastrointest Surg. (2016) 20:1966–74. doi: 10.1007/s11605-016-3243-7
19. Zhang J, Francois R, Iyer R, Seshadri M, Zajac-Kaye M, Hochwald SN, et al. Current understanding of the molecular biology of pancreatic neuroendocrine tumors. J Natl Cancer Inst. (2013) 105:1005–17. doi: 10.1093/jnci/djt135
20. National Comprehensive Cancer Network (NCCN). NCCN Clinical Practice Guidelines in Oncology. Neuroendocrine Tumors. Version 3.2018. Available online at: https://www.nccn.org/
21. Bilimoria KY, Bentrem DJ, Merkow RP, Tomlinson JS, Stewart AK, Ko CY, et al. Application of the pancreatic adenocarcinoma staging system to pancreatic neuroendocrine tumors. J Am Coll Surg. (2007) 205:558–63. doi: 10.1016/j.jamcollsurg.2007.05.009
22. Lee L, Igarashi H, Fujimori N, Hijioka M, Kawabe K, Oda Y, et al. Long-term outcomes and prognostic factors in 78 Japanese patients with advanced pancreatic neuroendocrine neoplasms: a single-center retrospective study. Jpn J Clin Oncol. (2015) 45:1131–8. doi: 10.1093/jjco/hyv143
23. Strosberg J, Gardner N, Kvols L. Survival and prognostic factor analysis in patients with metastatic pancreatic endocrine carcinomas. Pancreas. (2009) 38:255–8. doi: 10.1097/MPA.0b013e3181917e4e
24. Primavesi F, Klieser E, Cardini B, Marsoner K, Fröschl U, Thalhammer S, et al. Exploring the surgical landscape of pancreatic neuroendocrine neoplasia in Austria: results from the ASSO pNEN study group. Eur J Surg Oncol. (2019) 45:198–206. doi: 10.1016/j.ejso.2018.08.016
25. Fang C, Wang W, Zhang Y, Feng X, Sun J, Zeng Y, et al. Clinicopathologic characteristics and prognosis of gastroenteropancreatic neuroendocrine neoplasms: a multicenter study in South China. Chin J Cancer. (2017) 36:51. doi: 10.1186/s40880-017-0218-3
26. Xu Y, Li N, Xiang R, Sun P. Emerging roles of the p38 MAPK and PI3K/AKT/mTOR pathways in oncogene-induced senescence. Trends Biochem Sci. (2014) 39:268–76. doi: 10.1016/j.tibs.2014.04.004
27. Tchetina EV, Poole AR, Zaitseva EM, Sharapova EP, Kashevarova NG, Taskina EA, et al. Differences in Mammalian target of rapamycin gene expression in the peripheral blood and articular cartilages of osteoarthritic patients and disease activity. Arthritis. (2013) 2013:461486. doi: 10.1155/2013/461486
28. Jiao Y, Shi C, Edil BH, de Wilde RF, Klimstra DS, Maitra A, et al. DAXX/ATRX, MEN1, and mTOR pathway genes are frequently altered in pancreatic neuroendocrine tumors. Science. (2011) 331:1199–203. doi: 10.1126/science.1200609
29. Pavel ME, Hainsworth JD, Baudin E, Peeters M, Hörsch D, Winkler RE, et al. Everolimus plus octreotide long-acting repeatable for the treatment of advanced neuroendocrine tumours associated with carcinoid syndrome (RADIANT-2): a randomised, placebo-controlled, phase 3 study. Lancet. (2011) 378:2005–12. doi: 10.1016/S0140-6736(11)61742-X
Keywords: pancreatic neuroendocrine tumor, AJCC, overall survival, SEER, age
Citation: Wang Z, Jiang W, Zheng L, Yan J, Dai J, Huang C, Zhang Q, Yin Z, Gong X and Zhang Y (2019) Consideration of Age Is Necessary for Increasing the Accuracy of the AJCC TNM Staging System of Pancreatic Neuroendocrine Tumors. Front. Oncol. 9:906. doi: 10.3389/fonc.2019.00906
Received: 18 June 2019; Accepted: 02 September 2019;
Published: 19 September 2019.
Edited by:
Christian Macutkiewicz, Manchester University NHS Foundation Trust (MFT), United KingdomReviewed by:
Florian Primavesi, Innsbruck Medical University, AustriaZiv Radisavljevic, Brigham and Women's Hospital, Harvard Medical School, United States
Copyright © 2019 Wang, Jiang, Zheng, Yan, Dai, Huang, Zhang, Yin, Gong and Zhang. 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: Zhiqiang Yin, OTcyNjgzMDA0JiN4MDAwNDA7cXEuY29t; Xiangnan Gong, MTgxMTIwMDAzNjEmI3gwMDA0MDsxNjMuY29t; Yun Zhang, eXVuemhhbmcxMDI0JiN4MDAwNDA7MTI2LmNvbQ==
†These authors have contributed equally to this work