Fast fMRI Research Articles

The temporal specificity of BOLD fMRI is systematically related to anatomical and vascular features of the human brain.

The ability to detect fast responses with functional MRI depends on the speed of hemodynamic responses to neural activity, because hemodynamic responses act as a temporal low-pass filter which blurs rapid changes. However, the shape and timing of hemodynamic responses are highly variable across the brain and across stimuli. This heterogeneity of responses implies that the temporal specificity of fMRI signals, or the ability of fMRI to preserve fast information, could also vary substantially across the cortex. In this work we investigated how local differences in hemodynamic response timing affect the temporal specificity of fMRI. We used ultra-high field (7T) fMRI at high spatiotemporal resolution, studying the primary visual cortex (V1) as a model area for investigation. We used visual stimuli oscillating at slow and fast frequencies to probe the temporal specificity of individual voxels. As expected, we identified substantial variability in temporal specificity, with some voxels preserving their responses to fast neural activity more effectively than others. We investigated which voxels had the highest temporal specificity, and tested whether voxel timing was related to anatomical and vascular features. We found that low temporal specificity is only weakly explained by the presence of large veins or cerebral cortical depth. Notably, however, temporal specificity depended strongly on a voxel's position along the anterior-posterior anatomical axis of V1, with voxels within the calcarine sulcus being capable of preserving close to 25% of their amplitude as the frequency of stimulation increased from 0.05Hz to 0.20Hz, and voxels nearest to the occipital pole preserving less than 18%. These results indicate that detection biases in high-resolution fMRI will depend on the anatomical and vascular features of the area being imaged, and that these biases will differ depending on the timing of the underlying neuronal activity. While we attribute this variance primarily to hemodynamic effects, neuronal nonlinearities may also influence response timing. Importantly, this spatial heterogeneity of temporal specificity suggests that it could be exploited to achieve higher specificity in some locations, and that tailored data analysis strategies may help improve the detection and interpretation of fast fMRI responses.

Fast 3D fMRI acquisition with high spatial resolutions over a reduced FOV.

To develop and demonstrate a fast 3D fMRI acquisition technique with high spatial resolution over a reduced FOV, named k-t 3D reduced FOV imaging (3D-rFOVI). Based on 3D gradient-echo EPI, k-t 3D-rFOVI used a 2D RF pulse to reduce the FOV in the in-plane phase-encoding direction, boosting spatial resolution without increasing echo train length. For image acceleration, full sampling was applied in the central k-space region along the through-slab direction (kz) for all time frames, while randomized undersampling was used in outer kz regions at different time frames. Images were acquired at 3T and reconstructed using a method based on partial separability. fMRI detection sensitivity of k-t 3D-rFOVI was quantitively analyzed with simulation data. Human visual fMRI experiments were performed to evaluate k-t 3D-rFOVI and compare it with a commercial multiband EPI sequence. The simulation data showed that k-t 3D-rFOVI can detect 100% of fMRI activations with an acceleration factor (R) of 2 and ˜80% with R = 6. In the human fMRI data acquired with 1.5-mm spatial resolution and 800-ms volume TR (TRvol), k-t 3D-rFOVI with R = 4 detected 46% more activated voxels in the visual cortex than the multiband EPI. Additional fMRI experiments showed that k-t 3D-rFOVI can achieve TRvol of 480 ms with R = 6, while reliably detecting visual activation. k-t 3D-rFOVI can simultaneously achieve a high spatial resolution (1.5-mm isotropically) and short TRvol (480-ms) at 3T. It offers a robust acquisition technique for fast fMRI studies over a focused brain volume.

Association of body-mass index with physiological brain pulsations across adulthood – a fast fMRI study

Background/objectiveObesity is a risk factor for several brain-related health issues, and high body-mass index (BMI) is associated with an increased risk for several neurological conditions, including cognitive decline and dementia. Cardiovascular, respiratory, and vasomotor brain pulsations have each been shown to drive intracranial cerebrovascular fluid (CSF) flow, which is linked to the brain metabolite efflux that sustains homeostasis. While these three physiological pulsations are demonstrably altered in numerous brain diseases, there is no previous investigation of the association between physiological brain pulsations and BMI.Subjects/methodsWe measured the amplitudes of the physiological brain pulsations using amplitude of low frequency fluctation (ALFF) based method with resting-state functional magnetic resonance imaging via high temporal resolution whole-brain magnetic resonance encephalography (MREG) in 115 healthy subjects. We next undertook multiple linear regression to model the BMI effect voxel-wise whole-brain on very low frequency (VLF), respiration, cardiovascular, and respiratory induced modulation of cardiovascular pulsation amplitudes with age, pulse pressure, and gender as nuisance variables.ResultsIn our study population, BMI was positively associated with the amplitudes of vasomotor, respiratory, and respiratory induced modulations of cardiovascular pulsations (p < 0.05), while negatively associated with the amplitudes of cardiovascular pulsations (p < 0.05).ConclusionsThe findings suggest that BMI is a significant factor in alterations of cardiovascular pulsation of neurofluids. As physiological pulsations are the drivers of CSF flow and subsequent metabolite clearance, these results emphasize the need for further research into the mechanisms through which obesity affects brain clearance.

Functional magnetic resonance imaging signal has sub-second temporal accuracy

Neuronal activation sequence information is essential for understanding brain functions. Extracting such timing information from blood-oxygenation-level-dependent functional magnetic resonance imaging (fMRI) signals is confounded by local cerebral vascular reactivity (CVR), which varies across brain locations. Thus, detecting neuronal synchrony as well as inferring inter-regional causal modulation using fMRI signals can be biased. Here we used fast fMRI measurements sampled at 10 Hz to measure the fMRI latency difference between visual and sensorimotor areas when participants engaged in a visuomotor task. The regional fMRI timing was calibrated by subtracting the CVR latency measured by a breath-holding task. After CVR calibration, the fMRI signal at the lateral geniculate nucleus (LGN) preceded that at the visual cortex by 496 ms, followed by the fMRI signal at the sensorimotor cortex with a latency of 464 ms. Sequential LGN, visual, and sensorimotor cortex activations were found in each participant after the CVR calibration. These inter-regional fMRI timing differences across and within participants were more closely related to the reaction time after the CVR calibration. Our results suggested the feasibility of mapping brain activity using fMRI with accuracy in hundreds of milliseconds.

BOLD cardiorespiratory pulsatility in the brain: from noise to signal of interest.

Functional magnetic resonance imaging (fMRI) based on the Blood Oxygen Level Dependent (BOLD) contrast has been extensively used to map brain activity and connectivity in health and disease. Standard fMRI preprocessing includes different steps to remove confounds unrelated to neuronal activity. First, this narrative review explores how signal fluctuations due to cardiac and respiratory activity, usually considered as "physiological noise" and regressed out from fMRI time series. However, these signal components bear useful information about some mechanisms of brain functioning (e.g., glymphatic clearance) or cerebrovascular compliance in response to arterial pressure waves. Aging and chronic diseases can cause stiffening of the aorta and other main arteries, with a reduced dampening effect resulting in greater transmission of pressure impulses to the brain. Importantly, the continuous hammering of cardiac pulsations can produce local alterations of the mechanical properties of the small cerebral vessels, with a progressive deterioration that ultimately affects neuronal functionality. Second, the review emphasizes how fMRI can study the brain patterns most affected by cardiac pulsations in health and disease with high spatiotemporal resolution, offering the opportunity to identify much more specific risk markers than systemic factors based on measurements of the vascular compliance of large arteries or other global risk factors. In this regard, modern fast fMRI acquisition techniques allow a better characterization of these pulsatile signal components due to reduced aliasing effects, turning what has been traditionally considered as noise in a signal of interest that can be used to develop novel non-invasive biomarkers in different clinical contexts.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing.

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, because differences in fMRI frequency content can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

BOLD Response is more than just magnitude: Improving detection sensitivity through capturing hemodynamic profiles

Typical fMRI analyses often assume a canonical hemodynamic response function (HRF) that primarily focuses on the peak height of the overshoot, neglecting other morphological aspects. Consequently, reported analyses often reduce the overall response curve to a single scalar value. In this study, we take a data-driven approach to HRF estimation at the whole-brain voxel level, without assuming a response profile at the individual level. We then employ a roughness penalty at the population level to estimate the response curve, aiming to enhance predictive accuracy, inferential efficiency, and cross-study reproducibility. By examining a fast event-related FMRI dataset, we demonstrate the shortcomings and information loss associated with adopting the canonical approach. Furthermore, we address the following key questions: 1) To what extent does the HRF shape vary across different regions, conditions, and participant groups? 2) Does the data-driven approach improve detection sensitivity compared to the canonical approach? 3) Can analyzing the HRF shape help validate the presence of an effect in conjunction with statistical evidence? 4) Does analyzing the HRF shape offer evidence for whole-brain response during a simple task?

Brain state transition analysis using ultra-fast fMRI differentiates MCI from cognitively normal controls.

PurposeConventional resting-state fMRI studies indicate that many cortical and subcortical regions have altered function in Alzheimer’s disease (AD) but the nature of this alteration has remained unclear. Ultrafast fMRIs with sub-second acquisition times have the potential to improve signal contrast and enable advanced analyses to understand temporal interactions between brain regions as opposed to spatial interactions. In this work, we leverage such fast fMRI acquisitions from Alzheimer’s disease Neuroimaging Initiative to understand temporal differences in the interactions between resting-state networks in 55 older adults with mild cognitive impairment (MCI) and 50 cognitively normal healthy controls.MethodsWe used a sliding window approach followed by k-means clustering. At each window, we computed connectivity i.e., correlations within and across the regions of the default mode, salience, dorsal attention, and frontoparietal network. Visual and somatosensory networks were excluded due to their lack of association with AD. Using the Davies–Bouldin index, we identified clusters of windows with distinct connectivity patterns, also referred to as brain states. The fMRI time courses were converted into time courses depicting brain state transition. From these state time course, we calculated the dwell time for each state i.e., how long a participant spent in each state. We determined how likely a participant transitioned between brain states. Both metrics were compared between MCI participants and controls using a false discovery rate correction of multiple comparisons at a threshold of. 0.05.ResultsWe identified 8 distinct brain states representing connectivity within and between the resting state networks. We identified three transitions that were different between controls and MCI, all involving transitions in connectivity between frontoparietal, dorsal attention, and default mode networks (p<0.04).ConclusionWe show that ultra-fast fMRI paired with dynamic functional connectivity analysis allows us to capture temporal transitions between brain states. Most changes were associated with transitions between the frontoparietal and dorsal attention networks connectivity and their interaction with the default mode network. Although future work needs to validate these findings, the brain networks identified in our work are known to interact with each other and play an important role in cognitive function and memory impairment in AD.

A temporal sequence of thalamic activity unfolds at transitions in behavioral arousal state

Awakening from sleep reflects a profound transformation in neural activity and behavior. The thalamus is a key controller of arousal state, but whether its diverse nuclei exhibit coordinated or distinct activity at transitions in behavioral arousal state is unknown. Using fast fMRI at ultra-high field (7 Tesla), we measured sub-second activity across thalamocortical networks and within nine thalamic nuclei to delineate these dynamics during spontaneous transitions in behavioral arousal state. We discovered a stereotyped sequence of activity across thalamic nuclei and cingulate cortex that preceded behavioral arousal after a period of inactivity, followed by widespread deactivation. These thalamic dynamics were linked to whether participants subsequently fell back into unresponsiveness, with unified thalamic activation reflecting maintenance of behavior. These results provide an outline of the complex interactions across thalamocortical circuits that orchestrate behavioral arousal state transitions, and additionally, demonstrate that fast fMRI can resolve sub-second subcortical dynamics in the human brain.

Brain respiratory pulsatility of fast fMRI stabilizes during NREM Stage 2 sleep

Brain respiratory pulsatility of fast fMRI stabilizes during NREM Stage 2 sleep

Multi-echo balanced SSFP with a sequential phase-encoding order for functional MR imaging at 7T.

To develop a 2D multi-echo passband balanced SSFP (bSSFP) sequence using an echo-train readout with a sequential phase-encoding order (sequential multi-echo bSSFP), and evaluate its performance in fast functional brain imaging at 7 T. As images of sequential multi-echo bSSFP exhibit multiple ghosts due to periodic k-space modulations, a GRAPPA-based reconstruction method was proposed to eliminate ghosting artifacts. MRI experiments were performed to compare the image quality of multi-echo bSSFP and conventional single-echo bSSFP. Submillimeter-resolution fMRI using a checkerboard visual stimulus was conducted to compare the activation characteristics of multi-echo bSSFP, conventional single-echo bSSFP and standard gradient-echo EPI (GE-EPI). A higher mean structural similarity index was found between images of single-echo bSSFP and multi-echo bSSFP with a shorter echo train length (ETL). Multi-echo bSSFP (ETL=3) showed higher temporal SNR (tSNR) values than GRAPPA-accelerated single-echo bSSFP (R=2). In submillimeter-resolution fMRI experiments, multi-echo bSSFP (ETL=3) approached the imaging speed of GRAPPA-accelerated single-echo bSSFP (R=2), but without tSNR penalty and reduced activation due to acceleration. The median t-value and the number of significantly activated voxels were comparable between GE-EPI and multi-echo bSSFP (ETL=3) that provides virtually distortion-free functional images and inherits the activation patterns of conventional bSSFP. Sequential multi-echo bSSFP (ETL=3) is suitable for fast fMRI with submillimeter in-plane resolution, and offers an option to accelerate bSSFP imaging without tSNR penalty like parallel imaging.

Human NREM Sleep Promotes Brain-Wide Vasomotor and Respiratory Pulsations.

The physiological underpinnings of the necessity of sleep remain uncertain. Recent evidence suggests that sleep increases the convection of cerebrospinal fluid (CSF) and promotes the export of interstitial solutes, thus providing a framework to explain why all vertebrate species require sleep. Cardiovascular, respiratory and vasomotor brain pulsations have each been shown to drive CSF flow along perivascular spaces, yet it is unknown how such pulsations may change during sleep in humans. To investigate these pulsation phenomena in relation to sleep, we simultaneously recorded fast fMRI, magnetic resonance encephalography (MREG), and electroencephalography (EEG) signals in a group of healthy volunteers. We quantified sleep-related changes in the signal frequency distributions by spectral entropy analysis and calculated the strength of the physiological (vasomotor, respiratory, and cardiac) brain pulsations by power sum analysis in 15 subjects (age 26.5 ± 4.2 years, 6 females). Finally, we identified spatial similarities between EEG slow oscillation (0.2–2 Hz) power and MREG pulsations. Compared with wakefulness, nonrapid eye movement (NREM) sleep was characterized by reduced spectral entropy and increased brain pulsation intensity. These effects were most pronounced in posterior brain areas for very low-frequency (≤0.1 Hz) vasomotor pulsations but were also evident brain-wide for respiratory pulsations, and to a lesser extent for cardiac brain pulsations. There was increased EEG slow oscillation power in brain regions spatially overlapping with those showing sleep-related MREG pulsation changes. We suggest that reduced spectral entropy and enhanced pulsation intensity are characteristic of NREM sleep. With our findings of increased power of slow oscillation, the present results support the proposition that sleep promotes fluid transport in human brain.SIGNIFICANCE STATEMENT We report that the spectral power of physiological brain pulsation mechanisms driven by vasomotor, respiration, and cardiac rhythms in human brain increase during sleep, extending previous observations of their association with glymphatic brain clearance during sleep in rodents. The magnitudes of increased pulsations follow the rank order of vasomotor greater than respiratory greater than cardiac pulsations, with correspondingly declining spatial extents. Spectral entropy, previously known as vigilance and as an anesthesia metric, decreased during NREM sleep compared with the awake state in very low and respiratory frequencies, indicating reduced signal complexity. An EEG slow oscillation power increase occurring in the early sleep phase (NREM 1–2) spatially overlapped with pulsation changes, indicating reciprocal mechanisms between those measures.

An efficient functional magnetic resonance imaging data reduction strategy using neighborhood preserving embedding algorithm.

High dimensionality data have become common in neuroimaging fields, especially group‐level functional magnetic resonance imaging (fMRI) datasets. fMRI connectivity analysis is a widely used, powerful technique for studying functional brain networks to probe underlying mechanisms of brain function and neuropsychological disorders. However, data‐driven technique like independent components analysis (ICA), can yield unstable and inconsistent results, confounding the true effects of interest and hindering the understanding of brain functionality and connectivity. A key contributing factor to this instability is the information loss that occurs during fMRI data reduction. Data reduction of high dimensionality fMRI data in the temporal domain to identify the important information within group datasets is necessary for such analyses and is crucial to ensure the accuracy and stability of the outputs. In this study, we describe an fMRI data reduction strategy based on an adapted neighborhood preserving embedding (NPE) algorithm. Both simulated and real data results indicate that, compared with the widely used data reduction method, principal component analysis, the NPE‐based data reduction method (a) shows superior performance on efficient data reduction, while enhancing group‐level information, (b) develops a unique stratagem for selecting components based on an adjacency graph of eigenvectors, (c) generates more reliable and reproducible brain networks under different model orders when the outputs of NPE are used for ICA, (d) is more sensitive to revealing task‐evoked activation for task fMRI, and (e) is extremely attractive and powerful for the increasingly popular fast fMRI and very large datasets.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response.

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

Statistical power or more precise insights into neuro-temporal dynamics? Assessing the benefits of rapid temporal sampling in fMRI.

Functional magnetic resonance imaging (fMRI), a non-invasive and widely used human neuroimaging method, is most known for its spatial precision. However, there is a growing interest in its temporal sensitivity. This is despite the temporal blurring of neuronal events by the blood oxygen level dependent (BOLD) signal, the peak of which lags neuronal firing by 4-6 seconds. Given this, the goal of this review is to answer a seemingly simple question - "What are the benefits of increased temporal sampling for fMRI?". To answer this, we have combined fMRI data collected at multiple temporal scales, from 323 to 1000 milliseconds, with a review of both historical and contemporary temporal literature. After a brief discussion of technological developments that have rekindled interest in temporal research, we next consider the potential statistical and methodological benefits. Most importantly, we explore how fast fMRI can uncover previously unobserved neuro-temporal dynamics - effects that are entirely missed when sampling at conventional 1 to 2 second rates. With the intrinsic link between space and time in fMRI, this temporal renaissance also delivers improvements in spatial precision. Far from producing only statistical gains, the array of benefits suggest that the continued temporal work is worth the effort.

Investigating mechanisms of fast BOLD responses: The effects of stimulus intensity and of spatial heterogeneity of hemodynamics

Recent studies have demonstrated that fast fMRI can track neural activity well above the temporal limit predicted by the canonical hemodynamic response model. While these findings are promising, the biophysical mechanisms underlying these fast fMRI phenomena remain underexplored. In this study, we discuss two aspects of the hemodynamic response, complementary to several existing hypotheses, that can accommodate faster fMRI dynamics beyond those predicted by the canonical model. First, we demonstrate, using both visual and somatosensory paradigms, that the timing and shape of hemodynamic response functions (HRFs) vary across graded levels of stimulus intensity—with lower-intensity stimulation eliciting faster and narrower HRFs. Second, we show that as the spatial resolution of fMRI increases, voxel-wise HRFs begin to deviate from the canonical model, with a considerable portion of voxels exhibiting faster temporal dynamics than predicted by the canonical HRF. Collectively, both stimulus/task intensity and image resolution can affect the sensitivity of fMRI to fast brain activity, which may partly explain recent observations of fast fMRI signals. It is further noteworthy that, while the present investigations focus on fast neural responses, our findings suggest that a revised hemodynamic model may benefit the many fMRI studies using paradigms with wide ranges of contrast levels (e.g., resting or naturalistic conditions) or with modern, high-resolution MR acquisitions.

Functional MRI of Native and Non-native Speech Sound Production in Sequential German-English Bilinguals.

Bilingualism and multilingualism are highly prevalent. Non-invasive brain imaging has been used to study the neural correlates of native and non-native speech and language production, mainly on the lexical and syntactic level. Here, we acquired continuous fast event-related FMRI during visually cued overt production of exclusively German and English vowels and syllables. We analyzed data from 13 university students, native speakers of German and sequential English bilinguals. The production of non-native English sounds was associated with increased activity of the left primary sensorimotor cortex, bilateral cerebellar hemispheres (lobule VI), left inferior frontal gyrus, and left anterior insula compared to native German sounds. The contrast German > English sounds was not statistically significant. Our results emphasize that the production of non-native speech requires additional neural resources already on a basic phonological level in sequential bilinguals.

Frequency-specific task modulation of human brain functional networks: A fast fMRI study

How coherent neural oscillations are involved in task execution is a fundamental question in neuroscience. Although several electrophysiological studies have tackled this issue, the brain-wide task modulation of neural coherence remains uncharacterized. Here, with a fast fMRI technique, we studied shifts of brain-wide neural coherence across different task states in the ultraslow frequency range (0.01–0.7 Hz). First, we examined whether the shifts of the brain-wide neural coherence occur in a frequency-dependent manner. We quantified the shift of a region's average neural coherence by the inter-state variance of the mean coherence between the region and the rest of the brain. A clustering analysis based on the variance's spatial correlation between frequency components revealed four frequency bands (0.01–0.15 Hz, 0.15–0.37 Hz, 0.37–0.53 Hz, and 0.53–0.7 Hz) showing band-specific shifts of the brain-wide neural coherence. Next, we investigated the similarity of the inter-state variance's spectra between all pairs of regions. We found that regions showing similar spectra correspond to those forming functional modules of the brain network. Then, we investigated the relationship between identified frequency bands and modules’ inter-state variances. We found that modules showing the highest variance are those made up of parieto-occipital regions at 0.01–0.15 Hz, while it is replaced with another consisting of frontal regions above 0.15 Hz. Furthermore, these modules showed specific shifting patterns of the mean coherence across states at 0.01–0.15 Hz and above 0.15 Hz, suggesting that identified frequency bands differentially contribute to neural interactions during task execution. Our results highlight that usage of the fast fMRI enables brain-wide investigation of neural coherence up to 0.7 Hz, which opens a promising track for assessment of the large-scale neural interactions in the ultraslow frequency range.

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.

Concurrent electrophysiological and hemodynamic measurements of evoked neural oscillations in human visual cortex using sparsely interleaved fast fMRI and EEG

Electroencephalography (EEG) concurrently collected with functional magnetic resonance imaging (fMRI) is heavily distorted by the repetitive gradient coil switching during the fMRI acquisition. The performance of the typical template-based gradient artifact suppression method can be suboptimal because the artifact changes over time. Gradient artifact residuals also impede the subsequent suppression of ballistocardiography artifacts.Here we propose recording continuous EEG with temporally sparse fast fMRI (fast fMRI-EEG) to minimize the EEG artifacts caused by MRI gradient coil switching without significantly compromising the field-of-view and spatiotemporal resolution of fMRI. Using simultaneous multi-slice inverse imaging to achieve whole-brain fMRI with isotropic 5-mm resolution in 0.1 ​s, and performing these acquisitions once every 2 ​s, we have 95% of the duty cycle available to record EEG with substantially less gradient artifact. We found that the standard deviation of EEG signals over the entire acquisition period in fast fMRI-EEG was reduced to 54% of that in conventional concurrent echo-planar imaging (EPI) and EEG recordings (EPI-EEG) across participants. When measuring 15-Hz steady-state visual evoked potentials (SSVEPs), the baseline-normalized oscillatory neural response in fast fMRI-EEG was 2.5-fold of that in EPI-EEG. The functional MRI responses associated with the SSVEP delineated by EPI and fast fMRI were similar in the spatial distribution, the elicited waveform, and detection power. Sparsely interleaved fast fMRI-EEG provides high-quality EEG without substantially compromising the quality of fMRI in evoked response measurements, and has the potential utility for applications where the onset of the target stimulus cannot be precisely determined, such as epilepsy.