Detection of short-term activity avalanches in human brain default mode network with ultrafast MR encephalography

Inferring brain causal and temporal-lag networks for recognizing abnormal patterns of dementia

Brain functional network analysis has become a popular method to explore the laws of brain organization and identify biomarkers of neurological diseases. However, it is still a challenging task to construct an ideal brain network due to the limited understanding of the human brain. Existing methods often ignore the impact of temporal-lag on the results of brain network modeling, which may lead to some unreliable conclusions. To overcome this issue, we propose a novel brain functional network estimation method, which can simultaneously infer the causal mechanisms and temporal-lag values among brain regions. Specifically, our method converts the lag learning into an instantaneous effect estimation problem, and further embeds the search objectives into a deep neural network model as parameters to be learned. To verify the effectiveness of the proposed estimation method, we perform experiments on the Alzheimer's Disease Neuroimaging Initiative (ADNI) database by comparing the proposed model with several existing methods, including correlation-based and causality-based methods. The experimental results show that our brain networks constructed by the proposed estimation method can not only achieve promising classification performance, but also exhibit some characteristics of physiological mechanisms. Our approach provides a new perspective for understanding the pathogenesis of brain diseases. The source code is released at https://github.com/NJUSTxiazw/CTLN.

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.

Co-activation pattern alterations in autism spectrum disorder-A volume-wise hierarchical clustering fMRI study

INTRODUCTION

There has been a growing effort to characterize the time-varying functional connectivity of resting state (RS) fMRI brain networks (RSNs). Although voxel-wise connectivity studies have examined different sliding window lengths, nonsequential volume-wise approaches have been less common.

METHODS

Inspired by earlier co-activation pattern (CAP) studies, we applied hierarchical clustering (HC) to classify the image volumes of the RS-fMRI data on 28 adolescents with autism spectrum disorder (ASD) and their 27 typically developing (TD) controls. We compared the distribution of the ASD and TD groups' volumes in CAPs as well as their voxel-wise means. For simplification purposes, we conducted a group independent component analysis to extract 14 major RSNs. The RSNs' average z-scores enabled us to meaningfully regroup the RSNs and estimate the percentage of voxels within each RSN for which there was a significant group difference. These results were jointly interpreted to find global group-specific patterns.

RESULTS

We found similar brain state proportions in 58 CAPs (clustering interval from 2 to 30). However, in many CAPs, the voxel-wise means differed significantly within a matrix of 14 RSNs. The rest-activated default mode-positive and default mode-negative brain state properties vary considerably in both groups over time. This division was seen clearly when the volumes were partitioned into two CAPs and then further examined along the HC dendrogram of the diversifying brain CAPs. The ASD group network activations followed a more heterogeneous distribution and some networks maintained higher baselines; throughout the brain deactivation state, the ASD participants had reduced deactivation in 12/14 networks. During default mode-negative CAPs, the ASD group showed simultaneous visual network and either dorsal attention or default mode network overactivation.

CONCLUSION

Nonsequential volume gathering into CAPs and the comparison of voxel-wise signal changes provide a complementary perspective to connectivity and an alternative to sliding window analysis.

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.

Lag Analysis of Fast fMRI Reveals Delayed Information Flow Between the Default Mode and Other Networks in Narcolepsy

Narcolepsy is a chronic neurological disease characterized by dysfunction of the hypocretin system in brain causing disruption in the wake-promoting system. In addition to sleep attacks and cataplexy, patients with narcolepsy commonly report cognitive symptoms while objective deficits in sustained attention and executive function have been observed. Prior resting-state functional magnetic resonance imaging (fMRI) studies in narcolepsy have reported decreased inter/intranetwork connectivity regarding the default mode network (DMN). Recently developed fast fMRI data acquisition allows more precise detection of brain signal propagation with a novel dynamic lag analysis. In this study, we used fast fMRI data to analyze dynamics of inter resting-state network (RSN) information signaling between narcolepsy type 1 patients (NT1, n = 23) and age- and sex-matched healthy controls (HC, n = 23). We investigated dynamic connectivity properties between positive and negative peaks and, furthermore, their anticorrelative (pos-neg) counterparts. The lag distributions were significantly (P < 0.005, familywise error rate corrected) altered in 24 RSN pairs in NT1. The DMN was involved in 83% of the altered RSN pairs. We conclude that narcolepsy type 1 is characterized with delayed and monotonic inter-RSN information flow especially involving anticorrelations, which are known to be characteristic behavior of the DMN regarding neurocognition.

Heart Rate and Respiration Affect the Functional Connectivity of Default Mode Network in Resting-State Functional Magnetic Resonance Imaging

A growing number of brain imaging studies show functional connectivity (FC) between regions during emotional and cognitive tasks in humans. However, emotions are accompanied by changes in physiological parameters such as heart rate and respiration. These changes may affect blood oxygen level-dependent signals, as well as connectivity between brain areas. This study aimed to clarify the effects of physiological noise on the connectivity between areas related to the default mode network using resting-state functional magnetic resonance imaging (rs-fMRI). Healthy adult volunteers (age range: 19–51 years, mean age: 26.9 ± 9.1 years, 8 males and 8 females) underwent rs-fMRI for 10 min using a clinical 3T scanner (MAGNETOM Trio A Tim System, Siemens) with simultaneously recorded respiration and cardiac output. Physiological noise signals were subsequently removed from the acquired fMRI data using the DRIFTER toolbox. Image processing and analysis of the FC between areas related to the default mode network were performed using DPARSF. Network-Based Statistic (NBS) analysis of the functional connectome of the DMN and DMN-related area was used to perform three groups of comparison: without physiological noise correction, with cardiac noise correction, and with cardiac and respiratory noise correction. NBS analysis identified 36 networks with significant differences in three conditions in FC matrices. Post hoc comparison showed no differences between the three conditions, indicating that all three had the same networks. Among the 36 networks, strength of FC of 8 networks was modified under physiological noise correction. Connectivity between left and right anterior medial frontal regions increased strength of connectivity. These areas are located on the medial cerebral hemisphere, close to the sagittal sinus and arteries in the cerebral hemispheres, suggesting that medial frontal areas may be sensitive to cardiac rhythm close to arteries. The other networks observed temporal regions and showed a decrease in their connectivity strength by removing physiological noise, indicating that physiological noise, especially respiration, may be sensitive to BOLD signal in the temporal regions during resting state. Temporal lobe was highly correlated with anxiety-related respiration changes (Masaoka and Homma, 2000), speech processing, and respiratory sensation. These factors may affect the rs-fMRI signaling sensitivity.

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.

Cerebrovascular blood oxygenation level dependent pulsatility at baseline and following acute exercise among healthy adolescents

Arterial stiffness is linked to cerebral small vessel damage and neurodegeneration, but barriers to accessing deep cerebrovascular anatomy limit our ability to assess the brain. This study describes an adaptation of a cardiac-related scrubbing method as a means of generating blood oxygenation level-dependent pulsatility maps based on the cardiac cycle. We examine BOLD pulsatility at rest, based on the non-parametric deviation from null metric, as well as changes following acute physiological stress from 20 min of moderate-intensity cycling in 45 healthy adolescents. We evaluate the influence of repetition time (TR) and echo time (TE) using simulated and multi-echo empirical data, respectively. There were tissue-specific and voxel-wise BOLD pulsatility decreases 20 min following exercise cessation. BOLD pulsatility detection was comparable over a range of TR and TE values when scan volumes were kept constant; however, short TRs (≤500 ms) and TEs (∼14 ms) acquisitions would yield the most efficient detection. Results suggest cardiac-related BOLD pulsatility may represent a robust and easily adopted method of mapping cerebrovascular pulsatility with voxel-wise resolution.

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.

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.

Multimodal brain imaging with magnetoencephalography: A method for measuring blood pressure and cardiorespiratory oscillations

Studies with magnetoencephalography (MEG) are still quite rarely combined simultaneously with methods that can provide a metabolic dimension to MEG investigations. In addition, continuous blood pressure measurements which comply with MEG compatibility requirements are lacking. For instance, by combining methods reflecting neurovascular status one could obtain more information on low frequency fluctuations that have recently gained increasing interest as a mediator of functional connectivity within brain networks. This paper presents a multimodal brain imaging setup, capable to non-invasively and continuously measure cerebral hemodynamic, cardiorespiratory and blood pressure oscillations simultaneously with MEG. In the setup, all methods apart from MEG rely on the use of fibre optics. In particular, we present a method for measuring of blood pressure and cardiorespiratory oscillations continuously with MEG. The potential of this type of multimodal setup for brain research is demonstrated by our preliminary studies on human, showing effects of mild hypercapnia, gathered simultaneously with the presented modalities.

Dynamic Functional Brain Connectivity for Face Perception

Face perception is mediated by a distributed brain network comprised of the core system at occipito-temporal areas and the extended system at other relevant brain areas involving bilateral hemispheres. In this study we explored how the brain connectivity changes over the time for face-sensitive processing. We investigated the dynamic functional connectivity in face perception by analyzing time-dependent EEG phase synchronization in four different frequency bands: theta (4-7 Hz), alpha (8-14 Hz), beta (15-24 Hz), and gamma (25-45 Hz) bands in the early stages of face processing from 30 to 300 ms. High-density EEG were recorded from subjects who were passively viewing faces, buildings, and chairs. The dynamic connectivity within the core system and between the extended system were investigated. Significant differences between faces and non-faces mainly appear in theta band connectivity: (1) at the time segment of 90-120 ms between parietal area and occipito-temporal area in the right hemisphere, and (2) at the time segment of 150-180 ms between bilateral occipito-temporal areas. These results indicate (1) the importance of theta-band connectivity in the face-sensitive processing, and (2) that different parts of network are involved for the initial stage of face categorization and the stage of face structural encoding.