Design and validation of a patient-specific stereotactic frame for transcranial ultrasound therapy.

Transcranial focused ultrasound (tFUS) procedures such as neuromodulation and blood brain barrier opening require precise focus placement within the brain. MRI is currently the most reliable tool for focus localization but can be prohibitive for procedures requiring recurrent therapies. We designed, fabricated, and characterized a patient-specific, 3D-printed, stereotactic frame for repeated tFUS therapy. The frame is compact with minimal footprint, can be removed and re-secured between treatments while maintaining sub-mm accuracy and will allow for precise and repeatable transcranial FUS treatment without the need for MR-guidance following the initial calibration scan. Focus localization and repeatability were assessed via MR-thermometry and MR-ARFI on an ex vivo skull-phantom and in vivo non-human primates (NHP), respectively. Focal localization, registration, steering, and re-steering were accomplished during the initial MRI calibration scan session. Keeping steering coordinates fixed in subsequent therapy and imaging sessions, we found good agreement between steered foci and intended target, with target registration error of 1.2 ± 0.3 (n = 4, ex vivo) and 1.0 ± 0.5 (n = 3, in vivo) mm. Focus position (steered and non-steered) was consistent, with sub-mm variation in each dimension between studies. Our 3D-printed, patient-specific stereotactic frame can reliably position and orient the ultrasound transducer for repeated targeting of brain regions using a single MR-based calibration. The compact frame allows for high-precision tFUS to be carried out outside the magnet, and could help reduce the cost of tFUS treatments where repeated application of an ultrasound focus is required with high precision.

Practical targeting errors during optically tracked transcranial focused ultrasound using MR-ARFI and array-based steering.

OBJECTIVE: Transcranial focused ultrasound (tFUS) is being explored for neuroscience research and clinical applications due to its ability to affect precise brain regions noninvasively. The ability to target specific brain regions and localize the beam during these procedures is important for these applications to avoid damage and minimize off-target effects. Here, we present a method to combine optical tracking with magnetic resonance (MR) acoustic radiation force imaging to achieve targeting and localizing of the tFUS beam. This combined method provides steering coordinates to target brain regions within a clinically practical time frame. METHODS: Using an optically tracked hydrophone and bias correction with MR imaging we transformed the FUS focus coordinates into the MR space for targeting and error correction. We validated this method in vivo in 18 macaque FUS studies. RESULTS: Across these in vivo studies a single localization scan allowed for the average targeting error to be reduced from 4.8 mm to 1.4 mm and for multiple brain regions to be targeted with one transducer position. CONCLUSIONS: By reducing targeting error and providing the means to target multiple brain regions within a single session with high accuracy this method will allow further study of the effects of tFUS neuromodulation with more advanced approaches such as simultaneous dual or multi-site brain stimulation.

Small volume blood-brain barrier opening in macaques with a 1 MHz ultrasound phased array.

The use of focused ultrasound to open the blood-brain barrier (BBB) has the potential to deliver drugs to specific regions of the brain. The size of the BBB opening and ability to localize the opening determines the spatial extent and is a limiting factor in many applications of BBB opening where targeting a small brain region is desired. Here we evaluate the performance of a system designed for small opening volumes and highlight the unique challenges associated with pushing the spatial precision of this technique. To achieve small volume openings in cortical regions of the macaque brain, we tested a custom 1 MHz array transducer integrated into a magnetic resonance image-guided focused ultrasound system. Using real-time cavitation monitoring, we demonstrated twelve instances of single sonication, small volume BBB opening with average volumes of 59 ± 37 mm3 and 184 ± 2 mm3 in cortical and subcortical targets, respectively. We found high correlation between subject-specific acoustic simulations and observed openings when incorporating grey matter segmentation (R2 = 0.8577), and the threshold for BBB opening based on simulations was 0.53 MPa. Analysis of MRI-based safety assessment and cavitation signals indicate a safe pressure range for 1 MHz BBB opening and suggest that our system can be used to deliver drugs and gene therapy to small brain regions.

Disrupting nociceptive information processing flow through transcranial focused ultrasound neuromodulation of thalamic nuclei.

BACKGROUND: MRI-guided transcranial focused ultrasound (MRgFUS) as a next-generation neuromodulation tool can precisely target and stimulate deep brain regions with high spatial selectivity. Combined with MR-ARFI (acoustic radiation force imaging) and using fMRI BOLD signal as functional readouts, our previous studies have shown that low-intensity FUS can excite or suppress neural activity in the somatosensory cortex. OBJECTIVE: To investigate whether low-intensity FUS can suppress nociceptive heat stimulation-induced responses in thalamic nuclei during hand stimulation, and to determine how this suppression influences the information processing flow within nociception networks. FINDINGS: BOLD fMRI activations evoked by 47.5 °C heat stimulation of hand were detected in 24 cortical regions, which belong to sensory, affective, and cognitive nociceptive networks. Concurrent delivery of low-intensity FUS pulses (650 kHz, 550 kPa) to the predefined heat nociceptive stimulus-responsive thalamic centromedial_parafascicular (CM_para), mediodorsal (MD), ventral_lateral (VL_ and ventral_lateral_posteroventral (VLpv) nuclei suppressed their heat responses. Off-target cortical areas exhibited reduced, enhanced, or no significant fMRI signal changes, depending on the specific areas. Differentiable thalamocortical information flow during the processing of nociceptive heat input was observed, as indicated by the time to reach 10% or 30% of the heat-evoked BOLD signal peak. Suppression of thalamic heat responses significantly altered nociceptive processing flow and direction between the thalamus and cortical areas. Modulation of contralateral versus ipsilateral areas by unilateral thalamic activity differed. Signals detected in high-order cortical areas, such as dorsal frontal (DFC) and ventrolateral prefrontal (vlPFC) cortices, exhibited faster response latencies than sensory areas. CONCLUSIONS: The concurrent delivery of FUS suppressed nociceptive heat response in thalamic nuclei and disrupted the nociceptive network. This study offers new insights into the causal functional connections within the thalamocortical networks and demonstrates the modulatory effects of low-intensity FUS on nociceptive information processing.

Longitudinal multiparametric MRI of traumatic spinal cord injury in animal models.

Multi-parametric MRI (mpMRI) technology enables non-invasive and quantitative assessments of the structural, molecular, and functional characteristics of various neurological diseases. Despite the recognized importance of studying spinal cord pathology, mpMRI applications in spinal cord research have been somewhat limited, partly due to technical challenges associated with spine imaging. However, advances in imaging techniques and improved image quality now allow longitudinal investigations of a comprehensive range of spinal cord pathological features by exploiting different endogenous MRI contrasts. This review summarizes the use of mpMRI techniques including blood oxygenation level-dependent (BOLD) functional MRI (fMRI), diffusion tensor imaging (DTI), quantitative magnetization transfer (qMT), and chemical exchange saturation transfer (CEST) MRI in monitoring different aspects of spinal cord pathology. These aspects include cyst formation and axonal disruption, demyelination and remyelination, changes in the excitability of spinal grey matter and the integrity of intrinsic functional circuits, and non-specific molecular changes associated with secondary injury and neuroinflammation. These approaches are illustrated with reference to a nonhuman primate (NHP) model of traumatic cervical spinal cord injuries (SCI). We highlight the benefits of using NHP SCI models to guide future studies of human spinal cord pathology, and demonstrate how mpMRI can capture distinctive features of spinal cord pathology that were previously inaccessible. Furthermore, the development of mechanism-based MRI biomarkers from mpMRI studies can provide clinically useful imaging indices for understanding the mechanisms by which injured spinal cords progress and repair. These biomarkers can assist in the diagnosis, prognosis, and evaluation of therapies for SCI patients, potentially leading to improved outcomes.

Small volume blood-brain barrier opening in macaques with a 1 MHz ultrasound phased array.

Focused ultrasound blood-brain barrier (BBB) opening is a promising tool for targeted delivery of therapeutic agents into the brain. The volume of opening determines the extent of therapeutic administration and sets a lower bound on the size of targets which can be selectively treated. We tested a custom 1 MHz array transducer optimized for cortical regions in the macaque brain with the goal of achieving small volume openings. We integrated this device into a magnetic resonance image guided focused ultrasound system and demonstrated twelve instances of small volume BBB opening with average opening volumes of 59 ± 37 mm 3 and 184 ± 2 mm 3 in cortical and subcortical targets, respectively. We developed real-time cavitation monitoring using a passive cavitation detector embedded in the array and characterized its performance on a bench-top flow phantom mimicking transcranial BBB opening procedures. We monitored cavitation during in-vivo procedures and compared cavitation metrics against opening volumes and safety outcomes measured with FLAIR and susceptibility weighted MR imaging. Our findings show small BBB opening at cortical targets in macaques and characterize the safe pressure range for 1 MHz BBB opening. Additionally, we used subject-specific simulations to investigate variance in measured opening volumes and found high correlation (R 2 = 0.8577) between simulation predictions and observed measurements. Simulations suggest the threshold for 1 MHz BBB opening was 0.53 MPa. This system enables BBB opening for drug delivery and gene therapy to be targeted to more specific brain regions.

Detection of laser-associated heating in the brain during simultaneous fMRI and optogenetic stimulation.

PURPOSE: To calculate temperatures from T2 *-weighted images collected during optogenetic fMRI based on proton resonance frequency (PRF) shift thermometry, to monitor confounding heating effects and determine appropriate light parameters for optogenetic stimulation. METHODS: fMRI is mainly based on long-TE gradient-recalled echo acquisitions that are also suitable for measuring small temperature changes via the PRF shift. A motion- and respiration-robust processing pipeline was developed to calculate temperature changes based on the PRF shift directly from the T2 *-weighted images collected for fMRI with a two-shot 2D gradient-recalled echo-EPI sequence at 9.4T. Optogenetic fMRI protocols which differed in stimulation durations (3, 6 and 9 s) within a total block duration of 30 s were applied in a squirrel monkey to validate the methods with blue and green light (20 Hz, 30 mW) delivery interleaved between periods. General linear modeling was performed on the resulting time series temperature maps to verify if light delivery with each protocol resulted in significant heating in the brain around the optical fiber. RESULTS: The temperature SD was 0.05°C with the proposed imaging protocol and processing. Statistical analysis showed that the optogenetic stimulation protocol with a 3 s stimulation duration did not result in significant temperature rises. Significant temperature rises up to 0.13°C (p < 0. 05) were observed with 6 and 9 s stimulation durations for blue and green light. CONCLUSION: The proposed processing pipeline can be useful for the design of optogenetic stimulation protocols and for monitoring confounding heating effects.

Machine-learning-based video analysis of grasping behavior during recovery from cervical spinal cord injury.

Comprehensive characterizations of hand grasping behaviors after cervical spinal cord injuries are fundamental for developing rehabilitation strategies to promote recovery in spinal-cord-injured primates. We used the machine-learning-based video analysis software, DeepLabCut, to sensitively quantify kinematic aspects of grasping behavioral deficits in squirrel monkeys with C5-level spinal cord injuries. Three squirrel monkeys were trained to grasp sugar pellets from wells of varying depths before and after a left unilateral lesion of the cervical dorsal column. Using DeepLabCut, we identified post-lesion deficits in kinematic grasping behavior that included changes in digit orientation, increased variance in vertical and horizontal digit movement, and longer time to complete the task. While video-based analyses of grasping behavior demonstrated deficits in fine-scale digit function that persisted through at least 14 weeks post-injury, traditional end-point behavioral analyses showed a recovery of global hand function as evidenced by recovery of the proportion of successful retrievals by approximately 14 weeks post-injury. The combination of traditional end-point and video-based kinematic analyses provides a more comprehensive characterization of grasping behavior and highlights that global grasping performance may recover despite persistent fine-scale kinematic deficits in digit function. Machine-learning-based video analysis of kinematic digit function, in conjunction with traditional end-point behavioral analyses of grasping behavior, provide sensitive and specific indices for monitoring recovery of fine-grained hand sensorimotor behavior after spinal cord injury that can aid future studies that seek to develop targeted therapeutic interventions for improving behavioral outcomes.

Differential dose responses of transcranial focused ultrasound at brain regions indicate causal interactions.

We have previously shown that focused ultrasound (FUS) pulses in low pressure range exerted bidirectional and brain state-dependent neuromodulation in the nonhuman primate somatosensory cortices by fMRI. Here we aim to gain insights about the proposed neuron selective modulation of FUS and probe feedforward versus feedback interactions by simultaneously quantifying the stimulus (FUS pressures: 925, 425, 250 kPa) and response (% BOLD fMRI changes) function at the targeted area 3a/3b and off-target cortical areas at 7T. In resting-state, lowered intensities of FUS resulted in decreased fMRI signal changes at the target area 3a/3b and off-target area 1/2, S2, MCC, insula and auditory cortex, and no signal difference in thalamic VPL and MD nuclei. In activated states, concurrent high-intensity FUS significantly enhanced touch-evoked signals in area 1/2. Medium- and low-intensity FUS significantly suppressed touch-evoked BOLD signals in all areas except in the auditory cortex, VPL and MD thalamic nuclei. Distinct state dependent and dose-response curves led us to hypothesize that FUS's neuromodulatory effects may be mediated through preferential activation of different populations of neurons. Area 3a/3b may have distinct causal feedforward and feedback interactions with Area 1/2, S2, MCC, insula, and VPL. FUS offers a noninvasive neural stimulation tool for dissecting brain circuits and probing causal functional connections.

Differential Recovery of Submodality Touch Neurons and Interareal Communication in Sensory Input-Deprived Area 3b and S2 Cortices.

Cortical reactivation and regain of interareal functional connections have been linked to the recovery of hand grasping behavior after loss of sensory inputs in primates. We investigated contributions of neurons in two hierarchically organized somatosensory areas, 3b and S2, by characterizing local field potential (LFP) and multiunit spiking activity in five states (rest, stimulus-on, sustained, stimulus-off, and induced) and interareal communication after grasping behavior of dorsal column lesioned male squirrel monkeys had mostly recovered. Compared with normal cortex, fMRI, LFP, and spiking response magnitudes to step indentations were significantly weaker. The sustained component of the spiking recovered much better than the stimulus-off response. Correlation between overall spiking and γ LFP remained strong within each recovered areas 3b and S2. The interareal correlations of γ LFP were severely disrupted, except in the resting and stimulus-on periods. Interareal correlation of spiking was disrupted in the stimulus-off period only. In summary, submodality of low threshold mechanoreceptive neurons recovered differentially in input-deprived area 3b and S2 when impaired global hand grasping behavior returned. Slow-adapting-like neurons recovered, whereas rapid-adapting-like neurons did not. Interareal communications were also severely compromised. We propose that slow-adapting-like neurons and afferents in recovered area 3b and S2 mediate recovery of impaired grasping behavior after dorsal column tract lesion.SIGNIFICANCE STATEMENT Sensory feedback is essential for execution of hand grasping behavior in primates. Reactivations of somatosensory cortices have been attributed to recovery of such behavior after loss of sensory inputs via largely unknown mechanisms. In input-deprived area 3b and S2 cortex, after hand grasping behavior mostly recovered, we found slow-adapting-like neurons were greatly recovered, whereas rapid-adapting-like neurons did not. Communications between area 3b and S2 neurons were severely compromised. We suggest that recovery of slow-adapting-like neurons in input-deprived area 3b and S2 may mediate the recovery of hand grasping behavior.

RF shielding designs for birdcage coils for preclinical MRI at 9.4 T.

Birdcage coils are widely used in preclinical MRI as they perform well, allow for quadrature drive, and can provide a homogeneous transmit field. Unlike in larger bore scanners, an RF shield is essential to avoid strong cross-talk with gradient coils that are in close proximity. However, gradient switching induces eddy currents that heat the shield and coil and impair the temporal signal-to-noise ratio (tSNR). The motivation of this study is to investigate the performance of different designs of RF shields on a birdcage coil used for high resolution functional MRI of small primates at 9.4 T. We found the choice of materials for RF shields significantly affected ghosting and tSNR in fMRI scans. Both ultrathin foils and a slotted pattern reduce eddy currents and improve imaging quality. Our results also demonstrate that a 9-um-thick copper foil is sufficiently thin to reduce the eddy current effects for high-resolution fMRI scans and there is no need for high-cost 4-um-thick foil. For slotted shields, our results demonstrate that the number of slots should be carefully considered, and an excessive number of slots can lead to a lower SNR and tSNR. We believe the results from this study can be used as a reference to design future RF coil shields selection for preclinical scanners.

A reduced aperture allows for transcranial focus localization at lower pressure.

Localizing the focus during transcranial focused ultrasound procedures is important to ensure accurate targeting of specific brain regions and interpretation of results. Magnetic resonance acoustic radiation force imaging uses the displacement induced by the ultrasound focus in the brain to localize the beam, but the high pressure required to displace brain tissue may cause damage or confounds during subsequent neuromodulatory experiments. Here, reduced apertures were applied to a phased array transducer to generate comparable displacement to the full aperture but with 20% lower free field pressure.

Graph theory analysis identified two hubs that connect sensorimotor and cognitive and cortical and subcortical nociceptive networks in the non-human primate.

Pain perception involves multiple brain regions and networks. Understanding how these brain networks work together is fundamental for appreciating network-wise changes reported in patients with chronic pain disorders. Parcellating pain related networks and understanding their causal relationships is the first step to understand how painful information is processed, integrated, and modulated, and it requires direct manipulation of specific brain regions. Nonhuman primates (NHP) offer an ideal model system to achieve these goals because cortical and subcortical regions in the NHP brain are established based on a variety of different types of data collected in a way that is not feasible or, at least, extremely difficult in humans (i.e., histology data, tract-tracing, intracerebral recordings). In addition, different methodological techniques can also help characterize and further understand these brain cortical and subcortical regions over the course of development. Here we used a heat nociceptive stimulation that is proven to elicit activity of nociceptive neurons in the cortex to refine and parcellate the whole brain nociceptive functional networks, to identify key network hubs, and to characterize network-wise temporal dynamic signatures using high-resolution fMRI. We first functionally localized 24 cortical and subcortical regions that responded to heat nociceptive stimuli (somatosensory area 1/2, area 3a/3b, S2, posterior insula (pIns), anterior insula, area 7b, posterior parietal cortex, anterior cingulate cortex (ACC), prefrontal cortex, caudate, and mediodorsal (MD) and ventral posterior lateral (VPL) thalamic nuclei) and used them as seeds in resting state fMRI (rsfMRI) data analysis. We applied both hierarchical clustering and graph-theory analyses of the pairwise rsfMRI correlation metrics and identified five cortical and one subcortical sub-networks: strong resting state functional connectivity (rsFC) between ACC and prefrontal regions, parietal cortex and area 7b, S2 and posterior insula, areas 3a/3b and 1/2 within the S1 cortex, and thalamic MD and caudate nuclei. The rsFC strengths between cortical areas within each subnetwork were significantly stronger than those between subcortical regions. Regions within each sub-network also exhibited highly correlated temporal dynamics at rest, but the overall dynamic patterns varied drastically across sub-networks. Graph-theory analysis identified the MD nucleus as a hub that connects subcortical and cortical nociceptive sub-networks. The S2-pIns connection joins the sensory and affective/cognitive sub-networks.

Functional networks in non-human primate spinal cord and the effects of injury.

Spontaneous fluctuations of Blood Oxygenation-Level Dependent (BOLD) MRI signal in a resting state have previously been detected and analyzed to describe intrinsic functional networks in the spinal cord of rodents, non-human primates and human subjects. In this study we combined high resolution imaging at high field with data-driven Independent Component Analysis (ICA) to i) delineate fine-scale functional networks within and between segments of the cervical spinal cord of monkeys, and also to ii) characterize the longitudinal effects of a unilateral dorsal column injury on these networks. Seven distinct functional hubs were revealed within each spinal segment, with new hubs detected at bilateral intermediate and gray commissure regions in addition to the bilateral dorsal and ventral horns previously reported. Pair-wise correlations revealed significantly stronger connections between hubs on the dominant hand side. Unilateral dorsal-column injuries disrupted predominantly inter-segmental rather than intra-segmental functional connectivities as revealed by correlation strengths and graph-theory based community structures. The effects of injury on inter-segmental connectivity were evident along the length of the cord both below and above the lesion region. Connectivity strengths recovered over time and there was revival of inter-segmental communities as animals recovered function. BOLD signals of frequency 0.01-0.033 Hz were found to be most affected by injury. The results in this study provide new insights into the intrinsic functional architecture of spinal cord and underscore the potential of functional connectivity measures to characterize changes in networks after an injury and during recovery.

A thermal nociceptive patch in the S2 cortex of nonhuman primates: a combined functional magnetic resonance imaging and electrophysiology study.

ABSTRACT: Human functional magnetic resonance imaging (fMRI) and behavioral studies have established the roles of cortical areas along the Sylvian fissure in sensing subjective pain. Yet, little is known about how sensory aspects of painful information are represented and processed by neurons in these regions and how their electrophysiological activities are related to fMRI signals. The current study aims to partially address this critical knowledge gap by performing fMRI-guided microelectrode mapping and recording studies in the homologous region of the parietal operculum in squirrel monkeys under light anesthesia. In each animal studied (n = 8), we detected mesoscale mini-networks for heat nociception in cortical regions around the lateral sulcus. Within the network, we discovered a ∼1.5 × 1.5-mm2-sized cortical patch that solely contained heat nociceptive neurons that aligned with the heat fMRI activation locus. These neurons responded slowly to thermal (heat and cold) nociceptive stimuli exclusively, continued firing for several seconds after the succession of stimulation, and exhibited multidigit receptive fields and high spontaneous firing rates. Similar to the fMRI responses, increasing temperatures in the nociceptive range led to a nonlinear increase in firing rates. The finding of a clustering of heat nociceptive neurons provides novel insights into the unique functional organization of thermal nociception in the S2 subregion of the primate brain. With fMRI, it supports the existence of a modality-preferred heat nociceptive patch that is spatially separated and intermingled with touch patches containing neurons with comparable receptive fields and the presence of functionally distinct mini-networks in primate opercular cortex.

Sensitivity and specificity of CEST and NOE MRI in injured spinal cord in monkeys.

PURPOSE: The sensitivity and accuracy of chemical exchange saturation transfer (CEST) and nuclear Overhauser enhancement (NOE) effects for assessing injury-associated changes in cervical spinal cords were evaluated in squirrel monkeys. Multiple interacting pools of protons, including one identified by an NOE at -1.6 ppm relative to water (NOE(-1.6)), were derived and quantified from fitting proton Z-spectra. The effects of down-sampled data acquisitions and corrections for non-specific factors including T1, semi-solid magnetization transfer, and direct saturation of free water (DS), were investigated. The overall goal is to develop a protocol for rapid data acquisition for assessing the molecular signatures of the injured spinal cord and its surrounding regions. METHODS: MRI scans were recorded of anesthetized squirrel monkeys at 9.4 T, before and after a unilateral dorsal column sectioning of the cervical spinal cord. Z-spectral images at 51 different RF offsets were acquired. The amplitudes of CEST and NOE effects from multiple proton pools were quantified using a six-pool Lorenzian fitting of each Z-spectrum (MTRmfit). In addition, down-sampled data using reduced selections of RF offsets were analyzed and compared. An apparent exchange-dependent relaxation (AREXmfit) method was also used to correct for non-specific factors in quantifying regional spectra around lesion sites. RESULTS: The parametric maps from multi-pool fitting using the complete sampling data (P51e) detected unilateral changes at and around the injury. The maps derived from selected twofold down-sampled data with appropriate interpolation (P26sI51) revealed quite similar spatial distributions of different pools as those obtained using P51e at each resonance shift. Across 10 subjects, both data acquisition schemes detected significant decreases in NOE(-3.5) and NOE(-1.6) and increases in DS(0.0) and CEST(3.5) at the lesion site relative to measures of the normal tissues before injury. AREXmfit of cysts and other abnormal tissues at and around the lesion site also exhibited significant changes, especially at 3.5, -1.6 and -3.5 ppm RF offsets. CONCLUSION: These results confirm that a reduced set of RF offsets and down sampling are adequate for CEST imaging of injured spinal cord and allow shorter imaging times and/or permit additional signal averaging. AREXmfit correction improved the accuracy of CEST and NOE measures. The results provide a rapid (~13 mins), sensitive, and accurate protocol for deriving multiple NOE and CEST effects simultaneously in spinal cord imaging at high field.

Bidirectional and state-dependent modulation of brain activity by transcranial focused ultrasound in non-human primates.

Transcranial focused ultrasound (FUS) stimulation under MRI guidance, coupled with functional MRI (fMRI) monitoring of effects, offers a precise, noninvasive technology to dissect functional brain circuits and to modulate altered brain functional networks in neurological and psychiatric disorders. Here we show that ultrasound at moderate intensities modulated neural activity bi-directionally. Concurrent sonication of somatosensory areas 3a/3b with 250 kHz FUS suppressed the fMRI signals produced there by peripheral tactile stimulation, while at the same time eliciting fMRI activation at inter-connected, off-target brain regions. Direct FUS stimulation of the cortex resulted in different degrees of BOLD signal changes across all five off-target regions, indicating that its modulatory effects on active and resting neurons differed. This is the first demonstration of the dual suppressive and excitative modulations of FUS on a specific functional circuit and of ability of concurrent FUS and MRI to evaluate causal interactions between functional circuits with neuron-class selectivity.

Considerations for ultrasound exposure during transcranial MR acoustic radiation force imaging.

The aim of this study was to improve the sensitivity of magnetic resonance-acoustic radiation force imaging (MR-ARFI) to minimize pressures required to localize focused ultrasound (FUS) beams, and to establish safe FUS localization parameters for ongoing ultrasound neuromodulation experiments in living non-human primates. We developed an optical tracking method to ensure that the MR-ARFI motion-encoding gradients (MEGs) were aligned with a single-element FUS transducer and that the imaged slice was prescribed at the optically tracked location of the acoustic focus. This method was validated in phantoms, which showed that MR-ARFI-derived displacement sensitivity is maximized when the MR-ARFI MEGs were maximally aligned with the FUS propagation direction. The method was then applied in vivo to acquire displacement images in two healthy macaque monkeys (M fascicularis) which showed the FUS beam within the brain. Temperature images were acquired using MR thermometry to provide an estimate of in vivo brain temperature changes during MR-ARFI, and pressure and thermal simulations of the acoustic pulses were performed using the k-Wave package which showed no significant heating at the focus of the FUS beam. The methods presented here will benefit the multitude of transcranial FUS applications as well as future human applications.

On the Relationship between MRI and Local Field Potential Measurements of Spatial and Temporal Variations in Functional Connectivity.

Correlations between fluctuations in resting state BOLD fMRI signals are interpreted as measures of functional connectivity (FC), but the neural basis of their origins and their relationships to specific features of underlying electrophysiologic activity, have not been fully established. In particular, the dependence of FC metrics on different frequency bands of local field potentials (LFPs), and the relationship of dynamic changes in BOLD FC to underlying temporal variations of LFP correlations, are not known. We compared the spatial profiles of resting state coherences of different frequency bands of LFP signals, with high resolution resting state BOLD FC measurements. We also compared the probability distributions of temporal variations of connectivity in both modalities using a Markov chain model-based approach. We analyzed data obtained from the primary somatosensory (S1) cortex of monkeys. We found that in areas 3b and 1 of S1 cortex, low frequency LFP signal fluctuations were the main contributions to resting state LFP coherence. Additionally, the dynamic changes of BOLD FC behaved most similarly to the LFP low frequency signal coherence. These results indicate that, within the S1 cortex meso-scale circuit studied, resting state FC measures from BOLD fMRI mainly reflect contributions from low frequency LFP signals and their dynamic changes.

Intrinsic functional architecture of the non-human primate spinal cord derived from fMRI and electrophysiology.

Resting-state functional MRI (rsfMRI) has recently revealed correlated signals in the spinal cord horns of monkeys and humans. However, the interpretation of these rsfMRI correlations as indicators of functional connectivity in the spinal cord remains unclear. Here, we recorded stimulus-evoked and spontaneous spiking activity and local field potentials (LFPs) from monkey spinal cord in order to validate fMRI measures. We found that both BOLD and electrophysiological signals elicited by tactile stimulation co-localized to the ipsilateral dorsal horn. Temporal profiles of stimulus-evoked BOLD signals covaried with LFP and multiunit spiking in a similar way to those observed in the brain. Functional connectivity of dorsal horns exhibited a U-shaped profile along the dorsal-intermediate-ventral axis. Overall, these results suggest that there is an intrinsic functional architecture within the gray matter of a single spinal segment, and that rsfMRI signals at high field directly reflect this underlying spontaneous neuronal activity.