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

Correspondence between thalamic injury-induced changes in resting-state fMRI of monkeys and their sensorimotor behaviors and neural activities

Resting state functional MRI (rsfMRI) exploits variations in blood-oxygenation-level-dependent (BOLD) signals to infer resting state functional connectivity (FC) within and between brain networks. However, there have been few reports quantifying and validating the results of rsfMRI analyses with other metrics of brain circuits. We measured longitudinal changes in FC both within and between brain networks in three squirrel monkeys after focal lesions of the thalamic ventroposterior lateral nucleus (VPL) that were intended to disrupt the input to somatosensory cortex and impair manual dexterity. Local field potential signals were recorded to assess electrophysiological changes during each animal's recovery, and behavioral performances were measured longitudinally using a sugar-pellet grasping task. Finally, end-point histological evaluations were performed on brain tissue slices to quantify the VPL damage. The rsfMRI data analysis showed significant decrease in FC measures both within and between networks immediately post-injury, which started to recover at different time-points for each animal. The trajectories of FC recovery for each animal mirrored their individual behavioral recovery time-courses. Electrophysiological measurements of inter-electrode coherences and end-point histological measures also aligned well with the graded injury effects measured using rsfMRI-based FC. A simple algorithm employing FC measures from the somatosensory network could accurately predict each monkeys' behavioral recovery timeframe after four weeks post-injury. Whole brain between-network FC measures further revealed that the injury effects were not limited to thalamocortical connections but were rather more widespread. Overall, this study provides evidence of the validity of rsfMRI based FC measures as indicators of the functional integrity and behavioral relevance following an injury to a specific brain circuit.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and the underlying patterns of neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the amygdala of two male macaques. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5 Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-modal approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the macaque amygdala and activated them with a highly selective and potent DREADD agonist, deschloroclozapine. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Interestingly, activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-disciplinary approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

Detection and characterization of resting state functional networks in squirrel monkey brain

Resting-state fMRI based on analyzing BOLD signals is widely used to derive functional networks in the brain and how they alter during disease or injury conditions. Resting-state networks can also be used to study brain functional connectomes across species, which provides insights into brain evolution. The squirrel monkey (SM) is a non-human primate (NHP) that is widely used as a preclinical model for experimental manipulations to understand the organization and functioning of the brain. We derived resting-state networks from the whole brain of anesthetized SMs using Independent Component Analysis of BOLD acquisitions. We detected 15 anatomically constrained resting-state networks localized in the cortical and subcortical regions as well as in the white-matter. Networks encompassing visual, somatosensory, executive control, sensorimotor, salience and default mode regions, and subcortical networks including the Hippocampus-Amygdala, thalamus, basal-ganglia and brainstem region correspond well with previously detected networks in humans and NHPs. The connectivity pattern between the networks also agrees well with previously reported seed-based resting-state connectivity of SM brain. This study demonstrates that SMs share remarkable homologous network organization with humans and other NHPs, thereby providing strong support for their suitability as a translational animal model for research and additional insight into brain evolution across species.

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.

Localizing Spectral Interactions in the Resting State Network Using the Hilbert-Huang Transform

Brain synchronizations are orchestrated from neuronal oscillations through frequency interactions, such as the alpha rhythm during relaxation. Nevertheless, how the intrinsic interaction forges functional integrity across brain segregations remains elusive, thereby motivating recent studies to localize frequency interactions of resting-state fMRI (rs-fMRI). To this point, we aim to unveil the fMRI-based spectral interactions using the time-frequency (TF) analysis; however, Fourier-based TF analyses impose restrictions on revealing frequency interactions given the limited time points in fMRI signals. Instead of using the Fourier-based wavelet analysis to identify the fMRI frequency of interests, we employed the Hilbert-Huang transform (HHT) for probing the specific frequency contribution to the functional integration, called ensemble spectral interaction (ESI). By simulating data with time-variant frequency changes, we demonstrated the Hilbert TF maps with high spectro-temporal resolution and full accessibility in comparison with the wavelet TF maps. By detecting amplitude-to-amplitude frequency couplings (AAC) across brain regions, we elucidated the ESI disparity between the eye-closed (EC) and eye-open (EO) conditions in rs-fMRI. In the visual network, the strength of the spectral interaction within 0.03-0.04 Hz was amplified in EC compared with that in EO condition, whereas a canonical connectivity analysis did not present differences between conditions. Collectively, leveraging from the instantaneous frequency of HHT, we firstly addressed the ESI technique to map the fMRI-based functional connectivity in a brand-new AAC perspective. The ESI possesses potential in elucidating the functional connectivity at specific frequency bins, thereby providing additional diagnostic merits for future clinical neuroscience.

Using non-invasive neuroimaging to enhance the care, well-being and experimental outcomes of laboratory non-human primates (monkeys)

Over the past 10-20 years, neuroscience witnessed an explosion in the use of non-invasive imaging methods, particularly magnetic resonance imaging (MRI), to study brain structure and function. Simultaneously, with access to MRI in many research institutions, MRI has become an indispensable tool for researchers and veterinarians to guide improvements in surgical procedures and implants and thus, experimental as well as clinical outcomes, given that access to MRI also allows for improved diagnosis and monitoring for brain disease. As part of the PRIMEatE Data Exchange, we gathered expert scientists, veterinarians, and clinicians who treat humans, to provide an overview of the use of non-invasive imaging tools, primarily MRI, to enhance experimental and welfare outcomes for laboratory non-human primates engaged in neuroscientific experiments. We aimed to provide guidance for other researchers, scientists and veterinarians in the use of this powerful imaging technology as well as to foster a larger conversation and community of scientists and veterinarians with a shared goal of improving the well-being and experimental outcomes for laboratory animals.

Oscillations as a window into neuronal mechanisms underlying dorsal anterior cingulate cortex function

The function of dorsal Anterior Cingulate Cortex (dACC) remains poorly understood. While many methods, spanning bottom-up and top-down approaches, have been deployed, the view they offer is often conflicting. Integrating bottom-up and top-down approaches requires an intermediary with sufficient explanatory power, theoretical development, and empirical support. Oscillations in the local field potential (LFP) provide such a link. LFP oscillations arise from empirically well-characterized neuronal circuit motifs. Synchronizing the firing of individual units has appealing properties to bind disparate brain regions and propagate information, including gating, routing, and coding. Moreover, the LFP, rather than single unit activity, more closely relates to macro-scale recordings, such as the electroencephalogram and functional magnetic resonance imaging. Thus, LFP oscillations are a critical link that allow for the inference of neuronal micro-circuitry underlying macroscopic brain recordings.

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

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

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

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

Functional connectivity with cortical depth assessed by resting state fMRI of subregions of S1 in squirrel monkeys

Whereas resting state blood oxygenation-level dependent (BOLD) functional MRI has been widely used to assess functional connectivity between cortical regions, the laminar specificity of such measures is poorly understood. This study aims to determine: (a) whether the resting state functional connectivity (rsFC) between two functionally related cortical regions varies with cortical depth, (b) the relationship between layer-resolved tactile stimulus-evoked activation pattern and interlayer rsFC pattern between two functionally distinct but related somatosensory areas 3b and 1, and (c) the effects of spatial resolution on rsFC measures. We examined the interlayer rsFC between areas 3b and 1 of squirrel monkeys under anesthesia using tactile stimulus-driven and resting state BOLD acquisitions at submillimeter resolution. Consistent with previous observations in the areas 3b and 1, we detected robust stimulus-evoked BOLD activations with foci were confined mainly to the upper layers (centered at 21% of the cortical depth). By carefully placing seeds in upper, middle, and lower layers of areas 3b and 1, we observed strong rsFC between upper and middle layers of these two areas. The layer-resolved activation patterns in areas 3b and 1 agree with their interlayer rsFC patterns, and are consistent with the known anatomical connections between layers. In summary, using BOLD rsFC pattern, we identified an interlayer interareal microcircuit that shows strong intrinsic functional connections between upper and middle layer areas 3b and 1. RsFC can be used as a robust invasive tool to probe interlayer corticocortical microcircuits.

Resting-state white matter-cortical connectivity in non-human primate brain

Numerous studies have used functional magnetic resonance imaging (fMRI) to characterize functional connectivity between cortical regions by analyzing correlations in blood oxygenation level dependent (BOLD) signals in a resting state. However, to date, there have been only a handful of studies reporting resting state BOLD signals in white matter. Nonetheless, a growing number of reports has emerged in recent years suggesting white matter BOLD signals can be reliably detected, though their biophysical origins remain unclear. Moreover, recent studies have identified robust correlations in a resting state between signals from cortex and specific white matter tracts. In order to further validate and interpret these findings, we studied a non-human primate model to investigate resting-state connectivity patterns between parcellated cortical volumes and specific white matter bundles. Our results show that resting-state connectivity patterns between white and gray matter structures are not randomly distributed but share notable similarities with diffusion- and histology-derived anatomic connectivities. This suggests that resting-state BOLD correlations between white matter fiber tracts and the gray matter regions to which they connect are directly related to the anatomic arrangement and density of WM fibers. We also measured how different levels of baseline neural activity, induced by varying levels of anesthesia, modulate these patterns. As anesthesia levels were raised, we observed weakened correlation coefficients between specific white matter tracts and gray matter regions while key features of the connectivity pattern remained similar. Overall, results from this study provide further evidence that neural activity is detectable by BOLD fMRI in both gray and white matter throughout the resting brain. The combined use of gray and white matter functional connectivity could also offer refined full-scale functional parcellation of the entire brain to characterize its functional architecture.

Dissociated resting-state functional networks between the dream recall frequency and REM sleep percentage

Rapid eye movement (REM) sleep has been frequently associated with dreaming. However, mounting evidence obtained from behavioral, pharmacological, and brain imaging studies suggests that REM sleep is not indicative of the dream report and may originate from diverse neural substrates in brain functionality. The aim of the current study was to investigate the functional systems associated with inter-individual differences in dream recall and REM sleep through assessments of the resting-state functional connectivity. We collected resting-state functional magnetic resonance imaging (fMRI) data for functional connectivity evaluations from 43 healthy adult volunteers (23 men) before and after sleep. For assessment of the dream recall frequency, a 2-week sleep diary was maintained by all volunteers. In addition, whole-night polysomnography was performed for measuring the REM sleep percentage. Voxel-wise correlation analyses of 12 functional connectivity networks of interest with the dream recall frequency and REM sleep percentage were conducted using general linear model analysis. Both the dream recall frequency and REM sleep percentage showed negative associations with multiple brain functional networks. However, the dream recall frequency was mainly related to functional connectivity within the lateral visual network and thalamus, whereas the REM sleep percentage was mainly associated with connectivity within the frontoparietal networks and cerebellum. In addition, the dream recall frequency showed stronger coupling with the lateral visual network connectivity at night, whereas the coupling between the REM sleep percentage and cerebellum was higher in the morning. This indicated a significant time of day effect. Our results provide neuroimaging evidence that the functional system associated with the dream recall frequency is different from that associated with the REM sleep percentage.

The VALiDATe29 MRI Based Multi-Channel Atlas of the Squirrel Monkey Brain

We describe the development of the first digital atlas of the normal squirrel monkey brain and present the resulting product, VALiDATe29. The VALiDATe29 atlas is based on multiple types of magnetic resonance imaging (MRI) contrast acquired on 29 squirrel monkeys, and is created using unbiased, nonlinear registration techniques, resulting in a population-averaged stereotaxic coordinate system. The atlas consists of multiple anatomical templates (proton density, T1, and T2* weighted), diffusion MRI templates (fractional anisotropy and mean diffusivity), and ex vivo templates (fractional anisotropy and a structural MRI). In addition, the templates are combined with histologically defined cortical labels, and diffusion tractography defined white matter labels. The combination of intensity templates and image segmentations make this atlas suitable for the fundamental atlas applications of spatial normalization and label propagation. Together, this atlas facilitates 3D anatomical localization and region of interest delineation, and enables comparisons of experimental data across different subjects or across different experimental conditions. This article describes the atlas creation and its contents, and demonstrates the use of the VALiDATe29 atlas in typical applications. The atlas is freely available to the scientific community.

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.

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.

Biophysical and neural basis of resting state functional connectivity: Evidence from non-human primates

Functional MRI (fMRI) has evolved from simple observations of regional changes in MRI signals caused by cortical activity induced by a task or stimulus, to task-free acquisitions of images in a resting state. Such resting state signals contain low frequency fluctuations which may be correlated between voxels, and strongly correlated regions are deemed to reflect functional connectivity within synchronized circuits. Resting state functional connectivity (rsFC) measures have been widely adopted by the neuroscience community, and are being used and interpreted as indicators of intrinsic neural circuits and their functional states in a broad range of applications, both basic and clinical. However, there has been relatively little work reported that validates whether inter-regional correlations in resting state fluctuations of fMRI (rsfMRI) signals actually measure functional connectivity between brain regions, or to establish how MRI data correlate with other metrics of functional connectivity. In this mini-review, we summarize recent studies of rsFC within mesoscopic scale cortical networks (100μm-10mm) within a well defined functional region of primary somatosensory cortex (S1), as well as spinal cord and brain white matter in non-human primates, in which we have measured spatial patterns of resting state correlations and validated their interpretation with electrophysiological signals and anatomic connections. Moreover, we emphasize that low frequency correlations are a general feature of neural systems, as evidenced by their presence in the spinal cord as well as white matter. These studies demonstrate the valuable role of high field MRI and invasive measurements in an animal model to inform the interpretation of human imaging studies.

Effects of isoflurane anesthesia on resting-state fMRI signals and functional connectivity within primary somatosensory cortex of monkeys

INTRODUCTION

Correlated low-frequency fluctuations of resting-state functional magnetic resonance imaging (rsfMRI) signals have been widely used for inferring intrinsic brain functional connectivity (FC). In animal studies, accurate estimate of anesthetic effects on rsfMRI signals is demanded for reliable interpretations of FC changes. We have previously shown that inter-regional FC can reliably delineate local millimeter-scale circuits within digit representations of primary somatosensory cortex (S1) subregions (areas 3a, 3b, and 1) in monkeys under isoflurane anesthesia. The goals of this study are to determine (1) the general effects of isoflurane on rsfMRI signals in the S1 circuit and (2) whether the effects are functional- and regional- dependent, by quantifying the relationships between isoflurane levels, power and inter-regional correlation coefficients in digit and face regions of distinct S1 subregions.

METHODS

Functional MRI data were collected from male adult squirrel monkeys at three different isoflurane levels (1.25%, 0.875%, and 0.5%). All scans were acquired on a 9.4T magnet with a 3-cm-diameter surface transmit-receive coil centered over the S1 cortex. Power and seed-based inter-regional functional connectivity analyses were subsequently performed.

RESULTS

As anesthesia level increased, we observed (1) diminishing amplitudes of signal fluctuations, (2) reduced power of fluctuations in the low-frequency band used for connectivity measurements, (3) decreased inter-voxel connectivity around seed regions, and (4) weakened inter-regional FC across all pairs of regions of interest (digit-to-digit). The low-frequency power measures derived from rsfMRI signals from control muscle regions, however, did not exhibit any isoflurane level-related changes. Within the isoflurane dosage range we tested, the inter-regional functional connectivity differences were still detectable, and the effects of isoflurane did not differ across region-of-interest (ROI) pairs.

CONCLUSION

Our data demonstrate that isoflurane induced similar dose-dependent suppressive effects on the power of rsfMRI signals and local fine-scale FC across functionally related but distinct S1 subregions.