key: cord-0077444-fwu2qa4b authors: Du, Houwei; Wu, Sangru; Lei, Hanhan; Ambler, Gareth; Werring, David J.; Li, Hangfeng; Chen, Yangui; Chen, Dongping; Yuan, Qilin; Fang, Shuangfang; Chen, Ronghua; Zhang, Yixian; Wei, Jin; Chen, Guangliang; Chen, Jianhua; Liu, Nan title: Total Cerebral Small Vessel Disease Score and Cerebral Bleeding Risk in Patients With Acute Stroke Treated With Intravenous Thrombolysis date: 2022-04-11 journal: Front Aging Neurosci DOI: 10.3389/fnagi.2022.790262 sha: cf00fce78b77907575b455427b3f5c0ae3e6034c doc_id: 77444 cord_uid: fwu2qa4b OBJECTIVE: The aim of this study was to investigate the association of total cerebral small vessel disease (cSVD) score with the risk of intracerebral hemorrhage (ICH) in patients with acute ischemic stroke who received intravenous thrombolysis (IVT) using recombinant tissue-plasminogen activator (rt-PA). METHODS: We retrospectively reviewed clinical data from two stroke registries of patients with acute ischemic stroke treated with IVT. We assessed the baseline magnetic resonance (MR) visible cSVD markers and total cSVD score (ranging from 0 to 4) between patients with and without ICH after IVT. Logistic regression analysis was used to determine the association of total cSVD score with the risk of ICH after IVT, adjusted for cofounders selected by least absolute shrinkage and selection operator (LASSO). We additionally performed an E-value analysis to fully explain away a specific exposure-outcome association. Receiver operating characteristic (ROC) curve analysis was used to quantify the predictive potential of the total cSVD score for any ICH after IVT. RESULTS: Among 271 eligible patients, 55 (20.3%) patients experienced any ICH, 16 (5.9%) patients experienced a symptomatic ICH (sICH), and 5 (1.85%) patients had remote intracranial parenchymal hemorrhage (rPH). Logistic regression analysis showed that the risk of any ICH increased with increasing cSVD score [per unit increase, adjusted odds ratio (OR) 2.03, 95% CI 1.22–3.41, P = 0.007]. Sensitivity analyses using E-value revealed that it would need moderately robust unobserved confounding to render the exposure-outcome (cSVD-any ICH) association null. ROC analysis showed that compared with the National Institutes of Health Stroke Scale (NIHSS) score alone, a combination of cSVD and NIHSS score had a larger area under the curve for any ICH (0.811, 95% CI 0.756–0.866 vs. 0.784, 95% CI 0.723–0.846, P = 0.0004). CONCLUSION: The total cSVD score is associated with an increased risk of any ICH after IVT and improves prediction for any ICH compared with NIHSS alone. Current stroke guidelines recommend intravenous thrombolysis (IVT) using recombinant tissue-type plasminogen activator (rt-PA) as the first-line treatment for acute ischemic stroke (Chinese Society of Neurology and Chinese Stroke Society, 2018; Powers et al., 2018; Berge et al., 2021) . However, approximately 7.0-27.0% of patients treated with IVT using rt-PA experienced intracerebral hemorrhage (ICH) (Hacke et al., 2008; Dorňák et al., 2018; Anderson et al., 2019) , which might be related to unfavorable clinical outcomes. Previous efforts have been made to identify several risk factors for developing ICH after IVT, including aging, hypertension, atrial fibrillation, and stroke severity (Mazya et al., 2012; Whiteley et al., 2012) . Cerebral small vessel disease (cSVD) refers to a group of age-related (but not clearly age-derived) pathological processes with various etiologies that affect the microvasculature in the brain and might relate to intracranial bleeding risk (Pantoni et al., 2014; Charidimou et al., 2016) . Several studies showed that cSVD markers including cerebral microbleeds (CMBs) and white matter hyperintensities (WMH) were predictors for any ICH after IVT (Prats-Sánchez et al., 2016; Charidimou et al., 2017; Capuana et al., 2020) . However, a retrospective study of 434 patients with ischemic stroke failed to demonstrate the association of white matter disease with symptomatic ICH (sICH) after IVT (Capuana et al., 2020) . Recently, the role of enlarged perivascular spaces at basal ganglia (BG-EPVS) or centrum semiovale (CSO-EPVS) in ICH risk has been explored (Arba et al., 2017; Duperron et al., 2019) . To date, the effect of total cSVD burden on the risk of ICH after IVT in patients with acute ischemic stroke undergoing rt-PA thrombolysis remains unclear. Therefore, we aimed to investigate the relationship between total cSVD score and the risk of ICH after IVT in a retrospective observational study. This study is a retrospective observational study performed at two teaching hospitals (Union Hospital and Longyan First Hospital) of Fujian Medical University. We included the adult patients with confirmed acute ischemic stroke who were consecutively admitted between 1 January 2013 and 31 December 2020 based on the following criteria: (1) age > 18 years; (2) had no contraindication to IVT; and (3) receiving IVT using rt-PA within 4.5 h after the onset of symptoms. Patients were excluded if they met the following exclusion criteria: (1) did not undergo brain MRI detection 48 h within symptom onset; (2) insufficient MRI data due to poor image quality; and (3) missing data regarding baseline stroke severity assessed using the National Institutes of Health Stroke Scale (NIHSS) (Kwah and Diong, 2014) . Two authors blindly extracted the epidemiological, demographic, clinical, and laboratory data on admission and treatment data using a standardized data collection form. Hypertension was defined as documented systolic blood pressure ≥ 140 mmHg and/or diastolic blood pressure ≥ 90 mmHg, or by patients' self-reported diagnosis of hypertension and/or by the treatment of antihypertensives (Sacco et al., 2011; Du et al., 2020b) . Diabetes was defined as the current use of antidiabetic agents, having a fasting glucose level ≥ 7 mmol/L, having a non-fasting glucose level ≥ 11.1 mmol/L, or having a history of diabetes. Patients reporting current use of any type of tobacco during the past year were defined as current smokers (Sacco et al., 2011) . Those who consumed any dose of alcohol at least one time per week during the past year were defined as regular alcohol users (Du et al., 2020b) . We classified patients as having atrial fibrillation (AF) if they had been diagnosed with AF by electrocardiogram at admission or during the hospital stay or had a preadmission history (Emdin et al., 2017; Du et al., 2020b) . We assessed medical history of TIA, ischemic stroke, chronic heart failure, and ischemic heart disease by preadmission history and medical records. We obtained the information of prior use of medications as a chronic intake until the admission date. Eligible patients were treated with IVT using rt-PA (Boehringer Ingelheim Pharma GmbH) with and without bridging endovascular therapy based on the China stroke guideline (Chinese Society of Neurology and Chinese Stroke Society, 2018) and the best practice of the treating clinicians. The current China stroke guideline recommends rt-PA thrombolysis with a dose of 0.9 mg/kg. However, in clinical practice, some clinicians prefer to give a dose of 0.6 mg/kg to those with a higher ICH risk. Patients routinely underwent head computed tomography (CT) 24 h within IVT and 7 days after IVT to assess the occurrence of intracranial hemorrhage (ICH). Patients routinely underwent MRI within 48 h after onset to assess the location and size of infarction. Our main outcomes include any ICH, sICH, and remote parenchymal hemorrhages (rPH) after IVT, determined by a trained radiologist (JW) and a clinician (QY) using structured interviews blinded to baseline characteristics. sICH after IVT was defined as any subsequent visible hemorrhage associated with any neurological deterioration within 36 h of IVT, based on the National Institute of Neurological Disorders and Stroke (NINDS) criteria (Nisar et al., 2019) . rPH was defined as single or multiple ICHs without visible ischemic damage detected by cranial CT, remote from the infarction area (Prats-Sánchez et al., 2016) . Any disagreement was settled by referring to one senior clinician (NL) and radiologist (GC). All participants underwent brain MRI using 1. Glucose before AIS ‡ , mmol/L (IQR) 6.7 (5.6-9.5) 6.7 (5.9-8.4) Fasting blood-glucose § , mmol/L (IQR) 5.2 (4.5-6.6) 5.5 (4.9-6.6) Low density lipoprotein cholesterin # , (mmol/L), mean ± SD 2.68 ± 1.31 3.04 ± 1.00 Triglyceride ¶ , mmol/L (IQR) 1.2 (0.7-2.9) 1.2 (0.9-1.7) Cholesterol**, mmol/L, mean ± SD 4.37 ± 2.10 4.57 ± 1.12 SVD score (0-4), n (%) (DWI) sequence: TR, 3,570 ms; TE, 67 ms; FOV, 235 × 235 mm 2 . We determined well-established baseline imaging cSVD markers according to the previously published method (Wardlaw et al., 2013; Capuana et al., 2020; Du et al., 2020a) . Lacunes of presumed vascular origin were round or oval, located in the subcortex, similar to the signal of cerebrospinal fluid between 3 and 15 mm in diameter at fluid-attenuated inversion recovery (T2-FLAIR) and usually with low intensity of central cerebrospinal fluid and high-intensity margin around it (Wardlaw et al., 2013) . CMBs were a round or oval low-signal intensity focus of cerebral parenchyma on T2-weighted gradient-recalled echo (T2 * -GRE) or susceptibility-weighted imaging (SWI) with a diameter of 2-10 mm (Wardlaw et al., 2013; Capuana et al., 2020) . Enlarged perivascular spaces at basal ganglia (BG-EPVS) and centrum semiovale (CSO-EPVS) were defined on T2-weighted using a validated five-point ordinal scale as follows: 0 = no EPVS, 1 = 1-10 EPVS, 2 = 11-20 EPVS, 3 = 21-40 EPVS, and 4 = > 40 EPVS (Wardlaw et al., 2013; Du et al., 2020a) . We graded WMH in periventricular and deep white matter on T2-FLAIR sequences according to the Fazekas criteria (Staals et al., 2014) . Moderate to severe WMH was defined as a Fazekas score of 3 in periventricular white matter and/or 2-3 in deep white matter. Based on a well-validated method, we calculated a total cSVD score (0-4 points) to reflect the total cSVD burden, in which one score was defined as any of the following: ≥ 1 lacune; ≥ 1 CMB; ≥ 11 BG-EPVS; moderate to severe WMH (Staals et al., 2014; Du et al., 2021) . Categorical variables are summarized as absolute counts with percentages. Continuous variables were expressed as means (SD) if normally distributed or medians (interquartile range, IQR) if not normally distributed. We used the coefficient kappa test to evaluate the inter-rater and intra-rater reliabilities regarding MRI-visible cSVD assessment. The t-test or Wilcoxon test, where appropriate, were used to compare the difference in continuous variables, and the chi-square test or Fisher's exact test were used to compare the differences in categorical variables. We modeled the cSVD score on an ordinal scale continuously. Similar to several previous studies (Yilmaz et al., 2018; de Heus et al., 2020) , cSVD scores of 3 and 4 were combined due to a very small number of cases (n = 17 and n = 3, respectively). Multivariate logistic regression models were used to evaluate the association between the risk of any ICH after IVT and total cSVD score (0 vs. 1, 2, 3, 4). Other variables in the multivariable logistic regression analysis were first screened using univariable analysis with P < 0.2 and then using a penalized regression method known as the least absolute shrinkage and selection operator (LASSO) (Tibshirani, 1996) . An optimal value of the penalization parameter (lambda) for this model was obtained by 1-SE criteria using 10fold cross-validation (Tian et al., 2020) . We additionally performed an E-value analysis to fully explain away a specific exposure-outcome association. The E-value is defined as the minimum strength of association between an unmeasured confounder and an exposure/outcome to explain the observed association as non-causal (VanderWeele and Ding, 2017). We calculated the E-value using an online calculator. 1 The predictive performance of a combination of NIHSS and cSVD score for any ICH after IVT was evaluated by calculating the area under the receiver operating characteristic curve (AUROC) and its 95% confidence intervals (95% CI). All statistical analyses were performed using SPSS for Windows (version 25.0, IBM), STATA (version 15.0, Stata Corp., LP), and R version 4.5 (R Foundation, Vienna, Austria) with the "glmnet" package. 2 A P-value < 0.05 was considered statistically significant. A total of 785 patients with acute ischemic stroke underwent IVT between 1 January and 31 December 2020. Table 1 summarizes the prevalence of individual SVD markers. Regarding the total cSVD score, 81 (29.9%) patients had a score of 0, 99 (36.5%) patients had a score of 1, 71 (26.2%) patients had a score of 2, and 20 patients (7.4%) had a score of 3-4. Supplementary Table 2 shows that the intra-rater and inter-rater reliability for individual cSVD markers were satisfactory. Notably, 55 (20.3%) out of 271 patients experienced any ICH after IVT, 16 (5.9%) patients experienced a sICH, and 5 (1.9%) patients had rPH. (2) Moderate to severe white matter hyperintensities; (3) ≥ 1 CMBs; (4) CMB location (strict lobe, strict deep, and mixed); (5) cSVD score (0 vs. 1, 2, 3, 4); (6) NIHSS; (7) endovascular therapy; (8) atrial fibrillation; (9) age (per year increase); (10) current smoker; (11) hypertension; (12) ischemic heart disease; (13) prior use of antiplatelet treatment; (14) DWI mode 1 (single, scattered, and multiple); (15) DWI mode 2 (anterior circulation stroke, posterior circulation stroke, and mixed); (16) triglyceride; (17) ASCO stroke subtype. LASSO, least absolute shrinkage and selection operator; ASCO, atherothrombosis, small vessel disease, cardioembolism, other; NIHSS, National Institutes of Health Stroke Scale. , and a higher cSVD score (P < 0.001). We found similar ≥ 1 lacune presence between patients with and without any ICH after IVT [8 (14.5%) vs. 31 (14.3%)]. Table 2 shows that patients with rPH were more likely to have lacunes [3 (60.0%) vs. 36 (13.5%)], moderate to severe WMH [2 (40.0%) vs. 25 (9.3%)], and a higher cSVD score. There were no significant differences in individual SVD markers or total cSVD score among those with and without sICH ( Table 3 ). The LASSO model filtered out four variables (i.e., NIHSS, endovascular therapy, AF, and total cSVD score) from 17 risk characteristic individual factors for any ICH in the final multivariable regression model. Figure 1 shows how the regression coefficients change as a function of the penalization parameter [lambda (λ)]. On the far left, we saw the coefficients from standard logistic regression. As we go from left to right, the coefficients are "shrunk" toward zero. The coefficient estimated from the optimal model is indicated by the dashed line, which shows that only NIHSS, endovascular therapy, AF, and the total cSVD score were kept in the model. Multivariable logistic regression analysis showed that the total cSVD score (per unit increase, OR 2.03 95% 1.22-3.41, P = 0.007) was significantly associated with a higher risk of any ICH after IVT. The risk of any ICH increased with increasing cSVD score (cSVD score 1, OR 1.29, 95% CI 0.44-3.94, cSVD score 2, OR 4.90, 95% CI 1. OR 12.40, , compared with a score of 0, P < 0.001, AUC 0.698, 95% CI 0.620-0.777, FIGURE 2 | Bias plots assessed by E-value analysis (outcome is any ICH). Bias plots represent varying levels of unobserved confounder/exposure association and unobserved confounder/outcome association necessary to render observed exposure/outcome association null for total cSVD score per point increase. cSVD, cerebral small vessel disease; NIHSS, National Institutes of Health Stroke Scale; ROC, receiver operating characteristic curve. *P = 0.0004 (comparison of AUROC of NIHSS score and a combined NIHSS score and total SVD score using Z statistics). FIGURE 3 | Predictive performances of NIHSS score, SVD burden and a combined NIHSS score, and SVD burden using ROC analysis. ROC, receiver operating characteristic; cSVD, cerebral small vessel disease. Table 4 ). The very small number of sICH (n = 16) and rPH events (n = 5) did not allow us to perform the logistic regression analysis using the LASSO approach. Sensitivity analyses using the E-value revealed that it would need moderately robust unobserved confounding to render the exposure-outcome (cSVD-any ICH) association null (Figure 2 ). Based on a previous meta-analysis and our LASSO selection results, the NIHSS score was considered a robust predictor for any ICH after IVT (Whiteley et al., 2012) . Therefore, we compared the predictive performance of a combination of NIHSS and cSVD score with that of the NIHSS score alone. The results showed that combining cSVD and NIHSS score had a better predictive performance than NIHSS score alone (AUC 0.811, 95% CI 0.756-0.866 vs. 0.784, 95% CI 0.723-0.846, P = 0.0004, Table 5 and Figure 3 ). This study shows that a higher SVD burden is associated with increased risk of any ICH after IVT, suggesting that the risk of any ICH in patients with acute ischemic stroke who received IVT with rt-PA might be related to SVD severity. Moreover, the cSVD score adds to the NIHSS score alone in predicting any ICH after IVT. In our cohort, patients with any ICH were more likely to have ≥ 1 CMB presence. However, CMB presence was not screened out using the LASSO model. The limited sample size of our study does not enable us to perform the further analysis stratified by CMB burden. Our findings were supported by several studies that showed the CMB presence was not associated with an increased risk of any ICH (Kakuda et al., 2005; Kim et al., 2006; Fiehler et al., 2007) or sICH (Charidimou et al., 2017; Capuana et al., 2020) after IVT. Current stroke guidelines do not recommend IVT in patients with > 10 CMBs at baseline (Chinese Society of Neurology and Chinese Stroke Society, 2018; Powers et al., 2018; Berge et al., 2021) . A recent study using a multistep algorithm populated with data showed that in patients with > 10 CMBs, IVT is associated with higher mortality and net harm defined by age, stroke symptom severity, and treatment delay (Schlemm et al., 2020) . To date, there are no studies reporting directly on the risk of sICH in patients without IVT stratified by CMB burden. Future prospective studies are needed to define the CMB burden threshold in individual risk stratification predicting the bleeding risk after IVT and functional outcome. Evaluating the total cSVD burden may allow more accurate estimates of size effects on ICH risk after IVT than individual cSVD markers. In our cohort, the total SVD score showed a dose-effect response, reinforcing the validity of the findings. Similarly, a small sample-sized retrospective study showed that severe cSVD burden (cSVD score of 3-4) was associated with an increased risk of hemorrhagic transformation after IVT (OR = 8.429, 95% CI 1.643-43.227, P = 0.011) (Liu et al., 2019) . Our study adds to the aforementioned study by investigating the association between total cSVD score (0 vs. 1, 2, 3, 4) and the bleeding risk after IVT, using variables that were selected by LASSO. Moreover, our data supported the possible relationship between the risk of rPH and the total cSVD score. However, the small number of rPH events did not enable us to perform the regression analysis. A previous single-center observational study showed that patients with stroke with r-PH after IVT were more likely to have strictly lobar CMBs, suggesting amyloid angiopathy (Drelon et al., 2020) . Since the occurrence of r-PH may be related to worse outcomes, the association between r-PH and total cSVD score needs further investigation. Possible mechanisms that account for the role of cSVD in ICH after IVT include endothelial dysfunction, arterial stiffening (arterial stiffening might increase ICH risk), and increased blood-brain barrier permeability (Campbell and Khatri, 2020) . Consistent with a previous meta-analysis (Whiteley et al., 2012) , our data support that the NIHSS score is a robust predictor for ICH risk after IVT. The larger AUC of a combined NIHSS and cSVD score than an NIHSS score alone suggests that adding cSVD as neuroimaging biomarkers to a widely validated factor (NIHSS) might improve specificity and sensitivity in predicting ICH after IVT. Since most of our cohort had mild to moderate acute ischemic stroke [median NIHSS score 5 (IQR 3-10)], our results might not be generalized to the patients with severe stroke. Future large-scale studies are required to refine and validate robust risk prediction tools. There are some limitations of this study. First, this is a retrospective study with a limited sample size. Second, excluding patients who did not undergo MRI within 48 h after onset or insufficient MR image might introduce selection bias since most severe patients were not included. Third, we did not analyze the functional outcome data after discharge in this study; a longer follow-up period data will better determine the significance of our findings. Fourth, very few patients underwent MR before IVT; therefore, we could not exclude the new cSVD markers (particularly new CMBs) after IVT (Braemswig et al., 2019) . However, the cSVD findings did not influence the IVT decisions. Fifth, due to the small number of outcome events, we did not assess the predicting performance of cSVD for sICH and rPH, which is the most clinically relevant. Finally, due to the data unavailability, we were unable to calculate the estimated effect size adjusted for some unmeasured confounders. However, we used the E-value sensitivity analysis to quantify the potential implications of unmeasured confounders and found that it was unlikely to explain the entire outcome. This study shows that the total cSVD score is associated with a higher risk of any ICH after IVT. However, our study was not powered to assess the association between cSVD score and sICH risk. Further studies are warranted to determine whether rt-PA can be safely initiated in patients with acute ischemic stroke who have coexistent cSVD. The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s. The studies involving human participants were reviewed and approved by the Fujian Medical University Union Hospital Ethics Committees (2019KY076). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements. HD and NL: concept, design, and full access to all of the data in this study and take responsibility for the integrity of the data and the accuracy of the data analysis. HD, SW, HhL, and NL: acquisition, analysis, and interpretation of data. HD, SW, and HfL: drafting of the manuscript. GA, DW, HfL, YC, DC, QY, SF, RC, YZ, JW, GC, and JC: critical revision of the manuscript for important intellectual content. HD, HhL, and GA: statistical analysis. All authors contributed to the article and approved the submitted version. This study was supported by the Fujian Provincial Special Foundation for Natural Science Innovation Project (2016B014). Intensive blood pressure reduction with intravenous thrombolysis therapy for acute ischaemic stroke (enchanted): an international, randomised, open-label, blinded-endpoint, phase 3 trial Small vessel disease and clinical outcomes after iv rt-pa treatment European stroke organisation (eso) guidelines on intravenous thrombolysis for acute ischaemic stroke Predictors of new remote cerebral microbleeds after iv thrombolysis for ischemic stroke Hemorrhagic risk after intravenous thrombolysis for ischemic stroke in patients with cerebral microbleeds and white matter disease Leukoaraiosis, cerebral hemorrhage, and outcome after intravenous thrombolysis for acute ischemic stroke Microbleeds, cerebral hemorrhage, and functional outcome after stroke thrombolysis Chinese guidelines for diagnosis and treatment of acute ischemic stroke 2018. CHIN Day-to-day home blood pressure variability is associated with cerebral small vessel disease burden in a memory clinic population Predictors for intracranial hemorrhage following intravenous thrombolysis in posterior circulation stroke Remote brain hemorrhage after iv thrombolysis The effect of vascular risk factor burden on the severity of COVID-19 illness, a retrospective cohort study Association between steno-occlusive middle cerebral artery and basal ganglia perivascular spaces Small vessel disease and ischemic stroke risk during anticoagulation for atrial fibrillation after cerebral ischemia High dilated perivascular space burden: a new mri marker for risk of intracerebral hemorrhage Usual blood pressure, atrial fibrillation and vascular risk: evidence from 4.3 million adults Bleeding risk analysis in stroke imaging before thrombolysis (brasil): pooled analysis of t2 * -weighted magnetic resonance imaging data from 570 patients Thrombolysis with alteplase 3 to 4.5 hours after acute ischemic stroke Clinical importance of microbleeds in patients receiving iv thrombolysis Multiple cerebral microbleeds in hyperacute ischemic stroke: impact on prevalence and severity of early hemorrhagic transformation after thrombolytic treatment National institutes of health stroke scale (nihss) Relationship between total cerebral small vessel disease burden and hemorrhagic transformation of acute ischemic stroke patients after intravenous thrombolysis. CHIN Predicting the risk of symptomatic intracerebral hemorrhage in ischemic stroke treated with intravenous alteplase Symptomatic intracerebral hemorrhage after intravenous thrombolysis: predictive factors and validation of prediction models Thrombolysis in acute stroke patients with cerebral small vessel disease 2018 guidelines for the early management of patients with acute ischemic stroke: a guideline for healthcare professionals from the american heart association/american stroke association Remote intracerebral hemorrhage after intravenous thrombolysis Effect of prior medical treatments on ischemic stroke severity and outcome Benefit of intravenous thrombolysis in acute ischemic stroke patients with high cerebral microbleed burden Stroke subtype, vascular risk factors, and total mri brain small-vessel disease burden Prediction of overall survival in resectable intrahepatic cholangiocarcinoma: Isicc-applied prediction model Regression shrinkage and selection via the lasso Sensitivity analysis in observational research: introducing the e-value Neuroimaging standards for research into small vessel disease and its contribution to ageing and neurodegeneration Risk factors for intracranial hemorrhage in acute ischemic stroke patients treated with recombinant tissue plasminogen activator: a systematic review and meta-analysis of 55 studies Practical small vessel disease score relates to stroke, dementia, and death The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi. 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.Copyright © 2022 Du, Wu, Lei, Ambler, Werring, Li, Chen, Chen, Yuan, Fang, Chen, Zhang, Wei, Chen, Chen 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.