Neural defensive circuits underlie helping under threat in humans

  1. Joana B Vieira  Is a corresponding author
  2. Andreas Olsson
  1. Department of Psychology, Faculty of Health and Life Sciences, University of Exeter, United Kingdom
  2. Department of Clinical Neuroscience, Karolinska Institutet, Sweden

Abstract

Empathy for others’ distress has long been considered the driving force of helping. However, when deciding to help others in danger, one must consider not only their distress, but also the risk to oneself. Whereas the role of self-defense in helping has been overlooked in human research, studies in other animals indicate defensive responses are necessary for the protection of conspecifics. In this pre-registered study (N=49), we demonstrate that human defensive neural circuits are implicated in helping others under threat. Participants underwent fMRI scanning while deciding whether to help another participant avoid aversive electrical shocks, at the risk of also being shocked. We found that higher engagement of neural circuits that coordinate fast escape from self-directed danger (including the insula, PAG, and ACC) facilitated decisions to help others. Importantly, using representational similarity analysis, we found that the strength with which the amygdala and insula uniquely represented the threat to oneself (and not the other’s distress) predicted helping. Our findings indicate that in humans, as other mammals, defensive mechanisms play a greater role in helping behavior than previously understood.

Editor's evaluation

This work fills an important theoretical gap regarding the role of potential threats to the self in prosocial helping. Much of the prevailing knowledge about the motivations for prosocial behavior focuses on the distress of the conspecific-in-need. Leveraging animal research, the authors hypothesize that defensive neural circuitry may stimulate prosocial helping under threat, specifically that "cognitive" fear circuits would respond to more distal threats whereas "reactive" fear circuits would respond to imminent threats. In addition to examining helping behavior under conditions of threat to self, the authors included representational similarity analyses of neuroimaging data, examining how overlapping representations of self and other distress related to helping behavior. The potential to challenge existing empathy accounts of prosocial helping is intriguing and worthy of interrogation.

https://doi.org/10.7554/eLife.78162.sa0

Introduction

Helping someone in danger (e.g., by saving a person who fell on the train tracks, or running into a building in flames to rescue someone inside) may expose oneself to health and life-threatening risks. Nevertheless, such helping behaviors are observed across species (Hernandez-Lallement et al., 2020; Oliner, 2004; Preston, 2013; Rand and Epstein, 2014; de Waal and Lanting, 1997). Risky helping differs from other altruistic actions in that it occurs in the simultaneous presence of two highly salient cues: the distress of a conspecific in need, and a potential threat to the self. In humans, a wealth of research has been dedicated to the former, explaining how perceiving distress in others may trigger the motivation to help (Batson et al., 1987; de Waal and Preston, 2017), particularly if the helper is not under threat themselves (Preston, 2013). But virtually nothing is known about how, in a threatening situation, one’s own responses to the threat may drive decisions to help. More so, animal research suggests defensive brain mechanisms may in fact be implicated in aiding or protecting conspecifics (Ben-Ami Bartal et al., 2016; Bosch, 2013; Rickenbacher et al., 2017). Understanding the neurocognitive processes underlying the motivation to both safeguard oneself and helping others is critical to explain inter-individual behavior in dangerous contexts. The overarching goal of our study was thus to determine how one’s own defensive responses to threat guide decisions to help others in dangerous situations.

In humans, defensive responses to threat are graded as a function of the proximity or imminence of the threatening stimulus, paralleling predatory avoidance responses in other mammals (Fanselow and Lester, 1988; Mobbs et al., 2020). Distal and unpredictable threats are typically associated with risk assessment and intermittent anxiety, allowing for slower and more flexible escape decisions. As threat imminence increases and an attack becomes more likely, fixed and species-specific responses are triggered, such as freezing or, if immediate avoidance is necessary, fight-or-flight. Some behavioral reports indicate that different states along the defensive continuum may have dissociable effects on prosocial behavior. For example, following acute social stress, participants behave more prosocially in economic games (Tomova et al., 2017; von Dawans et al., 2012; von Dawans et al., 2019), make more moral decisions (Singer et al., 2017), and show greater empathy for others (Tomova et al., 2017). Importantly, it has been shown that individuals were more likely to help a co-participant avoid aversive electrical shocks when the threat of shock was imminent rather than distal (Vieira et al., 2020). This behavioral pattern was accompanied by faster reaction times and heart rate during imminent compared to distal threats, paralleling what has been found in response to imminent self-directed threats (Hashemi et al., 2019; Roelofs, 2017). Consistent with these laboratorial studies, higher danger in real-life situations (captured via public surveillance footage) has been associated with higher likelihood of bystander intervention (Lindegaard et al., 2022). Taken together, these findings suggest that defensive states triggered by high threat imminence may not only enable fast avoidance of self-directed threats, but also motivate helping when others are under threat. Yet, the neural basis of these effects is unclear. Specifically, it is unknown how the activation of specific sub-circuits underlying different defensive states (e.g., freezing vs. fight-or-flight) impacts decisions toward others in a threatening context.

We aimed to characterize the involvement of different defensive neural responses on helping under threat. To do so, we used a paradigm adapted from Vieira et al., 2020, in which participants make helping decisions at different stages of threat imminence (details in Figure 1A). Briefly, a participant is asked to decide whether or not to help a co-participant (in reality, a confederate) avoid aversive electrical shocks. In each trial, the participant watches a supposedly live video-feed of the co-participant, and a visual cue signaling an upcoming shock. The participant is asked to decide whether to help the co-participant avoid the shock at the risk of receiving a shock from her/himself. These decisions are prompted in some trials in the beginning of the trial (distal threat), and in others immediately prior to the shock delivery (imminent threat). If the participant decides not to help, the co-participant always receives a shock; if the participant decides to help, there is a fixed probability both participant and co-participant will receive a shock.

Outline of the experimental tasks.

(A) fMRI helping under threat task. Participants saw the co-participant on the screen, together with a visual cue signaling threat (an upcoming shock). There were three threat levels: safe (0 shocks, green circle), moderate threat (1 shock, yellow circle), and high threat (2 shocks, red circle). In each trial of the task, the circle started static on the left (4 s), and then moved to the right (4 s). Participants were prompted to decide whether they wanted to help the co-participant or not (1.25–1.75 s) either in the beginning of trial (distal) or right before the moment of shock delivery (imminent). Therefore, the available time to make a decision was identical in distal and imminent threats. If participants decided to help, there was a 70% chance both themselves and co-participant would receive shocks; if they decided not to help, the co-participant would always receive a shock, and the participant would not. Decisions prompted on safe trials were to arbitrarily choose to press 1 or 2, since no shocks would be administered. (B) After the fMRI task, outside the scanner, participants re-watched clips of the co-participant presented during the scan, and were asked to rate how much ‘discomfort, anxiety or uneasiness’ he was experiencing in each clip on a 9-point scale. They also presented images of the threat cues and asked to rate, on the same scale, how threatened they felt themselves when they saw those images during the scan.

Our hypotheses were guided by previous work on neural responses to the imminence of self-directed threats. It has been shown that the response to distal threats (i.e., unpredictable, spatially distant, retreating, slow moving) is coordinated by so-called ‘cognitive fear’ circuits, which include the ventromedial prefrontal cortex (vmPFC) and hippocampus (Mobbs et al., 2020; Qi et al., 2018). Conversely, imminent threats (i.e., predictable, spatially close, looming, and fast moving) predominantly engage ‘reactive’ fear circuits, which include the dorsal anterior cingulate (dACC), insula, and periaqueductal gray (PAG) (Mobbs et al., 2020; Qi et al., 2018). The amygdala plays a central role in both circuits, namely by coordinating adaptive switches between defensive states as a function of threat imminence (e.g., from freezing to fight-or-flight), through oxytocin-mediated communication between its central (CeA) and basolateral nucleus (BLA) (Terburg et al., 2018; Tovote et al., 2016). Based on these findings, we expected neural activation within the full defensive circuitry to respond to the threat level of the trial (i.e., safe, 1 shock, and 2 shocks), with higher engagement of brain regions previously included in cognitive fear circuits (i.e., vmPFC and hippocampus) in response to distal threats, and higher engagement of regions included in reactive fear circuits (i.e., insula, dACC, and PAG) in response to imminent threats. At the behavioral level, prior findings (Vieira et al., 2020) showed that helping decisions were more frequent under imminent than distal trials, suggesting that the activation of reactive fear circuits would facilitate helping behavior. We thus predicted that helping decisions would be associated with higher engagement of brain regions included in reactive fear circuits (i.e., amygdala, insula, ACC, and PAG).

One important aspect of our paradigm was that, as in most real-life dangerous situations, the threat and the conspecific in need were simultaneously presented. To dissociate the role of representations of threat and of other’s distress on helping behavior, after the scan we asked participants to rate the degree of distress experienced by the co-participant in each clip showed during the scan; also, participants rated how threatened they felt themselves when they saw the visual threat cues during the scan. These ratings were used as behavioral models in a representational similarity analysis (RSA; see Materials and methods) that identified neural representations of other’s distress and of threat to the self, and determined their association with helping behavior. It should be noted that our goal was specifically related to the link between neural representations of distress and threat, and helping behavior, and not to the dissociation of representations of distress and threat in the brain. This would be a highly interesting question to examine in future research.

The demonstration that neural representations of other’s distress are positively associated with helping decisions would support existing empathy-based explanations of helping. Indeed, it has been proposed that helping a conspecific in danger results primarily from an evolutionarily preserved motivation to care for offspring in mammals, which is triggered by signals of distress and vulnerability, and is especially likely to occur if the helper is not under threat themselves (Preston, 2013). However, evidence in rodents indicates that, rather than conflicting, defensive responses may be required for helping and caregiving: for example, anxious rat mothers display enhanced maternal behavior after pharmacological activation of defensive brain circuits (Bosch, 2013; Bosch, 2011; Bosch et al., 2005), whereas mice bred to have low anxiety display significant defects in maternal behaviors (Sheleg et al., 2017); and helping behavior in rats is compromised following treatment with anxiolytic drugs that suppress defensive circuits (Ben-Ami Bartal et al., 2016). According to these animal findings, an alternative prediction to empathy-based accounts is that the neural representation of threat to the self would promote helping of others.

Results

Helping decisions did not vary based on threat imminence or threat level

Using Generalized Linear Mixed Models (GLMMs), we found no significant effect of either threat imminence (β=0.017, se=0.014, t=1.24, p=0.218), threat level (β=−0.068, se=0.053, t=−1.282, p=0.205), nor a threat imminence*level interaction (β=−0.099, se=0.019, t=−0.497, p=0.622) on the percentage of helping decisions throughout the task (Figure 2A). Despite the lack of group-level effects of threat imminence (which were predicted based on previous work; Vieira et al., 2020), individual data showed that the number of participants helping more during imminent than distal threats was objectively higher (n=27) than those helping more during distal (n=11), or the same amount during imminent and distal (n=11) (Figure 2B). Note that these differences are descriptive, and no statistical inference was performed.

Behavioural results.

(A). There was no evidence of differential helping during imminent and distal threats, nor during 1-shock and 2-shock trials. (B) Difference between proportion of helping in imminent and distal trials (y axis) across subjects (x axis); 27 participants helped more during imminent than distal threats, 11 helped more during distal, and 11 helped the same amount. (C) Responses were faster during imminent than distal trials across threat levels. (D) Participants rated the co-participant’s distress as higher during imminent than distal trials, and as progressively higher across the three threat levels. (E) Participants reported feeling more threatened when watching 2-shock cues (red circle), followed by 1-shock cues (yellow circle), and safe cues (green circle), and when watching cues signaling imminent than distal threat.

Previous work also suggested the impact of threat imminence on helping might vary based on empathic tendencies (Vieira et al., 2020). We thus checked whether there was an interaction between the empathic concern scale of the Interpersonal Reactivity Inventory (Davis, 1983) and threat imminence on helping. When empathic concern was added in the model, results showed no significant effects of threat imminence (β=0.018, se=0.014, t=1.249, p=0.215) and threat level (β=−0.050, se=0.051, t=−0.985, p=0.330), and no significant interaction between imminence and empathic concern (β=0.018, se=0.014, t=1.292, p=0.20). However, a significant association emerged between empathic concern and helping behavior (β=−0.101, se=0.041, t=−2.473, p=0.02), indicating those higher in empathic concern displayed less frequent helping behavior. It should be noted here that the impact of low statistical power cannot be discounted. Our previous study (Vieira et al., 2020) had a larger sample, and to detect a correlation between helping and empathic concern of the magnitude of that reported in that study (with 80% power) we would need at least 123 participants (based on calculations carried on in GPower 3.1.9.2).

Also, to account for the possibility that decisions varied throughout the experiment (e.g., participants helped more in the beginning than toward the end), we also performed a mixed effects logistic regression on single trial dichotomous responses (help or no help), including the trial number as a fixed effect. This analysis revealed no significant effects, indicating the individuals did not respond differently as time passed.

Finally, in line with previous work (Vieira et al., 2020), analysis of reaction times showed individuals made faster decisions during imminent versus distal trials (β=−173.30, se=12.09, t=−14.33, p<0.0001), and for shock versus safe trials (β[1 shock]=−43.26, se=12.76, t=−3.39, p=0.001; β[2 shock]=−53.35, se=13.43, t=−3.97, p=0.0002), with no significant threat imminence*level interaction (Figure 2C).

Participants were sensitive to variations in the co-participant’s distress, threat imminence, and threat level (manipulation check)

After the scan, participants were asked to re-watch all clips of the co-participant during the scan, and rate the level of ‘discomfort, anxiety or uneasiness’ they thought he was experiencing in each clip. Of note, these clips were shown without the threat cues (see Figure 1B) that were also present during the scan, in order to isolate the response to distress and threat. Results showed participants rated the distress of the co-participant being significantly higher during imminent than distal clips (β=0.418, se=0.089, t=4.676, p<0.0001), and progressively higher across the three levels of threat (1 Sh: β=0.279, se=0.077, t=3.611, p=0.0005; 2 Sh: β=0.559, se=0.081, t=6.853, p<0.0001; reference class was safe). No significant threat level*imminence interaction was found (Figure 2F). These results suggest the video clips used in the scan successfully portrayed subtle variations in cues of distress by the co-participant.

Participants also presented isolated images of the threat cues used during the scanning task (namely, the green, yellow, and red circles, both in the distal and imminent positions; Figure 1B), and asked to rate how threatened they felt themselves when they saw those visual cues during the scan. Results showed participants rated threat stimuli as more threatening (1 Sh: β=2.255, se=0.297, t=7.598, p<0.0001; 2 Sh: β=3.92, se=0.297, t=13.89, p<0.0001; reference class was safe; Figure 2G). Additionally, imminent cues were rated as more threatening, but only for 1 shock and 2 shocks, and not for safe trials (imminence*1 Sh: β=0.894, se=0.420, t=2.129, p=0.034.; imminence*2 Sh: β=0.936, se=0.420, t=2.23, p=0.027).

Neural responses

Analysis of brain responses focused on a liberally defined set of brain regions that integrate the brain’s defensive system, namely the vmPFC and vlPFC/IFG (Mobbs et al., 2009; Mobbs et al., 2010; Wendt et al., 2017), the hippocampus (Qi et al., 2018), the insula (Mobbs et al., 2010; Wendt et al., 2017), the ACC (Mobbs et al., 2009; Mobbs et al., 2010), the amygdala (Terburg et al., 2018; Mobbs et al., 2009), and the midbrain (Qi et al., 2018; Mobbs et al., 2009; Mobbs et al., 2010; Wendt et al., 2017; Mobbs et al., 2007; Figure 3). The full size of the brain mask used included 23,269 voxels (186,152 mm3).

Combined ROI mask including bilateral ventral and lateral medial frontal cortex, dorsal ACC, insula, hippocampus, amygdala, and midbrain.

ACC, anterior cingulate cortex; ROI, region of interest.

Multivariate and univariate differentiation of threat imminence and level in the defensive circuitry

We performed a support vector machine regression to identify sites in which activation patterns were linearly associated with increasing threat level (i.e., from safe, to 1 shock and 2 shocks). As predicted, results showed that throughout the defensive circuitry (i.e., amygdala, insula, ACC, hipoccampus, and regions within the orbitofrontal cortex) multivariate activation tracked with threat level (FWE<0.05, k>10; Figure 4B).

Multivariate and univariate fMRI results.

(A) Local multivoxel activation patterns (identified by searchlight analysis) in the insula and dmPFC were distinguishable between distal and imminent threats, irrespective of helping decisions. (B) Local multivoxel activation patterns (identified by support vector regression) in the amygdala, insula, OFC/IFG, vmPFC, and ACC were linearly associated with varying threat level. (C) Local multivoxel activation patterns (identified by searchlight) in the insula, IFG, hippocampus, and ACC were distinguishable when making helping decisions under distal and imminent threat. (D) Clusters in the insula, IFG, OFC, and ACC displayed a significant threat imminence*threat level interaction. (E) Clusters in the insula, ACC, IFG/OFC, vmPFC, hippocampus, and PAG displayed a significant type of decision*threat imminence interaction. ACC, anterior cingulate cortex; AMY, amygdala; dmPFC, dorsomedial prefrontal cortex; Hipp, hippocampus; IFG, inferior frontal gyrus; Ins, insula; OFC, orbitofrontal cortex; PAG, periaqueductal gray; vmPFC, ventromedial prefrontal cortex. *p<0.05, **p<0.01, ***p<0.001.

Additionally, to identify brain regions with differential average activation to the imminence of threatening stimuli specifically, we ran a univariate threat imminence (distal, imminent) by threat level (safe, 1 shock, and 2 shocks) ANOVA. Brain regions displaying a significant threat imminence by threat level interaction included the bilateral insula, OFC and IFG and ACC (Table 1; Figure 4D). In shock trials (both 1 shock and 2 shocks), the bilateral insula and ACC presented higher activation for distal compared to imminent threats, whereas the bilateral IFG and right OFC showed higher activation in imminent compared to distal threats (full ANOVA results in Appendix 1—table 1). These results were opposite to our predictions that brain regions previously implicated in reactive fear circuits (i.e., insula and ACC) would be more active during imminent threats, and regions implicated in cognitive fear circuits (i.e., IFG and OFC) during distal threats. Finally, using a Searchlight cross-classification algorithm (12-mm-radius sphere), we also identified brain sites in which multivariate patterns were distinguishable between distal and imminent threats. We found that multivariate patterns in the bilateral insula and dorsal medial prefrontal gyrus dissociated between distal and imminent threats (FWE<0.05, k>10; Figure 4A).

Table 1
Multivariate results based on threat imminence (distal—imminent) and level (safe, 1 shock, 2 shocks; FWE<0.05).
Searchlight distal versus Imminent
R/Lkx, y, zTBA
Insula, superior temporal gyrusR4046, –4, –87.4322
InsulaL17–46, 2, –46.6413
18–42, 4, 87.15
Medial prefrontal cortexR2010, 48, 287.139
SVM regression SF – 1 SH – 2SH
HippocampusR1224, –10, –207.56
InsulaR2542, –8, –88.3413
Rolandic operculumR13952, 0, 07.1222, 47
Superior temporal gyrus, amygdalaR1138, 2, –247.73
AmygdalaL37–24, 4, –247.96
RectusR126, 30, –247.0711
Anterior cingulateR138, 38, 87.02
Middle orbital frontal gyrusR4232, 44, –147.3211
Middle orbital frontal gyrusL11–4, 48, –106.7511
44–28, 54, –107.59
4140, 56, –47.18
Superior frontal orbital gyrusR1916, 54, –146.7810
Medial frontal gyrusR306, 58, –86.9710
Searchlight help during distal versus imminent threats
R2936, –8, –168.04
HippocampusR2150, 4, –27.0122
InsulaR2332, 14, 147.1813
Inferior frontal gyrusL207–48, 20, –68.8147, 38, 22, 13, 45
Mid cingulate, dorsal anterior cingulateR432, 22, 307.2532, 9, 24, 6
InsulaL45–32, 26, 07.8347, 13, 45
Anterior cingulateR5310, 40, 87.9232, 10

Greater engagement of reactive fear circuits led to helping

To test our prediction that higher engagement of reactive fear circuits would lead to helping, we performed a decision type (help, not help) by threat imminence (distal, imminent) ANOVA, and focused on brain regions displaying a significant interaction between the two (which would indicate activation differences when making decisions under distal and imminent threat). Of note, due to the reduced number of not helping trials, ‘no help’ decisions in this analysis included not only threat trials in which participants did not help, but also responses made in safe trials (details in Materials and methods). Results showed a significant interaction in the midbrain PAG, bilateral insula, right hippocampus, dorsal ACC, OFC, and vmPFC, which was driven by the distal condition (Table 2; Figure 4E; Appendix 1—table 2). Indeed, during distal threats, higher activation in the hippocampus, vmPFC, and OFC was followed by decisions not to help, whereas higher activation in the dACC, PAG, and insula led to helping decisions.

Table 2
Results of the univariate ANOVAS (FWE<0.05).
Univariate threat level*imminence interaction
R/Lkx, y, zFBA
InsulaR76838, 20, 628.1413
InsulaL521–30, 24, –420.5813
IFG, OFCL593–40, 32, –1623.7211
OFCR4078, 32, –2023.4311
ACCR2784, 32, 2018.2924
OFC, IFGR5234, 40, –1213.1811, 47
Type of decision*imminence interaction
R/Lkx, y, zFBA
Midbrain602, –32, –416.85
InsulaR19238, –16, –220.6213, 47, 22, 44, 6, 45, 21
HippocampusR3830, –14, –2015.03
InsulaL595−36, –8, –418.6813, 22, 44, 6, 47, 45
Dorsal anterior cingulateL151–2, 14, 3016.1324, 6, 32, 5, 4
InsulaR74830, 20, –818.1847, 13, 22, 44, 6, 45, 21
Inferior frontal/orbital gyrusL46–34, 36, –1011.4711
Rectus, ventral medial frontal gyrusL2380, 46, –2019.9411, 10, 25

One limitation of this analysis is that it conflated responses made during safe trials and no help responses. It could be argued that these represent fundamentally different types of decisions. To overcome this limitation, we also performed a parametric modulator analysis, in which we used a GLM that included a parametric modulator for decisions made in shock trials (0 if no help and 1 if helped). This analysis allowed us to model all threat imminence and level conditions, but only assign parametric modulators to those in which a subsequent help or no help decision was made (i.e., distal 1 shock, distal 2 shocks, imminent 1 shock, and imminent 2 shocks). The disadvantage of this approach is that only participants with at least 1 ‘no help trial’ trial per condition were included (N=28). Results of the parametric modulator analysis were consistent with those of the ANOVA (Table 3), in that we only found significant modulation of brain activation by subsequent decision during distal threats. Activation of bilateral insula was increased before helping decisions (distal 2 shocks), and activation of the vmPFC was increased before not helping decisions (distal 1 shock and distal 2 shocks).

Table 3
Results of the parametric modulation analysis.
Distal 1 shock
R/Lkx, y, zTBA
(Neg) vmPFC, rectusL125–4, 30, –226.5511
Distal 2 shocks
R/Lkx, y, zTBA
(Pos) InsulaR11044, 2, 45.7613
(Pos) InsulaL150–38, 0, 85.5613
(Neg) vmPFC, medial frontal orbL64–4, 54, –125.1311

Following reviewer suggestions, we also estimated a first-level model that separated help, no help and safe decisions as a function of threat imminence and level, using a Bayesian approach (details in Materials and methods). The Bayesian first-level analysis was followed by model comparison (using random effects) on the target regions of interest (ROIs) that revealed significant results across the two frequentist analyses (i.e., insula and vmPFC). Using a probability threshold of 0.75 (corresponding to a BF of around 8) (Rosa et al., 2010), in the bilateral insula results indicated stronger evidence for the distal ‘help’ models (distal 1 shock and distal 2 shocks). In the bilateral vmPFC, results indicated stronger evidence across ‘no help’ models (distal 1 shock, distal 2 shocks, imminent 1 shock, and imminent 2 shocks) (Figure 5).

In the left and middle panels, comparison of ANOVA and parametric modulation activation maps.

Since results of these two analyses suggested the effects were driven by the distal condition, here we selected the distal help versus no help (no help+safe) contrast from the ANOVA (left) and the help and no help maps from the parametric modulator regressors. Commonalities were found in the insula (activation associated with helping decisions) and in the vmPFC (associated with not helping decisions). Red denotes higher activation during help decisions, and blue denotes higher activation during not help decisions. In the right panel, results of the Bayesian model selection (BMS; following Bayesian first-level analysis). Resulting model evidence maps were thresholded at 0.75 (BF of approx. 8). ROI masks were then applied for model comparison. Results showed stronger evidence for help models for the insula, and no help models for the vmPFC, in line with the frequentist analyses. Note that results from the parametric modulation and Bayesian analysis are inherently noisier, given the smaller number of participants and trials. ROI, region of interest; vmPFC, ventromedial prefrontal cortex.

We additionally ran a Searchlight analysis to localize dissociable neural patterns guiding decisions under distal and imminent threats. Results showed that, prior to helping decisions, patterns of activation in the insula, hippocampus and dorsal cingulate were distinguishable between distal and imminent threats (FWE<0.05, k>10; Figure 4C).

Neural representations of threat promoted helping

One of our goals was to determine whether helping decisions were predominantly driven by the response to another person’s distress, by one’s own defensive state, or by both. To this end, we performed an ROI-based RSA (Diedrichsen and Kriegeskorte, 2017). This analysis was done separately for imminent and distal trials, and comprised three steps (see Figure 5 and Materials and methods for details). First, we computed neural representational dissimilarity matrices (RDMs) that reflected trial-by-trial variation in activation patterns throughout the scanning task. Second, we used post-scan ratings to construct behavioral RDMs that reflected, respectively, between-trial differences in perceived distress experienced by the co-participant, and between-trial differences in how threatened the participant felt themselves during the scan. Finally, in the third step, we estimated the second-order similarity between neural and behavioral RDMs. This similarity metric allowed us to assess, for each ROI, whether trial-by-trial multivoxel patterns during the scan primarily represented the co-participant’s distress, or the threat to oneself. Importantly, it allowed us to determine whether and to what degree neural representations of other’s distress and of threat to the self were associated with helping behavior. Based on the second-order similarity between neural and behavioral RDMs, we found no evidence for any of the ROIs that neural activity predominantly represented other’s distress or threat to self (Appendix 1—table 5). However, results showed that, regardless of threat imminence, the similarity between neural and threat RDMs in the left amygdala (β=4.41, se=1.35, t=3.27, p=0.006) and left insula (β=2.46, se=0.97, t=2.53, p=0.047) was positively associated with helping behavior (Figure 6; Appendix 1—table 3). In other words, the more strongly these brain regions, especially the amygdala, represented the threat to oneself, the more frequently the participant decided to help.

Regardless of threat imminence, the similarity between neural and threat RDMs in the left amygdala and insula predicted higher frequency of helping decisions.

*p=0.047; **p=0.006. RDM, representational dissimilarity matrix.

Discussion

Our overarching goal with this study was to determine how one’s own defensive responses influence decisions to help others under threat. Our findings strongly suggested that neural circuits that coordinate fast avoidance responses from self-directed threats (reactive fear circuits) (Mobbs et al., 2020) also underlie the protection of others in dangerous situations. More, the extent to which key defensive regions (amygdala and insula) represent the threat to the self (and not others’ distress) predicts more frequent helping decisions.

Reactive fear circuits promote both self- and -other defence

Previous work has shown that increased imminence of an other-directed threat facilitates helping behavior, suggesting that the same way threat imminence triggers active avoidance (i.e., fight or flight) from self-directed threats, it may promote defensive helping when others are under threat (Vieira et al., 2020). Here, we examined the neural basis of this effect. We found that multivoxel activation within the defensive circuitry (i.e., amygdala, ACC, insula, hippocampus, and OFC) tracked with the level of threat. Also, we found dissociable local patterns of activation prior to helping as a function of threat imminence within the insula, hippocampus, and dorsal ACC. These findings are consistent with previous work highlighting the role of the hippocampus and dACC as key regions optimizing escape decisions within cognitive and reactive fear circuits, respectively (Qi et al., 2018). Importantly, we found that average activation within several regions of the defensive circuitry, namely the in the PAG, insula, hippocampus, dACC, OFC, and vmPFC, differed between trials leading up to helping or not helping decisions. Specifically, consistent with our predictions, activation in the insula, ACC, and PAG was higher before decisions to help, whereas activation in the hippocampus, vmPFC, and OFC was higher prior to decisions not to help. Therefore, although overall we did not find higher frequency of helping decisions under imminent than distal threat, our results suggest that greater engagement of reactive fear circuits facilitated helping behavior. Our failure to replicate prior behavioral effects of threat imminence could have been due to a difference in how participants were instructed (here, we instructed individuals to balance helping and not helping decisions) and/or statistical power (the previous study had a larger sample defined based on power calculations, N=98). The task instructions, in particular, used may have induced specific metacognitive processes that directly affected the behavioral decisions, and are virtually impossible to quantify. In any case, the present results are in line with prior findings in that they suggest that acute defensive states coordinated by reactive fear circuits promote helping under threat (Vieira et al., 2020). Of note, the differential activation of cognitive and reactive fear circuits based on the subsequent decision (help vs. not help) was only found when decisions were made under distal threat (i.e., when responses were prompted in the beginning of the trial). Also, contrary to our predictions, when comparing the average activation during distal and imminent threats independently of the decision, we found that the bilateral insula and ACC were more active during distal relative to imminent threats, whereas the bilateral IFG and OFC were more active during imminent threats. These results are opposite to those from prior studies that manipulated the imminence of self-directed threats (Mobbs et al., 2009; Wendt et al., 2017; Mobbs et al., 2007; Meyer et al., 2019). It is possible, however, that the disparity between our and previous findings is due to methodological reasons and does not reflect true differences in the processing of self- versus other-directed threats. For instance, in previous work only avoidance responses were made, whereas in our paradigm both help and no help decisions were possible (perhaps more analogous to an approach and avoidance option). Importantly, in our paradigm, every trial started with a distal threat (the circle was static on the left) that always evolved into an imminent threat (the circle would move to the right, and the shock would be administered at the end of the trial). This may have made the threat overly predictable (especially given the high number of trials required for fMRI), engaging reactive fear circuits to a greater extent during the distal phase of the trial.

Nonetheless, our results obtained in a paradigm wherein the threat was directed at another person (the co-participant) were overall consistent with previous fMRI research investigating neural responses to self-directed threats as a function of imminence, suggesting a parallel between responses to self- and other-directed threats. This self-other parallel has been demonstrated in other threat-related processes, such as learning (Olsson et al., 2020). Importantly, in line with previous demonstrations that acute defensive states may promote prosocial outcomes (Vieira et al., 2020), we found that greater engagement of reactive fear circuits may facilitate helping of others in a threatening situation.

The neural representation of threat to the self predicts helping

To decide whether to help another person in a dangerous situation, one must consider not only their distress, but also the threat in the environment. Here, we determined how the representations of another’s distress and of threat guide behavior. To do so, we collected ratings of both the co-participant’s distress, and of how threatened the participant felt during the scan. Crucially, these ratings were obtained after the scan, allowing us to obtain a behavioral metric of how participants independently represented another person’s distress and the threat value of the situation. Although collecting these ratings during the scan would have provided a more direct measure of the neural representation of other’s distress and threat, it would have also have increased scanning time considerably and potentially compromised data quality as a result of fatigue. Performing these ratings after the scan was consistent with previous approaches (Parkinson et al., 2014) and enabled us to avoid explicitly priming participants to consider those distress and threat cues during the scan, which could have influenced their behavior and neural responses. Using behavioral representations of distress and threat, we assessed the extent to which each brain region in the defensive circuitry represented those cues, and its association with helping decisions. We found no indication that any of our ROIs predominantly represented the other’s distress or the threat to own self. However, our results showed that, regardless of threat imminence, the more the left amygdala and left insula represented the threat to oneself, the more participants decided to help.

The association between representation of threat and helping was particularly strong in the left amygdala. The amygdala has long been known to have a pivotal role in the acquisition and expression of defensive responses in mammals (Tovote et al., 2015). For instance, in both humans and rodents, it has been shown to coordinate switches between defensive states across through the communication between its basolateral nucleus (BLA) and oxytocin(OT)-sensitive neurons in the central amygdala (CeA) (Terburg et al., 2018). Our present results suggest that the amygdala’s role in defensive responding may also be relevant for helping behavior. This is consistent with previous work in animals. In rodents, CeA activation by OT not only enables the transition from freezing to fight-or-flight, but has also been shown to trigger offspring care behaviors in females (Rickenbacher et al., 2017), and to enhance maternal aggression (Bosch, 2013). It has additionally been demonstrated that the administration of benzodiazepines, drugs with a known effect on CeA (Griessner et al., 2021), impairs helping behavior in rats (Ben-Ami Bartal et al., 2016). Taken together with these reports, our findings in humans suggest that amygdala-mediated defensive processes may also enable the provision of care to others, here in the form of helping. Of note, our effects were restricted to the left amygdala, which is consistent with several others demonstrations of hemispheric specialization in amygdala function in emotional and pain processing (Allen et al., 2021; Baas et al., 2004; Blair et al., 2005; Butler et al., 2018; Young and Williams, 2010).

We found no evidence that the representation of other’s distress in any ROI was associated with helping, including in brain regions that have previously linked with empathy for pain and distress states (ACC and insula) (Bernhardt and Singer, 2012; Lockwood, 2016). Relatedly, our behavioral analysis revealed that, contrary to previous research (FeldmanHall et al., 2015; Hein et al., 2016), empathic concern was negatively associated with the percentage of helping decisions. These findings are at odds with prior accounts of altruistic responding, which suggest empathy for distress is the key proximal mechanism driving helping behavior (Batson et al., 1987; de Waal and Preston, 2017). Several explanations may account for lack of evidence of an association between neural representation of distress and helping in our data. One is that the behavioral ratings, which were collected after the scan, were not sensitive enough to reflect variations in perceived distress between clips of the confederate. However, this seems unlikely, given that we found significant effects of imminence and threat level on ratings of distress. Another possibility is that distress is mainly represented in brain regions not included in our analysis. The amygdala, insula and ACC would be prime candidate regions to represent another individual’s distress, in light of previous research on empathy (Fallon et al., 2020; Lamm et al., 2011; Marsh et al., 2014; OConnell et al., 2019), but we did not detect an association between the degree to which these regions represented distress and helping behavior. Other potential regions would be those previously implicated in mentalizing, such as the temporo-parietal junction (TPJ) (Lamm et al., 2011; Patel et al., 2019; Schurz et al., 2014). At an exploratory level, we repeated the RSA analysis within anatomical masks of the left and right TPJ, but did not find evidence that the representation of distress in these regions were related to helping behavior (see Appendix 1—table 4). Additional research at the whole brain level is necessary to further assess the representation of distress in the brain, and its impact on helping behavior under threat. In any case, our present results suggest that, even if perceiving distress/need in others to some extent triggers altruistic motivation (Preston, 2013), the ability to provide help may ultimately rely on the activation of circuitry implicated in self-defense.

In summary, our results point to a parallel between responses to self- and other-directed threats, and suggest that the engagement of reactive fear circuits facilitates helping of others. Importantly, we showed that the extent to which the amygdala represents the threat to self (and not other’s distress) predicts helping decisions. These results challenge the idea that empathy for distress is the only proximal mechanism motivating helping decisions, and that overriding self-defensive responses is necessary to help others under threat. Rather, in dangerous situations, one’s own response to the threatening event may enable defensive helping of others, possibly through the activation of neural mechanisms subserving both individual defense and offspring care in mammals.

Materials and methods

This study has been pre-registered (https://osf.io/yvufn) and any deviations from the pre-registration are justified in the Appendix. Data and code will be made available at the OSF project page (https://osf.io/9cuva).

Participants

Forty-nine healthy volunteers (M=24.29, SD=4.78) participated in the experiment. Participants were recruited via flyers posted on- and off-campus, and local online recruitment systems. All participants were right-handed, had normal or corrected-to-normal vision, and were screened for history of psychiatric or neurological diagnoses, current medication, brain injuries, and substance abuse. Participants provided informed consent prior to the experiment, and were compensated for their participation. This work was approved by the Swedish Ethics Board (Etikprövningsmyndigheten).

fMRI helping under threat task

Request a detailed protocol

In each testing session, a participant and a confederate (henceforth, co-participant) were informed the experiment comprised two parts (only one of those parts involved an MRI scan), which would be randomly assigned to each one by flipping a coin. Participant and co-participant were then accompanied to separate testing rooms (the actual participant was taken to the MRI area) and did not interact again during the experiment (details about testing procedures, post-task questionnaires, and debriefing are available in the Appendix).

In the MRI, participants performed a task modified from previous work (Vieira et al., 2020) wherein they made trial-by-trial decisions about whether or not to help the co-participant avoid aversive electrical shocks to the wrist, at the risk of also being shocked (Figure 1A). Threat imminence was manipulated by varying the spatial position and movement of a visual cue signaling varying levels of threat on a computer screen. Respectively, a green circle signaled no threat (no shocks), a yellow circle signaled moderate threat (1 upcoming shock), and the red circle signaled high threat (2 upcoming shocks). In addition, a webcam feed of the co-participant was presented on the screen throughout the task. Unbeknownst to the participant, the video feed was in fact pre-recorded, and edited to select unique clips for each trial of the task.

Participants were informed that, throughout the experiment, they and the co-participant would see the same screen. Each trial started with a static cue on the left side of the screen (4 s), which then moved to the right (4 s). In shock trials, the co-participant would be administered an aversive shock to the wrist when the cue reached the right end of the screen, unless participants decided to help him. To decide whether they wanted to help the co-participant avoid the upcoming shock, participants made forced-choice responses by pressing 1 (Help) or 2 (Do not help) on an MRI-compatible button box as soon as the response slide was displayed (1.25–1.75 s). Responses were prompted sometimes in the beginning of the trial, when the visual cue was static on the left side (distal threat), and other times at the end of the trial, after the visual cue had moved to an endpoint on the right, and thus immediately before shock delivery (imminent threat). Of note, the time available to make a response was identical in distal and imminent conditions. Although in naturalistic settings higher threat imminence generally coincides with less time to make a decision, here we opted to dissociate imminence from decision for the sake of experimental control. Outcomes of participants’ decisions were as follows: if they chose not to help, the co-participant would always receive a shock; if they chose to help, there would be around 70% probability of both participant and co-participant receiving 1 shock (in moderate threat trials) or 2 shocks (high threat trials). Shocks were administered on the left ankle. Participants were instructed they should respond as quickly as possible. Also, to discourage missed responses, they were informed that a shock would be delivered to both participants (with 100% chance) whenever a response was not detected. Finally, to balance the number of helping and non-helping trials during the scanning session, participants were informed that they would have a pre-set number of times they could help on each run, and thus they should try to balance, per run, the number of times they helped and not helped. In reality, participants could help on as many trials as they wished. These task instructions may have affected overall helping performance as participants were asked to consider the overall balance of helping and not helping choices whenever making a new decision. Nonetheless, participants were still left with the decision whether to help or not on any given trial, allowing us to examine the corresponding brain activation. It is also important to point out that these instructions may have introduced additional metacognitive demands on the task. Yet, given the within-subject design, these demands are not expected to have introduced a systematic bias in the data. Shock administration always happened at the end of the trial, and participants were able to see the outcome of their decisions on the screen (i.e., the co-participant receiving or not receiving a shock; 4 s).

Safe trials followed an identical structure, with response slides presented at distal or imminent stages in relation to the end of the trial. However, participants were instructed that no shocks would be given and they should arbitrarily choose to press 1 or 2 when the response slide was displayed. It was made clear to them that their choice would have no consequences for them or the co-participant.

The task included 144 trials split into 8 functional runs (approx. 8 min). Each run comprised 18 trials, 9 distal and 9 imminent, and 6 of each threat level (resulting in 24 trials per condition, in total). Distal and imminent trials were presented in blocks, and the order of blocks was in each run. Within each block, safe, moderate threat, and high threat trials were randomized. The order of functional runs was randomized across participants. The task was programmed and delivered using E-prime 3.0 (Psychology Software Tools, Inc, https://www.pstnet.com).

Ratings

Request a detailed protocol

After the scan, participants were taken to a different testing room and asked to complete a follow-up task. Here, all video clips showed during the scanning task were presented to participants, in random order. Participants were informed that these had been recorded during the scan, and that their task now was to, for each clip, rate the level of distress, anxiety or concern they perceived in the co-participant, on a 9-point scale. Participants also presented images of the visual cues at distal (left side of the screen) and imminent positions (right side of the screen), and asked to rate on a 9-point scale how threatened they felt during the scan, whenever they saw those images (Figure 1B). Ratings of distress and threat were presented in separate blocks, and the order was randomized.

fMRI acquisition and preprocessing

Request a detailed protocol

Participants were scanned in a single session at the Stockholm University Brain Imaging Center (SUBIC), using a 3T Siemens scanner with a 64-channel head coil. First, a high-resolution T1-weighted anatomical scan was obtained (TR=2300 ms, TE=2.98 ms; FoV=256 mm, flip angle=9°, and 192 axial slices of 1 mm isovoxels), followed by 8 functional runs, of about 8 min each. Functional images were acquired with an echo-planar T2*-weighted imaging sequence with whole-brain coverage while participants performed the fMRI task (TR=1920 ms, TE=30 ms, FoV=192 mm, flip angle=70°, 62 interleaved slices of 2 mm isovoxels, and acceleration factor of 2).

Preprocessing of fMRI data was done using SPM12 (Wellcome Trust Centre for Neuroimaging; https://www.fil.ion.ucl.ac.uk), and included slice timing correction, realignment to the volume acquired immediately before the anatomical scan (i.e., the first image of the first functional sequence) using six-parameter rigid-body transformations (translation M=0.02 mm, min=−0.4, max=0.6; rotation M=0.0006 mm, min=0.02, max=0.008), coregistration with the structural data, normalization to standard space using the Montreal Neurological Institute (MNI) template with a voxel size of 2×2×2 mm3, and smoothing using a Gaussian kernel with an isotropic full-width-half-maximum of 4 mm (Gardumi et al., 2016; Hendriks et al., 2017). Finally, a high-pass filter with a cutoff of 128 s was applied to remove slow signal drifts.

Statistical analysis

Behavioral data

Request a detailed protocol

Behavioral data was in general analyzed using GLMMs, an approach that accounts for variation in the dependent variable that is explained by random sampling of, for instance, participant or trial number (random effects), in addition to the independent variables (fixed effects). Mixed-effects approaches have further been proposed to increase the generalizability of research findings to other individuals and stimuli (Yarkoni, 2020).

Our main behavioral variable was helping behavior, which was operationalized as the percentage of helping responses throughout the task. We modeled helping percentage using GLMMs as a function of threat imminence, threat level, and imminence by threat interaction (fixed effects). The subject was added as a random effect, with random intercept and slope per threat imminence and level. In a separate model, we also added a threat imminence*level*empathic concern interaction, following previous indications that threat imminence may affect helping behavior more strongly in individuals with higher caregiving tendencies. To account for the possibility that behavior varied throughout the experiment, we also performed a mixed effects logistic regression on single trial dichotomous responses (help or no help), including the trial number as a fixed effect in addition to the other fixed effects. Finally, following recent recommendations to consider within-individual effect sizes (Grice et al., 2020), we also calculated the difference between number of helping decisions under imminence and distal threat, per individual.

Reaction times were averaged per condition, and analyzed in a GLMM with threat imminence, threat level, and imminence by level interaction as fixed effects, and the subject as a random effect (intercept and slope). Similarly, post-task ratings of other’s distress and threat to self were analyzed in a GLMM with threat imminence, threat level, and imminence by level interaction as fixed effects, and the subject as a random effect (intercept and slope).

Imaging data

Request a detailed protocol
First-level analysis
Request a detailed protocol

First-level analysis was performed in SPM12 and was based on the general linear model. Time-series of each voxel were normalized by dividing the signal intensity of a given voxel at each point by the mean signal intensity of that voxel for each run and multiplying it by 100. Resulting regression coefficients thus represent a percent signal change from the mean. Regressors were created by convolving the train of stimulus events with a canonical hemodynamic response function. Three different GLMs were estimated based on the goal of the analysis. For assessing differences based on threat imminence and level, six regressors of interest were modeled corresponding to the time window of the visual threat cue (distal safe, distal 1 shock, distal 2 shocks, imminent safe, imminent 1 shock, and imminent 2 shocks). These regressors were defined based on the position of the threat cue on the screen (static=distal; approaching=imminent), and not relative to when the participant made a decision. In addition, eight regressors of no interest were added in the model, corresponding to the time window of the response and the outcome, plus the six motion parameters estimated during realignment.

To assess differences in neural response based on the type of decision, another model was estimated with six regressors of interest: help distal 1 shock, help imminent 1 shock, help distal 2 shocks, help imminent 2 shocks, no help distal, and no help imminent. Here, distal and imminent refer specifically to when the decision was prompted in the trial. Of note, due to the reduced number of no helping trials for some participants, the ‘no help’ regressor included both not help decisions, and decisions made during safe trials, wherein no shocks were given. In addition, eight regressors of no interest (response, outcome, and six motion parameters) were added. Because of potential concerns in modelling no help and safe trials together, we additionally created a model that included regressors for distal and imminent events (per threat level) with a subsequent decision, and for each shock regressor (1 shock and 2 shocks) we added a parametric modulator to reflect the response subsequently made in that trial (0=no help, 1=help). The limitation of this approach is that only participants with at least 1 ‘no help trial’ trial per condition were included (N=28). Following reviewer advice, we also implemented a Bayesian approach (using Variational Bayes in SPM12) to obtain a first-level model that produced separate estimates for help, no help and safe trials, based on threat imminence and level. We then performed model comparison on ROI evidence maps thresholded at p>0.75 (following Rosa et al., 2010).

Finally, a fourth model was created to enable subsequent trial-by-trial RSA, wherein one regressor was estimated for each individual trial, modelling the time window of the threat cue.

Regions of interest
Request a detailed protocol

Given our focus on defensive brain circuitry, our analysis targeted pre-specified ROIs that were anatomically defined, including the left and right amygdala, left and right hippocampus, left and right insula, midbrain, left and right ACC, left and right vmPFC, and left and right vlPFC (Figure 3). ROIs were defined on the Wake Forest University (WFU) Pickatlas toolbox (http://www.fmri.wfubmc.edu/cms/software; Maldjian et al., 2003).

Multivoxel pattern analysis and support vector regression
Request a detailed protocol

Beta values derived from first-level analyses were used in multivariate analyses, including multivoxel pattern analysis (MVPA) searchlight, support vector regression (SVR), and RSA. Spatially distributed patterns of activation across voxels can reveal distinguishable neural responses between experimental conditions even in the absence of significant average activation differences in single voxels, making multivariate approaches more sensitive than conventional univariate analysis (Formisano et al., 2008).

To identify local activation patterns that distinguish between distal and imminent threats, we used an MVPA Searchlight implemented in The Decoding Toolbox (TDT). A spherical searchlight (radius 12 mm) was moved throughout each participant’s data and, at each searchlight center, a support vector machine algorithm was trained to discriminate activation patterns in response to distal and imminent threats. Training was done iteratively on each 7 functional runs and tested on the 8th (leave-one-out cross-validation). Resulting percentage score at each voxel for participant was calculated and displayed in individual accuracy maps. Accuracy maps were then analyzed at the group level in a one-sample t test implemented in SPM12. A similar approach was taken to the identification of local activation patterns that discriminated helping decisions during distal and imminent threats. Additionally, an SVR was used to identify local patterns that showed a continuous linear association with threat level (safe, 1 shock, and 2 shocks). As for the classification searchlight, individual accuracy maps were analyzed at the group level in a one-sample t test. Group results were thresholded at voxelwise FWE<0.05 and only clusters with more than 10 voxels were further considered.

Univariate analysis
Request a detailed protocol

We also performed GLM-based univariate analyses. We analyzed BOLD signal to the threat cues, regardless of decision, in a threat imminence (distal and imminent) by threat level (safe, 1 shock, and 2 shocks) repeated-measures ANOVA. We also analyzed BOLD signal to threat cues prior to the decision in a decision type (help and not help) by threat imminence (distal and imminent) ANOVA. Full results for these analyses are available in . Finally, we performed a One-sample t test on parametric modulator contrast images to identify brain regions wherein average activation during the threat was modulated by subsequent decision. Univariate analysis results were first thresholded at p<0.001 uncorrected. With this threshold, clusters with more than 10 voxels were significant with FWE-corrected p<0.05.

Representational similarity analysis
Request a detailed protocol

One of the goals of the study was to characterize representations of other’s distress and of threat to oneself by defensive regions, and determine its relation to helping behavior. RSA was done separately for imminent and distal trials, and comprised three steps (Figure 7). On the first step, we modeled each trial in the first-level analysis (in SPM12), in order to estimate one beta coefficient per trial. Then, for each ROI, we extracted betas from each voxel in each trial to estimate the correlation of beta-values (expressed in r values) between all voxels per trial. Resulting r values were used to construct a representational similarity matrix across all trials that reflects the correlation between all voxels in each trial. This matrix was then transformed (1−r) to reflect dissimilarity instead of similarity (RDM). On the second step, distress ratings provided after the scan on the unique video clips shown in each trial were used to construct a dissimilarity matrix that reflects the difference in perceived distress of the co-participant between trials (expressed in Euclidean distances). Post-scan threat level ratings were used in an identical manner to construct a dissimilarity matrix that reflects between-trial differences in how threatened the participant felt during the scan. Finally, on the third step, we estimated the second-order similarity (kendall’s τ) between neural and behavioral RDMs. In a nutshell, this similarity metric allowed us to assess, for each ROI, whether trial-by-trial multivoxel patterns during the scan represented the co-participant’s distress and the threat to oneself. Importantly, it allowed us to determine whether neural representations of other’s distress and of threat to oneself were associated with helping behavior. To do so, second-order similarity values were entered in a linear model predicting average helping percentage during the scan. Predictors in this model were the similarity between neural and threat RDM, the similarity between neural and distress RDM, as well as threat imminence. Thirteen linear models were estimated, one for each ROI (i.e., left and right amygdala, left and right hippocampus, left and right insula, midbrain, left and right ACC, left and right vmPFC, and left and right vlPFC). False discovery rate (FDR) correction was applied to adjust the p value of all coefficient estimates, across all 13 models. FDR-corrected p values below α=0.05 were considered significant. Beta value extraction was performed in Matlab, and all remaining steps and analyses of the RSA were performed through custom-made scripts in R (code available at https://osf.io/nb6cf/).

Schematic of the RSA pipeline.

On step 1, we extracted the vector of trial-by-trial betas for each voxel in a given ROI. We then calculated the correlation (Pearson r) between all trial pairs. These correlation values were inverted (1−r) and used to create a trial-by-trial matrix, wherein each cell represents how correlated activation across all voxels of the ROI was in each trial pair (neural representational dissimilarity matrix, RDM). On step 2, post-scan ratings of the co-participant’s distress in each unique clip were used to construct a trial-by-trial matrix, wherein each cell contained the Euclidean distance between the rating of each pair of clips (distress RDM). A similar method was used with the ratings of threat to the participant (threat RDM). On step 3, the second-order similarity between the neural RDM and distress RDM, and between the neural RDM and threat RDM were calculated using a ranked correlation method (Kendall’s tau). ROI, region of interest; RSA, representational similarity analysis.

Appendix 1

Experimental procedures

The experiment involved one session. Upon arrival, both participant and co-participant (i.e., confederate) were welcomed into the lab and told that the experiment comprised two parts and each one of them would only perform one. The parts were assigned through a fake coin flip to give the impression of randomization. Participant and co-participant were then accompanied to separate rooms to receive more detailed instructions.

Prior to entering the MRI, written informed consent was obtained, and electrodes for administration of electrical shocks were placed on participants’ left ankle. Participants were then given a short practice block on a computer, to allow them to familiarize with the timing of responses. After the practice, participants were taken into the MRI scanner room.

Once in the scanner, shock intensity was individually calibrated using a standard work-up procedure. Participants were asked to select an intensity level that was “not painful, but very uncomfortable” and that it should be something that “if they could avoid, they would rather avoid”. It was emphasized to participants that, during the experiment, they would be able to avoid the shock if they so desired.

Participants were then given written instructions for the first block of the task. They were informed they would see the co-participant via webcam and they would both be presented the same stimuli on the screen. The first block consisted of 3 trials (safe, 1 shock and 2 shocks), wherein the participant was simply required to pay attention to the screen; they would not be asked to make any responses and would not be given any shocks. The goal of this block was to allow participants to see the consequences of safe and threat trials to the co-participant. This was meant to discourage them from making “test responses” in the first trials of the actual task (for instance, decide not to help just to see whether the co-participant would indeed receive a shock).

Thereafter, participants were given written instructions for the rest of the task. They were explained that, on each trial, they would be asked whether they wanted to help the co-participant avoid the shock(s) or not, and what would be the outcomes of their decisions. Participants were also informed that: 1. the co-participant was not aware that shock administration was decided by them, 2. they should try to balance out the number of help and not help decisions, since there was a preset number of times they could help in each run, 3. they would not swap places with the co-participant afterwards, 4. they would not meet the co-participant again after finishing the experiment, and 5. their behavior during the task would not be observed nor filmed.

After the scan, participants were taken to another testing room to perform the ratings task. At the end, they were asked to fill out questionnaires, including post-tasks questions designed to assess whether the instructions and cover-story were believable (e.g., To what extent do you think the shocks were controlled by you? How authentic did the co-participant seem to you?). Lastly, they were fully debriefed.

Post-task and questionnaire measures

After the fMRI helping under threat task and the ratings task, participants completed a series of questionnaires to assess individual differences in empathy and threat sensitivity, including the Interpersonal Reactivity Inventory (Davis, 1983), the Trait Fear Questionnaire (Kramer et al., 2020), the, and the Triarchic Psychopathy Measure (Patrick et al., 2009). Additionally, participants completed a series of post-task questions designed to assess the believability of the experiment.

Univariate analysis

Appendix 1—table 1
Threat imminence X Threat level Anova.
Main effect of threat imminence
R/Lkx, y, zFBA
InsulaR3642, 2,–612.7613
InsulaL113–28, 12,–2016.2947
IFGR141050, 20,–1025.4347
vmPFCL334–8, 32,–1617.4311, 25, 10
OFCL37–18, 56,–812.6211
Main effect of threat level
R/Lkx, y, zFBA
IFGR7236, 24, 820.2513, 45, 47, 44,
vmPFCL4430, 34,–1627.3510, 11, 32, 25, 9
Threat imminence*level
R/Lkx, y, zFBA
InsulaR76838, 20, 628.1413
InsulaL521–30, 24,–420.5813
IFG, OFCL593–40, 32,–1623.7211
OFCR4078, 32,–2023.4311
ACCR2784, 32, 2018.2924
OFC, IFGR5234, 40,–1213.1811, 47
Appendix 1—table 2
Type of decision X Threat imminence Anova.
Main effect of type of decision
R/Lkx, y, zFBA
HippocampusL128−32,–18, –1417.85
HippocampusR5828,–16, –1616.01
InsulaL638–40, 8, 624.1013, 47, 44, 45, 22
InsulaR94946, 14,–230.4313, 47, 44, 45, 22
ACCR2182, 22, 2824.3332, 6, 24, 8
IFGL104–38, 28,–1412.3821, 38, 47, 22
vmPFCL7680, 44,–1824.7411, 25, 10, 32
OFCR12936, 54,–1414.5011, 10
Main effect of threat imminence
R/Lkx, y, zFBA
InsulaR4036,–20, 821.0713
5638, 8, 417.9313
4532, 20, 1417.0913, 45
IFGL211–50, 42,–829.3747, 10
vmPFCL447–2, 54,–1436.3911
Type of decision*threat imminence
R/Lkx, y, zFBA
Midbrain602,–32, –416.85
InsulaR19238,–16, –220.6213, 47, 22, 44, 6, 45, 21
HippocampusR3830,–14, –2015.03
InsulaL595−36,–8, –418.6813, 22, 44, 6, 47, 45
Dorsal anterior cingulateL151–2, 14, 3016.1324, 6, 32, 5, 4
InsulaR74830, 20,–818.1847, 13, 22, 44, 6, 45, 21
IFG/orbital gyrusL46–34, 36,–1011.4711
Ventral med frontal gyrusL2380, 46,–2019.9411, 10, 25

RSA

Appendix 1—table 3
Estimates from models for each ROI, predicting helping percentage throughout the scan as a function of neural-distress similarity, neural-threat similarity, and threat imminence (all p values across all models were FDR-corrected).
EstimateStd. Errort valuep (FDR-corrected)
L amygdala
Intercept0.5360.0549.823<.00001
Threat4.4071.3483.269.006*
Distress–1.2922.049–0.630.887
Imminence0.0400.0550.731.837
L insula
Intercept0.5310.0598.883<.00001
Threat2.4610.9742.525.047*
Distress1.2861.5520.828.837
Imminence0.01290.0520.246.887
L ACC
Intercept0.6490.05711.257<.00001
Threat0.2160.9400.230.887
Distress–1.3391.304–1.026.742
Imminence0.0180.0540.348.887
L hippocampus
Intercept0.6200.0679.233<.00001
Threat0.0421.5640.027.978
Distress0.1041.6420.063.968
Imminence0.0170.0540.324.887
Midbrain
Intercept0.7220.06610.936<.00001
Threat–2.1961.694–1.296.573
Distress–2.6221.753–1.495.451
Imminence0.0400.0550.736.837
L vmPFC
Intercept0.5890.0599.842<.00001
Threat0.7521.0100.744.837
Distress0.6271.4950.419.887
Imminence0.0170.0540.316.887
L vlPFC
Intercept0.5890.0599.842<.00001
Threat0.7521.0100.744.837
Distress0.6271.4950.419.887
Imminence0.0170.0540.316.887
R amygdala
Intercept0.6950.06211.170<.00001
Threat–1.6511.584–1.042.742
Distress–1.3191.411–0.935.797
Imminence0.0170.0530.332.887
R insula
Intercept0.6030.06010.015<.00001
Threat0.4360.9760.447.887
Distress0.1891.6760.113.946
Imminence0.0200.0540.376.887
R ACC
Intercept0.5410.0608.883<.00001
Threat1.7811.2791.392.512
Distress1.6171.4561.110.739
Imminence0.0140.0530.276.887
R hippocampus
Intercept0.6770.06310.677<.00001
Threat–1.4511.432–1.012.742
Distress–0.7041.350–0.521.887
Imminence0.0170.0540.326.887
R vmPFC
Intercept0.6120.05810.490<.00001
Threat0.9331.0550.884.820
Distress–0.4671.345–0.347.887
Imminence0.0120.0470.272.887
R vlPFC
Intercept0.6030.0669.134<.00001
Threat–0.2421.405–0.172.916
Distress1.1181.6700.669.876
Imminence0.02070.0540.381.887
Appendix 1—table 4
Exploratory examination of association between neural-distress and neural-threat similarity, and helping in the temporo-parietal junction (TPJ) (anatomically defined ROI based on aal, including supramarginal and angula gyri).

Estimates from linear models for left and right TPJ.

EstimateStd. Errort valuep
L TPJ
Intercept0.5520.05110.756<.00001
Threat1.6690.9671.7250.088
Distress1.1031.4920.7390.462
Imminence0.0140.0530.2690.789
R TPJ
Intercept0.5590.0579.859<.00001
Threat1.4070.8421.6710.098
Distress0.5851.4910.3920.696
Imminence0.0180.0540.3390.735
Appendix 1—figure 1
Example RDMs (Neural, Distress and Threat) sorted by threat level.
Appendix 1—table 5
Comparison of Neural-Distress similarity, and Neural-Threat similarity, for distal and imminent threats, in each ROI.
DistalImminent
L Accz=–0.039, P=0.969z=–0.028, P=0.978
L Amyz=–0.047, P=0.963z=0.014, P=0.989
L Insulaz=–0.003, P=0.997z=0.0064 p=0.995
L Vlpfcz=–0.073, P=0.941z=–0.036, P=0.971
L Vmpfcz=–0.029, P=0.977z=–0.036, P=0.971
L Hippocampusz=0.001, P=0.999z=0.002, P=0.998
Midbrainz=–0.004, P=0.997z=0.026, P=0.979
R Accz=–0.030, P=0.976z=–0.005, P=0.996
R Amyz=0.004, P=0.997z=0.004, P=0.997
R Insulaz=–0.061, P=0.951z=–0.039, P=0.969
R Vlpfcz=–0.005, P=0.996z=–0.027, P=0.978
R Vmpfcz=0.002, P=0.998z=–0.028, P=0.978
R Hippocampusz=–0.010, P=0.992z=0.010, P=0.992

Response distributions per fMRI run

Appendix 1—figure 2
Response distributions per fMRI run.

Cross-validation

Sample size was limited based on resources and no power analysis was performed. The number of trials was selected to guarantee a balance between having a reasonable number of data points per condition, and keeping the length of the scan within limits that were manageable for participants. The number of trials we used (144) was in fact higher than that of previous work that implemented similar RSA approaches with neural and behavioural data (e.g., Parkinson et al., 2014, J Neurosci; Kim et al., 2022, NeuroImage). To maximize generalizability, we have nonetheless followed the suggestion of one of the manuscript reviewers and implemented a cross-validation procedure for the RSA analysis.

Our RSA analysis comprised several steps:

  1. Computing the initial voxel by trial matrix (which contains beta values for each voxel in each trial of the experiment).

  2. Calculating the neural RDM (a trial x trial dissimilarity matrix in which each cell represents the correlation across all voxels of the ROI between pairs of trials).

  3. Calculating second-order correlations between neural RDMs and behavioral RDMs, in which assessed the association between the neural activation pattern in each ROI and two behavioral models (provided by the distress and threat ratings),

  4. And finally, predicting helping behavior in a linear model based on two predictors: the correlation between neural and distress RDM and the correlation between neural and threat RDM.

    The main goal of cross-validation is to guide the selection of the best fitting model. We thus used cross-validation on point 4 above, since in the previous points no model estimation or comparison was performed. For each ROI, we ran a 10-fold cross validation (between 5 and 10 folds are typically recommended) on the linear model predicting percentage of helping responses as a function of neural similarity with distress ratings, neural similarity with threat ratings, and threat imminence.

Resulting R squares per ROI are shown in the graph below, followed by a table with the full results (root mean squared error, RMSE, mean absolute error, MAE, and corresponding SDs):

Appendix 1—figure 3
Cross-validation.

K-fold cross validation results

ROIRsquaredRsquaredSDRMSERMSESDMAEMAESD
L Amy0.1630.2000.2370.0470.1780.047
L ins0.2040.2060.2450.0360.1870.027
L ACC0.0770.0830.2460.0580.1880.046
Midbrain0.1070.0900.2460.0640.1950.046
L Hipp0.2850.2000.2530.0620.1920.042
L vmpfc0.1660.1950.2460.0560.1900.039
L vlpfc0.0280.0350.2520.0460.1910.031
R Amy0.0610.0530.2470.0550.1860.040
R ins0.2150.1640.2540.0510.1930.034
R ACC0.1330.1630.2480.0420.1870.034
R Hipp0.0980.1140.2500.0670.1910.047
R vmpfc0.0470.0620.2510.0360.1890.034
R vlpfc0.1250.1270.2490.0620.1930.050

It should be noted that cross validation is especially valuable to guide subsequent model comparison and selection. This was not the goal of our linear model, which aimed at performing statistical inference on the predictive power of each of the predictors. Therefore, in the present case, the interpretability of the cross-validation procedure is limited.

Deviations from pre-registration

A pre-registration for this study can be found at. Some changes were implemented in relation to the pre-registration, namely:

  1. We restricted our analysis to ROIs for which we had predictions, instead of reporting whole-brain analysis. Our data, including the whole-brain maps, will be available on OSF.

  2. Correction for multiple comparisons in the analysis of fMRI data was done via Family-Wise Error (FWE) and not False-Discovery Rate (FDR). FWE was deemed a more conservative approach, and adequate to our focus on regions-of-interest.

  3. Group results for MVPA searchlight analyses were done with t-test and not permutations. This alteration was in line with recommended procedures (Hebart et al., 2014). Median maps were assessed to confirm they were in agreement with the parametric results.

  4. Out of 49 participants, only one participant stated categorically not believing in the cover story. We deemed it to be a more conservative approach to keep this participant in the analysis, especially since it was unclear at what point of the experiment they started to have doubts.

  5. We pre-registered an exploratory whole-brain RSA Searchlight analysis that is still ongoing, and is thus not included in the present manuscript.

Data availability

This study has been pre-registered (https://osf.io/yvufn). Anonymized data, code, and materials used in the study are available on the OSF project page (osf.io/9cuva).

The following data sets were generated
    1. Vieira JB
    (2020) OSF
    ID 9cuva. Helping under threat - Part 2.

References

    1. Bosch OJ
    (2013) Maternal aggression in rodents: brain oxytocin and vasopressin mediate pup defence
    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 368:20130085.
    https://doi.org/10.1098/rstb.2013.0085
  1. Book
    1. Fanselow MS
    2. Lester L
    (1988)
    A functional behavioristic approach to aversively motivated behavior: predatory imminence as a determinant of the topography of defensive behavior
    In: Fanselow MS, editors. Evolution and Learning. Psychology Press. pp. 185–212.
  2. Book
    1. Oliner SP
    (2004)
    Do Unto Others: Extraordinary Acts Of Ordinary People
    Basic Books.
    1. Roelofs K
    (2017) Freeze for action: neurobiological mechanisms in animal and human freezing
    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 372:20160206.
    https://doi.org/10.1098/rstb.2016.0206

Article and author information

Author details

  1. Joana B Vieira

    1. Department of Psychology, Faculty of Health and Life Sciences, University of Exeter, Exeter, United Kingdom
    2. Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden
    Contribution
    Conceptualization, Data curation, Formal analysis, Supervision, Investigation, Methodology, Writing – original draft, Writing – review and editing
    For correspondence
    j.m.b.vieira@exeter.ac.uk
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7335-4588
  2. Andreas Olsson

    Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden
    Contribution
    Funding acquisition, Writing – review and editing
    Competing interests
    No competing interests declared

Funding

Vetenskapsrådet (2018-00877)

  • Andreas Olsson

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors thank Dean Mobbs, India Morrison, Björn Lindström, Artin Arshamian, and Henrik Ehrsson for their helpful and insightful comments on our manuscript, and Rita Almeida for support with analysis. This study was supported by a Consolidator Grant (2018-00877) from the Swedish Research Council (Vetenskapsrådet) to AO.

Ethics

This work was approved by the Swedish Ethics Board (ETIKPRÖVNINGSMYNDIGHETEN) in Sweden (Dnr 2018/2200-31 (original); Dnr 2019-04335 (amendment)). All participants provided written informed consent to participate in the experiment, and for the data to be published.

Version history

  1. Preprint posted: December 18, 2021 (view preprint)
  2. Received: February 25, 2022
  3. Accepted: September 29, 2022
  4. Version of Record published: October 25, 2022 (version 1)

Copyright

© 2022, Vieira and Olsson

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,745
    Page views
  • 264
    Downloads
  • 1
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Joana B Vieira
  2. Andreas Olsson
(2022)
Neural defensive circuits underlie helping under threat in humans
eLife 11:e78162.
https://doi.org/10.7554/eLife.78162

Share this article

https://doi.org/10.7554/eLife.78162

Further reading

    1. Cell Biology
    2. Evolutionary Biology
    Jonathan E Phillips, Duojia Pan
    Research Advance

    The genomes of close unicellular relatives of animals encode orthologs of many genes that regulate animal development. However, little is known about the function of such genes in unicellular organisms or the evolutionary process by which these genes came to function in multicellular development. The Hippo pathway, which regulates cell proliferation and tissue size in animals, is present in some of the closest unicellular relatives of animals, including the amoeboid organism Capsaspora owczarzaki. We previously showed that the Capsaspora ortholog of the Hippo pathway nuclear effector Yorkie/YAP/TAZ (coYki) regulates actin dynamics and the three-dimensional morphology of Capsaspora cell aggregates, but is dispensable for cell proliferation control (Phillips et al., 2022). However, the function of upstream Hippo pathway components, and whether and how they regulate coYki in Capsaspora, remained unknown. Here, we analyze the function of the upstream Hippo pathway kinases coHpo and coWts in Capsaspora by generating mutant lines for each gene. Loss of either kinase results in increased nuclear localization of coYki, indicating an ancient, premetazoan origin of this Hippo pathway regulatory mechanism. Strikingly, we find that loss of either kinase causes a contractile cell behavior and increased density of cell packing within Capsaspora aggregates. We further show that this increased cell density is not due to differences in proliferation, but rather actomyosin-dependent changes in the multicellular architecture of aggregates. Given its well-established role in cell density-regulated proliferation in animals, the increased density of cell packing in coHpo and coWts mutants suggests a shared and possibly ancient and conserved function of the Hippo pathway in cell density control. Together, these results implicate cytoskeletal regulation but not proliferation as an ancestral function of the Hippo pathway kinase cascade and uncover a novel role for Hippo signaling in regulating cell density in a proliferation-independent manner.

    1. Evolutionary Biology
    2. Immunology and Inflammation
    Zachary Paul Billman, Stephen Bela Kovacs ... Edward A Miao
    Research Article

    Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles, and birds gave rise to GSDMA–D in mammals. Uniquely in our tree, amphibian, reptile, and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.