Dissociation of Centrally and Peripherally Induced Transcranial Magnetic Stimulation Effects in Nonhuman Primates.

Transcranial magnetic stimulation (TMS) is a noninvasive brain stimulation method that is rapidly growing in popularity for studying causal brain-behavior relationships. However, its dose-dependent centrally induced neural mechanisms and peripherally induced sensory costimulation effects remain debated. Understanding how TMS stimulation parameters affect brain responses is vital for the rational design of TMS protocols. Studying these mechanisms in humans is challenging because of the limited spatiotemporal resolution of available noninvasive neuroimaging methods. Here, we leverage invasive recordings of local field potentials in a male and a female nonhuman primate (rhesus macaque) to study TMS mesoscale responses. We demonstrate that early TMS-evoked potentials show a sigmoidal dose-response curve with stimulation intensity. We further show that stimulation responses are spatially specific. We use several control conditions to dissociate centrally induced neural responses from auditory and somatosensory coactivation. These results provide crucial evidence regarding TMS neural effects at the brain circuit level. Our findings are highly relevant for interpreting human TMS studies and biomarker developments for TMS target engagement in clinical applications.SIGNIFICANCE STATEMENT Transcranial magnetic stimulation (TMS) is a widely used noninvasive brain stimulation method to stimulate the human brain. To advance its utility for clinical applications, a clear understanding of its underlying physiological mechanisms is crucial. Here, we perform invasive electrophysiological recordings in the nonhuman primate brain during TMS, achieving a spatiotemporal precision not available in human EEG experiments. We find that evoked potentials are dose dependent and spatially specific, and can be separated from peripheral stimulation effects. This means that TMS-evoked responses can indicate a direct physiological stimulation response. Our work has important implications for the interpretation of human TMS-EEG recordings and biomarker development.

Anatomical details affect electric field predictions for non-invasive brain stimulation in non-human primates.

Non-human primates (NHPs) have become key for translational research in noninvasive brain stimulation (NIBS). However, in order to create comparable stimulation conditions for humans it is vital to study the accuracy of current modeling practices across species. Numerical models to simulate electric fields are an important tool for experimental planning in NHPs and translation to human studies. It is thus essential whether and to what extent the anatomical details of NHP models agree with current modeling practices when calculating NIBS electric fields. Here, we create highly accurate head models of two non-human primates (NHP) MR data. We evaluate how muscle tissue and head field of view (depending on MRI parameters) affect simulation results in transcranial electric and magnetic stimulation (TES and TMS). Our findings indicate that the inclusion of anisotropic muscle can affect TES electric field strength up to 22% while TMS is largely unaffected. Additionally, comparing a full head model to a cropped head model illustrates the impact of head field of view on electric fields for both TES and TMS. We find opposing effects between TES and TMS with an increase up to 24.8% for TES and a decrease up to 24.6% for TMS for the cropped head model compared to the full head model. Our results provide important insights into the level of anatomical detail needed for NHP head models and can inform future translational efforts for NIBS studies.

Top-down, contextual entrainment of neuronal oscillations in the auditory thalamocortical circuit.

Prior studies have shown that repetitive presentation of acoustic stimuli results in an alignment of ongoing neuronal oscillations to the sequence rhythm via oscillatory entrainment by external cues. Our study aimed to explore the neural correlates of the perceptual parsing and grouping of complex repeating auditory patterns that occur based solely on statistical regularities, or context. Human psychophysical studies suggest that the recognition of novel auditory patterns amid a continuous auditory stimulus sequence occurs automatically halfway through the first repetition. We hypothesized that once repeating patterns were detected by the brain, internal rhythms would become entrained, demarcating the temporal structure of these repetitions despite lacking external cues defining pattern on- or offsets. To examine the neural correlates of pattern perception, neuroelectric activity of primary auditory cortex (A1) and thalamic nuclei was recorded while nonhuman primates passively listened to streams of rapidly presented pure tones and bandpass noise bursts. At arbitrary intervals, random acoustic patterns composed of 11 stimuli were repeated five times without any perturbance of the constant stimulus flow. We found significant delta entrainment by these patterns in the A1, medial geniculate body, and medial pulvinar. In A1 and pulvinar, we observed a statistically significant, pattern structure-aligned modulation of neuronal firing that occurred earliest in the pulvinar, supporting the idea that grouping and detecting complex auditory patterns is a top-down, context-driven process. Besides electrophysiological measures, a pattern-related modulation of pupil diameter verified that, like humans, nonhuman primates consciously detect complex repetitive patterns that lack physical boundaries.

Limitations of ex vivo measurements for in vivo neuroscience.

A long history of postmortem studies has provided significant insight into human brain structure and organization. Cadavers have also proven instrumental for the measurement of artifacts and nonneural effects in functional imaging, and more recently, the study of biophysical properties critical to brain stimulation. However, death produces significant changes in the biophysical properties of brain tissues, making an ex vivo to in vivo comparison complex, and even questionable. This study directly compares biophysical properties of electric fields arising from transcranial electric stimulation (TES) in a nonhuman primate brain pre- and postmortem. We show that pre- vs. postmortem, TES-induced intracranial electric fields differ significantly in both strength and frequency response dynamics, even while controlling for confounding factors such as body temperature. Our results clearly indicate that ex vivo cadaver and in vivo measurements are not easily equitable. In vivo examinations remain essential to establishing an adequate understanding of even basic biophysical phenomena in vivo.

Auditory properties in the parabelt regions of the superior temporal gyrus in the awake macaque monkey: an initial survey.

The superior temporal gyrus (STG) is on the inferior-lateral brain surface near the external ear. In macaques, 2/3 of the STG is occupied by an auditory cortical region, the "parabelt," which is part of a network of inferior temporal areas subserving communication and social cognition as well as object recognition and other functions. However, due to its location beneath the squamous temporal bone and temporalis muscle, the STG, like other inferior temporal regions, has been a challenging target for physiological studies in awake-behaving macaques. We designed a new procedure for implanting recording chambers to provide direct access to the STG, allowing us to evaluate neuronal properties and their topography across the full extent of the STG in awake-behaving macaques. Initial surveys of the STG have yielded several new findings. Unexpectedly, STG sites in monkeys that were listening passively responded to tones with magnitudes comparable to those of responses to 1/3 octave band-pass noise. Mapping results showed longer response latencies in more rostral sites and possible tonotopic patterns parallel to core and belt areas, suggesting the reversal of gradients between caudal and rostral parabelt areas. These results will help further exploration of parabelt areas.

Experimental validation of computational models for the prediction of phase distribution during multi-channel transcranial alternating current stimulation.

Transcranial alternating current stimulation (tACS) is a widely used noninvasive brain stimulation (NIBS) technique to affect neural activity. Neural oscillations exhibit phase-dependent associations with cognitive functions, and tools to manipulate local oscillatory phases can affect communication across remote brain regions. A recent study demonstrated that multi-channel tACS can generate electric fields with a phase gradient or traveling waves in the brain. Computational simulations using phasor algebra can predict the phase distribution inside the brain and aid in informing parameters in tACS experiments. However, experimental validation of computational models for multi-phase tACS is still lacking. Here, we develop such a framework for phasor simulation and evaluate its accuracy using in vivo recordings in nonhuman primates. We extract the phase and amplitude of electric fields from intracranial recordings in two monkeys during multi-channel tACS and compare them to those calculated by phasor analysis using finite element models. Our findings demonstrate that simulated phases correspond well to measured phases (r = 0.9). Further, we systematically evaluated the impact of accurate electrode placement on modeling and data agreement. Finally, our framework can predict the amplitude distribution in measurements given calibrated tissues’ conductivity. Our validated general framework for simulating multi-phase, multi-electrode tACS provides a streamlined tool for principled planning of multi-channel tACS experiments.

Predicting the phase distribution during multi-channel transcranial alternating current stimulation in silico and in vivo.

BACKGROUND: Transcranial alternating current stimulation (tACS) is a widely used noninvasive brain stimulation (NIBS) technique to affect neural activity. TACS experiments have been coupled with computational simulations to predict the electromagnetic fields within the brain. However, existing simulations are focused on the magnitude of the field. As the possibility of inducing the phase gradient in the brain using multiple tACS electrodes arises, a simulation framework is necessary to investigate and predict the phase gradient of electric fields during multi-channel tACS. OBJECTIVE: Here, we develop such a framework for phasor simulation using phasor algebra and evaluate its accuracy using in vivo recordings in monkeys. METHODS: We extract the phase and amplitude of electric fields from intracranial recordings in two monkeys during multi-channel tACS and compare them to those calculated by phasor analysis using finite element models. RESULTS: Our findings demonstrate that simulated phases correspond well to measured phases (r = 0.9). Further, we systematically evaluated the impact of accurate electrode placement on modeling and data agreement. Finally, our framework can predict the amplitude distribution in measurements given calibrated tissues' conductivity. CONCLUSIONS: Our validated general framework for simulating multi-phase, multi-electrode tACS provides a streamlined tool for principled planning of multi-channel tACS experiments.

An open-access dataset of naturalistic viewing using simultaneous EEG-fMRI.

In this work, we present a dataset that combines functional magnetic imaging (fMRI) and electroencephalography (EEG) to use as a resource for understanding human brain function in these two imaging modalities. The dataset can also be used for optimizing preprocessing methods for simultaneously collected imaging data. The dataset includes simultaneously collected recordings from 22 individuals (ages: 23-51) across various visual and naturalistic stimuli. In addition, physiological, eye tracking, electrocardiography, and cognitive and behavioral data were collected along with this neuroimaging data. Visual tasks include a flickering checkerboard collected outside and inside the MRI scanner (EEG-only) and simultaneous EEG-fMRI recordings. Simultaneous recordings include rest, the visual paradigm Inscapes, and several short video movies representing naturalistic stimuli. Raw and preprocessed data are openly available to download. We present this dataset as part of an effort to provide open-access data to increase the opportunity for discoveries and understanding of the human brain and evaluate the correlation between electrical brain activity and blood oxygen level-dependent (BOLD) signals.

Projection from visual areas V2 and prostriata to caudal auditory cortex in the monkey.

Studies in humans and monkeys report widespread multisensory interactions at or near primary visual and auditory areas of neocortex. The range and scale of these effects has prompted increased interest in interconnectivity between the putatively "unisensory" cortices at lower hierarchical levels. Recent anatomical tract-tracing studies have revealed direct projections from auditory cortex to primary visual area (V1) and secondary visual area (V2) that could serve as a substrate for auditory influences over low-level visual processing. To better understand the significance of these connections, we looked for reciprocal projections from visual cortex to caudal auditory cortical areas in macaque monkeys. We found direct projections from area prostriata and the peripheral visual representations of area V2. Projections were more abundant after injections of temporoparietal area and caudal parabelt than after injections of caudal medial belt and the contiguous areas near the fundus of the lateral sulcus. Only one injection was confined to primary auditory cortex (area A1) and did not demonstrate visual connections. The projections from visual areas originated mainly from infragranular layers, suggestive of a "feedback"-type projection. The selective localization of these connections to peripheral visual areas and caudal auditory cortex suggests that they are involved in spatial localization.

Dissociation of broadband high-frequency activity and neuronal firing in the neocortex.

Broadband high-frequency activity (BHA; 70 to 150 Hz), also known as "high gamma," a key analytic signal in human intracranial (electrocorticographic) recordings, is often assumed to reflect local neural firing [multiunit activity (MUA)]. As the precise physiological substrates of BHA are unknown, this assumption remains controversial. Our analysis of laminar multielectrode data from V1 and A1 in monkeys outlines two components of stimulus-evoked BHA distributed across the cortical layers: an "early-deep" and "late-superficial" response. Early-deep BHA has a clear spatial and temporal overlap with MUA. Late-superficial BHA was more prominent and accounted for more of the BHA signal measured near the cortical pial surface. However, its association with local MUA is weak and often undetectable, consistent with the view that it reflects dendritic processes separable from local neuronal firing.

Gradients of connectivity distance in the cerebral cortex of the macaque monkey.

Cortical connectivity conforms to a series of organizing principles that are common across species. Spatial proximity, similar cortical type, and similar connectional profile all constitute factors for determining the connectivity between cortical regions. We previously demonstrated another principle of connectivity that is closely related to the spatial layout of the cerebral cortex. Using functional connectivity from resting-state fMRI in the human cortex, we found that the further a region is located from primary cortex, the more distant are its functional connections with the other areas of the cortex. However, it remains unknown whether this relationship between cortical layout and connectivity extends to other primate species. Here, we investigated this relationship using both resting-state functional connectivity as well as gold-standard tract-tracing connectivity in the macaque monkey cortex. For both measures of connectivity, we found a gradient of connectivity distance extending between primary and frontoparietal regions. In the human cortex, the further a region is located from primary areas, the stronger its connections to distant portions of the cortex, with connectivity distance highest in frontal and parietal regions. The similarity between the human and macaque findings provides evidence for a phylogenetically conserved relationship between the spatial layout of cortical areas and connectivity.

Electric field dynamics in the brain during multi-electrode transcranial electric stimulation.

Neural oscillations play a crucial role in communication between remote brain areas. Transcranial electric stimulation with alternating currents (TACS) can manipulate these brain oscillations in a non-invasive manner. Recently, TACS using multiple electrodes with phase shifted stimulation currents were developed to alter long-range connectivity. Typically, an increase in coordination between two areas is assumed when they experience an in-phase stimulation and a disorganization through an anti-phase stimulation. However, the underlying biophysics of multi-electrode TACS has not been studied in detail. Here, we leverage direct invasive recordings from two non-human primates during multi-electrode TACS to characterize electric field magnitude and phase as a function of the phase of stimulation currents. Further, we report a novel "traveling wave" stimulation where the location of the electric field maximum changes over the stimulation cycle. Our results provide a mechanistic understanding of the biophysics of multi-electrode TACS and enable future developments of novel stimulation protocols.

Delineating the Macroscale Areal Organization of the Macaque Cortex In Vivo.

Complementing long-standing traditions centered on histology, fMRI approaches are rapidly maturing in delineating brain areal organization at the macroscale. The non-human primate (NHP) provides the opportunity to overcome critical barriers in translational research. Here, we establish the data requirements for achieving reproducible and internally valid parcellations in individuals. We demonstrate that functional boundaries serve as a functional fingerprint of the individual animals and can be achieved under anesthesia or awake conditions (rest, naturalistic viewing), though differences between awake and anesthetized states precluded the detection of individual differences across states. Comparison of awake and anesthetized states suggested a more nuanced picture of changes in connectivity for higher-order association areas, as well as visual and motor cortex. These results establish feasibility and data requirements for the generation of reproducible individual-specific parcellations in NHPs, provide insights into the impact of scan state, and motivate efforts toward harmonizing protocols.

Spatiotemporal structure of intracranial electric fields induced by transcranial electric stimulation in humans and nonhuman primates.

Transcranial electric stimulation (TES) is an emerging technique, developed to non-invasively modulate brain function. However, the spatiotemporal distribution of the intracranial electric fields induced by TES remains poorly understood. In particular, it is unclear how much current actually reaches the brain, and how it distributes across the brain. Lack of this basic information precludes a firm mechanistic understanding of TES effects. In this study we directly measure the spatial and temporal characteristics of the electric field generated by TES using stereotactic EEG (s-EEG) electrode arrays implanted in cebus monkeys and surgical epilepsy patients. We found a small frequency dependent decrease (10%) in magnitudes of TES induced potentials and negligible phase shifts over space. Electric field strengths were strongest in superficial brain regions with maximum values of about 0.5 mV/mm. Our results provide crucial information of the underlying biophysics in TES applications in humans and the optimization and design of TES stimulation protocols. In addition, our findings have broad implications concerning electric field propagation in non-invasive recording techniques such as EEG/MEG.

Anatomical evidence of multimodal integration in primate striate cortex.

The primary visual cortex (area 17 or V1) is not thought to receive input from nonvisual extrastriate cortical areas. However, this has yet to be shown to be the case using sensitive tracers in the part of area 17 subserving the peripheral visual field. Here we show using retrograde tracers that peripheral area 17 subserving the visual field at an eccentricity of 10-20 degrees receives projections from the core and parabelt areas of the auditory cortex as well as from the polysensory area of the temporal lobe (STP). The relative strength of these projections was calculated for each injection by computing the proportions of retrogradely labeled neurons located in the auditory and STP areas with respect to number of labeled neurons constituting the established projection from the superior temporal sulci (STS) motion complex (middle temporal area, medial superior temporal, fundus of the superior temporal area). In peripheral area V1 the projection from auditory cortex corresponds to 9.5% of that of the STS motion complex and STP to 35% of that from the STS motion complex. Compared to peripheral area 17, central and paracentral area 17 showed considerably weaker inputs from auditory cortex (0.2-0.8%) but slightly more from STP cortex (3.5-6.1%). The present results show that the connectivity of area 17 is eccentricity dependent. Direct projections from auditory and STP cortex to peripheral area 17 have important consequences for higher visual functions of area 17, including multimodal integration at early stages of the visual cortical pathway.

Quantitative analysis of connectivity in the visual cortex: extracting function from structure.

It is generally agreed that information flow through the cortex is constrained by a hierarchical architecture. Lack of precise data on areal connectivity leads to indeterminacy of existing models. The authors introduce two quantitative parameters (SLN and FLN) that hold the promise of resolving such indeterminacy. In the visual system, using a very incomplete database, provisional hierarchies are in line with the recent proposal of higher functions of area V1 and suggest a hitherto unsuspected central function of the frontal eye field.

Feedforward and feedback projections of caudal belt and parabelt areas of auditory cortex: refining the hierarchical model.

Our working model of the primate auditory cortex recognizes three major regions (core, belt, parabelt), subdivided into thirteen areas. The connections between areas are topographically ordered in a manner consistent with information flow along two major anatomical axes: core-belt-parabelt and caudal-rostral. Remarkably, most of the connections supporting this model were revealed using retrograde tracing techniques. Little is known about laminar circuitry, as anterograde tracing of axon terminations has rarely been used. The purpose of the present study was to examine the laminar projections of three areas of auditory cortex, pursuant to analysis of all areas. The selected areas were: middle lateral belt (ML); caudomedial belt (CM); and caudal parabelt (CPB). Injections of anterograde tracers yielded data consistent with major features of our model, and also new findings that compel modifications. Results supporting the model were: (1) feedforward projection from ML and CM terminated in CPB; (2) feedforward projections from ML and CPB terminated in rostral areas of the belt and parabelt; and (3) feedback projections typified inputs to the core region from belt and parabelt. At odds with the model was the convergence of feedforward inputs into rostral medial belt from ML and CPB. This was unexpected since CPB is at a higher stage of the processing hierarchy, with mainly feedback projections to all other belt areas. Lastly, extending the model, feedforward projections from CM, ML, and CPB overlapped in the temporal parietal occipital area (TPO) in the superior temporal sulcus, indicating significant auditory influence on sensory processing in this region. The combined results refine our working model and highlight the need to complete studies of the laminar inputs to all areas of auditory cortex. Their documentation is essential for developing informed hypotheses about the neurophysiological influences of inputs to each layer and area.

Anatomy of hierarchy: feedforward and feedback pathways in macaque visual cortex.

The laminar location of the cell bodies and terminals of interareal connections determines the hierarchical structural organization of the cortex and has been intensively studied. However, we still have only a rudimentary understanding of the connectional principles of feedforward (FF) and feedback (FB) pathways. Quantitative analysis of retrograde tracers was used to extend the notion that the laminar distribution of neurons interconnecting visual areas provides an index of hierarchical distance (percentage of supragranular labeled neurons [SLN]). We show that: 1) SLN values constrain models of cortical hierarchy, revealing previously unsuspected areal relations; 2) SLN reflects the operation of a combinatorial distance rule acting differentially on sets of connections between areas; 3) Supragranular layers contain highly segregated bottom-up and top-down streams, both of which exhibit point-to-point connectivity. This contrasts with the infragranular layers, which contain diffuse bottom-up and top-down streams; 4) Cell filling of the parent neurons of FF and FB pathways provides further evidence of compartmentalization; 5) FF pathways have higher weights, cross fewer hierarchical levels, and are less numerous than FB pathways. Taken together, the present results suggest that cortical hierarchies are built from supra- and infragranular counterstreams. This compartmentalized dual counterstream organization allows point-to-point connectivity in both bottom-up and top-down directions.

The spectrotemporal filter mechanism of auditory selective attention.

Although we have convincing evidence that attention to auditory stimuli modulates neuronal responses at or before the level of primary auditory cortex (A1), the underlying physiological mechanisms are unknown. We found that attending to rhythmic auditory streams resulted in the entrainment of ongoing oscillatory activity reflecting rhythmic excitability fluctuations in A1. Strikingly, although the rhythm of the entrained oscillations in A1 neuronal ensembles reflected the temporal structure of the attended stream, the phase depended on the attended frequency content. Counter-phase entrainment across differently tuned A1 regions resulted in both the amplification and sharpening of responses at attended time points, in essence acting as a spectrotemporal filter mechanism. Our data suggest that selective attention generates a dynamically evolving model of attended auditory stimulus streams in the form of modulatory subthreshold oscillations across tonotopically organized neuronal ensembles in A1 that enhances the representation of attended stimuli.

Dual mechanism of neuronal ensemble inhibition in primary auditory cortex.

Inhibition plays an essential role in shaping and refining the brain's representation of sensory stimulus attributes. In primary auditory cortex (A1), so-called "sideband" inhibition helps to sharpen the tuning of local neuronal responses. Several distinct types of anatomical circuitry could underlie sideband inhibition, including direct thalamocortical (TC) afferents, as well as indirect intracortical mechanisms. The goal of the present study was to characterize sideband inhibition in A1 and to determine its mechanism by analyzing laminar profiles of neuronal ensemble activity. Our results indicate that both lemniscal and nonlemniscal TC afferents play a role in inhibitory responses via feedforward inhibition and oscillatory phase reset, respectively. We propose that the dynamic modulation of excitability in A1 due to the phase reset of ongoing oscillations may alter the tuning of local neuronal ensembles and can be regarded as a flexible overlay on the more obligatory system of lemniscal feedforward type responses.