Purpose
Slice‐wise shimming can improve field homogeneity, but suffers from large noise propagation in the shim calculation. Here, we propose a robust shim current optimization for higher‐order ...dynamic shim updating, based on Tikhonov regularization with a variable regularization parameter, λ.
Theory and Methods
λ was selected for each slice separately in a fully automatic procedure based on a combination of boundary constraints and an L‐curve search algorithm. Shimming performance was evaluated for second order slice‐wise shimming of the brain at 7T, by simulation on a database of field maps from 143 subjects, and by direct measurements in 8 subjects.
Results
Simulations yielded on average 36% reduction in the shim current norm for just 0.4 Hz increase in residual field SD as compared to unconstrained unregularized optimization. In vivo results yielded on average 34.0 Hz residual field SD as compared to 34.3 Hz with a constrained unregularized optimization, while simultaneously reducing the shim current norm to 2.8 A from 3.9 A. The proposed regularization also reduced the average step in the shim current between slices.
Conclusion
Slice‐wise variable Tikhonov regularization yielded reduced current norm and current steps to a negligible cost in field inhomogeneity. The method holds promise to increase the robustness, and thereby the utility, of higher‐order shim updating.
Magnetic resonance imaging and spectroscopy of the spinal cord stand to benefit greatly from the increased signal-to-noise ratio of ultra-high field. However, ultra-high field also poses considerable ...technical challenges, especially related to static and dynamic B
fields. Breathing causes the field to fluctuate with the respiratory cycle, giving rise to artifacts such as ghosting and apparent motion in images. We here investigated the spatial and temporal characteristics of breathing-induced B
fields in the cervical spinal cord at 7T. We analyzed the magnitude and spatial profile of breathing-induced fields during breath-holds in an expired and inspired breathing state. We also measured the temporal field evolution during free breathing by acquiring a time series of fast phase images, and a principal component analysis was performed on the measured field evolution. In all subjects, the field shift was largest around the vertebral level of C7 and lowest at the top of the spinal cord. At C7, we measured peak-to-peak field fluctuations of 36 Hz on average during normal free breathing; increasing to on average 113 Hz during deep breathing. The first principal component could explain more than 90% of the field variations along the foot-head axis inside the spinal cord in all subjects. We further implemented a proof-of-principle shim correction, demonstrating the feasibility of using the shim system to compensate for the breathing-induced fields inside the spinal cord. Effective correction strategies will be crucial to unlock the full potential of ultra-high field for spinal cord imaging.
Studies of human primary somatosensory cortex (S1) have placed a strong emphasis on the cortical representation of the hand and the propensity for plasticity therein. Despite many reports of group ...differences and experience-dependent changes in cortical digit somatotopy, relatively little work has considered the variability of these maps across individuals and to what extent this detailed functional architecture is dynamic over time. With the advent of 7 T fMRI, it is increasingly feasible to map such detailed organization noninvasively in individual human participants. Here, we extend the ability of ultra-high-field imaging beyond a technological proof of principle to investigate the intersubject variability of digit somatotopy across participants and the stability of this organization across a range of intervals. Using a well validated phase-encoding paradigm and an active task, we demonstrate the presence of highly reproducible maps of individual digits in S1, sharply contrasted by a striking degree of intersubject variability in the shape, extent, and relative position of individual digit representations. Our results demonstrate the presence of very stable fine-grain somatotopy of the digits in human S1 and raise the issue of population variability in such detailed functional architecture of the human brain. These findings have implications for the study of detailed sensorimotor plasticity in the context of both learning and pathological dysfunction. The simple task and 10 min scan required to derive these maps also raises the potential for this paradigm as a tool in the clinical setting.
We applied ultra-high-resolution fMRI at 7 T to map sensory digit representations in the human primary somatosensory cortex (S1) at the level of individual participants across multiple time points. The resulting fine-grain maps of individual digits in S1 reveal the stability in this fine-grain functional organization over time, contrasted with the variability in these maps across individuals.
Purpose
Spinal cord MRI at ultrahigh field is hampered by time‐varying magnetic fields associated with the breathing cycle, giving rise to ghosting artifacts in multi‐shot acquisitions. Here, we ...suggest a correction approach based on linking the signal from a respiratory bellows to field changes inside the spinal cord. The information is used to correct the data at the image reconstruction level.
Methods
The correction was demonstrated in the context of multi‐shot T2*‐weighted imaging of the cervical spinal cord at 7T. A respiratory trace was acquired during a high‐resolution multi‐echo gradient‐echo sequence, used for structural imaging and quantitative T2* mapping, and a multi‐shot EPI time series, as would be suitable for fMRI. The coupling between the trace and the breathing‐induced fields was determined by a short calibration scan in each individual. Images were reconstructed with and without trace‐based correction.
Results
In the multi‐echo acquisition, breathing‐induced fields caused severe ghosting in images with long TE, which led to a systematic underestimation of T2* in the spinal cord. The trace‐based correction reduced the ghosting and increased the estimated T2* values. Breathing‐related ghosting was also observed in the multi‐shot EPI images. The correction largely removed the ghosting, thereby improving the temporal signal‐to‐noise ratio of the time series.
Conclusions
Trace‐based retrospective correction of breathing‐induced field variations can reduce ghosting and improve quantitative metrics in multi‐shot structural and functional T2*‐weighted imaging of the spinal cord. The method is straightforward to implement and does not rely on sequence modifications or additional hardware beyond a respiratory bellows.
Limited understanding of infant pain has led to its lack of recognition in clinical practice. While the network of brain regions that encode the affective and sensory aspects of adult pain are well ...described, the brain structures involved in infant nociceptive processing are completely unknown, meaning we cannot infer anything about the nature of the infant pain experience. Using fMRI we identified the network of brain regions that are active following acute noxious stimulation in newborn infants, and compared the activity to that observed in adults. Significant infant brain activity was observed in 18 of the 20 active adult brain regions but not in the infant amygdala or orbitofrontal cortex. Brain regions that encode sensory and affective components of pain are active in infants, suggesting that the infant pain experience closely resembles that seen in adults. This highlights the importance of developing effective pain management strategies in this vulnerable population.
Current clinical and experimental literature strongly supports the phenomenon of reduced pain perception whilst attention is distracted away from noxious stimuli. This study used functional MRI to ...elucidate the underlying neural systems and mechanisms involved. An analogue of the Stroop task, the counting Stroop, was used as a cognitive distraction task whilst subjects received intermittent painful thermal stimuli. Pain intensity scores were significantly reduced when subjects took part in the more cognitively demanding interference task of the counting Stroop than in the less demanding neutral task. When subjects were distracted during painful stimulation, brain areas associated with the affective division of the anterior cingulate cortex (ACC) and orbitofrontal regions showed increased activation. In contrast, many areas of the pain matrix (i.e. thalamus, insula, cognitive division of the ACC) displayed reduced activation, supporting the behavioural results of reduced pain perception.
Experience-dependent reorganisation of functional maps in the cerebral cortex is well described in the primary sensory cortices. However, there is relatively little evidence for such cortical ...reorganisation over the short-term. Using human somatosensory cortex as a model, we investigated the effects of a 24 hr gluing manipulation in which the right index and right middle fingers (digits 2 and 3) were adjoined with surgical glue. Somatotopic representations, assessed with two 7 tesla fMRI protocols, revealed rapid off-target reorganisation in the non-manipulated fingers following gluing, with the representation of the ring finger (digit 4) shifted towards the little finger (digit 5) and away from the middle finger (digit 3). These shifts were also evident in two behavioural tasks conducted in an independent cohort, showing reduced sensitivity for discriminating the temporal order of stimuli to the ring and little fingers, and increased substitution errors across this pair on a speeded reaction time task.
We present the reliability of ultra-high field T2* MRI at 7T, as part of the UK7T Network's “Travelling Heads” study. T2*-weighted MRI images can be processed to produce quantitative susceptibility ...maps (QSM) and R2* maps. These reflect iron and myelin concentrations, which are altered in many pathophysiological processes. The relaxation parameters of human brain tissue are such that R2* mapping and QSM show particularly strong gains in contrast-to-noise ratio at ultra-high field (7T) vs clinical field strengths (1.5–3T). We aimed to determine the inter-subject and inter-site reproducibility of QSM and R2* mapping at 7T, in readiness for future multi-site clinical studies.
Ten healthy volunteers were scanned with harmonised single- and multi-echo T2*-weighted gradient echo pulse sequences. Participants were scanned five times at each “home” site and once at each of four other sites. The five sites had 1× Philips, 2× Siemens Magnetom, and 2× Siemens Terra scanners. QSM and R2* maps were computed with the Multi-Scale Dipole Inversion (MSDI) algorithm (https://github.com/fil-physics/Publication-Code). Results were assessed in relevant subcortical and cortical regions of interest (ROIs) defined manually or by the MNI152 standard space.
Mean susceptibility (χ) and R2* values agreed broadly with literature values in all ROIs. The inter-site within-subject standard deviation was 0.001–0.005 ppm (χ) and 0.0005–0.001 ms−1 (R2*). For χ this is 2.1–4.8 fold better than 3T reports, and 1.1–3.4 fold better for R2*. The median ICC from within- and cross-site R2* data was 0.98 and 0.91, respectively. Multi-echo QSM had greater variability vs single-echo QSM especially in areas with large B0 inhomogeneity such as the inferior frontal cortex. Across sites, R2* values were more consistent than QSM in subcortical structures due to differences in B0-shimming. On a between-subject level, our measured χ and R2* cross-site variance is comparable to within-site variance in the literature, suggesting that it is reasonable to pool data across sites using our harmonised protocol.
The harmonized UK7T protocol and pipeline delivers on average a 3-fold improvement in the coefficient of reproducibility for QSM and R2* at 7T compared to previous reports of multi-site reproducibility at 3T. These protocols are ready for use in multi-site clinical studies at 7T.