Design and application of a multimodality-compatible 1Tx/6Rx RF coil for monkey brain MRI at 7T.

OBJECTIVE: Blood-oxygen-level-dependent functional MRI allows to investigte neural activities and connectivity. While the non-human primate plays an essential role in neuroscience research, multimodal methods combining functional MRI with other neuroimaging and neuromodulation enable us to understand the brain network at multiple scales. APPROACH: In this study, a tight-fitting helmet-shape receive array with a single transmit loop for anesthetized macaque brain MRI at 7T was fabricated with four openings constructed in the coil housing to accommodate multimodal devices, and the coil performance was quantitatively evaluated and compared to a commercial knee coil. In addition, experiments over three macaques with infrared neural stimulation (INS), focused ultrasound stimulation (FUS), and transcranial direct current stimulation (tDCS) were conducted. MAIN RESULTS: The RF coil showed higher transmit efficiency, comparable homogeneity, improved SNR and enlarged signal coverage over the macaque brain. Infrared neural stimulation was applied to the amygdala in deep brain region, and activations in stimulation sites and connected sites were detected, with the connectivity consistent with anatomical information. Focused ultrasound stimulation was applied to the left visual cortex, and activations were acquired along the ultrasound traveling path, with all time course curves consistent with pre-designed paradigms. The existence of transcranial direct current stimulation electrodes brought no interference to the RF system, as evidenced through high-resolution MPRAGE structure images. SIGNIFICANCE: This pilot study reveals the feasibility for brain investigation at multiple spatiotemporal scales, which may advance our understanding in dynamic brain networks.

Critical factors in achieving fine-scale functional MRI: Removing sources of inadvertent spatial smoothing.

Ultra-high Field (≥7T) functional magnetic resonance imaging (UHF-fMRI) provides opportunities to resolve fine-scale features of functional architecture such as cerebral cortical columns and layers, in vivo. While the nominal resolution of modern fMRI acquisitions may appear to be sufficient to resolve these features, several common data preprocessing steps can introduce unwanted spatial blurring, especially those that require interpolation of the data. These resolution losses can impede the detection of the fine-scale features of interest. To examine quantitatively and systematically the sources of spatial resolution losses occurring during preprocessing, we used synthetic fMRI data and real fMRI data from the human visual cortex-the spatially interdigitated human V2 "thin" and "thick" stripes. The pattern of these cortical columns lies along the cortical surface and thus can be best appreciated using surface-based fMRI analysis. We used this as a testbed for evaluating strategies that can reduce spatial blurring of fMRI data. Our results show that resolution losses can be mitigated at multiple points in preprocessing pathway. We show that unwanted blur is introduced at each step of volume transformation and surface projection, and can be ameliorated by replacing multi-step transformations with equivalent single-step transformations. Surprisingly, the simple approaches of volume upsampling and of cortical mesh refinement also helped to reduce resolution losses caused by interpolation. Volume upsampling also serves to improve motion estimation accuracy, which helps to reduce blur. Moreover, we demonstrate that the level of spatial blurring is nonuniform over the brain-knowledge which is critical for interpreting data in high-resolution fMRI studies. Importantly, our study provides recommendations for reducing unwanted blurring during preprocessing as well as methods that enable quantitative comparisons between preprocessing strategies. These findings highlight several underappreciated sources of a spatial blur. Individually, the factors that contribute to spatial blur may appear to be minor, but in combination, the cumulative effects can hinder the interpretation of fine-scale fMRI and the detectability of these fine-scale features of functional architecture.

Infrared neural stimulation with 7T fMRI: A rapid in vivo method for mapping cortical connections of primate amygdala.

We have previously shown that INS-fMRI is a rapid method for mapping mesoscale brain networks in the macaque monkey brain. Focal stimulation of single cortical sites led to the activation of connected cortical locations, resulting in a global connectivity map. Here, we have extended this method for mapping brainwide networks following stimulation of single subcortical sites. As a testbed, we focused on the basal nucleus of the amygdala in the macaque monkey. We describe methods to target basal nucleus locations with submillimeter precision, pulse train stimulation methods, and statistical tests for assessing non-random nature of activations. Using these methods, we report that stimulation of precisely targeted loci in the basal nucleus produced sparse and specific activations in the brain. Activations were observed in the insular and sensory association cortices as well as activations in the cingulate cortex, consistent with known anatomical connections. What is new here is that the activations were focal and, in some cases, exhibited shifting topography with millimeter shifts in stimulation site. The precision of the method enables networks mapped from different nearby sites in the basal nucleus to be distinguished. While further investigation is needed to improve the sensitivity of this method, our analyses do support the reproducibility and non-random nature of some of the activations. We suggest that INS-fMRI is a promising method for mapping large-scale cortical and subcortical networks at high spatial resolution.

Combining brain perturbation and neuroimaging in non-human primates.

Brain perturbation studies allow detailed causal inferences of behavioral and neural processes. Because the combination of brain perturbation methods and neural measurement techniques is inherently challenging, research in humans has predominantly focused on non-invasive, indirect brain perturbations, or neurological lesion studies. Non-human primates have been indispensable as a neurobiological system that is highly similar to humans while simultaneously being more experimentally tractable, allowing visualization of the functional and structural impact of systematic brain perturbation. This review considers the state of the art in non-human primate brain perturbation with a focus on approaches that can be combined with neuroimaging. We consider both non-reversible (lesions) and reversible or temporary perturbations such as electrical, pharmacological, optical, optogenetic, chemogenetic, pathway-selective, and ultrasound based interference methods. Method-specific considerations from the research and development community are offered to facilitate research in this field and support further innovations. We conclude by identifying novel avenues for further research and innovation and by highlighting the clinical translational potential of the methods.

Focal Electrical Stimulation of Cortical Functional Networks.

Traditional electrical stimulation of brain tissue typically affects relatively large volumes of tissue spanning multiple millimeters. This low spatial resolution stimulation results in nonspecific functional effects. In addition, a primary shortcoming of these designs was the failure to take advantage of inherent functional organization in the cerebral cortex. Here, we describe a new method to electrically stimulate the brain which achieves selective targeting of single feature-specific domains in visual cortex. We provide evidence that this paradigm achieves mesoscale, functional network-specificity, and intensity dependence in a way that mimics visual stimulation. Application of this approach to known feature domains (such as color, orientation, motion, and depth) in visual cortex may lead to important functional improvements in the specificity and sophistication of brain stimulation methods and has implications for visual cortical prosthetic design.

Functionally specific optogenetic modulation in primate visual cortex.

In primates, visual perception is mediated by brain circuits composed of submillimeter nodes linked together in specific networks that process different types of information, such as eye specificity and contour orientation. We hypothesized that optogenetic stimulation targeted to cortical nodes could selectively activate such cortical networks. We used viral transfection methods to confer light sensitivity to neurons in monkey primary visual cortex. Using intrinsic signal optical imaging and single-unit electrophysiology to assess effects of targeted optogenetic stimulation, we found that (i) optogenetic stimulation of single ocular dominance columns (eye-specific nodes) revealed preferential activation of nearby same-eye columns but not opposite-eye columns, and (ii) optogenetic stimulation of single orientation domains increased visual response of matching orientation domains and relatively suppressed nonmatching orientation selectivity. These findings demonstrate that optical stimulation of single nodes leads to modulation of functionally specific cortical networks related to underlying neural architecture.

Optical imaging reveals functional domains in primate sensorimotor cortex.

Motor cortex (M1) and somatosensory cortex (S1) are central to arm and hand control. Efforts to understand encoding in M1 and S1 have focused on temporal relationships between neural activity and movement features. However, it remains unclear how the neural activity is spatially organized within M1 and S1. Optical imaging methods are well-suited for revealing the spatio-temporal organization of cortical activity, but their application is sparse in monkey sensorimotor cortex. Here, we investigate the effectiveness of intrinsic signal optical imaging (ISOI) for measuring cortical activity that supports arm and hand control in a macaque monkey. ISOI revealed spatial domains that were active in M1 and S1 in response to instructed reaching and grasping. The lateral M1 domains overlapped the hand representation and contained a population of neurons with peak firing during grasping. In contrast, the medial M1 domain overlapped the arm representation and a population of neurons with peak firing during reaching. The S1 domain overlapped the hand representations of areas 1 and 2 and a population of neurons with peak firing upon hand contact with the target. Our single unit recordings indicate that ISOI domains report the locations of spatial clusters of functionally related neurons. ISOI is therefore an effective tool for surveilling the neocortex for "hot zones" of activity that supports movement. Combining the strengths of ISOI with other imaging modalities (e.g., fMRI, 2-photon) and with electrophysiological methods can open new frontiers in understanding the spatio-temporal organization of cortical signals involved in movement control.

A 16-channel AC/DC array coil for anesthetized monkey whole-brain imaging at 7T.

Functional magnetic resonance imaging (fMRI) in monkeys is important for bridging the gap between invasive animal brain studies and non-invasive human brain studies. To resolve the finer functional structure of the monkey brain, ultra-high-field (UHF) MR is essential, and high-performance, close-fitting RF receive coils are typically desired to fully leverage the intrinsic gains provided by UHF MRI. Moreover, static field (B0) inhomogeneity arising from the tissue susceptibility interface is more severe at UHF, presenting an obstacle to achieving high-resolution fMRI. B0 shim of the monkey head is challenging due to its smaller size and more complex sources of B0 offsets in multi-modal imaging tasks. In the present work, we have customized an array coil for lightly-anesthetized monkey fMRI in the 7T human scanner that combines RF and multi-coil (MC) B0 shim functionality (also referred to as AC/DC coils) to provide high imaging SNR and high-spatial-order, rapidly switchable B0-shim capability. Additional space was retained on the coil to render it compatible with monkey multi-modal imaging studies. Both MC global (whole-volume) and dynamic (slice-optimized) shim methods were tested and evaluated, and the benefits of MC shim for fMRI experiments was also studied. A minor reduction in RF coil performance was found after introducing additional B0 shim circuitry. However, the proposed RF coil provided higher image SNR and more uniform contrast compared to a commercially available coil for human knee imaging. Compared with static 2nd-order shim, the B0 inhomogeneity was reduced by 56.8%, and 95-percentile B0 offset was reduced to within 28.2 Hz through MC shim, versus 68.7 Hz with 2nd-order static shim. As a result, functional image quality could be improved, and brain activation can be better detected using the proposed AC/DC monkey coil.

TNAP activity is localized at critical sites of retinal neurotransmission across various vertebrate species.

Evidence is emerging with regard to the role of tissue non-specific alkaline phosphatase (TNAP) in neural functions. As an ectophosphatase, this enzyme might influence neural activity and synaptic transmission in diverse ways. The localization of the enzyme in known neural circuits, such as the retina, might significantly advance an understanding of its role in normal and pathological functioning. However, the presence of TNAP in the retina is scarcely investigated. Our multispecies comparative study (zebrafish, cichlid, frog, chicken, mouse, rat, golden hamster, guinea pig, rabbit, sheep, cat, dog, ferret, squirrel monkey, human) using enzyme histochemistry and Western blots has shown the presence of TNAP activity in the retina of several mammalian species, including humans. Although the TNAP activity pattern varies across species, we have observed the following trends: (1) in all investigated species (except golden hamster), retinal vessels display TNAP activity; (2) TNAP activity consistently occurs in the photoreceptor layer; (3) in majority of the investigated species, marked TNAP activity is present in the outer and inner plexiform layers. In zebrafish, frog, chicken, guinea pig, and rat, TNAP histochemistry has revealed several sublayers of the inner plexiform layer. Frog, golden hamster, guinea pig, mouse, and human retinas possess a subpopulation of amacrine cells positively staining for TNAP activity. The expression of TNAP in critical sites of retinal signal transmission across a wide range of species suggests its fundamental, evolutionally conserved role in vision.

Quantifying tissue temperature changes induced by infrared neural stimulation: numerical simulation and MR thermometry.

Infrared neural stimulation (INS) delivered via short pulse trains is an innovative tool that has potential for us use for studying brain function and circuitry, brain machine interface, and clinical use. The prevailing mechanism for INS involves the conversion of light energy into thermal transients, leading to neuronal membrane depolarization. Due to the potential risks of thermal damage, it is crucial to ensure that the resulting local temperature increases are within non-damaging limits for brain tissues. Previous studies have estimated damage thresholds using histological methods and have modeled thermal effects based on peripheral nerves. However, additional quantitative measurements and modeling studies are needed for the central nervous system. Here, we performed 7 T MRI thermometry on ex vivo rat brains following the delivery of infrared pulse trains at five different intensities from 0.1-1.0 J/cm2 (each pulse train 1,875 nm, 25 us/pulse, 200 Hz, 0.5 s duration, delivered through 200 µm fiber). Additionally, we utilized the General BioHeat Transfer Model (GBHTM) to simulate local temperature changes in perfused brain tissues while delivering these laser energies to tissue (with optical parameters of human skin) via three different sizes of optical fibers at five energy intensities. The simulation results clearly demonstrate that a 0.5 second INS pulse train induces an increase followed by an immediate drop in temperature at stimulation offset. The delivery of multiple pulse trains with 2.5 s interstimulus interval (ISI) leads to rising temperatures that plateau. Both thermometry and modeling results show that, using parameters that are commonly used in biological applications (200 µm diameter fiber, 0.1-1.0 J/cm2), the final temperature increase at the end of the 60 sec stimuli duration does not exceed 1°C with stimulation values of 0.1-0.5 J/cm2 and does not exceed 2°C with stimulation values of up to 1.0 J/cm2. Thus, the maximum temperature rise is consistent with the thermal damage threshold reported in previous studies. This study provides a quantitative evaluation of the temperature changes induced by INS, suggesting that existing practices pose minimal major safety concerns for biological tissues.

A novel interface for cortical columnar neuromodulation with multipoint infrared neural stimulation.

Cutting edge advances in electrical visual cortical prosthetics have evoked perception of shapes, motion, and letters in the blind. Here, we present an alternative optical approach using pulsed infrared neural stimulation. To interface with dense arrays of cortical columns with submillimeter spatial precision, both linear array and 100-fiber bundle array optical fiber interfaces were devised. We deliver infrared stimulation through these arrays in anesthetized cat visual cortex and monitor effects by optical imaging in contralateral visual cortex. Infrared neural stimulation modulation of response to ongoing visual oriented gratings produce enhanced responses in orientation-matched domains and suppressed responses in non-matched domains, consistent with a known higher order integration mediated by callosal inputs. Controls include dynamically applied speeds, directions and patterns of multipoint stimulation. This provides groundwork for a distinct type of prosthetic targeted to maps of visual cortical columns.

Mesoscale organization of ventral and dorsal visual pathways in macaque monkey revealed by 7T fMRI.

In human and nonhuman primate brains, columnar (mesoscale) organization has been demonstrated to underlie both lower and higher order aspects of visual information processing. Previous studies have focused on identifying functional preferences of mesoscale domains in specific areas; but there has been little understanding of how mesoscale domains may cooperatively respond to single visual stimuli across dorsal and ventral pathways. Here, we have developed ultrahigh-field 7 T fMRI methods to enable simultaneous mapping, in individual macaque monkeys, of response in both dorsal and ventral pathways to single simple color and motion stimuli. We provide the first evidence that anatomical V2 cytochrome oxidase-stained stripes are well aligned with fMRI maps of V2 stripes, settling a long-standing controversy. In the ventral pathway, a systematic array of paired color and luminance processing domains across V4 was revealed, suggesting a novel organization for surface information processing. In the dorsal pathway, in addition to high quality motion direction maps of MT, MST and V3A, alternating color and motion direction domains in V3 are revealed. As well, submillimeter motion domains were observed in peripheral LIPd and LIPv. In sum, our study provides a novel global snapshot of how mesoscale networks in the ventral and dorsal visual pathways form the organizational basis of visual objection recognition and vision for action.

Orientation selectivity mapping in the visual cortex.

The orientation map is one of the most well-studied functional maps of the visual cortex. However, results from the literature are of different qualities. Clear boundaries among different orientation domains and blurred uncertain distinctions were shown in different studies. These unclear imaging results will lead to an inaccuracy in depicting cortical structures, and the lack of consideration in experimental design will also lead to biased depictions of the cortical features. How we accurately define orientation domains will impact the entire field of research. In this study, we test how spatial frequency (SF), stimulus size, location, chromatic, and data processing methods affect the orientation functional maps (including a large area of dorsal V4, and parts of dorsal V1) acquired by intrinsic signal optical imaging. Our results indicate that, for large imaging fields, large grating stimuli with mixed SF components should be considered to acquire the orientation map. A diffusion model image enhancement based on the difference map could further improve the map quality. In addition, the similar outcomes of achromatic and chromatic gratings indicate two alternative types of afferents from LGN, pooling in V1 to generate cue-invariant orientation selectivity.

Two-photon imaging of excitatory and inhibitory neural response to infrared neural stimulation.

Significance: Pulsed infrared neural stimulation (INS, 1875 nm) is an emerging neurostimulation technology that delivers focal pulsed heat to activate functionally specific mesoscale networks and holds promise for clinical application. However, little is known about its effect on excitatory and inhibitory cell types in cerebral cortex. Aim: Estimates of summed population neuronal response time courses provide a potential basis for neural and hemodynamic signals described in other studies. Approach: Using two-photon calcium imaging in mouse somatosensory cortex, we have examined the effect of INS pulse train application on hSyn neurons and mDlx neurons tagged with GCaMP6s. Results: We find that, in anesthetized mice, each INS pulse train reliably induces robust response in hSyn neurons exhibiting positive going responses. Surprisingly, mDlx neurons exhibit negative going responses. Quantification using the index of correlation illustrates responses are reproducible, intensity-dependent, and focal. Also, a contralateral activation is observed when INS applied. Conclusions: In sum, the population of neurons stimulated by INS includes both hSyn and mDlx neurons; within a range of stimulation intensities, this leads to overall excitation in the stimulated population, leading to the previously observed activations at distant post-synaptic sites.

Brainwide mesoscale functional networks revealed by focal infrared neural stimulation of the amygdala.

The primate amygdala serves to evaluate emotional content of sensory inputs and modulate emotional and social behaviors; it modulates cognitive, multisensory and autonomic circuits predominantly via the basal (BA), lateral (LA), and central (CeA) nuclei, respectively. Based on recent electrophysiological evidence suggesting mesoscale (millimeters-scale) nature of intra-amygdala functional organization, we have investigated the connectivity of these nuclei using Infrared Neural Stimulation of single mesoscale sites coupled with mapping in ultrahigh field 7T functional Magnetic Resonance Imaging (INS-fMRI). Stimulation of multiple sites within amygdala of single individuals evoked 'mesoscale functional connectivity maps', allowing comparison of BA, LA and CeA connected brainwide networks. This revealed a mesoscale nature of connected sites, complementary spatial patterns of functional connectivity, and topographic relationships of nucleus-specific connections. Our data reveal a functional architecture of systematically organized brainwide networks mediating sensory, cognitive, and autonomic influences from the amygdala.

Motor versus Psychomotor? Deciphering the Neural Source of Psychomotor Retardation in Depression.

Major depressive disorder (MDD) is characterized by psychomotor retardation whose underlying neural source remains unclear. Psychomotor retardation may either be related to a motor source like the motor cortex or, alternatively, to a psychomotor source with neural changes outside motor regions, like input regions such as visual cortex. These two alternative hypotheses in main (n = 41) and replication (n = 18) MDD samples using 7 Tesla MRI are investigated. Analyzing both global and local connectivity in primary motor cortex (BA4), motor network and middle temporal visual cortex complex (MT+), the main findings in MDD are: 1) Reduced local and global synchronization and increased local-to-global output in motor regions, which do not correlate with psychomotor retardation, though. 2) Reduced local-to-local BA4 - MT+ functional connectivity (FC) which correlates with psychomotor retardation. 3) Reduced global synchronization and increased local-to-global output in MT+ which relate to psychomotor retardation. 4) Reduced variability in the psychophysical measures of MT+ based motion perception which relates to psychomotor retardation. Together, it is shown that visual cortex MT+ and its relation to motor cortex play a key role in mediating psychomotor retardation. This supports psychomotor over motor hypothesis about the neural source of psychomotor retardation in MDD.

BMI 2.0: Toward a technological interface with brainwide networks.

The field of brain machine interface has long sought a technology for brainwide interaction. In this issue of Neuron, Kim et al.1 present a novel method for dynamic, patterned, and precise optogenetic stimulation of mouse cortex in ultra-high-field MRI that portends such an interface.

Spatial frequency representation in V2 and V4 of macaque monkey.

Spatial frequency (SF) is an important attribute in the visual scene and is a defining feature of visual processing channels. However, there remain many unsolved questions about how extrastriate areas in primate visual cortex code this fundamental information. Here, using intrinsic signal optical imaging in visual areas of V2 and V4 of macaque monkeys, we quantify the relationship between SF maps and (1) visual topography and (2) color and orientation maps. We find that in orientation regions, low to high SF is mapped orthogonally to orientation; in color regions, which are reported to contain orthogonal axes of color and lightness, low SFs tend to be represented more frequently than high SFs. This supports a population-based SF fluctuation related to the 'color/orientation' organizations. We propose a generalized hypercolumn model across cortical areas, comprised of two orthogonal parameters with additional parameters.

Targeted Optical Neural Stimulation: A New Era for Personalized Medicine.

Targeted optical neural stimulation comprises infrared neural stimulation and optogenetics, which affect the nervous system through induced thermal transients and activation of light-sensitive proteins, respectively. The main advantage of this pair of optical tools is high functional selectivity, which conventional electrical stimulation lacks. Over the past 15 years, the mechanism, safety, and feasibility of optical stimulation techniques have undergone continuous investigation and development. When combined with other methods like optical imaging and high-field functional magnetic resonance imaging, the translation of optical stimulation to clinical practice adds high value. We review the theoretical foundations and current state of optical stimulation, with a particular focus on infrared neural stimulation as a potential bridge linking optical stimulation to personalized medicine.

Infrared neural stimulation in human cerebral cortex.

BACKGROUND: Modulation of brain circuits by electrical stimulation has led to exciting and powerful therapies for diseases such as Parkinson's. Because human brain organization is based in mesoscale (millimeter-scale) functional nodes, having a method that can selectively target such nodes could enable more precise, functionally specific stimulation therapies. Infrared Neural Stimulation (INS) is an emerging stimulation technology that stimulates neural tissue via delivery of tiny heat pulses. In nonhuman primates, this optical method provides focal intensity-dependent stimulation of the brain without tissue damage. However, whether INS application to the human central nervous system (CNS) is similarly effective is unknown. OBJECTIVE: To examine the effectiveness of INS on human cerebral cortex in intraoperative setting and to evaluate INS damage threshholds. METHODS: Five epileptic subjects undergoing standard lobectomy for epilepsy consented to this study. Cortical response to INS was assessed by intrinsic signal optical imaging (OI, a method that detects changes in tissue reflectance due to neuronal activity). A custom integrated INS and OI system was developed specifically for short-duration INS and OI acquisition during surgical procedures. Single pulse trains of INS with intensities from 0.2 to 0.8 J/cm2 were delivered to the somatosensory cortex and responses were recorded via optical imaging. Following tissue resection, histological analysis was conducted to evaluate damage threshholds. RESULTS: As assessed by OI, and similar to results in monkeys, INS induced responses in human cortex were highly focal (millimeter sized) and led to relative suppression of nearby cortical sites. Intensity dependence was observed at both stimulated and functionally connected sites. Histological analysis of INS-stimulated human cortical tissue provided damage threshold estimates. CONCLUSION: This is the first study demonstrating application of INS to human CNS and shows feasibility for stimulating single cortical nodes and associated sites and provided INS damage threshold estimates for cortical tissue. Our results suggest that INS is a promising tool for stimulation of functionally selective mesoscale circuits in the human brain, and may lead to advances in the future of precision medicine.