MRI of spontaneous fluctuations after acute cerebral ischemia in nonhuman primates

A novel approach for assessing hypoperfusion in stroke using spatial independent component analysis of resting-state fMRI

Individualized treatment of acute stroke depends on the timely detection of ischemia and potentially salvageable tissue in the brain. Using functional MRI (fMRI), it is possible to characterize cerebral blood flow from blood‐oxygen‐level‐dependent (BOLD) signals without the administration of exogenous contrast agents. In this study, we applied spatial independent component analysis to resting‐state fMRI data of 37 stroke patients scanned within 24 hr of symptom onset, 17 of whom received follow‐up scans the next day. Our analysis revealed “Hypoperfusion spatially‐Independent Components” (HICs) whose spatial patterns of BOLD signal resembled regions of delayed perfusion depicted by dynamic susceptibility contrast MRI. These HICs were detected even in the presence of excessive patient motion, and disappeared following successful tissue reperfusion. The unique spatial and temporal features of HICs allowed them to be distinguished with high accuracy from other components in a user‐independent manner (area under the curve = 0.93, balanced accuracy = 0.90, sensitivity = 1.00, and specificity = 0.85). Our study therefore presents a new, noninvasive method for assessing blood flow in acute stroke that minimizes interpretative subjectivity and is robust to severe patient motion.

The Effect of Scan Length on the Assessment of BOLD Delay in Ischemic Stroke

Objectives: To evaluate the impact of resting-state functional MRI scan length on the diagnostic accuracy, image quality and lesion volume estimation of BOLD delay maps used for brain perfusion assessment in acute ischemic stroke.

Methods: Sixty-three acute ischemic stroke patients received a 340 s resting-state functional MRI within 24 h of stroke symptom onset. BOLD delay maps were calculated from the full scan and four shortened versions (68 s, 136 s, 204 s, 272 s). The BOLD delay lesions on these maps were compared in terms of spatial overlap and volumetric agreement with the lesions derived from the full scans and with time-to-maximum (Tmax) lesions derived from DSC-MRI in a subset of patients (n = 10). In addition, the interpretability and quality of these maps were compared across different scan lengths using mixed models.

Results: Shortened BOLD delay scans showed a small volumetric bias (ranging from 0.05 to 5.3 mL; between a 0.13% volumetric underestimation and a 7.7% overestimation relative to the mean of the volumes, depending on scan length) compared to the full scan. Decreased scan length was associated with decreased spatial overlap with both the BOLD delay lesions derived from the full scans and with Tmax lesions. Only the two shortest scan lengths (68 and 136 s) were associated with substantially decreased interpretability, decreased structure clarity, and increased noisiness of BOLD delay maps.

Conclusions: BOLD delay maps derived from resting-state fMRI scans lasting 272 and 204 s provide sufficient diagnostic quality and adequate assessment of perfusion lesion volumes. Such shortened scans may be helpful in situations where quick clinical decisions need to be made.

Non-invasive monitoring of longitudinal changes in cerebral hemodynamics in acute ischemic stroke using BOLD signal delay

Relative delays in blood-oxygen-level-dependent (BOLD) signal oscillations can be used to assess cerebral perfusion without using contrast agents. However, little is currently known about the utility of this method in detecting clinically relevant perfusion changes over time. We investigated the relationship between longitudinal BOLD delay changes, vessel recanalization, and reperfusion in 15 acute stroke patients with vessel occlusion examined within 24 h of symptom onset (D0) and one day later (D1). We created BOLD delay maps using time shift analysis of resting-state functional MRI data and quantified perfusion lesion volume changes (using the D1/D0 volume ratio) and severity changes (using a linear mixed model) over time. Between baseline and follow-up, BOLD delay lesions shrank (median D1/D0 ratio = 0.2, IQR = 0.03-0.7) and BOLD delay severity decreased (b = -4.4 s) in patients with recanalization, whereas they grew (median D1/D0 ratio = 1.47, IQR = 1.1-1.7) and became more severe (b = 4.3 s) in patients with persistent vessel occlusion. Clinically relevant changes in cerebral perfusion in early stroke can be detected using BOLD delay, making this non-invasive method a promising option for detecting tissue at risk of infarction and monitoring stroke patients following recanalization therapy.

Dynamic Alterations in Spontaneous Neural Activity in Multiple Brain Networks in Subacute Stroke Patients: A Resting-State fMRI Study

Objective: To examine whether subacute stroke patients would exhibit abnormal dynamic characteristics of brain activity relative to healthy controls (HC) and to investigate whether the altered dynamic regional indexes were associated with clinical behavior in stroke patients. Methods: The dynamic amplitude of low-frequency fluctuations (dALFF) and dynamic regional homogeneity (dReHo) in 42 subacute stroke patients and 55 healthy controls were compared. Correlation analyses between dALFF and dReHo in regions showing significant intergroup differences and clinical scores (i.e., the National Institutes of Health Stroke Scale, Fugl-Meyer assessment and lesion volume size) were conducted in stroke patients. Receiver operating characteristic (ROC) curve analysis was used to determine the potential value of altered dynamic regional indexes to identify stroke patients. Results: Significantly dALFF in the bilateral cerebellum posterior lobe (CPL), ipsilesional superior parietal lobe, ipsilesional inferior temporal gyrus (ITG), the midline supplementary motor area (SMA), ipsilesional putamen and lentiform nucleus were detected in stroke patients compared to HC. Relative to the HC group, the stroke patients showed significant differences in dReHo in the contralesional rectal gyrus, contralesional ITG, contralesional pons, ipsilesional middle frontal gyrus (MFG). Significant correlations between dALFF variability in midline SMA and Fugl-Meyer assessment (FMA) scores or between dReHo variability in the ipsilesional MFG and FMA scores were detected in stroke patients. Furthermore, the ROC curve revealed that dynamic ALFF at SMA and ReHo at ipsilesional MFG might have the potential to distinguish stroke patients. Conclusion: The pattern of intrinsic brain activity variability is altered in stroke patients compared with HC, and dynamic ALFF/ReHo might be potential tools to assess stroke patients' motor function.

Power Spectral Density Analysis of Long-Term Motor Recovery in Patients With Subacute Stroke

BACKGROUND

Prognostic measures of long-term motor recovery are important in patients with stroke presenting with severe hemiplegia.

OBJECTIVE

We aimed to investigate whether initial power spectral density (PSD) analysis of resting-state functional magnetic resonance (fMRI) data can provide a sensitive prognostic predictor in patients with subacute stroke with severe hand disability.

METHODS

Twelve patients with good recovery, 14 patients with poor recovery, and 12 healthy subjects were included. PSD analysis was performed using resting-state fMRI data. Contralesional and ipsilesional PSD in the motor cortex were measured. Pearson correlation analysis was performed to assess a possible association between the difference in ipsilesional versus contralesional PSD and motor outcomes. A receiver operating characteristic (ROC) curve was constructed to estimate the discriminative value of the difference between the ipsilesional PSD and the contralesional PSD for good versus poor recovery.

RESULTS

There were no differences in PSD between the contralesional and ipsilesional hemispheres in the good recovery group ( P = .77). In contrast, there were significant differences in PSD between the 2 hemispheres in the poor recovery group ( P = .07). The difference in PSD between the 2 hemispheres had a positive correlation with post Brunnstrom stage scores. ROC analysis showed that the difference in PSD between the 2 hemispheres was sensitive in discriminating good versus poor recovery.

CONCLUSION

The present study suggests that PSD in the motor cortex may be a sensitive predictor of late-onset motor recovery following stroke.

Assessment and modulation of resting-state neural networks after stroke

PURPOSE OF REVIEW

Stroke is a major cause of disability; however, most patients experience spontaneous partial recovery of functions in subacute to chronic phases. Poststroke loss and recovery of functions have been increasingly correlated with brain-wide alterations in the connectivity of neural networks, which is described in this review. Elucidation of the mechanisms of functional brain remodeling could reveal targets and strategies for more effective neurorehabilitation.

RECENT FINDINGS

Data from recent resting-state functional MRI, electroencephalography, magnetoencephalography, and optical imaging studies in patients and animal models have demonstrated that loss of function after stroke is closely associated with disrupted connectivity in large-scale networks beyond the lesion territory. Restoration of functional connectivity in the surviving networks appears to be critical for functional recovery, and this may be promoted with specific therapeutic strategies, such as robot-assisted training and noninvasive brain stimulation. The adaptability of functional networks relies on the structural integrity of neuronal pathways, but the relationship between the two remains incompletely understood. Furthermore, disturbed neurovascular coupling after stroke can confound hemodynamically based measurements of functional connectivity.

SUMMARY

Identification of key network processes in adaptive brain plasticity can aid in the prediction of functional outcome and the development of therapeutic interventions to support and promote recovery after stroke.

Altered resting-state FMRI signals in acute stroke patients with ischemic penumbra

Background

Identifying the ischemic penumbra in acute stroke subjects is important for the clinical decision making process. The aim of this study was to use resting-state functional magnetic resonance singal (fMRI) to investigate the change in the amplitude of low-frequency fluctuations (ALFF) of these subjects in three different subsections of acute stroke regions: the infarct core tissue, the penumbra tissue, and the normal brain tissue. Another aim of this study was to test the feasilbility of consistently detecting the penumbra region of the brain through ALFF analysis.

Methods

Sixteen subjects with first-ever acute ischemic stroke were scanned within 27 hours of the onset of stroke using magnetic resonance imaging. The core of infarct regions and penumbra regions were determined by diffusion and perfusion-weighted imaging respectively. The ALFF were measured from resting-state blood oxygen level dependent (BOLD) fMRI scans. The averaged relative ALFF value of each regions were correlated with the time after the onset of stroke.

Results

Relative ALFF values were significantly different in the infarct core tissue, penumbra tissue and normal brain tissue. The locations of lesions in the ALFF maps did not match perfectly with diffusion and perfusion-weighted imagings; however, these maps provide a contrast that can be used to differentiate between penumbra brain tissue and normal brain tissue. Significant correlations between time after stroke onset and the relative ALFF values were present in the penumbra tissue but not in the infarct core and normal brain tissue.

Conclusion

Preliminary results from this study suggest that the ALFF reflects the underlying neurovascular activity and has a great potential to estimate the brain tissue viability after ischemia. Results also show that the ALFF may contribute to acute stroke imaging for thrombolytic or neuroprotective therapies.

Spontaneous BOLD signal fluctuations in young healthy subjects and elderly patients with chronic kidney disease

Spontaneous fluctuations in blood oxygenation level-dependent (BOLD) images are the basis of resting-state fMRI and frequently used for functional connectivity studies. However, there may be intrinsic information in the amplitudes of these fluctuations. We investigated the possibility of using the amplitude of spontaneous BOLD signal fluctuations as a biomarker for cerebral vasomotor reactivity.

We compared the coefficient of variation (CV) of the time series (defined as the temporal standard deviation of the time series divided by the mean signal intensity) in two populations: 1) Ten young healthy adults and 2) Ten hypertensive elderly subjects with chronic kidney disease (CKD).

We found a statistically significant increase (P<0.01) in the CV values for the CKD patients compared with the young healthy adults in both gray matter (GM) and white matter (WM). The difference was independent of the exact segmentation method, became more significant after correcting for physiological signals using RETROICOR, and mainly arose from very low frequency components of the BOLD signal fluctuation (f<0.025 Hz). Furthermore, there was a strong relationship between WM and GM signal fluctuation CV's (R² = 0.87) in individuals, with a ratio of about 1∶3.

These results suggest that amplitude of the spontaneous BOLD signal fluctuations may be used to assess the cerebrovascular reactivity mechanisms and provide valuable information about variations with age and different disease states.

Noncontrast mapping of arterial delay and functional connectivity using resting-state functional MRI: a study in Moyamoya patients

PURPOSE

To investigate if delays in resting-state spontaneous fluctuations of the BOLD (sfBOLD) signal can be used to create maps similar to time-to-maximum of the residue function (Tmax) in Moyamoya patients and to determine whether sfBOLD delays affect the results of brain connectivity mapping.

MATERIALS AND METHODS

Ten patients were scanned at 3 Tesla using a gradient-echo echo planar imaging sequence for sfBOLD imaging. Cross correlation analysis was performed between each brain voxel signal and a reference signal comprised of either the superior sagittal sinus (SSS) or whole brain (WB) average time course. sfBOLD delay maps were created based on the time shift necessary to maximize the correlation coefficient, and compared with dynamic susceptibility contrast Tmax maps. Standard and time-shifted resting-state BOLD connectivity analyses of the default mode network were compared.

RESULTS

Good linear correlations were found between sfBOLD delays and Tmax using the SSS as reference (r(2)  = 0.8, slope = 1.4, intercept = -4.6) or WB (r(2)  = 0.7, slope = 0.8, intercept = -3.2). New nodes of connectivity were found in delayed regions when accounting for delays in the analysis.

CONCLUSION

Resting-state sfBOLD imaging can create delay maps similar to Tmax maps without the use of contrast agents in Moyamoya patients. Accounting for these delays may affect the results of functional connectivity maps.

A review of current imaging methods used in stroke research

Stroke is a serious healthcare problem with high mortality and long-term disability. However, to date, our ability to prevent and cure stroke remains limited. One important goal in stroke research is to identify the extent and location of lesion for treatment. In addition, accurately differentiating salvageable tissue from infarct and evaluating therapeutic efficacies are indispensible. These objectives could potentially be met with the assistance of modern neuroimaging techniques. This paper reviews current imaging methods commonly used in ischemic stroke research. These methods include positron emission tomography, computed tomography, T1 MRI, T2 MRI, diffusion and perfusion MRI, diffusion tensor imaging, blood-brain barrier permeability MRI, pH-weighted MRI, and functional MRI.

Multimodal MRI of nonhuman primate stroke

Stroke is the fourth leading cause of death. Despite decades of research, no neuroprotective drug has proven to be effective clinically. One widely accepted view to account for this negative outcome is that the rodent stroke model simply does not adequately reflect the complexity of human stroke. Recent failures of several high-profile neuroprotective drugs for stroke treatment in phase III clinical trials further underscore the importance of developing adequate animal models for stroke research. The brain organization and vascular circuitry of nonhuman primates (NHPs) are more homologous with humans than the widely used rodent for stroke modeling. The Stroke Therapy Academic Industry Roundtable, a national committee commissioned by the American Heart Association, recommended that clinically relevant NHP stroke models be established for developing and assessing neuroprotective drugs. The aim of this article is to review the challenges and applications of magnetic resonance imaging studies of NHP stroke models.

Imaging Stroke Evolution after Middle Cerebral Artery Occlusion in Non-human Primates

This article reviews imaging approaches applied to the study of stroke in nonhuman primates. We briefly survey the various surgical and minimally invasive experimental stroke models in nonhuman primates, followed by a summary of studies using computed tomography, positron emission tomography and magnetic resonance imaging and spectroscopy to monitor stroke from the hyperacute phase (within minutes of the onset of cerebral ischemia) to the chronic phase (1 month and beyond).

Brain metabolites B1-corrected proton T1 mapping in the rhesus macaque at 3 T

The accuracy of metabolic quantification in MR spectroscopy is limited by the unknown radiofrequency field and T(1). To address both issues in proton ((1)H) MR spectroscopy, we obtained radiofrequency field-corrected T(1) maps of N-acetylaspartate, choline, and creatine in five healthy rhesus macaques at 3 T. For efficient use of the 4 hour experiment, we used a new three-point protocol that optimizes the precision of T(1) in three-dimensional (1)H-MR spectroscopy localization for extensive, approximately 30%, brain coverage at 0.6 x 0.6 x 0.5 cm(3) = 180-microL spatial resolution. The resulting mean T(1)s in 700 voxels were N-acetylaspartate = 1232 +/- 44, creatine = 1238 +/- 23 and choline = 1107 +/- 56 ms (mean +/- standard error of the mean). Their histograms from all 140 voxels in each animal were similar in position and shape, characterized by standard errors of the mean of the full width at half maximum divided by their means of better than 8%. Regional gray matter N-acetylaspartate, choline, and creatine T(1)s (1333 +/- 43, 1265 +/- 52, and 1131 +/- 28 ms) were 5-10% longer than white matter: 1188 +/- 34, 1201 +/- 24, and 1082 +/- 50 ms (statistically significant for the N-acetylaspartate only), all within 10% of the corresponding published values in the human brain.

Metabolite proton T(2) mapping in the healthy rhesus macaque brain at 3 T

The structure and metabolism of the rhesus macaque brain, an advanced model for neurologic diseases and their treatment response, is often studied noninvasively with MRI and (1)H-MR spectroscopy. Due to the shorter transverse relaxation time (T(2)) at the higher magnetic fields these studies favor, the echo times used in (1)H-MR spectroscopy subject the metabolites to unknown T(2) weighting, decreasing the accuracy of quantification which is key for inter- and intra-animal comparisons. To establish the "baseline" (healthy animal) T(2) values, we mapped them for the three main metabolites' T(2)s at 3 T in four healthy rhesus macaques and tested the hypotheses that their mean values are similar (i) among animals; and (ii) to analogs regions in the human brain. This was done with three-dimensional multivoxel (1)H-MR spectroscopy at (0.6 x 0.6 x 0.5 cm)(3) = 180 microL spatial resolution over a 4.2 x 3.0 x 2.0 = 25 cm(3) ( approximately 30%) of the macaque brain in a two-point protocol that optimizes T(2) precision per unit time. The estimated T(2)s in several gray and white matter regions are all within 10% of those reported in the human brain (mean +/- standard error of the mean): N-acetylaspartate = 316 +/- 7, creatine = 177 +/- 3, and choline = 264 +/- 9 ms, with no statistically significant gray versus white matter differences.

Mapping Transient Hyperventilation Induced Alterations with Estimates of the Multi-Scale Dynamics of BOLD Signal

Temporal blood oxygen level dependent (BOLD) contrast signals in functional MRI during rest may be characterized by power spectral distribution (PSD) trends of the form 1/f(alpha). Trends with 1/f characteristics comprise fractal properties with repeating oscillation patterns in multiple time scales. Estimates of the fractal properties enable the quantification of phenomena that may otherwise be difficult to measure, such as transient, non-linear changes. In this study it was hypothesized that the fractal metrics of 1/f BOLD signal trends can map changes related to dynamic, multi-scale alterations in cerebral blood flow (CBF) after a transient hyperventilation challenge. Twenty-three normal adults were imaged in a resting-state before and after hyperventilation. Different variables (1/f trend constant alpha, fractal dimension D(f), and, Hurst exponent H) characterizing the trends were measured from BOLD signals. The results show that fractal metrics of the BOLD signal follow the fractional Gaussian noise model, even during the dynamic CBF change that follows hyperventilation. The most dominant effect on the fractal metrics was detected in grey matter, in line with previous hyperventilation vaso-reactivity studies. The alpha was able to differentiate also blood vessels from grey matter changes. D(f) was most sensitive to grey matter. H correlated with default mode network areas before hyperventilation but this pattern vanished after hyperventilation due to a global increase in H. In the future, resting-state fMRI combined with fractal metrics of the BOLD signal may be used for analyzing multi-scale alterations of cerebral blood flow.

Spontaneous low-frequency blood oxygenation level-dependent fluctuations and functional connectivity analysis of the 'resting' brain

Functional magnetic resonance imaging techniques using the blood oxygenation level-dependent (BOLD) contrast are widely used to map human brain function by relating local hemodynamic responses to neuronal stimuli compared to control conditions. There is increasing interest in spontaneous cerebral BOLD fluctuations that are prominent in the low-frequency range (<0.1 Hz) and show intriguing spatio-temporal correlations in functional networks. The nature of these signal fluctuations remains unclear, but there is accumulating evidence for a neural basis opening exciting new avenues to study human brain function and its connectivity at rest. Moreover, an increasing number of patient studies report disease-dependent variation in the amplitude and spatial coherence of low-frequency BOLD fluctuations (LFBF) that may afford greater diagnostic sensitivity and easier clinical applicability than standard fMRI. The main disadvantage of this emerging tool relates to physiological (respiratory, cardiac and vasomotion) and motion confounds that are challenging to disentangle requiring thorough preprocessing. Technical aspects of functional connectivity fMRI analysis and the neuroscientific potential of spontaneous LFBF in the default mode and other resting-state networks have been recently reviewed. This review will give an update on the current knowledge of the nature of LFBF, their relation to physiological confounds and potential for clinical diagnostic and pharmacological studies.