Analysis and correction of field fluctuations in fMRI data using field monitoring

The possible influence of third-order shim coils on gradient-magnet interactions: an inter-field and inter-site study

OBJECTIVE

To assess the possible influence of third-order shim coils on the behavior of the gradient field and in gradient-magnet interactions at 7 T and above.

MATERIALS AND METHODS

Gradient impulse response function measurements were performed at 5 sites spanning field strengths from 7 to 11.7 T, all of them sharing the same exact whole-body gradient coil design. Mechanical fixation and boundary conditions of the gradient coil were altered in several ways at one site to study the impact of mechanical coupling with the magnet on the field perturbations. Vibrations, power deposition in the He bath, and field dynamics were characterized at 11.7 T with the third-order shim coils connected and disconnected inside the Faraday cage.

RESULTS

For the same whole-body gradient coil design, all measurements differed greatly based on the third-order shim coil configuration (connected or not). Vibrations and gradient transfer function peaks could be affected by a factor of 2 or more, depending on the resonances. Disconnecting the third-order shim coils at 11.7 T also suppressed almost completely power deposition peaks at some frequencies.

DISCUSSION

Third-order shim coil configurations can have major impact in gradient-magnet interactions with consequences on potential hardware damage, magnet heating, and image quality going beyond EPI acquisitions.

Impact of B 0 $$ {\mathrm{B}}_0 $$ field imperfections correction on BOLD sensitivity in 3D-SPARKLING fMRI data

PURPOSE

Static and dynamicB 0 $$ {\mathrm{B}}_0 $$ field imperfections are detrimental to functional MRI (fMRI) applications, especially at ultra-high magnetic fields (UHF). In this work, a field camera is used to assess the benefits of retrospectively correctingB 0 $$ {\mathrm{B}}_0 $$ field perturbations on Blood Oxygen Level Dependent (BOLD) sensitivity in non-Cartesian three-dimensional (3D)-SPARKLING fMRI acquisitions.

METHODS

fMRI data were acquired at 1 mm3 $$ {}^3 $$ and for a 2.4s-TR while concurrently monitoring in real-time field perturbations using a Skope Clip-on field camera in a novel experimental setting involving a shorter TR than the required minimal TR of the field probes. Measurements of the dynamic field deviations were used along with a staticΔ B 0 $$ \Delta {\mathrm{B}}_0 $$ map to retrospectively correct static and dynamic field imperfections, respectively. In order to evaluate the impact of such a correction on fMRI volumes, a comparative study was conducted on healthy volunteers.

RESULTS

Correction ofB 0 $$ {\mathrm{B}}_0 $$ deviations improved image quality and yielded between 20% and 30% increase in median temporal signal-to-noise ratio (tSNR).Using fMRI data collected during a retinotopic mapping experiment, we demonstrated a significant increase in sensitivity to the BOLD contrast and improved accuracy of the BOLD phase maps: 44% (resp., 159%) more activated voxels were retrieved when using a significance control level based on a p-value of 0.001 without correcting for multiple comparisons (resp., 0.05 with a false discovery rate correction).

CONCLUSION

3D-SPARKLING fMRI hugely benefits from static and dynamicB 0 $$ {\mathrm{B}}_0 $$ imperfections correction. However, the proposed experimental protocol is flexible enough to be deployed on a large spectrum of encoding schemes, including arbitrary non-Cartesian readouts.

NMR Magnetometer Based on Dynamic Nuclear-Polarization for Low-Strength Magnetic Field Measurement

Nuclear magnetic resonance (NMR) magnetometers are considered due to their ability to map magnetic fields with high precision and calibrate other magnetic field measurement devices. However, the low signal-to-noise ratio of low-strength magnetic fields limits the precision when measuring magnetic fields below 40 mT. Therefore, we developed a new NMR magnetometer that combines the dynamic nuclear polarization (DNP) technique with pulsed NMR. The dynamic pre-polarization technique enhances the SNR under a low magnetic field. Pulsed NMR was used in conjunction with DNP to improve measurement accuracy and speed. The efficacy of this approach was validated through simulation and analysis of the measurement process. Next, a complete set of equipment was constructed, and we successfully measured magnetic fields of 30 mT and 8 mT with an accuracy of only 0.5 Hz (11 nT) at 30 mT (0.4 ppm) and 1 Hz (22 nT) at 8mT (3 ppm).

Real-time shimming with FID navigators

PURPOSE

To implement a method for real-time field control using rapid FID navigator (FIDnav) measurements and evaluate the efficacy of the proposed approach for mitigating dynamic field perturbations and improvingT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality.

METHODS

FIDnavs were embedded in a gradient echo sequence and a subject-specific linear calibration model was generated on the scanner to facilitate rapid shim updates in response to measured FIDnav signals. To confirm the accuracy of FID-navigated field updates, phantom and volunteer scans were performed with online updates of the scanner B0 shim settings. To evaluate improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality with real-time shimming, 10 volunteers were scanned at 3T while performing deep-breathing and nose-touching tasks designed to modulate the B0 field. Quantitative image quality metrics were compared with and without FID-navigated field control. An additional volunteer was scanned at 7T to evaluate performance at ultra-high field.

RESULTS

Applying measured FIDnav shim updates successfully compensated for applied global and linear field offsets in phantoms and across all volunteers. FID-navigated real-time shimming led to a substantial reduction in field fluctuations and a consequent improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality in volunteers performing deep-breathing and nose-touching tasks, with 7.57% ± 6.01% and 8.21% ± 10.90% improvement in peak SNR and structural similarity, respectively.

CONCLUSION

FIDnavs facilitate rapid measurement and application of field coefficients for slice-wise B0 shimming. The proposed approach can successfully counteract spatiotemporal field perturbations and substantially improvesT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality, which is important for a variety of clinical and research applications, particularly at ultra-high field.

Integration of an RF coil and commercial field camera for ultrahigh-field MRI

PURPOSE

To develop an RF coil with an integrated commercial field camera for ultrahigh field (7T) neuroimaging. The RF coil would operate within a head-only gradient coil and be subject to the corresponding design constraints. The RF coil can thereafter be used for subject-specific correction of k-space trajectories-notably in gradient-sensitive sequences such as single-shot spiral imaging.

METHODS

The transmit and receive performance was evaluated before and after the integration of field probes, whereas field probes were evaluated when in an optimal configuration external to the coil and after their integration. Diffusion-weighted EPI and single-shot spiral acquisitions were employed to evaluate the efficacy of correcting higher order field perturbations and the consequent effect on image quality.

RESULTS

Field probes had a negligible effect on RF-coil performance, including the transmit efficiency, transmit uniformity, and mean SNR over the brain. Modest reductions in field-probe signal lifetimes were observed, caused primarily by nonidealities in the gradient and shim fields of the head-only gradient coil at the probe positions. The field-monitoring system could correct up to second-order field perturbations in single-shot spiral imaging.

CONCLUSION

The integrated RF coil and field camera was capable of concurrent-field monitoring within a 7T head-only scanner and facilitated the subsequent correction of k-space trajectories during spiral imaging.

Thermal variation in gradient response: measurement and modeling

PURPOSE

Many aspects and imperfections of gradient dynamics in MRI have been successfully captured by linear time-invariant (LTI) models. Changes in gradient behavior due to heating, however, violate time invariance. The goal of this work is to study such changes at the level of transfer functions and model them by thermal extension of the LTI framework.

METHODS

To study the impact of gradient heating on transfer functions, a clinical MR system was heated using a range of high-amplitude DC and AC waveforms, each followed by measuring transfer functions in rapid succession while the system cooled down. Simultaneously, gradient temperature was monitored with an array of temperature sensors positioned according to initial infrared recordings of the gradient tube. The relation between temperatures and transfer functions is cast into local and global linear models. The models are analysed in terms of self-consistency, conditioning, and prediction performance.

RESULTS

Pronounced thermal effects are observed in the time resolved transfer functions, largely attributable to in-coil eddy currents and mechanical resonances. Thermal modeling is found to capture these effects well. The keys to good model performance are well-placed temperature sensors and suitable training data.

CONCLUSION

Heating changes gradient response, violating time invariance. The utility of LTI modeling can nevertheless be recovered by a linear thermal extension, relying on temperature sensing and adequate one-time training.

The backbone network of dynamic functional connectivity

Temporal networks have become increasingly pervasive in many real-world applications, including the functional connectivity analysis of spatially separated regions of the brain. A major challenge in analysis of such networks is the identification of noise confounds, which introduce temporal ties that are nonessential, or links that are formed by chance due to local properties of the nodes. Several approaches have been suggested in the past for static networks or temporal networks with binary weights for extracting significant ties whose likelihood cannot be reduced to the local properties of the nodes. In this work, we propose a data-driven procedure to reveal the irreducible ties in dynamic functional connectivity of resting-state fMRI data with continuous weights. This framework includes a null model that estimates the latent characteristics of the distributions of temporal links through optimization, followed by a statistical test to filter the links whose formation can be reduced to the activities and local properties of their interacting nodes. We demonstrate the benefits of this approach by applying it to a resting-state fMRI dataset, and provide further discussion on various aspects and advantages of it.

Advances in spiral fMRI: A high-resolution study with single-shot acquisition

Spiral fMRI has been put forward as a viable alternative to rectilinear echo-planar imaging, in particular due to its enhanced average k-space speed and thus high acquisition efficiency. This renders spirals attractive for contemporary fMRI applications that require high spatiotemporal resolution, such as laminar or columnar fMRI. However, in practice, spiral fMRI is typically hampered by its reduced robustness and ensuing blurring artifacts, which arise from imperfections in both static and dynamic magnetic fields. Recently, these limitations have been overcome by the concerted application of an expanded signal model that accounts for such field imperfections, and its inversion by iterative image reconstruction. In the challenging ultra-high field environment of 7 Tesla, where field inhomogeneity effects are aggravated, both multi-shot and single-shot 2D spiral imaging at sub-millimeter resolution was demonstrated with high depiction quality and anatomical congruency. In this work, we further these advances towards a time series application of spiral readouts, namely, single-shot spiral BOLD fMRI at 0.8 mm in-plane resolution. We demonstrate that high-resolution spiral fMRI at 7 T is not only feasible, but delivers both excellent image quality, BOLD sensitivity, and spatial specificity of the activation maps, with little artifactual blurring. Furthermore, we show the versatility of the approach with a combined in/out spiral readout at a more typical resolution (1.5 mm), where the high acquisition efficiency allows to acquire two images per shot for improved sensitivity by echo combination.

Reliability of task-specific neuronal activation assessed with functional PET, ASL and BOLD imaging

Mapping the neuronal response during cognitive processing is of crucial importance to gain new insights into human brain function. BOLD imaging and ASL are established MRI methods in this endeavor. Recently, the novel approach of functional PET (fPET) was introduced, enabling absolute quantification of glucose metabolism at rest and during task execution in a single measurement. Here, we report test-retest reliability of fPET in direct comparison to BOLD imaging and ASL. Twenty healthy subjects underwent two PET/MRI measurements, providing estimates of glucose metabolism, cerebral blood flow (CBF) and blood oxygenation. A cognitive task was employed with different levels of difficulty requiring visual-motor coordination. Task-specific neuronal activation was robustly detected with all three imaging approaches. The highest reliability was obtained for glucose metabolism at rest. Although this dropped during task performance it was still comparable to that of CBF. In contrast, BOLD imaging yielded high performance only for qualitative spatial overlap of task effects but not for quantitative comparison. Hence, the combined assessment of fPET and ASL offers reliable and simultaneous absolute quantification of glucose metabolism and CBF at rest and task.

Feasibility of spiral fMRI based on an LTI gradient model

Spiral imaging is very well suited for functional MRI, however its use has been limited by the fact that artifacts caused by gradient imperfections and B₀ inhomogeneity are more difficult to correct compared to EPI. Effective correction requires accurate knowledge of the traversed k-space trajectory. With the goal of making spiral fMRI more accessible, we have evaluated image reconstruction using trajectories predicted by the gradient impulse response function (GIRF), which can be determined in a one-time calibration step. GIRF-predicted reconstruction was tested for high-resolution (0.8 mm) fMRI at 7T. Image quality and functional results of the reconstructions using GIRF-prediction were compared to reconstructions using the nominal trajectory and concurrent field monitoring. The reconstructions using nominal spiral trajectories contain substantial artifacts and the activation maps contain misplaced activation. Image artifacts are substantially reduced when using the GIRF-predicted reconstruction, and the activation maps for the GIRF-predicted and monitored reconstructions largely overlap. The GIRF reconstruction provides a large increase in the spatial specificity of the activation compared to the nominal reconstruction. The GIRF-reconstruction generates image quality and fMRI results similar to using a concurrently monitored trajectory. The presented approach does not prolong or complicate the fMRI acquisition. Using GIRF-predicted trajectories has the potential to enable high-quality spiral fMRI in situations where concurrent trajectory monitoring is not available.

High-resolution fMRI at 7 Tesla: challenges, promises and recent developments for individual-focused fMRI studies

Limited detection power has been a bottleneck for subject-specific functional MRI (fMRI) studies, however the higher signal-to-noise ratio afforded by ultra-high magnetic fields (≥ 7 Tesla) provides levels of sensitivity and resolution needed to study individual subjects. What may be surprising is that higher imaging resolution may provide both higher specificity and sensitivity due to reductions in partial volume effects and reduced physiological noise. However, challenges remain to ensure high data quality and to reduce variability in ultra-high field fMRI. We discuss session-specific biases including those caused by factors related to instrumentation, anatomy, and physiology-which can translate into variability across sessions-and how to minimize these to help ultra-high field fMRI reach its full potential for individual-focused studies.

FMRI based on transition-band balanced SSFP in comparison with EPI on a high-performance 0.55 T scanner

PURPOSE

Low-field (<1 tesla) MRI scanners allow more widespread diagnostic use for a range of cardiac, musculoskeletal, and neurological applications. However, the feasibility of performing robust fMRI at low field has yet to be fully demonstrated. To address this gap, we investigated task-based fMRI using a highly sensitive transition-band balanced steady-state free precession approach and standard EPI on a 0.55 tesla scanner equipped with modern high-performance gradient coils and a receive array.

METHODS

TR and flip-angle of transition-band steady-state free precession were optimized for 0.55 tesla by simulations. Static shimming was employed to compensate for concomitant field effects. Visual task-based fMRI data were acquired from 8 healthy volunteers. For comparison, standard EPI data were also acquired with TE = T 2 ∗ . Retrospective image-based correction for physiological effects (RETROICOR) was used to quantify physiological noise effects.

RESULTS

Activation was robustly detected using both methods in a 4-min scan time. Transition-band steady-state free precession was found to be sensitive to interference from subtle spatial and temporal (field drift, respiration) variations in the magnetic field, counteracting potential advantages of the reduced magnetic susceptibility effects compared to its utilization at high field. These adverse effects could be partially remedied with static shimming and postprocessing approaches. Standard EPI proved more robust against the sources of interference.

CONCLUSION

BOLD contrast is sufficiently large at 0.55 tesla for robust detection of brain activation and may be employed to broaden the spectrum of applications of low-field MRI. Standard EPI outperforms transition-band steady-state free precession in terms of signal stability.

Ground-truth "resting-state" signal provides data-driven estimation and correction for scanner distortion of fMRI time-series dynamics

The fMRI community has made great strides in decoupling neuronal activity from other physiologically induced T₂* changes, using sensors that provide a ground-truth with respect to cardiac, respiratory, and head movement dynamics. However, blood oxygenation level-dependent (BOLD) time-series dynamics are also confounded by scanner artifacts, in complex ways that can vary not only between scanners but even, for the same scanner, between sessions. Unfortunately, the lack of an equivalent ground truth for BOLD time-series has thus far stymied the development of reliable methods for identification and removal of scanner-induced noise, a problem that we have previously shown to severely impact detection sensitivity of resting-state brain networks. To address this problem, we first designed and built a phantom capable of providing dynamic signals equivalent to that of the resting-state brain. Using the dynamic phantom, we then compared the ground-truth time-series with its measured fMRI data. Using these, we introduce data-quality metrics: Standardized Signal-to-Noise Ratio (ST-SNR) and Dynamic Fidelity that, unlike currently used measures such as temporal SNR (tSNR), can be directly compared across scanners. Dynamic phantom data acquired from four "best-case" scenarios: high-performance scanners with MR-physicist-optimized acquisition protocols, still showed scanner instability/multiplicative noise contributions of about 6-18% of the total noise. We further measured strong non-linearity in the fMRI response for all scanners, ranging between 8-19% of total voxels. To correct scanner distortion of fMRI time-series dynamics at a single-subject level, we trained a convolutional neural network (CNN) on paired sets of measured vs. ground-truth data. The CNN learned the unique features of each session's noise, providing a customized temporal filter. Tests on dynamic phantom time-series showed a 4- to 7-fold increase in ST-SNR and about 40-70% increase in Dynamic Fidelity after denoising, with CNN denoising outperforming both the temporal bandpass filtering and denoising using Marchenko-Pastur principal component analysis. Critically, we observed that the CNN temporal denoising pushes ST-SNR to a regime where signal power is higher than that of noise (ST-SNR > 1). Denoising human-data with ground-truth-trained CNN, in turn, showed markedly increased detection sensitivity of resting-state networks. These were visible even at the level of the single-subject, as required for clinical applications of fMRI.

Effect of Physiological Noise on Thoracolumbar Spinal Cord Functional Magnetic Resonance Imaging in 3T Magnetic Field

Introduction:

Functional Magnetic Resonance Imaging (fMRI) methods have been used to study sensorimotor processing in the spinal cord. However, these techniques confront unwanted noises to the measured signal from the physiological fluctuations. In the spinal cord imaging, most of the challenges are consequences of cardiac and respiratory movement artifacts that are considered as significant sources of noise, especially in the thoracolumbar region. In this study, we investigated the effect of each source of physiological noise and their contribution to the outcome of the analysis of the blood-oxygen-level-dependent signal in the human thoracolumbar spinal cord.

Methods:

Fifteen young healthy male volunteers participated in the study, and pain stimuli were delivered on the L5 dermatome between the two malleoli. Respiratory and cardiac signals were recorded during the imaging session, and the generated respiration and cardiac regressors were included in the general linear model for quantification of the effect of each of them on the task-analysis results. The sum of active voxels of the clusters was calculated in the spinal cord in three correction states (respiration correction only, cardiac correction only, and respiration and cardiac noise corrections) and analyzed with analysis of variance statistical test and receiver operating characteristic curve.

Results:

The results illustrated that cardiac noise correction had an effective role in increasing the active voxels (Mean±SD = 23.46±9.46) compared to other noise correction methods. Cardiac effects were higher than other physiological noise sources

Conclusion:

In summary, our results indicate great respiration effects on the lumbar and thoracolumbar spinal cord fMRI, and its contribution to the heartbeat effect can be a significant variable in the individual fMRI data analysis. Displacement of the spinal cord and the effects of this noise in the thoracolumbar and lumbar spinal cord fMRI results are significant and cannot be ignored.

Dynamic distortion correction for functional MRI using FID navigators

PURPOSE

To develop a method for slice-wise dynamic distortion correction for EPI using rapid spatiotemporal B₀ field measurements from FID navigators (FIDnavs) and to evaluate the efficacy of this new approach relative to an established data-driven technique.

METHODS

A low-resolution reference image was used to create a forward model of FIDnav signal changes to enable estimation of spatiotemporal B₀ inhomogeneity variations up to second order from measured FIDnavs. Five volunteers were scanned at 3 T using a 64-channel coil with FID-navigated EPI. The accuracy of voxel shift measurements and geometric distortion correction was assessed for experimentally induced magnetic field perturbations. The temporal SNR was evaluated in EPI time-series acquired at rest and with a continuous nose-touching action, before and after image realignment.

RESULTS

Field inhomogeneity coefficients and voxel shift maps measured using FIDnavs were in excellent agreement with multi-echo EPI measurements. The FID-navigated distortion correction accurately corrected image geometry in the presence of induced magnetic field perturbations, outperforming the data-driven approach in regions with large field offsets. In functional MRI scans with nose touching, FIDnav-based correction yielded temporal SNR gains of 30% in gray matter. Following image realignment, which accounted for global image shifts, temporal SNR gains of 3% were achieved.

CONCLUSIONS

Our proposed application of FIDnavs enables slice-wise dynamic distortion correction with high temporal efficiency. We achieved improved signal stability by leveraging the encoding information from multichannel coils. This approach can be easily adapted to other EPI-based sequences to improve temporal SNR for a variety of clinical and research applications.

A novel sequence to improve auditory functional MRI with variable silent delays

PURPOSE

Auditory functional MRI (fMRI) often uses silent inter-volume delays for stimulus presentation. However, maintaining the steady-state of the magnetization usually requires constant delays. Here, a novel acquisition scheme dubbed "pre-Saturated EPI using Multiple delays in Steady-state" (SEPIMS) is proposed, using spin saturation at a fixed delay before each volume to maintain steady-state conditions, independent of previous spin history. This concept allows for variable inter-volume delays and thus for flexible stimulus design in auditory fMRI. The purpose was to compare the signal stability of SEPIMS and conventional sparse EPI (CS-EPI).

METHODS

The saturation module comprises two non-selective adiabatic saturation pulses. The efficiency of the saturation and its effect on the SEPIMS signal stability is tested in vitro and in vivo.

RESULTS

Data show that SEPIMS yields the same signal stability as CS-EPI, even for extreme variations between inter-volume delay durations. However, dual saturation pulses are required to achieve sufficiently high saturation efficiency in compartments with long T₁ values. Importantly, spoiler gradient pulses after the EPI readout have to be optimized to avoid eddy-current-induced image distortions.

CONCLUSION

The proposed SEPIMS sequence maintains high signal stability in the presence of variable inter-volume durations, thus allowing for flexible stimulus design.

Gradient Response Harvesting for Continuous System Characterization During MR Sequences

MRI gradient systems are required to generate magnetic field gradient waveforms with very high fidelity. This is commonly implemented by gradient system calibration and pre-emphasis. However, a number of mechanisms, particularly thermal changes, cause variation in the gradient response over time, which cannot be addressed by calibration approaches. To overcome this limitation, we present a novel method termed gradient response harvesting, where the gradient response is continuously characterized during the course of a normal MR sequence. Snippets of field measurements are repeatedly acquired during an MR sequence, and from these multiple field measurements and the known nominal MR sequence gradients, the gradient response and gradient/field offsets are calculated. The calculation is implemented in a model-based and a model-free variant. The method is demonstrated for EPI with high gradient duty-cycle, where the continuous gradient characterization is used to obtain k-space trajectory estimates that are employed in the subsequent image reconstruction. During the course of the MR sequence, changes in both the envelope and the phase of the gradient response functions were observed, including shifts of mechanical resonances. The gradient response changes were also reflected in the calculated uninterrupted gradient waveforms and thus in the k-space trajectories. Using the updated encoding information in the image reconstruction removed ghosting artifacts, that otherwise impaired the image quality. We introduced the concept of gradient response harvesting and demonstrated its feasibility. The obtained gradient response functions may be used for quality assurance/preventive maintenance, real-time adaptation of gradient pre-emphasis or to calculate uninterrupted gradient field evolutions.

Characteristics of respiratory measures in young adults scanned at rest, including systematic changes and "missed" deep breaths

Breathing rate and depth influence the concentration of carbon dioxide in the blood, altering cerebral blood flow and thus functional magnetic resonance imaging (fMRI) signals. Such respiratory fluctuations can have substantial influence in studies of fMRI signal covariance in subjects at rest, the so-called "resting state functional connectivity" technique. If respiration is monitored during fMRI scanning, it is typically done using a belt about the subject's abdomen to record abdominal circumference. Several measures have been derived from these belt records, including the windowed envelope of the waveform (ENV), the windowed variance in the waveform (respiration variation, RV), and a measure of the amplitude of each breath divided by the cycle time of the breath (respiration volume per time, RVT). Any attempt to gauge respiratory contributions to fMRI signals requires a respiratory measure, but little is known about how these measures compare to each other, or how they perform beyond the small studies in which they were initially proposed. Here, we examine the properties of these measures in hundreds of healthy young adults scanned for an hour each at rest, a subset of the Human Connectome Project chosen for having high-quality physiological records. We find: 1) ENV, RV, and RVT are all correlated, and ENV and RV are more highly correlated to each other than to RVT; 2) respiratory events like deep breaths exhibit characteristic heart rate elevations, fMRI signal changes, head motions, and image quality abnormalities time-locked to large deflections in the belt traces; 3) all measures can "miss" deep breaths; 4) RVT "misses" deep breaths more than ENV or RV; 5) all respiratory measures change systematically over the course of a 14.4-min scan. We discuss the implications of these findings for the literature and ways to move forward in modeling respiratory influences on fMRI scans.

Identification of physiological response functions to correct for fluctuations in resting-state fMRI related to heart rate and respiration

Functional magnetic resonance imaging (fMRI) is widely viewed as the gold standard for studying brain function due to its high spatial resolution and non-invasive nature. However, it is well established that changes in breathing patterns and heart rate strongly influence the blood oxygen-level dependent (BOLD) fMRI signal and this, in turn, can have considerable effects on fMRI studies, particularly resting-state studies. The dynamic effects of physiological processes are often quantified by using convolution models along with simultaneously recorded physiological data. In this context, physiological response function (PRF) curves (cardiac and respiratory response functions), which are convolved with the corresponding physiological fluctuations, are commonly employed. While it has often been suggested that the PRF curves may be region- or subject-specific, it is still an open question whether this is the case. In the present study, we propose a novel framework for the robust estimation of PRF curves and use this framework to rigorously examine the implications of using population-, subject-, session- and scan-specific PRF curves. The proposed framework was tested on resting-state fMRI and physiological data from the Human Connectome Project. Our results suggest that PRF curves vary significantly across subjects and, to a lesser extent, across sessions from the same subject. These differences can be partly attributed to physiological variables such as the mean and variance of the heart rate during the scan. The proposed methodological framework can be used to obtain robust scan-specific PRF curves from data records with duration longer than 5 min, exhibiting significantly improved performance compared to previously defined canonical cardiac and respiration response functions. Besides removing physiological confounds from the BOLD signal, accurate modeling of subject- (or session-/scan-) specific PRF curves is of importance in studies that involve populations with altered vascular responses, such as aging subjects.

A single-chip integrated transceiver for high field NMR magnetometry

We present the design and performance of a broad-band single-chip integrated transceiver specifically conceived for nuclear magnetic resonance magnetometry. The single-chip transceiver is realized using a standard silicon complementary metal-oxide-semiconductor integrated circuit technology. A radio-frequency (RF) transmit amplifier, a transmit/receive switch, a low noise RF receive amplifier, a quadrature (IQ)-mixer, and two intermediate frequency amplifiers are integrated on a single silicon chip of 1.8 mm². The advantages and problematic aspects with respect to conventional discrete electronic approaches are discussed. We show the results of magnetic field measurements performed at 1.4 and 7.05 T, using solid and liquid samples having volumes from 40 μl down to 100 pl. Particular attention is devoted to the comparison of the experimentally measured magnetic field standard deviation with respect to the Cramer-Rao lower bound value. With a sample of distilled water (T₁ ≅ T₂ ≅ 3 s, T₂ ^*≅ 20 ms) having a volume of 40 μl, a standard deviation of 2.5 nT at 7.05 T (i.e., 0.5 ppb) in 1 s of averaging time is achieved, with a projected Cramer-Rao lower bond of 8 pT (i.e., 1.1 ppt).

Integrated fMRI Preprocessing Framework Using Extended Kalman Filter for Estimation of Slice-Wise Motion

Functional MRI acquisition is sensitive to subjects' motion that cannot be fully constrained. Therefore, signal corrections have to be applied a posteriori in order to mitigate the complex interactions between changing tissue localization and magnetic fields, gradients and readouts. To circumvent current preprocessing strategies limitations, we developed an integrated method that correct motion and spatial low-frequency intensity fluctuations at the level of each slice in order to better fit the acquisition processes. The registration of single or multiple simultaneously acquired slices is achieved online by an Iterated Extended Kalman Filter, favoring the robust estimation of continuous motion, while an intensity bias field is non-parametrically fitted. The proposed extraction of gray-matter BOLD activity from the acquisition space to an anatomical group template space, taking into account distortions, better preserves fine-scale patterns of activity. Importantly, the proposed unified framework generalizes to high-resolution multi-slice techniques. When tested on simulated and real data the latter shows a reduction of motion explained variance and signal variability when compared to the conventional preprocessing approach. These improvements provide more stable patterns of activity, facilitating investigation of cerebral information representation in healthy and/or clinical populations where motion is known to impact fine-scale data.

Serial correlations in single-subject fMRI with sub-second TR

When performing statistical analysis of single-subject fMRI data, serial correlations need to be taken into account to allow for valid inference. Otherwise, the variability in the parameter estimates might be under-estimated resulting in increased false-positive rates. Serial correlations in fMRI data are commonly characterized in terms of a first-order autoregressive (AR) process and then removed via pre-whitening. The required noise model for the pre-whitening depends on a number of parameters, particularly the repetition time (TR). Here we investigate how the sub-second temporal resolution provided by simultaneous multislice (SMS) imaging changes the noise structure in fMRI time series. We fit a higher-order AR model and then estimate the optimal AR model order for a sequence with a TR of less than 600 ms providing whole brain coverage. We show that physiological noise modelling successfully reduces the required AR model order, but remaining serial correlations necessitate an advanced noise model. We conclude that commonly used noise models, such as the AR(1) model, are inadequate for modelling serial correlations in fMRI using sub-second TRs. Rather, physiological noise modelling in combination with advanced pre-whitening schemes enable valid inference in single-subject analysis using fast fMRI sequences.