Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 25 August 2021
Sec. Neuroinflammation and Neuropathy

Prognostic Value of a New Integrated Parameter—Both Collateral Circulation and Permeability Surface—in Hemorrhagic Transformation of Middle Cerebral Artery Occlusion Acute Ischemic Stroke: Retrospective Cohort Study

\r\nChanchan LiChanchan Li1Xiaozhu HaoXiaozhu Hao1Luyi LinLuyi Lin1Chengfeng SunChengfeng Sun1Hai YuHai Yu2Zhenwei YaoZhenwei Yao1Xiaoyuan FengXiaoyuan Feng1Yanmei Yang*Yanmei Yang1*
  • 1Department of Radiology, Huashan Hospital, Fudan University, Shanghai, China
  • 2Department of Neurology, Huashan Hospital, Fudan University, Shanghai, China

Background: Multimodal CT, including CT angiography (CTA) and CT perfusion (CTP), was increasingly used in stroke triage. This study was to determine the relationship between a new integrated parameter—both collateral circulation and relative permeability surface (PS)—and the hemorrhagic transformation (HT) in acute ischemic stroke (AIS) with middle cerebral artery occlusion (MCAO).

Methods: We retrospectively reviewed consecutive AIS patients with MCAO who underwent baseline CTA/CTP within 4 h of symptom onset and follow-up susceptibility-weighted imaging (SWI) within 3 weeks. Collateral circulation was assessed on the baseline CTA. Baseline CTP data were postprocessed to generate PS parameter. The patients with poor collateral circulation and at the same time with high relative PS were classified as the group of both poor collateral circulation and high relative PS. HT was defined according to European Cooperative Acute Stroke Study II criteria on follow-up SWI imaging. Multivariate logistic regression analysis was performed using HT as an outcome variable.

Results: The group of patients with both poor collateral circulation and high relative PS was thirteen and thirty-three (52%) developed HT of the final cohort sixty-three AIS patients with MCAO. Multivariate logistic analysis revealed the new integrated parameter—both collateral circulation and relative PS (odds ratio, 16.59; 95% confidence interval, 13.09–19.10; P < 0.001) was independent predictor of HT. The area under the curve was 0.85 (95% confidence interval, 0.81–0.89). The sensitivity was 57%, specificity 97% and positive predictive value 92%, negative predictive value 58%.

Conclusions: For AIS patients with MCAO, these with poor collateral circulation on CTA and at the same time with high relative PS on CTP were at high risk for HT.

Introduction

Hemorrhagic transformation (HT) is an unwanted complication of acute ischemic stroke (AIS) that may severely worsen the prognosis (Knight et al., 2019). Multimodal CT, including CT angiography (CTA) and CT perfusion (CTP), is a part of routine acute stroke care in many stroke centers. Numerous studies have tried to determine clinical and imaging parameters associated with HT so as to identify patients at highest risk for thrombolytic or endovascular therapies. Indeed, risks factors for HT have been investigated and included higher age, higher stroke severity, large vessel occlusion and collateral score on CTA and permeability surface (PS) on CTP et al. (Hom et al., 2011; Bivard et al., 2016; Marsh et al., 2016; Wang et al., 2016; Kalinin et al., 2017; Li Q. et al., 2017; Li Y. et al., 2017; Liu et al., 2017; Puig et al., 2017). However, many studies taken it for granted yet could not confirm the cause-and-effect relationship between collateral circulation and PS for these two parameters were simultaneously acquired (Li Q. et al., 2017; Puig et al., 2017; Horsch et al., 2018). Therefore, it was inappropriate to draw a conclusion that inadequate collateral circulation contributed to the increased PS or to identify collateral circulation as a confounding factor in the HT prediction model.

In this study, we integrated the collateral circulation and PS into a new parameter—both collateral circulation and relative PS. We hypothesized that this new parameter contributes to HT and we aimed to determine the relationship between this new parameter and HT in AIS with middle cerebral artery occlusion (MCAO).

Materials and Methods

Study Population and Design

A single-center retrospective cohort study was conducted in consecutive AIS patients, according to Guidelines for Diagnosis and Treatment of Acute Ischemic Stroke in China 2014, presenting to the Stroke Unit of Huashan Hospital, Fudan University between January 2016 and December 2019. Patients were included if the following inclusion criteria were met: (1) age > 18 years and < 85 years; (2) National Institutes of Health Stroke Scale (NIHSS) was 4 to 22 at admission; (3) time from symptom onset to no-contrast CT (NCCT)/CTP/CTA was within 4 h and before recanalization treatment; (4) AIS with middle cerebral artery (MCA)-M1 occlusion documented on CTA; (5) the admission NCCT showed no evidence of intracerebral hemorrhage; (6) time from symptom onset to follow-up susceptibility-weighted imaging (SWI) was within 3 weeks. Patients with one or more of the following conditions were excluded: (1) history of brain tumor or brain surgery; (2) ischemic stroke in the bilateral cerebral hemispheres; (3) low imaging quality; (4) abnormal endovascular procedure including inappropriate thrombectomy and impaling vessels.

This study was approved by the Ethics Committee of Huashan Hospital, Fudan University (2020-939) and was registered at the Chinese Clinical Trial Registry (ChiCTR2000035575). Informed consent was waived for collection and analysis of preexisting data. This study was carried out in accordance with the Declaration of Helsinki.

Clinical Variables

The following clinical variables were collected for each patient: sex, age, admission NIHSS, risk factors (hypertension, hyperlipidemia, diabetes mellitus, atrial fibrillation, previous stroke, and current smoking), recanalization therapy (intravenous thrombolysis and endovascular thrombectomy) and time to recanalization therapy, time to imaging.

Patients without contraindication to thrombolysis received intravenous recombinant tissue plasminogen activator (rt-PA) treatment. Eligible patients without extensive early ischemic signs on pre-interventional scans received endovascular treatment. Patients with contraindications to thrombolysis received supportive care.

Image Acquisition

All patients underwent baseline NCCT/CTP/CTA within 4 h of symptom onset and follow-up SWI within 3 weeks. Whole-brain NCCT, CTP, and CTA were simultaneously performed on a 256-section CT scanner (Brilliance iCT, Phillips Medical Systems, Cleveland, OH, United States) as follows: Jog mode, 120 kVp, 150 mAs, field of view 220 mm, matrix 512 × 512. A dual-head power injector (Stellate Injection System, Indianola, PA, United States) was used to inject 40 ml of non-ionic contrast medium (Ultravist, iodine 370 mg I/ml; Bayer Healthcare, Berlin, Germany) at 5 ml/s followed by 20 ml saline with 5-s (CTP) and 8-s (CTA) delay. CTP included a 50-s scan reconstructed at 0.4-s intervals to produce a series of 312 sequential images for 13 sections, covering a total of 120 mm from the foramen magnum to the lateral ventricles. SWI was performed on a 3.0T superconducting magnetic resonance imaging (MRI) scanner (Verio, Siemens AG, Erlangen, Germany).

Image Analysis

Baseline CTA data were processed using the Phillips Brilliance Workspace portal software (Vision 5.0.2), including 10-mm axial and multi-planar maximum intensity projection reconstructions. The leptomeningeal collateral circulation was graded according to a previously reported scoring system on a scale from 0 to 3: 0 = absent collateral, 1 = collateral filling < 50% of the occluded territory, 2 = collateral filling > 50% < 100% of the occluded territory, and 3 = collateral filling 100% of the occluded territory (Figure 1). The results were dichotomized into poor collateral circulation (a score of 0 or 1) versus good collateral circulation (a score of 2 or 3) (Tan et al., 2009). MCA-M1 occlusion was defined as a main MCA trunk occlusion before the bifurcation, with or without ipsilateral internal carotid artery occlusion.

FIGURE 1
www.frontiersin.org

Figure 1. Examples of collateral circulation from 0 to 3 (white arrow). (A) 0 = absent collateral, (B) 1 = collateral filling < 50% of the occluded territory, (C) 2 = collateral filling > 50% < 100% of the occluded territory, and (D) 3 = collateral filling 100% of the occluded territory. All four patients were acute ischemic stroke (AIS) with right middle cerebral artery (MCA)-M1 occlusion. AIS, acute ischemic stroke; MCA, middle cerebral artery.

Postprocessing of baseline CTP source images was performed using Phillips brain perfusion software (Vision 5.0.2) to generate parametric maps of PS (Figure 2A). The adiabatic approximation tissue homogeneity model and a deconvolution technique were used to calculate CTP parameters. HT was defined as a hypointense susceptibility effect-induced area in the ischemic region on follow-up SWI (within 3 weeks) (Figure 2B). If present, HT was classified according to the European Cooperative Acute Stroke Study II classification into four subtypes, that is, hemorrhagic infarction (HI) 1 and 2 and parenchymal hematoma (PH)1 and 2. Symptomatic intracerebral hemorrhage (SICH) was identified as PH with worsening of neurological deficit ≥ 4 points on the NIHSS that was temporally related to its imaging appearance (Larrue et al., 2001). Baseline PS maps and follow-up SWI images were co-registered automatically using OpenCV 2.4.9 (Figure 2C). Then the regions of interest (ROIs) were drawn within the HT region and the corresponding contralateral region. This method made the ROI on the PS maps on the affected side were in the exact location where the HT happened. Relative PS were calculated by normalizing the values in the HT region with those in the contralateral mirror region. The results were dichotomized into high relative PS (≥2.89) versus low relative PS (< 2.89) (Li Q. et al., 2017).

FIGURE 2
www.frontiersin.org

Figure 2. (A) CT perfusion permeability surface (PS) map, follow-up (B) susceptibility-weighted imaging (SWI) in the same patient as in Figure 1B and (C) co-registered map. Obviously increased (A) PS on the affected side compared with the contralateral normal side. Follow-up (B) SWI showed hemorrhagic transformation (HT) (white arrow) which belonged to parenchymal hematoma (PH) 2. PS map and follow-up SWI image were co-registered automatically using OpenCV 2.4.9. Then the regions of interest (ROIs) were drawn within the HT region and the corresponding contralateral region. Relative PS were calculated by normalizing the values in the HT region with those in the contralateral mirror region. PS, permeability surface; SWI, susceptibility-weighted imaging; HT, hemorrhagic transformation; PH, parenchymal hematoma; ROIs, regions of interest.

Two radiologists (CL and XH) who were blinded to the study groups independently assessed the MCA-M1 status and collateral circulation. Two radiologists (LL and CS) assessed the follow-up SWI for HT and relative PS in selected ROIs. According to the new parameter—both collateral circulation and relative PS, the patients with poor collateral circulation and at the same time with high relative PS were classified as the group of both poor collateral circulation and high relative PS.

Statistical Analysis

Statistical analyses were performed using SPSS 19.0 software (SPSS, Chicago, IL, United States). Quantitative variable was expressed as mean ± standard deviation and categorical variable was expressed as proportion. Differences were assessed using the χ2 test (categorical variables), Student’s-t test (variables with normal distribution), or Mann–Whitney rank test (variables without a normal distribution). Univariate and multivariate logistic regression analysis were performed using HT as an outcome variable. Because we had only a limited number of outcomes and only three variables could be selected maximally to have at least 10 outcomes per variable. Variables analyzed were age, admission NIHSS, the new parameter—both collateral circulation and relative PS (Horsch et al., 2018). Interrater reliability was evaluated using weighted k values for ordered categorical variables. All P < 0.05 was considered significant.

Results

Sixty-three patients were included in this study from January 2016 to December 2019. Study flow chart was shown in Supplementary Figure I. The baseline characteristics had no significant difference between the patients involved and not involved in this study (Supplementary Table I). Of the final cohort comprising sixty-three AIS patients with M1-MCA occlusion, thirteen was the group of both poor collateral circulation and high relative PS and fifty belonged to the other. Thirty-three (52%) developed HT. ECASS II types were as follows: fifteen HI1(45%), eleven HI2(33%), four PH1(12%), three PH2(9%), and four SICH (12%). The mean time from symptom onset to HT was 5.21 ± 4.92 days. Interrater Reliability was excellent for collateral score (k = 0.95, P = 0.014) and relative PS (k = 0.91, P = 0.031).

Clinical characteristics were comparable between the patients with both poor collateral circulation and high relative PS (n = 13) and the other ones (n = 50) (Table 1). In the univariate logistic regression analysis, age, admission NIHSS, the new parameter—both collateral circulation and relative PS were all associated with HT. In the multivariate logistic regression analysis with three selected variables, only the new parameter—both collateral circulation and relative PS (odds ratio, 16.59; 95% confidence interval, 13.09–19.10; P < 0.001) was independent predictor of HT (Table 2). The area under the curve was 0.85 (95% confidence interval, 0.81–0.89). The sensitivity was 57%, specificity 97% and positive predictive value 92%, negative predictive value 58%.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of patients.

TABLE 2
www.frontiersin.org

Table 2. Univariate and multivariate logistic regression analysis for HT.

Discussion

Our results indicated that the new integrated parameter—both collateral circulation and relative PS was strongly associated with HT after AIS. This study highlighted the promise of using multimodal CT parameter to predict HT in AIS patients with MCA-M1 occlusion.

Although purely clinical factors are useful in the decision-making process before rt-PA administration, in practice, most of these are insufficiently powerful in isolation to alter decision-making regarding treatment with thrombolysis, apart from absolute contraindications. This has led to the recent development of various HT prediction scoring systems using a variety of clinical and imaging-based predictors to estimate the pretreatment risk of HT (Hom et al., 2011; Bivard et al., 2016; Marsh et al., 2016; Wang et al., 2016; Kalinin et al., 2017; Li Q. et al., 2017; Li Y. et al., 2017; Liu et al., 2017; Puig et al., 2017). MRI is not routinely available in the acute setting, whereas multimodal CT is widely and rapidly accessible in most stroke centers.

Hemorrhagic transformation is a dynamic and complex phenomenon, and its pathophysiology is still not clear. The underlying mechanism may relate to impairment of the autoregulatory capacity of the cerebral vasculature and disruption of the blood-brain barrier (BBB) (Hom et al., 2011; Wang et al., 2016; Li Q. et al., 2017; Li Y. et al., 2017; Puig et al., 2017). Our study showed the poor collateral circulation and at the same time the high relative PS provided a key link to a higher incidence of HT. To our knowledge, this integrated parameter—both collateral circulation and relative PS has not been reported previously. Our results supported the view that the extent of hemodynamic derangement can help identify those at highest risk for HT and HT was the result of severe ischemic damage to the vessel walls that leads to BBB disruption (Hom et al., 2011; Wang et al., 2016; Li Q. et al., 2017; Li Y. et al., 2017; Puig et al., 2017). This new integrated parameter showed 97% specificity and 92% positive predictive value in predicting HT of AIS patients with MCA-M1 occlusion. For AIS patients with MCAO, these with poor collateral circulation on CTA and at the same time with high relative PS on CTP were at high risk for HT. It is the strength of our study. This may help clinicians to tighten current rt-PA treatment indications to exclude patients who may previously have been included.

CT angiography maximum intensity projection reconstruction is effective to evaluate leptomeningeal collateral flow. After MCA-M1 occlusion, leptomeningeal collateral is rapidly recruited to supply the ischemic area (Li Q. et al., 2017). PS map on CTP is intended to reflect the degree of BBB disruption, although it is potentially fraught with substantial variability due to differences in acquisition and postprocessing techniques, including the size and molecular charge of the contrast agent used. To a certain extent, our findings prompted that collateral blood flow may be a potential therapeutic target for reducing the incidence of HT. However, there was one uncertainty that should be taken into consideration that our results could not confirm that poor collateral circulation contributed to the increased PS or good collateral circulation contributed to the decreased PS for these two parameters were simultaneously acquired in this study.

This study has many limitations. First, this is a single-center study with a small sample size, that is, an insufficient number of patients with HT to stratify into symptomatic versus asymptomatic groups, or PH versus HI groups. The parameters threshold may not be accurate. Second, this study has the limitations and potential biases inherent in a retrospective study. The duration of the CTP acquisition is relatively short (50 s), which may make it questionable if it is truly measuring BBB permeability because our study is retrospectively conducted with imaging performed for clinical purposes rather than imaging optimized for research purposes. Therefore, further prospective research is required.

Conclusion

In conclusion, the new integrated parameter—both collateral circulation and relative PS predicts HT in AIS patients with MCAO. The incorporation of readily available multimodal CT parameter into the clinical decision-making algorithm may allow for a greater appreciation of the risk-benefit profile before decisions about thrombolytic therapy in AIS. Further analysis with a larger number of patients and clinical end points in an independent data set is required to optimize and validate the methodology and potentially make this multimodal image protocol more applicable in the clinical setting.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Huashan Hospital, Fudan University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

CL, ZY, XF, and YY were responsible for study concept and the design. CL, XH, LL, CS, and HY collected the data. CL, XH, LL, and CS analyzed the data and conducted the statistical analysis. CL wrote the manuscript. ZY, XF, and YY supervised the manuscript and provided technical or information support. All authors approved the final version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Grant numbers: 81771788 and 81801660). None of these funding sources were involved in the study design, data collection, data analysis or interpretation, or the writing of the report and the decision to submit the article for publication.

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.

Acknowledgments

We would like to thank all the subjects for participating in this study and the support by the Huashan Hospital, Fudan University. We would also like to thank Shengjie Zhang for the help of image analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2021.703734/full#supplementary-material

Abbreviations

AIS, acute ischemic stroke; BBB, blood-brain barrier; CTA, CT angiography; CTP, CT perfusion; HI, hemorrhagic infarction; HT, hemorrhagic transformation; MCA, middle cerebral artery; MCAO, middle cerebral artery occlusion; MRI, magnetic resonance imaging; NCCT, no-contrast CT; NIHSS, National Institutes of Health Stroke Scale; PH, parenchymal hematoma; PS, permeability surface; ROI, regions of interest; rt-PA, recombinant tissue plasminogen activator; SICH, symptomatic intra cerebral hemorrhage; SWI, susceptibility-weighted imaging.

References

Bivard, A., Cheng, X., Lin, L. T., Levi, C., Spratt, N., Kleinig, T., et al. (2016). Global white matter hypoperfusion on CT predicts larger infarcts and hemorrhagic transformation after acute ischemic. CNS Neurosci. Ther. 22, 238–243. doi: 10.1111/cns.12491

PubMed Abstract | CrossRef Full Text | Google Scholar

Hom, J., Dankbaar, J. W., Soares, B. P., Schneider, T., Cheng, S. C., Bredno, J., et al. (2011). Blood-brain barrier permeability assessed by perfusion CT predicts symptomatic hemorrhagic transformation and malignant edema in acute ischemic stroke. Am. J. Neuroradiol. 32, 41–48. doi: 10.3174/ajnr.A2244

PubMed Abstract | CrossRef Full Text | Google Scholar

Horsch, A. D., Bennink, E., van Seeters, T., Kappelle, L. J., Graaf, Y. V. D., Mali, W. P. T. M., et al. (2018). Computed tomography perfusion derived blood-brain barrier permeability does not yet improve prediction of hemorrhagic transformation. Cerebrovasc. Dis. 45, 26–32. doi: 10.1159/000485043

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalinin, M. N., Khasanova, D. R., and Ibatullin, M. M. (2017). The hemorrhagic transformation index score: a prediction tool in middle cerebral artery ischemic stroke. BMC Neurol. 17:177. doi: 10.1186/s12883-017-0958-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Knight, G. A., Nario, J. J. Q., and Gupta, A. (2019). Cause of acute stroke: a patterned approach. Radiol. Clin. North Am. 57, 1093–1108. doi: 10.1016/j.rcl.2019.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Larrue, V., von Kummer, R., Muller, A., and Bluhmki, E. (2001). Risk factors for severe hemorrhagic transformation in ischemic stroke patients treated with recombinant tissue plasminogen activator: a secondary analysis of the European-Australasian acute stroke study (ECASS II). Stroke 32, 438–441. doi: 10.1161/01.str.32.2.438

CrossRef Full Text | Google Scholar

Li, Q., Gao, X. Y., Yao, Z. W., Feng, X. Y., He, H. J., Xue, J., et al. (2017). Permeability surface of deep middle cerebral artery territory on computed tomographic perfusion predicts hemorrhagic transformation after stroke. Stroke 48, 2412–2418. doi: 10.1161/STROKEAHA.117.017486

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Xia, Y., Liu, N., Jackson, A., Wintermark, M., Zhang, Y., et al. (2017). Focal low and global high permeability predict the possibility, risk, and location of hemorrhagic transformation following intra-arterial thrombolysis therapy in acute stroke. Am. J. Neuroradiol. 38, 1730–1736. doi: 10.3174/ajnr.A5287

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Wu, B., Zhao, J. L., Cao, Y. Y., Dedhia, N., Caplan, L. R., et al. (2017). Computed tomography perfusion alberta stroke program early computed tomography score in associated with hemorrhagic transformation after acute cardioembolic stroke. Front. Neurol. 8:591. doi: 10.3389/fneur.2017.00591

PubMed Abstract | CrossRef Full Text | Google Scholar

Marsh, E. B., Llinas, R. H., Schneider, A. L. C., Hillis, A. E. H., Lawrence, E., Dziedzic, P., et al. (2016). Predicting hemorrhagic transformation of acute ischemic stroke: prospective validation of the HeRS score. Medicine 95:e2430. doi: 10.1097/MD.0000000000002430

PubMed Abstract | CrossRef Full Text | Google Scholar

Puig, J., Blasco, G., Daunis, P., van Eendendburg, C., Carrillo, M., Aboud, C., et al. (2017). High-permeability region size on perfusion CT predicts hemorrhagic transformation after intravenous thrombolysis in stroke. PLoS One. 12:e0188238. doi: 10.1371/journal.pone.0188238

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, I. Y. L., Demchuk, A. M., Hopyan, J., Zhang, J., Gladstone, D., Wong, K., et al. (2009). CT angiography clot burden score and collateral score: correlation with clinical and radiologic outcomes in acute middle cerebral artery infarct. Am. J. Neuroradiol. 30, 525–531. doi: 10.3174/ajnr.A1408

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Li, M. C., Chen, Q. X., and Wang, J. (2016). Hemorrhagic transformation after tissue plasminogen activator reperfusion therapy for ischemic stroke: mechanisms, models, and biomarkers. Cell. Mol. Neurobiol. 52, 1572–1579. doi: 10.1007/s12035-014-8952-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: collateral circulation, stroke, angiography, perfusion imaging, cerebral hemorrhage

Citation: Li C, Hao X, Lin L, Sun C, Yu H, Yao Z, Feng X and Yang Y (2021) Prognostic Value of a New Integrated Parameter—Both Collateral Circulation and Permeability Surface—in Hemorrhagic Transformation of Middle Cerebral Artery Occlusion Acute Ischemic Stroke: Retrospective Cohort Study. Front. Aging Neurosci. 13:703734. doi: 10.3389/fnagi.2021.703734

Received: 30 April 2021; Accepted: 26 July 2021;
Published: 25 August 2021.

Edited by:

Rubem C. A. Guedes, Federal University of Pernambuco, Brazil

Reviewed by:

Norbert Nighoghossian, Université Claude Bernard Lyon 1, France
Chushuang Chen, The University of Newcastle, Australia

Copyright © 2021 Li, Hao, Lin, Sun, Yu, Yao, Feng and Yang. 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: Yanmei Yang, yym9876@sohu.com

Disclaimer: 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.