Reducing task-based fMRI scanning time using simultaneous multislice echo planar imaging

Shorter TR combined with finer atlas positively modulate topological organization of brain network: A resting state fMRI study

BACKGROUND

The use of shorter TR and finer atlases in rs-fMRI can provide greater detail on brain function and anatomy. However, there is limited understanding of the effect of this combination on brain network properties.

METHODS

A study was conducted with 20 healthy young volunteers who underwent rs-fMRI scans with both shorter (0.5s) and long (2s) TR. Two atlases with different degrees of granularity (90 vs 200 regions) were used to extract rs-fMRI signals. Several network metrics, including small-worldness, Cp, Lp, Eloc, and Eg, were calculated. Two-factor ANOVA and two-sample t-tests were conducted for both the single spectrum and five sub-frequency bands.

RESULTS

The network constructed using the combination of shorter TR and finer atlas showed significant enhancements in Cp, Eloc, and Eg, as well as reductions in Lp and γ in both the single spectrum and subspectrum (p < 0.05, Bonferroni correction). Network properties in the 0.082-0.1 Hz frequency range were weaker than those in the 0.01-0.082 Hz range.

CONCLUSION

Our findings suggest that the use of shorter TR and finer atlas can positively affect the topological characteristics of brain networks. These insights can inform the development of brain network construction methods.

Do mean values tell the full story? Cardiac cycle and biological sex comparisons in temporally derived neurovascular coupling metrics

Previous reports have noted cerebrovascular regulation differs across the cardiac cycle, with greater regulation occurring within systole. However, this methodological notion has not been meticulously scrutinized during temporally deduced neurovascular coupling (NVC) metrics with additional respect to biological sex. Analyses of 111 healthy individuals (40 females/71 males) were performed where participants engaged in the "Where's Waldo?" paradigm. All NVC parameters were quantified in the posterior and middle cerebral arteries at 310 unique timepoints. Several individuals completed repeat testing which enabled for between-day (3 timepoints) and within-day (7 timepoints) reliability comparisons in 17 and 11 individuals, respectively. One-way analysis of variance compared NVC metrics between diastole, mean, and systole values, as well as differences between biological sexes. Greater absolute cerebral blood velocity (CBv; baseline and peak) and total activation (area under the curve) were noted within systole for both posterior cerebral artery (PCA; P < 0.001) and middle cerebral artery (MCA; P < 0.001) values; however, the relative percent increase in CBv was greater within diastole (P < 0.001). Females had an elevated diastolic and mean CBv and a greater diastolic cerebrovascular conductance (P < 0.050). No sex differences were present for systolic CBv measures and within parameters quantifying the NVC response (area under the curve/relative CBv increase) across the cardiac cycle (P > 0.072). Future investigations seeking to differentiate cerebral regulatory mechanisms between clinical populations may benefit by performing their analyses across the cardiac cycle, as certain pathogenesis may affect one aspect of the cardiac cycle independently. Minimal differences were noted between females and males for metrics characterizing the NVC response across the cardiac cycle.NEW & NOTEWORTHY Neurovascular coupling (NVC) studies commonly assess the mean cerebral hemodynamic response with little consideration for diastole, systole, and biological sex. Greater total activation expressed as the area under the curve was seen within systole compared with mean and diastole. Resting cerebral blood velocity sex differences were more prevalent during diastole when the cerebrovasculature was pressure-passive. Future studies should assess the NVC response across the cardiac cycle as it may help delineate the underlying pathophysiology of various clinical populations.

Improved Activation and Hemodynamic Response Function of Olfactory fMRI Using Simultaneous Multislice with Reduced TR Acquisition

Objectives

The respiration could decrease the time synchronization between odor stimulation and data acquisition, consequently deteriorating the functional activation and hemodynamic response function (HRF) in olfactory functional magnetic resonance imaging (fMRI) with a conventional repetition time (TR). In this study, we aimed to investigate whether simultaneous multislice (SMS) technology with reduced TR could improve the blood oxygen level-dependent (BOLD) activation and optimize HRF modeling in olfactory fMRI.

Methods

Sixteen young healthy subjects with normal olfaction underwent olfactory fMRI on a 3T MRI scanner using a 64 channel head coil. FMRI data were acquired using SMS acceleration at three different TRs: 3000 ms, 1000 ms, and 500 ms. Both metrics of BOLD activation (activated voxels, mean, and maximum t-scores) and the HRF modeling (response height and time to peak) were calculated in the bilateral amygdalae, hippocampi, and insulae.

Results

The 500 ms and 1000 ms TRs both significantly improved the number of activated voxels, mean, and maximum t-score in the amygdalae and insulae, compared with a 3000 ms TR (all P < 0.05). But the increase of these metrics in the hippocampi did not reach a statistical significance (all P > 0.05). No significant difference in any BOLD activation metrics between TR 500 ms and 1000 ms was observed in all regions of interest (ROIs) (all P > 0.05). The HRF curves showed that higher response height and shorter time to peak in all ROIs were obtained at 500 ms and 1000 ms TRs compared to 3000 ms TR. TR 500 ms had a more significant effect on response height than TR 1000 ms in the amygdalae (P = 0.017), and there was no significant difference in time to peak between TR 500 ms and 1000 ms in all ROIs (all P > 0.05).

Conclusions

The fast image acquisition technique of SMS with reduced TR could significantly improve the functional activation and HRF curve in olfactory fMRI.

Insufficient sampling frequencies skew heart rate variability estimates: Implications for extracting heart rate metrics from neuroimaging and physiological data

BACKGROUND

While cardiac pulsations are widely present within physiological and neuroimaging data, it is unknown the extent this information can provide valid and reliable heart rate and heart rate variability (HRV) estimates. The objective of this study was to demonstrate how a slight temporal shift due to an insufficient sampling frequency can impact the validity/accuracy of deriving cardiac metrics.

METHODS

Twenty-two participants were instrumented with valid/reliable industry-standard or open-source electrocardiograms. Five-minute lead II recordings were collected at 1000 Hz in an upright orthostatic position. Following artifact removal, the 1000 Hz recording for each participant was downsampled to frequencies ranging 2-500 Hz. The validity of each participant's downsampled recording was compared against their 1000 Hz recording ("reference-standard") using Bland-Altman plots with 95 % limits of agreement (LOA), coefficient of variation (CoV), intraclass correlation coefficients, and adjusted r-squared values.

RESULTS

Downsampled frequencies of ≥ 50 and ≥ 90 Hz produced highly robust measures with narrow log-transformed 95 % LOA (<±0.01) and low CoV values (≤3.5 %) for heart rate and HRV metrics, respectively. Below these thresholds, the log-transformed 95 % LOA became wider (LOA range: ±0.1-1.9) and more variable (CoV range: 1.5-111.6 %).

CONCLUSION

These results provide an important consideration for obtaining cardiac information from physiological data. Compared to the "reference-standard" ECG, a seemingly negligible temporal shift of the systolic contraction (R wave) greater than 11-milliseconds (90 Hz) away from its true value, lessened the validity of the HRV. Further research is warranted to determine the minimum sampling frequency required to obtain valid heart rate/HRV metrics from pulsatile waveforms.

An empirical investigation of the benefit of increasing the temporal resolution of task-evoked fMRI data with multi-band imaging

Objective

There is a tendency for reducing TR in MRI experiments with multi-band imaging. We empirically investigate its benefit for the group-level statistical outcome in task-evoked fMRI.

Methods

Three visual fMRI data sets were collected from 17 healthy adult participants. Multi-band acquisition helped vary the TR (2000/1000/410 ms, respectively). Because these data sets capture different temporal aspects of the haemodynamic response (HRF), we tested several HRF models. We computed a composite descriptive statistic, H, from β’s of each first-level model fit and carried it to the group-level analysis. The number of activated voxels and the t value of the group-level analysis as well as a goodness-of-fit measure were used as surrogate markers of data quality for comparison.

Results

Increasing the temporal sampling rate did not provide a universal improvement in the group-level statistical outcome. Rather, both the voxel-wise and ROI-averaged group-level results varied widely with anatomical location, choice of HRF and the setting of the TR. Correspondingly, the goodness-of-fit of HRFs became worse with increasing the sampling frequency.

Conclusion

Rather than universally increasing the temporal sampling rate in cognitive fMRI experiments, these results advocate the performance of a pilot study for the specific ROIs of interest to identify the appropriate temporal sampling rate for the acquisition and the correspondingly suitable HRF for the analysis of the data.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10334-021-00918-z.

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.

The Impact of Echo Time Shifts and Temporal Signal Fluctuations on BOLD Sensitivity in Presurgical Planning at 7 T

OBJECTIVES

Gradients in the static magnetic field caused by tissues with differing magnetic susceptibilities lead to regional variations in the effective echo time, which modifies both image signal and BOLD sensitivity. Local echo time changes are not considered in the most commonly used metric for BOLD sensitivity, temporal signal-to-noise ratio (tSNR), but may be significant, particularly at ultrahigh field close to air cavities (such as the sinuses and ear canals) and near gross brain pathologies and postoperative sites.

MATERIALS AND METHODS

We have studied the effect of local variations in echo time and tSNR on BOLD sensitivity in 3 healthy volunteers and 11 patients with tumors, postoperative cavities, and venous malformations at 7 T. Temporal signal-to-noise ratio was estimated from a 5-minute run of resting state echo planar imaging with a nominal echo time of 22 milliseconds. Maps of local echo time were derived from the phase of a multiecho GE scan. One healthy volunteer performed 10 runs of a breath-hold task. The t-map from this experiment served as a criterion standard BOLD sensitivity measure. Two runs of a less demanding breath-hold paradigm were used for patients.

RESULTS

In all subjects, a strong reduction in the echo time (from 22 milliseconds to around 11 milliseconds) was found close to the ear canals and sinuses. These regions were characterized by high tSNR but low t-values in breath-hold t-maps. In some patients, regions of particular interest in presurgical planning were affected by reductions in the echo time to approximately 13-15 milliseconds. These included the primary motor cortex, Broca's area, and auditory cortex. These regions were characterized by high tSNR values (70 and above). Breath-hold results were corrupted by strong motion artifacts in all patients.

CONCLUSIONS

Criterion standard BOLD sensitivity estimation using hypercapnic experiments is challenging, especially in patient populations. Taking into consideration the tSNR, commonly used for BOLD sensitivity estimation, but ignoring local reductions in the echo time (eg, from 22 to 11 milliseconds), would erroneously suggest functional sensitivity sufficient to map BOLD signal changes. It is therefore important to consider both local variations in the echo time and temporal variations in signal, using the product metric of these two indices for instance. This should ensure a reliable estimation of BOLD sensitivity and to facilitate the identification of potential false-negative results. This is particularly true at high fields, such as 7 T and in patients with large pathologies and postoperative cavities.