Relationship Between Changes in the Temporal Dynamics of the Blood-Oxygen-Level-Dependent Signal and Hypoperfusion in Acute Ischemic Stroke

Functional connectivity-hemodynamic (un)coupling changes in chronic mild brain injury are associated with mental health and neurocognitive indices: a resting state fMRI study

PURPOSE

To examine hemodynamic and functional connectivity alterations and their association with neurocognitive and mental health indices in patients with chronic mild traumatic brain injury (mTBI).

METHODS

Resting-state functional MRI (rs-fMRI) and neuropsychological assessment of 37 patients with chronic mTBI were performed. Intrinsic connectivity contrast (ICC) and time-shift analysis (TSA) of the rs-fMRI data allowed the assessment of regional hemodynamic and functional connectivity disturbances and their coupling (or uncoupling). Thirty-nine healthy age- and gender-matched participants were also examined.

RESULTS

Patients with chronic mTBI displayed hypoconnectivity in bilateral hippocampi and parahippocampal gyri and increased connectivity in parietal areas (right angular gyrus and left superior parietal lobule (SPL)). Slower perfusion (hemodynamic lag) in the left anterior hippocampus was associated with higher self-reported symptoms of depression (r =  - 0.53, p = .0006) and anxiety (r =  - 0.484, p = .002), while faster perfusion (hemodynamic lead) in the left SPL was associated with lower semantic fluency (r =  - 0.474, p = .002). Finally, functional coupling (high connectivity and hemodynamic lead) in the right anterior cingulate cortex (ACC)) was associated with lower performance on attention and visuomotor coordination (r =  - 0.50, p = .001), while dysfunctional coupling (low connectivity and hemodynamic lag) in the left ventral posterior cingulate cortex (PCC) and right SPL was associated with lower scores on immediate passage memory (r =  - 0.52, p = .001; r =  - 0.53, p = .0006, respectively). Uncoupling in the right extrastriate visual cortex and posterior middle temporal gyrus was negatively associated with cognitive flexibility (r =  - 0.50, p = .001).

CONCLUSION

Hemodynamic and functional connectivity differences, indicating neurovascular (un)coupling, may be linked to mental health and neurocognitive indices in patients with chronic mTBI.

Resting-State fMRI: Emerging Concepts for Future Clinical Application

Resting-state functional magnetic resonance imaging (rsfMRI) has been developed as a method of investigating spontaneous neural activity. Based on its low-frequency signal synchronization, rsfMRI has made it possible to identify multiple macroscopic structures termed resting-state networks (RSNs) on a single scan of less than 10 minutes. It is easy to implement even in clinical practice, in which assigning tasks to patients can be challenging. These advantages have accelerated the adoption and growth of rsfMRI. Recently, studies on the global rsfMRI signal have attracted increasing attention. Because it primarily arises from physiological events, less attention has hitherto been paid to the global signal than to the local network (i.e., RSN) component. However, the global signal is not a mere nuisance or a subsidiary component. On the contrary, it is quantitatively the dominant component that accounts for most of the variance in the rsfMRI signal throughout the brain and provides rich information on local hemodynamics that can serve as an individual-level diagnostic biomarker. Moreover, spatiotemporal analyses of the global signal have revealed that it is closely and fundamentally associated with the organization of RSNs, thus challenging the basic assumptions made in conventional rsfMRI analyses and views on RSNs. This review introduces new concepts emerging from rsfMRI spatiotemporal analyses focusing on the global signal and discusses how they may contribute to future clinical medicine. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

Resting-state networks representation of the global phenomena

Resting-state functional magnetic resonance imaging (rsfMRI) has been widely applied to investigate spontaneous neural activity, often based on its macroscopic organization that is termed resting-state networks (RSNs). Although the neurophysiological mechanisms underlying the RSN organization remain largely unknown, accumulating evidence points to a substantial contribution from the global signals to their structured synchronization. This study further explored the phenomenon by taking advantage of the inter- and intra-subject variations of the time delay and correlation coefficient of the signal timeseries in each region using the global mean signal as the reference signal. Consistent with the hypothesis based on the empirical and theoretical findings, the time lag and correlation, which have consistently been proven to represent local hemodynamic status, were shown to organize networks equivalent to RSNs. The results not only provide further evidence that the local hemodynamic status could be the direct source of the RSNs' spatial patterns but also explain how the regional variations in the hemodynamics, combined with the changes in the global events' power spectrum, lead to the observations. While the findings pose challenges to interpretations of rsfMRI studies, they further support the view that rsfMRI can offer detailed information related to global neurophysiological phenomena as well as local hemodynamics that would have great potential as biomarkers.

Distinct hemodynamic and functional connectivity features of fatigue in clinically isolated syndrome and multiple sclerosis: accounting for the confounding effect of concurrent depression symptoms

PURPOSE

This study aims to identify common and distinct hemodynamic and functional connectivity (FC) features for self-rated fatigue and depression symptoms in patients with clinically isolated syndrome (CIS) and relapsing-remitting multiple sclerosis (RR-MS).

METHODS

Twenty-four CIS, 29 RR-MS patients, and 39 healthy volunteers were examined using resting-state fMRI (rs-fMRI) to obtain whole-brain maps of (i) hemodynamic response patterns (through time shift analysis), (ii) FC (via intrinsic connectivity contrast maps), and (iii) coupling between hemodynamic response patterns and FC. Each regional map was correlated with fatigue scores, controlling for depression, and with depression scores, controlling for fatigue.

RESULTS

In CIS patients, the severity of fatigue was associated with accelerated hemodynamic response in the insula, hyperconnectivity of the superior frontal gyrus, and evidence of reduced hemodynamics-FC coupling in the left amygdala. In contrast, depression severity was associated with accelerated hemodynamic response in the right limbic temporal pole, hypoconnectivity of the anterior cingulate gyrus, and increased hemodynamics-FC coupling in the left amygdala. In RR-MS patients, fatigue was associated with accelerated hemodynamic response in the insula and medial superior frontal cortex, increased functional role of the left amygdala, and hypoconnectivity of the dorsal orbitofrontal cortex, while depression symptom severity was linked to delayed hemodynamic response in the medial superior frontal gyrus; hypoconnectivity of the insula, ventromedial thalamus, dorsolateral prefrontal cortex, and posterior cingulate; and decreased hemodynamics-FC coupling of the medial orbitofrontal cortex.

CONCLUSION

There are distinct FC and hemodynamic responses, as well as different magnitude and topography of hemodynamic connectivity coupling, associated with fatigue and depression in early and later stages of MS.

Cerebrovascular Reactivity Has Negligible Contribution to Hemodynamic Lag After Stroke: Implications for Functional Magnetic Resonance Imaging Studies

BACKGROUND

Functional magnetic resonance imaging (fMRI) is ubiquitously used to study poststroke recovery. However, the fMRI-derived hemodynamic responses are vulnerable to vascular insult which can result in reduced magnitude and temporal delays (lag) in the hemodynamic response function (HRF). The cause of HRF lag remains controversial, and a better understanding of it is required to ensure accurate interpretation of poststroke fMRI studies. In this longitudinal study, we investigate the relationship between hemodynamic lag and cerebrovascular reactivity (CVR) following stroke.

METHODS

Voxel-wise lag maps were calculated relative to a mean gray matter reference signal for 27 healthy controls and 59 patients with stroke across 2 time points (≈2 weeks and ≈4 months poststroke) and 2 conditions: resting-state and breath-holding. The breath-holding condition was additionally used to calculate CVR in response to hypercapnia. HRF lag was computed for both conditions across tissue compartments: lesion, perilesional tissue, unaffected tissue of the lesioned hemisphere, and their homolog regions in the unaffected hemisphere. CVR and lag maps were correlated. Group, condition, and time effects were assessed using ANOVA analyses.

RESULTS

Compared with the average gray matter signal, a relative hemodynamic lead was observed in the primary sensorimotor cortices in resting-state and bilateral inferior parietal cortices in the breath-holding condition. Whole-brain hemodynamic lag was significantly correlated across conditions irrespective of group, with regional differences across conditions suggestive of a neural network pattern. Patients showed relative lag in the lesioned hemisphere which significantly reduced over time. Breath-hold derived lag and CVR had no significant voxel-wise correlation in controls, or patients within the lesioned hemisphere or the homologous regions of the lesion and perilesional tissue in the right hemisphere (mean r<0.1).

CONCLUSIONS

The contribution of altered CVR to HRF lag was negligible. We suggest that HRF lag is largely independent of CVR, and could partly reflect intrinsic neural network dynamics among other factors.

Reproducibility of cerebral perfusion measurements using BOLD delay

BOLD delay is an emerging, noninvasive method for assessing cerebral perfusion that does not require the use of intravenous contrast agents and is thus particularly suited for longitudinal monitoring. In this study, we assess the reproducibility of BOLD delay using data from 136 subjects with normal cerebral perfusion scanned on two separate occasions with scanners, sequence parameters, and intervals between scans varying between subjects. The effects of various factors on the reproducibility of BOLD delay, defined here as the differences in BOLD delay values between the scanning sessions, were investigated using a linear mixed model. Reproducibility was additionally assessed using the intraclass correlation coefficient of BOLD delay between sessions. Reproducibility was highest in the posterior cerebral artery territory. The mean BOLD delay test-retest difference after accounting for the aforementioned factors was 1.2 s (95% CI = 1.0 to 1.4 s). Overall, BOLD delay shows good reproducibility, but care should be taken when interpreting longitudinal BOLD delay changes that are either very small or are located in certain brain regions.

Reliability and Sensitivity to Alterered Hemodynamics Measured with Resting-state fMRI Metrics: Comparison with 123I-IMP SPECT

Blood oxygenation level-dependent (BOLD) contrast is sensitive to local hemodynamic changes and thus is applicable to imaging perfusion or vascular reactivity. However, knowledge about its measurement characteristics compared to reference standard perfusion imaging is limited. This study longitudinally evaluated perfusion in patients with steno-occlusive disease using resting-state functional MRI (rsfMRI) acquired before and within nine days of anterior circulation revascularization in patients with large cerebral artery steno-occlusive diseases. The reliability and sensitivity to longitudinal changes of rsfMRI temporal correlation (Rc) and time delay (TDc) relative to the cerebellar signal were examined voxel-wise in comparison with single-photon emission CT (SPECT) cerebral blood flow (CBF) using the within-subject standard deviation (Sw) and intraclass correlation coefficients (ICCs). For statistical comparisons, the standard deviation (SD) of longitudinal changes within the cerebellum, the number of voxels with significant changes in the left middle cerebral artery territory ipsilateral to surgery, and their average changes relative to the cerebellar SD were evaluated. The test-retest reliability of the fMRI metrics was also similarly evaluated using the human connectome project (HCP) healthy young adult dataset. The test-retest time interval was 31 ± 18 days. Test-retest reliability was significantly higher for SPECT (cerebellar SD: -2.59 ± 0.20) than for fMRI metrics (cerebellar SD: Rc, -2.34 ± 0.24, p = 0.04; TDc, -2.19 ± 0.21, p = 0.003). Sensitivity to postoperative changes, which was evaluated as the number of voxels, was significantly higher for fMRI TDc (8.78 ± 0.72) than for Rc (7.42 ± 1.48, p = 0.03) or SPECT CBF (6.88 ± 0.67, p < 0.001). The ratio between the average Rc, TDc, and SPECT CBF changes within the left MCA target region and cerebellar SD was also significantly higher for fMRI TDc (1.21 ± 0.79) than Rc (0.48 ± 0.94, p = 0.006) or SPECT CBF (0.23 ± 0.57, p = 0.001). The measurement variability of time delay was also larger than that of temporal correlation in HCP data within the cerebellum (t = -8.7, p < 0.001) or in the whole-brain (t = -27.4, p < 0.001) gray matter. These data suggest that fMRI time delay is more sensitive to the hemodynamic changes than SPECT CBF, although the reliability is lower. The implication for fMRI connectivity studies is that temporal correlation can be significantly decreased due to altered hemodynamics, even in cases with normal CBF.

Toward an integrative neurovascular framework for studying brain networks

Brain functional connectivity based on the measure of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals has become one of the most widely used measurements in human neuroimaging. However, the nature of the functional networks revealed by BOLD fMRI can be ambiguous, as highlighted by a recent series of experiments that have suggested that typical resting-state networks can be replicated from purely vascular or physiologically driven BOLD signals. After going through a brief review of the key concepts of brain network analysis, we explore how the vascular and neuronal systems interact to give rise to the brain functional networks measured with BOLD fMRI. This leads us to emphasize a view of the vascular network not only as a confounding element in fMRI but also as a functionally relevant system that is entangled with the neuronal network. To study the vascular and neuronal underpinnings of BOLD functional connectivity, we consider a combination of methodological avenues based on multiscale and multimodal optical imaging in mice, used in combination with computational models that allow the integration of vascular information to explain functional connectivity.

Physiological instability is linked to mortality in primary central nervous system lymphoma: A case-control fMRI study

Primary central nervous system lymphoma (PCNSL) is an aggressive brain disease where lymphocytes invade along perivascular spaces of arteries and veins. The invasion markedly changes (peri)vascular structures but its effect on physiological brain pulsations has not been previously studied. Using physiological magnetic resonance encephalography (MREG_BOLD ) scanning, this study aims to quantify the extent to which (peri)vascular PCNSL involvement alters the stability of physiological brain pulsations mediated by cerebral vasculature. Clinical implications and relevance were explored. In this study, 21 PCNSL patients (median 67y; 38% females) and 30 healthy age-matched controls (median 63y; 73% females) were scanned for MREG_BOLD signal during 2018-2021. Motion effects were removed. Voxel-by-voxel Coefficient of Variation (CV) maps of MREG_BOLD signal was calculated to examine the stability of physiological brain pulsations. Group-level differences in CV were examined using nonparametric covariate-adjusted tests. Subject-level CV alterations were examined against control population Z-score maps wherein clusters of increased CV values were detected. Spatial distributions of clusters and findings from routine clinical neuroimaging were compared [contrast-enhanced, diffusion-weighted, fluid-attenuated inversion recovery (FLAIR) data]. Whole-brain mean CV was linked to short-term mortality with 100% sensitivity and 100% specificity, as all deceased patients revealed higher values (n = 5, median 0.055) than surviving patients (n = 16, median 0.028) (p < .0001). After adjusting for medication, head motion, and age, patients revealed higher CV values (group median 0.035) than healthy controls (group median 0.024) around arterial territories (p ≤ .001). Abnormal clusters (median 1.10 × 10⁵ mm³ ) extended spatially beyond FLAIR lesions (median 0.62 × 10⁵ mm³ ) with differences in volumes (p = .0055).

A noninvasive and comprehensive method for continuous assessment of cerebral blood flow pulsation based on magnetic induction phase shift

Cerebral blood flow (CBF) monitoring is of great significance for treating and preventing strokes. However, there has not been a fully accepted method targeting continuous assessment in clinical practice. In this work, we built a noninvasive continuous assessment system for cerebral blood flow pulsation (CBFP) that is based on magnetic induction phase shift (MIPS) technology and designed a physical model of the middle cerebral artery (MCA). Physical experiments were carried out through different simulations of CBF states. Four healthy volunteers were enrolled to perform the MIPS and ECG synchronously monitoring trials. Then, the components of MIPS related to the blood supply level and CBFP were investigated by signal analysis in time and frequency domain, wavelet decomposition and band-pass filtering. The results show that the time-domain baseline of MIPS increases with blood supply level. A pulse signal was identified in the spectrum (0.2-2 Hz in 200-2,000 ml/h groups, respectively) of MIPS when the simulated blood flow rate was not zero. The pulsation frequency with different simulated blood flow rates is the same as the squeezing frequency of the feeding pump. Similar to pulse waves, the MIPS signals on four healthy volunteers all had periodic change trends with obvious peaks and valleys. Its frequency is close to that of the ECG signal and there is a certain time delay between them. These results indicate that the CBFP component can effectively be extracted from MIPS, through which different blood supply levels can be distinguished. This method has the potential to become a new solution for non-invasive and comprehensive monitoring of CBFP.

Effects of Thresholding on Voxel-Wise Correspondence of Breath-Hold and Resting-State Maps of Cerebrovascular Reactivity

Functional magnetic resonance imaging for presurgical brain mapping enables neurosurgeons to identify viable tissue near a site of operable pathology which might be at risk of surgery-induced damage. However, focal brain pathology (e.g., tumors) may selectively disrupt neurovascular coupling while leaving the underlying neurons functionally intact. Such neurovascular uncoupling can result in false negatives on brain activation maps thereby compromising their use for surgical planning. One way to detect potential neurovascular uncoupling is to map cerebrovascular reactivity using either an active breath-hold challenge or a passive resting-state scan. The equivalence of these two methods has yet to be fully established, especially at a voxel level of resolution. To quantitatively compare breath-hold and resting-state maps of cerebrovascular reactivity, we first identified threshold settings that optimized coverage of gray matter while minimizing false responses in white matter. When so optimized, the resting-state metric had moderately better gray matter coverage and specificity. We then assessed the spatial correspondence between the two metrics within cortical gray matter, again, across a wide range of thresholds. Optimal spatial correspondence was strongly dependent on threshold settings which if improperly set tended to produce statistically biased maps. When optimized, the two CVR maps did have moderately good correspondence with each other (mean accuracy of 73.6%). Our results show that while the breath-hold and resting-state maps may appear qualitatively similar they are not quantitatively identical at a voxel level of resolution.

Total perfusion-diffusion mismatch detected using resting-state functional MRI

Total perfusion-diffusion mismatch is a well-recognised phenomenon in patients with acute ischaemic stroke. We describe a case of total perfusion-diffusion mismatch detected using an emerging contrast-agent-free perfusion imaging technique in a young patient with acute cerebellar stroke.

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.

Anxiety and depression severity in neuropsychiatric SLE are associated with perfusion and functional connectivity changes of the frontolimbic neural circuit: a resting-state f(unctional) MRI study

OBJECTIVE

To examine the hypothesis that perfusion and functional connectivity disturbances in brain areas implicated in emotional processing are linked to emotion-related symptoms in neuropsychiatric SLE (NPSLE).

METHODS

Resting-state fMRI (rs-fMRI) was performed and anxiety and/or depression symptoms were assessed in 32 patients with NPSLE and 18 healthy controls (HC). Whole-brain time-shift analysis (TSA) maps, voxel-wise global connectivity (assessed through intrinsic connectivity contrast (ICC)) and within-network connectivity were estimated and submitted to one-sample t-tests. Subgroup differences (high vs low anxiety and high vs low depression symptoms) were assessed using independent-samples t-tests. In the total group, associations between anxiety (controlling for depression) or depression symptoms (controlling for anxiety) and regional TSA or ICC metrics were also assessed.

RESULTS

Elevated anxiety symptoms in patients with NPSLE were distinctly associated with relatively faster haemodynamic response (haemodynamic lead) in the right amygdala, relatively lower intrinsic connectivity of orbital dlPFC, and relatively lower bidirectional connectivity between dlPFC and vmPFC combined with relatively higher bidirectional connectivity between ACC and amygdala. Elevated depression symptoms in patients with NPSLE were distinctly associated with haemodynamic lead in vmPFC regions in both hemispheres (lateral and medial orbitofrontal cortex) combined with relatively lower intrinsic connectivity in the right medial orbitofrontal cortex. These measures failed to account for self-rated, milder depression symptoms in the HC group.

CONCLUSION

By using rs-fMRI, altered perfusion dynamics and functional connectivity was found in limbic and prefrontal brain regions in patients with NPSLE with severe anxiety and depression symptoms. Although these changes could not be directly attributed to NPSLE pathology, results offer new insights on the pathophysiological substrate of psychoemotional symptomatology in patients with lupus, which may assist its clinical diagnosis and treatment.

Tibetan medicine Duoxuekang ameliorates hypobaric hypoxia-induced brain injury in mice by restoration of cerebrovascular function

ETHNOPHARMACOLOGICAL RELEVANCE

Duoxuekang (DXK, ཁྲག་འཕེལ་བདེ་བྱེད།) is a clinical experience prescription of CuoRu-Cailang, a famous Tibetan medicine master, which has effective advantages in the treatment of hypobaric hypoxia (HH)-induced brain injury. However, its underlying mechanisms remain unclear.

AIM OF THE STUDY

The present study was designed to investigate the effects of DXK on cerebrovascular function of HH-induced brain injury in mice.

MATERIALS AND METHODS

DSC-MR imaging was used to evaluate the effect of DXK on the brain blood perfusion of patients with hypoxic brain injury. HPLC analysis was used to detect the content of salidroside, gallic acid, tyrosol, corilagin, ellagic acid, isorhamnetin, quercetin and gingerol in DXK. The model of HH-induced brain injury in mice was established by an animal hypobaric and hypoxic chamber. The BABL/c mice were randomly divided into six groups: control group, model group, Hongjingtian oral liquid group (HOL, 3.3 ml/kg) and DXK groups (0.9, 1.8 and 3.6 g/kg). All mice (except the control group) were intragastrically administrated for a continuous 7 days and put into the animal hypobaric and hypoxic chamber after the last intragastric administration. Hematoxylin-eosin staining was employed to evaluate the pathological changes of brain tissue. Masson and Weigert stainings were used to detect the content of collagen fibers and elastic fibers of brain, respectively. Routine blood test and biochemical kits were used to analyze hematological parameters and oxidative stress indices. Immunofluorescence staining was applied to detect the protein levels of VEGF, CD31/vWF and α-SMA.

RESULTS

The results of DSC-MR imaging confirmed that DXK can increased CBV in the left temporal lobe while decreased MTT in the right frontal lobe, right temporal lobe and right occipital lobe of the brain. DXK contains salidroside, gallic acid, tyrosol, corilagin, ellagic acid, isorhamnetin, quercetin and gingerol. Compared with the model group, DXK can ameliorate the atrophy and deformation, and increase the number of pyramidal neurons in hippocampal CA3 area and cortical neurocytes. Masson and Weigert stainings results revealed that DXK can significantly increase the content of collagen fibers and elastic fibers in brain. Routine blood test results demonstrated that DXK can dramatically decrease the levels of WBC, MCH and MCHC, while increase RBC, HGB, HCT, MCV and PLT in the blood samples. Biochemical results revealed that DXK can markedly increase SOD, CAT and GSH activities, while decrease MDA activity. Immunofluorescence revealed that DXK can notably increase the protein levels of VEGF, CD31/vWF and α-SMA.

CONCLUSIONS

In conclusion, this study proved that DXK can ameliorate HH-induced brain injury by improving brain blood perfusion, increasing the number of collagen and elastic fibers and inhibiting oxidative stress injury. The underlying mechanisms may be involved in maintaining the integrity of cerebrovascular endothelial cells and vascular function. However, further in vivo and in vitro investigations are still needed to elucidate the mechanisms of DXK on regulating cerebral blood vessels.

Evidence of Age-Related Hemodynamic and Functional Connectivity Impairment: A Resting State fMRI Study

Purpose: To assess age-related changes in intrinsic functional brain connectivity and hemodynamics during adulthood in the context of the retrogenesis hypothesis, which states that the rate of age-related changes is higher in late-myelinating (prefrontal, lateral-posterior temporal) cerebrocortical areas as compared to early myelinating (parietal, occipital) regions. In addition, to examine the dependence of age-related changes upon concurrent subclinical depression symptoms which are common even in healthy aging. Methods: Sixty-four healthy adults (28 men) aged 23-79 years (mean 45.0, SD = 18.8 years) were examined. Resting-state functional MRI (rs-fMRI) time series were used to compute voxel-wise intrinsic connectivity contrast (ICC) maps reflecting the strength of functional connectivity between each voxel and the rest of the brain. We further used Time Shift Analysis (TSA) to estimate voxel-wise hemodynamic lead or lag for each of 22 ROIs from the automated anatomical atlas (AAL). Results: Adjusted for depression symptoms, gender and education level, reduced ICC with age was found primarily in frontal, temporal regions, and putamen, whereas the opposite trend was noted in inferior occipital cortices (p < 0.002). With the same covariates, increased hemodynamic lead with advancing age was found in superior frontal cortex and thalamus, with the opposite trend in inferior occipital cortex (p < 0.002). There was also evidence of reduced coupling between voxel-wise intrinsic connectivity and hemodynamics in the inferior parietal cortex. Conclusion: Age-related intrinsic connectivity reductions and hemodynamic changes were demonstrated in several regions-most of them part of DMN and salience networks-while impaired neurovascular coupling was, also, found in parietal regions. Age-related reductions in intrinsic connectivity were greater in anterior as compared to posterior cortices, in line with implications derived from the retrogenesis hypothesis. These effects were affected by self-reported depression symptoms, which also increased with age.

15 Years MR-encephalography

Objective

This review article gives an account of the development of the MR-encephalography (MREG) method, which started as a mere ‘Gedankenexperiment’ in 2005 and gradually developed into a method for ultrafast measurement of physiological activities in the brain. After going through different approaches covering k-space with radial, rosette, and concentric shell trajectories we have settled on a stack-of-spiral trajectory, which allows full brain coverage with (nominal) 3 mm isotropic resolution in 100 ms. The very high acceleration factor is facilitated by the near-isotropic k-space coverage, which allows high acceleration in all three spatial dimensions.

Methods

The methodological section covers the basic sequence design as well as recent advances in image reconstruction including the targeted reconstruction, which allows real-time feedback applications, and—most recently—the time-domain principal component reconstruction (tPCR), which applies a principal component analysis of the acquired time domain data as a sparsifying transformation to improve reconstruction speed as well as quality.

Applications

Although the BOLD-response is rather slow, the high speed acquisition of MREG allows separation of BOLD-effects from cardiac and breathing related pulsatility. The increased sensitivity enables direct detection of the dynamic variability of resting state networks as well as localization of single interictal events in epilepsy patients. A separate and highly intriguing application is aimed at the investigation of the glymphatic system by assessment of the spatiotemporal patterns of cardiac and breathing related pulsatility.

Discussion

MREG has been developed to push the speed limits of fMRI. Compared to multiband-EPI this allows considerably faster acquisition at the cost of reduced image quality and spatial resolution.

Respiratory-related brain pulsations are increased in epilepsy-a two-centre functional MRI study

Resting-state functional MRI has shown potential for detecting changes in cerebral blood oxygen level-dependent signal in patients with epilepsy, even in the absence of epileptiform activity. Furthermore, it has been suggested that coefficient of variation mapping of fast functional MRI signal may provide a powerful tool for the identification of intrinsic brain pulsations in neurological diseases such as dementia, stroke and epilepsy. In this study, we used fast functional MRI sequence (magnetic resonance encephalography) to acquire ten whole-brain images per second. We used the functional MRI data to compare physiological brain pulsations between healthy controls (n = 102) and patients with epilepsy (n = 33) and furthermore to drug-naive seizure patients (n = 9). Analyses were performed by calculating coefficient of variation and spectral power in full band and filtered sub-bands. Brain pulsations in the respiratory-related frequency sub-band (0.11-0.51 Hz) were significantly (P < 0.05) increased in patients with epilepsy, with an increase in both signal variance and power. At the individual level, over 80% of medicated and drug-naive seizure patients exhibited areas of abnormal brain signal power that correlated well with the known clinical diagnosis, while none of the controls showed signs of abnormality with the same threshold. The differences were most apparent in the basal brain structures, respiratory centres of brain stem, midbrain and temporal lobes. Notably, full-band, very low frequency (0.01-0.1 Hz) and cardiovascular (0.8-1.76 Hz) brain pulses showed no differences between groups. This study extends and confirms our previous results of abnormal fast functional MRI signal variance in epilepsy patients. Only respiratory-related brain pulsations were clearly increased with no changes in either physiological cardiorespiratory rates or head motion between the subjects. The regional alterations in brain pulsations suggest that mechanisms driving the cerebrospinal fluid homeostasis may be altered in epilepsy. Magnetic resonance encephalography has both increased sensitivity and high specificity for detecting the increased brain pulsations, particularly in times when other tools for locating epileptogenic areas remain inconclusive.

Is the Blood Oxygenation Level-Dependent fMRI Response to Motor Tasks Altered in Children After Neonatal Stroke?

Functional MRI is increasingly being used in the assessment of brain activation and connectivity following stroke. Many of these studies rely on the Blood Oxygenation Level Dependent (BOLD) contrast. However, the stability, as well as the accuracy of the BOLD response to motor task in the ipsilesional hemisphere, remains ambiguous. In this work, the BOLD signal acquired from both healthy and affected hemispheres was analyzed in 7-year-old children who sustained a Neonatal Arterial Ischemic Stroke (NAIS). Accordingly, a repetitive motor task of the contralesional and the ipsilesional hands was performed by 33 patients with unilateral lesions. These patients were divided into two groups: those without cerebral palsy (NAIS), and those with cerebral palsy (CP). The BOLD signal time course was obtained from distinctly defined regions of interest (ROIs) extracted from the functional activation maps of 30 healthy controls with similar age and demographic characteristics as the patients. An ROI covering both the primary motor cortex (M1) and the primary somatosensory cortex (S1) was also tested. Compared with controls, NAIS patients without CP had similar BOLD amplitude variation for both the contralesional and the ipsilesional hand movements. However, in the case of NAIS patients with CP, a significant difference in the averaged BOLD amplitude was found between the healthy and affected hemisphere. In both cases, no progressive attenuation of the BOLD signal amplitude was observed throughout the task epochs. Besides, results also showed a correlation between the BOLD signal percentage variation of the lesioned hemisphere and the dexterity level. These findings suggest that for patients who sustained a NAIS with no extensive permanent motor impairment, BOLD signal-based data analysis can be a valuable tool for the evaluation of functional brain networks.

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.

Symptomatic psychosis risk and physiological fluctuation in functional MRI data

BACKGROUND

Physiological brain pulsations have been shown to play a critical role in maintaining interstitial homeostasis in the glymphatic brain clearance mechanism. We investigated whether psychotic symptomatology is related to the physiological variation of the human brain using fMRI.

METHODS

The participants (N = 277) were from the Northern Finland Birth Cohort 1986. Psychotic symptoms were evaluated with the Positive Symptoms Scale of the Structured Interview for Prodromal Syndromes (SIPS). We used the coefficient of variation of BOLD signal (CV_BOLD) as a proxy for physiological brain pulsatility. The CV_BOLD-analyses were controlled for motion, age, sex, and educational level. The results were also compared with fMRI and voxel-based morphometry (VBM) meta-analyses of schizophrenia patients (data from the Brainmap database).

RESULTS

At the global level, participants with psychotic-like symptoms had higher CV_BOLD in cerebrospinal fluid (CSF) and white matter (WM), when compared to participants with no psychotic symptoms. Voxel-wise analyses revealed that CV_BOLD was increased, especially in periventricular white matter, basal ganglia, cerebellum and parts of the cortical structures. Those brain regions, which included alterations of physiological fluctuation in symptomatic psychosis risk, overlapped <6% with the regions that were found to be affected in the meta-analyses of previous fMRI and VBM studies in schizophrenia patients. Motion did not vary as a function of SIPS.

CONCLUSIONS

Psychotic-like symptoms were associated with elevated CV_BOLD in a variety of brain regions. The CV_BOLD findings may produce new information about cerebral physiological fluctuations that have been out of reach in previous fMRI and VBM studies.

The relationship of genetic susceptibilities for psychosis with physiological fluctuation in functional MRI data

Previously, schizophrenia is found to be related to the variability of the functional magnetic resonance imaging (fMRI) signal in the white matter. However, evidence about the relationship between genetic vulnerabilities and physiological fluctuation in the brain is lacking. We investigated whether familial risk for psychosis (FR) and polygenic risk score for schizophrenia (PRS) are linked with physiological fluctuation in fMRI data. We used data from the Oulu Brain and Mind study (n = 140-149, aged 20-24 years) that is a substudy of the Northern Finland Birth Cohort 1986. The participants underwent a resting-state fMRI scan. Coefficient of variation (CV) of blood oxygen level dependent (BOLD) signal (CVBOLD) was used as a proxy of physiological fluctuation in the brain. Familial risk was defined to be present if at least one parent had been diagnosed with psychosis previously. PRS was computed based on the results of the prior GWAS by the Schizophrenia Working Group. FR or PRS were not associated with CVBOLD in cerebrospinal fluid, white matter, or grey matter. The findings did not provide evidence for the previous suggestions that genetic vulnerabilities for schizophrenia become apparent in alterations of the variation of the BOLD signal in the brain.

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.

Low Frequency Systemic Hemodynamic "Noise" in Resting State BOLD fMRI: Characteristics, Causes, Implications, Mitigation Strategies, and Applications

Advances in functional magnetic resonance imaging (fMRI) acquisition have improved signal to noise to the point where the physiology of the subject is the dominant noise source in resting state fMRI data (rsfMRI). Among these systemic, non-neuronal physiological signals, respiration and to some degree cardiac fluctuations can be removed through modeling, or in the case of newer, faster acquisitions such as simultaneous multislice acquisition, simple spectral filtering. However, significant low frequency physiological oscillation (∼0.01-0.15 Hz) remains in the signal. This is problematic, as it is the precise frequency band occupied by the neuronally modulated hemodynamic responses used to study brain connectivity, precluding its removal by spectral filtering. The source of this signal, and its method of production and propagation in the body, have not been conclusively determined. Here, we summarize the defining characteristics of the systemic low frequency noise signal, and review some current theories about the signal source and the evidence supporting them. The strength and distribution of the systemic LFO signal make characterizing and removing it essential for accurate quantification, especially for resting state connectivity, when no stimulation can be compared with the signal. Widespread correlated non-neuronal signals obscure and distort the more localized patterns of neuronal correlations between interacting brain regions; they may even cause apparent connectivity between regions with no neuronal interaction. Here, we discuss a simple method we have developed to parse the global, moving, blood-borne signal from the stationary, neuronal connectivity signals, substantially reducing the negative correlations that result from global signal regression. Finally, we will discuss some of the uses to which the moving systemic low frequency oscillation can be put if we consider it a "signal" carrying information, rather than simply "noise" complicating the interpretation of resting state connectivity. Properly utilizing this signal may offer insights into subtle hemodynamic alterations that can be used as early indicators of circulatory dysfunction in a number of neuropsychiatric conditions, such as prodromal stroke, moyamoya, and Alzheimer's disease.

The impact of ischemic stroke on connectivity gradients

The functional organization of the brain can be represented as a low-dimensional space that reflects its macroscale hierarchy. The dimensions of this space, described as connectivity gradients, capture the similarity of areas' connections along a continuous space. Studying how pathological perturbations with known effects on functional connectivity affect these connectivity gradients provides support for their biological relevance. Previous work has shown that localized lesions cause widespread functional connectivity alterations in structurally intact areas, affecting a network of interconnected regions. By using acute stroke as a model of the effects of focal lesions on the connectome, we apply the connectivity gradient framework to depict how functional reorganization occurs throughout the brain, unrestricted by traditional definitions of functional network boundaries. We define a three-dimensional connectivity space template based on functional connectivity data from healthy controls. By projecting lesion locations into this space, we demonstrate that ischemic strokes result in dimension-specific alterations in functional connectivity over the first week after symptom onset. Specifically, changes in functional connectivity were captured along connectivity Gradients 1 and 3. The degree of functional connectivity change was associated with the distance from the lesion along these connectivity gradients (a measure of functional similarity) regardless of the anatomical distance from the lesion. Together, these results provide support for the biological validity of connectivity gradients and suggest a novel framework to characterize connectivity alterations after stroke.

Spreading Depolarization Waves in Neurological Diseases: A Short Review about its Pathophysiology and Clinical Relevance

Lesion growth following acutely injured brain tissue after stroke, subarachnoid hemorrhage and traumatic brain injury is an important issue and a new target area for promising therapeutic interventions. Spreading depolarization or peri-lesion depolarization waves were demonstrated as one of the significant contributors of continued lesion growth. In this short review, we discuss the pathophysiology for SD forming events and try to list findings detected in neurological disorders like migraine, stroke, subarachnoid hemorrhage and traumatic brain injury in both human as well as experimental studies. Pharmacological and non-pharmacological treatment strategies are highlighted and future directions and research limitations are discussed.

One-step analysis of brain perfusion and function for acute stroke patients after reperfusion: A resting-state fMRI study

BACKGROUND

Resting-state functional magnetic resonance imaging (rs-fMRI) can noninvasively estimate the perfusion and function of the brain.

PURPOSE

To investigate the perfusion and functional status using rs-fMRI in acute ischemic stroke (AIS) patients after reperfusion therapy.

STUDY TYPE

Prospective.

SUBJECTS

Twenty-five AIS patients who underwent dynamic susceptibility contrast (DSC) upon hospital admission and both rs-fMRI and DSC scans at 24 hours after reperfusion therapy.

FIELD STRENGTH/SEQUENCE

3T; DSC, rs-fMRI.

ASSESSMENT

The time delay of the blood oxygenation level-dependent (BOLD) signal was calculated using time-shift-analysis (TSA) and compared with the time to peak (TTP) derived from the DSC. For patients who exhibited partial or complete reperfusion in the supratentorial hemisphere, we quantified the function of different regions (healthy tissue, reperfused tissue, not reperfused tissue) by using three rs-fMRI measurements (functional connectivity, the amplitude of low-frequency fluctuation [ALFF] and regional homogeneity [ReHo]). Correlations between the functional measurements and modified Rankin Scale (mRS) scores were calculated.

STATISTICAL TESTS

Dice coefficient (DC) analysis, two-sample t-tests, Pearson correlation coefficient.

RESULTS

Twelve patients who exhibited complete reperfusion on their TTP maps showed no time-delayed areas on the TSA maps. For the remaining 13 patients with partial reperfusion (5/13) or no reperfusion (8/13) on the TTP maps, the TSA detected comparable time-delayed areas. Eleven out of 13 patients showed moderate to good overlap (mean DC, 0.58 ± 0.1) between the TTP and TSA results. Fourteen patients were chosen for functional analyses and most patients (12/14) showed abnormal functional connectivity in the reperfused regions. The reperfused and not reperfused tissues had lower mean ReHo values than those of the healthy tissue (both P < 0.001). The mRS scores showed negative correlation with mean ReHo values of reperfused region (R = -0.523, P = 0.027). DATA CONCLUSION: rs-fMRI might be a useful way to estimate both the perfusion and functional status for AIS patients after reperfusion therapy.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:221-229.

Non-invasive evaluation of cerebral perfusion in patients with transient ischemic attack: an fMRI study

Detection of hypoperfused tissue due to the ischemia is considered to be important in understanding the cerebral perfusion status and may be helpful in guiding therapeutic decisions for patients with transient ischemic attack (TIA). We hypothesized that the combination of two non-invasive fMRI techniques: resting-state BOLD-fMRI time-shift analysis (TSA) approach and 3D ASL, could detect the cerebral hemodynamic status in TIA patients noninvasively. From April 2015 to June 2016, 51 TIA patients were recruited in this study. We calculated the time delay between the resting-state BOLD signal at each voxel and the whole-brain signal using TSA approach and compared the results to CBF map derived from ASL. Out of the 51 patients, 24 patients with normal arrival time and CBF were in Stage 0; 14 patients who showed delayed arrival time and normal CBF which indicated elevated CBV were in Stage I; the other 13 patients who had both delayed arrival time and decreased CBF were in Stage II, the group average spatial overlap, i.e., Dice coefficient, of the two measurements was 0.55. Four patients in Stage 0 (17.4%), three patients in Stage I (23.1%) and five patients in Stage II (45.5%) suffered ischemic stroke or TIA symptoms in 1 year after MRI scan. The patients in Stage II was at highest risk of subsequent events when compared to other two stages. The combination of resting-state BOLD-fMRI and ASL hold the potential to noninvasively identify the hemodynamic status in TIA patients and help predict the risk of subsequent events.

Detecting Perfusion Pattern Based on the Background Low-Frequency Fluctuation in Resting-State Functional Magnetic Resonance Imaging Data and Its Influence on Resting-State Networks: An Iterative Postprocessing Approach

Resting-state functional magnetic resonance imaging (RS-fMRI) is based on the assumption that the vascular response and the blood oxygenation level-dependent response are homogenous across the entire brain. However, this a priori hypothesis is not consistent with the well-known variability of cerebral vascular territories. To explore whether the RS networks are influenced by varied vascular speed in different vascular territories, we assessed the time-shift maps that give an estimate of the local timing of the vascular response and checked whether local differences in this timing have an impact on the estimates of RS networks. Two hundred seventeen elderly (≥60 years), healthy participants (73.74 ± 4.41 years, 143 females, 203 right handed) underwent one MRI examination, including an RS-fMRI session. After preprocessing, statistical analyses included time-shift analyses and RS-fMRI analyses using as regressor the delay maps obtained from the time-shift analyses. The functional connectivity map of default mode network (DMN) of each participant was then calculated by using the seed-to-voxel analysis in the REST toolbox. Faster cerebrovascular responses were notably present in the primary motor and somatosensory and peri-insular cortex, while slower responses were present in various regions, including notably the posterior cingulate cortex (PCC). Moreover, significant changes notably in the DMN, including medial prefrontal cortex (t = 11.95), PCC (t = 11.52), right middle temporal lobe (t = 10.72), and right angular gyrus (t = 10.88), were observed also taking into account the cerebrovascular delayed maps. As the most prominent example of the RS networks, DMN activation patterns change as a function of the cerebrovascular delay. These data suggest that a group correction for vascular maps in RS-fMRI measurements is essential to correctly depict functional differences and exclude potential confounding effects, notably in the elderly with increasing prevalence of vascular comorbidity.

Altered physiological brain variation in drug-resistant epilepsy

INTRODUCTION

Functional magnetic resonance imaging (fMRI) combined with simultaneous electroencephalography (EEG-fMRI) has become a major tool in mapping epilepsy sources. In the absence of detectable epileptiform activity, the resting state fMRI may still detect changes in the blood oxygen level-dependent signal, suggesting intrinsic alterations in the underlying brain physiology.

METHODS

In this study, we used coefficient of variation (CV) of critically sampled 10 Hz ultra-fast fMRI (magnetoencephalography, MREG) signal to compare physiological variance between healthy controls (n = 10) and patients (n = 10) with drug-resistant epilepsy (DRE).

RESULTS

We showed highly significant voxel-level (p < 0.01, TFCE-corrected) increase in the physiological variance in DRE patients. At individual level, the elevations range over three standard deviations (σ) above the control mean (μ) CV_MREG values solely in DRE patients, enabling patient-specific mapping of elevated physiological variance. The most apparent differences in group-level analysis are found on white matter, brainstem, and cerebellum. Respiratory (0.12-0.4 Hz) and very-low-frequency (VLF = 0.009-0.1 Hz) signal variances were most affected.

CONCLUSIONS

The CV_MREG increase was not explained by head motion or physiological cardiorespiratory activity, that is, it seems to be linked to intrinsic physiological pulsations. We suggest that intrinsic brain pulsations play a role in DRE and that critically sampled fMRI may provide a powerful tool for their identification.

Cerebral Hypoperfusion and Other Shared Brain Pathologies in Ischemic Stroke and Alzheimer's Disease

Newly emerged evidence reveals that ischemic stroke and Alzheimer's disease (AD) share pathophysiological changes in brain tissue including hypoperfusion, oxidative stress, immune exhaustion, and inflammation. A mechanistic link between hypoperfusion and amyloid β accumulation can lead to cell damage as well as to motor and cognitive deficits. This review will discuss decreased cerebral perfusion and other related pathophysiological changes common to both ischemic stroke and AD, such as vascular damages, cerebral blood flow alteration, abnormal expression of amyloid β and tau proteins, as well as behavioral and cognitive deficits. Furthermore, this review highlights current treatment options and potential therapeutic targets that warrant further investigation.

Brain perfusion alterations in tick-borne encephalitis-preliminary report

BACKGROUND

Magnetic resonance imaging (MRI) changes in tick-borne encephalitis (TBE) are non-specific and the pathophysiological mechanisms leading to their formation remain unclear. This study investigated brain perfusion in TBE patients using dynamic susceptibility-weighted contrast-enhanced magnetic resonance perfusion imaging (DSC-MRI perfusion).

METHODS

MRI scans were performed for 12 patients in the acute phase, 3-5days after the diagnosis of TBE. Conventional MRI and DSC-MRI perfusion studies were performed. Cerebral blood flow (CBF), cerebral blood volume (CBV), mean transit time (MTT), and time to peak (TTP) parametric maps were created. The bilateral frontal, parietal, and temporal subcortical regions and thalamus were selected as regions of interest. Perfusion parameters of TBE patients were compared to those of a control group.

RESULTS

There was a slight increase in CBF and CBV, with significant prolongation of TTP in subcortical areas in the study subjects, while MTT values were comparable to those of the control group. A significant increase in thalamic CBF (p<0.001) and increased CBV (p<0.05) were observed. Increased TTP and a slight reduction in MTT were also observed within this area.

CONCLUSIONS

The DSC-MRI perfusion study showed that TBE patients had brain perfusion disturbances, expressed mainly in the thalami. These results suggest that DSC-MRI perfusion may provide important information regarding the areas affected in TBE patients.