Skip to main content

ORIGINAL RESEARCH article

Front. Neurosci., 02 March 2022
Sec. Brain Imaging Methods

Intravoxel Incoherent Motion Magnetic Resonance Imaging Used in Preoperative Screening of High-Risk Patients With Moyamoya Disease Who May Develop Postoperative Cerebral Hyperperfusion Syndrome

\r\nFeng Gao*&#x;Feng Gao1*†Wei Zhao&#x;Wei Zhao2†Yu ZhengYu Zheng3Yu DuanYu Duan4Ming JiMing Ji1Guangwu LinGuangwu Lin1Zhenfang Zhu*Zhenfang Zhu1*
  • 1Department of Radiology, Huadong Hospital Fudan University, Shanghai, China
  • 2Department of Radiology, The Second Xiangya Hospital of Central South University, Changsha, China
  • 3Department of Radiology, Chengdu Second People’s Hospital, Chengdu, China
  • 4Department of Neurosurgery, Huadong Hospital Fudan University, Shanghai, China

Objective: This study aimed to investigate the feasibility of preoperative intravoxel incoherent motion (IVIM) MRI for the screening of high-risk patients with moyamoya disease (MMD) who may develop postoperative cerebral hyperperfusion syndrome (CHS).

Methods: This study composed of two parts. In the first part 24 MMD patients and 24 control volunteers were enrolled. IVIM-MRI was performed. The relative pseudo-diffusion coefficient, perfusion fraction, apparent diffusion coefficient, and diffusion coefficient (rD*, rf, rADC, and rD) values of the IVIM sequence were compared according to hemispheres between MMD patient and healthy control groups. In the second part, 98 adult patients (124 operated hemispheres) with MMD who underwent surgery were included. Preoperative IVIM-MRI was performed. The rD*, rf, rADC, rD, and rfD* values of the IVIM sequence were calculated and analyzed. Operated hemispheres were divided into CHS and non-CHS groups. Patients’ age, sex, Matsushima type, Suzuki stage, and IVIM-MRI examination results were compared between CHS and non-CHS groups.

Results: Only the rf value was significantly higher in the healthy control group than in the MMD group (P < 0.05). Out of 124 operated hemispheres, 27 were assigned to the CHS group. Patients with clinical presentation of Matsushima types I–V were more likely to develop CHS after surgery (P < 0.05). The rf values of the ipsilateral hemisphere were significantly higher in the CHS group than in the non-CHS group (P < 0.05). The rfD* values of the ACA and MCA supply areas of the ipsilateral hemisphere were significantly higher in the CHS group than in the non-CHS group (P < 0.05). Only the rf value of the anterior cerebral artery supply area in the contralateral hemisphere was higher in the CHS group than in the non-CHS group (P < 0.05). The rf values of the middle and posterior cerebral artery supply areas and the rD, rD*, and rADC values of the both hemispheres were not significantly different between the CHS and non-CHS groups (P > 0.05).

Conclusion: Preoperative non-invasive IVIM-MRI analysis, particularly the f-value of the ipsilateral hemisphere, may be helpful in predicting CHS in adult patients with MMD after surgery. MMD patients with ischemic onset symptoms are more likely to develop CHS after surgery.

Introduction

Moyamoya disease (MMD) is a common neurological disease characterized by progressive steno-occlusion of the distal internal carotid artery, leading to the compensatory development of collateral vessels (moyamoya vessels) at the base of the brain (Suzuki and Takaku, 1969). MMD is more prevalent in the Asian population and more frequently occurs in women. Revascularization is recommended for patients with MMD (Research Committee on the Pathology and Treatment of Spontaneous Occlusion of the Circle of Willis, and Health Labour Sciences Research Grant for Research on Measures for Infractable Diseases, 2012). Superficial temporal artery-middle cerebral artery (STA-MCA) bypass in combination with indirect procedures has been employed as treatment for patients with MMD (Kazumata et al., 2014). Although an STA-MCA bypass carries a relatively small amount of blood, cerebral hyperperfusion is frequently observed after revascularization surgery, particularly in adult patients (Zhao et al., 2013; Kazumata et al., 2014). Furthermore, patients with MMD are more likely to develop cerebral hyperperfusion syndrome (CHS) after revascularization surgery than those with other atherosclerotic occlusive cerebrovascular diseases (Fujimura et al., 2011). Postoperative complications in patients with MMD can cause not only postoperative neurological symptoms but also permanent neurological defects (Cho et al., 2014; Antonucci et al., 2016; Kim et al., 2016). Postoperative CHS in MMD can lead to an increase in intracranial vascular pressure and may induce intracranial hemorrhage and even hematogenous cerebral edema (Mesiwala et al., 2008; Park et al., 2018). Studies have shown that appropriately lowering the patients’ blood pressure during the perioperative period may reduce the incidence of postoperative CHS (Hayashi et al., 2012). In clinical practice, lowering the blood pressure of patients with MMD after surgery may also lead to cerebral vascular insufficiency and cerebral infarction (Lee et al., 2018). Identifying patients with MMD who are at high risk for CHS development after surgery and implementing appropriate perioperative management are important to prevent the occurrence of postoperative hyperperfusion and improve the prognosis of these patients (Ostergaard et al., 2014).

Impaired cerebral autoregulation has been implicated in the mechanism of cerebral hyperperfusion mediated by endothelial dysfunction. Various imaging modalities such as single-photon emission computed tomography (SPECT), computed tomography (CT), and magnetic resonance imaging (MRI) can determine the cerebral perfusion status. SPECT and CT perfusion are radiative, whereas MRI techniques, including dynamic contrast-enhanced MRI, dynamic susceptibility contrast-enhanced (DSC) MRI, arterial spin labeling (ASL) MRI, and methods used to analyze cerebral microcirculation, such as intravoxel incoherent motion (IVIM) MRI, and vascular space occupancy are non-radiative (Callewaert et al., 2021). Brain IVIM technique is advantageous over other technique, such as positron emission tomography, CT, DSC MRI, and ASL MRI, because it does not require contrast media injection or radiation exposure and is not affected by the arterial input (Hara et al., 2019). Microvessels are below the currently achievable spatial resolution of MRI and other clinically applicable imaging methods, making the direct non-invasive measurement of microvascular density impossible to be performed. Nonetheless, MRI can be applied to measure microvascular perfusion, providing indirect readouts of the underlying microvasculature alteration (Federau et al., 2012; Callewaert et al., 2021).

IVIM sequence, which was proposed by Le Bihan et al. (1986), is a diffusion MRI technique that can simultaneously provide information regarding both molecular diffusion and perfusion (Le Bihan, 2019). Microscopic movement in organisms includes the diffusion of water molecules and microcirculation of blood. IVIM is a double exponential model that divides biological tissues into two components: a slow-moving component in which water molecules diffuse based on Brownian motion and a fast-moving component in which water molecules move due to blood circulation (Federau et al., 2012). Four parameters, including apparent diffusion coefficient (ADC), diffusion coefficient (D), pseudo-diffusion coefficient (D*), and perfusion fraction (f), can be can be obtained from IVIM sequence. Promising results had been obtained when perfusion was measured using the IVIM technique in various diseases of different body organs, such as liver (Guiu et al., 2012), kidney (Thoeny et al., 2009), pancreas (Lemke et al., 2009), breast (Sigmund et al., 2011), prostate (Riches et al., 2009), and head and neck tumors (Sumi and Nakamura, 2013). The application of IVIM technique to the brain has increased over recent years (Paschoal et al., 2018). Previous study showed that f and fD* are more sensitive to detect changes during the stroke process, as these parameters are more related to cerebral blood flow (CBF) (Yao et al., 2016). IVIM also can be used to assessment of patients with cerebral small vessel disease (Wardlaw et al., 2013). Hara et al. (2019) used IVIM to evaluate the hemodynamic disturbance of MMD by comparison with the gold-standard 15O-gas positron emission tomography method. All those studies showed that IVIM may be used to non-invasively assess cerebral hemodynamic impairment in patients with MMD. To the best of our knowledge, thus far, no study has used preoperative IVIM-MRI to predict patients with MMD who may develop postoperative CHS. This study aimed to investigate potential risk factors for CHS after surgery in patients with MMD using the IVIM technique.

Materials and Methods

Clinical Data and Participants

This study was approved by the ethics committee of our hospital. Written informed consent was obtained from all patients. All experiments were performed in accordance with the relevant guidelines and regulations set by the ethics committee. This study was composed of two parts. First part compared parameters of MMD patient and healthy control groups. Second part only included MMD patients, preoperative parameters between CHS and non-CHS group were compared.

From January 2018 to July 2018, 24 adult patients with MMD (mean age, 41.67 ± 10.606; range, 19–65 years; female, 12; male, 12) with 48 brain hemispheres and 24 healthy adult volunteers (mean age, 43.67 ± 8.499; range, 19–65 years; female, 14; male, 10) with 48 brain hemispheres were included in this study. MMD was diagnosed using digital substraction angiography (DSA) according to the criteria of the Research Committee on Spontaneous Occlusion of the Circle of Willis (MMD) of the Ministry of Health and Welfare, Japan (Fukui, 1997). The inclusion criteria for the MMD group were as follows: (i) patients aged > 18 years and (ii) patients with no contraindications to MR examination. All healthy volunteers were from our hospital and underwent routine physical examinations. The inclusion criteria for the healthy control group were as follows: (i) patients aged > 18 years; (ii) no previous neurological history; (iii) MRI examination revealing no obvious hemorrhagic foci or encephalomalacia; (iv) no autoimmune disease; (v) magnetic resonance angiography examination showing no obvious stenosis and occlusion of cerebral arteries; and (vi) patients with no contraindications to MRI.

From August 2018 to May 2021, 98 adult patients with MMD (mean age, 42.87 ± 11.32 years; range, 19–65 years; female, 52; male, 46) with 124 operated hemispheres were included in this study (26 patients underwent bilateral revascularization surgeries, 72 patients underwent unilateral revascularization surgeries. Patients with bilateral cerebral hemispheres that both meet surgical indications, the side with more severe clinical symptoms or more worse perfusion status will be operated first.). All MR examinations were performed within 1 week before the revascularization surgery. MMD was diagnosed using DSA (Fukui, 1997). The inclusion criteria were as follows: (i) patients aged > 18 years; (ii) patients with no contraindications to MR examination; and (iii) patients who underwent STA-MCA bypass and encephalo-duro-myo-synangiosis (EDMS) combined surgery. In our hospital, most patients underwent combined surgery, although only a few patients underwent EDMS surgery. Therefore, we only included patients who underwent combined surgery in this study to avoid bias caused by different surgery types.

Magnetic Resonance Imaging Examination Protocol and Image Analysis

All preoperative MR examinations were performed with a 3-T whole-body MRI scanner (Siemens Skyra Freedom; Siemens Medical Solutions, Erlangen, Germany) using a head-neck coil. All patients were instructed to stay still and not to think about anything with their eyes closed during the MRI examination. T1-weighted imaging, T2-weighted imaging, fluid-attenuated inversion-recovery sequence, time-of-flight magnetic resonance angiography, and IVIM-MRI were performed in this study.

Intravoxel Incoherent Motion Scanning Protocol and Imaging Analysis

Multiple b-values were used (0, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 900, and 1,000 s/mm2) in the IVIM sequence. Detailed scanning parameters are listed in Table 1. Four parameters, including apparent diffusion coefficient (ADC), diffusion coefficient (D), pseudo-diffusion coefficient (D*), and perfusion fraction (f), were calculated and derived using a Siemens IVIM work-in-progress package (MR Body Diffusion Toolbox version 1.3.0) (Figure 1). Regions of interest (ROIs) with dimensions of approximately 200 mm2 were placed manually to avoid subarachnoid spaces, cerebral infarction, and cerebral hemorrhage areas. As ROIs in this study were placed manually, to ensure the accuracy of measurement, the ROIs were place three times for each patient, and the average value was considered as the final result. Each cerebral hemisphere was divided into three areas according to the different supply cerebral arteries [anterior cerebral artery (ACA) supply area, MCA supply area, and posterior cerebral artery (PCA) supply area]. A total of 22 ROIs were placed symmetrically on the peripheral white matter, basal ganglia, and cerebellar hemisphere levels (Togao et al., 2006; Zaharchuk et al., 2011). Moreover, the ipsilateral cerebellar hemisphere was used as a reference to standardize the evaluated parameters. Subsequently, rADC, rD, rD*, and rf were obtained [i.e., rADC = ADC (different blood artery supply area)/ADC (ipsilateral cerebellar hemisphere) (Figure 2)].

TABLE 1
www.frontiersin.org

Table 1. MRI scanning sequences and parameters.

FIGURE 1
www.frontiersin.org

Figure 1. Intravoxel incoherent motion MRI imaging obtained after post-processing (the columns show three different levels, peripheral white matter level, basal ganglia level, cerebellar hemisphere level; the line shows four different IVIM parameters, ADC, D, f, D*value). (A–C) Image of ADC value; (D–F) image of D-value; (G–I) image of f value; (J–L) image of D* value.

FIGURE 2
www.frontiersin.org

Figure 2. The placement of ROIs on IVIM images in three different level. (A) peripheral white matter level; (B) basal ganglia level; (C) cerebellar hemisphere level. A total of 22 ROIs were placed symmetrically on the three transverse images.

Diagnosis of Cerebral Hyperperfusion Syndrome and Group Analysis

In the first part, the participants were divided into an MMD patient group and a healthy control volunteer group. All parameters were compared between the two groups. The rD*, rf, rADC, and rD values of the IVIM sequence were compared according to hemisphere between the MMD and healthy control groups. The D*, f, ADC, and D values were compared according to cerebellum between the MMD and healthy control groups.

In the second part, symptomatic CHS was defined as the postoperative development of a severe headache, seizures, and new neurological deficits, with neither definite hematomas nor definite acute infarctions presenting or observed on brain CT images, diffusion magnetic resonance images, or both (van Mook et al., 2005; Iwata et al., 2014). A senior neurosurgeon blinded to the study confirmed the diagnosis of postoperative CHS. Operated hemispheres were divided into CHS and non-CHS groups based on whether the patient developed CHS after surgery. We then compared the parameters between the CHS and non-CHS groups, and each hemisphere was divided according to different supply areas: the ACA, MCA, and PCA supply areas. In this section, the rD*, rf, rADC, rD, and rfD* values of the IVIM sequence were compared according to different cerebral artery supply areas in the two hemispheres (ipsilateral/surgery side and contralateral/non-surgery side).

Statistical Analyses

The independent t-test or Mann–Whitney U-test, where appropriate, was used to assess between-group differences in continuous variables (age as well as the rD*, rf, rADC, and rD values of the IVIM sequence). The chi-square test was applied to evaluate between-group differences in categorical variables [sex, hypertension, diabetes, Suzuki stage (Suzuki and Takaku, 1969), and Matsushima type (Matsushima et al., 1990, Supplementary Table 1). The IVIM parameters showed statistically significant differences between the CHS and non-CHS groups; the receiver operating characteristic (ROC) curve was analyzed to obtain the diagnostic efficacy of these parameters. Statistical results were considered significant when the P-values were < 0.05. Statistical analyses were performed using a commercially available computer software program (IBM SPSS version 22).

Results

The rD, rD*, rADC, and rf values of the MMD and healthy control groups are shown in Table 2. Only the rf value was significantly higher in the healthy control group than in the MMD group (P < 0.05). No significant differences were found in rADC, rD, and rD* between the MMD and healthy control groups (P > 0.05). No significant differences were found in ADC, f, D, and D* of cerebellum between the MMD and healthy control groups (P > 0.05).

TABLE 2
www.frontiersin.org

Table 2. Examination results of the MMD and control groups.

There were 27 (27/124, 21.77%) operated hemispheres developed CHS after revascularization surgery. The main symptoms of CHS included severe headache (six cases), persistent vomiting (two cases), dysphagia (three cases), speech impairment (seven cases), face and eye pain (one case), muscle strength decrease (the upper limbs in two cases and the lower limbs in four cases), and seizure (two cases). These symptoms were observed 1–9 days (mean ± standard deviation: 3.444 ± 2.532 days) after combined surgery. All symptoms improved before discharge from the hospital.

The CHS and non-CHS groups exhibited a significant difference with respect to the Matsushima type (P < 0.05). Overall, 88.89% (24/27 operated hemispheres) of operative patients in the CHS group presented with Matsushima types I–V, whereas 58.76% (57/97 operated hemispheres) of operative patients in the non-CHS group presented with Matsushima types I–V. No statistically significant differences in age, sex, hypertension, diabetes, and Suzuki stage were observed between the CHS and non-CHS groups (P > 0.05; Table 3).

TABLE 3
www.frontiersin.org

Table 3. Clinical characteristics of the patients in the CHS and non-CHS groups.

The preoperative microvascular perfusion status of the CHS and non-CHS groups is shown in Table 4. The hemispheres were divided into ipsilateral (hemispheres performed with surgery) and contralateral hemispheres. The rf values of the ACA, MCA, and PCA supply areas of the ipsilateral hemisphere were significantly higher in the CHS group than in the non-CHS group (P < 0.05). Furthermore, rADC, rD, and rD* in all areas of the ipsilateral hemisphere showed no significant differences between the CHS and non-CHS groups (P > 0.05). Only the rf value of the ACA supply area in the contralateral hemisphere was significantly higher in the CHS group than in the non-CHS group (P < 0.05). The rf values of the MCA and PCA supply areas of the contralateral hemisphere were not significantly different between the CHS and non-CHS groups (P > 0.05). The values of rADC, rD, and rD* in all areas of the contralateral hemisphere were not significantly different between the CHS and non-CHS groups (P > 0.05). The rfD* values of the ACA and MCA supply areas of the ipsilateral hemisphere were significantly higher in the CHS group than in the non-CHS group (P < 0.05). Furthermore, rfD* value of PCA supply areas of ipsilateral hemisphere and in all areas of the contralateral hemisphere showed no significant differences between the CHS and non-CHS groups (P > 0.05).

TABLE 4
www.frontiersin.org

Table 4. Preoperative microvascular perfusion status of the CHS and non-CHS groups.

The parameters of the IVIM sequence, which showed statistically significant differences between the CHS and non-CHS groups, were used to generate ROC curves for the prediction of postoperative CHS. The results showed that the areas under the curve (AUCs) of rf in the ACA, MCA, and PCA supply areas of the ipsilateral hemisphere for predicting postoperative CHS were 0.700, 0.733, and 0.701, respectively. The AUC for rf in the ACA supply area in the contralateral hemisphere was 0.625 for predicting postoperative CHS. The AUC for rfD* in the ACA and MCA supply area in the ipsilateral hemisphere for predicting postoperative CHS were 0.654 and 0.680. The threshold, sensitivity, and specificity of each diagnostic indicator are shown in Table 5 and Figure 3.

TABLE 5
www.frontiersin.org

Table 5. Sensitivity and specificity of diagnostic thresholds for different diagnostic indicators.

FIGURE 3
www.frontiersin.org

Figure 3. ROC curve of rfD* in the ACA, MCA supply areas of the ipsilateral hemisphere, rf in the ACA, MCA, PCA supply areas of the ipsilateral hemisphere and in the ACA supply area of the contralateral hemisphere for predicting postoperative CHS. The AUCs of rf in the ACA, MCA, and PCA supply areas of the ipsilateral hemisphere for predicting postoperative CHS were 0.700, 0.733, and 0.701, respectively. The AUC for rf in the ACA supply area in the contralateral hemisphere was 0.625 for predicting postoperative CHS. The AUC for rfD* in the ACA and MCA supply area in the ipsilateral hemisphere for predicting postoperative CHS were 0.654 and 0.680.

Discussion

Our study used the IVIM technique to compare the microvascular perfusion status between patients with MMD and healthy volunteers and subsequently compared the preoperative microvascular perfusion status between CHS and non-CHS groups with MMD. The results indicated that the rf value was higher in the healthy control group than in the MMD group. Preoperative non-invasive IVIM-MRI analysis, particularly the f-value of the ipsilateral hemisphere, may be helpful in predicting CHS in patients with MMD after surgery. Patients with MMD with Matsushima types I–V are more likely to develop CHS after surgery.

In the IVIM sequence, the number of b-values and specific values are also related to the accuracy of the IVIM measurement parameters. In our study, referring to previous study results, a total of 16 b-values were set in the IVIM sequence (Federau et al., 2012). The b-values were 0, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 900, and 1,000 s/mm2, among which there were eight b-values lower than 200 s/mm2. The base for separate perfusion and diffusion parameters in IVIM is the signal-to-noise ratio (SNR). Wu et al. (2015) suggested that when b = 1,000 s/mm2, the minimum SNR is 30 in the brain and can obtain a reliable f-value to measure cerebral blood volume (CBV). In this study, the maximum b-value was 1,000 s/mm2.

D is the diffusion coefficient of free water, and ADC is the sum of the contributions of all diffusion coefficients related to the resulting motion (Le Bihan et al., 1986). Although D and ADC are different in mathematical theory, they have a strong correlation. The D-value, along with the ADC, is widely used in the diagnosis of acute stroke and has shown reliable diagnostic efficacy (Suo et al., 2016). In this study, we chose subjects without acute cerebral infarction to exclude different blood perfusion statuses that may influence the D and ADC values. The results showed no significant difference in rD and rADC values between the MMD patient and healthy control groups (P > 0.05). The IVIM microvascular perfusion parameters comprise a perfusion fraction (f), which describes the fraction of incoherent signal that arises from the vascular compartment, and a D*, which macroscopically describes the incoherent blood movement in the microvasculature compartment. A previous study showed that f represents the ratio of capillary volume in the voxel to the entire tissue volume, which can reflect the proportion of water molecules in the capillary network (Lee et al., 2014). In other words, f is related to the microvascular perfusion status of the tissues and is related to the capillary density score (Lee et al., 2014). Previous study showed that f is related to CBV, D* is related to 1/MTT, and fD* is related to CBF (Paschoal et al., 2018). In our study, the results showed that only the rf in the healthy control group was significantly higher than that in the MMD patient group (P < 0.05). This finding is also consistent with the hemodynamic change in MMD, which is characterized by varying degrees of steno-occlusion of the distal internal carotid artery and reduction of cerebral microvascular perfusion. The D* value is also related to the microvascular perfusion of the capillary network, although our study showed no statistically significant differences in rD* between the MMD and healthy control groups (P > 0.05). This result may be related to the stability of the D* value. A previous study showed that the repeatability of the D* value is poor, which may be due to the high sensitivity of this value; the D* value may be influenced by the partial volume effect caused by capillary blood flow and lacunae, which contain the cerebrospinal fluid (CSF) (Bisdas et al., 2013). In this study, patients with MMD with bilateral vascular involvement were enrolled; therefore, we chose a relatively normal ipsilateral cerebellar hemisphere rather than the other hemispheres to standardize the D, D*, ADC, and f values. The value of ADC, f, D, and D* value in cerebellum between MMD and healthy control groups were also compared, in order to evaluate the influence of crossed cerebellar diaschisis (CCD). The result showed that no significant differences were found of ADC, f, D, and D* value in cerebellum between the MMD and healthy control groups (P > 0.05, Table 2). However, MMD may also affect the PCAs, which may also cause some deviation in the measurements.

In patients with MMD, chronic ischemia leads to the development of new pathological vessels, with impaired cerebrovascular autonomic regulation and cerebrovascular reactivity, and the vessels are unable to control the increased cerebral blood flow after surgery. Therefore, CHS is considered to occur after surgery (Hayashi et al., 2012). The incidence of CHS in our study was 21.77% (27/124), which is consistent with a previous study result ranging from 6.7 to 38.2% (Fujimura et al., 2011, 2012). Although patient age, sex, hypertension, diabetes, and Suzuki stage showed no significant differences between the CHS and non-CHS groups in this study, the Matsushima type had a statistically significant difference between the two groups. A previous study demonstrated that adult patients with MMD commonly present with hemorrhage due to intimal tearing and subsequent blood vessel rupture (Dusick et al., 2011; Okada et al., 2012). Our study showed that approximately 34.68% (43/124 operated hemisphere) of the operative patients presented with Matsushima type VI (hemorrhage) symptoms, although among operative patients with postoperative CHS, approximately 88.89% (24/27 operated hemisphere) presented with Matsushima types I–V (ischemic) symptoms and 11.11% (3/27) presented with Matsushima type VI (hemorrhage) symptoms. The results were consistent with those of a previous study, which concluded that symptomatic hyperperfusion occurred more often in patients with ischemic onset than in those with hemorrhage (Hayashi et al., 2010).

Previous studies mainly used perfusion imaging to evaluate postoperative CHS and demonstrated that either an increase of more than 100% or a comparable postoperative increase in the regional cerebral blood flow after revascularization in patients with MMD may indicate CHS (Kaku et al., 2012; Fujimura et al., 2015). However, these radiologic techniques can only measure hemodynamics in a specific ROI. In this study, we measured the microvascular perfusion of both ipsilateral and contralateral surgery hemispheres according to different cerebral artery supply areas.

A previous study showed that although IVIM perfusion might be affected by physiological factors, such as the cardiac cycle and non-vascular components, there was a significant correlation between f and CBV, fD*, and CBF (Federau et al., 2014; Wu et al., 2015). Our study showed that in the CHS group, the rf values of the ACA, MCA, and PCA supply areas in the ipsilateral surgery hemisphere were significantly higher than those in the non-CHS group (P < 0.05). As mentioned above, f is related to the microvascular perfusion status of the tissues and is related to the capillary density score (Lee et al., 2014). Higher f-value may indicate a higher proportion of water molecules in the capillary network and may also indicate the higher density of micro collateral vessels. Consequently, areas with a higher density of micro collateral vessels may receive more blood volume after revascularization surgery. Our study also showed that the rfD* values of the ACA and MCA supply areas of the ipsilateral hemisphere were significantly higher in the CHS group than in the non-CHS group (P < 0.05). This may also reflect that in the CHS group, the preoperative CBF of ipsilateral hemisphere was higher than non-CHS group. Studies have shown that if revascularization surgery is performed when the cerebrovascular reserve is normal, surgery may increase the risk of postoperative CHS and may also affect the establishment of collateral circulation (Han et al., 2011). The results of this study are also consistent with those of previous studies; patients with MMD with higher preoperative f-values are more likely to develop postoperative CHS than patients with lower f-values. This finding also indicates that patients with MMD who develop postoperative CHS have better preoperative cerebral microvascular perfusion status than patients without postoperative CHS. The f-value for the ACA supply area in the contralateral hemisphere was significantly higher in the CHS group than in the non-CHS group (P < 0.05). The rf values of the MCA and PCA supply areas were not significantly different between the CHS and non-CHS groups (P > 0.05). We hypothesized that it might be related to the anterior communicating artery, which connects the bilateral ACA supply cerebral area. Therefore, the f-values of the MCA and PCA supply areas of the contralateral hemisphere were not significantly correlated with postoperative CHS.

In this study, the rADC, rD, and rD* values of all areas in both ipsilateral and contralateral hemispheres showed no significant difference between the CHS and non-CHS groups (P > 0.05). As mentioned above, the D and ADC values were not associated with microvascular perfusion. Each operative patient with MMD in both the CHS and non-CHS groups showed no acute infraction; thus, the D and ADC values had no significant difference between the two groups. Previous studies have shown that the D* value has poor stability and may be affected by the CSF (Wirestam et al., 2001). Intracranial ischemia and encephalomalacia are common in patients with no acute infraction. Although we tried our best to avoid cerebral sulci, intracerebral ischemia, and encephalomalacia when placing the ROIs, the CSF may affect the D* value more or less. This may explain why the D* value was not significantly different between the CHS and non-CHS groups.

We used the rf value of the ACA supply area in the ipsilateral and contralateral hemispheres to generate ROC curves for the prediction of postoperative cerebral hyperperfusion. The results indicated that the AUC for the rf values was higher than that for the rfD* values to predict postoperative CHS. And the AUC for the rf values of the ACA, MCA, and PCA supply areas in the ipsilateral hemisphere was higher than that for the rf values of the ACA supply area in the contralateral hemisphere to predict postoperative CHS. These results also indicated that the microvascular perfusion status of the ipsilateral hemisphere was more related to the postoperative CHS than that of the contralateral hemisphere. The blood flow of the anastomosed vessel increases immediately after revascularization surgery, which may directly affect the microvascular perfusion status of the ipsilateral surgery hemisphere. The microvascular perfusion status of the ACA supply area in the contralateral hemisphere may also be affected by the anterior communicating artery. With respect to the prediction of postoperative CHS, in the ipsilateral hemisphere, the AUC for the rf values of the MCA supply areas was higher than that for the rf values of the ACA and PCA supply areas. This was also consistent with the operation. In this study, all patients underwent STA-MCA bypass surgery; hence, the preoperative microvascular perfusion status of the MCA supply area may show the highest relation to postoperative CHS.

This study has some limitations. First, the sample size was small; therefore, identifying the independent predictors of CHS using multivariate logistic regression is not feasible. Second, CHS was diagnosed mainly based on the patient’s clinical symptoms, instead of the standard diagnostic method (van Mook et al., 2005; Kawamata et al., 2011). Third, this study lacked microvascular perfusion information of patients with MMD who developed postoperative CHS because most patients who developed postoperative CHS were not appropriate to perform the perfusion study. A previous study showed that f is assumed to be proportional to CBV and D* is assumed to be proportional to the reciprocal value of the mean transit time (Federau et al., 2014). This study lacked other perfusion techniques for comparison. This will be improved in future studies.

IVIM can be used to evaluate the microvascular perfusion status of patients with MMD. Preoperative non-invasive IVIM-MRI analysis, particularly the f-value of the ipsilateral hemisphere, may be helpful in predicting CHS in patients with MMD after surgery. Patients with MMD with ischemic onset symptoms are more likely to develop CHS after surgery.

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 Huadong Hospital. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

WZ, FG, and GL contributed to conception and design of the study. WZ, FG, YZ, and ZZ organized the database. WZ and FG performed the statistical analysis and wrote the first draft of the manuscript. YZ, YD, and ZZ wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

This study was funded by the National Natural Science Foundation Project (Grant No. 81771816).

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/fnins.2022.826021/full#supplementary-material

References

Antonucci, M. U., Burns, T. C., Pulling, T. M., Rosenberg, J., Marks, M. P., Steinberg, G. K., et al. (2016). Acute preoperative infarcts and poor cerebrovascular reserve are independent risk factors for severe ischemic complications following direct Extracranial-Intracranial bypass for moyamoya disease. AJNR Am. J. Neuroradiol. 37, 228–235. doi: 10.3174/ajnr.A4535

PubMed Abstract | CrossRef Full Text | Google Scholar

Bisdas, S., Koh, T. S., Roder, C., Braun, C., Schittenhelm, J., Ernemann, U., et al. (2013). Intravoxel incoherent motion diffusion-weighted MR imaging of gliomas: Feasibility of the method and initial results. Neuroradiology 55, 1189–1196. doi: 10.1007/s00234-013-1229-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Callewaert, B., Jones, E., Himmelreich, U., and Gsell, W. (2021). Non-Invasive evaluation of cerebral microvasculature using Pre-Clinical MRI: principles, advantages and limitations. Diagnostics (Basel) 11:926. doi: 10.3390/diagnostics11060926

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, W. S., Kim, J. E., Kim, C. H., Ban, S. P., Kang, H. S., Son, Y. J., et al. (2014). Long-term outcomes after combined revascularization surgery in adult moyamoya disease. Stroke 45, 3025–3031. doi: 10.1161/strokeaha.114.005624

CrossRef Full Text | Google Scholar

Dusick, J. R., Gonzalez, N. R., and Martin, N. A. (2011). Clinical and angiographic outcomes from indirect revascularization surgery for Moyamoya disease in adults and children: a review of 63 procedures. Neurosurgery 68, 34–43, discussion 43. doi: 10.1227/NEU.0b013e3181fc5ec2

PubMed Abstract | CrossRef Full Text | Google Scholar

Federau, C., Maeder, P., O’Brien, K., Browaeys, P., Meuli, R., and Hagmann, P. (2012). Quantitative measurement of brain perfusion with intravoxel incoherent motion MR imaging. Radiology 265, 874–881. doi: 10.1148/radiol.12120584

PubMed Abstract | CrossRef Full Text | Google Scholar

Federau, C., O’Brien, K., Meuli, R., Hagmann, P., and Maeder, P. (2014). Measuring brain perfusion with intravoxel incoherent motion (IVIM): initial clinical experience. J. Magn. Reson. Imaging 39, 624–632. doi: 10.1002/jmri.24195

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujimura, M., Inoue, T., Shimizu, H., Saito, A., Mugikura, S., and Tominaga, T. (2012). Efficacy of prophylactic blood pressure lowering according to a standardized postoperative management protocol to prevent symptomatic cerebral hyperperfusion after direct revascularization surgery for moyamoya disease. Cerebrovasc. Dis. 33, 436–445. doi: 10.1159/000336765

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujimura, M., Niizuma, K., Endo, H., Sato, K., Inoue, T., Shimizu, H., et al. (2015). Quantitative analysis of early postoperative cerebral blood flow contributes to the prediction and diagnosis of cerebral hyperperfusion syndrome after revascularization surgery for moyamoya disease. Neurol. Res. 37, 131–138. doi: 10.1179/1743132814Y.0000000432

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujimura, M., Shimizu, H., Inoue, T., Mugikura, S., Saito, A., and Tominaga, T. (2011). Significance of focal cerebral hyperperfusion as a cause of transient neurologic deterioration after extracranial-intracranial bypass for moyamoya disease: comparative study with non-moyamoya patients using N-isopropyl-p-[(123)I]iodoamphetamine single-photon emission computed tomography. Neurosurgery 68, 957–964, discussion 964–965.

Google Scholar

Fukui, M. (1997). Guidelines for the diagnosis and treatment of spontaneous occlusion of the circle of Willis (‘moyamoya’ disease). Research committee on spontaneous occlusion of the circle of Willis (moyamoya disease) of the ministry of health and welfare, Japan. Clin. Neurol. Neurosurg. 99(Suppl. 2), S238–S240.

PubMed Abstract | Google Scholar

Guiu, B., Petit, J. M., Capitan, V., Aho, S., Masson, D., Lefevre, P. H., et al. (2012). Intravoxel incoherent motion diffusion-weighted imaging in nonalcoholic fatty liver disease: a 3.0-T MR study. Radiology 265, 96–103. doi: 10.1148/radiol.12112478

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, J. S., Abou-Hamden, A., Mandell, D. M., Poublanc, J., Crawley, A. P., Fisher, J. A., et al. (2011). Impact of extracranial-intracranial bypass on cerebrovascular reactivity and clinical outcome in patients with symptomatic moyamoya vasculopathy. Stroke 42, 3047–3054. doi: 10.1161/STROKEAHA.111.615955

PubMed Abstract | CrossRef Full Text | Google Scholar

Hara, S., Hori, M., Ueda, R., Hagiwara, A., Hayashi, S., Inaji, M., et al. (2019). Intravoxel incoherent motion perfusion in patients with Moyamoya disease: comparison with (15)O-gas positron emission tomography. Acta Radiol. Open 8:1741009221. doi: 10.1177/2058460119846587

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayashi, K., Horie, N., Suyama, K., and Nagata, I. (2012). Incidence and clinical features of symptomatic cerebral hyperperfusion syndrome after vascular reconstruction. World Neurosurg. 78, 447–454. doi: 10.1016/j.wneu.2011.10.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayashi, T., Shirane, R., Fujimura, M., and Tominaga, T. (2010). Postoperative neurological deterioration in pediatric moyamoya disease: watershed shift and hyperperfusion. J. Neurosurg. Pediatr. 6, 73–81. doi: 10.3171/2010.4.PEDS09478

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwata, T., Mori, T., Miyazaki, Y., Tanno, Y., Kasakura, S., and Aoyagi, Y. (2014). Global oxygen extraction fraction by blood sampling to anticipate cerebral hyperperfusion phenomenon after carotid artery stenting. Neurosurgery 75, 546–551, discussion 551. doi: 10.1227/NEU.0000000000000485

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaku, Y., Iihara, K., Nakajima, N., Kataoka, H., Fukuda, K., Masuoka, J., et al. (2012). Cerebral blood flow and metabolism of hyperperfusion after cerebral revascularization in patients with moyamoya disease. J. Cereb. Blood Flow Metab. 32, 2066–2075. doi: 10.1038/jcbfm.2012.110

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawamata, T., Kawashima, A., Yamaguchi, K., Hori, T., and Okada, Y. (2011). Usefulness of intraoperative laser Doppler flowmetry and thermography to predict a risk of postoperative hyperperfusion after superficial temporal artery-middle cerebral artery bypass for moyamoya disease. Neurosurg. Rev. 34, 355–362, discussion 362. doi: 10.1007/s10143-011-0331-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazumata, K., Ito, M., Tokairin, K., Ito, Y., Houkin, K., Nakayama, N., et al. (2014). The frequency of postoperative stroke in moyamoya disease following combined revascularization: a single-university series and systematic review. J. Neurosurg. 121, 432–440. doi: 10.3171/2014.1.JNS13946

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, T., Oh, C. W., Bang, J. S., Kim, J. E., and Cho, W. S. (2016). Moyamoya disease: treatment and outcomes. J. Stroke 18, 21–30.

Google Scholar

Le Bihan, D. (2019). What can we see with IVIM MRI? Neuroimage 187, 56–67. doi: 10.1016/j.neuroimage.2017.12.062

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Bihan, D., Breton, E., Lallemand, D., Grenier, P., Cabanis, E., and Laval-Jeantet, M. (1986). MR imaging of intravoxel incoherent motions: Application to diffusion and perfusion in neurologic disorders. Radiology 161, 401–407. doi: 10.1148/radiology.161.2.3763909

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H. J., Rha, S. Y., Chung, Y. E., Shim, H. S., Kim, Y. J., Hur, J., et al. (2014). Tumor perfusion-related parameter of diffusion-weighted magnetic resonance imaging: correlation with histological microvessel density. Magn. Reson. Med. 71, 1554–1558. doi: 10.1002/mrm.24810

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, S. U., Oh, C. W., Kwon, O. K., Bang, J. S., Ban, S. P., Byoun, H. S., et al. (2018). Surgical treatment of adult moyamoya disease. Curr. Treat. Options Neurol. 20:22.

Google Scholar

Lemke, A., Laun, F. B., Klauss, M., Re, T. J., Simon, D., Delorme, S., et al. (2009). Differentiation of pancreas carcinoma from healthy pancreatic tissue using multiple b-values: comparison of apparent diffusion coefficient and intravoxel incoherent motion derived parameters. Investig. Radiol. 44, 769–775. doi: 10.1097/RLI.0b013e3181b62271

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsushima, Y., Aoyagi, M., Niimi, Y., Masaoka, H., and Ohno, K. (1990). Symptoms and their pattern of progression in childhood moyamoya disease. Brain Dev. 12, 784–789. doi: 10.1016/s0387-7604(12)80007-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Mesiwala, A. H., Sviri, G., Fatemi, N., Britz, G. W., and Newell, D. W. (2008). Long-term outcome of superficial temporal artery-middle cerebral artery bypass for patients with moyamoya disease in the US. Neurosurg. Focus 24:E15. doi: 10.3171/FOC/2008/24/2/E15

PubMed Abstract | CrossRef Full Text | Google Scholar

Okada, Y., Kawamata, T., Kawashima, A., Yamaguchi, K., Ono, Y., and Hori, T. (2012). The efficacy of superficial temporal artery-middle cerebral artery anastomosis in patients with moyamoya disease complaining of severe headache. J. Neurosurg. 116, 672–679. doi: 10.3171/2011.11.JNS11944

PubMed Abstract | CrossRef Full Text | Google Scholar

Ostergaard, L., Engedal, T. S., Aamand, R., Mikkelsen, R., Iversen, N. K., Anzabi, M., et al. (2014). Capillary transit time heterogeneity and flow-metabolism coupling after traumatic brain injury. J. Cereb. Blood Flow Metab. 34, 1585–1598. doi: 10.1038/jcbfm.2014.131

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, W., Park, E. S., Lee, S., Park, J. C., Chung, J., Lee, J. M., et al. (2018). Intracranial hemorrhage after superficial temporal Artery-middle cerebral artery direct anastomosis for adults with moyamoya disease. World Neurosurg. 119, e774–e782. doi: 10.1016/j.wneu.2018.07.266

PubMed Abstract | CrossRef Full Text | Google Scholar

Paschoal, A. M., Leoni, R. F., Dos, S. A., and Paiva, F. F. (2018). Intravoxel incoherent motion MRI in neurological and cerebrovascular diseases. Neuroimage Clin. 20, 705–714. doi: 10.1016/j.nicl.2018.08.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Research Committee on the Pathology and Treatment of Spontaneous Occlusion of the Circle of Willis, and Health Labour Sciences Research Grant for Research on Measures for Infractable Diseases (2012). Guidelines for diagnosis and treatment of moyamoya disease (spontaneous occlusion of the circle of Willis). Neurol. Med. Chir. (Tokyo) 52, 245–266. doi: 10.2176/nmc.52.245

PubMed Abstract | CrossRef Full Text | Google Scholar

Riches, S. F., Hawtin, K., Charles-Edwards, E. M., and de Souza, N. M. (2009). Diffusion-weighted imaging of the prostate and rectal wall: comparison of biexponential and monoexponential modelled diffusion and associated perfusion coefficients. NMR Biomed. 22, 318–325. doi: 10.1002/nbm.1328

PubMed Abstract | CrossRef Full Text | Google Scholar

Sigmund, E. E., Cho, G. Y., Kim, S., Finn, M., Moccaldi, M., Jensen, J. H., et al. (2011). Intravoxel incoherent motion imaging of tumor microenvironment in locally advanced breast cancer. Magn. Reson. Med. 65, 1437–1447. doi: 10.1002/mrm.22740

PubMed Abstract | CrossRef Full Text | Google Scholar

Sumi, M., and Nakamura, T. (2013). Head and neck tumors: assessment of perfusion-related parameters and diffusion coefficients based on the intravoxel incoherent motion model. AJNR Am. J. Neuroradiol. 34, 410–416. doi: 10.3174/ajnr.A3227

PubMed Abstract | CrossRef Full Text | Google Scholar

Suo, S., Cao, M., Zhu, W., Li, L., Li, J., Shen, F., et al. (2016). Stroke assessment with intravoxel incoherent motion diffusion-weighted MRI. NMR Biomed. 29, 320–328. doi: 10.1002/nbm.3467

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, J., and Takaku, A. (1969). Cerebrovascular “moyamoya” disease. Disease showing abnormal net-like vessels in base of brain. Arch. Neurol. 20, 288–299. doi: 10.1001/archneur.1969.00480090076012

PubMed Abstract | CrossRef Full Text | Google Scholar

Thoeny, H. C., Binser, T., Roth, B., Kessler, T. M., and Vermathen, P. (2009). Noninvasive assessment of acute ureteral obstruction with diffusion-weighted MR imaging: a prospective study. Radiology 252, 721–728. doi: 10.1148/radiol.2523082090

PubMed Abstract | CrossRef Full Text | Google Scholar

Togao, O., Mihara, F., Yoshiura, T., Tanaka, A., Noguchi, T., Kuwabara, Y., et al. (2006). Cerebral hemodynamics in Moyamoya disease: correlation between perfusion-weighted MR imaging and cerebral angiography. AJNR Am. J. Neuroradiol. 27, 391–397.

PubMed Abstract | Google Scholar

van Mook, W. N., Rennenberg, R. J., Schurink, G. W., van Oostenbrugge, R. J., Mess, W. H., Hofman, P. A., et al. (2005). Cerebral hyperperfusion syndrome. Lancet Neurol. 4, 877–888.

Google Scholar

Wardlaw, J. M., Smith, E. E., Biessels, G. J., Cordonnier, C., Fazekas, F., Frayne, R., et al. (2013). Neuroimaging standards for research into small vessel disease and its contribution to ageing and neurodegeneration. Lancet Neurol. 12, 822–838. doi: 10.1016/S1474-4422(13)70124-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wirestam, R., Borg, M., Brockstedt, S., Lindgren, A., Holtas, S., and Stahlberg, F. (2001). Perfusion-related parameters in intravoxel incoherent motion MR imaging compared with CBV and CBF measured by dynamic susceptibility-contrast MR technique. Acta Radiol. 42, 123–128.

PubMed Abstract | Google Scholar

Wu, W. C., Chen, Y. F., Tseng, H. M., Yang, S. C., and My, P. C. (2015). Caveat of measuring perfusion indexes using intravoxel incoherent motion magnetic resonance imaging in the human brain. Eur. Radiol. 25, 2485–2492. doi: 10.1007/s00330-015-3655-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, Y., Zhang, S., Tang, X., Zhang, S., Shi, J., Zhu, W., et al. (2016). Intravoxel incoherent motion diffusion-weighted imaging in stroke patients: initial clinical experience. Clin. Radiol. 71, 911–938. doi: 10.1016/j.crad.2016.04.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaharchuk, G., Do, H. M., Marks, M. P., Rosenberg, J., Moseley, M. E., and Steinberg, G. K. (2011). Arterial spin-labeling MRI can identify the presence and intensity of collateral perfusion in patients with moyamoya disease. Stroke 42, 2485–2491. doi: 10.1161/STROKEAHA.111.616466

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, W. G., Luo, Q., Jia, J. B., and Yu, J. L. (2013). Cerebral hyperperfusion syndrome after revascularization surgery in patients with moyamoya disease. Br. J. Neurosurg. 27, 321–325.

Google Scholar

Keywords: moyamoya disease, perfusion imaging, magnetic resonance imaging, microvascular circulation, surgery

Citation: Gao F, Zhao W, Zheng Y, Duan Y, Ji M, Lin G and Zhu Z (2022) Intravoxel Incoherent Motion Magnetic Resonance Imaging Used in Preoperative Screening of High-Risk Patients With Moyamoya Disease Who May Develop Postoperative Cerebral Hyperperfusion Syndrome. Front. Neurosci. 16:826021. doi: 10.3389/fnins.2022.826021

Received: 30 November 2021; Accepted: 26 January 2022;
Published: 02 March 2022.

Edited by:

Kevin C. Chan, NYU Grossman School of Medicine, United States

Reviewed by:

Steve Huntz Fung, Houston Methodist Research Institute, United States
Matthew ManHin Cheung, CUHK Medical Centre, Hong Kong SAR, China

Copyright © 2022 Gao, Zhao, Zheng, Duan, Ji, Lin and Zhu. 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: Feng Gao, Z2FvZmVuZ2ExQDEyNi5jb20=; Zhenfang Zhu, aGRrYWl4aW4wMDFAc2luYS5jb20=

These authors have contributed equally to this work and share first authorship

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.