Abstract
Background Acute ischemic stroke caused by large-vessel occlusion is effectively treated by endovascular treatment (EVT). However, treatment could be further refined by improved understanding of the underlying pathophysiology, including dynamic cerebral autoregulation (dCA). Near-infrared spectroscopy (NIRS) requires virtually no setup time and enables dCA investigation during EVT by measuring dynamic concentration in cortical oxygenated hemoglobin (OxyHb) continuously.
Methods The aim was to investigate dCA during EVT, before and after recanalization, and at follow-up after 24 hours and 90 days. We applied interhemispheric transfer function analysis of low-frequency (LF) oscillations (0.07-0.2 Hz) in OxyHb that yields the dCA measures, gain and phase shift. Mixed-effect models were used to examine dCA over time. We then predicted outcome by dCA in multivariate logistic regression models.
Results We found higher LF gain in patients with favorable outcome. While dCA was stable for most patients during EVT, LF phase shift increased immediately after recanalization for patients with milder symptom severity and at 24-hour follow-up for patients with more severe symptoms. Adjusted for age, infarct size before EVT, and recanalization success, average LF gain predicted independent functional outcome, symptom severity and mortality at 90-day follow-up.
Conclusions Interhemispheric transfer function analysis based on NIRS provides insights into dCA changes during EVT that could potentially be applied in the development of individualized treatment. The method relies on intact contralateral dCA which needs further investigations. The potential of dCA to strengthen prediction models also warrant additional examination.
Clinical Trial Registration ClinicalTrial.gov: NCT03738644. https://clinicaltrials.gov/study/NCT03738644.
Introduction
Stroke remains a worldwide leading cause of disability and death despite the emergence of effective recanalization therapies.1 Patients with a large-vessel occlusion (LVO) face the worst prognosis and constitute up to 38% of patients with acute ischemic stroke (AIS).2 Endovascular treatment (EVT) drastically improves functional outcome after anterior circulation LVO when performed within 6 hours3, or within 24 hours in cases with mismatch between neurologic deficit and infarct core.4 Posterior circulation EVT has shown similar advantages compared to best medical treatment.5 However, procedure-related complications and futile recanalization still restrain the overall efficacy of EVT6,7. Impairment of cerebral autoregulation is one of the proposed mechanisms for unexplained futile recanalization.7,8
Cerebral autoregulation is the maintenance of suitable blood flow despite changes in cerebral perfusion pressure.9 Dynamic cerebral autoregulation (dCA) determines how cerebral blood flow is regulated during rapid changes in perfusion pressure.10 Prevailing methods indicate varying degrees of dCA impairment in patients with AIS including LVO8,11 which could endanger vulnerable brain tissue by hypo-and hyperperfusion during and after EVT. Post-thrombectomy studies have shown dCA is an independent predictor of functional outcome12, and that it can be used to individualize blood pressure management.13 However, dCA has never been examined before or during EVT.
Most dCA investigations rely on transcranial Doppler sonography examining blood flow velocity in the middle cerebral artery (VMCA) which requires expertise and setup time.14 Assessment of dCA before EVT would compromise patient safety when every minute counts. Monitoring dCA during EVT also require modalities that do not interfere with the digital subtraction angiography. Near-infrared spectroscopy (NIRS) is an optical method which continuously measures dynamic changes in cortical hemoglobin concentrations without needing significant setup time or causing artifacts on the digital subtraction angiography.15,16
Low-frequency oscillations (LFO, approximately 0.1 Hz) are prevalent in both systemic (i.e., arterial blood pressure, ABP) and cerebral circulation (e.g., VMCA, NIRS).15,17,18 Quantifying the changes between them in the frequency domain by transfer function analysis (TFA) is one of the most utilized and standardized methods of assessing dCA.19,20 Applying TFA for comparison of the NIRS signal between the unaffected and the ischemic hemisphere, we aimed to investigate dCA before, during and after EVT.
Methods
This prospective observational study (ClinicalTrials.gov: NCT03738644) was conducted between November 2018 and November 2020 after approval from the Scientific Ethics Committees for the Capital Region of Denmark (H-18028704) and in accordance with the World Medical Association Declaration of Helsinki. All participants or their proxy signed written informed consent. Reporting complies with STrengthening the Reporting of OBservational studies in Epidemiology (STROBE) guidelines for observational studies.21 Corresponding author had full access to all the data in the study and take responsibility for data integrity and analysis. Further methodical details are provided in supplementary materials.
All AIS patients admitted to Rigshospitalet with large-vessel occlusion receiving EVT were screened for eligibility. The exclusion criteria are shown in the enrollment flowchart (Figure 1). Patients with bilateral infarction were excluded as no contralateral hemisphere could be defined (see interhemispheric transfer function analysis section).
Patients were diagnosed and treated as per standard-of-care. Large-vessel occlusion was confirmed and Alberta Stroke Program Early CT Score (ASPECTS)22 or Posterior circulation-ASPECTS (PC-ASPECTS)23 were determined by staff neuroradiologists assisted by RAPID AI software (iSchemaView, Menlo Park, California, USA) and dichotomized to favorable (ASPECTS/PC-ASPECTS ≥ 6) or unfavorable (ASPECTS/PC-ASPECTS < 6). Patients were transferred immediately to the angiosuite where symptom severity was re-assessed using the National Institutes of Health Stroke Scale (NIHSS) before EVT was performed under general anesthesia. In addition to conventional monitoring,24 NIRS was applied throughout the procedure (described below). Reperfusion success was assessed by the treating interventional neuroradiologist (modified Treatment In Cerebral Infarction, mTICI; Grade 0-2a unsuccessful achieving no reperfusion to reperfusion in less than half of the occluded artery territory; Grade 2b-3 successful reperfusion in more than half of the occluded artery territory or complete reperfusion).25 Post-procedure, the patients were extubated as soon as possible and underwent dual-energy computed tomography (CT) or magnetic resonance imaging (MRI) after 24 hours to detect intracranial hemorrhage or other complications. Patients underwent further examinations as per standard-of-care. Stroke etiology was classified based on the Causative Classification System for Ischemic Stroke.26
Follow-up
Patients had two follow-up (FU) examinations with NIRS in the supine position for 20 minutes after at least 5 minutes of rest. The first FU was performed 24 hours after EVT (+/-6 hours from end of endovascular reperfusion efforts) in addition to NIHSS assessment. The second FU was completed after 90 days (+/- 14 days) with functional outcome scored modified Rankin Scale (mRS) and independency (defined as mRS of 0-2) as well as NIHSS. Occurrence of re-hospitalizations, new vascular events (recurrent stroke, myocardial infarction, or surgery for peripheral artery disease) and all-cause mortality were registered. Patients unable to attend 90-day FU, were offered home visits within a 2-hour transportation radius. Patients who did not have in-person FU, were followed up by phone interview and using the electronic health record system.
NIRS examination
The NIRS monitoring was performed using a continuous wave-NIRS (CW-NIRS) system (Octamon, Artinis Medical Systems, Elst, the Netherlands) with three long-distance channels (35 mm) per hemisphere, examining dynamic hemoglobin concentrations in the prefrontal cortex in the border zones between territories of the middle cerebral (MCA) and the anterior cerebral arteries (ACA) (two channels) and ACA territories exclusively (one channel per side).27 A fourth channel with short diode-receiver distance (10 mm) mainly examined the extracerebral tissues on both sides.
Time segment selection
Per recommendation from the Cerebrovascular Research Network only steady state data segments of five minutes were analyzed.20 Thus, data segments were excluded within 5 minutes of inducing general anesthesia and within 2 minutes of changes in anesthetics, opioids, or vasopressors.28 Very noisy data segments were also excluded from the analysis. Individual NIRS signal and all vitals were evaluated for steady state with an allowed variation of 10%. The first data segments free of these exclusionary criteria were chosen.
The first segment (PRE) occurred after sedation but before any attempted revascularization. The second data segment (POST) was just after the final reperfusion status was achieved, including the abandonment of reperfusion efforts. Third (24-hour) and fourth (90-day) data segments were selected based on the extent of motion artifacts. Patients who did not have complete PRE, POST and 24-hour segments were excluded from further analysis.
Interhemispheric transfer function analysis
The Cerebrovascular Research Network guidelines for TFA20 were adapted with minor exceptions described throughout. TFA was only performed on oxygenated hemoglobin (OxyHb) concentrations as they are more reliable29 and better reflect the arterial compartment.15 Input to TFA was OxyHb from the contralateral hemisphere, while the corresponding contralateral channel from the ischemic hemisphere served as output. Fourier transformation was performed, and gain (amplitude ratio) and phase shift were then calculated in three frequency intervals: High-frequency (HF, 0.2-0.5 Hz), low-frequency (LF, 0.07-0.2 Hz) and very low-frequency (VLF, 0.02-0.07 Hz). Power spectral density (PSD) specified amplitudes ((µM*mm)2/Hz) for each hemisphere.
Positive phase shifts denoted LFOs in the contralateral hemisphere occurring before LFOs in the ischemic hemisphere and vice versa. Because dyssynchronization of phase shifts can occur in either direction, absolute phase shift was calculated. Similarly, gain of more than 1 resulted from amplitudes in the ischemic hemisphere exceeding amplitudes in the contralateral hemisphere and vice versa. Theoretically, dCA would be deemed intact with no substantial phase shift and gain equal to one. This interpretation assumes dCA in the contralateral hemisphere is intact.
Any differences between input and output signals were assumed to result from intracerebral changes as extracerebral circulation is not under the influence of cerebral autoregulation. We tested this assumption by performing a TFA between short-distance channels from each side with the contralateral side as input and the affected side as output at 24-hour FU, under which the clinical condition was most stable and less prone to artifacts. Segment selection, data preparation and TFA methodology were equal to the processing for long-distance channels described above. Unfortunately, short-distance channels are more prone to oversaturation, reducing the number of channels with adequate data quality, especially in acute settings. Thus, incorporation of short-distance channel TFA in statistical models was not possible.
Statistics
Statistical analysis was performed in two subsets of participants: 1) All participants with sufficient data at PRE, POST and 24-hour FU (FU24-group) 2) All patients with sufficient data at PRE, POST, 24-hour and 90-day FU (FU90-group) excluding patients with no in-person FU and patients who died. Baseline data from subjects without in-person FU was also analyzed to identify differences in subsets.
Normally distributed data are presented as mean and standard deviation (SD), while non-normal data are presented as median and interquartile range (IQR). Two-sample testing was performed as appropriate. False discovery rate was applied when performing multiple comparisons.
Linear mixed-effect models were fitted with LF gain or LF absolute phase shift as outcome, subjects as random effects, and fixed effects comprised of time segment in addition to and different patient characteristics (e.g., stroke etiology, age, favorable ASPECTS before EVT, anterior or posterior circulation stroke), treatment (e.g., anesthesia, vasopressor, successful recanalization) or outcome (e.g., NIHSS, mRS, mortality). In case of significant fixed effects, mixed-effect models were fitted with LF PSD from each hemisphere as an outcome to examine the LFO amplitude in each hemisphere. Interaction between time segment and other fixed effects were assessed to determine different temporal developments between groups.
Significant fixed effects from mixed-effect models were used as outcome in logistic regression models with corresponding TFA parameter as univariate predictor. Two multivariate predictions were then performed. The acute prediction model was adjusted for age, recanalization success, and favorable ASPECTS/PC-ASPECTS.30 NIHSS was not used as co-variate as some patients were intubated before arrival in the angiosuite and did not have a NIHSS before EVT. The 24-hour prediction was modelled with 24-hour NIHSS and age as co-variates. For dichotomous outcomes, receiver operating characteristics curves were generated, and area under the receiver operating characteristic curve (AUROC) was generated to assess the accuracy of the predictor in both univariate and multivariate models without validation subsets.
Models estimates (β) are presented with 95% confidence intervals (CI). Significance level was chosen at 5%.
Results
Baseline characteristics
Baseline information and medical history of all enrolled subjects and grouped into subjects with and without in-person 90-day FU are listed in Table 1 and Table S1. No differences between the FU and non-FU groups were detected.
Index stroke characteristics, treatment and outcome are listed in Table 2 and Table S2. Most of the living patients had in-person FU (82.8%). There were no significant differences concerning etiology, onset, occluded artery, procedural technique, complications, anesthetics, or vasopressors. Subjects who did not complete the 90 day-FU had higher NIHSS before EVT, lower proportion of successful reperfusion and longer time from last-known-well to reperfusion possibly leading to lower 24-hour ASPECTS/PC-ASPECTS as well as higher 24-hour NIHSS, 90-day mRS, and all cause-mortality. Non-FU subjects also experienced a numerically higher occurrence of complications and vascular events. Removing subjects with fatal outcome within 90 days from the non-FU group, did not alter the differences between groups (data not shown).
Short distance channels TFA
TFA of short distance channel at 24-hour FU is shown in Table S3. Across frequency ranges, gain did not differ from 1 and phase shift did not differ from 0. We found no side-to-side difference for mean OxyHb and PSD.
Average TFA results
The median time from last-known-well to PRE was 4.0 hours (2.5; 7.9). Mean OxyHb decreased at both POST and 24-hour FU similarly in both hemispheres. PSD increased bilaterally in both VLF and LF range after 24 hours. Overall, vitals were quite steady. TFA showed increased absolute phase difference after 24 hours in the LF range and a temporary reduction after recanalization in the VLF range normalizing after 24 hours. Mean OxyHb, vital parameters, PSD, and the results of the TFA for the FU24-group are listed in Table S3.
In all ranges, FU subjects showed a PSD reduction between 24 hours and 90 days. Coherence and gain were stable between 24 hours and 90 days. LF and VFL absolute phase differences showed significant reduction, while quite stable in other ranges. All results for the FU90-group are shown in Table S4.
Mixed-effect models
LF gain
In the FU24-group LF gain was significantly lower in subjects with increasing mRS (β=-0.04, CI=[-0.07; -0.02], p=0.002), dependent outcome (β=-0.13, CI=[-0.24; -0.02], p=0.019), increasing NIHSS (β=- 1.0%, CI=[-1.8%; -0.1%], p=0.029), and in subjects with fatal 90-day outcome (β=-0.21, CI=[-0.35;- 0.06], p=0.005) across PRE, POST and 24-hour data segments. The observed reduction in LF gain was the result of greater PSD reduction in the ischemic hemisphere.
Grouping patients by anesthetics, vasopressor, stroke etiology (large-artery atherosclerosis or cardioembolic), age, occlusion complication level (simple isolated thrombus or complicated ICA-top or tandem thrombus), favorable ASPECTS/PC-ASPECTS before EVT, anterior or posterior circulation stroke, and hypotension during recording (average mean arterial pressure (MAP) < 70 mmHg) had no effect on LF gain. Grouping patients by recanalization success (unsuccessful: β=-0.15, CI=[-0.33; 0.03], p=0.105) and hypertension (MAP< 130 mmHg before recanalization or 90 mmHg after recanalization) during recording (β=0.09, CI=[-0.01; 0.19], p=0.082) showed no statistical change in LF gain. There was no interaction effect between the time segment and any other fixed effects. Time segment had no effect on LF gain with PRE as index.
Associations between LF gain and outcome were not prevalent in the FU90 group. While LF gain was still numerically lower in subjects with worse outcomes, the difference did not reach a significance level. In contrast to FU24-group a significant interaction effect between favorable ASPECTS/PC-ASPECTS and the time segment was observed (p=0.035) as LF gain increased more after recanalization and declined more at 24-hour and 90-day FU in patients with non-favorable ASPECTS/PC-ASPECTS. The interaction effect is plotted in Figure S1. LF gain was numerically lower at 24-hour and 90-day FU with PRE as index and significantly lower compared to POST.
LF phase shift
Mixed-effect models showed no effect on absolute phase shift of grouping patients by mRS, mortality, vasopressor, stroke etiology, age, occlusion complication level, favorable ASPECTS/PC-ASPECTS before EVT, anterior or posterior circulation stroke, hypotension nor hypertension during recording. Compared to patients anesthetized with propofol, those receiving sevoflurane showed a trend towards numerically lower absolute phase shift (β=-46.0%, CI=[-71.9%; +3.9%], p=0.064). We found an interaction between NIHSS and time segment as subjects with milder stroke severity had an immediate increase in absolute phase shift after recanalization, as opposed to patients with greater stroke severity increasing in absolute phase shift between POST and 24-hour FU (Figure 2).
Results from 24-hour and 90-day subsets was equivalent. Absolute phase shift decreased between 24-hour and 90-day FU and was not different from index (PRE) at 90-day FU.
All results from linear mixed-effect models can be found in Table S6-S9.
Prediction models
LF gain showed associations to outcome and therefore, was used as co-variate in the following prediction models. LF gain was averaged across PRE, POST, and 24-hour FU because mixed-effects models showed no effect of time segment on LF gain. Thus, average LF gain was representable for LF gain in the first 24 hours. It could theoretically still be used to predict outcome in the large proportion of enrolled subjects excluded from analysis due to missing state period (16%). The full output from multivariate models is reported in supplementary materials.
Modified Rankin scale
Logistic regression models were used to predict 90-day independent (mRS 0-2) or dependent (mRS > 2) functional outcome with LF gain as predicting variable. Average LF gain was significant predictor of independent outcome (β=2.36, CI=[0.41; 4.53], p=0.023, AUROC=0.641) with an odds ratio (OR) of 1.60 (CI=[1.07;2.41]) for a 0.2 increase in LF gain.
In the acute predictions model (adjusted for age, recanalization, and favorable ASPECTS/PC-ASPECTS before EVT), average LF gain remained a significant predictor of independent outcome (β=2.60, CI=[0.32; 5.22], p=0.035, AUROC=0.802, Table S10) with an OR of 1.68 (CI=[1.04; 2.73]) for a 0.2 increase in LF gain. Age was the only other significant predictor.
The 24-hour model adjusted for age and trichotomized 24-hour NIHSS (0-5; 6-13; >14) and showed that LF gain (β=2.29, CI=[-0.14; 5.11], p=0.081, AUROC=0.876, Table S11) had an OR of 1.58 (CI=[0.94; 2.65] for 90-day independency with a 0.2 increment. Age and trichotomized 24-hour NIHSS were significant predictors.
NIHSS
Predicting 1-step changes in 90-day NIHSS is of minor clinical significance prompting NIHSS categorization (Near-remission 0-1; Mild 2-3; Moderate 4-6; Severe 7-18; Fatal outcome). Ordinal logistic regression model was used to predict 90-day categorized NIHSS with average LF gain as a significant predictor (β=2.14, CI=[0.42; 3.94], p=0.016). With an increase of 0.2 in LF gain, the common OR was 1.53 (CI=[1.08;2.17]) for better categorical 90-day NIHSS. Predicted probabilities for the univariate model is depicted in Figure 3.
In multivariate acute prediction, LF gain remained a significant predictor (β=2.26, CI=[0.45; 4.19], p=0.017, Table S12) with a common OR of 1.57 (CI=[1.08; 2.28]) for better categorical 90-day NIHSS. Favorable ASPECTS/PC-ASPECTS before EVT was the only other significant predictor.
The 24-hour model showed average LF gain (β=3.11, CI=[-0.46; 3.66], p=0.072, AUROC=0.876, Table S13) had an OR of 1.37 (CI=[0.91; 2.07] for better categorical 90-day NIHSS with a 0.2 increment.
Mortality
Logistic regression models were used to predict 90-day survival. In the univariate model, average LF gain was a significant predictor of 90-day survival (β=3.73, CI=[1.04; 6.85], p=0.011, AUROC=0.716) with an OR of 2.11 (CI=[1.19; 3.73]) for a0.2 increase in LF gain. In the acute model, average LF gain remained a significant predictor of 90-day survival (β=3.78, CI=[0.77; 7.35], p=0.022, AUROC=0.863, Table S14) with an OR of 2.13 (CI=[1.12; 4.06]) for a 0.2 increase in LF gain. The 24-hour model (Table S15) was strongly affected by collinearity between trichotomized NIHSS and LF gain.
Discussion
Dynamic CA (dCA) in AIS is of immense interest due to its pathophysiological role and consistent association with long-term outcome.8,31 Further, dCA has shown promising potential for individualized blood pressure management after thrombectomy.13 In order to investigate the potential for individualizing intraprocedural treatment it is important to study dCA in the hyperacute phase of stroke and large-vessel occlusion before, during and after revascularization. By our account, this study is the first to examine dCA not only during EVT but also before and after recanalization and possibly the earliest dCA assessment after ischemic stroke. Here, we have shown the time course of interhemispheric TFA based on cortical OxyHb across hyperacute, acute, and chronic phases of large-vessel occlusion and that LF gain was associated with and thus had predictive capabilities of 90-day symptom severity, functional outcome, and mortality.
Phase shift
Studies of dCA in AIS patients have predominantly been based on transcranial Doppler sonography and TFA between ABP and VMCA. Associations between impaired dCA in the ischemic hemisphere and poor outcome have been shown in numerous studies regarding infarct size, edema, hemorrhagic transformation, symptom severity at discharge and long-term functional outcome (mRS) but mostly with LF phase shift as covariate.8,32
While we did not find associations between interhemispheric phase shift and functional outcome, our data did show an interaction between NIHSS and time segments. Absolute phase shift increased immediately after recanalization in patients with milder symptom severity whereas patients with higher NIHSS exhibited an increase at 24-hour FU. The differences in temporal development of phase shift dyssynchronization could be related to the penumbra as there were no effect of ASPECTS/PC-ASPECTS or the recanalization success. Dichotomizing recanalization to complete (mTICI 3) or incomplete (mTICI 2a or 2b), Sheriff et al. also showed an interaction effect when examining ABP-VMCA phase shift over time. Such interactions emphasize the importance of examination time and the need to account for clinical characteristics when studying dCA in AIS patients. Phase shift between ABP and VMCA have shown no interhemispheric difference within 6 hours of onset with moderate to severe stroke symptoms in the MCA territory.33 Later, reductions compared to healthy controls and contralateral hemisphere are consistently seen from 20 hours and throughout 7 days from onset in patients with MCA occlusion,12,34,35 while perhaps already normalizing after 10 days.34
Using interhemispheric TFA, we found an absolute phase shift relatively close to 0, indicating equal dCA across hemispheres. Studies of contralateral phase shift dyssynchronization have shown inconsistent results which could result from vasoactive substances accumulating especially in LVO and affecting dCA even in the contralateral hemisphere.7,8 Contralateral dyssynchronization could explain why we found no association to functional outcome that is otherwise reported quite consistently in TFA of ABP-VMCA in the ischemic hemisphere. Examining bilateral ABP-VMCA phase shift and interhemispheric phase shift concurrently could reveal this in the future. In this study, simultaneous interhemispheric TFA and bilateral ABP-OxyHb TFA would have diminished the sample size greatly as many patients did not have sufficient ABP segments before recanalization.
While TFA of ABP and VMCA have been used in many studies of AIS patients, only one study has previously used OxyHb as both input and output.36 This study was performed 1-2 days after onset and showed increased LF absolute phase shift in non-thrombolyzed patients compared to thrombolyzed patients. The thrombolyzed patients had a higher NIHSS at admission, but a higher rate of total remission and conceivably less infarction volume, perhaps resulting in equal dCA across both hemispheres. The current study found no such difference as effect of recanalization success probably because of significant differences concerning study design, cohort and TFA methodology.
Gain
Similar to the current study, Castro et al. showed that LF gain between ABP and VMCA in the affected hemisphere of moderate to severe MCA strokes within 6 hours of onset was a predictor of 90-day dependency.37 Gain remained stable between hyperacute (<6 hours), acute (24 hours) and chronic (90-day) stroke stages,37 while inconsistent findings are reported in the subacute stage (3 to 7 days after onset).38,39 However, gain in the ischemic hemisphere increased with outcome severity, interpreted as an increasingly impaired dCA with a diminished ability to dampen the amplitude of LFOs from systemic perfusion.40 Conversely, a recent meta-analysis found lower gain and thus more intact dCA across all AIS patients compared to healthy controls.11 The inconsistencies concerning ABP-VMCA gain could be explained by heterogeneities between stroke populations and methodical variability including scaling methods that is required when using different modalities in TFA. LF gain in patients with favorable outcome from the current cohort was close to one, indicating equal dCA between hemispheres and intact dCA in the ischemic hemisphere under the assumption of intact contralateral dCA. With increasing outcome severity, LFO amplitude was relatively stable in the contralateral hemisphere while LFO amplitude in the ischemic hemisphere declined resulting in lower gain likely due to affected dCA. Thus, direct comparisons of LF gain between conventional TFA setup (ABP-VMCA) and interhemispheric TFA seems implausible perhaps due to differences in the examined vasculature. NIRS should be sensitive to all microvascular regulations, while that may not be true for VMCA, which is based on upstream arteriolar resistance.41 Multiple TFAs based on simultaneous measurements of ABP, VMCA and OxyHb would be needed in the future to determine how different TFA models interact.
We can only compare amplitude and gain directly to subacute patients examined by Phillip et al.42 Numerically amplitude in the ischemic hemisphere was higher resulting in higher gain in thrombolyzed patients while neither reached statistical significance.36 This could be consistent with the current study as thrombolyzed patients had higher remission rate than non-thrombolyzed patients, but NIHSS at time of examination and long-term outcome is unknown.
LF gain was a significant predictor of 90-day categorized NIHSS, functional independency and mortality after adjusting for age, recanalization success, and favorable ASPECTS/PC-ASPECTS before EVT. However, with the added knowledge at 24-hour FU, trichotomized NIHSS seemed to confound the effect of average LF gain on predictions of independency and categorized NIHSS at 90-day FU. Average LF gain could still be impactful in the 24-hour prediction model especially when 24-hour NIHSS is unknown (e.g., due to prolonged intubation).
Strengths and limitations
In this study, we applied one of the most frequently used and standardized methods to evaluate dCA, though with the novelty adjustment of using NIRS signals as both input and output. There are certain advantages to this approach. First, the setup of NIRS does not require expertise and takes less than 60 seconds enabling dCA studies in hyperacute settings before recanalization therapies without compromising patient safety. Second, only very few patients cannot be examined with NIRS due to extracranial hematomas or poor signal quality. Third, no scaling is required to adjust for different measurement modalities, removing any uncertainty concerning the scaling method. Lastly, while extracranial tissue contributes to the long-distance signal, the applied model accounts for that because LFO in the extracerebral tissue can reasonably be assumed to be static and without regional differences. This assumption was upheld in the short-distance channel TFA. Also, there is absolute certainty of intracerebral signal component due to contrast artifacts in all long-distance channels. Thus, we are confident that any changes in LF gain or phase shift are due to intracerebral changes, i.e. dynamic autoregulation.
Other approaches to study dCA exclusively by NIRS have been applied previously. Interestingly, Becker et al. interestingly proposed to use short-distance channels as both regressors to the long-distance channel and as input to the TFA as surrogate of systemic perfusion.43 This approach shares some of the advantages of our method while also enabling regional dCA examination in both hemispheres. However, it also requires sufficient data quality of short-distance channels, which are more prone to detector saturation and ultimately hinder the incorporation of results from short-distance channel TFA as variables in our statistical models. Short-distance channels also have a small intracerebral signal component and cannot be thought of as purely representative of systemic perfusion.44 This was also evident in our data as contrast artifacts appeared concurrently in some short-distance channels, albeit at a much lower intensity than in long-distance channels. The intracerebral component with dCA changes could also explain some of the variance in our short-distance channel TFA. Although regression of short-distance channels improves intracerebral sensitivity,45 regression methods have substantial heterogeneity and are tainted by the trade-off between sacrificing intracerebral signals and accepting extracerebral contamination. Even NIRS oximeters with multi-distance channels (spatially resolved spectroscopy-NIRS) commercially available for clinical use have demonstrated extracerebral contamination and have the added uncertainties of propriety processing resulting in different values of cortical hemoglobin saturation46,47 which renders reproducibility across NIRS devices impossible. Also, NIRS oximeters could easily have overlooked contrast artifacts, which ensured an intracerebral signal component but should be eliminated from TFA segments.
Another NIRS-only approach to examine phase shift was shown by Elting et al. who applied TFA between OxyHb and deoxygenated hemoglobin (DeoxyHb) from single channels measurements with spatially resolved spectroscopy-NIRS. Accounting for transit time (HF phase shift), this approach produced similar phase shift to TFA between ABP and VMCA.48 This method was also feasible in ICU patients with neurologic disease,49 but has never been applied in acute settings nor in AIS patients alone. Whether such an approach can be reproduced with CW-NIRS remains unanswered. Here, we choose to analyze OxyHb alone as the signal quality is better than DeoxyHb especially in hyperacute settings.
Various limitations restrict the conclusions from this study. Interhemispheric TFA examines dCA relative to the contralateral hemisphere. While amplitude can be assessed per hemisphere (PSD) phase shift cannot. This leaves uncertainty as to whether phase shift is equally impaired or equally intact in patients with large-vessel occlusions when not doing simultaneous TFA between ABP and OxyHb on each hemisphere. Also, we had no control group of healthy subjects to compare with.
AIS patients including those receiving EVT are inherently heterogenous concerning stroke etiology, co-morbidities, imaging parameters, and treatment. In this study we chose to prioritize everyday generalizability and not exclude patients based on heterogeneities, but account for as many of them as possible in the statistical models instead. However, accounting for all heterogeneities is nearly impossible and requires a very large sample size especially for less common occurrences (e.g., etiologies, infarct locations, co-morbidities, concomitant treatment). One such example could be patients receiving sevoflurane to maintain general anesthesia which could affect dCA differently from patients receiving propofol.50 LFO amplitude was lower during EVT regardless of anesthetics, but gain was unaltered. Our results show uncertainty concerning sevoflurane effect on interhemispheric phase shift, and thus sevoflurane subjects remained in all analysis.
The main exclusion criteria from patients examined with NIRS was the absence of a steady state period mainly before recanalization. Reasons for this were a short procedural time to access the occlusion and attempt recanalization within exclusion period after intubation, difficulties to establish steady state, e.g. due to cardiovascular comorbidity as well as high or low sensitivity to anesthetics or vasopressors. Therefore, we believe no selection bias has been introduced by data segment selection.
We did not evaluate collateral circulation due to differences in standard-of-care imaging (CT-angiography, MRI time-of-flight angiography) and missing contralateral series from the digital subtraction angiography. Collateral circulation could certainly have an impact on both outcome and dCA.51,52 However, collaterals are better developed in patients with large-artery stenosis compared to patients with cardioembolic stroke etiology and we did not find any effect of etiology.
Results from the FU90-group should be interpreted with caution as non-FU patients did not attend mainly due to very poor functional or fatal outcome. Efforts were made to minimize this trend by offering FU visits at the rehabilitation center or at home visits despite COVID-19 restrictions. The exclusion of a significant portion of patients with the worst outcome may explain why associations between LF gain and outcome were not significant in the FU90-group.
In conclusion, interhemispheric TFA based on NIRS is a feasible method for investigating dCA in AIS patients before, during and after EVT. LF gain remained stable during EVT up to 90 days post-treatment and had associations to as well as predictive capabilities of long-term symptom severity, functional outcome, and mortality. Development of LF phase shifts depended on symptom severity and increased immediately after recanalization in AIS patient with lower NIHSS. Further studies involving concurrent dCA investigations using different modalities are warranted to assess dCA in the contralateral hemisphere and to validate interhemispheric TFA but will also advance our understanding of the relations between different TFA models for studying dCA.
Declaration of interests
The authors declare that there is no conflict of interest.
Funding
The study was supported by Simon Fougner Hartmanns Family Foundation, Gangsted Foundation, Sophus Jacobsen Foundation and Rigshospitalet’s Research Foundation.
Supplemental Material
Supplemental Methods
Figure S1
Tables S1–S15
References no. 53-63
Author contributions
AVH and HKI conceived and designed the study with contributions from TCT, HGB, GB, CS, and HWS. AVH acquired the data with major contributions from TGL as well as TCT, HGB, GB, CS and KH. AVH performed data analysis with contributions from TGL and HKI. All authors contributed to interpretation of results. AVH drafted the article. All authors revised the article critically and approved the final version.
Data availability statement
All processing code is open source available. NIRS data segments can be provided upon reasonable request. Individual patient data are identifiable, and availability requires approval of both centers as well as local ethics committees.
Acknowledgements
We gratefully appreciate all the patients for their participation and thank all members of staff for their help in patient recruitment and executing the study.
Non-standard Abbreviations and Acronyms
- ABP
- Arterial Blood Pressure
- ACA
- Anterior Cerebral Artery
- AIS
- Acute Ischemic Stroke
- ASPECTS
- Alberta Stroke Program Early CT Score
- BA
- Basilar artery.
- CW-NIRS
- Continuous Wave Near-Infrared Spectroscopy
- dCA
- Dynamic Cerebral Autoregulation
- DeoxyHb
- Deoxygenated Hemoglobin Concentration
- EVT
- Endovascular Treatment
- FU
- Follow-up
- FU24-group
- All participants with sufficient data at PRE, POST and 24-hour FU.
- FU90-group
- All patients with sufficient data at PRE, POST, 24-hour and 90-day FU.
- HF
- High-frequency (0.2-0.5 Hz)
- ICA
- Internal carotid artery.
- IVT
- Intravenous thrombolysis.
- LF
- Low-frequency (0.07-0.2 Hz)
- LFO
- Low-frequency oscillations
- LVO
- Large-vessel Occlusion
- M1
- First part of middle cerebral artery.
- M2
- Second part of middle cerebral artery.
- MCA
- Middle Cerebral Artery
- mRS
- modified Rankin Scale
- mTICI
- Modified Treatment In Cerebral Infarction
- NIHSS
- National Institutes of Health Stroke Scale
- NIRS
- Near-Infrared Spectroscopy
- OxyHb
- Oxygenated Hemoglobin Concentration
- PC-ASPECTS
- Posterior circulation Alberta Stroke Program Early CT Score
- PCA
- Posterior cerebral artery.
- POST
- Second time segment (after achieving final reperfusion status)
- PRE
- First time segment (after sedation but before any attempted revascularization)
- PSD
- Power spectral
- STROBE
- STrengthening the Reporting of OBservational studies in Epidemiology
- TFA
- Transfer Function Analysis
- TIA
- Transient ischemic attack
- VLF
- Very low-frequency (0.02-0.07 Hz).
- VMCA
- Flow velocity in the Middle Cerebral Artery