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.

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.

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.

High spatial correspondence at a columnar level between activation and resting state fMRI signals and local field potentials.

Although blood oxygenation level-dependent (BOLD) fMRI has been widely used to map brain responses to external stimuli and to delineate functional circuits at rest, the extent to which BOLD signals correlate spatially with underlying neuronal activity, the spatial relationships between stimulus-evoked BOLD activations and local correlations of BOLD signals in a resting state, and whether these spatial relationships vary across functionally distinct cortical areas are not known. To address these critical questions, we directly compared the spatial extents of stimulated activations and the local profiles of intervoxel resting state correlations for both high-resolution BOLD at 9.4 T and local field potentials (LFPs), using 98-channel microelectrode arrays, in functionally distinct primary somatosensory areas 3b and 1 in nonhuman primates. Anatomic images of LFP and BOLD were coregistered within 0.10 mm accuracy. We found that the point spread functions (PSFs) of BOLD and LFP responses were comparable in the stimulus condition, and both estimates of activations were slightly more spatially constrained than local correlations at rest. The magnitudes of stimulus responses in area 3b were stronger than those in area 1 and extended in a medial to lateral direction. In addition, the reproducibility and stability of stimulus-evoked activation locations within and across both modalities were robust. Our work suggests that the intrinsic resolution of BOLD is not a limiting feature in practice and approaches the intrinsic precision achievable by multielectrode electrophysiology.

Discrete Modules and Mesoscale Functional Circuits for Thermal Nociception within Primate S1 Cortex.

This study addresses one long-standing question of whether functional separations are preserved for somatosensory modalities of touch, heat, and cold nociception within primate primary somatosensory (S1) cortex. This information is critical for understanding how the nature of pain is represented in the primate brain. Using a combination of submillimeter-resolution fMRI and microelectrode local field potential (LFP) and spike recordings, we identified spatially segregated cortical zones for processing touch and nociceptive heat and cold stimuli in somatotopically appropriate areas 3a, 3b, 1, and 2 of S1 in male monkeys. The distances between zones were comparable (∼3.4 mm) across stimulus modalities (heat, cold, and tactile), indicating the existence of uniform, modality-specific modules. Stimulus-evoked LFP maps validated the fMRI maps in areas 3b and 1. Isolation of heat and cold nociceptive neurons from the fMRI zones confirmed the validity of using fMRI to probe nociceptive regions and circuits. Resting-state fMRI analysis revealed distinct intrinsic functional circuits among functionally related zones. We discovered distinct modular structures and networks for thermal nociception within S1 cortex, a finding that has significant implications for studying chronic pain syndromes and guiding the selection of neuromodulation targets for chronic pain management.SIGNIFICANCE STATEMENT Primate S1 subregions contain discrete heat and cold nociceptive modules. Modules with the same properties exhibit strong functional connection. Nociceptive fMRI response coincides with LFP and spike activities of nociceptive neurons. Functional separation of heat and cold pain is retained within primate S1 cortex.

Correlated Disruption of Resting-State fMRI, LFP, and Spike Connectivity between Area 3b and S2 following Spinal Cord Injury in Monkeys.

This study aims to understand how functional connectivity (FC) between areas 3b and S2 alters following input deprivation and the neuronal basis of disrupted FC of resting-state fMRI signals. We combined submillimeter fMRI with microelectrode recordings to localize the deafferented digit regions in areas 3b and S2 by mapping tactile stimulus-evoked fMRI activations before and after cervical dorsal column lesion in each male monkey. An average afferent disruption of 97% significantly reduced fMRI, local field potential (LFP), and spike responses to stimuli in both areas. Analysis of resting-state fMRI signal correlation, LFP coherence, and spike cross-correlation revealed significantly reduced functional connectivity between deafferented areas 3b and S2. The degrees of reductions in stimulus responsiveness and FC after deafferentation differed across fMRI, LFP, and spiking signals. The reduction of FC was much weaker than that of stimulus-evoked responses. Whereas the largest stimulus-evoked signal drop (∼80%) was observed in LFP signals, the greatest FC reduction was detected in the spiking activity (∼30%). fMRI signals showed mild reductions in stimulus responsiveness (∼25%) and FC (∼20%). The overall deafferentation-induced changes were quite similar in areas 3b and S2 across signals. Here we demonstrated that FC strength between areas 3b and S2 was much weakened by dorsal column lesion, and stimulus response reduction and FC disruption in fMRI covary with those of LFP and spiking signals in deafferented areas 3b and S2. These findings have important implications for fMRI studies aiming to probe FC alterations in pathological conditions involving deafferentation in humans.SIGNIFICANCE STATEMENT By directly comparing fMRI, local field potential, and spike signals in both tactile stimulation and resting states before and after severe disruption of dorsal column afferent, we demonstrated that reduction in fMRI responses to stimuli is accompanied by weakened resting-state fMRI functional connectivity (FC) in input-deprived and reorganized digit regions in area 3b of the S1 and S2. Concurrent reductions in local field potential and spike FC validated the use of resting-state fMRI signals for probing neural intrinsic FC alterations in pathological deafferented cortex, and indicated that disrupted FC between mesoscale functionally highly related regions may contribute to the behavioral impairments.

Injury alters intrinsic functional connectivity within the primate spinal cord.

Recent demonstrations of correlated low-frequency MRI signal variations between subregions of the spinal cord at rest in humans, similar to those found in the brain, suggest that such resting-state functional connectivity constitutes a common feature of the intrinsic organization of the entire central nervous system. We report our detection of functional connectivity within the spinal cords of anesthetized squirrel monkeys at rest and show that the strength of connectivity within these networks is altered by the effects of injuries. By quantifying the low-frequency MRI signal correlations between different horns within spinal cord gray matter, we found distinct functional connectivity relationships between the different sensory and motor horns, a pattern that was similar to activation patterns evoked by nociceptive heat or tactile stimulation of digits. All horns within a single spinal segment were functionally connected, with the strongest connectivity occurring between ipsilateral dorsal and ventral horns. Each horn was strongly connected to the same horn on neighboring segments, but this connectivity reduced drastically along the spinal cord. Unilateral injury to the spinal cord significantly weakened the strength of the intrasegment horn-to-horn connectivity only on the injury side and in slices below the lesion. These findings suggest resting-state functional connectivity may be a useful biomarker of functional integrity in injured and recovering spinal cords.

High-resolution functional MRI identified distinct global intrinsic functional networks of nociceptive posterior insula and S2 regions in squirrel monkey brain.

Numerous functional imaging and electrophysiological studies in humans and animals indicate that the two contiguous areas of secondary somatosensory cortex (S2) and posterior insula (pIns) are core regions in nociceptive processing and pain perception. In this study, we tested the hypothesis that the S2-pIns connection serves as a hub for connecting distinct sensory and affective nociceptive processing networks in the squirrel monkey brain. At 9.4T, we first mapped the brain regions that respond to nociceptive heat stimuli with high-resolution fMRI, and then used seed-based resting-state fMRI (rsfMRI) analysis to delineate and refine the global intrinsic functional connectivity circuits of the proximal S2 and pIns regions. In each subject, nociceptive (47.5°C) heat-evoked fMRI activations were detected in many brain regions, including primary somatosensory (S1), S2, pIns, area 7b, anterior cingulate cortex (ACC), primary motor cortex, prefrontal cortex, supplementary motor area, thalamus, and caudate. Using the heat-evoked fMRI activation foci in S2 and pIns as the seeds, voxel-wise whole-brain resting-state functional connectivity (rsFC) analysis revealed strong functional connections between contralateral S2 and pIns, as well as their corresponding regions in the ipsilateral hemisphere. Spatial similarity and overlap analysis identified each region as part of two distinct intrinsic functional networks with 7% overlap: sensory S2-S1-area 7b and affective pIns-ACC-PCC networks. Moreover, a high degree of overlap was observed between the combined rsFC maps of nociceptive S2 and pIns regions and the nociceptive heat-evoked activation map. In summary, our study provides evidence for the existence of two distinct intrinsic functional networks for S2 and pIns nociceptive regions, and these two networks are joined via the S2-pIns connection. Brain regions that are involved in processing nociceptive inputs are also highly interconnected at rest. The presence of robust and distinct S1-S2-area 7b and pIns-ACC-PCC rsFC networks under anesthesia underscores their fundamental roles in processing nociceptive information.

Differential fMRI Activation Patterns to Noxious Heat and Tactile Stimuli in the Primate Spinal Cord.

Mesoscale local functional organizations of the primate spinal cord are largely unknown. Using high-resolution fMRI at 9.4 T, we identified distinct interhorn and intersegment fMRI activation patterns to tactile versus nociceptive heat stimulation of digits in lightly anesthetized monkeys. Within a spinal segment, 8 Hz vibrotactile stimuli elicited predominantly fMRI activations in the middle part of ipsilateral dorsal horn (iDH), along with significantly weaker activations in ipsilateral (iVH) and contralateral (cVH) ventral horns. In contrast, nociceptive heat stimuli evoked widespread strong activations in the superficial part of iDH, as well as in iVH and contralateral dorsal (cDH) horns. As controls, only weak signal fluctuations were detected in the white matter. The iDH responded most strongly to both tactile and heat stimuli, whereas the cVH and cDH responded selectively to tactile versus nociceptive heat, respectively. Across spinal segments, iDH activations were detected in three consecutive segments in both tactile and heat conditions. Heat responses, however, were more extensive along the cord, with strong activations in iVH and cDH in two consecutive segments. Subsequent subunit B of cholera toxin tracer histology confirmed that the spinal segments showing fMRI activations indeed received afferent inputs from the stimulated digits. Comparisons of the fMRI signal time courses in early somatosensory area 3b and iDH revealed very similar hemodynamic stimulus-response functions. In summary, we identified with fMRI distinct segmental networks for the processing of tactile and nociceptive heat stimuli in the cervical spinal cord of nonhuman primates. Significance statement: This is the first fMRI demonstration of distinct intrasegmental and intersegmental nociceptive heat and touch processing circuits in the spinal cord of nonhuman primates. This study provides novel insights into the local functional organizations of the primate spinal cord for pain and touch, information that will be valuable for designing and optimizing therapeutic interventions for chronic pain management.

Correlated inter-regional variations in low frequency local field potentials and resting state BOLD signals within S1 cortex of monkeys.

The hypothesis that specific frequency components of the spontaneous local field potentials (LFPs) underlie low frequency fluctuations of resting state fMRI (rsfMRI) signals was tested. The previous analyses of rsfMRI signals revealed differential inter-regional correlations among areas 3a, 3b, and 1 of primary somatosensory cortex (S1) in anesthetized monkeys (Wang et al. [2013]: Neuron 78:1116-1126). Here LFP band(s) which correlated between S1 regions, and how these inter-regional correlation differences covaried with rsfMRI signals were examined. LFP signals were filtered into seven bands (delta, theta, alpha, beta, gamma low, gamma high, and gamma very high), and then a Hilbert transformation was applied to obtain measures of instantaneous amplitudes and temporal lags between regions of interest (ROI) digit-digit pairs (areas 3b-area 1, area 3a-area 1, area 3a-area 3b) and digit-face pairs (area 3b-face, area 1-face, and area 3a-face). It was found that variations in the inter-regional correlation strengths between digit-digit and digit-face pairs in the delta (1-4 Hz), alpha (9-14 Hz), beta (15-30 Hz), and gamma (31-50 Hz) bands parallel those of rsfMRI signals to varying degrees. Temporal lags between digit-digit area pairs varied across LFP bands, with area 3a mostly leading areas 1/2 and 3b. In summary, the data demonstrates that the low and middle frequency range (1-50 Hz) of spontaneous LFP signals similarly covary with the low frequency fluctuations of rsfMRI signals within local circuits of S1, supporting a neuronal electrophysiological basis of rsfMRI signals. Inter-areal LFP temporal lag differences provided novel insights into the directionality of information flow among S1 areas at rest. Hum Brain Mapp 37:2755-2766, 2016. © 2016 Wiley Periodicals, Inc.

Parallel functional reorganizations of somatosensory areas 3b and 1, and S2 following spinal cord injury in squirrel monkeys.

Multiple somatosensory cortices of adult primates reorganize following spinal cord injury, but little is known about the temporal dynamics and inter-areal differences of the reorganization. Using longitudinal high-resolution fMRI in combination with microelectrode recordings and tracer histology, we previously illustrated a two-phase dynamic spatial reorganization of digit representations in area 3b within weeks after a unilateral lesion of the dorsal column in squirrel monkeys (Chen et al., 2012). Here we report that higher-order area 1 and secondary somatosensory cortex (S2) underwent similar spatial reorganizations, which were characterized by shifted and expanded digit activations at week 4 after lesion, which then shifted back and contracted by week 8. In addition, the responsiveness of areas 3b and 1, and S2, as measured by the magnitude of the BOLD signal change to tactile stimuli, was reduced markedly at 4 weeks and then recovered to ~50% of the prelesion level at 8 weeks, a time when behavioral recovery was complete, as assessed by successful food retrieval rates. Across animals, the extents of spatial reorganizations and changes in cortical responsiveness and activation sizes in all three areas were correlated with the degree of afferent disruption. In summary, our data show that more severe afferent disruption was associated with greater cortical plasticity and behavioral impairment. Reorganization that occurred in area 3b, area 1, and S2 were similar across most measures.

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.

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.

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.

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.

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.

Active coping with stress suppresses glucose metabolism in the rat hypothalamus.

We used 18F-fluorodeoxyglucose small-animal positron-emission tomography to determine whether different styles of coping with stress are associated with different patterns of neuronal activity in the hypothalamus. Adult rats were subjected to immobilization (IMO)-stress or to a non-immobilized condition for 30 min, in random order on separate days, each of which was followed by brain-scanning. Some rats in the immobilized condition were allowed to actively cope with the stress by chewing a wooden stick during IMO, while the other immobilized rats were given nothing to chew on. Voxel-based statistical analysis of the brain imaging data shows that chewing counteracted the stress-induced increased glucose uptake in the hypothalamus to the level of the non-immobilized condition. Region-of-interest analysis of the glucose uptake values further showed that chewing significantly suppressed stress-induced increased glucose uptake in the paraventricular hypothalamic nucleus and the anterior hypothalamic area but not in the lateral hypothalamus. Together with the finding that the mean plasma corticosterone concentration at the termination of the IMO was also significantly suppressed when rats had an opportunity to chew a wooden stick, our results showed that active coping by chewing inhibited the activation of the hypothalamic-pituitary-adrenal axis to reduce the endocrine stress response.

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.