High-speed real-time resting-state FMRI using multi-slab echo-volumar imaging

Real-time fMRI using multi-band echo-volumar imaging with millimeter spatial resolution and sub-second temporal resolution at 3 tesla

Purpose

In this study we develop undersampled echo-volumar imaging (EVI) using multi-band/simultaneous multi-slab encoding in conjunction with multi-shot slab-segmentation to accelerate 3D encoding and to reduce the duration of EVI encoding within slabs. This approach combines the sampling efficiency of single-shot 3D encoding with the sensitivity advantage of multi-echo acquisition. We describe the pulse sequence development and characterize the spatial-temporal resolution limits and BOLD sensitivity of this approach for high-speed task-based and resting-state fMRI at 3 T. We study the feasibility of further acceleration using compressed sensing (CS) and assess compatibility with NORDIC denoising.

Methods

Multi-band echo volumar imaging (MB-EVI) combines multi-band encoding of up to 6 slabs with CAIPI shifting, accelerated EVI encoding within slabs using up to 4-fold GRAPPA accelerations, 2-shot kz-segmentation and partial Fourier acquisitions along the two phase-encoding dimensions. Task-based and resting-state fMRI at 3 Tesla was performed across a range of voxel sizes (between 1 and 3 mm isotropic), repetition times (118-650 ms), and number of slabs (up to 12). MB-EVI was compared with multi-slab EVI (MS-EVI) and multi-band-EPI (MB-EPI).

Results

Image quality and temporal SNR of MB-EVI was comparable to MS-EVI when using 2-3 mm spatial resolution. High sensitivity for mapping task-based activation and resting-state connectivity at short TR was measured. Online deconvolution of T2* signal decay markedly reduced spatial blurring and improved image contrast. The high temporal resolution of MB-EVI enabled sensitive mapping of high-frequency resting-state connectivity above 0.3 Hz with 3 mm isotropic voxel size (TR: 163 ms). Detection of task-based activation with 1 mm isotropic voxel size was feasible in scan times as short as 1 min 13 s. Compressed sensing with up to 2.4-fold retrospective undersampling showed negligible loss in image quality and moderate region-specific losses in BOLD sensitivity. NORDIC denoising significantly enhanced fMRI sensitivity without introducing image blurring.

Conclusion

Combining MS-EVI with multi-band encoding enables high overall acceleration factors and provides flexibility for maximizing spatial-temporal resolution and volume coverage. The high BOLD sensitivity of this hybrid MB-EVI approach and its compatibility with online image reconstruction enables high spatial-temporal resolution real-time task-based and resting state fMRI.

Characterization of pulsations in the brain and cerebrospinal fluid using ultra-high field magnetic resonance imaging

The development of innovative non-invasive neuroimaging methods and biomarkers is critical for studying brain disease. Imaging of cerebrospinal fluid (CSF) pulsatility may inform the brain fluid dynamics involved in clearance of cerebral metabolic waste. In this work, we developed a methodology to characterize the frequency and spatial localization of whole brain CSF pulsations in humans. Using 7 Tesla (T) human magnetic resonance imaging (MRI) and ultrafast echo-planar imaging (EPI), in-vivo images were obtained to capture pulsations of the CSF signal. Physiological data were simultaneously collected and compared with the 7 T MR data. The primary components of signal pulsations were identified using spectral analysis, with the most evident frequency bands identified around 0.3, 1.2, and 2.4 Hz. These pulsations were mapped spatially and temporally onto the MR image domain and temporally onto the physiological measures of electrocardiogram and respiration. We identified peaks in CSF pulsations that were distinct from peaks in grey matter and white matter regions. This methodology may provide novel in vivo biomarkers of disrupted brain fluid dynamics.

Effect of sleep deprivation and NREM sleep stage on physiological brain pulsations

Introduction

Sleep increases brain fluid transport and the power of pulsations driving the fluids. We investigated how sleep deprivation or electrophysiologically different stages of non-rapid-eye-movement (NREM) sleep affect the human brain pulsations.

Methods

Fast functional magnetic resonance imaging (fMRI) was performed in healthy subjects (n = 23) with synchronous electroencephalography (EEG), that was used to verify arousal states (awake, N1 and N2 sleep). Cardiorespiratory rates were verified with physiological monitoring. Spectral power analysis assessed the strength, and spectral entropy assessed the stability of the pulsations.

Results

In N1 sleep, the power of vasomotor (VLF < 0.1 Hz), but not cardiorespiratory pulsations, intensified after sleep deprived vs. non-sleep deprived subjects. The power of all three pulsations increased as a function of arousal state (N2 > N1 > awake) encompassing brain tissue in both sleep stages, but extra-axial CSF spaces only in N2 sleep. Spectral entropy of full band and respiratory pulsations decreased most in N2 sleep stage, while cardiac spectral entropy increased in ventricles.

Discussion

In summary, the sleep deprivation and sleep depth, both increase the power and harmonize the spectral content of human brain pulsations.

Measuring brain beats: Cardiac-aligned fast functional magnetic resonance imaging signals

Blood and cerebrospinal fluid (CSF) pulse and flow throughout the brain, driven by the cardiac cycle. These fluid dynamics, which are essential to healthy brain function, are characterized by several noninvasive magnetic resonance imaging (MRI) methods. Recent developments in fast MRI, specifically simultaneous multislice acquisition methods, provide a new opportunity to rapidly and broadly assess cardiac-driven flow, including CSF spaces, surface vessels and parenchymal vessels. We use these techniques to assess blood and CSF flow dynamics in brief (3.5 min) scans on a conventional 3 T MRI scanner in five subjects. Cardiac pulses are measured with a photoplethysmography (PPG) on the index finger, along with functional MRI (fMRI) signals in the brain. We, retrospectively, align the fMRI signals to the heartbeat. Highly reliable cardiac-gated fMRI temporal signals are observed in CSF and blood on the timescale of one heartbeat (test-retest reliability within subjects R²  > 50%). In blood vessels, a local minimum is observed following systole. In CSF spaces, the ventricles and subarachnoid spaces have a local maximum following systole instead. Slower resting-state scans with slice timing, retrospectively, aligned to the cardiac pulse, reveal similar cardiac-gated responses. The cardiac-gated measurements estimate the amplitude and phase of fMRI pulsations in the CSF relative to those in the arteries, an estimate of the local intracranial impedance. Cardiac aligned fMRI signals can provide new insights about fluid dynamics or diagnostics for diseases where these dynamics are important.

The glymphatic system: Current understanding and modeling

We review theoretical and numerical models of the glymphatic system, which circulates cerebrospinal fluid and interstitial fluid around the brain, facilitating solute transport. Models enable hypothesis development and predictions of transport, with clinical applications including drug delivery, stroke, cardiac arrest, and neurodegenerative disorders like Alzheimer's disease. We sort existing models into broad categories by anatomical function: Perivascular flow, transport in brain parenchyma, interfaces to perivascular spaces, efflux routes, and links to neuronal activity. Needs and opportunities for future work are highlighted wherever possible; new models, expanded models, and novel experiments to inform models could all have tremendous value for advancing the field.

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).

Cardiovascular Pulsatility Increases in Visual Cortex Before Blood Oxygen Level Dependent Response During Stimulus

The physiological pulsations that drive tissue fluid homeostasis are not well characterized during brain activation. Therefore, we used fast magnetic resonance encephalography (MREG) fMRI to measure full band (0–5 Hz) blood oxygen level-dependent (BOLD_FB) signals during a dynamic visual task in 23 subjects. This revealed brain activity in the very low frequency (BOLD_VLF) as well as in cardiac and respiratory bands. The cardiovascular hemodynamic envelope (CHe) signal correlated significantly with the visual BOLD_VLF response, considered as an independent signal source in the V1-V2 visual cortices. The CHe preceded the canonical BOLD_VLF response by an average of 1.3 (± 2.2) s. Physiologically, the observed CHe signal could mark increased regional cardiovascular pulsatility following vasodilation.

Case Series of Applications of Resting State Functional MRI in Brain Tumor Surgery: A Novel Technique

Background The extent of resection for brain tumors is a critical factor in determining the oncologic outcome for a patient. However, a balance between preservation of neurological function and maximal resection is essential for true benefit.

Functional magnetic resonance imaging (fMRI) is one of the approaches that augments the neurosurgeon's ability to attain maximal safe resection by providing preoperative mapping. It may not be possible to perform awake craniotomy with intraoperative localization by direct cortical stimulation in all patients, such as children and those with neurocognitive impairment. Task-based fMRI may have limited value in these cases due to low patient cooperability.

Methods In this article we present in a case-based format, the various clinical scenarios where resting state fMRI (rs-fMRI) can be helpful in guiding neurosurgical resection. rs-fMRI of the patients has been acquired on Philips 1.5 T system. Seed voxel method has been used for processing and analysis.

Conclusion rs-fMRI does not require active patient cooperation to generate useful information and thus can be a promising tool in patients unable to cooperate for task-based studies.

Slab boundary artifact correction in multislab imaging using convolutional-neural-network-enabled inversion for slab profile encoding

PURPOSE

This study aims to propose a novel algorithm for slab boundary artifact correction in both single-band multislab imaging and simultaneous multislab (SMSlab) imaging.

THEORY AND METHODS

In image domain, the formation of slab boundary artifacts can be regarded as modulating the artifact-free images using the slab profiles and introducing aliasing along the slice direction. Slab boundary artifact correction is the inverse problem of this process. An iterative algorithm based on convolutional neural networks (CNNs) is proposed to solve the problem, termed CNN-enabled inversion for slab profile encoding (CPEN). Diffusion-weighted SMSlab images and reference images without slab boundary artifacts were acquired in 7 healthy subjects for training. Images of 5 healthy subjects were acquired for testing, including single-band multislab and SMSlab images with 1.3-mm or 1-mm isotropic resolution. CNN-enabled inversion for slab profile encoding was compared with a previously reported method (i.e., nonlinear inversion for slab profile encoding [NPEN]).

RESULTS

CNN-enabled inversion for slab profile encoding reduces the slab boundary artifacts in both single-band multislab and SMSlab images. It also suppresses the slab boundary artifacts in the diffusion metric maps. Compared with NPEN, CPEN shows fewer residual artifacts in different acquisition protocols and more significant improvements in quantitative assessment, and it also accelerates the computation by more than 35 times.

CONCLUSION

CNN-enabled inversion for slab profile encoding can reduce the slab boundary artifacts in multislab acquisitions. It shows better slab boundary artifact correction capacity, higher robustness, and computation efficiency when compared with NPEN. It has the potential to improve the accuracy of multislab acquisitions in high-resolution DWI and functional MRI.

Spectral analysis of physiological brain pulsations affecting the BOLD signal

Physiological pulsations have been shown to affect the global blood oxygen level dependent (BOLD) signal in human brain. While these pulsations have previously been regarded as noise, recent studies show their potential as biomarkers of brain pathology. We used the extended 5 Hz spectral range of magnetic resonance encephalography (MREG) data to investigate spatial and frequency distributions of physiological BOLD signal sources. Amplitude spectra of the global image signals revealed cardiorespiratory envelope modulation (CREM) peaks, in addition to the previously known very low frequency (VLF) and cardiorespiratory pulsations. We then proceeded to extend the amplitude of low frequency fluctuations (ALFF) method to each of these pulsations. The respiratory pulsations were spatially dominating over most brain structures. The VLF pulsations overcame the respiratory pulsations in frontal and parietal gray matter, whereas cardiac and CREM pulsations had this effect in central cerebrospinal fluid (CSF) spaces and major blood vessels. A quasi‐periodic pattern (QPP) analysis showed that the CREM pulsations propagated as waves, with a spatiotemporal pattern differing from that of respiratory pulsations, indicating them to be distinct intracranial physiological phenomenon. In conclusion, the respiration has a dominant effect on the global BOLD signal and directly modulates cardiovascular brain pulsations.

Cardiovascular brain impulses in Alzheimer's disease

Accumulation of amyloid-β is a key neuropathological feature in brain of Alzheimer’s disease patients. Alterations in cerebral haemodynamics, such as arterial impulse propagation driving the (peri)vascular CSF flux, predict future Alzheimer’s disease progression. We now present a non-invasive method to quantify the three-dimensional propagation of cardiovascular impulses in human brain using ultrafast 10 Hz magnetic resonance encephalography. This technique revealed spatio-temporal abnormalities in impulse propagation in Alzheimer’s disease. The arrival latency and propagation speed both differed in patients with Alzheimer’s disease. Our mapping of arterial territories revealed Alzheimer’s disease-specific modifications, including reversed impulse propagation around the hippocampi and in parietal cortical areas. The findings imply that pervasive abnormality in (peri)vascular CSF impulse propagation compromises vascular impulse propagation and subsequently glymphatic brain clearance of amyloid-β in Alzheimer’s disease.

Trading off spatio-temporal properties in 3D high-speed fMRI using interleaved stack-of-spirals trajectories

PURPOSE

Highly undersampled acquisitions have been proposed to push the limits of temporal resolution in functional MRI. This contribution is aimed at identifying parameter sets that let the user trade-off between ultra-high temporal resolution and spatial signal quality by varying the sampling densities. The proposed method maintains the synergies of a temporal resolution that enables direct filtering of physiological artifacts for highest statistical power, and 3D read-outs with optimal use of encoding capabilities of multi-coil arrays for efficient sampling and high signal-to-noise ratio (SNR).

METHODS

One- to four-shot interleaved spherical stack-of-spiral trajectories with repetition times from 96 to 352 ms at a nominal resolution of 3 mm using different sampling densities were compared for image quality and temporal SNR (tSNR). The one- and three-shot trajectories were employed in a resting state study for functional characterization.

RESULTS

Compared to a previously described single-shot trajectory, denser sampled trajectories of the same type are shown to be less prone to blurring and off-resonance vulnerability that appear in addition to the variable density artifacts of the point spread function. While the multi-shot trajectories lead to a decrease in tSNR efficiency, the high SNR due to the 3D read-out, combined with notable increases in image quality, leads to superior overall results of the three-shot interleaved stack of spirals. A resting state analysis of 15 subjects shows significantly improved functional sensitivity in areas of high off-resonance gradients.

CONCLUSION

Mild variable-density sampling leads to excellent tSNR behavior and no increased off-resonance vulnerability, and is suggested unless maximum temporal resolution is sought.

Self-Positivity or Self-Negativity as a Function of the Medial Prefrontal Cortex

Self and emotions are key motivational factors of a person strivings for health and well-being. Understanding neural mechanisms supporting the relationship between these factors bear far-reaching implications for mental health disorders. Recent work indicates a substantial overlap between self-relevant and emotion information processing and has proposed the medial prefrontal cortex (MPFC) as one shared neural signature. However, the precise cognitive and neural mechanisms represented by the MPFC in investigations of self- and emotion-related processing are largely unknown. Here we examined whether the neural underpinnings of self-related processing in the MPFC link to positive or negative emotions. We collected fMRI data to test the distinct and shared neural circuits of self- and emotion-related processing while participants performed personal (self, friend, or stranger) and emotion (happy, sad, or neutral) associative matching tasks. By exploiting tight control over the factors that determine the effects of self-relevance and emotions (positive: Happy vs. neutral; negative: Sad vs. neutral), our univariate analysis revealed that the ventral part of the MPFC (vmPFC), which has established involvement in self-prioritisation effects, was not recruited in the negative emotion prioritisation effect. In contrast, there were no differences in brain activity between the effects of positive emotion- and self-prioritisation. These results were replicated by both region of interest (ROI)-based analysis in the vmPFC and the seed- to voxel functional connectivity analysis between the MPFC and the rest of the brain. The results suggest that the prioritisation effects for self and positive emotions are tightly linked together, and the MPFC plays a large role in discriminating between positive and negative emotions in relation to self-relevance.

The variability of functional MRI brain signal increases in Alzheimer's disease at cardiorespiratory frequencies

Biomarkers sensitive to prodromal or early pathophysiological changes in Alzheimer's disease (AD) symptoms could improve disease detection and enable timely interventions. Changes in brain hemodynamics may be associated with the main clinical AD symptoms. To test this possibility, we measured the variability of blood oxygen level-dependent (BOLD) signal in individuals from three independent datasets (totaling 80 AD patients and 90 controls). We detected a replicable increase in brain BOLD signal variability in the AD populations, which constituted a robust biomarker for clearly differentiating AD cases from controls. Fast BOLD scans showed that the elevated BOLD signal variability in AD arises mainly from cardiovascular brain pulsations. Manifesting in abnormal cerebral perfusion and cerebrospinal fluid convection, present observation presents a mechanism explaining earlier observations of impaired glymphatic clearance associated with AD in humans.

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.

Real-Time Resting-State Functional Magnetic Resonance Imaging Using Averaged Sliding Windows with Partial Correlations and Regression of Confounding Signals

Background/Introduction: There is considerable interest in using real-time functional magnetic resonance imaging (fMRI) for monitoring functional connectivity dynamics. To date, the majority of real-time resting-state fMRI studies have examined a limited number of brain regions. This is, in part, due to the computational demands of traditional seed- and independent component analysis-based methods, in particular when using increasingly available high-speed fMRI methods. Methods: This study describes a computationally efficient, real-time, seed-based, resting-state fMRI analysis pipeline using moving averaged sliding-windows (ASW) with partial correlations and regression of motion parameters and signals from white matter and cerebrospinal fluid. Results: Analytical and numerical analyses of ASW correlation and sliding-window regression as a function of window width show selectable bandpass filter characteristics and effective suppression of artifactual correlations resulting from signal drifts and transients. The analysis pipeline is compatible with multislab echo-volumar imaging and simultaneous multislice echo-planar imaging with repetition times as short as 136 msec. High-speed, resting-state fMRI data in healthy controls demonstrate the effectiveness of this approach for minimizing artifactual correlations in white and gray matter, which was comparable to conventional regression across the entire scan. Integrating sliding-window averaging (width: W₁) within a second-level sliding-window (width: W₂) enabled monitoring of intra- and internetwork correlation dynamics of up to 12 resting-state networks with bandpass filter characteristics determined by the first-level sliding-window and temporal resolution W₁ + W₂. Conclusions: The computational performance and confound tolerance make this seed-based, resting-state fMRI approach suitable for real-time monitoring of data quality and resting-state connectivity dynamics in neuroscience and clinical research studies. Impact statement Using averaged sliding-windows for seed-based correlation and regression of confounding signals provides a powerful model-free approach to increase tolerance to artifactual signal transients in resting-state analysis. The algorithmic efficiency of this sliding-window approach enables real-time, seed-based, resting-state functional magnetic resonance imaging (fMRI) of multiple networks with computation of connectivity matrices and online monitoring of data quality. Integration of a second-level sliding-window enables mapping of resting-state connectivity dynamics. Sensitivity and tolerance to confounding signals compare favorably with conventional correlation and confound regression across the entire scan. This methodological advance has the potential to enhance the clinical utility of resting-state fMRI and facilitate neuroscience applications.

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.

Transient Arousal Modulations Contribute to Resting-State Functional Connectivity Changes Associated with Head Motion Parameters

Correlations of resting-state functional magnetic resonance imaging (rsfMRI) signals are being widely used for assessing the functional brain connectivity in health and disease. However, an association was recently observed between rsfMRI connectivity modulations and the head motion parameters and regarded as a causal relationship, which has raised serious concerns about the validity of many rsfMRI findings. Here, we studied the origin of this rsfMRI-motion association and its relationship to arousal modulations. By using a template-matching method to locate arousal-related fMRI changes, we showed that the effects of high motion time points on rsfMRI connectivity are largely due to their significant overlap with arousal-affected time points. The finding suggests that the association between rsfMRI connectivity and the head motion parameters arises from their comodulations at transient arousal modulations, and this information is critical not only for proper interpretation of motion-associated rsfMRI connectivity changes, but also for controlling the potential confounding effects of arousal modulation on rsfMRI metrics.

From Light to Sound: Photoacoustic and Ultrasound Imaging in Fundamental Research of Alzheimer's Disease

Alzheimer's disease (AD) causes severe cognitive dysfunction and has long been studied for the underlining physiological and pathological mechanisms. Several biomedical imaging modalities have been applied, including MRI, PET, and high-resolution optical microscopy, for research purposes. However, there is still a strong need for imaging tools that can provide high spatiotemporal resolutions with relatively deep penetration to enhance our understanding of AD pathology and monitor treatment progress in fundamental research. Photoacoustic (PA) imaging and ultrasound (US) imaging can potentially address these unmet needs in AD research. PA imaging provides functional information with endogenous and/or exogenous contrast, while US imaging provides structural information. Recent studies have demonstrated the ability to monitor physiological parameters in small-animal brains with PA and US imaging as well as the feasibility of using US imaging as a therapeutic tool for AD. This concise review aims to introduce recent advances in AD research using PA and US imaging, provide the fundamentals, and discuss the potentials and challenges for future advances.

Presurgical Localization of the Primary Sensorimotor Cortex in Gliomas : When is Resting State FMRI Beneficial and Sufficient?

PURPOSE

Functional magnetic resonance imaging (fMRI) has an established role in neurosurgical planning; however, ambiguity surrounds the comparative value of resting and task-based fMRI relative to anatomical localization of the sensorimotor cortex. This study was carried out to determine: 1) how often fMRI adds to prediction of motor risks beyond expert neuroradiological review, 2) success rates of presurgical resting and task-based sensorimotor mapping, and 3) the impact of accelerated resting fMRI acquisitions on network detectability.

METHODS

Data were collected at 2 centers from 71 patients with a primary brain tumor (31 women; mean age 41.9 ± 13.9 years) and 14 healthy individuals (6 women; mean age 37.9 ± 12.7 years). Preoperative 3T MRI included anatomical scans and resting fMRI using unaccelerated (TR = 3.5 s), intermediate (TR = 1.56 s) or high temporal resolution (TR = 0.72 s) sequences. Task fMRI finger tapping data were acquired in 45 patients. Group differences in fMRI reproducibility, spatial overlap and success frequencies were assessed with t‑tests and χ²-tests.

RESULTS

Radiological review identified the central sulcus in 98.6% (70/71) patients. Task-fMRI succeeded in 100% (45/45). Resting fMRI failed to identify a sensorimotor network in up to 10 patients; it succeeded in 97.9% (47/48) of accelerated fMRIs, compared to only 60.9% (14/23) of unaccelerated fMRIs ([Formula: see text](2) = 17.84, p < 0.001). Of the patients 12 experienced postoperative deterioration, largely predicted by anatomical proximity to the central sulcus.

CONCLUSION

The use of fMRI in patients with residual or intact presurgical motor function added value to uncertain anatomical localization in just a single peri-Rolandic glioma case. Resting fMRI showed high correspondence to task localization when acquired with accelerated sequences but offered limited success at standard acquisitions.

Breath hold effect on cardiovascular brain pulsations - A multimodal magnetic resonance encephalography study

Ultra-fast functional magnetic resonance encephalography (MREG) enables separate assessment of cardiovascular, respiratory, and vasomotor waves from brain pulsations without temporal aliasing. We examined effects of breath hold- (BH) related changes on cardiovascular brain pulsations using MREG to study the physiological nature of cerebrovascular reactivity. We used alternating 32 s BH and 88 s resting normoventilation (NV) to change brain pulsations during MREG combined with simultaneously measured respiration, continuous non-invasive blood pressure, and cortical near-infrared spectroscopy (NIRS) in healthy volunteers. Changes in classical resting-state network BOLD-like signal and cortical blood oxygenation were reproduced based on MREG and NIRS signals. Cardiovascular pulsation amplitudes of MREG signal from anterior cerebral artery, oxygenated hemoglobin concentration in frontal cortex, and blood pressure decreased after BH. MREG cardiovascular pulse amplitudes in cortical areas and sagittal sinus increased, while cerebrospinal fluid and white matter remained unchanged. Respiratory centers in the brainstem - hypothalamus - thalamus - amygdala network showed strongest increases in cardiovascular pulsation amplitude. The spatial propagation of averaged cardiovascular impulses altered as a function of successive BH runs. The spread of cardiovascular pulse cycles exhibited a decreasing spatial similarity over time. MREG portrayed spatiotemporally accurate respiratory network activity and cardiovascular pulsation dynamics related to BH challenges at an unpreceded high temporal resolution.

Real-time presurgical resting-state fMRI in patients with brain tumors: Quality control and comparison with task-fMRI and intraoperative mapping

Resting-state functional magnetic resonance imaging (rsfMRI) is a promising task-free functional imaging approach, which may complement or replace task-based fMRI (tfMRI) in patients who have difficulties performing required tasks. However, rsfMRI is highly sensitive to head movement and physiological noise, and validation relative to tfMRI and intraoperative electrocortical mapping is still necessary. In this study, we investigate (a) the feasibility of real-time rsfMRI for presurgical mapping of eloquent networks with monitoring of data quality in patients with brain tumors and (b) rsfMRI localization of eloquent cortex compared with tfMRI and intraoperative electrocortical stimulation (ECS) in retrospective analysis. Five brain tumor patients were studied with rsfMRI and tfMRI on a clinical 3T scanner using MultiBand(8)-echo planar imaging (EPI) with repetition time: 400 ms. Moving-averaged sliding-window correlation analysis with regression of motion parameters and signals from white matter and cerebrospinal fluid was used to map sensorimotor and language resting-state networks. Data quality monitoring enabled rapid optimization of scan protocols, early identification of task noncompliance, and head movement-related false-positive connectivity to determine scan continuation or repetition. Sensorimotor and language resting-state networks were identifiable within 1 min of scan time. The Euclidean distance between ECS and rsfMRI connectivity and task-activation in motor cortex, Broca's, and Wernicke's areas was 5-10 mm, with the exception of discordant rsfMRI and ECS localization of Wernicke's area in one patient due to possible cortical reorganization and/or altered neurovascular coupling. This study demonstrates the potential of real-time high-speed rsfMRI for presurgical mapping of eloquent cortex with real-time data quality control, and clinically acceptable concordance of rsfMRI with tfMRI and ECS localization.

Dynamic lag analysis reveals atypical brain information flow in autism spectrum disorder

This study investigated whole‐brain dynamic lag pattern variations between neurotypical (NT) individuals and individuals with autism spectrum disorder (ASD) by applying a novel technique called dynamic lag analysis (DLA). The use of 3D magnetic resonance encephalography data with repetition time = 100 msec enables highly accurate analysis of the spread of activity between brain networks. Sixteen resting‐state networks (RSNs) with the highest spatial correlation between NT individuals (n = 20) and individuals with ASD (n = 20) were analyzed. The dynamic lag pattern variation between each RSN pair was investigated using DLA, which measures time lag variation between each RSN pair combination and statistically defines how these lag patterns are altered between ASD and NT groups. DLA analyses indicated that 10.8% of the 120 RSN pairs had statistically significant (P‐value <0.003) dynamic lag pattern differences that survived correction with surrogate data thresholding. Alterations in lag patterns were concentrated in salience, executive, visual, and default‐mode networks, supporting earlier findings of impaired brain connectivity in these regions in ASD. 92.3% and 84.6% of the significant RSN pairs revealed shorter mean and median temporal lags in ASD versus NT, respectively. Taken together, these results suggest that altered lag patterns indicating atypical spread of activity between large‐scale functional brain networks may contribute to the ASD phenotype. Autism Res 2020, 13: 244–258. © 2019 The Authors. Autism Research published by International Society for Autism Research published by Wiley Periodicals, Inc.

Lay Summary

Autism spectrum disorder (ASD) is characterized by atypical neurodevelopment. Using an ultra‐fast neuroimaging procedure, we investigated communication across brain regions in adults with ASD compared with neurotypical (NT) individuals. We found that ASD individuals had altered information flow patterns across brain regions. Atypical patterns were concentrated in salience, executive, visual, and default‐mode network areas of the brain that have previously been implicated in the pathophysiology of the disorder.

A 3D k-space Fourier encoding and reconstruction framework for simultaneous multi-slab acquisition

PURPOSE

To propose a novel 3D k-space Fourier encoding and reconstruction framework for simultaneous multi-slab (SMSlab) acquisition and demonstrate its efficacy in high-resolution imaging.

METHODS

First, it is illustrated in theory how the inter-slab gap interferes with the formation of the SMSlab 3D k-space. Then, joint RF and gradient encoding are applied to remove the inter-slab gap interference and form a SMSlab 3D k-space. In vivo experiments are performed to validate the proposed theory. Acceleration in the proposed SMSlab 3D k-space is also evaluated.

RESULTS

High-resolution (1.0 mm isotropic) images can be reconstructed using the proposed SMSlab 3D framework. Controlled aliasing in parallel imaging sampling and 2D GRAPPA reconstruction can also be applied in the SMSlab 3D k-space. Compared with conventional multi-slab acquisition, SMSlab exhibits better SNR maintainability (such as lower g-factors), especially at high acceleration factors.

CONCLUSION

It is demonstrated that the joint application of RF and gradient encoding enables SMSlab within a 3D Fourier encoding framework. Images with high isotropic resolution can be reconstructed, and further acceleration is also applicable. The proposed SMSlab 3D k-space can be valuable for both high-resolution and high-efficiency diffusion and functional MRI.

3D Multi-Resolution Optical Flow Analysis of Cardiovascular Pulse Propagation in Human Brain

The brain is cleaned from waste by glymphatic clearance serving a similar purpose as the lymphatic system in the rest of the body. Impairment of the glymphatic brain clearance precedes protein accumulation and reduced cognitive function in Alzheimer's disease (AD). Cardiovascular pulsations are a primary driving force of the glymphatic brain clearance. We developed a method to quantify cardiovascular pulse propagation in the human brain with magnetic resonance encephalography (MREG). We extended a standard optical flow estimation method to three spatial dimensions, with a multi-resolution processing scheme. We added application-specific criteria for discarding inaccurate results. With the proposed method, it is now possible to estimate the propagation of cardiovascular pulse wavefronts from the whole brain MREG data sampled at 10 Hz. The results show that on average the cardiovascular pulse propagates from major arteries via cerebral spinal fluid spaces into all tissue compartments in the brain. We present an example, that cardiovascular pulsations are significantly altered in AD: coefficient of variation and sample entropy of the pulse propagation speed in the lateral ventricles change in AD. These changes are in line with the theory of glymphatic clearance impairment in AD. The proposed non-invasive method can assess a performance indicator related to the glymphatic clearance in the human brain.

Beyond T2 and 3T: New MRI techniques for clinicians

Technological advances in Magnetic Resonance Imaging (MRI) in terms of field strength and hybrid MR systems have led to improvements in tumor imaging in terms of anatomy and functionality. This review paper discusses the applications of such advances in the field of radiation oncology with regards to treatment planning, therapy guidance and monitoring tumor response and predicting outcome.

Functional Connectivity in Adult Brain Tumor Patients: A Systematic Review

Brain tumor (BT) patients often experience reduced cognitive abilities and disrupted adaptive functioning before and after treatment. An innovative approach to understanding the underlying brain networks associated with these outcomes has been to study the brain's functional connectivity (FC), the spatially distributed and temporally correlated activity throughout the brain, and how it can be affected by a tumor. The present review synthesized the extant BT FC literature that utilizes functional magnetic resonance imaging to study FC strength of commonly observed networks during rest and task. A systematic review of English articles using PubMed was conducted. Search terms included brain tumor OR glioma AND functional connectivity, independent component analysis, ICA, psychophysiological interaction, OR PPI. Studies in which participants were diagnosed with BTs as adults that evaluated specific networks of interest using independent component analysis or seed-based component analysis were included. Twenty-five studies met inclusion criteria. BT patients often presented with decreases in FC strength within well-established networks and increases in atypical FC patterns. Network differences were tumor adjacent and distal, and left hemisphere tumors generally had a greater impact on FC. FC alterations often correlated with behavioral or cognitive outcomes when assessed. Overall, BTs appear to lead to various alterations in FC across different functional networks, and the most common change is a decrease in expected FC strength. More longitudinal studies are needed to determine the time course of network alterations across treatment and recovery, the role of medical treatments in BT survivors' FC, and the potential of FC patterns as biomarkers of cognitive outcomes.

Characterizing Rapid Fluctuations of Resting State Functional Connectivity in Demyelinating, Neurodegenerative, and Psychiatric Conditions: From Static to Time-Varying Analysis

Functional magnetic resonance imaging (fMRI) at resting state (RS) has been widely used to characterize the main brain networks. Functional connectivity (FC) has been mostly assessed assuming that FC is static across the whole fMRI examination. However, FC is highly variable at a very fast time-scale, as demonstrated by neurophysiological techniques. Time-varying functional connectivity (TVC) is a novel approach that allows capturing reoccurring patterns of interaction among functional brain networks. Aim of this review is to provide a description of the methods currently used to assess TVC on RS fMRI data, and to summarize the main results of studies applying TVC in healthy controls and patients with multiple sclerosis (MS). An overview of the main results obtained in neurodegenerative and psychiatric conditions is also provided. The most popular TVC approach is based on the so-called “sliding windows,” in which the RS fMRI acquisition is divided in small temporal segments (windows). A window of fixed length is shifted over RS fMRI time courses, and data within each window are used to calculate FC and its variability over time. Sliding windows can be combined with clustering techniques to identify recurring FC states or used to assess global TVC properties of large-scale functional networks or specific brain regions. TVC studies have used heterogeneous methodologies so far. Despite this, similar results have been obtained across investigations. In healthy subjects, the default-mode network (DMN) exhibited the highest degree of connectivity dynamism. In MS patients, abnormal global TVC properties and TVC strengths were found mainly in sensorimotor, DMN and salience networks, and were associated with more severe structural MRI damage and with more severe physical and cognitive disability. Conversely, abnormal TVC measures of the temporal network were correlated with better cognitive performances and less severe fatigue. In patients with neurodegenerative and psychiatric conditions, TVC abnormalities of the DMN, attention and executive networks were associated to more severe clinical manifestations. TVC helps to provide novel insights into fundamental properties of functional networks, and improves the understanding of brain reorganization mechanisms. Future technical advances might help to clarify TVC association with disease prognosis and response to treatment.

The potential of MR-Encephalography for BCI/Neurofeedback applications with high temporal resolution

Real-time functional magnetic resonance imaging (rt-fMRI) enables the update of various brain-activity measures during an ongoing experiment as soon as a new brain volume is acquired. However, the recorded Blood-oxygen-level dependent (BOLD) signal also contains physiological artifacts such as breathing and heartbeat, which potentially cause misleading false positive effects especially problematic in brain-computer interface (BCI) and neurofeedback (NF) setups. The low temporal resolution of echo planar imaging (EPI) sequences (which is in the range of seconds) prevents a proper separation of these artifacts from the BOLD signal. MR-Encephalography (MREG) has been shown to provide the high temporal resolution required to unalias and correct for physiological fluctuations and leads to increased specificity and sensitivity for mapping task-based activation and functional connectivity as well as for detecting dynamic changes in connectivity over time. By comparing a simultaneous multislice echo planar imaging (SMS-EPI) sequence and an MREG sequence using the same nominal spatial resolution in an offline analysis for three different experimental fMRI paradigms (perception of house and face stimuli, motor imagery, Stroop task), the potential of this novel technique for future BCI and NF applications was investigated. First, adapted general linear model pre-whitening which accounts for the high temporal resolution in MREG was implemented to calculate proper statistical results and be able to compare these with the SMS-EPI sequence. Furthermore, the respiration- and cardiac pulsation-related signals were successfully separated from the MREG signal using independent component analysis which were then included as regressors for a GLM analysis. Only the MREG sequence allowed to clearly separate cardiac pulsation and respiration components from the signal time course. It could be shown that these components highly correlate with the recorded respiration and cardiac pulsation signals using a respiratory belt and fingertip pulse plethysmograph. Temporal signal-to-noise ratios of SMS-EPI and MREG were comparable. Functional connectivity analysis using partial correlation showed a reduced standard error in MREG compared to SMS-EPI. Also, direct time course comparisons by down-sampling the MREG signal to the SMS-EPI temporal resolution showed lower variance in MREG. In general, we show that the higher temporal resolution is beneficial for fMRI time course modeling and this aspect can be exploited in offline application but also, is especially attractive, for real-time BCI and NF applications.

Sampling Rate Effects on Resting State fMRI Metrics

Low image sampling rates used in resting state functional magnetic resonance imaging (rs-fMRI) may cause aliasing of the cardiorespiratory pulsations over the very low frequency (VLF) BOLD signal fluctuations which reflects to functional connectivity (FC). In this study, we examine the effect of sampling rate on currently used rs-fMRI FC metrics. Ultra-fast fMRI magnetic resonance encephalography (MREG) data, sampled with TR 0.1 s, was downsampled to different subsampled repetition times (sTR, range 0.3–3 s) for comparisons. Echo planar k-space sampling (TR 2.15 s) and interleaved slice collection schemes were also compared against the 3D single shot trajectory at 2.2 s sTR. The quantified connectivity metrics included stationary spatial, time, and frequency domains, as well as dynamic analyses. Time domain methods included analyses of seed-based functional connectivity, regional homogeneity (ReHo), coefficient of variation, and spatial domain group level probabilistic independent component analysis (ICA). In frequency domain analyses, we examined fractional and amplitude of low frequency fluctuations. Aliasing effects were spatially and spectrally analyzed by comparing VLF (0.01–0.1 Hz), respiratory (0.12–0.35 Hz) and cardiac power (0.9–1.3 Hz) FFT maps at different sTRs. Quasi-periodic pattern (QPP) of VLF events were analyzed for effects on dynamic FC methods. The results in conventional time and spatial domain analyses remained virtually unchanged by the different sampling rates. In frequency domain, the aliasing occurred mainly in higher sTR (1–2 s) where cardiac power aliases over respiratory power. The VLF power maps suffered minimally from increasing sTRs. Interleaved data reconstruction induced lower ReHo compared to 3D sampling (p < 0.001). Gradient recalled echo-planar imaging (EPI BOLD) data produced both better and worse metrics. In QPP analyses, the repeatability of the VLF pulse detection becomes linearly reduced with increasing sTR. In conclusion, the conventional resting state metrics (e.g., FC, ICA) were not markedly affected by different TRs (0.1–3 s). However, cardiorespiratory signals showed strongest aliasing in central brain regions in sTR 1–2 s. Pulsatile QPP and other dynamic analyses benefit linearly from short TR scanning.

FMRI correlates of olfactory processing in typically-developing school-aged children

Human olfactory processing is understudied relative to other sensory modalities, despite its links to neurodevelopmental and neurodegenerative disorders. To address this limitation, we developed a fast, robust fMRI odor paradigm that is appropriate for all ages and levels of cognitive functioning. To test this approach, thirty-four typically developing children aged 7-12 underwent fMRI during brief, repeated exposure to phenylethyl alcohol, a flower-scented odor. Prior to fMRI scanning, olfactory testing (odor detection and identification) was conducted. During fMRI stimulus presentation, odorant release was synchronized to each participant's inspiratory phase to ensure participants were inhaling during the odorant exposure. Between group differences and correlations between activation and odor detection threshold scores were tested using the FMRIB Software Library. Results demonstrated that our 2-min paradigm significantly activated primary and secondary olfactory regions. In addition, a significant relationship between odor detection threshold and higher activation in the right amygdala and lower activation in the left frontal, insular, occipital, and cerebellar regions was observed, suggesting that this approach is sensitive to individual differences in olfactory processing. These findings demonstrate the feasibility of studying olfactory function in children using brain imaging techniques.

Optimization of rs-fMRI parameters in the Seed Correlation Analysis (SCA) in DPARSF toolbox: A preliminary study

There are a number of various methods of resting-state functional magnetic resonance imaging (rs-fMRI) analysis such as independent component analysis, multivariate autoregressive models, or seed correlation analysis however their results depend on arbitrary choice of parameters. Therefore, the aim of this work was to optimize the parameters in the seed correlation analysis using the Data Processing Assistant for Resting-State fMRI (DPARSF) toolbox for rs-fMRI data received from a Siemens Magnetom Skyra 3-Tesla scanner using a whole-brain, gradient-echo echo planar sequence with a 32-channel head coil. Different ranges of the following parameters: amplitude of low-frequency fluctuation (ALFF), Gaussian kernel at FWHM and radius of spherical ROI for 109 regions were tested for 20 healthy volunteers. The highest values of functional connectivity (FC) correlations were found for ALFF 0.01-0.08, spherical ROIs with the 8-mm radius and Gaussian kernel 8 mm at FWHM in all the studied areas that is, Auditory, Sensimotor, Visual, and Default Mode Network. The dominating influence of ALFF and smoothing on values of FC correlations was noted.

Capturing Dynamic Connectivity from Resting State fMRI using Time-Varying Graphical Lasso

Functional connectivity (FC) within the human brain evaluated through functional magnetic resonance imaging (fMRI) data has attracted increasing attention and has been employed to study the development of the brain or health conditions of the brain. Many different approaches have been proposed to estimate FC from fMRI data, whereas many of them rely on an implicit assumption that functional connectivity should be static throughout the fMRI scan session. Recently, the fMRI community has realized the limitation of assuming static connectivity and dynamic approaches are more prominent in the resting state fMRI (rs-fMRI) analysis. The sliding window technique has been widely used in many studies to capture network dynamics, but has a number of limitations. In this study, we apply a time-varying graphical lasso (TVGL) model, an extension from the traditional graphical lasso, to address the challenge, which can greatly improve the estimation of FC. The performance of estimating dynamic FC is evaluated with the TVGL through both simulated experiments and real rs-fMRI data from the Philadelphia Neurodevelopmental Cohort (PNC) project. Improved performance is achieved over the sliding window technique. In particular, group differences and transition behaviours between young adults and children are investigated using the estimated dynamic connectivity networks, which help us to better unveil the mechanisms underlying the evolution of the brain over time.

Identifying Respiration-Related Aliasing Artifacts in the Rodent Resting-State fMRI

Resting-state functional magnetic resonance imaging (rs-fMRI) combined with optogenetics and electrophysiological/calcium recordings in animal models is becoming a popular platform to investigate brain dynamics under specific neurological states. Physiological noise originating from the cardiac and respiration signal is the dominant interference in human rs-fMRI and extensive efforts have been made to reduce these artifacts from the human data. In animal fMRI studies, physiological noise sources including the respiratory and cardiorespiratory artifacts to the rs-fMRI signal fluctuation have typically been less investigated. In this article, we demonstrate evidence of aliasing effects into the low-frequency rs-fMRI signal fluctuation mainly due to respiration-induced B0 offsets in anesthetized rats. This aliased signal was examined by systematically altering the fMRI sampling rate, i.e., the time of repetition (TR), in free-breathing conditions and by adjusting the rate of ventilation. Anesthetized rats under ventilation showed a significantly narrower frequency bandwidth of the aliasing effect than free-breathing animals. It was found that the aliasing effect could be further reduced in ventilated animals with a muscle relaxant. This work elucidates the respiration-related aliasing effects on the rs-fMRI signal fluctuation from anesthetized rats, indicating non-negligible physiological noise needed to be taken care of in both awake and anesthetized animal rs-fMRI studies.

Targeted partial reconstruction for real-time fMRI with arbitrary trajectories

PURPOSE

A partial image reconstruction formalism is introduced for the targeted extraction of real-time feedback from arbitrary trajectories when full image reconstruction in real time is computationally too demanding.

METHODS

Explicit calculation and storage of linear combinations of lines of the reconstruction matrix by an incomplete basis change in spatial coordinates lead to translation of the expensive full reconstruction from a frame-wise application to a region of interest (ROI)-wise application. This step is independent from signal data and can be executed before the experiment. Subsequently, the results of the sum over fully reconstructed voxels can be evaluated directly. Data from a high-speed fMRI acquisition was used to investigate the targeted partial reconstruction of a functional ROI atlas, incorporating an intravolume dephasing correction. The same data and ROIs were used for a comparison of the time series obtained with those obtained from already existing methods for compartment-wise reconstruction. To examine real-time feasibility, the reconstruction was implemented and tested for online reconstruction performance.

RESULTS

The reconstruction yields results that are virtually identical to the standard reconstruction (i.e., the magnitude sums over the ROIs), with negligible discrepancies even after termination of the conjugate gradient algorithm at a feasible number of iterations. Notably, more discrepancies arise with existing compartment-wise reconstructions. The online real-time implementation evaluated 1 ROI within 2.8 ms in the case of a highly parallel 3D whole brain acquisition.

CONCLUSION

The high reconstruction fidelity and speed are satisfying for the exemplary application of real-time functional feedback using a highly parallel 3D whole brain acquisition.

Fluctuations of the EEG-fMRI correlation reflect intrinsic strength of functional connectivity in default mode network

Both functional magnetic resonance imaging (fMRI) and electrophysiological recordings have revealed that resting-state functional connectivity is temporally variable in human brain. Combined full-band electroencephalography-fMRI (fbEEG-fMRI) studies have shown that infraslow (<.1 Hz) fluctuations in EEG scalp potential are correlated with the blood-oxygen-level-dependent (BOLD) fMRI signals and that also this correlation appears variable over time. Here, we used simultaneous fbEEG-fMRI to test the hypothesis that correlation dynamics between BOLD and fbEEG signals could be explained by fluctuations in the activation properties of resting-state networks (RSNs) such as the extent or strength of their activation. We used ultrafast magnetic resonance encephalography (MREG) fMRI to enable temporally accurate and statistically robust short-time-window comparisons of infra-slow fbEEG and BOLD signals. We found that the temporal fluctuations in the fbEEG-BOLD correlation were dependent on RSN connectivity strength, but not on the mean signal level or magnitude of RSN activation or motion during scanning. Moreover, the EEG-fMRI correlations were strongest when the intrinsic RSN connectivity was strong and close to the pial surface. Conversely, weak fbEEG-BOLD correlations were attributable to periods of less coherent or spatially more scattered intrinsic RSN connectivity, or RSN activation in deeper cerebral structures. The results thus show that the on-average low correlations between infra-slow EEG and BOLD signals are, in fact, governed by the momentary coherence and depth of the underlying RSN activation, and may reach systematically high values with appropriate source activities. These findings further consolidate the notion of slow scalp potentials being directly coupled to hemodynamic fluctuations.

Estimation of Dynamic Sparse Connectivity Patterns From Resting State fMRI

Functional connectivity (FC) estimated from functional magnetic resonance imaging (fMRI) time series, especially during resting state periods, provides a powerful tool to assess human brain functional architecture in health, disease, and developmental states. Recently, the focus of connectivity analysis has shifted toward the subnetworks of the brain, which reveals co-activating patterns over time. Most prior works produced a dense set of high-dimensional vectors, which are hard to interpret. In addition, their estimations to a large extent were based on an implicit assumption of spatial and temporal stationarity throughout the fMRI scanning session. In this paper, we propose an approach called dynamic sparse connectivity patterns (dSCPs), which takes advantage of both matrix factorization and time-varying fMRI time series to improve the estimation power of FC. The feasibility of analyzing dynamic FC with our model is first validated through simulated experiments. Then, we use our framework to measure the difference between young adults and children with real fMRI data set from the Philadelphia Neurodevelopmental Cohort (PNC). The results from the PNC data set showed significant FC differences between young adults and children in four different states. For instance, young adults had reduced connectivity between the default mode network and other subnetworks, as well as hyperconnectivity within the visual system in states 1 and 3, and hypoconnectivity in state 2. Meanwhile, they exhibited temporal correlation patterns that changed over time within functional subnetworks. In addition, the dSCPs model indicated that older people tend to spend more time within a relatively connected FC pattern. Overall, the proposed method provides a valid means to assess dynamic FC, which could facilitate the study of brain networks.

Time efficient whole-brain coverage with MR Fingerprinting using slice-interleaved echo-planar-imaging

Magnetic resonance fingerprinting (MRF) is a promising method for fast simultaneous quantification of multiple tissue parameters. The objective of this study is to improve the coverage of MRF based on echo-planar imaging (MRF-EPI) by using a slice-interleaved acquisition scheme. For this, the MRF-EPI is modified to acquire several slices in a randomized interleaved manner, increasing the effective repetition time of the spoiled gradient echo readout acquisition in each slice. Per-slice matching of the signal-trace to a precomputed dictionary allows the generation of T1 and T2* maps with integrated B1+ correction. Subsequent compensation for the coil sensitivity profile and normalization to the cerebrospinal fluid additionally allows for quantitative proton density (PD) mapping. Numerical simulations are performed to optimize the number of interleaved slices. Quantification accuracy is validated in phantom scans and feasibility is demonstrated in-vivo. Numerical simulations suggest the acquisition of four slices as a trade-off between quantification precision and scan-time. Phantom results indicate good agreement with reference measurements (Difference T1: -2.4 ± 1.1%, T2*: -0.5 ± 2.5%, PD: -0.5 ± 7.2%). In-vivo whole-brain coverage of T1, T2* and PD with 32 slices was acquired within 3:36 minutes, resulting in parameter maps of high visual quality and comparable performance with single-slice MRF-EPI at 4-fold scan-time reduction.

Influences of Age, Sex, and Moderate Alcohol Drinking on the Intrinsic Functional Architecture of Adolescent Brains

The transition from adolescent to adult cognition and emotional control requires neurodevelopmental maturation likely involving intrinsic functional networks (IFNs). Normal neurodevelopment may be vulnerable to disruption from environmental insult such as alcohol consumption commonly initiated during adolescence. To test potential disruption to IFN maturation, we used resting-state functional magnetic resonance imaging (rs-fMRI) in 581 no-to-low alcohol-consuming and 117 moderate-to-high-drinking youth. Functional seed-to-voxel connectivity analysis assessed age, sex, and moderate alcohol drinking on default-mode, executive-control, salience, reward, and emotion networks and tested cognitive and motor coordination correlates of network connectivity. Among no-to-low alcohol-consuming adolescents, executive-control frontolimbicstriatal connectivity was stronger in older than younger adolescents, particularly boys, and predicted better ability in balance, memory, and impulse control. Connectivity patterns in moderate-to-high-drinking youth were tested mainly in late adolescence when drinking was initiated. Implicated was the emotion network with attenuated connectivity to default-mode network regions. Our cross-sectional rs-fMRI findings from this large cohort of adolescents show sexual dimorphism in connectivity and suggest neurodevelopmental rewiring toward stronger and spatially more distributed executive-control networking in older than younger adolescents. Functional network rewiring in moderate-to-high-drinking adolescents may impede maturation of affective and self-reflection systems and obscure maturation of complex social and emotional behaviors.

Inner-volume echo volumar imaging (IVEVI) for robust fetal brain imaging

PURPOSE

Fetal functional MRI studies using conventional 2-dimensional single-shot echo-planar imaging sequences may require discarding a large data fraction as a result of fetal and maternal motion. Increasing the temporal resolution using echo volumar imaging (EVI) could provide an effective alternative strategy. Echo volumar imaging was combined with inner volume (IV) imaging (IVEVI) to locally excite the fetal brain and acquire full 3-dimensional images, fast enough to freeze most fetal head motion.

METHODS

IVEVI was implemented by modifying a standard multi-echo echo-planar imaging sequence. A spin echo with orthogonal excitation and refocusing ensured localized excitation. To introduce T2* weighting and to save time, the k-space center was shifted relative to the spin echo. Both single and multi-shot variants were tested. Acoustic noise was controlled by adjusting the amplitude and switching frequency of the readout gradient. Image-based shimming was used to minimize B₀ inhomogeneities within the fetal brain.

RESULTS

The sequence was first validated in an adult. Eight fetuses were scanned using single-shot IVEVI at a 3.5 × 3.5 × 5.0 mm³ resolution with a readout duration of 383 ms. Multishot IVEVI showed reduced geometric distortions along the second phase-encode direction.

CONCLUSIONS

Fetal EVI remains challenging. Although effective echo times comparable to the T2* values of fetal cortical gray matter at 3 T could be achieved, controlling acoustic noise required longer readouts, leading to substantial distortions in single-shot images. Although multishot variants enabled us to reduce susceptibility-induced geometric distortions, sensitivity to motion was increased. Future studies should therefore focus on improvements to multishot variants. Magn Reson Med 80:279-285, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Safety and EEG data quality of concurrent high-density EEG and high-speed fMRI at 3 Tesla

PURPOSE

Concurrent EEG and fMRI is increasingly used to characterize the spatial-temporal dynamics of brain activity. However, most studies to date have been limited to conventional echo-planar imaging (EPI). There is considerable interest in integrating recently developed high-speed fMRI methods with high-density EEG to increase temporal resolution and sensitivity for task-based and resting state fMRI, and for detecting interictal spikes in epilepsy. In the present study using concurrent high-density EEG and recently developed high-speed fMRI methods, we investigate safety of radiofrequency (RF) related heating, the effect of EEG on cortical signal-to-noise ratio (SNR) in fMRI, and assess EEG data quality.

MATERIALS AND METHODS

The study compared EPI, multi-echo EPI, multi-band EPI and multi-slab echo-volumar imaging pulse sequences, using clinical 3 Tesla MR scanners from two different vendors that were equipped with 64- and 256-channel MR-compatible EEG systems, respectively, and receive only array head coils. Data were collected in 11 healthy controls (3 males, age range 18-70 years) and 13 patients with epilepsy (8 males, age range 21-67 years). Three of the healthy controls were scanned with the 256-channel EEG system, the other subjects were scanned with the 64-channel EEG system. Scalp surface temperature, SNR in occipital cortex and head movement were measured with and without the EEG cap. The degree of artifacts and the ability to identify background activity was assessed by visual analysis by a trained expert in the 64 channel EEG data (7 healthy controls, 13 patients).

RESULTS

RF induced heating at the surface of the EEG electrodes during a 30-minute scan period with stable temperature prior to scanning did not exceed 1.0° C with either EEG system and any of the pulse sequences used in this study. There was no significant decrease in cortical SNR due to the presence of the EEG cap (p > 0.05). No significant differences in the visually analyzed EEG data quality were found between EEG recorded during high-speed fMRI and during conventional EPI (p = 0.78). Residual ballistocardiographic artifacts resulted in 58% of EEG data being rated as poor quality.

CONCLUSION

This study demonstrates that high-density EEG can be safely implemented in conjunction with high-speed fMRI and that high-speed fMRI does not adversely affect EEG data quality. However, the deterioration of the EEG quality due to residual ballistocardiographic artifacts remains a significant constraint for routine clinical applications of concurrent EEG-fMRI.

On the detection of high frequency correlations in resting state fMRI

Current studies of resting-state connectivity rely on coherent signal fluctuations at frequencies below 0.1 Hz, however, recent studies using high-speed fMRI have shown that fluctuations above 0.5 Hz may exist. This study replicates the feasibility of measuring high frequency (HF) correlations in six healthy controls and a patient with a brain tumor while analyzing non-physiological signal sources via simulation. Resting-state data were acquired using a high-speed multi-slab echo-volumar imaging pulse sequence with 136 ms temporal resolution. Bandpass frequency filtering in combination with sliding window seed-based connectivity analysis using running mean of the correlation maps was employed to map HF correlations up to 3.7 Hz. Computer simulations of Rician noise and the underlying point spread function were analyzed to estimate baseline spatial autocorrelation levels in four major networks (auditory, sensorimotor, visual, and default-mode). Using seed regions based on Brodmann areas, the auditory and default-mode networks were observed to have significant frequency band dependent HF correlations above baseline spatial autocorrelation levels. Correlations in the sensorimotor network were at trend level. The auditory network was still observed using a unilateral single voxel seed. In the patient, HF auditory correlations showed a spatial displacement near the tumor consistent with the displacement seen at low frequencies. In conclusion, our data suggest that HF connectivity in the human brain may be observable with high-speed fMRI, however, the detection sensitivity may depend on the network observed, data acquisition technique, and analysis method.

Biomarkers, designs, and interpretations of resting-state fMRI in translational pharmacological research: A review of state-of-the-Art, challenges, and opportunities for studying brain chemistry

A decade of research and development in resting-state functional MRI (RSfMRI) has opened new translational and clinical research frontiers. This review aims to bridge between technical and clinical researchers who seek reliable neuroimaging biomarkers for studying drug interactions with the brain. About 85 pharma-RSfMRI studies using BOLD signal (75% of all) or arterial spin labeling (ASL) were surveyed to investigate the acute effects of psychoactive drugs. Experimental designs and objectives include drug fingerprinting dose-response evaluation, biomarker validation and calibration, and translational studies. Common biomarkers in these studies include functional connectivity, graph metrics, cerebral blood flow and the amplitude and spectrum of BOLD fluctuations. Overall, RSfMRI-derived biomarkers seem to be sensitive to spatiotemporal dynamics of drug interactions with the brain. However, drugs cause both central and peripheral effects, thus exacerbate difficulties related to biological confounds, structured noise from motion and physiological confounds, as well as modeling and inference testing. Currently, these issues are not well explored, and heterogeneities in experimental design, data acquisition and preprocessing make comparative or meta-analysis of existing reports impossible. A unifying collaborative framework for data-sharing and data-mining is thus necessary for investigating the commonalities and differences in biomarker sensitivity and specificity, and establishing guidelines. Multimodal datasets including sham-placebo or active control sessions and repeated measurements of various psychometric, physiological, metabolic and neuroimaging phenotypes are essential for pharmacokinetic/pharmacodynamic modeling and interpretation of the findings. We provide a list of basic minimum and advanced options that can be considered in design and analyses of future pharma-RSfMRI studies. Hum Brain Mapp 38:2276-2325, 2017. © 2017 Wiley Periodicals, Inc.

Enhanced subject-specific resting-state network detection and extraction with fast fMRI

Resting-state networks have become an important tool for the study of brain function. An ultra-fast imaging technique that allows to measure brain function, called Magnetic Resonance Encephalography (MREG), achieves an order of magnitude higher temporal resolution than standard echo-planar imaging (EPI). This new sequence helps to correct physiological artifacts and improves the sensitivity of the fMRI analysis. In this study, EPI is compared with MREG in terms of capability to extract resting-state networks. Healthy controls underwent two consecutive resting-state scans, one with EPI and the other with MREG. Subject-level independent component analyses (ICA) were performed separately for each of the two datasets. Using Stanford FIND atlas parcels as network templates, the presence of ICA maps corresponding to each network was quantified in each subject. The number of detected individual networks was significantly higher in the MREG data set than for EPI. Moreover, using short time segments of MREG data, such as 50 seconds, one can still detect and track consistent networks. Fast fMRI thus results in an increased capability to extract distinct functional regions at the individual subject level for the same scan times, and also allow the extraction of consistent networks within shorter time intervals than when using EPI, which is notably relevant for the analysis of dynamic functional connectivity fluctuations. Hum Brain Mapp 38:817-830, 2017. © 2016 Wiley Periodicals, Inc.

Ultra-fast magnetic resonance encephalography of physiological brain activity - Glymphatic pulsation mechanisms?

The theory on the glymphatic convection mechanism of cerebrospinal fluid holds that cardiac pulsations in part pump cerebrospinal fluid from the peri-arterial spaces through the extracellular tissue into the peri-venous spaces facilitated by aquaporin water channels. Since cardiac pulses cannot be the sole mechanism of glymphatic propulsion, we searched for additional cerebrospinal fluid pulsations in the human brain with ultra-fast magnetic resonance encephalography. We detected three types of physiological mechanisms affecting cerebral cerebrospinal fluid pulsations: cardiac, respiratory, and very low frequency pulsations. The cardiac pulsations induce a negative magnetic resonance encephalography signal change in peri-arterial regions that extends centrifugally and covers the brain in ≈1 Hz cycles. The respiratory ≈0.3 Hz pulsations are centripetal periodical pulses that occur dominantly in peri-venous areas. The third type of pulsation was very low frequency (VLF 0.001-0.023 Hz) and low frequency (LF 0.023-0.73 Hz) waves that both propagate with unique spatiotemporal patterns. Our findings using critically sampled magnetic resonance encephalography open a new view into cerebral fluid dynamics. Since glymphatic system failure may precede protein accumulations in diseases such as Alzheimer's dementia, this methodological advance offers a novel approach to image brain fluid dynamics that potentially can enable early detection and intervention in neurodegenerative diseases.

Advantages and disadvantages of a fast fMRI sequence in the context of EEG-fMRI investigation of epilepsy patients: A realistic simulation study

EEG-fMRI is an established technique to allow mapping BOLD changes in response to interictal discharges recorded in the EEG of epilepsy patients. Traditional fMRI experiments rely on an echo planar imaging (EPI) sequence with a time resolution given by its time-to-repetition (TR) of ~2 s. Recently, multiple fast fMRI sequences have been developed to get around the limited temporal resolution of the EPI sequence, and achieved a TR in the 100 ms range or lower. One such sequence is called Magnetic Resonance EncephaloGraphy (MREG). Its high temporal resolution should offer increased detection sensitivity and statistical power in the context of epilepsy studies and in fMRI experiments in general. The aim of this work was to investigate the advantages and disadvantages offered by MREG. This was done by superimposing artificial event-related BOLD responses on EPI and MREG background signals, from 5 epileptic patients, that were free of epileptic discharges (spikes) on simultaneously recorded EEG. These functional datasets simulated different spiking rates and hemodynamic response amplitudes, and were analyzed with the commonly used General Linear Model (GLM) with the canonical hemodynamic response function (HRF) as a fixed model of the response. Robustness to violation of the assumptions of the GLM was additionally assessed with similar simulations using variable spike-to-spike response amplitudes and 8 non-canonical HRFs. Consistent with previous work, MREG yields higher maximum statistical t-values than EPI, but our simulations showed these statistics to be inflated, as the false positive rate at a standard threshold was high. At thresholds set to appropriately control specificity, EPI showed better true positive rate and larger cluster size than MREG. However, the lack of an appropriate calibration of the amplitude of the responses across the sequences precludes definitive judgment on their relative sensitivity. In addition, we show that a mismatch between the assumed and actual HRF impairs more MREG detection performance, but that EPI is more affected by non-modeled spike-to-spike variations of response amplitude. Filtering-out physiological noise, which is not aliased at the fast sampling rate of MREG, and the modeling of temporal autocorrelation are advantageous in increasing the detection power of MREG. This simulation study 1) warrants care when interpreting statistical t-values from fast fMRI sequences, 2) proposes thresholds for valid inferences and processing methods for maximal sensitivities, and 3) demonstrates the relative robustness/susceptibility of MREG and EPI to violation of the GLM's assumptions.

A proof-of-principle study of multi-site real-time functional imaging at 3T and 7T: Implementation and validation

Real-time functional Magnetic Resonance Imaging (rtfMRI) is used mainly for neurofeedback or for brain-computer interfaces (BCI). But multi-site rtfMRI could in fact help in the application of new interactive paradigms such as the monitoring of mutual information flow or the controlling of objects in shared virtual environments. For that reason, a previously developed framework that provided an integrated control and data analysis of rtfMRI experiments was extended to enable multi-site rtfMRI. Important new components included a data exchange platform for analyzing the data of both MR scanners independently and/or jointly. Information related to brain activation can be displayed separately or in a shared view. However, a signal calibration procedure had to be developed and integrated in order to permit the connecting of sites that had different hardware and to account for different inter-individual brain activation levels. The framework was successfully validated in a proof-of-principle study with twelve volunteers. Thus the overall concept, the calibration of grossly differing signals, and BCI functionality on each site proved to work as required. To model interactions between brains in real-time, more complex rules utilizing mutual activation patterns could easily be implemented to allow for new kinds of social fMRI experiments.

Windowed correlation: a suitable tool for providing dynamic fMRI-based functional connectivity neurofeedback on task difficulty

The goal of neurofeedback training is to provide participants with relevant information on their ongoing brain processes in order to enable them to change these processes in a meaningful way. Under the assumption of an intrinsic brain-behavior link, neurofeedback can be a tool to guide a participant towards a desired behavioral state, such as a healthier state in the case of patients. Current research in clinical neuroscience regarding the most robust indicators of pathological brain processes in psychiatric and neurological disorders indicates that fMRI-based functional connectivity measures may be among the most important biomarkers of disease. The present study therefore investigated the general potential of providing fMRI neurofeedback based on functional correlations, computed from short-window time course data at the level of single task periods. The ability to detect subtle changes in task performance with block-wise functional connectivity measures was evaluated based on imaging data from healthy participants performing a simple motor task, which was systematically varied along two task dimensions representing two different aspects of task difficulty. The results demonstrate that fMRI-based functional connectivity measures may provide a better indicator for an increase in overall (motor) task difficulty than activation level-based measures. Windowed functional correlations thus seem to provide relevant and unique information regarding ongoing brain processes, which is not captured equally well by standard activation level-based neurofeedback measures. Functional connectivity markers, therefore, may indeed provide a valuable tool to enhance and monitor learning within an fMRI neurofeedback setup.