Scanning fast and slow: current limitations of 3 Tesla functional MRI and future potential

Alterations in Functional Network Topology Within Normal Hemispheres Contralateral to Anterior Circulation Steno-Occlusive Disease: A Resting-State BOLD Study

The purpose of this study was to assess spatially remote effects of hemodynamic impairment on functional network topology contralateral to unilateral anterior circulation steno-occlusive disease (SOD) using resting-state blood oxygen level-dependent (BOLD) imaging, and to investigate the relationships between network connectivity and cerebrovascular reactivity (CVR), a measure of hemodynamic stress. Twenty patients with unilateral, chronic anterior circulation SOD and 20 age-matched healthy controls underwent resting-state BOLD imaging. Five-minute standardized baseline BOLD acquisition was followed by acetazolamide infusion to measure CVR. The BOLD baseline was used to analyze network connectivity contralateral to the diseased hemispheres of SOD patients. Compared to healthy controls, reduced network degree (z-score = −1.158 ± 1.217, P < 0.001, false discovery rate (FDR) corrected), local efficiency (z-score = −1.213 ± 1.120, P < 0.001, FDR corrected), global efficiency (z-score = −1.346 ± 1.119, P < 0.001, FDR corrected), and enhanced modularity (z-score = 1.000 ± 1.205, P = 0.002, FDR corrected) were observed in the contralateral, normal hemispheres of SOD patients. Network degree (P = 0.089, FDR corrected; P = 0.027, uncorrected) and nodal efficiency (P = 0.089, FDR corrected; P = 0.045, uncorrected) showed a trend toward a positive association with CVR. The results indicate remote abnormalities in functional connectivity contralateral to the diseased hemispheres in patients with unilateral SOD, despite the absence of macrovascular disease or demonstrable hemodynamic impairment. The clinical impact of remote functional disruptions requires dedicated investigation but may portend far reaching consequence for even putatively unilateral cerebrovascular disease.

Association between tissue hypoxia, perfusion restrictions, and microvascular architecture alterations with lesion-induced impairment of neurovascular coupling

Functional magnetic resonance imaging (fMRI) has been mainly utilized for the preoperative localization of eloquent cortical areas. However, lesion-induced impairment of neurovascular coupling (NVC) in the lesion border zone may lead to false-negative fMRI results. The purpose of this study was to determine physiological factors impacting the NVC. Twenty patients suffering from brain lesions were preoperatively examined using multimodal neuroimaging including fMRI, magnetoencephalography (MEG) during language or sensorimotor tasks (depending on lesion location), and a novel physiologic MRI approach for the combined quantification of oxygen metabolism, perfusion state, and microvascular architecture. Congruence of brain activity patterns between fMRI and MEG were found in 13 patients. In contrast, we observed missing fMRI activity in perilesional cortex that demonstrated MEG activity in seven patients, which was interpreted as lesion-induced impairment of NVC. In these brain regions with impaired NVC, physiologic MRI revealed significant brain tissue hypoxia, as well as significantly decreased macro- and microvascular perfusion and microvascular architecture. We demonstrated that perilesional hypoxia with reduced vascular perfusion and architecture is associated with lesion-induced impairment of NVC. Our physiologic MRI approach is a clinically applicable method for preoperative risk assessment for the presence of false-negative fMRI results and may prevent severe postoperative functional deficits.

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.

Discriminating stress from rest based on resting-state connectivity of the human brain: A supervised machine learning study

Acute stress induces large-scale neural reorganization with relevance to stress-related psychopathology. Here, we applied a novel supervised machine learning method, combining the strengths of a priori theoretical insights with a data-driven approach, to identify which connectivity changes are most prominently associated with a state of acute stress and individual differences therein. Resting-state functional magnetic resonance imaging scans were taken from 334 healthy participants (79 females) before and after a formal stress induction. For each individual scan, mean time-series were extracted from 46 functional parcels of three major brain networks previously shown to be potentially sensitive to stress effects (default mode network (DMN), salience network (SN), and executive control networks). A data-driven approach was then used to obtain discriminative spatial linear filters that classified the pre- and post-stress scans. To assess potential relevance for understanding individual differences, probability of classification using the most discriminative filters was linked to individual cortisol stress responses. Our model correctly classified pre- versus post-stress states with highly significant accuracy (above 75%; leave-one-out validation relative to chance performance). Discrimination between pre- and post-stress states was mainly based on connectivity changes in regions from the SN and DMN, including the dorsal anterior cingulate cortex, amygdala, posterior cingulate cortex, and precuneus. Interestingly, the probability of classification using these connectivity changes were associated with individual cortisol increases. Our results confirm the involvement of DMN and SN using a data-driven approach, and specifically single out key regions that might receive additional attention in future studies for their relevance also for individual differences.

Does higher sampling rate (multiband + SENSE) improve group statistics - An example from social neuroscience block design at 3T

Multiband (MB) or Simultaneous multi-slice (SMS) acquisition schemes allow the acquisition of MRI signals from more than one spatial coordinate at a time. Commercial availability has brought this technique within the reach of many neuroscientists and psychologists. Most early evaluation of the performance of MB acquisition employed resting state fMRI or the most basic tasks. In this study, we tested whether the advantages of using MB acquisition schemes generalize to group analyses using a cognitive task more representative of typical cognitive neuroscience applications. Twenty-three subjects were scanned on a Philips 3 ​T scanner using five sequences, up to eight-fold acceleration with MB-factors 1 to 4, SENSE factors up to 2 and corresponding TRs of 2.45s down to 0.63s, while they viewed (i) movie blocks showing complex actions with hand object interactions and (ii) control movie blocks without hand object interaction. Data were processed using a widely used analysis pipeline implemented in SPM12 including the unified segmentation and canonical HRF modelling. Using random effects group-level, voxel-wise analysis we found that all sequences were able to detect the basic action observation network known to be recruited by our task. The highest t-values were found for sequences with MB4 acceleration. For the MB1 sequence, a 50% bigger voxel volume was needed to reach comparable t-statistics. The group-level t-values for resting state networks (RSNs) were also highest for MB4 sequences. Here the MB1 sequence with larger voxel size did not perform comparable to the MB4 sequence. Altogether, we can thus recommend the use of MB4 (and SENSE 1.5 or 2) on a Philips scanner when aiming to perform group-level analyses using cognitive block design fMRI tasks and voxel sizes in the range of cortical thickness (e.g. 2.7 ​mm isotropic). While results will not be dramatically changed by the use of multiband, our results suggest that MB will bring a moderate but significant benefit.

Amalgamating evidence of dynamics

Many approaches to evidence amalgamation focus on relatively static information or evidence: the data to be amalgamated involve different variables, contexts, or experiments, but not measurements over extended periods of time. However, much of scientific inquiry focuses on dynamical systems; the system's behavior over time is critical. Moreover, novel problems of evidence amalgamation arise in these contexts. First, data can be collected at different measurement timescales, where potentially none of them correspond to the underlying system's causal timescale. Second, missing variables have a significantly different impact on time series measurements than they do in the traditional static setting; in particular, they make causal and structural inference much more difficult. In this paper, we argue that amalgamation should proceed by integrating causal knowledge, rather than at the level of "raw" evidence. We defend this claim by first outlining both of these problems, and then showing that they can be solved only if we operate on causal structures. We therefore must use causal discovery methods that are reliable given these problems. Such methods do exist, but their successful application requires careful consideration of the problems that we highlight.

Acute stress alters the 'default' brain processing

Active adaptation to acute stress is essential for coping with daily life challenges. The stress hormone cortisol, as well as large scale re-allocations of brain resources have been implicated in this adaptation. Stress-induced shifts between large-scale brain networks, including salience (SN), central executive (CEN) and default mode networks (DMN), have however been demonstrated mainly under task-conditions. It remains unclear whether such network shifts also occur in the absence of ongoing task-demands, and most critically, whether these network shifts are predictive of individual variation in the magnitude of cortisol stress-responses. In a sample of 335 healthy participants, we investigated stress-induced functional connectivity changes (delta-FC) of the SN, CEN and DMN, using resting-state fMRI data acquired before and after a socially evaluated cold-pressor test and a mental arithmetic task. To investigate which network changes are associated with acute stress, we evaluated the association between cortisol increase and delta-FC of each network. Stress-induced cortisol increase was associated with increased connectivity within the SN, but with decreased coupling of DMN at both local (within network) and global (synchronization with brain regions also outside the network) levels. These findings indicate that acute stress prompts immediate connectivity changes in large-scale resting-state networks, including the SN and DMN in the absence of explicit ongoing task-demands. Most interestingly, this brain reorganization is coupled with individuals' cortisol stress-responsiveness. These results suggest that the observed stress-induced network reorganization might function as a neural mechanism determining individual stress reactivity and, therefore, it could serve as a promising marker for future studies on stress resilience and vulnerability.

Functional neuroimaging of the central autonomic network: recent developments and clinical implications

Purpose

The central autonomic network (CAN) is an intricate system of brainstem, subcortical, and cortical structures that play key roles in the function of the autonomic nervous system. Prior to the advent of functional neuroimaging, in vivo studies of the human CAN were limited. The purpose of this review is to highlight the contribution of functional neuroimaging, specifically functional magnetic resonance imaging (fMRI), to the study of the CAN, and to discuss recent advances in this area. Additionally, we aim to emphasize exciting areas for future research.

Methods

We reviewed the existing literature in functional neuroimaging of the CAN. Here, we focus on fMRI research conducted in healthy human subjects, as well as research that has been done in disease states, to understand CAN function. To minimize confounding, papers examining CAN function in the context of cognition, emotion, pain, and affective disorders were excluded.

Results

fMRI has led to significant advances in the understanding of human CAN function. The CAN is composed of widespread brainstem and forebrain structures that are intricately connected and play key roles in reflexive and modulatory control of autonomic function.

Conclusions

fMRI technology has contributed extensively to current knowledge of CAN function. It holds promise to serve as a biomarker in disease states. With ongoing advancements in fMRI technology, there is great opportunity and need for future research involving the CAN.

Flexible 23-channel coil array for high-resolution magnetic resonance imaging at 3 Tesla

PURPOSE

The purpose of this work is the design, implementation and evaluation of a mechanically flexible receive-only coil array for magnetic resonance imaging (MRI) at 3 T that can be applied to various target organs and provides high parallel imaging performance.

METHODS

A 23-channel array was designed based on a rigid-flex printed circuit board (PCB). The flexible multi-layer part contains the copper traces forming the coil elements. The rigid part of the PCB houses the solder joints and lumped elements. The coil housing consists of rigid caps mounted above the rigid parts. Adhesive PTFE sheets cover all flexible parts. The developed array was tested on the bench as well as in phantom and in vivo MRI experiments employing parallel imaging acceleration factors up to six.

RESULTS

Efficient mutual decoupling between receive elements and detuning between receive array and body coil was achieved. An increased signal-to-noise ratio in comparison to commercial reference coils is demonstrated, especially in regions close to the developed array and for high parallel imaging acceleration factors. Exemplary in vivo images of head, ankle, knee, shoulder and hand are presented.

CONCLUSION

Based on high sensitivity close to the array and low g-factors, this flexible coil is well suited for studies of occipital and temporal cortex, as well as musculoskeletal targets like knee, ankle, elbow and wrist.

Simultaneous perfusion and permeability assessments using multiband multi-echo EPI (M2-EPI) in brain tumors

PURPOSE

To study a multiband multi-echo EPI (M2-EPI) sequence for dynamic susceptibility contrast (DSC) perfusion imaging with leakage correction and vascular permeability measurements, and to evaluate the benefits of increased temporal resolution provided by this acquisition strategy on the accuracy of perfusion and permeability estimations.

METHODS

A novel M2-EPI sequence was developed, and a pharmacokinetic model accounting for contrast agent extravasation was used to produce perfusion maps and additional vascular permeability maps. The advantage of M2-EPI for DSC perfusion imaging was demonstrated in vivo in 5 patients with brain tumors, and numerical simulations were performed to evaluate the advantage of improved temporal resolution afforded by the technique.

RESULTS

In contrast to underestimations of cerebral blood volume (CBV) in tumors using the single-echo acquisition strategy, M2-EPI provided more plausible estimates of CBV. A quantitative evaluation showed higher estimated values of CBV and mean transit time in tumor tissues using M2-EPI (CBV: 3.08 ± 0.78 mL/100 g versus 1.56 ± 1.38 mL/100 g [P = .006]; mean transit time: 4.94 ± 1.17 seconds versus 1.83 ± 2.06 seconds [P = 0.033]). Numerical simulations showed that higher temporal resolution provided by M2-EPI was associated with more accurate estimates of cerebral blood flow, CBV, and permeability parameters.

CONCLUSION

The novel M2-EPI acquisition strategy for DSC imaging facilitates leakage-corrected perfusion measurements with additional permeability assessments and more accurate estimates of perfusion/permeability parameters, and may be used as a quantitative tool for the diagnosis, prognosis, and treatment monitoring of brain tumors.

Spectral Diversity in Default Mode Network Connectivity Reflects Behavioral State

Default mode network (DMN) functional connectivity is thought to occur primarily in low frequencies (<0.1 Hz), resulting in most studies removing high frequencies during data preprocessing. In contrast, subtractive task analyses include high frequencies, as these are thought to be task relevant. An emerging line of research explores resting fMRI data at higher-frequency bands, examining the possibility that functional connectivity is a multiband phenomenon. Furthermore, recent studies suggest DMN involvement in cognitive processing; however, without a systematic investigation of DMN connectivity during tasks, its functional contribution to cognition cannot be fully understood. We bridged these concurrent lines of research by examining the contribution of high frequencies in the relationship between DMN and dorsal attention network at rest and during task execution. Our findings revealed that the inclusion of high frequencies alters between network connectivity, resulting in reduced anticorrelation and increased positive connectivity between DMN and dorsal attention network. Critically, increased positive connectivity was observed only during tasks, suggesting an important role for high-frequency fluctuations in functional integration. Moreover, within-DMN connectivity during task execution correlated with RT only when high frequencies were included. These results show that DMN does not simply deactivate during task execution and suggest active recruitment while performing cognitively demanding paradigms.

fMRI measurements of amygdala activation are confounded by stimulus correlated signal fluctuation in nearby veins draining distant brain regions

Imaging the amygdala with functional MRI is confounded by multiple averse factors, notably signal dropouts due to magnetic inhomogeneity and low signal-to-noise ratio, making it difficult to obtain consistent activation patterns in this region. However, even when consistent signal changes are identified, they are likely to be due to nearby vessels, most notably the basal vein of rosenthal (BVR). Using an accelerated fMRI sequence with a high temporal resolution (TR = 333 ms) combined with susceptibility-weighted imaging, we show how signal changes in the amygdala region can be related to a venous origin. This finding is confirmed here in both a conventional fMRI dataset (TR = 2000 ms) as well as in information of meta-analyses, implying that “amygdala activations” reported in typical fMRI studies are likely confounded by signals originating in the BVR rather than in the amygdala itself, thus raising concerns about many conclusions on the functioning of the amygdala that rely on fMRI evidence alone.

On the generalizability of resting-state fMRI machine learning classifiers

Machine learning classifiers have become increasingly popular tools to generate single-subject inferences from fMRI data. With this transition from the traditional group level difference investigations to single-subject inference, the application of machine learning methods can be seen as a considerable step forward. Existing studies, however, have given scarce or no information on the generalizability to other subject samples, limiting the use of such published classifiers in other research projects. We conducted a simulation study using publicly available resting-state fMRI data from the 1000 Functional Connectomes and COBRE projects to examine the generalizability of classifiers based on regional homogeneity of resting-state time series. While classification accuracies of up to 0.8 (using sex as the target variable) could be achieved on test datasets drawn from the same study as the training dataset, the generalizability of classifiers to different study samples proved to be limited albeit above chance. This shows that on the one hand a certain amount of generalizability can robustly be expected, but on the other hand this generalizability should not be overestimated. Indeed, this study substantiates the need to include data from several sites in a study investigating machine learning classifiers with the aim of generalizability.

The spectral diversity of resting-state fluctuations in the human brain

In order to assess whole-brain resting-state fluctuations at a wide range of frequencies, resting-state fMRI data of 20 healthy subjects were acquired using a multiband EPI sequence with a low TR (354 ms) and compared to 20 resting-state datasets from standard, high-TR (1800 ms) EPI scans. The spatial distribution of fluctuations in various frequency ranges are analyzed along with the spectra of the time-series in voxels from different regions of interest. Functional connectivity specific to different frequency ranges (<0.1 Hz; 0.1–0.25 Hz; 0.25–0.75 Hz; 0.75–1.4 Hz) was computed for both the low-TR and (for the two lower-frequency ranges) the high-TR datasets using bandpass filters. In the low-TR data, cortical regions exhibited highest contribution of low-frequency fluctuations and the most marked low-frequency peak in the spectrum, while the time courses in subcortical grey matter regions as well as the insula were strongly contaminated by high-frequency signals. White matter and CSF regions had highest contribution of high-frequency fluctuations and a mostly flat power spectrum. In the high-TR data, the basic patterns of the low-TR data can be recognized, but the high-frequency proportions of the signal fluctuations are folded into the low frequency range, thus obfuscating the low-frequency dynamics. Regions with higher proportion of high-frequency oscillations in the low-TR data showed flatter power spectra in the high-TR data due to aliasing of the high-frequency signal components, leading to loss of specificity in the signal from these regions in high-TR data. Functional connectivity analyses showed that there are correlations between resting-state signal fluctuations of distant brain regions even at high frequencies, which can be measured using low-TR fMRI. On the other hand, in the high-TR data, loss of specificity of measured fluctuations leads to lower sensitivity in detecting functional connectivity. This underlines the advantages of low-TR EPI sequences for resting-state and potentially also task-related fMRI experiments.