Preferential Neuronal Responses to Snakes in the Monkey Medial Prefrontal Cortex Support an Evolutionary Origin for Ophidiophobia

Abstract

Ophidiophobia (snake phobia) is one of the most common specific phobias. It has been proposed that specific phobia may have evolutionary origin, and that attentional bias to specific items may promote onset of phobia. Noninvasive imaging studies of patients with specific phobia reported that the medial prefrontal cortex (mPFC), especially the rostral part of the anterior cingulate cortex (rACC), and amygdala are activated during presentation of phobogenic stimuli. We propose that the mPFC-amygdala circuit may be involved in pathogenesis of phobia. The mPFC receives inputs from the phylogenically old subcortical visual pathway including the superior colliculus, pulvinar and amygdala, while mPFC neurons are highly sensitive to snakes that are the first modern predator of primates, and discriminate snakes with striking postures from those with non-striking postures. Furthermore, the mPFC has been implicated in attentional allocation, and promotes amygdala-dependent aversive conditioning. These findings suggest that the rACC focuses attention on snakes, and promotes aversive conditioning to snakes, which may lead to anxiety and ophidiophobia.

Introduction

Specific phobia is a type of anxiety disorder with a lifetime prevalence ranging from 3% to 15% across various countries (Eaton et al., 2018). Patients with specific phobia display immediate excessive fear or anxiety in response to a specific object or situation despite the absence of actual danger imposed by the specific object or situation, and the disturbances by specific phobia impose a significant impact on patients’ daily activity (American Psychiatric Association, 2013). Specific phobia consists of five types, and ophidiophobia (snake phobia) belongs to an animal type of specific phobia. Psychological studies reported that 53.3% of survey participants felt anxiety in response to snakes (Davey, 1994), and 2–3% of the participants were ophidiophobia-like (Klorman et al., 1974; Klieger, 1987; Polák et al., 2016), and that ophidiophobia is one of the most common specific phobias (Fredrikson et al., 1996).

It has been proposed that this specific phobia may have an evolutionary origin (Mineka and Öhman, 2002; Russell et al., 2015; Rakison, 2018). The fear module of primates, including the superior colliculus, pulvinar, and amygdala, may have evolved particularly in response to the threat posed by snakes (Isbell, 1994, 2006, 2009; Öhman and Mineka, 2003; Soares et al., 2017). Consistently, monkey pulvinar neuronal activity is more sensitive to snake images than other control images (Le et al., 2013, 2014, 2016). Furthermore, human and monkey behavioral studies support the existence of an innate fear module sensitive to snake images. Human infants (5–6 months old) who likely had less experience of snakes looked at snake images longer than control images and displayed stronger sympathetic responses to snake images (Hoehl et al., 2017; Rakison, 2018), while young children and adults detected snakes faster than control images (LoBue and DeLoache, 2008; Masataka et al., 2010). In monkeys, isolation-reared (snake-naïve) mouse lemurs and pig-tailed macaques avoided a snake model and odor of snakes (Nelson et al., 2003; Weiss et al., 2015). Also, humans are predisposed to efficiently learn aversive conditioning between snakes and unconditioned shock (Seligman, 1970; Öhman and Mineka, 2003). These findings suggest that the innate fear module for snake detection and learning may be involved in snake phobia.

It has been proposed that attentional bias to specific animals may promote anxiety (i.e., phobia; Matthews and MacLeod, 2002; Heeren et al., 2013; LoBue and Rakison, 2013): images of specific animals are automatically processed in a fear module regardless of attention (Öhman and Soares, 1994), compete with cognitive activity, and when the activity for the specific animals exceeds cognitive activity, capture attention for awareness to induce anxiety. Consistent with this hypothesis, snake images, which are prone to induce phobia (see above), captured more attention than other control images in monkeys and humans (Shibasaki and Kawai, 2009; Soares et al., 2014; Kawai and Koda, 2016; Gomes et al., 2017). It is suggested that the anterior cingulate cortex (ACC) is involved in attentional allocation to emotional stimuli (Niedenthal and Kitayama, 1994) and that the affective (rostral) part of the ACC (rACC) is involved in salience evaluation of emotional stimuli (Bush et al., 2000). Furthermore, phobogenic stimuli were detected faster than control stimuli in subjects with strong fear comparable to specific phobia and patients with phobia (Öhman et al., 2001; Bar-Haim et al., 2007), while the rACC was activated in response to phobogenic stimuli including snakes in patients with specific phobia (Rauch et al., 1995; Pissiota et al., 2003; Carlsson et al., 2004; Amir et al., 2005; Straube et al., 2007; Britton et al., 2009). Thus, the currently available data suggest that the ACC, especially rACC, plays an important role in the pathogenesis of ophidiophobia.

Based on evolutional theories for the origin of phobias (see above), we hypothesized that the ACC, which receives projections from the pulvinar (Porrino et al., 1981; Romanski et al., 1997), is also sensitive to snakes. Consistently, neurons in the medial prefrontal cortex (mPFC) including the ACC responded faster and stronger to snakes in monkeys (Dinh et al., 2018). Furthermore, behavioral studies reported that monkeys, as well as humans, are more sensitive to striking postures in snake detection and threat assessment (Masataka et al., 2010; Etting et al., 2014). Based on these behavioral findings and evolutionary reasoning, here we tested whether mPFC neurons are more responsive to images of snakes with striking postures in non-human primates.

The Behavioral Task for Neurophysiological Recording from The Monkey mPFC

To test this prediction, we recorded and analyzed monkey mPFC neuronal responses to snake images with striking and non-striking postures. The recording was performed using two adult macaque monkeys (Macaca fuscata: 1 female and 1 male) weighing 7.1–8.6 kg, which was the same monkeys used in the previous experiment (Dinh et al., 2018). In the previous study (Dinh et al., 2018), the responses of mPFC neurons to snake images regardless of their postures were compared with other categories of images (e.g., raptors, carnivores, humans, and monkey faces, et cetera). In this study, the neuronal responses to snake images with striking and non-striking postures were analyzed. While the monkeys were individually housed with food available ad libitum, the monkeys’ water intake was restricted in their home cages and they obtained juice rewards during training as well as recording experiments (Dinh et al., 2018). The monkeys were treated according to the United States Public Health Service Policy on Human Care and Use of Laboratory Animals, the National Institutes of Health Guide for the Care and Use of Laboratory Animals, and the Guidelines for the Care and Use of Animals at the University of Toyama. Experimental procedures in this study have been approved by the Ethical Committee for Animal Experiments at the University of Toyama (Dinh et al., 2018).

The detailed procedures are described in Dinh et al. (2018). Briefly, the monkey sat in a monkey chair in front of a computer display. Eye-movements were monitored by a CCD camera (Matsuda, 1996). In the present study, six snake photos consisting of three snake images with striking postures and another three snake images with non-striking postures were used (Figure 1). The present study used the same visual stimuli in the previous study in the monkey pulvinar (Le et al., 2014). The monkeys discriminated against the snake images in a sequential delayed nonmatching-to-sample task (DNMS; Dinh et al., 2018). Briefly, the task started with a buzzer tone, followed by a fixation cross for 1.5 s. Then, a sample stimulus was presented for 500 ms (sample phase) between 1 and 4 times with a 1.5-s interval. After a 1.5-s interval from the last presentation of the sample stimulus, a new stimulus (target phase) was presented. The monkey could obtain a juice reward if the monkey pressed a button in the chair within 2 s.

FIGURE 1

Figure 1. Visual stimuli used in a neurophysiological recording from the monkey medial prefrontal cortex (mPFC). (A) Three photos of snakes with striking postures. (B) Three photos of snakes with non-striking postures. The same visual stimuli used in the study by Le et al. (2014) were also used in the present study.

Upon reaching a 96% correct-response criterion in the DNMS task, a U-shaped epoxy-resin plate was attached to the skull under anesthesia (Nishijo et al., 1988; Tazumi et al., 2010; Dinh et al., 2018). The monkeys were retrained with the DNMS task two weeks after the surgery, while the head was painlessly fixed using a head-restraining device for the U-shaped plate. After re-training of the monkeys with the head fixed, neuronal activity was recorded from the mPFC during the DNMS task. The detailed procedures for electrophysiological recording were reported previously (Dinh et al., 2018). Briefly, a glass-coated tungsten microelectrode was stereotaxically inserted into the mPFC to record neuronal activity. All the data, including analog signals of neuronal activities, were stored in a computer off-line analysis.

Characteristics of mPFC Neuronal Responses to Snake Postures

Only neuronal responses to sample stimuli were analyzed (Dinh et al., 2018). Briefly, significant responses to each snake image were determined by a Wilcoxon signed–rank test (p < 0.05), in which neuronal activity during 500-ms after (post) stimulus onset was compared with that during 100-ms before (pre) stimulus onset. The mPFC neurons were defined as snake-responsive neurons if a given mPFC neuron showed significant responses to at least one of the six snake images. The 235 mPFC neurons were tested with all six snake images in the DNMS task, and 95 neurons responded to one or more snake images (snake-responsive neurons). The numbers of mPFC neurons and snake-responsive neurons in the individual monkeys are indicated in Supplementary Table 1. There were no significant differences in the proportion of snake-responsive neurons among the mPFC neurons between the two monkeys (χ²-test; χ2(1) χ(1)2 = 1.147, p = 0.284). Figure 2 shows a representative snake-responsive mPFC neuron. The neuron responded stronger to snakes with striking postures (Figure 2A). For each snake-responsive neuron, response magnitudes to each stimulus were computed [i.e., (mean firing rates during the post 500 ms period) minus (mean firing rates in the pre 100 ms period)]. There was a significant difference in the response magnitudes among the six snake images in this neuron (one-way ANOVA: F_(5,63) = 6.215, p < 0.0001; Figure 2B). Subsequent post hoc tests indicated that response magnitudes to the three snake images with striking postures were larger than those to the three snake images with non-striking postures (Tukey test, p < 0.05).

FIGURE 2

Figure 2. A representative mPFC neuron sensitive to snake postures. (A) Peri-event histograms of the mPFC neuronal activity in response to each snake photo. Each raster display above each histogram indicates neuronal activity. The red horizontal bars above the raster display indicate the stimulus presentation period (500 ms). Zero in the abscissas indicates the stimulus onset. Calibration at the right bottom of the figure indicates the number of spikes per trial in each bin. Bin width = 50 ms. (B) Response magnitudes of this neuron to the six snake images. Histograms indicate mean ± SEM. *p < 0.05.

Snake-responsive neurons were further categorized based on comparison of response magnitudes to the snake images with striking postures with those to non-striking postures. For this categorization, two peri-event histograms, one for the snakes with striking postures and the other for the snakes with non-striking postures, were constructed. The mPFC neurons were categorized as striking-selective neurons if the response magnitudes to the three snake images with striking postures were larger than those to the three snake images with non-striking postures (unpaired t-test, p < 0.05), while the mPFC neurons were categorized as non-striking-selective neurons if the response magnitudes to the three snake images with non-striking postures were larger than those to the three snake images with striking postures (unpaired t-test, p < 0.05). The remaining mPFC neurons were defined as nonselective neurons. Of the 95 snake-responsive mPFC neurons, 64 neurons responded stronger to snakes with striking postures (striking-selective), 20 neurons responded stronger to snakes with non-striking postures (non-striking-selective), and 11 neurons showed no difference in response magnitudes between these two categories of the snake images (non-selective; Figure 3A). Since there were no significant differences in the proportion of striking-selective neurons among the mPFC neurons between the two monkeys (χ²-test; χ2(1) χ(1)2 = 0.095, p = 0.758), the combined data from the two monkeys were further analyzed. The ratio of the striking-selective neurons was significantly greater than that of the non-striking-selective neurons (χ²-test; χ2(1) χ(1)2 = 20.51, p < 0.0001). Furthermore, averaged response magnitudes of all snake-responsive neurons to the snakes with striking postures and those to the snakes with non-striking postures were compared by paired t-tests (p < 0.05). The results indicated that the averaged response magnitudes to the snakes with striking postures were significantly larger than those to the snakes with non-striking postures (paired t-test: t₍₉₄₎ = 5.459, p < 0.0001; Figure 3B).

FIGURE 3

Figure 3. Classification (A) and response characteristics (B,C) of snake-responsive neurons. (A) The ratios of three types of mPFC neurons to the snake postures. ****Significant difference from non-striking-selective mPFC neurons (p < 0.0001). (B) Comparison of mean response magnitudes to snakes with striking postures and those to non-striking postures (n = 95). ****Significant difference from snakes with non-striking postures (p < 0.0001). (C) Comparison of mean response latencies to snakes with striking postures and those to non-striking postures (n = 87). ***Significant difference from snakes with non-striking postures (p < 0.001).

Response latencies to the snake images were also analyzed. For each neuron, the same two data sets used to analyze response magnitudes for snakes with striking and non-striking postures (see above) were also used to analyze response latencies. Response latency of the mPFC neuron was measured as the period between the stimulus onset and the time point at which the neuronal firing rate exceeded the mean ± 2 SD of the baseline firing rate. The mean response latencies of all striking-selective neurons were compared to non-striking-selective neurons by paired t-tests (p < 0.05; n = 87 for latency measures). The results indicated that averaged response latencies to the snakes with striking postures were significantly shorter than those to the snakes with non-striking postures (paired t-test: t₍₈₆₎ = 3.804, p = 0.0003; Figure 3C).

Population Coding of Snake Postures in The Monkey mPFC

To investigate temporal changes in population coding of snake postures in the mPFC, the initial 200-ms after stimulus onset was divided into four 50-ms epochs: epoch 1 (0–50 ms), epoch 2 (50–100 ms), epoch 3 (100–150 ms) and epoch 4 (150–200 ms). The mean response magnitude in each epoch was similarly calculated, as follows: (mean firing rate in each epoch) minus (mean firing rate during the 100-ms pre-period). These four data sets were separately analyzed by multidimensional scaling (MDS). MDS can convert the relationships (Euclidean distances as dissimilarity in this study) between all possible pairs of stimuli to a geometric representation of the stimuli in a space (see Young and Hamer, 1987 for more details). In this study, data matrices of response magnitudes in a 95 × 6 array (95 snake-responsive neurons × 6 stimuli) were generated in each epoch and were subjected to MDS analysis. Then, six visual stimuli were positioned in the two-dimensional space so that the geometrical relationships among the stimuli represented the original relationships (MDS, PROXSCAL procedure, SPSS statistical package, version 16; Shepard, 1962; Kruskal, 1964).

In the 2-D MDS spaces, r² values of epochs 1, 2, 3 and 4 were 0.671, 0.819, 0.953, and 0.762, respectively. The distributions of the six snake images in 2D space derived from each epoch data are shown in Figure 4. In epoch 2 and 3, two clusters, one for snakes with striking postures and another for snakes with non-striking postures, were recognized. Discriminant analyses of the MDS data indicated that the correct percent of discrimination between the two categories of the snake images was 97.32% and 100% in epochs 2 and 3, respectively (p = 0.016 and p = 0.012, respectively; Table 1). In epochs 1 and 4, the two categories of the snake images were not significantly separated.

FIGURE 4

Figure 4. Distributions of the six snake photos in a 2-D space derived from multidimensional scaling (MDS) analysis of response magnitudes of the 95 mPFC neurons to these photos. (A–D) Response magnitudes were analyzed by MDS epoch 1 (A), epoch 2 (B), epoch 3 (C), and epoch 4 (D). In epochs 2 and 3 (B,C), the snakes with striking postures were separated from the snakes with non-striking postures (p < 0.05).

TABLE 1

Table 1. Discriminant analyses of the multidimensional scaling (MDS) results in the medial prefrontal cortex (mPFC).

A previous study reported that response magnitudes and latencies in the mPFC were positively correlated with the corresponding data in the pulvinar (Dinh et al., 2018). Furthermore, in the previous study in which responses to the same snake images with striking and non-striking postures were analyzed in the pulvinar, snakes with striking postures were significantly separated from those with non-striking postures in epoch 2 in the MDS analysis (Le et al., 2014). This study investigated how the MDS data for the snakes with striking and non-striking postures in the mPFC (this study) were related to those in the pulvinar (Le et al., 2014). To investigate the relationships between the two MDS results, we analyzed the relationships between the Euclidean distances between all stimulus pairs in MDS spaces in epochs 2 and 3 in the mPFC and those in the pulvinar data in epoch 2 (Le et al., 2014) by linear regression analysis. The results indicated that the distances between all possible stimulus pairs in the mPFC MDS in epoch 2 were significantly and positively correlated with the corresponding data in the pulvinar MDS space in epoch 2 (r² = 0.717; F_(1,13) = 32.97, p < 0.0001), and that the distances in the mPFC MDS in epoch 3 were also significantly and positively correlated with the corresponding data in the pulvinar MDS space in epoch 2 (r² = 0.754; F_(1,13) = 39.74, p < 0.0001). These findings suggest that the mPFC receives information from the pulvinar to discriminate snake postures.

Locations of The mPFC Neurons

The histological procedures to identify recording sites in the mPFC were reported previously (Dinh et al., 2018). Briefly, after the last recording session, the monkeys were perfused under deep anesthesia, and brain sections were processed with Cresyl violet. Then, recording sites of mPFC neurons were stereotaxically determined based on the stereotaxic brain atlas of monkeys (Macaca fuscata; Kusama and Mabuchi, 1970). The mPFC was divided into three subareas (Dinh et al., 2018); the area anterior to the cingulate sulcus (anterior mPFC), area dorsal to the fundus of the cingulate sulcus (dorsal mPFC), and area ventral to the fundus of the cingulate sulcus (ventral mPFC; Figure 5). Consistent with the previous study in the mPFC (Dinh et al., 2018), the snake-responsive neurons were located mainly in the ventral mPFC (Figure 5). Especially, the striking-selective neurons were located in the ventral mPFC corresponding to the pregenual part of the ACC (i.e., rACC).

FIGURE 5

Figure 5. Locations of the three types of snake-responsive neurons in the mPFC. (A–H) Locations of mPFC neurons plotted in coronal sections at different A-P levels. The number under each section, distance (mm) anteriorly from the interaural line. Open circles, striking-selective neurons (n = 64); open squares, non-striking-selective neurons (n = 20); open triangles, nonselective neurons (n = 11); dots, nonresponsive neurons (n = 140). The mPFC was divided into three subareas: an anterior part of the mPFC (white areas in A,B), a dorsal part of the mPFC (dark gray areas in C–H), and a ventral part of the mPFC (light gray areas in C–H). PS, principal sulcus; LOS, lateral orbital sulcus; MOS, medial orbital sulcus; CIS, cingulate sulcus.

Discussion

In the present study, we showed that the monkey mPFC neurons responded stronger and faster to snakes with striking postures compared with snakes with non-striking postures. Furthermore, population activity patterns of the mPFC neurons discriminated against snakes with striking postures from those with non-striking postures in epochs 2 and 3. Ethological studies reported that the striking speed of snakes is fast enough to preclude motor escape responses from mammal species (Penning et al., 2016). Given that biting strikes are generally preceded by striking postures (Arnold and Bennett, 1984; Greene, 1988), fast discrimination of snake posture may be critical to surviving agonistic encounters with snakes. Consistently, monkeys respond strongly to snake models with striking postures (Etting et al., 2014). We previously reported that mPFC and pulvinar neurons responded stronger and quicker to snakes compared with other animals (Le et al., 2013; Dinh et al., 2018), and that pulvinar neurons also discriminated against snakes with striking postures from those with non-striking postures (Le et al., 2014). Furthermore, MDS locations of snakes with and without striking postures in the mPFC were significantly correlated with those in the pulvinar, suggesting that the mPFC receives inputs directly from the pulvinar or indirectly via the amygdala since the pulvinar also sends projections to the amygdala (Rafal et al., 2015; Diano et al., 2017; Elorette et al., 2018), and the pulvinar and amygdala send projections to the mPFC (Porrino et al., 1981; Romanski et al., 1997). These results suggest that the mPFC and subcortical visual pathways are particularly sensitive to evolutionary threats (striking postures of snakes), which is consistent with the Snake Detection Theory, in which predations of primates by snakes are important selective forces to shape the primate brain to be sensitive to snakes (Isbell, 2006, 2009). These findings support the idea of the evolutionary origin of specific phobia, especially ophidiophobia (see “Introduction” section).

The striking-selective neurons were located mostly in the mPFC area roughly corresponding to the pregenual part of the ACC. Human imaging studies reported that phobogenic stimuli elicited activation in the rACC that roughly corresponds to the pregenual part of the ACC (see “Introduction” section), and that functional connectivity between the rACC and amygdala increased (Gold et al., 2015) and activity in the rACC covaried with activity in the amygdala (Carlson et al., 2009) in presence of threat stimuli. Previous electrophysiological studies reported that rACC neurons showed excitatory responses to aversive cues and inhibitory responses to rewarding cues (Takenouchi et al., 1999; Amemori and Graybiel, 2012) and that microstimulation of the rACC increased bias for avoidance behavior (Amemori and Graybiel, 2012), suggesting that activity in the rACC may be anxiogenic. Consistently, a meta-analysis of imaging studies of patients and at-risk groups suggests that the pregenual ACC is associated with the expression of emotional psychopathology (Marusak et al., 2016). Furthermore, cortical thickness in the pregenual ACC was increased in patients with animal phobia including ophidiophobia compared with healthy controls (Rauch et al., 2004). These findings suggest that the rACC, along with the amygdala, is an important brain region for the development of ophidiophbia.

Attentional bias to specific animals promotes anxiety and phobia (see “Introduction” section). The mPFC is reported to be involved in attentional allocation to various salient visual stimuli (Dalley et al., 2004; Guillem et al., 2011; Kaping et al., 2011; Mao et al., 2017). These findings suggest that neurons in the rACC, which receive inputs from the amygdala and pulvinar, may guide attention to focus on snakes. Consistent with this idea, gray matter volume in the rACC was positively correlated with covert attention to threatening stimuli (Carlson et al., 2012). The increases in attention bias to snakes might be associated with the pathogenesis of ophidiophobia. Human behavioral studies reported that training anxious subjects to guide attention to non-threatening faces from threatening faces reduced anxiety (Heeren et al., 2013). The same method may be applied to patients to treat ophidiophobia. On the other hand, the pregenual ACC (rACC) has bidirectional (reciprocal) connections with the amygdala (Calderazzo et al., 2020), and these bidirectional connections are excitatory (Laviolette et al., 2005; Bissière et al., 2008). It is reported that activity in the rACC is required to develop aversive conditioning in the amygdala (Laviolette et al., 2005; Bissière et al., 2008). These findings suggest that rACC neurons responsive to snakes may guide aversive conditioning in the amygdala. A human behavioral study using aversive conditioning reported that the association between snakes and unconditioned shock was more strongly formed than the association between other control objects and shock (Öhman and Mineka, 2003). These findings suggest that the rACC might initially and preferentially respond to snakes through direct inputs from the pulvinar. The rACC is also likely to facilitate amygdala-dependent conditioning processes, which in turn might further facilitate rACC responses to snakes. The increased activity in the rACC may promote attention to snakes, which may increase the likelihood of developing ophidiophobia. Taken together, these findings highlight the role of the rACC-amygdala circuits in the pathogenesis of specific phobia.

Limitation

First, the responses to snakes could depend on low-grade visual features of the snake stimuli, but not on their postures. However, our previous study reported that the responses to snakes in the mPFC were decreased by scrambling of the stimuli using two different methods (quadrate-scrambling, in which snake images were cut into small fragments, and the fragments were randomly reassembled, and Fourier-scrambling, in which shapes were degraded while the global low-level visual features of the original snake images were maintained) (Dinh et al., 2018). These findings suggest that mPFC neurons responded to snake shapes (i.e., snake postures), but not to low-grade visual features of snake images.

Second, the mPFC neurons were recorded from the two monkeys. Since the proportions of snake-responsive and striking-selective neurons among the mPFC neurons were similar between the two monkeys (Supplementary Table 1), the data from the two monkeys were combined for further analyses. Although this way of analysis (analyses of combined data of neuronal responses recorded from a few animals) is relatively common in monkey neurophysiological studies (e.g., Tsao et al., 2006; Tsujimoto et al., 2010; Cavanaugh et al., 2020), further studies using more animals are required to test the present hypothesis.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The animal study was reviewed and approved by the Ethical Committee for Animal Experiments in University of Toyama.

Author Contributions

HNishij conceived the study and designed the experiment. HD and QL performed the experiment. QL and HNishij analyzed the data and wrote the article. HNishij, HNishim, JM, TS, RM, CT, and TO revised the article. All authors contributed to the article and approved the submitted version.

Funding

The study was supported partly by Grant-in-Aid for Scientific Research (B) (16H04652 and 20H03417) from the Japan Society for Promotion of Science (JSPS), Takeda Science Foundation, and research grants from the University of Toyama.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor declared a past co-authorship with one of the authors RM.

Acknowledgments

Animals were provided by NBRP “Japanese Monkeys” through the National BioResource Project (NBRP) of the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnbeh.2021.653250/full#supplementary-material.

References

American Psychiatric Association. (2013). Diagnostic and Statistical Manual of Mental Disorders, 5th edition. Arlington: American Psychiatric Association.

Amemori, K., and Graybiel, A. M. (2012). Localized microstimulation of primate pregenual cingulate cortex induces negative decision-making. Nat. Neurosci. 15, 776–785. doi: 10.1038/nn.3088

Amir, N., Klumpp, H., Elias, J., Bedwell, J. S., Yanasak, N., and Miller, L. S. (2005). Increased activation of the anterior cingulate cortex during processing of disgust faces in individuals with social phobia. Biol. Psychiatry 57, 975–981. doi: 10.1016/j.biopsych.2005.01.044

Arnold, S. J., and Bennett, A. F. (1984). Behavioral variation in natural populations III: antipredator display in the garter snake Thamnophis radix. Anim. Behav. 32, 1108–1118. doi: 10.1016/S0003-3472(84)80227-4

Bar-Haim, Y., Lamy, D., Pergamin, L., Bakermans-Kranenburg, M. J., and van Ijzendoorn, M. H. (2007). Threat-related attentional bias in anxious and nonanxious individuals: a meta-analytic study. Psychol. Bull. 133, 1–24. doi: 10.1037/0033-2909.133.1.1

Bissière, S., Plachta, N., Hoyer, D., McAllister, K. H., Olpe, H. R., Grace, A. A., et al. (2008). The rostral anterior cingulate cortex modulates the efficiency of amygdala-dependent fear learning. Biol. Psychiatry 63, 821–831. doi: 10.1016/j.biopsych.2007.10.022

Britton, J. C., Gold, A. L., Deckersbach, T., and Rauch, S. L. (2009). Functional MRI study of specific animal phobia using an event-related emotional counting stroop paradigm. Depress. Anxiety. 26, 796–805. doi: 10.1002/da.20569

Bush, G., Luu, P., and Posner, M. I. (2000). Cognitive and emotional influences in anterior cingulate cortex. Trends Cogn. Sci. 4, 215–222. doi: 10.1016/s1364-6613(00)01483-2

Calderazzo, S. M., Busch, S. E., Moore, T. L., Rosene, D. L., and Medalla, M. (2020). Distribution and overlap of entorhinal, premotor and amygdalar connections in the monkey anterior cingulate cortex. J. Comp. Neurol. 529, 885–904. doi: 10.1002/cne.24986

Carlson, J. M., Beacher, F., Reinke, K. S., Habib, R., Harmon-Jones, E., Mujica-Parodi, L. R., et al (2012). Nonconscious attention bias to threat is correlated with anterior cingulate cortex gray matter volume: a voxel-based morphometry result and replication. NeuroImage 59, 1713–1718. doi: 10.1016/j.neuroimage.2011.09.040

Carlson, J. M., Reinke, K. S., and Habib, R. (2009). A left amygdala mediated network for rapid orienting to masked fearful faces. Neuropsychologia 47, 1386–1389. doi: 10.1016/j.neuropsychologia.2009.01.026

Carlsson, K., Petersson, K. M., Lundqvist, D., Karlsson, A., Ingvar, M., and Ohman, A. (2004). Fear and the amygdala: manipulation of awareness generates differential cerebral responses to phobic and fear-relevant (but nonfeared) stimuli. Emotion 4, 340–353. doi: 10.1037/1528-3542.4.4.340

Cavanaugh, J., McAlonan, K., and Wurtz, R. H. (2020). Organization of corollary discharge neurons in monkey medial dorsal thalamus. J. Neurosci. 40, 6367–6378. doi: 10.1523/JNEUROSCI.2344-19.2020

Dalley, J. W., Theobald, D. E., Bouger, P., Chudasama, Y., and Cardinal, R. N. (2004). Cortical cholinergic function and deficits in visual attentional performance in rats following 192 IgG-saporin-induced lesions of the medial prefrontal cortex. Cereb. Cortex 14, 922–932. doi: 10.1093/cercor/bhh052

Davey, G. C. (1994). Self-reported fears to common indigenous animals in an adult UK population: the role of disgust sensitivity. Br. J. Psychol. 85, 541–554. doi: 10.1111/j.2044-8295.1994.tb02540.x

Diano, M., Celeghin, A., Bagnis, A., and Tamietto, M. (2017). Amygdala response to emotional stimuli without awareness: facts and interpretations. Front. Psychol. 7:2029. doi: 10.3389/fpsyg.2016.02029

Dinh, H. T., Nishimaru, H., Matsumoto, J., Takamura, Y., Le, Q. V., Hori, E., et al. (2018). Superior neuronal detection of snakes and conspecific faces in the macaque medial prefrontal cortex. Cereb. Cortex 28, 2131–2145. doi: 10.1093/cercor/bhx118

Eaton, W. W., Bienvenu, O. J., and Miloyan, B. (2018). Specific phobias. Lancet Psychiatry 5, 678–686. doi: 10.1016/S2215-0366(18)30169-X

Elorette, C., Forcelli, P. A., Saunders, R. C., and Malkova, L. (2018). Colocalization of tectal inputs with amygdala-projecting neurons in the macaque pulvinar. Front. Neural Circuits 12:91. doi: 10.3389/fncir.2018.00091

Etting, S. F., Isbell, L. A., and Grote, M. N. (2014). Factors increasing snake detection and perceived threat in captive rhesus macaques (Macaca mulatta). Am. J. Primatol. 76, 135–145. doi: 10.1002/ajp.22216

Fredrikson, M., Annas, P., Fischer, H., and Wik, G. (1996). Gender and age differences in the prevalence of specific fears and phobias. Behav. Res. Ther. 34, 33–39. doi: 10.1016/0005-7967(95)00048-3

Gold, A. L., Morey, R. A., and McCarthy, G. (2015). Amygdala-prefrontal cortex functional connectivity during threat-induced anxiety and goal distraction. Biol. Psychiatry 77, 394–403. doi: 10.1016/j.biopsych.2014.03.030

Gomes, N., Silva, S., Silva, C. F., and Soares, S. C. (2017). Beware the serpent: the advantage of ecologically-relevant stimuli in accessing visual awareness. Evol. Hum. Behav. 38, 227–234. doi: 10.1016/j.evolhumbehav.2016.10.004

Greene, H. W. (1988). “Antipredator mechanisms in reptiles,” in Biology of the Reptilia (vol. 16), Ecology B, Defense and Life History, eds C. Gans and R. B. Huey (New York, NY: Alan Liss), 1–152.

Guillem, K., Bloem, B., Poorthuis, R. B., Loos, M., Smit, A. B., Maskos, U., et al. (2011). Nicotinic acetylcholine receptor β2 subunits in the medial prefrontal cortex control attention. Science 333, 888–891. doi: 10.1126/science.1207079

Heeren, A., De Raedt, R., Koster, E. H., and Philippot, P. (2013). The (neuro)cognitive mechanisms behind attention bias modification in anxiety: proposals based on theoretical accounts of attentional bias. Front. Hum. Neurosci. 7:119. doi: 10.3389/fnhum.2013.00119

Hoehl, S., Hellmer, K., Johansson, M., and Gredebäck, G. (2017). Itsy bitsy spider…: infants react with increased arousal to spiders and snakes. Front. Psychol. 8:1710. doi: 10.3389/fpsyg.2017.01710

Isbell, L. A. (1994). Predation on primates: ecological patterns and evolutionary consequences. Evol. Anthropol. 3, 61–71. doi: 10.1002/evan.1360030207

Isbell, L. A. (2006). Snakes as agents of evolutionary change in primate brains. J. Hum. Evol. 51, 1–35. doi: 10.1016/j.jhevol.2005.12.012

Isbell, L. A. (2009). The Fruit, the Tree and the Serpent: Why We See So Well. Cambridge, MA: Harvard University Press.

Kaping, D., Vinck, M., Hutchison, R. M., Everling, S., and Womelsdorf, T. (2011). Specific contributions of ventromedial, anterior cingulate and lateral prefrontal cortex for attentional selection and stimulus valuation. PLoS Biol. 9:e1001224. doi: 10.1371/journal.pbio.1001224

Kawai, N., and Koda, H. (2016). Japanese monkeys (Macaca fuscata) quickly detect snakes but not spiders: evolutionary origins of fear-relevant animals. J. Comp. Psychol. 130, 299–303. doi: 10.1037/com0000032

Klieger, D. M. (1987). The snake anxiety questionnaire as a measure of ophidophobia. Educ. Psychol. Meas. 47, 449–459. doi: 10.1177/0013164487472017

Klorman, R., Weerts, T. C., Hastings, J. C., Melamed, B. G., and Lang, P. J. (1974). Psychometric description of some specific-fear questionnaires. Behav. Ther. 5, 401–409. doi: 10.1016/s0005-7894(74)80008-0

Kruskal, J. B. (1964). Multidimensional scaling by optimizing goodness of fit to a nonmetric hypothesis. Psychometrika 29, 1–27. doi: 10.1007/bf02289565

Kusama, T., and Mabuchi, M. (1970). Stereotaxic Atlas of the Brain of Macaca fuscata. Tokyo: Tokyo University Press.

Laviolette, S. R., Lipski, W. J., and Grace, A. A. (2005). A subpopulation of neurons in the medial prefrontal cortex encodes emotional learning with burst and frequency codes through a dopamine D4 receptor-dependent basolateral amygdala input. J. Neurosci. 25, 6066–6075. doi: 10.1523/JNEUROSCI.1168-05.2005

Le, Q. V., Isbell, L. A., Matsumoto, J., Le, V. Q., Nishimaru, H., Hori, E., et al. (2016). Snakes elicit earlier and monkey faces, later, gamma oscillations in macaque pulvinar neurons. Sci. Rep. 6:20595. doi: 10.1038/srep20595

Le, Q. V., Isbell, L. A., Matsumoto, J., Quang, L. V., Tran, A. H., Maior, R. S., et al. (2014). Monkey pulvinar neurons’ differential firing to snake postures. PLoS One 9:e114258. doi: 10.1371/journal.pone.0114258

Le, Q. V., Isbell, L. A., Nguyen, M. N., Matsumoto, J., Hori, E., Maior, R. S., et al. (2013). Pulvinar neurons reveal neurobiological evidence of past selection for rapid detection of snakes. Proc. Natl. Acad. Sci. U S A. 110, 19000–19005. doi: 10.1073/pnas.1312648110

LoBue, V., and DeLoache, J. S. (2008). Detecting the snake in the grass: attention to fear-relevant stimuli by adults and young children. Psychol. Sci. 19, 284–289. doi: 10.1111/j.1467-9280.2008.02081.x

LoBue, V., and Rakison, D. H. (2013). What we fear most: A developmental advantage for threat-relevant stimuli. Dev. Rev. 33, 285–303. doi: 10.1016/j.dr.2013.07.005

Mao, C. V., Araujo, M. F., Nishimaru, H., Matsumoto, J., Tran, A. H., Hori, E., et al. (2017). Pregenual anterior cingulate gyrus involvement in spontaneous social interactions in primates-evidence from behavioral, pharmacological, neuropsychiatric and neurophysiological findings. Front. Neurosci. 11:34. doi: 10.3389/fnins.2017.00034

Masataka, N., Hayakawa, S., and Kawai, N. (2010). Human young children as well as adults demonstrate “superior” rapid snake detection when typical striking posture is displayed by the snake. PLoS One 5:e15122. doi: 10.1371/journal.pone.0015122

Marusak, H. A., Thomason, M. E., Peters, C., Zundel, C., Elrahal, F., and Rabinak, C. A. (2016). You say ‘prefrontal cortex’ and I say ‘anterior cingulate’: meta-analysis of spatial overlap in amygdala-to-prefrontal connectivity and internalizing symptomology. Trans. Psychiatry 6:e944. doi: 10.1038/tp.2016.218

Matsuda, K. (1996). Measurement system of eye positions by using oval fitting of a pupil. Neurosci. Res. Suppl. 25:S270. doi: 10.1016/0168-0102(96)89315-1

Matthews, A., and MacLeod, C. (2002). Induced processing biases have causal effects on anxiety. Cogn. Emot. 16, 331–354. doi: 10.1080/02699930143000518

Mineka, S., and Öhman, A. (2002). Phobias and preparedness: The selective, automatic and encapsulated nature of fear. Biol. Psychiatry 52, 927–937. doi: 10.1016/s0006-3223(02)01669-4

Nelson, E. E., Shelton, S. E., and Kalin, N. H. (2003). Individual differences in the responses of naïve rhesus monkeys to snakes. Emotion 3, 3–11. doi: 10.1037/1528-3542.3.1.3

Niedenthal, P. M., and Kitayama, S. (1994). The Heart’s Eye-Emotional Influences in Perception and Attention. New York, NY: Academic Press.

Nishijo, H., Ono, T., and Nishino, H. (1988). Single neuron responses in amygdala of alert monkey during complex sensory stimulation with affective significance. J. Neurosci. 8, 3570–3583. doi: 10.1523/JNEUROSCI.08-10-03570.1988

Öhman, A., Flykt, A., and Esteves, F. (2001). Emotion drives attention: detecting the snake in the grass. J. Exp. Psychol. Gen. 130, 466–478. doi: 10.1037/0096-3445.130.3.466

Öhman, A., and Mineka, S. (2003). The malicious serpent: snakes as a prototypical stimulus for an evolved module of fear. Curr. Dir. Psychol. Sci. 12, 5–9. doi: 10.1111/1467-8721.01211

Öhman, A., and Soares, J. F. F. (1994). “Unconscious anxiety”: phobic responses to masked stimuli. J. Abnorm. Psychol. 103, 231–240. doi: 10.1037/0021-843x.103.2.231

Penning, D. A., Sawvel, B., and Moon, B. R. (2016). Debunking the viper’s strike: harmless snakes kill a common assumption. Biol. Lett. 12:20160011. doi: 10.1098/rsbl.2016.0011

Pissiota, A., Frans, O., Michelgård, A., Appel, L., Långström, B., Flaten, M. A., et al. (2003). Amygdala and anterior cingulate cortex activation during affective startle modulation: a PET study of fear. Eur. J. Neurosci. 18, 1325–1331. doi: 10.1046/j.1460-9568.2003.02855.x

Polák, J., Sedláčková, K., Nácar, D., Landová, E., and Frynta, D. (2016). Fear the serpent: a psychometric study of snake phobia. Psychiatry Res. 242, 163–168. doi: 10.1016/j.psychres.2016.05.024

Porrino, L. J., Crane, A. M., and Goldman-Rakic, P. S. (1981). Direct and indirect pathways from the amygdala to the frontal lobe in rhesus monkeys. J. Comp. Neurol. 198, 121–136. doi: 10.1002/cne.901980111

Rafal, R. D., Koller, K., Bultitude, J. H., Mullins, P., Ward, R., Mitchell, A. S., et al. (2015). Connectivity between the superior colliculus and the amygdala in humans and macaque monkeys: virtual dissection with probabilistic DTI tractography. J. Neurophysiol. 114, 1947–1962. doi: 10.1152/jn.01016.2014

Rakison, D. H. (2018). Do 5-month-old infants possess an evolved detection mechanism for snakes, sharks and rodents? J. Cogn. Dev. 19, 456–476. doi: 10.1080/15248372.2018.1488717

Rauch, S. L., Savage, C. R., Alpert, N. M., Miguel, E. C., Baer, L., Breiter, H. C., et al. (1995). A positron emission tomographic study of simple phobic symptom provocation. Arch. Gen. Psychiatry 52, 20–28. doi: 10.1001/archpsyc.1995.03950130020003

Rauch, S. L., Wright, C. I., Martis, B., Busa, E., McMullin, K. G., Shin, L. M., et al. (2004). A magnetic resonance imaging study of cortical thickness in animal phobia. Biol. Psychiatry 55, 946–952. doi: 10.1016/j.biopsych.2003.12.022

Romanski, L. M., Giguere, M., Bates, J. F., and Goldman-Rakic, P. S. (1997). Topographic organization of medial pulvinar connections with the prefrontal cortex in the rhesus monkey. J. Comp. Neurol. 379, 313–332.

Russell, M. J., and Maslej, M. M., andrews, P. W. (2015). “Role of evolution,” in Phobias: The Psychology of Irrational Fear, an Encyclopedia, eds I. Milosevic and R. E. McCabe (Goleta, CA: ABC-CLIO), 130–133

Shepard, R. N. (1962). The analysis of proximities: multidimensional scaling with an unknown distance function. Psychometrika 27, 219–246. doi: 10.1007/bf02289621

Shibasaki, M., and Kawai, N. (2009). Rapid detection of snakes by Japanese monkeys (Macaca fuscata): an evolutionarily predisposed visual system. J. Comp. Psychol. 123, 131–135. doi: 10.1037/a0015095

Seligman, M. E. P. (1970). On the generality of the laws of learning. Psychol. Rev. 77, 406–418. doi: 10.1037/h0029790

Soares, S. C., Lindström, B., Esteves, F., and Öhman, A. (2014). The hidden snake in the grass: superior detection of snakes in challenging attentional conditions. PLoS One 9:e114724. doi: 10.1371/journal.pone.0114724

Soares, S. C., Maior, R. S., Isbell, L. A., Tomaz, C., and Nishijo, H. (2017). Fast detector/first responder: interactions between the superior colliculus-pulvinar pathway and stimuli relevant to primates. Front. Neurosci. 11:67. doi: 10.3389/fnins.2017.00067

Straube, T., Mentzel, H. J., and Miltner, W. H. (2007). Waiting for spiders: brain activation during anticipatory anxiety in spider phobics. NeuroImage 37, 1427–1436. doi: 10.1016/j.neuroimage.2007.06.023

Takenouchi, K., Nishijo, H., Uwano, T., Tamura, R., Takigawa, M., and Ono, T. (1999). Emotional and behavioral correlates of the anterior cingulate cortex during associative learning in rats. Neurosci. 93, 1271–1287. doi: 10.1016/s0306-4522(99)00216-x

Tazumi, T., Hori, E., Maior, R. S., Ono, T., and Nishijo, H. (2010). Neural correlates to seen gaze-direction and head orientation in the macaque monkey amygdala. Neurosci. 169, 287–301. doi: 10.1016/j.neuroscience.2010.04.028

Tsao, D. Y., Freiwald, W. A., Tootell, R. B., and Livingstone, M. S. (2006). A cortical region consisting entirely of face-selective cells. Science 311, 670–674. doi: 10.1126/science.1119983

Tsujimoto, S., Genovesio, A., and Wise, S. P. (2010). Evaluating self-generated decisions in frontal pole cortex of monkeys. Nat. Neurosci. 13, 120–126. doi: 10.1038/nn.2453

Weiss, L., Brandl, P., and Frynta, D. (2015). Fear reactions to snakes in naïve mouse lemurs and pig-tailed macaques. Primates 56, 279–284. doi: 10.1007/s10329-015-0473-3

Young, F. W., and Hamer, R. M. (1987). Multidimensional Scaling: History, Theory and Applications. Hillsdale, NJ: Lawrence Erlbaum Associates.