Macaques Research Articles

Invited Session III: Diversity in chromatic processing across the animal kingdom: Comparative connectomics reveals noncanonical wiring for color vision in human foveal retina.

The Old World macaque monkey and New World common marmoset provide fundamental models for human visual processing, yet the human ancestral lineage diverged from these monkey lineages over 25 Mya. We therefore asked whether fine-scale synaptic wiring in the nervous system is preserved across these three primate families, despite long periods of independent evolution. We applied connectomic electron microscopy to the specialized foveal retina where circuits for highest acuity and color vision reside. Synaptic motifs arising from the cone photoreceptor type sensitive to short (S) wavelengths and associated with "blue-yellow" (S-ON and S-OFF) color-coding circuitry were reconstructed. We found that distinctive circuitry arises from S cones for each of the three species. The S cones contacted neighboring L and M (long- and middle-wavelength sensitive) cones in humans, but such contacts were rare or absent in macaques and marmosets. We discovered a major S-OFF pathway in the human retina and established its absence in marmosets. Further, the S-ON and S-OFF chromatic pathways make excitatory-type synaptic contacts with L and M cone types in humans, but not in macaques or marmosets. Our results predict that early-stage chromatic signals are distinct in the human retina and imply that solving the human connectome at the nanoscale level of synaptic wiring will be critical for fully understanding the neural basis of human color vision.

Decoding motor plans using a closed-loop ultrasonic brain–machine interface

Brain–machine interfaces (BMIs) enable people living with chronic paralysis to control computers, robots and more with nothing but thought. Existing BMIs have trade-offs across invasiveness, performance, spatial coverage and spatiotemporal resolution. Functional ultrasound (fUS) neuroimaging is an emerging technology that balances these attributes and may complement existing BMI recording technologies. In this study, we use fUS to demonstrate a successful implementation of a closed-loop ultrasonic BMI. We streamed fUS data from the posterior parietal cortex of two rhesus macaque monkeys while they performed eye and hand movements. After training, the monkeys controlled up to eight movement directions using the BMI. We also developed a method for pretraining the BMI using data from previous sessions. This enabled immediate control on subsequent days, even those that occurred months apart, without requiring extensive recalibration. These findings establish the feasibility of ultrasonic BMIs, paving the way for a new class of less-invasive (epidural) interfaces that generalize across extended time periods and promise to restore function to people with neurological impairments.

Cross-species amplification from non-human primate DNAs in commercial human DNA assay kits

Species specificity of commercial human DNA quantification kits and short tandem repeat (STR) profiling kits was examined using primate DNA samples. These samples comprised 33 individuals from eight primate species, each with gender and kinship data, including human (Homo sapiens), chimpanzee (Pan troglodytes), gorilla (Gorilla gorilla), and orangutan (Pongo pygmaeus) of Hominidae family, and Japanese macaque (Macaca fuscata), long-tailed macaque (Macaca fascicularis), hamadryas baboon (Papio hamadryas), and savannah monkey (Chlorocebus sp.) of Cercopithecidae family. The findings revealed varying levels of cross-species amplifications in all non-human DNA samples that correlated with their evolutionary proximity to humans, both kit types. Moreover, cross-species amplification, including female DNA samples, was observed in a Y-chromosomal STR profiling kit. Additionally, species specificity differed among the commercial kits examined. The cross-species amplification data presented in this study offer valuable assistance in interpreting the results of individual human identification in forensic cases involving non-human primates.

Beat processing in newborn infants cannot be explained by statistical learning based on transition probabilities

Newborn infants have been shown to extract temporal regularities from sound sequences, both in the form of learning regular sequential properties, and extracting periodicity in the input, commonly referred to as a regular pulse or the ‘beat’. However, these two types of regularities are often indistinguishable in isochronous sequences, as both statistical learning and beat perception can be elicited by the regular alternation of accented and unaccented sounds. Here, we manipulated the isochrony of sound sequences in order to disentangle statistical learning from beat perception in sleeping newborn infants in an EEG experiment, as previously done in adults and macaque monkeys. We used a binary accented sequence that induces a beat when presented with isochronous timing, but not when presented with randomly jittered timing. We compared mismatch responses to infrequent deviants falling on either accented or unaccented (i.e., odd and even) positions. Results showed a clear difference between metrical positions in the isochronous sequence, but not in the equivalent jittered sequence. This suggests that beat processing is present in newborns. Despite previous evidence for statistical learning in newborns the effects of this ability were not detected in the jittered condition. These results show that statistical learning by itself does not fully explain beat processing in newborn infants.

A proposed structural connectivity matrices approach for the superior fronto-occipital fascicle in the Harvard-Oxford Atlas comparative framework following the Pandya comparative extrapolation principle.

A key set of connections necessary for the most complex brain functions are the long association cortico-cortical fiber tracts. These pathways have been described by the Dejerines and others using post mortem histological or brain dissection techniques. Given methodological limitations, these fiber connections have not been delineated completely in humans. Although the stem portions of fiber tracts have been identified in humans, their precise origins and terminations remain to be determined. By contrast, the origins and terminations as well as the stems of long cortico-cortical association fiber pathways in monkeys have been detailed in the macaque monkey brain using experimental tract tracing methods. Deepak Pandya made major contributions to the delineation of fiber tracts in the monkey brain. In the early 1990s, he compared his observations in monkeys with the original descriptions in humans by the Dejerines. With the advent of diffusion-weighted imaging, Dr. Pandya extended this line of investigation to the human brain with Dr. Nikos Makris. In this translational analysis of long association cortico-cortical fiber tracts, they applied a principle of extrapolation from monkey to human. In the present study, we addressed the reasoning and the complex methodology in translating brain structural connectivity from monkey to human in one cortico-cortical fiber tract, namely the superior fronto-occipital fascicle, which was delineated in both species by Dr. Pandya and colleagues. Furthermore, we represented this information in the form of connectional matrices in the context of the HOA2.0-ComPaRe framework, a homological monkey-to-human translational system used in neuroimaging.

Inserting a Neuropixels probe into awake monkey cortex: two probes, two methods

BackgroundNeuropixels probes have revolutionized neurophysiological studies in the rodent, but inserting these probes through the much thicker primate dura remains a challenge. New methodsHere we describe two methods we have developed for the insertion of two types of Neuropixels probes acutely into the awake macaque monkey cortex. For the fine rodent probe (Neuropixels 1.0, IMEC), which is unable to pierce native primate dura, we developed a dural-eyelet method to insert the probe repeatedly without breakage. For the thicker short NHP probe (Neuropixels NP1010), we developed an artificial dura system to insert the probe. Results and comparison with existing methodsWe have now conducted successful experiments in 3 animals across 7 recording chambers with the procedures described here and have achieved recordings with similar yields over several months in each case. ConclusionWe hope that our hardware, surgical preparation, methods for insertion and methods for removal of broken probe parts are of value to primate physiologists everywhere.

Opportunistic predation and predation-related events on long-tailed macaque and proboscis monkey in Kinabatangan, Sabah, Malaysia

We report two cases of predation and predation-related events on long-tailed macaques and proboscis monkeys living in sympatry in the lower Kianabatangan, Sabah, Malaysia.

Investigating the dietary niches of fossil Plio-Pleistocene European macaques: The case of Macaca majori Azzaroli, 1946 from Sardinia

The genus Macaca includes medium- to large-bodied monkeys and represents one of the most diverse primate genera, also having a very large geographic range. Nowadays, wild macaque populations are found in Asia and Africa, inhabiting a wide array of habitats. Fossil macaques were also present in Europe from the Late Miocene until the Late Pleistocene. Macaques are considered ecologically flexible monkeys that exhibit highly opportunistic dietary strategies, which may have been critical to their evolutionary success. Nevertheless, available ecological information regarding fossil European species is very sparse, limiting our knowledge of their evolutionary history in this geographic area. To further our understanding of fossil European macaque ecology, we investigated the dietary ecology of Macaca majori, an insular endemic species from Sardinia. In particular, we characterized the dental capabilities and potential dietary adaptations of M. majori through dental topographic and enamel thickness analyses of two M2s from the Early Pleistocene site of Capo Figari (1.8 Ma). We also assessed its diet through dental microwear texture analysis, while the microwear texture of M. majori was also compared with microwear textures from other European fossil macaques from mainland Europe. The dental topographic and enamel thickness analyses suggest that M. majori frequently consumes hard/mechanically challenging and/or abrasive foods. The results of the dental microwear analysis are consistent with this interpretation and further suggest that M. majori probably exhibited more durophagous dietary habits than mainland Plio-Pleistocene macaques. Overall, our results indicate that M. majori probably occupied a different dietary niche compared to its mainland fossil relatives, which suggests that they may have inhabited different paleoenvironments.

Cellular and laminar architecture: A short history and commentary.

The feedforward/feedback classification, as originally stated in relation to early visual areas in the macaque monkey, has had a significant influence on ideas of laminar interactions, area reciprocity, and cortical hierarchical organization. In some contrast with this macroscale "laminar connectomics," a more cellular approach to cortical connections, as briefly surveyed here, points to a still underappreciated heterogeneity of neuronal subtypes and complex microcircuitries. From the perspective of heterogeneities, the question of how brain regions interact and influence each other quickly leads to discussions about concurrent hierarchical and nonhierarchical cortical features, brain organization as a multiscale system forming nested groups and hierarchies, connectomes annotated by multiple biological attributes, and interleaved and overlapping scales of organization.

Local field potential responses associated with perceptual filling-in at blind spot in macaque primary visual cortex.

To understand the neural mechanisms of perceptual filling-in at the blind spot (BS), we analyzed neural activity in the region representing the visual field corresponding to the BS (BS region) in the primary visual cortex (V1) of the macaque monkey. We inserted a linear array electrode into the BS region or surrounding region and recorded the multiunit activities (MUAs) and local field potential (LFP). We examined the responses of MUAs and LFP to a large visual stimulus that entirely covered the BS (surface stimuli) while the monkey performed a visual fixation task in either the monocular condition without receiving direct retinal input or the binocular condition receiving retinal information. We observed clear MUA responses in the deep layers within the BS region under monocular conditions, confirming previous reports that V1 neurons in the BS region are activated when perceptual filling-in occurs. Current source density analysis using LFP showed that MUA responses were mainly observed in layer 5. Although LFP responses were generally stronger in the binocular condition than in the monocular condition, a notable exception was observed in the BS region. LFP responses in the low-beta band in the superficial layers were stronger in the monocular condition than in the binocular condition. These results suggest that low-beta activity in the superficial layer is related to the occurrence of perceptual filling-in in the BS. The origin of this activity is considered to be feedback signals from the extrastriate areas to the V1.NEW & NOTEWORTHY Two characteristic activities were induced in the blind spot (BS) region in response to the stimulus, causing perceptual filling-in: 1) beta-band LFP responses in the superficial layers and 2) neuronal responses in the deep layers, mainly in layer 5. These data suggest that the feedback signal from the extrastriate areas to the BS region in V1 is involved in perceptual filling-in.

Quantification of human mature frataxin protein expression in nonhuman primate hearts after gene therapy

Deficiency in human mature frataxin (hFXN-M) protein is responsible for the devastating neurodegenerative and cardiodegenerative disease of Friedreich’s ataxia (FRDA). It results primarily through epigenetic silencing of the FXN gene by GAA triplet repeats on intron 1 of both alleles. GAA repeat lengths are most commonly between 600 and 1200 but can reach 1700. A subset of approximately 3% of FRDA patients have GAA repeats on one allele and a mutation on the other. FRDA patients die most commonly in their 30s from heart disease. Therefore, increasing expression of heart hFXN-M using gene therapy offers a way to prevent early mortality in FRDA. We used rhesus macaque monkeys to test the pharmacology of an adeno-associated virus (AAV)hu68.CB7.hFXN therapy. The advantage of using non-human primates for hFXN-M gene therapy studies is that hFXN-M and monkey FXN-M (mFXN-M) are 98.5% identical, which limits potential immunologic side-effects. However, this presented a formidable bioanalytical challenge in quantification of proteins with almost identical sequences. This could be overcome by the development of a species-specific quantitative mass spectrometry-based method, which has revealed for the first time, robust transgene-specific human protein expression in monkey heart tissue. The dose response is non-linear resulting in a ten-fold increase in monkey heart hFXN-M protein expression with only a three-fold increase in dose of the vector.

Motor actions are spatially organized in motor and dorsal premotor cortex.

Frontal motor areas are central to controlling voluntary movements. In non-human primates, the motor areas contain independent, somatotopic, representations of the forelimb (i.e., motor maps). But are the neural codes for actions spatially organized within those forelimb representations? Addressing this question would provide insight into the poorly understood structure-function relationships of the cortical motor system. Here, we tackle the problem using high resolution optical imaging and motor mapping in motor (M1) and dorsal premotor (PMd) cortex. Two macaque monkeys performed an instructed reach-to-grasp task while cortical activity was recorded with intrinsic signal optical imaging (ISOI). The spatial extent of activity in M1 and PMd was then quantified in relation to the forelimb motor maps, which we obtained from the same hemisphere with intracortical microstimulation. ISOI showed that task-related activity was concentrated in patches that collectively overlapped <40% of the M1 and PMd forelimb representations. The spatial organization of the patches was consistent across task conditions despite small variations in forelimb use. Nevertheless, the largest condition differences in forelimb use were reflected in the magnitude of cortical activity. Distinct time course profiles from patches in arm zones and patches in hand zones suggest functional differences within the forelimb representations. The results collectively support an organizational framework wherein the forelimb representations contain subzones enriched with neurons tuned for specific actions. Thus, the often-overlooked spatial dimension of neural activity appears to be an important organizing feature of the neural code in frontal motor areas.

A competition framework for fixation-preference in strabismus.

Strabismic subjects often develop the ability to fixate on a target with either eye. Previous studies have shown that fixation-preference behavior varies systematically depending on spatial location of the target. We hypothesized that, when an eccentric target is presented, oculomotor fixation-preference in strabismus may be accounted for in a competitive decision framework wherein the brain must choose between two possible retinal errors to prepare a conjugate saccade that results in one of the eyes acquiring the eccentric target. We tested this framework by recording from visuo-motor neurons in the superior colliculus (SC) of two strabismic rhesus macaque monkeys as they performed a delayed saccade task under binocular viewing conditions. In one experiment, visual targets were presented at one of two locations corresponding to the neuronal receptive field location with respect to either the viewing or the deviated eye. Robust visual sensory responses were observed when targets were presented at either location indicating the presence of competing sensory signals for eye-choice. In a second experiment, a single visual target was placed at the neuronal receptive field location where the animal switched fixation on some trials and did not on other trials. At such target locations where either eye could acquire the target, both visual and build-up activity was greater in trials when the saccade encoded by the neuron "won." These findings provide evidence for the influence of visual suppression within SC sensory activity and support the possible utilization of a competition framework, one that has been previously described for when a binocularly aligned animal chooses from among multiple targets, to drive fixation-preference behavior in strabismus.

No Evidence for Cross-Modal fMRI Adaptation in Macaque Parieto-Premotor Mirror Neuron Regions.

To probe the presence of mirror neurons in the human brain, cross-modal fMRI adaptation has been suggested as a suitable technique. The rationale behind this suggestion is that this technique allows making more accurate inferences about neural response properties underlying fMRI voxel activations, beyond merely showing shared voxels that are active during both action observation and execution. However, the validity of using cross-modal fMRI adaptation to demonstrate the presence of mirror neurons in parietal and premotor brain regions has been questioned given the inconsistent and weak results obtained in human studies. A better understanding of cross-modal fMRI adaptation effects in the macaque brain is required as the rationale for using this approach is based on several assumptions related to macaque mirror neuron response properties that still need validation. Here, we conducted a cross-modal fMRI adaptation study in macaque monkeys, using the same action execution and action observation tasks that successfully yielded mirror neuron region cross-modal action decoding in a previous monkey MVPA study. We scanned two male rhesus monkeys while they first executed a sequence of either reach-and-grasp or reach-and-touch hand actions and then observed a video of a human actor performing these motor acts. Both whole-brain and region-of-interest analyses failed to demonstrate cross-modal fMRI adaptation effects in parietal and premotor mirror neuron regions. Our results, in line with previous findings in non-human primates, show that cross-modal motor-to-visual fMRI adaptation is not easily detected in monkey brain regions known to house mirror neurons. Thus, our results advocate caution in using cross-modal fMRI adaptation as a method to infer whether mirror neurons can be found in the primate brain.

Activity in primate visual cortex is minimally driven by spontaneous movements

Organisms process sensory information in the context of their own moving bodies, an idea referred to as embodiment. This idea is important for developmental neuroscience, robotics and systems neuroscience. The mechanisms supporting embodiment are unknown, but a manifestation could be the observation in mice of brain-wide neuromodulation, including in the primary visual cortex, driven by task-irrelevant spontaneous body movements. We tested this hypothesis in macaque monkeys (Macaca mulatta), a primate model for human vision, by simultaneously recording visual cortex activity and facial and body movements. We also sought a direct comparison using an analogous approach to those used in mouse studies. Here we found that activity in the primate visual cortex (V1, V2 and V3/V3A) was associated with the animals’ own movements, but this modulation was largely explained by the impact of the movements on the retinal image, that is, by changes in visual input. These results indicate that visual cortex in primates is minimally driven by spontaneous movements and may reflect species-specific sensorimotor strategies.

DNA methylation and transcriptome analysis reveal epigenomic differences among three macaque species.

Macaques (genus Macaca) are the most widely distributed non-human primates, and their evolutionary history, gene expression profiles, and genetic differences have been extensively studied. However, the DNA methylomes of macaque species are not available in public databases, which hampers understanding of epigenetic differences among macaque species. Epigenetic modifications can potentially affect development, physiology, behavior, and evolution. Here, we investigated the methylation patterns of the Tibetan macaque (M. thibetana; TM), Chinese rhesus macaque (M. mulatta lasiota; CR), and crab-eating macaque (M. fascicularis; CE) through whole-genome bisulfite sequencing from peripheral blood. We compared genome-wide methylation site information for the three species. We identified 12,128 (CR vs. CE), 59,165 (CR vs. TM), and 39,751 (CE vs. TM) differentially methylated regions (DMRs) in the three macaques. Furthermore, we obtained the differentially expressed genes (DEGs) among the three macaque species. The differences between CR and CE were smaller at both the methylome and transcriptome levels than compared with TM (CR vs. TM and CE vs. TM). We also found a change in the density of single nucleotide mutations in DMRs relative to their flanking regions, indicating a potential mechanism through which genomic alterations may modulate methylation landscapes, thereby influencing the transcriptome. Functional enrichment analyses showed the DMR-related genes were enriched in developmental processes and neurological functions, such as the growth hormone-related pathway, insulin secretion pathway, thyroid hormone synthesis pathway, morphine addiction, and GABAergic synapses. These differences may be associated with variations in physiology and habitat among the macaques. Our study provides one of the first genome-wide comparisons of genetic, gene expression, and epigenetic variations across different macaques. Our results should facilitate further research on comparative genomic and genetic differences in macaque species.

Phylogenomics Reveals High Levels of Incomplete Lineage Sorting at the Ancestral Nodes of the Macaque Radiation.

The genus Macaca includes 23 species assigned into 4 to 7 groups. It exhibits the largest geographic range and represents the most successful example of adaptive radiation of nonhuman primates. However, intrageneric phylogenetic relationships among species remain controversial and have not been resolved so far. In this study, we conducted a phylogenomic analysis on 16 newly generated and 8 published macaque genomes. We found strong evidence supporting the division of this genus into 7 species groups. Incomplete lineage sorting (ILS) was the primary factor contributing to the discordance observed among gene trees; however, we also found evidence of hybridization events, specifically between the ancestral arctoides/sinica and silenus/nigra lineages that resulted in the hybrid formation of the fascicularis/mulatta group. Combined with fossil data, our phylogenomic data were used to establish a scenario for macaque radiation. These findings provide insights into ILS and potential ancient introgression events that were involved in the radiation of macaques, which will lead to a better understanding of the rapid speciation occurring in nonhuman primates.

Visual feature tuning properties of stimulus-driven saccadic inhibition in macaque monkeys.

Saccadic inhibition refers to a short-latency transient cessation of saccade generation after visual sensory transients. This oculomotor phenomenon occurs with a latency that is consistent with a rapid influence of sensory responses, such as stimulus-induced visual bursts, on oculomotor control circuitry. However, the neural mechanisms underlying saccadic inhibition are not well understood. Here, we exploited the fact that macaque monkeys experience robust saccadic inhibition to test the hypothesis that inhibition time and strength exhibit systematic visual feature tuning properties to a multitude of visual feature dimensions commonly used in vision science. We measured saccades in three monkeys actively controlling their gaze on a target, and we presented visual onset events at random times. Across seven experiments, the visual onsets tested size, spatial frequency, contrast, orientation, motion direction, and motion speed dependencies of saccadic inhibition. We also investigated how inhibition might depend on the behavioral relevance of the appearing stimuli. We found that saccadic inhibition starts earlier, and is stronger, for large stimuli of low spatial frequencies and high contrasts. Moreover, saccadic inhibition timing depends on motion direction and orientation, with earlier inhibition systematically occurring for horizontally drifting vertical gratings. On the other hand, saccadic inhibition is stronger for faster motions and when the appearing stimuli are subsequently foveated. Besides documenting a range of feature tuning dimensions of saccadic inhibition to the properties of exogenous visual stimuli, our results establish macaque monkeys as an ideal model system for unraveling the neural mechanisms underlying a ubiquitous oculomotor phenomenon in visual neuroscience.NEW & NOTEWORTHY Visual onsets dramatically reduce saccade generation likelihood with very short latencies. Such latencies suggest that stimulus-induced visual responses, normally jump-starting perceptual and scene analysis processes, can also directly impact the decision of whether to generate saccades or not, causing saccadic inhibition. Consistent with this, we found that changing the appearance of the visual onsets systematically alters the properties of saccadic inhibition. These results constrain neurally inspired models of coordination between saccade generation and exogenous sensory stimulation.

Safety, immunogenicity and efficacy of Relcovax®, a dual receptor binding domain (RBD) and nucleocapsid (N) subunit protein vaccine candidate against SARS-CoV-2 virus

SARS-CoV-2, severe acute respiratory syndrome coronavirus-2, causes coronavirus disease- 2019 (COVID-19). Mostly, COVID-19 causes respiratory symptoms that can resemble those of a cold, the flu, or pneumonia. COVID-19 may harm more than just lungs and respiratory systems. It may also have an impact on other parts of the body and debilitating effects on humans, necessitating the development of vaccines at an unprecedented rate in order to protect humans from infections. In response to SARS-CoV-2 infection, mRNA, viral vector-based carrier and inactivated virus-based vaccines, as well as subunit vaccines, have recently been developed. We developed Relcovax®, a dual antigen (Receptor binding domain (RBD) and Nucleocapsid (N) proteins) subunit protein vaccine candidate. Preliminary mouse preclinical studies revealed that Relcovax® stimulates cell-mediated immunity and provides broader protection against two SARS-CoV-2 variants, including the delta strain. Before conducting human studies, detailed preclinical safety assessments are required, so Relcovax® was tested for safety, and immunogenicity in 28-day repeated dose toxicity studies in rats and rabbits. In the toxicity studies, there were no mortality or morbidity, abnormal clinical signs, abnormalities in a battery of neurobehavioral observations, abnormalities in detailed clinical and ophthalmological examinations, or changes in body weights or feed consumption. In any of the studies, no abnormal changes in organ weights, haematology, clinical chemistry, urinalysis parameters, or pathological findings were observed. Immunogenicity tests on rats and rabbits revealed 100 % seroconversion. Relcovax® was therefore found to be safe in animals, with a No Observed Adverse Effect Level (NOAEL) of 20 µg/protein in rats and rabbits. In efficacy studies, Relcovax® immunised hamsters demonstrated dose-dependent protection against SARS-CoV-2 infection, with a high dose (20 µg/protein) being the most protective, while in cynomolgus macaque monkey study, lowest dose 5 µg/protein had the highest efficacy. In conclusion, Relcovax® was found to be safe, immunogenic, and efficacious in in vivo studies.

Variable cardiac responses in rhesus macaque monkeys after discrete mediodorsal thalamus manipulations

The control of some physiological parameters, such as the heart rate, is known to have a role in cognitive and emotional processes. Cardiac changes are also linked to mental health issues and neurodegeneration. Thus, it is not surprising that many of the brain structures typically associated with cognition and emotion also comprise a circuit—the central automatic network—responsible for the modulation of cardiovascular output. The mediodorsal thalamus (MD) is involved in higher cognitive processes and is also known to be connected to some of the key neural structures that regulate cardiovascular function. However, it is unclear whether the MD has any role in this circuitry. Here, we show that discrete manipulations (microstimulation during anaesthetized functional neuroimaging or localized cytotoxin infusions) to either the magnocellular or the parvocellular MD subdivisions led to observable and variable changes in the heart rate of female and male rhesus macaque monkeys. Considering the central positions that these two MD subdivisions have in frontal cortico-thalamocortical circuits, our findings suggest that MD contributions to autonomic regulation may interact with its identified role in higher cognitive processes, representing an important physiological link between cognition and emotion.