Model-based navigation of transcranial focused ultrasound neuromodulation in humans: Application to targeting the amygdala and thalamus.

BACKGROUND: Transcranial focused ultrasound (tFUS) neuromodulation has shown promise in animals but is challenging to translate to humans because of the thicker skull that heavily scatters ultrasound waves. OBJECTIVE: We develop and disseminate a model-based navigation (MBN) tool for acoustic dose delivery in the presence of skull aberrations that is easy to use by non-specialists. METHODS: We pre-compute acoustic beams for thousands of virtual transducer locations on the scalp of the subject under study. We use the hybrid angular spectrum solver mSOUND, which runs in ∼4 s per solve per CPU yielding pre-computation times under 1 h for scalp meshes with up to 4000 faces and a parallelization factor of 5. We combine this pre-computed set of beam solutions with optical tracking, thus allowing real-time display of the tFUS beam as the operator freely navigates the transducer around the subject' scalp. We assess the impact of MBN versus line-of-sight targeting (LOST) positioning in simulations of 13 subjects. RESULTS: Our navigation tool has a display refresh rate of ∼10 Hz. In our simulations, MBN increased the acoustic dose in the thalamus and amygdala by 8-67 % compared to LOST and avoided complete target misses that affected 10-20 % of LOST cases. MBN also yielded a lower variability of the deposited dose across subjects than LOST. CONCLUSIONS: MBN may yield greater and more consistent (less variable) ultrasound dose deposition than transducer placement with line-of-sight targeting, and thus could become a helpful tool to improve the efficacy of tFUS neuromodulation.

Robust deep learning estimation of cortical bone porosity from MR T1-weighted images for individualized transcranial focused ultrasound planning.

Objective: Transcranial focused ultrasound (tFUS) is an emerging neuromodulation approach that has been demonstrated in animals but is difficult to translate to humans because of acoustic attenuation and scattering in the skull. Optimal dose delivery requires subject-specific skull porosity estimates which has traditionally been done using CT. We propose a deep learning (DL) estimation of skull porosity from T1-weighted MRI images which removes the need for radiation-inducing CT scans. Approach: We evaluate the impact of different DL approaches, including network architecture, input size and dimensionality, multichannel inputs, data augmentation, and loss functions. We also propose back-propagation in the mask (BIM), a method whereby only voxels inside the skull mask contribute to training. We evaluate the robustness of the best model to input image noise and MRI acquisition parameters and propagate porosity estimation errors in thousands of beam propagation scenarios. Main results: Our best performing model is a cGAN with a ResNet-9 generator with 3D 64×64×64 inputs trained with L1 and L2 losses. The model achieved a mean absolute error of 6.9% in the test set, compared to 9.5% with the pseudo-CT of Izquierdo et al. (38% improvement) and 9.4% with the generic pixel-to-pixel image translation cGAN pix2pix (36% improvement). Acoustic dose distributions in the thalamus were more accurate with our approach than with the pseudo-CT approach of both Burgos et al. and Izquierdo et al, resulting in near-optimal treatment planning and dose estimation at all frequencies compared to CT (reference). Significance: Our DL approach porosity estimates with ~7% error, is robust to input image noise and MRI acquisition parameters (sequence, coils, field strength) and yields near-optimal treatment planning and dose estimates for both central (thalamus) and lateral brain targets (amygdala) in the 200-1000 kHz frequency range.

Multiphoton parallel transmission (MP-pTx): Pulse design methods and numerical validation.

PURPOSE: We propose and evaluate multiphoton parallel transmission (MP-pTx) to mitigate flip angle inhomogeneities in high-field MRI. MP-pTx is an excitation method that utilizes a single, conventional birdcage coil supplemented with low-frequency (kHz) irradiation from a multichannel shim array and/or gradient channels. SAR analysis is simplified to that of a conventional birdcage coil, because only the radiofrequency (RF) field from the birdcage coil produces significant SAR. METHODS: MP-pTx employs an off-resonance RF pulse from a conventional birdcage coil supplemented with oscillating z $$ z $$ -directed fields from a multichannel shim array and/or the gradient coils. We simulate the ability of MP-pTx to create uniform nonselective brain excitations at 7 T using realistic B 1 + $$ {\mathrm{B}}_1^{+} $$ and Δ B 0 $$ \Delta {\mathrm{B}}_0 $$ field maps. The RF, shim array, and gradient waveform's amplitudes and phases are optimized using a genetic algorithm followed by sequential quadratic programming. RESULTS: A 1 ms MP-pTx excitation using a 32-channel shim array with current constrained to less than 50 Amp-turns reduced the transverse magnetization's normalized root-mean-squared error from 29% for a conventional birdcage excitation to 6.6% and was nearly 40% better than a 1 ms birdcage coil 5 kT-point excitation with optimized kT-point locations and comparable pulse power. CONCLUSION: The MP-pTx method resembles conventional pTx in its goals and approach but replaces the parallel RF channels with cheaper, low-frequency shim channels. The method mitigates high-field flip angle inhomogeneities to a level better than 3 T CP-mode and comparable to 7 T pTx while retaining the straightforward SAR characteristics of conventional birdcage excitations, as low-frequency shim array fields produce negligible SAR.

Experimental validation of a PNS-optimized whole-body gradient coil.

PURPOSE: Peripheral nerve stimulation (PNS) limits the usability of state-of-the-art whole-body and head-only MRI gradient coils. We used detailed electromagnetic and neurodynamic modeling to set an explicit PNS constraint during the design of a whole-body gradient coil and constructed it to compare the predicted and experimentally measured PNS thresholds to those of a matched design without PNS constraints. METHODS: We designed, constructed, and tested two actively shielded whole-body Y-axis gradient coil winding patterns: YG1 is a conventional symmetric design without PNS-optimization, whereas YG2's design used an additional constraint on the allowable PNS threshold in the head-imaging landmark, yielding an asymmetric winding pattern. We measured PNS thresholds in 18 healthy subjects at five landmark positions (head, cardiac, abdominal, pelvic, and knee). RESULTS: The PNS-optimized design YG2 achieved 46% higher average experimental thresholds for a head-imaging landmark than YG1 while incurring a 15% inductance penalty. For cardiac, pelvic, and knee imaging landmarks, the PNS thresholds increased between +22% and +35%. For abdominal imaging, PNS thresholds did not change significantly between YG1 and YG2 (-3.6%). The agreement between predicted and experimental PNS thresholds was within 11.4% normalized root mean square error for both coils and all landmarks. The PNS model also produced plausible predictions of the stimulation sites when compared to the sites of perception reported by the subjects. CONCLUSION: The PNS-optimization improved the PNS thresholds for the target scan landmark as well as most other studied landmarks, potentially yielding a significant improvement in image encoding performance that can be safely used in humans.

Influence of peripheral axon geometry and local anatomy on magnetostimulation chronaxie.

Objective.Rapid switching of magnetic resonance imaging (MRI) gradient fields induces electric fields that can cause peripheral nerve stimulation (PNS) and so accurate characterization of PNS is required to maintain patient safety and comfort while maximizing MRI performance. The minimum magnetic gradient amplitude that causes stimulation, the PNS threshold, depends on intrinsic axon properties and the spatial and temporal properties of the induced electric field. The PNS strength-duration curve is widely used to characterize simulation thresholds for periodic waveforms and is parameterized by the chronaxie and rheobase. Safety limits to avoid unwanted PNS in MRI rely on a single chronaxie value to characterize the response of all nerves. However, experimental magnetostimulation peripheral nerve chronaxie values vary by an order of magnitude. Given the diverse range of chronaxies observed and the importance of this number in MRI safety models, we seek a deeper understanding of the mechanisms contributing to chronaxie variability.Approach.We use a coupled electromagnetic-neurodynamic PNS model to assess geometric sources of chronaxie variability. We study the impact of the position of the stimulating magnetic field coil relative to the body, along with the effect of local anatomical features and nerve trajectories on the driving function and the resulting chronaxie.Main results.We find realistic variation of local axon and tissue geometry can modulate a given axon's chronaxie by up to two-fold. Our results identify the temporal rate of charge redistribution as the underlying determinant of the chronaxie.Significance.This charge distribution is a function of both intrinsic axon properties and the spatial stimulus along the nerve; thus, examination of the local tissue topology, which shapes the electric fields, as well as the nerve trajectory, are critical for better understanding chronaxie variations and defining more biologically informed MRI safety guidelines.

Transcranial focused ultrasound of the amygdala modulates fear network activation and connectivity.

BACKGROUND: Current noninvasive brain stimulation methods are incapable of directly modulating subcortical brain regions critically involved in psychiatric disorders. Transcranial Focused Ultrasound (tFUS) is a newer form of noninvasive stimulation that could modulate the amygdala, a subcortical region implicated in fear. OBJECTIVE: We investigated the effects of active and sham tFUS of the amygdala on fear circuit activation, skin conductance responses (SCR), and self-reported anxiety during a fear-inducing task. We also investigated amygdala tFUS' effects on amygdala-fear circuit resting-state functional connectivity. METHODS: Thirty healthy individuals were randomized in this double-blinded study to active or sham tFUS of the left amygdala. We collected fMRI scans, SCR, and self-reported anxiety during a fear-inducing task (participants viewed red or green circles which indicated the risk of receiving an aversive stimulus), as well as resting-state scans, before and after tFUS. RESULTS: Compared to sham tFUS, active tFUS was associated with decreased (pre to post tFUS) blood-oxygen-level-dependent fMRI activation in the amygdala (F(1,25) = 4.86, p = 0.04, η2 = 0.16) during the fear task, and lower hippocampal (F(1,27) = 4.41, p = 0.05, η2 = 0.14), and dorsal anterior cingulate cortex (F(1,27) = 6.26, p = 0.02; η2 = 0.19) activation during the post tFUS fear task. The decrease in amygdala activation was correlated with decreased subjective anxiety (r = 0.62, p = 0.03). There was no group effect in SCR changes from pre to post tFUS (F(1,23) = 0.85, p = 0.37). The active tFUS group also showed decreased amygdala-insula (F(1,28) = 4.98, p = 0.03) and amygdala-hippocampal (F(1,28) = 7.14, p = 0.01) rsFC, and increased amygdala-ventromedial prefrontal cortex (F(1,28) = 3.52, p = 0.05) resting-state functional connectivity. CONCLUSIONS: tFUS can change functional connectivity and brain region activation associated with decreased anxiety. Future studies should investigate tFUS' therapeutic potential for individuals with clinical levels of anxiety.

Comparison of tight-fitting 7T parallel-transmit head array designs using excitation uniformity and local specific absorption rate metrics.

PURPOSE: We model the performance of parallel transmission (pTx) arrays with 8, 16, 24, and 32 channels and varying loop sizes built on a close-fitting helmet for brain imaging at 7 T and compare their local specific absorption rate (SAR) and flip-angle performances to that of birdcage coil (used as a baseline) and cylindrical 8-channel and 16-channel pTx coils (single-row and dual-row). METHODS: We use the co-simulation approach along with MATLAB scripting for batch-mode simulation of the coils. For each coil, we extracted B1 + maps and SAR matrices, which we compressed using the virtual observation points algorithm, and designed slice-selective RF shimming pTx pulses with multiple local SAR and peak power constraints to generate L-curves in the transverse, coronal, and sagittal orientations. RESULTS: Helmet designs outperformed cylindrical pTx arrays at a constant number of channels in the flip-angle uniformity at a constant local SAR metric: up to 29% for 8-channel arrays, and up to 34% for 16-channel arrays, depending on the slice orientation. For all helmet arrays, increasing the loop diameter led to better local SAR versus flip-angle uniformity tradeoffs, although this effect was more pronounced for the 8-channel and 16-channel systems than the 24-channel and 32-channel systems, as the former have more limited degrees of freedom and therefore benefit more from loop-size optimization. CONCLUSION: Helmet pTx arrays significantly outperformed cylindrical arrays with the same number of channels in local SAR and flip-angle uniformity metrics. This improvement was especially pronounced for non-transverse slice excitations. Loop diameter optimization for helmets appears to favor large loops, compatible with nearest-neighbor decoupling by overlap.

An 8-channel Tx dipole and 20-channel Rx loop coil array for MRI of the cervical spinal cord at 7 Tesla.

The quality of cervical spinal cord images can be improved by the use of tailored radiofrequency (RF) coil solutions for ultrahigh field imaging; however, very few commercial and research 7-T RF coils currently exist for the spinal cord, and in particular, those with parallel transmission (pTx) capabilities. This work presents the design, testing, and validation of a pTx/Rx coil for the human neck and cervical/upper thoracic spinal cord. The pTx portion is composed of eight dipoles to ensure high homogeneity over this large region of the spinal cord. The Rx portion is made up of twenty semiadaptable overlapping loops to produce high signal-to-noise ratio (SNR) across the patient population. The coil housing is designed to facilitate patient positioning and comfort, while also being tight fitting to ensure high sensitivity. We demonstrate RF shimming capabilities to optimize B1 + uniformity, power efficiency, and/or specific absorption rate efficiency. B1 + homogeneity, SNR, and g-factor were evaluated in adult volunteers and demonstrated excellent performance from the occipital lobe down to the T4-T5 level. We compared the proposed coil with two state-of-the-art head and head/neck coils, confirming its superiority in the cervical and upper thoracic regions of the spinal cord. This coil solution therefore provides a convincing platform for producing the high image quality necessary for clinical and research scanning of the upper spinal cord.

Prediction of experimental cardiac magnetostimulation thresholds using pig-specific body models.

PURPOSE: Modern high-amplitude gradient systems can be limited by the International Electrotechnical Commission 60601-2-33 cardiac stimulation (CS) limit, which was set in a conservative manner based on electrode experiments and E-field simulations in uniform ellipsoidal body models. Here, we show that coupled electromagnetic-electrophysiological modeling in detailed body and heart models can predict CS thresholds, suggesting that such modeling might lead to more detailed threshold estimates in humans. Specifically, we compare measured and predicted CS thresholds in eight pigs. METHODS: We created individualized porcine body models using MRI (Dixon for the whole body, CINE for the heart) that replicate the anatomy and posture of the animals used in our previous experimental CS study. We model the electric fields induced along cardiac Purkinje and ventricular muscle fibers and predict the electrophysiological response of these fibers, yielding CS threshold predictions in absolute units for each animal. Additionally, we assess the total modeling uncertainty through a variability analysis of the 25 main model parameters. RESULTS: Predicted and experimental CS thresholds agree within 19% on average (normalized RMS error), which is smaller than the 27% modeling uncertainty. No significant difference was found between the modeling predictions and experiments (p < 0.05, paired t-test). CONCLUSION: Predicted thresholds matched the experimental data within the modeling uncertainty, supporting the model validity. We believe that our modeling approach can be applied to study CS thresholds in humans for various gradient coils, body shapes/postures, and waveforms, which is difficult to do experimentally.

Peripheral nerve stimulation informed design of a high-performance asymmetric head gradient coil.

PURPOSE: Peripheral nerve stimulation (PNS) limits the image encoding performance of both body gradient coils and the latest generation of head gradients. We analyze a variety of head gradient design aspects using a detailed PNS model to guide the design process of a new high-performance asymmetric head gradient to raise PNS thresholds and maximize the usable image-encoding performance. METHODS: A novel three-layer coil design underwent PNS optimization involving PNS predictions of a series of candidate designs. The PNS-informed design process sought to maximize the usable parameter space of a coil with <10% nonlinearity in a 22 cm region of linearity, a relatively large inner diameter (44 cm), maximum gradient amplitude of 200 mT/m, and a high slew rate of 900 T/m/s. PNS modeling allowed identification and iterative adjustment of coil features with beneficial impact on PNS such as the number of winding layers, shoulder accommodation strategy, and level of asymmetry. PNS predictions for the final design were compared to measured thresholds in a constructed prototype. RESULTS: The final head gradient achieved up to 2-fold higher PNS thresholds than the initial design without PNS optimization and compared to existing head gradients with similar design characteristics. The inclusion of a third intermediate winding layer provided the additional degrees of freedom necessary to improve PNS thresholds without significant sacrifices to the other design metrics. CONCLUSION: Augmenting the design phase of a new high-performance head gradient coil by PNS modeling dramatically improved the usable image-encoding performance by raising PNS thresholds.

8-channel Tx dipole and 20-channel Rx loop coil array for MRI of the cervical spinal cord at 7 Tesla.

The quality of cervical spinal cord images can be improved by the use of tailored radiofrequency coil solutions for ultra-high field imaging; however, very few commercial and research 7 Tesla radiofrequency coils currently exist for the spinal cord, and in particular those with parallel transmit capabilities. This work presents the design, testing and validation of a pTx/Rx coil for the human neck and cervical/upper-thoracic spinal cord. The pTx portion is composed of 8 dipoles to ensure high homogeneity over this large region of the spinal cord. The Rx portion is made of 20 semi-adaptable overlapping loops to produce high Signal-to-noise ratio (SNR) across the patient population. The coil housing is designed to facilitate patient positioning and comfort, while being tight fitting to ensure high sensitivity. We demonstrate RF shimming capabilities to optimize B 1 + uniformity, power efficiency and/or specific absorption rate (SAR) efficiency. B 1 + homogeneity, SNR and g-factor was evaluated in adult volunteers and demonstrated excellent performance from the occipital lobe down to the T4-T5 level. We compared the proposed coil with two state-of-the-art head and head/neck coils, confirming its superiority in the cervical and upper-thoracic regions of the spinal cord. This coil solution therefore provides a convincing platform for producing the high image quality necessary for clinical and research scanning of the upper spinal cord.

Parallel transmission 2D RARE imaging at 7T with transmit field inhomogeneity mitigation and local SAR control.

PURPOSE: We develop and test a parallel transmit (pTx) pulse design framework to mitigate transmit field inhomogeneity with control of local specific absorption rate (SAR) in 2D rapid acquisition with relaxation enhancement (RARE) imaging at 7T. METHODS: We design large flip angle RF pulses with explicit local SAR constraints by numerical simulation of the Bloch equations. Parallel computation and analytical expressions for the Jacobian and the Hessian matrices are employed to reduce pulse design time. The refocusing-excitation "spokes" pulse pairs are designed to satisfy the Carr-Purcell-Meiboom-Gill (CPMG) condition using a combined magnitude least squares-least squares approach. RESULTS: In a simulated dataset, the proposed approach reduced peak local SAR by up to 56% for the same level of refocusing uniformity error and reduced refocusing uniformity error by up to 59% (from 32% to 7%) for the same level of peak local SAR compared to the circularly polarized birdcage mode of the pTx array. Using explicit local SAR constraints also reduced peak local SAR by up to 46% compared to an RF peak power constrained design. The excitation and refocusing uniformity error were reduced from 20%-33% to 4%-6% in single slice phantom experiments. Phantom experiments demonstrated good agreement between the simulated excitation and refocusing uniformity profiles and experimental image shading. CONCLUSION: PTx-designed excitation and refocusing CPMG pulse pairs can mitigate transmit field inhomogeneity in the 2D RARE sequence. Moreover, local SAR can be decreased significantly using pTx, potentially leading to better slice coverage, enabling larger flip angles or faster imaging.

Measurement of magnetostimulation thresholds in the porcine heart.

PURPOSE: Powerful MRI gradient systems can surpass the International Electrotechnical Commission (IEC) 60601-2-33 limit for cardiac stimulation (CS), which was determined by simple electromagnetic simulations and electrode stimulation experiments. Only a few canine studies measured magnetically induced CS thresholds in vivo and extrapolating them to human safety limits can be challenging. METHODS: We measured cardiac magnetostimulation thresholds in 10 healthy, anesthetized pigs using capacitors discharged into a flat spiral coil to produce damped sinusoidal waveforms with effective stimulus duration ts,eff  = 0.45 ms. Electrocardiography (ECG), blood pressure, and peripheral oximetry signals were recorded to determine threshold coil currents yielding cardiac capture. Dixon and CINE MR volumes from each animal were segmented to generate porcine-specific electromagnetic models to calculate dB/dt and E-field values in the porcine heart at threshold. For comparison, we also simulated maximum dB/dt and E-field values created by three MRI gradient systems in the heart of a human body model. RESULTS: The average dB/dt threshold estimated in the porcine heart was 1.66 ± 0.23 kT/s, which is 11-fold greater than the IEC dB/dt limit at ts,eff  = 0.45 ms, and 31-fold greater than the maximum value created by the investigated MRI gradients in the human heart. The average E-field threshold estimated in the porcine heart was 92.9 ± 13.5 V/m, which is 6-fold greater than the IEC E-field limit at ts,eff  = 0.45 ms and 37-fold greater than the maximum gradient-induced E-field in the human heart. CONCLUSION: This first measurement of cardiac magnetostimulation thresholds in pigs indicates that the IEC cardiac safety limit is conservative for the investigated stimulus duration (ts,eff  = 0.45 ms).

A patient-friendly 16-channel transmit/64-channel receive coil array for combined head-neck MRI at 7 Tesla.

PURPOSE: To extend the coverage of brain coil arrays to the neck and cervical-spine region to enable combined head and neck imaging at 7 Tesla (T) ultra-high field MRI. METHODS: The coil array structures of a 64-channel receive coil and a 16-channel transmit coil were merged into one anatomically shaped close-fitting housing. Transmit characteristics were evaluated in a B1+ -field mapping study and an electromagnetic model. Receive SNR and the encoding capability for accelerated imaging were evaluated and compared with a commercially available 7 T brain array coil. The performance of the head-neck array coil was demonstrated in human volunteers using high-resolution accelerated imaging. RESULTS: In the brain, the SNR matches the commercially available 32-channel brain array and showed improvements in accelerated imaging capabilities. More importantly, the constructed coil array improved the SNR in the face area, neck area, and cervical spine by a factor of 1.5, 3.4, and 5.2, respectively, in regions not covered by 32-channel brain arrays at 7 T. The interelement coupling of the 16-channel transmit coil ranged from -14 to -44 dB (mean = -19 dB, adjacent elements <-18 dB). The parallel 16-channel transmit coil greatly facilitates B1+ field shaping required for large FOV neuroimaging at 7 T. CONCLUSION: This new head-neck array coil is the first demonstration of a device of this nature used for combined full-brain, head-neck, and cervical-spine imaging at 7 T. The array coil is well suited to provide large FOV images, which potentially improves ultrahigh field neuroimaging applications for clinical settings.

A 31-channel integrated "AC/DC" B0 shim and radiofrequency receive array coil for improved 7T MRI.

PURPOSE: To test an integrated "AC/DC" array approach at 7T, where B0 inhomogeneity poses an obstacle for functional imaging, diffusion-weighted MRI, MR spectroscopy, and other applications. METHODS: A close-fitting 7T 31-channel (31-ch) brain array was constructed and tested using combined Rx and ΔB0 shim channels driven by a set of rapidly switchable current amplifiers. The coil was compared to a shape-matched 31-ch reference receive-only array for RF safety, signal-to-noise ratio (SNR), and inter-element noise correlation. We characterize the coil array's ability to provide global and dynamic (slice-optimized) shimming using ΔB0 field maps and echo planar imaging (EPI) acquisitions. RESULTS: The SNR and average noise correlation were similar to the 31-ch reference array. Global and slice-optimized shimming provide 11% and 40% improvements respectively compared to baseline second-order spherical harmonic shimming. Birdcage transmit coil efficiency was similar for the reference and AC/DC array setups. CONCLUSION: Adding ΔB0 shim capability to a 31-ch 7T receive array can significantly boost 7T brain B0 homogeneity without sacrificing the array's rdiofrequency performance, potentially improving ultra-high field neuroimaging applications that are vulnerable to off-resonance effects.

A Huygens' surface approach to rapid characterization of peripheral nerve stimulation.

PURPOSE: Peripheral nerve stimulation (PNS) modeling has a potential role in designing and operating MRI gradient coils but requires computationally demanding simulations of electromagnetic fields and neural responses. We demonstrate compression of an electromagnetic and neurodynamic model into a single versatile PNS matrix (P-matrix) defined on an intermediary Huygens' surface to allow fast PNS characterization of arbitrary coil geometries and body positions. METHODS: The Huygens' surface approach divides PNS prediction into an extensive pre-computation phase of the electromagnetic and neurodynamic responses, which is independent of coil geometry and patient position, and a fast coil-specific linear projection step connecting this information to a specific coil geometry. We validate the Huygens' approach by performing PNS characterizations for 21 body and head gradients and comparing them with full electromagnetic-neurodynamic modeling. We demonstrate the value of Huygens' surface-based PNS modeling by characterizing PNS-optimized coil windings for a wide range of patient positions and poses in two body models. RESULTS: The PNS prediction using the Huygens' P-matrix takes less than a minute (instead of hours to days) without compromising numerical accuracy (error ≤ 0.1%) compared to the full simulation. Using this tool, we demonstrate that coils optimized for PNS at the brain landmark using a male model can also improve PNS for other imaging applications (cardiac, abdominal, pelvic, and knee imaging) in both male and female models. CONCLUSION: Representing PNS information on a Huygens' surface extended the approach's ability to assess PNS across body positions and models and test the robustness of PNS optimization in gradient design.

Mapping the subcortical connectivity of the human default mode network.

The default mode network (DMN) mediates self-awareness and introspection, core components of human consciousness. Therapies to restore consciousness in patients with severe brain injuries have historically targeted subcortical sites in the brainstem, thalamus, hypothalamus, basal forebrain, and basal ganglia, with the goal of reactivating cortical DMN nodes. However, the subcortical connectivity of the DMN has not been fully mapped, and optimal subcortical targets for therapeutic neuromodulation of consciousness have not been identified. In this work, we created a comprehensive map of DMN subcortical connectivity by combining high-resolution functional and structural datasets with advanced signal processing methods. We analyzed 7 Tesla resting-state functional MRI (rs-fMRI) data from 168 healthy volunteers acquired in the Human Connectome Project. The rs-fMRI blood-oxygen-level-dependent (BOLD) data were temporally synchronized across subjects using the BrainSync algorithm. Cortical and subcortical DMN nodes were jointly analyzed and identified at the group level by applying a novel Nadam-Accelerated SCAlable and Robust (NASCAR) tensor decomposition method to the synchronized dataset. The subcortical connectivity map was then overlaid on a 7 Tesla 100 µm ex vivo MRI dataset for neuroanatomic analysis using automated segmentation of nuclei within the brainstem, thalamus, hypothalamus, basal forebrain, and basal ganglia. We further compared the NASCAR subcortical connectivity map with its counterpart generated from canonical seed-based correlation analyses. The NASCAR method revealed that BOLD signal in the central lateral nucleus of the thalamus and ventral tegmental area of the midbrain is strongly correlated with that of the DMN. In an exploratory analysis, additional subcortical sites in the median and dorsal raphe, lateral hypothalamus, and caudate nuclei were correlated with the cortical DMN. We also found that the putamen and globus pallidus are negatively correlated (i.e., anti-correlated) with the DMN, providing rs-fMRI evidence for the mesocircuit hypothesis of human consciousness, whereby a striatopallidal feedback system modulates anterior forebrain function via disinhibition of the central thalamus. Seed-based analyses yielded similar subcortical DMN connectivity, but the NASCAR result showed stronger contrast and better spatial alignment with dopamine immunostaining data. The DMN subcortical connectivity map identified here advances understanding of the subcortical regions that contribute to human consciousness and can be used to inform the selection of therapeutic targets in clinical trials for patients with disorders of consciousness.

Safety and image quality at 7T MRI for deep brain stimulation systems: Ex vivo study with lead-only and full-systems.

Ultra-high field MRI at 7 T can produce much better visualization of sub-cortical structures compared to lower field, which can greatly help target verification as well as overall treatment monitoring for patients with deep brain stimulation (DBS) implants. However, use of 7 T MRI for such patients is currently contra-indicated by guidelines from the device manufacturers due to the safety issues. The aim of this study was to provide an assessment of safety and image quality of ultra-high field magnetic resonance imaging at 7 T in patients with deep brain stimulation implants. We performed experiments with both lead-only and complete DBS systems implanted in anthropomorphic phantoms. RF heating was measured for 43 unique patient-derived device configurations. Magnetic force measurements were performed according to ASTM F2052 test method, and device integrity was assessed before and after experiments. Finally, we assessed electrode artifact in a cadaveric brain implanted with an isolated DBS lead. RF heating remained below 2°C, similar to a fever, with the 95% confidence interval between 0.38°C-0.52°C. Magnetic forces were well below forces imposed by gravity, and thus not a source of concern. No device malfunctioning was observed due to interference from MRI fields. Electrode artifact was most noticeable on MPRAGE and T2*GRE sequences, while it was minimized on T2-TSE images. Our work provides the safety assessment of ultra-high field MRI at 7 T in patients with DBS implants. Our results suggest that 7 T MRI may be performed safely in patients with DBS implants for specific implant models and MRI hardware.

Safety and imaging performance of two-channel RF shimming for fetal MRI at 3T.

PURPOSE: This study investigates whether two-channel radiofrequency (RF) shimming can improve imaging without increasing specific absorption rate (SAR) for fetal MRI at 3T. METHODS: Transmit field ( B1+ ) average and variation in the fetus was simulated in seven numerical pregnant body models. Safety was quantified by maternal and fetal peak local SAR and fetal average SAR. The shim parameter space was divided into improved B1+ (magnitude and homogeneity) and improved SAR regions, and an overlap where RF shimming improved both classes of metrics compared with birdcage mode was assessed. Additionally, the effect of fetal position, tissue detail, and dielectric properties on transmit field and SAR was studied. RESULTS: A region of subject-specific RF shim parameter space improving both B1+ and SAR metrics was found for five of the seven models. Optimizing only B1+ metrics improved B1+ efficiency across models by 15% on average and 28% for the best-case model. B1+ variation improved by 26% on average and 49% for the best case. However, for these shim settings, fetal SAR increased by up to 106%. The overlap region, where both B1+ and SAR metrics improve, showed an average B1+ efficiency improvement of 6% on average across models and 19% for the best-case model. B1+ variation improved by 13% on average and 40% for the best case. RFS could also decrease maternal/fetal SAR by up to 49%/58%. CONCLUSION: RF shimming can improve imaging compared with birdcage mode without increasing fetal and maternal SAR when a patient-specific SAR model is incorporated into the shimming procedure.