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

Resting state networks of awake adolescent and adult squirrel monkeys using ultra-high field (9.4T) functional magnetic resonance imaging

Resting state networks (RSNs) are increasingly forwarded as candidate biomarkers for neuropsychiatric disorders. Such biomarkers may provide objective measures for evaluating novel therapeutic interventions in nonhuman primates often used in translational neuroimaging research. This study aimed to characterize the RSNs of awake squirrel monkeys and compare the characteristics of those networks in adolescent and adult subjects. Twenty-seven squirrel monkeys (n=12 adolescents [6 male/6 female] ∼2.5 years and n=15 adults [7 male/8 female] ∼9.5 years) were gradually acclimated to awake scanning procedures; whole-brain fMRI images were acquired with a 9.4 Tesla scanner. Group level independent component (ICA) analysis (30 ICs) with dual regression was used to detect and compare RSNs. Twenty ICs corresponding to physiologically meaningful networks representing a range of neural functions, including motor, sensory, reward, and cognitive processes were identified in both adolescent and adult monkeys. The reproducibility of these RSNs was evaluated across several ICA model orders. Adults showed a trend for greater connectivity compared to adolescent subjects in two of the networks of interest: (1) in the right occipital region with the OFC network and (2) in the left temporal cortex, bilateral occipital cortex, and cerebellum with the posterior cingulate network. However, when age was entered into the above model, this trend for significance was lost. These results demonstrate that squirrel monkey RSNs are stable and consistent with RSNs previously identified in humans, rodents, and other nonhuman primate species. These data also identify several networks in adolescence that are conserved and others that may change into adulthood.Significance Statement Functional magnetic resonance imaging procedures have revealed important information about how the brain is modified by experimental manipulations, disease states, and aging throughout the lifespan. Preclinical neuroimaging, especially in nonhuman primates, has become a frequently used means to answer targeted questions related to brain resting-state functional connectivity. The present study characterized resting state networks (RSNs) in adult and adolescent squirrel monkeys; twenty RSNs corresponding to networks representing a range of neural functions were identified. The RSNs identified here can be utilized in future studies examining the effects of experimental manipulations on brain connectivity in squirrel monkeys. These data also may be useful for comparative analysis with other primate species to provide an evolutionary perspective for understanding brain function and organization.

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.

Estimating and mitigating the effects of systemic low frequency oscillations (sLFO) on resting state networks in awake non-human primates using time lag dependent methodology

AIM

Resting-state fMRI (rs-fMRI) is often used to infer regional brain interactions from the degree of temporal correlation between spontaneous low-frequency fluctuations, thought to reflect local changes in the BOLD signal due to neuronal activity. One complication in the analysis and interpretation of rs-fMRI data is the existence of non-neuronal low frequency physiological noise (systemic low frequency oscillations; sLFOs) which occurs within the same low frequency band as the signal used to compute functional connectivity. Here, we demonstrate the use of a time lag mapping technique to estimate and mitigate the effects of the sLFO signal on resting state functional connectivity of awake squirrel monkeys.

METHODS

Twelve squirrel monkeys (6 male/6 female) were acclimated to awake scanning procedures; whole-brain fMRI images were acquired with a 9.4 Tesla scanner. Rs-fMRI data was preprocessed using an in-house pipeline and sLFOs were detected using a seed regressor generated by averaging BOLD signal across all voxels in the brain, which was then refined recursively within a time window of -16-12 s. The refined regressor was then used to estimate the voxel-wise sLFOs; these regressors were subsequently included in the general linear model to remove these moving hemodynamic components from the rs-fMRI data using general linear model filtering. Group level independent component analysis (ICA) with dual regression was used to detect resting-state networks and compare networks before and after sLFO denoising.

RESULTS

Results show sLFOs constitute ~64% of the low frequency fMRI signal in squirrel monkey gray matter; they arrive earlier in regions in proximity to the middle cerebral arteries (e.g., somatosensory cortex) and later in regions close to draining vessels (e.g., cerebellum). Dual regression results showed that the physiological noise was significantly reduced after removing sLFOs and the extent of reduction was determined by the brain region contained in the resting-state network.

CONCLUSION

These results highlight the need to estimate and remove sLFOs from fMRI data before further analysis.

Resting state networks of awake adolescent and adult squirrel monkeys using ultra-high field (9.4T) functional magnetic resonance imaging

Resting state networks (RSNs) are increasingly forwarded as candidate biomarkers for neuropsychiatric disorders. Such biomarkers may provide objective measures for evaluating novel therapeutic interventions in nonhuman primates often used in translational neuroimaging research. This study aimed to characterize the RSNs of awake squirrel monkeys and compare the characteristics of those networks in adolescent and adult subjects. Twenty-seven squirrel monkeys ( n =12 adolescents [6 male/6 female] ∼2.5 years and n =15 adults [7 male/8 female] ∼9.5 years) were gradually acclimated to awake scanning procedures; whole-brain fMRI images were acquired with a 9.4 Tesla scanner. Group level independent component (IC) analysis (30 ICs) with dual regression was used to detect and compare RSNs. Twenty ICs corresponding to physiologically meaningful networks representing a range of neural functions, including motor, sensory, reward (e.g., basal ganglia), and cognitive processes were identified in both adolescent and adult monkeys. Significant age-related differences between the adult and adolescent subjects (adult > adolescent) were found in two networks of interest: (1) the right upper occipital region with an OFC IC and (2) the left temporal cortex, bilateral visual areas, and cerebellum with the cingulate IC. These results demonstrate that squirrel monkey RSNs are stable and consistent with RSNs previously identified in humans, rodents, and other nonhuman primate species. These data also identify several networks in adolescence that are conserved and others that may change into adulthood.

Significance Statement

Functional magnetic resonance imaging procedures have revealed important information about how the brain is modified by experimental manipulations, disease states, and aging throughout the lifespan. Preclinical neuroimaging, especially in nonhuman primates, has become a frequently used means to answer targeted questions related to brain resting-state functional connectivity. The present study characterized resting state networks (RSNs) in adult and adolescent squirrel monkeys; twenty RSNs corresponding to networks representing a range of neural functions were identified. The RSNs identified here can be utilized in future studies examining the effects of experimental manipulations on brain connectivity in squirrel monkeys. These data also may be useful for comparative analysis with other primate species to provide an evolutionary perspective for understanding brain function and organization.

In vivo symmetric multi-contrast MRI brain templates and atlas for spontaneously hypertensive rats

Spontaneously hypertensive rats (SHRs) are a valuable animal model of essential hypertension. The increasing use of SHRs in neuroimaging has generated an urgent demand for a template set that provides a reference for advanced data analysis. Structural T2-weighted magnetic resonance imaging (MRI), diffusion tensor imaging (DTI) and functional MRI scans that were used to build the template set were obtained from 8 SHRs longitudinally scanned in vivo at 10, 24 and 52 weeks of age. These symmetric multi-contrast templates were constructed by iterative registration and averaging. The cortical atlas was derived from the Tohoku atlas, and the subcortical regions were manually delineated based on the templates. A set of SHR brain images named the Hebei Medical University rat brain template set (HRT) comprised 3D symmetric T2WI, raw T2-weighted signal with no added diffusion weighting (B0), fractional anisotropy (FA), mean diffusivity (MD) and blood oxygen level-dependent (BOLD) templates; tissue probability maps (TPMs) of gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF); and a whole-brain atlas with 163 labels. We quantitatively validated the template and characterized the longitudinal changes in brain morphology in different brain tissues as SHRs aged. To our knowledge, the HRT is the first MRI template set for SHRs. We believe that the HRT can serve as a beneficial tool for precise analysis of the SHR brain using structural and functional MRI, which can promote neuroimaging studies on essential hypertension.

Cortical Morphology and White Matter Tractography of Three Phylogenetically Distant Primates: Evidence for a Simian Elaboration

Comparative neuroimaging has been used to identify changes in white matter architecture across primate species phylogenetically close to humans, but few have compared the phylogenetically distant species. Here, we acquired postmortem diffusion imaging data from ring-tailed lemurs (Lemur catta), black-capped squirrel monkeys (Saimiri boliviensis), and rhesus macaques (Macaca mulatta). We were able to establish templates and surfaces allowing us to investigate sulcal, cortical, and white matter anatomy. The results demonstrate an expansion of the frontal projections of the superior longitudinal fasciculus complex in squirrel monkeys and rhesus macaques compared to ring-tailed lemurs, which correlates with sulcal anatomy and the lemur’s smaller prefrontal granular cortex. The connectivity of the ventral pathway in the parietal region is also comparatively reduced in ring-tailed lemurs, with the posterior projections of the inferior longitudinal fasciculus not extending toward parietal cortical areas as in the other species. In the squirrel monkeys we note a very specific occipito-parietal anatomy that is apparent in their surface anatomy and the expansion of the posterior projections of the optical radiation. Our study supports the hypothesis that the connectivity of the prefrontal-parietal regions became relatively elaborated in the simian lineage after divergence from the prosimian lineage.

The squirrel monkey model in clinical neuroscience

Clinical neuroscience research relying on animal models brought valuable translational insights into the function and pathologies of the human brain. The anatomical, physiological, and behavioural similarities between humans and mammals have prompted researchers to study cerebral mechanisms at different levels to develop and test new treatments. The vast majority of biomedical research uses rodent models, which are easily manipulable and have a broadly resembling organisation to the human nervous system but cannot satisfactorily mimic some disorders. For these disorders, macaque monkeys have been used as they have a more comparable central nervous system. Still, this research has been hampered by limitations, including high costs and reduced samples. This review argues that a squirrel monkey model might bridge the gap by complementing translational research from rodents, macaque, and humans. With the advent of promising new methods such as ultrasound imaging, tool miniaturisation, and a shift towards open science, the squirrel monkey model represents a window of opportunity that will potentially fuel new translational discoveries in the diagnosis and treatment of brain pathologies.

Long-term effects of intermittent early life stress on primate prefrontal-subcortical functional connectivity

Correlational studies of humans suggest that exposure to early life stress has long-term effects on neural circuits involved in vulnerability and resilience to mental health disorders. Stress-related mental health disorders are more prevalent in women than in men. Here, female squirrel monkeys are randomized to intermittently stressful (IS) social separations or a non-separated (NS) control condition conducted from 17 to 27 weeks of age. Nine years later in mid-life adulthood, resting-state functional magnetic resonance imaging was employed to parcellate prefrontal cortex (PFC). Resulting subdivisions were then used to characterize functional connectivity within PFC, and between PFC subdivisions and subcortical regions that are known to be altered by stress. Extensive hyper-connectivity of medial and orbitofrontal PFC with amygdala, hippocampus, and striatum was observed in IS compared to NS monkeys. Functional hyper-connectivity in IS monkeys was associated with previously reported indications of diminished anxiety-like behavior induced by prepubertal stress. Hyper-connectivity of PFC with amygdala and with hippocampus was also associated with increased ventral striatal dopamine D2 and/or D3 receptor (DRD2/3) availability assessed with positron emission tomography (PET) of [¹¹C]raclopride binding in adulthood. Ventral striatal DRD2/3 availability has been linked to cognitive control, which plays a key role in stress coping as an aspect of emotion regulation. These findings provide causal support for enduring neurobiological effects of early life stress and suggest novel targets for new treatments of stress-related mental health disorders.

Evolution of Brain Connections: Integrating Diffusion MR Tractography With Gene Expression Highlights Increased Corticocortical Projections in Primates

Diffusion MR tractography permits investigating the 3D structure of cortical pathways as interwoven paths across the entire brain. We use high-resolution scans from diffusion spectrum imaging and high angular resolution diffusion imaging to investigate the evolution of cortical pathways within the euarchontoglire (i.e., primates, rodents) lineage. More specifically, we compare cortical fiber pathways between macaques (Macaca mulatta), marmosets (Callithrix jachus), and rodents (mice, Mus musculus). We integrate these observations with comparative analyses of Neurofilament heavy polypeptide (NEFH) expression across the cortex of mice and primates. We chose these species because their phylogenetic position serves to trace the early evolutionary history of the human brain. Our comparative analysis from diffusion MR tractography, cortical white matter scaling, and NEFH expression demonstrates that the examined primates deviate from mice in possessing increased long-range cross-cortical projections, many of which course across the anterior to posterior axis of the cortex. Our study shows that integrating gene expression data with diffusion MR data is an effective approach in identifying variation in connectivity patterns between species. The expansion of corticocortical pathways and increased anterior to posterior cortical integration can be traced back to an extension of neurogenetic schedules during development in primates.

High Angular Resolution Diffusion MRI Reveals Conserved and Deviant Programs in the Paths that Guide Human Cortical Circuitry

Diffusion magnetic resonance (MR) tractography represents a novel opportunity to investigate conserved and deviant developmental programs between humans and other species such as mice. To that end, we acquired high angular resolution diffusion MR scans of mice [embryonic day (E) 10.5 to postnatal week 4] and human brains [gestational week (GW) 17-30] at successive stages of fetal development to investigate potential evolutionary changes in radial organization and emerging pathways between humans and mice. We compare radial glial development as well as commissural development (e.g., corpus callosum), primarily because our findings can be integrated with previous work. We also compare corpus callosal growth trajectories across primates (i.e., humans and rhesus macaques) and rodents (i.e., mice). One major finding is that the developing cortex of humans is predominated by pathways likely associated with a radial glial organization at GW 17-20, which is not as evident in age-matched mice (E 16.5, 17.5). Another finding is that, early in development, the corpus callosum follows a similar developmental timetable in primates (i.e., macaques and humans) as in mice. However, the corpus callosum grows for an extended period of time in primates compared with rodents. Taken together, these findings highlight deviant developmental programs underlying the emergence of cortical pathways in the human brain.

Large-scale comparative neuroimaging: Where are we and what do we need?

Neuroimaging has a lot to offer comparative neuroscience. Although invasive "gold standard" techniques have a better spatial resolution, neuroimaging allows fast, whole-brain, repeatable, and multi-modal measurements of structure and function in living animals and post-mortem tissue. In the past years, comparative neuroimaging has increased in popularity. However, we argue that its most significant potential lies in its ability to collect large-scale datasets of many species to investigate principles of variability in brain organisation across whole orders of species-an ambition that is presently unfulfilled but achievable. We briefly review the current state of the field and explore what the current obstacles to such an approach are. We propose some calls to action.

A Web-Based Atlas Combining MRI and Histology of the Squirrel Monkey Brain

The squirrel monkey (Saimiri sciureus) is a commonly-used surrogate for humans in biomedical research. In the neuroimaging community, MRI and histological atlases serve as valuable resources for anatomical, physiological, and functional studies of the brain; however, no digital MRI/histology atlas is currently available for the squirrel monkey. This paper describes the construction of a web-based multi-modal atlas of the squirrel monkey brain. The MRI-derived information includes anatomical MRI contrast (i.e., T2-weighted and proton-density-weighted) and diffusion MRI metrics (i.e., fractional anisotropy and mean diffusivity) from data acquired both in vivo and ex vivo on a 9.4 Tesla scanner. The histological images include Nissl and myelin stains, co-registered to the corresponding MRI, allowing identification of cyto- and myelo-architecture. In addition, a bidirectional neuronal tracer, biotinylated dextran amine (BDA) was injected into the primary motor cortex, enabling highly specific identification of regions connected to the injection location. The atlas integrates the results of common image analysis methods including diffusion tensor imaging glyphs, labels of 57 white-matter tracts identified using DTI-tractography, and 18 cortical regions of interest identified from Nissl-revealed cyto-architecture. All data are presented in a common space, and all image types are accessible through a web-based atlas viewer, which allows visualization and interaction of user-selectable contrasts and varying resolutions. By providing an easy to use reference system of anatomical information, our web-accessible multi-contrast atlas forms a rich and convenient resource for comparisons of brain findings across subjects or modalities. The atlas is called the Combined Histology-MRI Integrated Atlas of the Squirrel Monkey (CHIASM). All images are accessible through our web-based viewer ( https://chiasm.vuse.vanderbilt.edu /), and data are available for download at ( https://www.nitrc.org/projects/smatlas/ ).

Histologically derived fiber response functions for diffusion MRI vary across white matter fibers-An ex vivo validation study in the squirrel monkey brain

Understanding the relationship between the diffusion-weighted MRI signal and the arrangement of white matter fibers is fundamental for accurate voxel-wise reconstruction of the fiber orientation distribution (FOD) and subsequent fiber tractography. Spherical deconvolution reconstruction techniques model the diffusion signal as the convolution of the FOD with a response function that represents the signal profile of a single fiber orientation. Thus, given the signal and a fiber response function, the FOD can be estimated in every imaging voxel by deconvolution. However, the selection of the appropriate response function remains relatively under-studied, and requires further validation. In this work, using 3D histologically defined FODs and the corresponding diffusion signal from three ex vivo squirrel monkey brains, we derive the ground truth response functions. We find that the histologically derived response functions differ from those conventionally used. Next, we find that response functions statistically vary across brain regions, which suggests that the practice of using the same kernel throughout the brain is not optimal. We show that different kernels lead to different FOD reconstructions, which in turn can lead to different tractography results depending on algorithmic parameters, with large variations in the accuracy of resulting reconstructions. Together, these results suggest there is room for improvement in estimating and understanding the relationship between the diffusion signal and the underlying FOD.

Discrimination of the hierarchical structure of cortical layers in 2-photon microscopy data by combined unsupervised and supervised machine learning

The laminar organization of the cerebral cortex is a fundamental characteristic of the brain, with essential implications for cortical function. Due to the rapidly growing amount of high-resolution brain imaging data, a great demand arises for automated and flexible methods for discriminating the laminar texture of the cortex. Here, we propose a combined approach of unsupervised and supervised machine learning to discriminate the hierarchical cortical laminar organization in high-resolution 2-photon microscopic neural image data of mouse brain without observer bias, that is, without the prerequisite of manually labeled training data. For local cortical foci, we modify an unsupervised clustering approach to identify and represent the laminar cortical structure. Subsequently, supervised machine learning is applied to transfer the resulting layer labels across different locations and image data, to ensure the existence of a consistent layer label system. By using neurobiologically meaningful features, the discrimination results are shown to be consistent with the layer classification of the classical Brodmann scheme, and provide additional insight into the structure of the cerebral cortex and its hierarchical organization. Thus, our work paves a new way for studying the anatomical organization of the cerebral cortex, and potentially its functional organization.

Challenges in diffusion MRI tractography - Lessons learned from international benchmark competitions

Diffusion MRI (dMRI) fiber tractography has become a pillar of the neuroimaging community due to its ability to noninvasively map the structural connectivity of the brain. Despite widespread use in clinical and research domains, these methods suffer from several potential drawbacks or limitations. Thus, validating the accuracy and reproducibility of techniques is critical for sound scientific conclusions and effective clinical outcomes. Towards this end, a number of international benchmark competitions, or "challenges", has been organized by the diffusion MRI community in order to investigate the reliability of the tractography process by providing a platform to compare algorithms and results in a fair manner, and evaluate common and emerging algorithms in an effort to advance the state of the field. In this paper, we summarize the lessons from a decade of challenges in tractography, and give perspective on the past, present, and future "challenges" that the field of diffusion tractography faces.

A fiber coherence index for quality control of B-table orientation in diffusion MRI scans

PURPOSE

The diffusion MRI "b-vector" table describing the diffusion sensitization direction can be flipped and permuted in dimension due to different orientation conventions used in scanners and incorrect or improperly utilized file formats. This can lead to incorrect fiber orientation estimates and subsequent tractography failure. Here, we present an automated quality control procedure to detect when the b-table is flipped and/or permuted incorrectly.

METHODS

We define a "fiber coherence index" to describe how well fibers are connected to each other, and use it to automatically detect the correct configuration of b-vectors. We examined the performance on 3981 research subject scans (Baltimore Longitudinal Study of Aging), 1065 normal subject scans of high image quality (Human Connectome Project), and 202 patient scans (Vanderbilt University Medical Center), as well as 9 in-vivo and 9 ex-vivo animal data.

RESULTS

The coherence index resulted in a 99.9% (3979/3981) and 100% (1065/1065) success rate in normal subject scans, 98% (198/202) in patient scans, and 100% (18/18) in both in-vivo and ex-vivo animal data in detecting the correct gradient table in datasets without severe image artifacts. The four failing cases (4/202) in patient scans, and two failures in healthy subject scans (2/3981), all showed prominent motion or signal dropout artifacts.

CONCLUSIONS

The fiber coherence measure can be used as an automatic quality assurance check in any diffusion analysis pipeline. Additionally, the success of this fiber coherence measure suggests potential broader applications, including evaluating data quality, or even providing diagnostic value as a biomarker of white matter integrity.

Limits to anatomical accuracy of diffusion tractography using modern approaches

Diffusion MRI fiber tractography is widely used to probe the structural connectivity of the brain, with a range of applications in both clinical and basic neuroscience. Despite widespread use, tractography has well-known pitfalls that limits the anatomical accuracy of this technique. Numerous modern methods have been developed to address these shortcomings through advances in acquisition, modeling, and computation. To test whether these advances improve tractography accuracy, we organized the 3-D Validation of Tractography with Experimental MRI (3D-VoTEM) challenge at the ISBI 2018 conference. We made available three unique independent tractography validation datasets - a physical phantom and two ex vivo brain specimens - resulting in 176 distinct submissions from 9 research groups. By comparing results over a wide range of fiber complexities and algorithmic strategies, this challenge provides a more comprehensive assessment of tractography's inherent limitations than has been reported previously. The central results were consistent across all sub-challenges in that, despite advances in tractography methods, the anatomical accuracy of tractography has not dramatically improved in recent years. Taken together, our results independently confirm findings from decades of tractography validation studies, demonstrate inherent limitations in reconstructing white matter pathways using diffusion MRI data alone, and highlight the need for alternative or combinatorial strategies to accurately map the fiber pathways of the brain.

Anatomical accuracy of standard-practice tractography algorithms in the motor system - A histological validation in the squirrel monkey brain

For two decades diffusion fiber tractography has been used to probe both the spatial extent of white matter pathways and the region to region connectivity of the brain. In both cases, anatomical accuracy of tractography is critical for sound scientific conclusions. Here we assess and validate the algorithms and tractography implementations that have been most widely used - often because of ease of use, algorithm simplicity, or availability offered in open source software. Comparing forty tractography results to a ground truth defined by histological tracers in the primary motor cortex on the same squirrel monkey brains, we assess tract fidelity on the scale of voxels as well as over larger spatial domains or regional connectivity. No algorithms are successful in all metrics, and, in fact, some implementations fail to reconstruct large portions of pathways or identify major points of connectivity. The accuracy is most dependent on reconstruction method and tracking algorithm, as well as the seed region and how this region is utilized. We also note a tremendous variability in the results, even though the same MR images act as inputs to all algorithms. In addition, anatomical accuracy is significantly decreased at increased distances from the seed. An analysis of the spatial errors in tractography reveals that many techniques have trouble properly leaving the gray matter, and many only reveal connectivity to adjacent regions of interest. These results show that the most commonly implemented algorithms have several shortcomings and limitations, and choices in implementations lead to very different results. This study should provide guidance for algorithm choices based on study requirements for sensitivity, specificity, or the need to identify particular connections, and should serve as a heuristic for future developments in tractography.

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.