Accelerated multislice MRI with patterned excitation

PURPOSE

Accelerate multislice 2D MRI by using RF pulses that simultaneously act on different slices to combine contrast preparation and image acquisition.

THEORY AND METHODS

MRI applications often require the use of multiple RF pulses to generate desired contrast and prepare the signal for readout. Examples are the use of inversion prepulses to generate T1 contrast, or the use of spin-echo preparations to generate T2 or diffusion contrast. In multislice MRI, this separation of contrast preparation and readout can render scans time-inefficient and lengthy. We introduce a class of pulse sequences that overcomes this inefficiency by combining contrast preparation and signal readout. This is accomplished by using RF pulses that manipulate the magnetization of multiple slices simultaneously and a gradient crusher scheme that selects a target slice for readout.

RESULTS

Feasibility of the method was demonstrated for spin echo-based measurement of water diffusion and tissue pulsation in human brain at 3 T. Increases in time-efficiency and reductions in scan time were highly dependent on specific implementation and reached as high as 25% and 53%, respectively.

CONCLUSION

A novel approach to multislice MRI is demonstrated that reduces scan time for specific applications.

Cortical lesions uniquely predict motor disability accrual and form rarely in the absence of new white matter lesions in multiple sclerosis

Background and objectives

Cortical lesions (CL) are common in multiple sclerosis (MS) and associate with disability and progressive disease. We asked whether CL continue to form in people with stable white matter lesions (WML) and whether the association of CL with worsening disability relates to pre-existing or new CL.

Methods

A cohort of adults with MS were evaluated annually with 7 tesla (T) brain magnetic resonance imaging (MRI) and 3T brain and spine MRI for 2 years, and clinical assessments for 3 years. CL were identified on 7T images at each timepoint. WML and brain tissue segmentation were performed using 3T images at baseline and year 2.

Results

59 adults with MS had ≥1 7T follow-up visit (mean follow-up time 2±0.5 years). 9 had "active" relapsing-remitting MS (RRMS), defined as new WML in the year prior to enrollment. Of the remaining 50, 33 had "stable" RRMS, 14 secondary progressive MS (SPMS), and 3 primary progressive MS. 16 total new CL formed in the active RRMS group (median 1, range 0-10), 7 in the stable RRMS group (median 0, range 0-5), and 4 in the progressive MS group (median 0, range 0-1) (p=0.006, stable RR vs PMS p=0.88). New CL were not associated with greater change in any individual disability measure or in a composite measure of disability worsening (worsening Expanded Disability Status Scale or 9-hole peg test or 25-foot timed walk). Baseline CL volume was higher in people with worsening disability (median 1010μl, range 13-9888 vs median 267μl, range 0-3539, p=0.001, adjusted for age and sex) and in individuals with RRMS who subsequently transitioned to SPMS (median 2183μl, range 270-9888 vs median 321μl, range 0-6392 in those who remained RRMS, p=0.01, adjusted for age and sex). Baseline WML volume was not associated with worsening disability or transition from RRMS to SPMS.

Discussion

CL formation is rare in people with stable WML, even in those with worsening disability. CL but not WML burden predicts future worsening of disability, suggesting that the relationship between CL and disability progression is related to long-term effects of lesions that form in the earlier stages of disease, rather than to ongoing lesion formation.

Effect of motion, cortical orientation and spatial resolution on quantitative imaging of cortical R2* and magnetic susceptibility at 0.3 mm in-plane resolution at 7 T

MR images of the effective relaxation rate R2* and magnetic susceptibility χ derived from multi-echo T2*-weighted (T2*w) MRI can provide insight into iron and myelin distributions in the brain, with the potential of providing biomarkers for neurological disorders. Quantification of R2* and χ at submillimeter resolution in the cortex in vivo has been difficult because of challenges such as head motion, limited signal to noise ratio, long scan time, and motion related magnetic field fluctuations. This work aimed to improve the robustness for quantifying intracortical R2* and χ and analyze the effects from motion, spatial resolution, and cortical orientation. T2*w data was acquired with a spatial resolution of 0.3 × 0.3 × 0.4 mm3 at 7 T and downsampled to various lower resolutions. A combined correction for motion and B0 changes was deployed using volumetric navigators. Such correction improved the T2*w image quality rated by experienced image readers and test-retest reliability of R2* and χ quantification with reduced median inter-scan differences up to 10 s-1 and 5 ppb, respectively. R2* and χ near the line of Gennari, a cortical layer high in iron and myelin, were as much as 10 s-1 and 10 ppb higher than the region at adjacent cortical depth. In addition, a significant effect due to the cortical orientation relative to the static field (B0) was observed in χ with a peak-to-peak amplitude of about 17 ppb. In retrospectively downsampled data, the capability to distinguish different cortical depth regions based on R2* or χ contrast remained up to isotropic 0.5 mm resolution. This study highlights the unique characteristics of R2* and χ along the cortical depth at submillimeter resolution and the need for motion and B0 corrections for their robust quantification in vivo.

Cerebrovascular activity is a major factor in the cerebrospinal fluid flow dynamics

Cerebrospinal fluid (CSF) provides physical protection to the central nervous system as well as an essential homeostatic environment for the normal functioning of neurons. Additionally, it has been proposed that the pulsatile movement of CSF may assist in glymphatic clearance of brain metabolic waste products implicated in neurodegeneration. In awake humans, CSF flow dynamics are thought to be driven primarily by cerebral blood volume fluctuations resulting from a number of mechanisms, including a passive vascular response to blood pressure variations associated with cardiac and respiratory cycles. Recent research has shown that mechanisms that rely on the action of vascular smooth muscle cells ("cerebrovascular activity") such as neuronal activity, changes in intravascular CO2, and autonomic activation from the brainstem, may lead to CSF pulsations as well. Nevertheless, the relative contribution of these mechanisms to CSF flow remains unclear. To investigate this further, we developed an MRI approach capable of disentangling and quantifying CSF flow components of different time scales associated with these mechanisms. This approach was evaluated on human control subjects (n = 12) performing intermittent voluntary deep inspirations, by determining peak flow velocities and displaced volumes between these mechanisms in the fourth ventricle. We found that peak flow velocities were similar between the different mechanisms, while displaced volumes per cycle were about a magnitude larger for deep inspirations. CSF flow velocity peaked at around 10.4 s (range 7.1-14.8 s, n = 12) following deep inspiration, consistent with known cerebrovascular activation delays for this autonomic challenge. These findings point to an important role of cerebrovascular activity in the genesis of CSF pulsations. Other regulatory triggers for cerebral blood flow such as autonomic arousal and orthostatic challenges may create major CSF pulsatile movement as well. Future quantitative comparison of these and possibly additional types of CSF pulsations with the proposed approach may help clarify the conditions that affect CSF flow dynamics.

Multiple sclerosis cortical lesion detection with deep learning at ultra-high-field MRI

Manually segmenting multiple sclerosis (MS) cortical lesions (CLs) is extremely time consuming, and past studies have shown only moderate inter-rater reliability. To accelerate this task, we developed a deep-learning-based framework (CLAIMS: Cortical Lesion AI-Based Assessment in Multiple Sclerosis) for the automated detection and classification of MS CLs with 7 T MRI. Two 7 T datasets, acquired at different sites, were considered. The first consisted of 60 scans that include 0.5 mm isotropic MP2RAGE acquired four times (MP2RAGE×4), 0.7 mm MP2RAGE, 0.5 mm T₂ *-weighted GRE, and 0.5 mm T₂ *-weighted EPI. The second dataset consisted of 20 scans including only 0.75 × 0.75 × 0.9 mm³ MP2RAGE. CLAIMS was first evaluated using sixfold cross-validation with single and multi-contrast 0.5 mm MRI input. Second, the performance of the model was tested on 0.7 mm MP2RAGE images after training with either 0.5 mm MP2RAGE×4, 0.7 mm MP2RAGE, or alternating the two. Third, its generalizability was evaluated on the second external dataset and compared with a state-of-the-art technique based on partial volume estimation and topological constraints (MSLAST). CLAIMS trained only with MP2RAGE×4 achieved results comparable to those of the multi-contrast model, reaching a CL true positive rate of 74% with a false positive rate of 30%. Detection rate was excellent for leukocortical and subpial lesions (83%, and 70%, respectively), whereas it reached 53% for intracortical lesions. The correlation between disability measures and CL count was similar for manual and CLAIMS lesion counts. Applying a domain-scanner adaptation approach and testing CLAIMS on the second dataset, the performance was superior to MSLAST when considering a minimum lesion volume of 6 μL (lesion-wise detection rate of 71% versus 48%). The proposed framework outperforms previous state-of-the-art methods for automated CL detection across scanners and protocols. In the future, CLAIMS may be useful to support clinical decisions at 7 T MRI, especially in the field of diagnosis and differential diagnosis of MS patients.

Cortical lesion hotspots and association of subpial lesions with disability in multiple sclerosis

BACKGROUND

Dramatic improvements in visualization of cortical (especially subpial) multiple sclerosis (MS) lesions allow assessment of impact on clinical course.

OBJECTIVE

Characterize cortical lesions by 7 tesla (T) T₂^*-/T₁-weighted magnetic resonance imaging (MRI); determine relationship with other MS pathology and contribution to disability.

METHODS

Sixty-four adults with MS (45 relapsing-remitting/19 progressive) underwent 3 T brain/spine MRI, 7 T brain MRI, and clinical testing.

RESULTS

Cortical lesions were found in 94% (progressive: median 56/range 2-203; relapsing-remitting: 15/0-168; p = 0.004). Lesion distribution across 50 cortical regions was nonuniform (p = 0.006), with highest lesion burden in supplementary motor cortex and highest prevalence in superior frontal gyrus. Leukocortical and white matter lesion volumes were strongly correlated (r = 0.58, p < 0.0001), while subpial and white matter lesion volumes were moderately correlated (r = 0.30, p = 0.002). Leukocortical (p = 0.02) but not subpial lesions (p = 0.40) were correlated with paramagnetic rim lesions; both were correlated with spinal cord lesions (p = 0.01). Cortical lesion volumes (total and subtypes) were correlated with expanded disability status scale, 25-foot timed walk, nine-hole peg test, and symbol digit modality test scores.

CONCLUSION

Cortical lesions are highly prevalent and are associated with disability and progressive disease. Subpial lesion burden is not strongly correlated with white matter lesions, suggesting differences in inflammation and repair mechanisms.

Autonomic arousals contribute to brain fluid pulsations during sleep

During sleep, slow waves of neuro-electrical activity engulf the human brain and aid in the consolidation of memories. Recent research suggests that these slow waves may also promote brain health by facilitating the removal of metabolic waste, possibly by orchestrating the pulsatile flow of cerebro-spinal fluid (CSF) through local neural control over vascular tone. To investigate the role of slow waves in the generation of CSF pulsations, we analyzed functional MRI data obtained across the full sleep-wake cycle and during a respiratory task during wakefulness. This revealed a novel generating mechanism that relies on the autonomic regulation of cerebral vascular tone without requiring slow electrocortical activity or even sleep. Therefore, the role of CSF pulsations in brain waste clearance may, in part, depend on proper autoregulatory control of cerebral blood flow.

Sensitivity limitations of high-resolution perfusion-based human fMRI at 7 Tesla

The study of the brain's functional organization at laminar and columnar level of the cortex with blood oxygenation-level dependent (BOLD) functional MRI (fMRI) is affected by the contribution of large veins downstream from the microvascular response to brain activity. Blood volume- and especially perfusion-based techniques may reduce this problem because of their reduced sensitivity to venous effects, but may not allow the same spatial resolution because of smaller signal changes associated with brain activity. Here we investigated the practical resolution limits of perfusion-weighted fMRI in human visual stimulation experiments. For this purpose, we used a highly sensitive, single-shot perfusion labeling (SSPL) technique at 7 T and compared sensitivity to detect visual activation at low (2 mm, n = 10) and high (1 mm, n = 8) nominal isotropic spatial, and 3 s temporal, resolution with BOLD in 5½-minute-long experiments. Despite the smaller absolute signal change with activation, 2 mm resolution SSPL yielded comparable sensitivity to BOLD. This was attributed to a superior suppression of physiological noise with SSPL. However, at 1 mm nominal resolution, SSPL sensitivity fell on average at least 42% below that of BOLD, and detection of visual activation was compromised. This is explained by the fact that at high resolution, with both techniques, typically thermal noise rather than physiological noise dominates sensitivity. The observed sensitivity loss implies that to perform 1-mm resolution, perfusion weighted fMRI with a robustness similar to BOLD, scan times that are almost 3 times longer than the comparable BOLD experiment are required. This is in line with or slightly better than previous comparisons between perfusion-weighted fMRI and BOLD. The lower sensitivity has to be weighed against the spatial fidelity advantages of high-resolution perfusion-weighted fMRI.

Navigator-Guided Motion and B0 Correction of T2*-Weighted Magnetic Resonance Imaging Improves Multiple Sclerosis Cortical Lesion Detection

BACKGROUND

Cortical lesions are common in multiple sclerosis (MS). T2*-weighted (T2*w) imaging at 7 T is relatively sensitive for cortical lesions, but quality is often compromised by motion and main magnetic field (B0) fluctuations.

PURPOSE

The aim of this study was to determine whether motion and B0 correction with a navigator-guided gradient-recalled echo sequence can improve cortical lesion detection in T2*w magnetic resonance imaging.

MATERIALS AND METHODS

In this prospective study, a gradient-recalled echo sequence incorporating a navigator allowing for motion and B0 field correction was applied to collect T2*w images at 7 T from adults with MS between August 2019 and March 2020. T2*-weighted images were acquired in 1 to 3 partially overlapping scans per individual and were reconstructed using global average B0 correction ("uncorrected") or motion correction and spatially linear B0 correction ("corrected"). Image quality rating and manual segmentation of cortical lesions were performed on uncorrected and corrected images. Lesions seen on a single scan were retrospectively evaluated on the complementary scan. The association of cortical lesions with clinical disability was assessed. Mixed models were used to determine the effect of correction on lesion detection as well as on the relationship between disability and lesion count.

RESULTS

A total of 22 T2*w scans were performed on 11 adults with MS (mean [SD] age, 49 [11] years; 8 women). Quality improved for 20 of 22 scans (91%) after correction. A total of 69 cortical lesions were identified on uncorrected images (median per scan, 2; range, 0-11) versus 148 on corrected images (median per scan, 4.5; range, 0-25; rate ratio [RR], 2.1; P < 0.0001). For low-quality uncorrected scans with moderate to severe motion artifact (18/22, 82%), there was an improvement in cortical lesion detection with correction (RR, 2.5; P < 0.0001), whereas there was no significant change in cortical lesion detection for high-quality scans (RR, 1.3; P = 0.43).

CONCLUSIONS

Navigator-guided motion and B0 correction substantially improves the overall image quality of T2*w magnetic resonance imaging at 7 T and increases its sensitivity for cortical lesions.

fMRI-based detection of alertness predicts behavioral response variability

Levels of alertness are closely linked with human behavior and cognition. However, while functional magnetic resonance imaging (fMRI) allows for investigating whole-brain dynamics during behavior and task engagement, concurrent measures of alertness (such as EEG or pupillometry) are often unavailable. Here, we extract a continuous, time-resolved marker of alertness from fMRI data alone. We demonstrate that this fMRI alertness marker, calculated in a short pre-stimulus interval, captures trial-to-trial behavioral responses to incoming sensory stimuli. In addition, we find that the prediction of both EEG and behavioral responses during the task may be accomplished using only a small fraction of fMRI voxels. Furthermore, we observe that accounting for alertness appears to increase the statistical detection of task-activated brain areas. These findings have broad implications for augmenting a large body of existing datasets with information about ongoing arousal states, enriching fMRI studies of neural variability in health and disease.

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.

Eight-channel parallel transmit-receive system for 7 T MRI with optically controlled and monitored on-coil current-mode RF amplifiers

PURPOSE

To demonstrate a practical implementation of an eight-channel parallel-transmit system for brain imaging at 7 T based on on-coil amplifier technology.

METHODS

An eight-channel parallel transmit-receive system was built with optimized on-coil switch-mode current RF power amplifiers. The amplifiers were optically controlled from an eight-channel interface that was connected to a 7 T MRI scanner. The interface also optically received a down-converted version of the coil current sensed in each amplifier for monitoring and feedback adjustments.

RESULTS

Each on-coil amplifier delivered more than 100 W peak power and provided enough amplifier decoupling (<-15 dB) for the implemented eight-channel array configuration. Phantom and human images were acquired to demonstrate practical operation of this new technology in a 7 T MRI scanner.

CONCLUSION

Further development and improvement of previously demonstrated on-coil technology led to successful implementation of an eight-channel parallel-transmit system able to deliver strong B₁ fields for typical brain imaging applications. This is an important step forward toward implementation of on-coil RF amplification for high-field MRI.

B0-field dependence of MRI T1 relaxation in human brain

Tissue longitudinal relaxation characterized by recovery time T1 or rate R1 is a fundamental MRI contrast mechanism that is increasingly being used to study the brain's myelination patterns in both health and disease. Nevertheless, the quantitative relationship between T1 and myelination, and its dependence on B0 field strength, is still not well known. It has been theorized that in much of brain tissue, T1 field-dependence is driven by that of macromolecular protons (MP) through a mechanism called magnetization transfer (MT). Despite the explanatory power of this theory and substantial support from in-vitro experiments at low fields (<3 ​T), in-vivo evidence across clinically relevant field strengths is lacking. In this study, T1-weighted MRI was acquired in a group of eight healthy volunteers at four clinically relevant field strengths (0.55, 1.5, 3 and 7 ​T) using the same pulse sequence at a single site, and jointly analyzed based on the two-pool model of MT. MP fraction and free-water pool T1 were obtained in several brain structures at 3 and 7 ​T, which allowed distinguishing between contributions from macromolecular content and iron to tissue T1. Based on this, the T1 of MP in white matter, indirectly determined by assuming a field independent T1 of free water, was shown to increase approximately linearly with B0. This study advances our understanding of the T1 contrast mechanism and its relation to brain myelin content across the wide range of currently available MRI strengths, and it has the potential to inform design of T1 mapping methods for improved reproducibility in the human brain.

Reducing motion sensitivity in 3D high-resolution T2*-weighted MRI by navigator-based motion and nonlinear magnetic field correction

T₂*-weighted gradient echo (GRE) MRI at high field is uniquely sensitive to the magnetic properties of tissue and allows the study of brain and vascular anatomy at high spatial resolution. However, it is also sensitive to B₀ field changes induced by head motion and physiological processes such as the respiratory cycle. Conventional motion correction techniques do not take these field changes into account, and consequently do not fully recover image quality in T₂*-weighted MRI. Here, a novel approach was developed to address this by monitoring the B₀ field with a volumetric EPI phase navigator. The navigator was acquired at a shorter echo time than that of the (higher resolution) T₂*-weighted GRE imaging data and accelerated with parallel imaging for high temporal resolution. At 4 ​mm isotropic spatial resolution and 0.54 ​s temporal resolution, the accuracy for estimation of rotation and translation was better than 0.2° and 0.1 ​mm, respectively. The 10% and 90% percentiles of B₀ measurement error using the navigator were -1.8 and 1.5 Hz  at 7 T, respectively. A fast retrospective reconstruction algorithm correcting for both motion and nonlinear B₀ changes was also developed. The navigator and reconstruction algorithm were evaluated in correcting motion-corrupted high-resolution T₂*-weighted GRE MRI on healthy human subjects at 7 ​T. Excellent image quality was demonstrated with the proposed correction method.

All-night functional magnetic resonance imaging sleep studies

BACKGROUND

Previous functional magnetic resonance imaging (fMRI) sleep studies have been hampered by the difficulty of obtaining extended amounts of sleep in the sleep-adverse environment of the scanner and often have resorted to manipulations such as sleep depriving subjects before scanning. These manipulations limit the generalizability of the results.

NEW METHOD

The current study is a methodological validation of procedures aimed at obtaining all-night fMRI data in sleeping subjects with minimal exposure to experimentally induced sleep deprivation. Specifically, subjects slept in the scanner on two consecutive nights, allowing the first night to serve as an adaptation night.

RESULTS/COMPARISON WITH EXISTING METHOD(S)

Sleep scoring results from simultaneously acquired electroencephalography data on Night 2 indicate that subjects (n = 12) reached the full spectrum of sleep stages including slow-wave (M = 52.1 min, SD = 26.5 min) and rapid eye movement (REM, M = 45.2 min, SD = 27.9 min) sleep and exhibited a mean of 2.1 (SD = 1.1) nonREM-REM sleep cycles.

CONCLUSIONS

It was found that by diligently applying fundamental principles and methodologies of sleep and neuroimaging science, performing all-night fMRI sleep studies is feasible. However, because the two nights of the study were performed consecutively, some sleep deprivation from Night 1 as a cause of the Night 2 results is likely, so consideration should be given to replicating the current study with a washout period. It is envisioned that other laboratories can adopt the core features of this protocol to obtain similar results.

Effect of head motion on MRI B0 field distribution

PURPOSE

To identify and characterize the sources of B₀ field changes due to head motion, to reduce motion sensitivity in human brain MRI.

METHODS

B₀ fields were measured in 5 healthy human volunteers at various head poses. After measurement of the total field, the field originating from the subject was calculated by subtracting the external field generated by the magnet and shims. A subject-specific susceptibility model was created to quantify the contribution of the head and torso. The spatial complexity of the field changes was analyzed using spherical harmonic expansion.

RESULTS

Minor head pose changes can cause substantial and spatially complex field changes in the brain. For rotations and translations of approximately 5 º and 5 mm, respectively, at 7 T, the field change that is associated with the subject's magnetization generates a standard deviation (SD) of about 10 Hz over the brain. The stationary torso contributes to this subject-associated field change significantly with a SD of about 5 Hz. The subject-associated change leads to image-corrupting phase errors in multi-shot <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> -weighted acquisitions.

CONCLUSION

The B₀ field changes arising from head motion are problematic for multishot <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> -weighted imaging. Characterization of the underlying sources provides new insights into mitigation strategies, which may benefit from individualized predictive field models in addition to real-time field monitoring and correction strategies.

Synthesized tissue-equivalent dielectric phantoms using salt and polyvinylpyrrolidone solutions

PURPOSE

To explore the use of polyvinylpyrrolidone (PVP) for simulated materials with tissue-equivalent dielectric properties.

METHODS

PVP and salt were used to control, respectively, relative permittivity and electrical conductivity in a collection of 63 samples with a range of solute concentrations. Their dielectric properties were measured with a commercial probe and fitted to a 3D polynomial in order to establish an empirical recipe. The material's thermal properties and MR spectra were measured.

RESULTS

The empirical polynomial recipe (available at https://www.amri.ninds.nih.gov/cgi-bin/phantomrecipe) provides the PVP and salt concentrations required for dielectric materials with permittivity and electrical conductivity values between approximately 45 and 78, and 0.1 to 2 siemens per meter, respectively, from 50 MHz to 4.5 GHz. The second- (solute concentrations) and seventh- (frequency) order polynomial recipe provided less than 2.5% relative error between the measured and target properties. PVP side peaks in the spectra were minor and unaffected by temperature changes.

CONCLUSION

PVP-based phantoms are easy to prepare and nontoxic, and their semitransparency makes air bubbles easy to identify. The polymer can be used to create simulated material with a range of dielectric properties, negligible spectral side peaks, and long T₂ relaxation time, which are favorable in many MR applications. Magn Reson Med 80:413-419, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Optically controlled on-coil amplifier with RF monitoring feedback

PURPOSE

To develop a new optically controlled on-coil amplifier that facilitates safe use of multi-channel radiofrequency (RF) transmission in MRI by real-time monitoring of signal phase and amplitude.

METHODS

Monitoring was carried out with a 4-channel prototype system by sensing, down sampling, digitizing, and optically transmitting the RF transmit signal to a remote PC to control the amplifiers. Performance was evaluated with benchtop and 7 T MRI experiments.

RESULTS

Monitored amplitude and phase were stable across repetitions and had standard deviations of 0.061 μT and 0.0073 rad, respectively. The feedback system allowed inter-channel phase and B₁ amplitude to be adjusted within two iterations. MRI experiments demonstrated the feasibility of this approach to perform safe and accurate multi-channel RF transmission and monitoring at high field.

CONCLUSION

We demonstrated a 4-channel transceiver system based on optically controlled on-coil amplifiers with RF signal monitoring and feedback control. The approach allows the safe and precise control of RF transmission fields, required to achieve uniform excitation at high field. Magn Reson Med 79:2833-2841, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Linear Discriminant Analysis Achieves High Classification Accuracy for the BOLD fMRI Response to Naturalistic Movie Stimuli

Naturalistic stimuli like movies evoke complex perceptual processes, which are of great interest in the study of human cognition by functional MRI (fMRI). However, conventional fMRI analysis based on statistical parametric mapping (SPM) and the general linear model (GLM) is hampered by a lack of accurate parametric models of the BOLD response to complex stimuli. In this situation, statistical machine-learning methods, a.k.a. multivariate pattern analysis (MVPA), have received growing attention for their ability to generate stimulus response models in a data-driven fashion. However, machine-learning methods typically require large amounts of training data as well as computational resources. In the past, this has largely limited their application to fMRI experiments involving small sets of stimulus categories and small regions of interest in the brain. By contrast, the present study compares several classification algorithms known as Nearest Neighbor (NN), Gaussian Naïve Bayes (GNB), and (regularized) Linear Discriminant Analysis (LDA) in terms of their classification accuracy in discriminating the global fMRI response patterns evoked by a large number of naturalistic visual stimuli presented as a movie. Results show that LDA regularized by principal component analysis (PCA) achieved high classification accuracies, above 90% on average for single fMRI volumes acquired 2 s apart during a 300 s movie (chance level 0.7% = 2 s/300 s). The largest source of classification errors were autocorrelations in the BOLD signal compounded by the similarity of consecutive stimuli. All classifiers performed best when given input features from a large region of interest comprising around 25% of the voxels that responded significantly to the visual stimulus. Consistent with this, the most informative principal components represented widespread distributions of co-activated brain regions that were similar between subjects and may represent functional networks. In light of these results, the combination of naturalistic movie stimuli and classification analysis in fMRI experiments may prove to be a sensitive tool for the assessment of changes in natural cognitive processes under experimental manipulation.

Optically controlled switch-mode current-source amplifiers for on-coil implementation in high-field parallel transmission

PURPOSE

We tested the feasibility of implementing parallel transmission (pTX) for high-field MRI using a radiofrequency (RF) amplifier design to be located on or in the immediate vicinity of an RF transmit coil.

METHOD

We designed a current-source switch-mode amplifier based on miniaturized, nonmagnetic electronics. Optical RF carrier and envelope signals to control the amplifier were derived, through a custom-built interface, from the RF source accessible in the scanner control. Amplifier performance was tested by benchtop measurements as well as with imaging at 7T (300 MHz) and 11.7 T (500 MHz). The ability to perform pTX was evaluated by measuring interchannel coupling and phase adjustment in a two-channel setup.

RESULTS

The amplifier delivered in excess of 44 W RF power and caused minimal interference with MRI. The interface derived accurate optical control signals with carrier frequencies ranging from 64 to 750 MHz. Decoupling better than 14 dB was obtained between two coil loops separated by only 1 cm. Application to MRI was demonstrated by acquiring artifact-free images at 7 T and 11.7 T.

CONCLUSION

We propose an optically controlled miniaturized RF amplifier for on-coil implementation at high fields that should facilitate implementation of high-density pTX arrays. Magn Reson Med 76:340-349, 2016. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.

A 7T spine array based on electric dipole transmitters

PURPOSE

The goal of this study was to explore the feasibility of using an array of electric dipole antennas for RF transmission in spine MRI at high fields.

METHOD

A two-channel transmit array based on an electric dipole design was quantitatively optimized for 7T spine imaging and integrated with a receive array combining eight loop coils. Using B1+ mapping, the transmit efficiency of the dipole array was compared with a design using quadrature loop pairs. The radiofrequency energy deposition for each array was measured using a home-built dielectric phantom and MR thermometry. The performance of the proposed array was qualitatively demonstrated in human studies.

RESULTS

The results indicate dramatically improved transmit efficiency for the dipole design compared with the loop excitation. A gain of up to 76% was achieved within the spinal region.

CONCLUSION

For imaging of the spine, electric dipole-based transmitters provide an attractive alternative to the traditional loop-based design. Easy integration with existing receive array technology facilitates practical use at high fields.

Detection of demyelination in multiple sclerosis by analysis of [Formula: see text] relaxation at 7 T

Multiple sclerosis (MS) is a relatively common cause of inflammatory demyelinating lesions of the central nervous system. In an attempt to detect and characterize ongoing demyelination in MS patient brains, we used a novel magnetic resonance imaging (MRI) technique, involving the fitting of a three-component model to the [Formula: see text] relaxation behavior at high-field (7 T). This model allowed estimation of the amount of myelin water (and thus indirectly myelin content), axonal water, and interstitial water. In this study, 25 relapsing-remitting MS patients underwent a 7 T MRI from which 12 gadolinium-enhancing lesions, 61 non-enhancing lesions, and their corresponding contralateral normal appearing white matter (NAWM) regions were analyzed. In both enhancing and non-enhancing lesions, the amplitude of myelin water was significantly decreased, and interstitial and axonal water were increased relative to the contralateral NAWM. Longer relaxation time [Formula: see text] of interstitial and axonal water, and lower frequency shift of axonal water, were also observed in both enhancing and non-enhancing lesions when compared to the contralateral NAWM. No significant difference was found between enhancing lesions and non-enhancing lesions. These findings suggest that the fitting of a three-component model to the [Formula: see text] decay curve in MS lesions may help to quantify myelin loss.

A torque balance measurement of anisotropy of the magnetic susceptibility in white matter

PURPOSE

Recent MRI studies have suggested that the magnetic susceptibility of white matter (WM) in the human brain is anisotropic, providing a new contrast mechanism for the visualization of fiber bundles and allowing the extraction of cellular compartment-specific information. This study provides an independent confirmation and quantification of this anisotropy.

METHODS

Anisotropic magnetic susceptibility results in a torque exerted on WM when placed in a uniform magnetic field, tending to align the WM fibers with the field. To quantify the effect, excised spinal cord samples were placed in a torque balance inside the magnet of a 7 T MRI system and the magnetic torque was measured as function of orientation.

RESULTS

All tissue samples (n = 5) showed orienting effects, confirming the presence of anisotropic susceptibility. Analysis of the magnetic torque resulted in reproducible values for the WM volume anisotropy that ranged from 13.6 to 19.2 ppb.

CONCLUSION

The independently determined anisotropy values confirm estimates inferred from MRI experiments and validate the use of anisotropy to extract novel information about brain fiber structure and myelination.

Characterization of a dielectric phantom for high-field magnetic resonance imaging applications

PURPOSE

In this work, a generic recipe for an inexpensive and nontoxic phantom was developed within a range of biologically relevant dielectric properties from 150 MHz to 4.5 GHz.

METHODS

The recipe includes deionized water as the solvent, NaCl to primarily control conductivity, sucrose to primarily control permittivity, agar-agar to gel the solution and reduce heat diffusivity, and benzoic acid to preserve the gel. Two hundred and seventeen samples were prepared to cover the feasible range of NaCl and sucrose concentrations. Their dielectric properties were measured using a commercial dielectric probe and were fitted to a 3D polynomial to generate a recipe describing the properties as a function of NaCl concentration, sucrose concentration, and frequency.

RESULTS

Results indicated that the intuitive linear and independent relationships between NaCl and conductivity and between sucrose and permittivity are not valid. A generic polynomial recipe was developed to characterize the complex relationship between the solutes and the resulting dielectric values and has been made publicly available as a web application. In representative mixtures developed to mimic brain and muscle tissue, less than 2% difference was observed between the predicted and measured conductivity and permittivity values.

CONCLUSIONS

It is expected that the recipe will be useful for generating dielectric phantoms for general magnetic resonance imaging (MRI) coil development at high magnetic field strength, including coil safety evaluation as well as pulse sequence evaluation (including B₁(+) mapping, B₁(+) shimming, and selective excitation pulse design), and other non-MRI applications which require biologically equivalent dielectric properties.

Early anti-correlated BOLD signal changes of physiologic origin

Negative BOLD signals that are synchronous with resting state fluctuations have been observed in large vessels in the cortical sulci and surrounding the ventricles. In this study, we investigated the origin of these negative BOLD signals by applying a Cued Deep Breathing (CDB) task to create transient hypocapnia and a resultant global fMRI signal decrease. We hypothesized that a global stimulus would amplify the effect in large vessels and that using a global negative (vasoconstrictive) stimulus would test whether these voxels exhibit either inherently negative or simply anti-correlated BOLD responses. Significantly anti-correlated, but positive, BOLD signal changes during respiratory challenges were identified in voxels primarily located near edges of brain spaces containing CSF. These positive BOLD responses occurred earlier than the negative CDB response across most of gray matter voxels. These findings confirm earlier suggestions that in some brain regions, local, fractional changes in CSF volume may overwhelm BOLD-related signal changes, leading to signal anti-correlation. We show that regions with CDB anti-correlated signals coincide with most, but not all, of the regions with negative BOLD signal changes observed during a visual and motor stimulus task. Thus, the addition of a physiological challenge to fMRI experiments can help identify which negative BOLD signals are passive physiological anti-correlations and which may have a putative neuronal origin.

Independent sources of spontaneous BOLD fluctuation along the visual pathway

In resting-state functional magnetic resonance imaging (fMRI) experiments, correlation analysis can be used to identify clusters of cortical regions that may be functionally connected. Although such functional connectivity is often assumed to reflect cortico-cortical connections, a potential confound is the contribution of subcortical brain regions, many of which have strong anatomical connectivity to cortical regions and may also enable cortico-cortical interactions through trans-thalamic pathways. To investigate this, we performed resting state fMRI of the human visual system, including cortical regions and subcortical nuclei of the pulvinar and lateral geniculate. Regression analysis was used to investigate the dependence of the measured inter-regional correlations upon afferents from specific retinal, thalamic and cortical regions as well as systemic global signal fluctuation. A high level of inter-hemispheric correlation (cc = 0.95) was found in the visual cortex that could not be explained by activity in the subcortical nuclei investigated; in addition a relatively low level of inter-hemispheric correlation (cc = 0.39-0.42) was found in vision-related thalamic nuclei that could not be explained by direct anatomical connections or their cortical inputs. These findings suggest that spontaneous fMRI signal correlations within the human visual system originate from a mixture of independent signal sources that may be transmitted through thalamo-cortical, cortico-thalamic, and cortico-cortical connections either trans-callosal or trans-thalamic in origin. Our findings thus call for more cautious interpretation of resting state functional connectivity in terms of any single type of anatomical connectivity.

Neuroelectrical decomposition of spontaneous brain activity measured with functional magnetic resonance imaging

Spontaneous activity in the human brain occurs in complex spatiotemporal patterns that may reflect functionally specialized neural networks. Here, we propose a subspace analysis method to elucidate large-scale networks by the joint analysis of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data. The new approach is based on the notion that the neuroelectrical activity underlying the fMRI signal may have EEG spectral features that report on regional neuronal dynamics and interregional interactions. Applying this approach to resting healthy adults, we indeed found characteristic spectral signatures in the EEG correlates of spontaneous fMRI signals at individual brain regions as well as the temporal synchronization among widely distributed regions. These spectral signatures not only allowed us to parcel the brain into clusters that resembled the brain's established functional subdivision, but also offered important clues for disentangling the involvement of individual regions in fMRI network activity.

Micro-compartment specific T2* relaxation in the brain

MRI at high field can be sensitized to the magnetic properties of tissues, which introduces a signal dependence on the orientation of white matter (WM) fiber bundles relative to the magnetic field. In addition, study of the NMR relaxation properties of this signal has indicated contributions from compartmentalized water environments inside and outside the myelin sheath that may be separable. Here we further investigated the effects of water compartmentalization on the MRI signal with the goal of extracting compartment-specific information. By comparing MRI measurements of human and marmoset brain at 7T with magnetic field modeling, we show that: (1) water between the myelin lipid bilayers, in the axonal, and in the interstitial space each experience characteristic magnetic field effects that depend on fiber orientation (2) these field effects result in characteristic relaxation properties and frequency shifts for these compartments; and (3) compartmental contributions may be separated by multi-component fitting of the MRI signal relaxation (i.e. decay) curve. We further show the potential application of these findings to the direct mapping of myelin content and assessment of WM fiber integrity with high field MRI.

In vivo evaluation of the effect of stimulus distribution on FIR statistical efficiency in event-related fMRI

Technical developments in MRI have improved signal to noise, allowing use of analysis methods such as Finite impulse response (FIR) of rapid event related functional MRI (er-fMRI). FIR is one of the most informative analysis methods as it determines onset and full shape of the hemodynamic response function (HRF) without any a priori assumptions. FIR is however vulnerable to multicollinearity, which is directly related to the distribution of stimuli over time. Efficiency can be optimized by simplifying a design, and restricting stimuli distribution to specific sequences, while more design flexibility necessarily reduces efficiency. However, the actual effect of efficiency on fMRI results has never been tested in vivo. Thus, it is currently difficult to make an informed choice between protocol flexibility and statistical efficiency. The main goal of this study was to assign concrete fMRI signal to noise values to the abstract scale of FIR statistical efficiency. Ten subjects repeated a perception task with five random and m-sequence based protocol, with varying but, according to literature, acceptable levels of multicollinearity. Results indicated substantial differences in signal standard deviation, while the level was a function of multicollinearity. Experiment protocols varied up to 55.4% in standard deviation. Results confirm that quality of fMRI in an FIR analysis can significantly and substantially vary with statistical efficiency. Our in vivo measurements can be used to aid in making an informed decision between freedom in protocol design and statistical efficiency.

Finding thalamic BOLD correlates to posterior alpha EEG

Oscillatory electrical brain activity in the alpha (8-13 Hz) band is a prominent feature of human electroencephalography (EEG) during alert wakefulness, and is commonly thought to arise primarily from the occipital and parietal parts of the cortex. While the thalamus is considered to play a supportive role in the generation and modulation of cortical alpha rhythms, its precise function remains controversial and incompletely understood. To address this, we evaluated the correlation between the blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) signals in the thalamus and the spontaneous modulation of posterior alpha rhythms based on EEG-fMRI data acquired concurrently during an eyes-closed task-free condition. We observed both negative and positive correlations in the thalamus. The negative correlations were mostly seen within the visual thalamus, with a preference for the pulvinar over lateral geniculate nuclei. The positive correlations were found at the anterior and medial dorsal nuclei. Through functional connectivity analysis of the fMRI data, the pulvinar was found to be functionally associated with the same widespread cortical visual areas where the fMRI signals were negatively correlated with the posterior alpha modulation. In contrast, the dorsal nuclei were part of a distinct functional network that included brain stem, cingulate cortex and cerebellum. These observations are consistent with previous animal electrophysiology studies and the notion that the visual thalamus, and the pulvinar in particular, is intimately involved in the generation and spontaneous modulation of posterior alpha rhythms, facilitated by its reciprocal and widespread interaction with the cortical visual areas. We further postulate that the anterior and medial dorsal nuclei, being part of the ascending neuromodulatory system, may indirectly modulate cortical alpha rhythms by affecting vigilance and arousal levels.

Statistical feature extraction for artifact removal from concurrent fMRI-EEG recordings

We propose a set of algorithms for sequentially removing artifacts related to MRI gradient switching and cardiac pulsations from electroencephalography (EEG) data recorded during functional magnetic resonance imaging (fMRI). Special emphasis is directed upon the use of statistical metrics and methods for the extraction and selection of features that characterize gradient and pulse artifacts. To remove gradient artifacts, we use channel-wise filtering based on singular value decomposition (SVD). To remove pulse artifacts, we first decompose data into temporally independent components and then select a compact cluster of components that possess sustained high mutual information with the electrocardiogram (ECG). After the removal of these components, the time courses of remaining components are filtered by SVD to remove the temporal patterns phase-locked to the cardiac timing markers derived from the ECG. The filtered component time courses are then inversely transformed into multi-channel EEG time series free of pulse artifacts. Evaluation based on a large set of simultaneous EEG-fMRI data obtained during a variety of behavioral tasks, sensory stimulations and resting conditions showed excellent data quality and robust performance attainable with the proposed methods. These algorithms have been implemented as a Matlab-based toolbox made freely available for public access and research use.

Development and characterization of a dynamic lesion phantom for the quantitative evaluation of dynamic contrast-enhanced MRI

PURPOSE

To develop a dynamic lesion phantom that is capable of producing physiological kinetic curves representative of those seen in human dynamic contrast-enhanced MRI (DCE-MRI) data. The objective of this phantom is to provide a platform for the quantitative comparison of DCE-MRI protocols to aid in the standardization and optimization of breast DCE-MRI.

METHODS

The dynamic lesion consists of a hollow, plastic mold with inlet and outlet tubes to allow flow of a contrast agent solution through the lesion over time. Border shape of the lesion can be controlled using the lesion mold production method. The configuration of the inlet and outlet tubes was determined using fluid transfer simulations. The total fluid flow rate was determined using x-ray images of the lesion for four different flow rates (0.25, 0.5, 1.0, and 1.5 ml/s) to evaluate the resultant kinetic curve shape and homogeneity of the contrast agent distribution in the dynamic lesion. High spatial and temporal resolution x-ray measurements were used to estimate the true kinetic curve behavior in the dynamic lesion for benign and malignant example curves. DCE-MRI example data were acquired of the dynamic phantom using a clinical protocol.

RESULTS

The optimal inlet and outlet tube configuration for the lesion molds was two inlet molds separated by 30° and a single outlet tube directly between the two inlet tubes. X-ray measurements indicated that 1.0 ml/s was an appropriate total fluid flow rate and provided truth for comparison with MRI data of kinetic curves representative of benign and malignant lesions. DCE-MRI data demonstrated the ability of the phantom to produce realistic kinetic curves.

CONCLUSIONS

The authors have constructed a dynamic lesion phantom, demonstrated its ability to produce physiological kinetic curves, and provided estimations of its true kinetic curve behavior. This lesion phantom provides a tool for the quantitative evaluation of DCE-MRI protocols, which may lead to improved discrimination of breast cancer lesions.

Nonexponential T₂ decay in white matter

Visualizing myelin in human brain may help the study of diseases such as multiple sclerosis. Previous studies based on T(1) and T(2) relaxation contrast have suggested the presence of a distinct water pool that may report directly on local myelin content. Recent work indicates that T(2) contrast may offer particular advantages over T(1) and T(2) contrast, especially at high field. However, the complex mechanism underlying T(2) relaxation may render interpretation difficult. To address this issue, T(2) relaxation behavior in human brain was studied at 3 and 7 T. Multiple gradient echoes covering most of the decay curve were analyzed for deviations from mono-exponential behavior. The data confirm the previous finding of a distinct rapidly relaxing signal component (T(2) ∼ 6 ms), tentatively attributed to myelin water. However, in extension to previous findings, this rapidly relaxing component displayed a substantial resonance frequency shift, reaching 36 Hz in the corpus callosum at 7 T. The component's fractional amplitude and frequency shift appeared to depend on both field strength and fiber orientation, consistent with a mechanism originating from magnetic susceptibility effects. The findings suggest that T(2) contrast at high field may be uniquely sensitive to tissue myelin content and that proper interpretation will require modeling of susceptibility-induced resonance frequency shifts.

An anthropomorphic phantom for quantitative evaluation of breast MRI

PURPOSE

In this study, the authors aim to develop a physical, tissue-mimicking phantom for quantitative evaluation of breast MRI protocols. The objective of this phantom is to address the need for improved standardization in breast MRI and provide a platform for evaluating the influence of image protocol parameters on lesion detection and discrimination. Quantitative comparisons between patient and phantom image properties are presented.

METHODS

The phantom is constructed using a mixture of lard and egg whites, resulting in a random structure with separate adipose- and glandular-mimicking components. T1 and T2 relaxation times of the lard and egg components of the phantom were estimated at 1.5 T from inversion recovery and spin-echo scans, respectively, using maximum-likelihood methods. The image structure was examined quantitatively by calculating and comparing spatial covariance matrices of phantom and patient images. A static, enhancing lesion was introduced by creating a hollow mold with stereolithography and filling it with a gadolinium-doped water solution.

RESULTS

Measured phantom relaxation values fall within 2 standard errors of human values from the literature and are reasonably stable over 9 months of testing. Comparison of the covariance matrices of phantom and patient data demonstrates that the phantom and patient data have similar image structure. Their covariance matrices are the same to within error bars in the anterior-posterior direction and to within about two error bars in the right-left direction. The signal from the phantom's adipose-mimicking material can be suppressed using active fat-suppression protocols. A static, enhancing lesion can also be included with the ability to change morphology and contrast agent concentration.

CONCLUSIONS

The authors have constructed a phantom and demonstrated its ability to mimic human breast images in terms of key physical properties that are relevant to breast MRI. This phantom provides a platform for the optimization and standardization of breast MRI imaging protocols for lesion detection and characterization.

Negative BOLD-fMRI signals in large cerebral veins

Reductions in blood oxygenation level dependent (BOLD)-functional magnetic resonance imaging (fMRI) signals below baseline levels have been observed under several conditions as negative activation in task-activation studies or anticorrelation in resting-state experiments. Converging evidence suggests that negative BOLD signals (NBSs) can generally be explained by local reductions in neural activity. Here, we report on NBSs that accompany hemodynamic changes in regions devoid of neural tissue. The NBSs were investigated with high-resolution studies of the visual cortex (VC) at 7 T. Task-activation studies were performed to localize a task-positive area in the VC. During rest, robust negative correlation with the task-positive region was observed in focal regions near the ventricles and dispersed throughout the VC. Both positive and NBSs were dependent on behavioral condition. Comparison with high-resolution structural images showed that negatively correlated regions overlapped with larger pial and ependymal veins near sulcal and ventricular cerebrospinal fluid (CSF). Results from multiecho fMRI showed that NBSs were consistent with increases in local blood volume. These findings confirm theoretical predictions that tie neural activity to blood volume increases, which tend to counteract positive fMRI signal changes associated with increased blood oxygenation. This effect may be more salient in high-resolution studies, in which positive and NBS may be more often spatially distinct.

Magnetic susceptibility mapping of brain tissue in vivo using MRI phase data

Phase images in susceptibility-weighted MRI of brain provide excellent contrast. However, the phase is affected by tissue geometry and orientation relative to the main magnetic field (B(0)), and phase changes extend beyond areas of altered susceptibility. Magnetic susceptibility, on the other hand, is an intrinsic tissue property, closely reflecting tissue composition. Therefore, recently developed inverse Fourier-based methods were applied to calculate susceptibility maps from high-resolution phase images acquired at a single orientation at 7 T in the human brain (in vivo and fixed) and at 11.7 T in fixed marmoset brain. In susceptibility images, the contrast of cortical layers was more consistent than in phase images and was independent of the structures' orientation relative to B(0). The contrast of iron-rich deep-brain structures (red nucleus and substantia nigra) in susceptibility images agreed more closely with iron-dominated R(2) (*) images than the phase image contrast, which extended outside the structures. The mean susceptibility in these regions was significantly correlated with their estimated iron content. Susceptibility maps calculated using this method overcome the orientation-dependence and non-locality of phase image contrast and seem to reflect underlying tissue composition. Susceptibility images should be easier to interpret than phase images and could improve our understanding of the sources of susceptibility contrast.

Large-scale spontaneous fluctuations and correlations in brain electrical activity observed with magnetoencephalography

Knowledge about the intrinsic functional architecture of the human brain has been greatly expanded by the extensive use of resting-state functional magnetic resonance imaging (fMRI). However, the neurophysiological correlates and origins of spontaneous fMRI signal changes remain poorly understood. In the present study, we characterized the power modulations of spontaneous magnetoencephalography (MEG) rhythms recorded from human subjects during wakeful rest (with eyes open and eyes closed) and light sleep. Through spectral, correlation and coherence analyses, we found that resting-state MEG rhythms demonstrated ultraslow (<0.1 Hz) spontaneous power modulations that synchronized over a large spatial distance, especially between bilaterally homologous regions in opposite hemispheres. These observations are in line with the known spatio-temporal properties of spontaneous fMRI signals, and further suggest that the coherent power modulation of spontaneous rhythmic activity reflects the electrophysiological signature of the large-scale functional networks previously observed with fMRI in the resting brain.

Prospective head-movement correction for high-resolution MRI using an in-bore optical tracking system

In MRI of the human brain, subject motion is a major cause of magnetic resonance image quality degradation. To compensate for the effects of head motion during data acquisition, an in-bore optical motion tracking system is proposed. The system comprises two MR-compatible infrared cameras that are fixed on a holder right above and in front of the head coil. The resulting close proximity of the cameras to the object allows precise tracking of its movement. During image acquisition, the MRI scanner uses this tracking information to prospectively compensate for head motion by adjusting the gradient field direction and radio frequency (RF) phases and frequencies. Experiments performed on subjects demonstrate robust system performance with translation and rotation accuracies of 0.1 mm and 0.15 degrees, respectively.

Hemodynamic nonlinearities affect BOLD fMRI response timing and amplitude

The interpretation of functional magnetic resonance imaging (fMRI) studies based on blood oxygen-level dependent (BOLD) contrast generally relies on the assumption of a linear relationship between evoked neuronal activity and fMRI response. While nonlinearities in this relationship have been suggested by a number of studies, it remains unclear to what extent they relate to the neurovascular response and are therefore inherent to BOLD fMRI. Full characterization of potential vascular nonlinearities is required for accurate inferences about the neuronal system under study. To investigate the extent of vascular nonlinearities, evoked activity was studied in humans with BOLD fMRI (n=28) and magnetoencephalography (MEG) (n=5). Brief (600-800 ms) rapidly repeated (1 Hz) visual stimuli were delivered using a stimulation paradigm that minimized neuronal nonlinearities. Nevertheless, BOLD fMRI experiments showed substantial remaining nonlinearities. The smallest stimulus separation (200-400 ms) resulted in significant response broadening (15-20% amplitude decrease; 10-12% latency increase; 6-14% duration increase) with respect to a linear prediction. The substantial slowing and widening of the response in the presence of preceding stimuli suggest a vascular rather than neuronal origin to the observed nonlinearity. This was confirmed by the MEG data, which showed no significant neuro-electric nonlinear interactions between stimuli as little as 200 ms apart. The presence of substantial vascular nonlinearities has important implications for rapid event-related studies by fMRI and other imaging modalities that infer neuronal activity from hemodynamic parameters.

Sources of functional magnetic resonance imaging signal fluctuations in the human brain at rest: a 7 T study

Signal fluctuations in functional magnetic resonance imaging (fMRI) can result from a number of sources that may have a neuronal, physiologic or instrumental origin. To determine the relative contribution of these sources, we recorded physiological (respiration and cardiac) signals simultaneously with fMRI in human volunteers at rest with their eyes closed. State-of-the-art technology was used including high magnetic field (7 T), a multichannel detector array and high-resolution (3 mm(3)) echo-planar imaging. We investigated the relative contribution of thermal noise and other sources of variance to the observed fMRI signal fluctuations both in the visual cortex and in the whole brain gray matter. The following sources of variance were evaluated separately: low-frequency drifts due to scanner instability, effects correlated with respiratory and cardiac cycles, effects due to variability in the respiratory flow rate and cardiac rate, and other sources, tentatively attributed to spontaneous neuronal activity. We found that low-frequency drifts are the most significant source of fMRI signal fluctuations (3.0% signal change in the visual cortex, TE=32 ms), followed by spontaneous neuronal activity (2.9%), thermal noise (2.1%), effects due to variability in physiological rates (respiration 0.9%, heartbeat 0.9%), and correlated with physiological cycles (0.6%). We suggest the selection and use of four lagged physiological noise regressors as an effective model to explain the variance related to fluctuations in the rates of respiration volume change and cardiac pulsation. Our results also indicate that, compared to the whole brain gray matter, the visual cortex has higher sensitivity to changes in both the rate of respiration and the spontaneous resting-state activity. Under the conditions of this study, spontaneous neuronal activity is one of the major contributors to the measured fMRI signal fluctuations, increasing almost twofold relative to earlier experiments under similar conditions at 3 T.

From phonemes to articulatory codes: an fMRI study of the role of Broca's area in speech production

We used event-related functional magnetic resonance imaging to investigate the neuroanatomical substrates of phonetic encoding and the generation of articulatory codes from phonological representations. Our focus was on the role of the left inferior frontal gyrus (LIFG) and in particular whether the LIFG plays a role in sublexical phonological processing such as syllabification or whether it is directly involved in phonetic encoding and the generation of articulatory codes. To answer this question, we contrasted the brain activation patterns elicited by pseudowords with high- or low-sublexical frequency components, which we expected would reveal areas related to the generation of articulatory codes but not areas related to phonological encoding. We found significant activation of a premotor network consisting of the dorsal precentral gyrus, the inferior frontal gyrus bilaterally, and the supplementary motor area for low- versus high-sublexical frequency pseudowords. Based on our hypothesis, we concluded that these areas and in particular the LIFG are involved in phonetic and not phonological encoding. We further discuss our findings with respect to the mechanisms of phonetic encoding and provide evidence in support of a functional segregation of the posterior part of Broca's area, the pars opercularis.

Modulation of spontaneous fMRI activity in human visual cortex by behavioral state

The phenomenon of spontaneous fMRI activity is increasingly being exploited to investigate the connectivity of functional networks in human brain with high spatial-resolution. Although mounting evidence points towards a neuronal contribution to this activity, its functional role and dependence on behavioral state remain unclear. In this work, we used BOLD fMRI at 7 T to study the modulation of spontaneous activity in occipital areas by various behavioral conditions, including resting with eyes closed, eyes open with visual fixation, and eyes open with fixation and focal visual stimulation. Spontaneous activity was separated from evoked activity and from signal fluctuations related to cardiac and respiratory cycles. We found that spontaneous activity in visual areas was substantially reduced (amplitude (44%) and coherence (25%)) with the fixation conditions relative to the eyes-closed condition. No significant further modulation was observed when the visual stimulus was added. The observed dependence on behavioral condition suggests that part of spontaneous fMRI signal fluctuations represents neuronal activity. Possible mechanisms for the modulation of spontaneous activity by behavioral state are discussed. The observed linear superposition of spontaneous fMRI activity with focal evoked activity related to visual processing has important implications for fMRI studies, which ideally should take into account the effect of spontaneous activity to properly define brain activations during task conditions.

Low frequency BOLD fluctuations during resting wakefulness and light sleep: a simultaneous EEG-fMRI study

Recent blood oxygenation level dependent functional MRI (BOLD fMRI) studies of the human brain have shown that in the absence of external stimuli, activity persists in the form of distinct patterns of temporally correlated signal fluctuations. In this work, we investigated the spontaneous BOLD signal fluctuations during states of reduced consciousness such as drowsiness and sleep. For this purpose, we performed BOLD fMRI on normal subjects during varying levels of consciousness, from resting wakefulness to light (non-slow wave) sleep. Depth of sleep was determined based on concurrently acquired EEG data. During light sleep, significant increases in the fluctuation level of the BOLD signal were observed in several cortical areas, among which visual cortex was the most significant. Correlations among brain regions involved with the default-mode network persisted during light sleep. These results suggest that activity in areas such as the default-mode network and primary sensory cortex, as measured from BOLD fMRI fluctuations, does not require a level of consciousness typical of wakefulness.

Metabolic origin of BOLD signal fluctuations in the absence of stimuli

Blood oxygen level-dependent (BOLD) functional magnetic resonance imaging studies have shown the existence of ongoing blood flow fluctuations in the absence of stimuli. Although this so-called 'resting-state activity' appears to be correlated across brain regions with apparent functional relationship, its origin might be predominantly vascular and not directly representing neuronal signaling. To investigate this, we simultaneously measured BOLD and perfusion signals on healthy human subjects (n=11) and used their ratio (BOLD/perfusion ratio or BPR) as an indicator of metabolic demand. BPR during rest and sleep was compared with that during a visual task (VT) and a breath-holding task (BH), which are challenges with substantial and little metabolic involvement, respectively. Within the visual cortex, BPR was 3.76+/-1.23 during BH, which was significantly higher than during the VT (1.76+/-0.27) and rest (1.56+/-0.41). Meanwhile, BPR values during VT and rest were not significantly different, suggesting a similar metabolic involvement. Eight subjects showed stage 1 and 2 sleep, during which temporally correlated BOLD and perfusion activity continued. In these subjects, there was no significant difference in BPR between the sleep and waking conditions (1.79+/-0.54 and 1.66+/-0.67, respectively), but both were lower than the BPR during BH. These data suggest that resting-state activity, at least in part, represents a metabolic process.

Reducing correlated noise in fMRI data

The sensitivity of functional MRI (fMRI) in detecting neuronal activation is dependent on the relative levels of signal and noise in the time-series data. The temporal noise level within a single voxel is generally substantially higher than the intrinsic NMR (thermal) noise, and the noise is often correlated between voxels. This work introduces and evaluates a method that allows fMRI sensitivity improvement by reduction of these correlated noise sources. The method allows model-free estimation of the correlated noise from brain regions not activated by the functional paradigm using a short (1-2 min) reference scan. A single regressor representing this noise-source estimate is added to the design matrix used in the data analysis. Results obtained from five volunteers show an average t-score improvement of 11.3% and a 24.2% increase in the size of the activated area.

Receive coil arrays and parallel imaging for functional magnetic resonance imaging of the human brain

The use of multiple small receive coils has several advantages over a single larger (e.g. birdcage) coil. With an array of small receive coils, image signal-to-noise ratio (SNR) can be optimized throughout the field-of-view and the image acquisition process can be accelerated by use of parallel imaging (PI) techniques. In an accelerated PI experiment, data are undersampled during acquisition. Subsequently, artifact-free images are reconstructed based on the independently acquired signals from the elements of the receive coil array. PI techniques have recently been applied to functional MRI (fMRI) experiments of the human brain in order to improve the performance of commonly used single-shot techniques like echo-planar imaging (EPI). Potential benefits of PI-fMRI include the reduction of geometrical distortions due to off-resonance signals, the reduction of signal-loss in areas with substantial signal inhomogeneity, increases of the spatial and temporal resolution of the fMRI experiment and reduction of gradient acoustic noise. Although the loss in SNR, inherent to PI, can severely compromise MRI image quality, the effect on fMRI quality, which is governed by the temporal stability of the signal, is often not as severe. On the other hand, PI's potential in mitigating the often severe image artifacts present in single-shot fMRI render it an important tool, in particular with the recent surge in high field MRI applications.

Low-frequency fluctuations in the cardiac rate as a source of variance in the resting-state fMRI BOLD signal

Heart rate fluctuations occur in the low-frequency range (<0.1 Hz) probed in functional magnetic resonance imaging (fMRI) studies of resting-state functional connectivity and most fMRI block paradigms and may be related to low-frequency blood-oxygenation-level-dependent (BOLD) signal fluctuations. To investigate this hypothesis, temporal correlations between cardiac rate and resting-state fMRI signal timecourses were assessed at 3 T. Resting-state BOLD fMRI and accompanying physiological data were acquired and analyzed using cross-correlation and regression. Time-shifted cardiac rate timecourses were included as regressors in addition to established physiological regressors (RETROICOR (Glover, G.H., Li, T.Q., Ress, D., 2000. Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn Reson Med 44, 162-167) and respiration volume per unit time (Birn, R.M., Diamond, J.B., Smith, M.A., Bandettini, P.A., 2006b. Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. NeuroImage 31, 1536-1548). Significant correlations between the cardiac rate and BOLD signal timecourses were revealed, particularly negative correlations in gray matter at time shifts of 6-12 s and positive correlations at time shifts of 30-42 s (TR=6 s). Regressors consisting of cardiac rate timecourses shifted by delays of between 0 and 24 s explained an additional 1% of the BOLD signal variance on average over the whole brain across 9 subjects, a similar additional variance to that explained by respiration volume per unit time and RETROICOR regressors, even when used in combination with these other physiological regressors. This suggests that including such time-shifted cardiac rate regressors will be beneficial for explaining physiological noise variance and will thereby improve the statistical power in future task-based and resting-state fMRI studies.

High-field MRI of brain cortical substructure based on signal phase

The ability to detect brain anatomy and pathophysiology with MRI is limited by the contrast-to-noise ratio (CNR), which depends on the contrast mechanism used and the spatial resolution. In this work, we show that in MRI of the human brain, large improvements in contrast to noise in high-resolution images are possible by exploiting the MRI signal phase at high magnetic field strength. Using gradient-echo MRI at 7.0 tesla and a multichannel detector, a nominal voxel size of 0.24 x 0.24 x 1.0 mm3 (58 nl) was achieved. At this resolution, a strong phase contrast was observed both between as well as within gray matter (GM) and white matter (WM). In gradient-echo phase images obtained on normal volunteers at this high resolution, the CNR between GM and WM ranged from 3:1 to 20:1 over the cortex. This CNR is an almost 10-fold improvement over conventional MRI techniques that do not use image phase, and it is an approximately 100-fold improvement when including the gains in resolution from high-field and multichannel detection. Within WM, phase contrast appeared to be associated with the major fiber bundles, whereas contrast within GM was suggestive of the underlying layer structure. The observed phase contrast is attributed to local variations in magnetic susceptibility, which, at least in part, appeared to originate from iron stores. The ability to detect cortical substructure from MRI phase contrast at high field is expected to greatly enhance the study of human brain anatomy in vivo.

An adaptive filter for suppression of cardiac and respiratory noise in MRI time series data

The quality of MRI time series data, which allows the study of dynamic processes, is often affected by confounding sources of signal fluctuation, including the cardiac and respiratory cycle. An adaptive filter is described, reducing these signal fluctuations as long as they are repetitive and their timing is known. The filter, applied in image domain, does not require temporal oversampling of the artifact-related fluctuations. Performance is demonstrated for suppression of cardiac and respiratory artifacts in 10-minute brain scans on 6 normal volunteers. Experimental parameters resemble a typical fMRI experiment (17 slices; 1700 ms TR). A second dataset was acquired at a rate well above the Nyquist frequency for both cardiac and respiratory cycle (single slice; 100 ms TR), allowing identification of artifacts specific to the cardiac and respiratory cycles, aiding assessment of filtering performance. Results show significant reduction in temporal standard deviation (SD(t)) in all subjects. For all 6 datasets with 1700 ms TR combined, the filtering method resulted in an average reduction in SD(t) of 9.2% in 2046 voxels substantially affected by respiratory artifacts, and 12.5% for the 864 voxels containing substantial cardiac artifacts. The maximal SD(t) reduction achieved was 52.7% for respiratory and 55.3% for cardiac filtering. Performance was found to be at least equivalent to the previously published RETROICOR method. Furthermore, the interaction between the filter and fMRI activity detection was investigated using Monte Carlo simulations, demonstrating that filtering algorithms introduce a systematic error in the detected BOLD-related signal change if applied sequentially. It is demonstrated that this can be overcome by combining physiological artifact filtering and detection of BOLD-related signal changes simultaneously. Visual fMRI data from 6 volunteers were analyzed with and without the filter proposed here. Inclusion of the cardio-respiratory regressors in the design matrix yielded a 4.6% t-score increase and 4.0% increase in the number of significantly activated voxels.