- Research article
- Open Access
Emotion reactivity-related brain network analysis in generalized anxiety disorder: a task fMRI study
BMC Psychiatry volume 20, Article number: 429 (2020)
Generalized anxiety disorder (GAD) is closely associated with emotional dysregulation. Patients with GAD tend to overreact to emotional stimuli and are impaired in emotional regulation. Using emotional regulation task, studies have found hypo-activation in prefrontal cortex (PFC) of GAD patients and concluded with inadequate top-down control. However, results remain inconsistent concerning PFC and limbic area’s reactivity to emotional stimuli. What’s more, only a few studies aim to identify how limbic area interacts with PFC in GAD patients. The current study aims to identify the difference in PFC-limbic circuitry response to emotional stimuli between GAD patients and healthy controls (HCs) from the perspective of brain network. Through brain network analysis, it revealed the connectivity between limbic area and PFC, and moreover, the orientation of connectivity, all of which gave a better test of inadequate top-down control hypothesis.
During fMRI scanning, participants were required to complete an emotional face identification task (fearful, neutral, happy facial expression). 30 participants (16 GAD patients, 14 HCs) were included in the formal analysis. A Bayesian-network based method was used to identify the brain network consisting of several pre-hypothesized regions of interest (ROIs) under each condition (negative, positive, neutral). In total, six graphs were obtained. Each of them represented the brain network that was common to the group under corresponding condition.
Results revealed that GAD patients showed more bottom-up connection but less top-down connection regardless of condition, relative to HCs. Also, the insula was more connected but the amygdala was less connected regardless of condition, relative to HCs. the results also revealed a very different brain network response between GAD patients and HCs even under neutral condition.
More bottom-up connection but less top-down connection may indicate that GAD patients are insufficient in top-down control, in keeping with inadequate top-down control hypothesis. The more connected insula may indicate GAD patients’ abnormality in interoception processing. Relative to HCs, distinct brain network response pattern in GAD patients under neutral condition suggests GAD patients’ abnormality in distinguishing safety from threat and intolerance of uncertainty.
Generalized anxiety disorder (GAD) is thought to be one of the least successfully treated psychiatric disorders, which is largely attributable to its unclear neurobiological basis. Patients with GAD demonstrate pervasive, sustained, uncontrollable worry as their major symptomatic complaint [1, 2] and this is closely associated with emotional dysregulation [3, 4]. Neuroimaging and clinical studies in patients with GAD have contributed to the emotion dysregulation model . The emotion dysregulation model holds that patients with GAD are inclined to overreact to both negative and positive emotional stimuli and exhibit abnormalities in emotion regulation [3, 4].
fMRI studies using emotional task have consistently found that abnormalities in prefrontal cortex (PFC)-limbic area in patients with GAD implicate in GAD’s pathology [5, 6]. According to emotional dysregulation model, patients with GAD tend to overrespond to emotional stimuli, but particularly, negative stimuli (e.g. heightened emotional intensity) [3, 4]. Relevant studies have mainly been conducted by using emotion processing task. Many studies have indicated that compared with healthy controls (HCs), GAD patients tended to over-activate PFC and/or cingulate cortex while confronting angry faces [7, 8], rating fearful faces , processing emotional images  and fearful pictures . Though most of the extant studies using emotional processing task are consistent with this general picture, there are also some inconsistencies. For example, study has found that GAD patients were attenuated in PFC and anterior cingulate cortex (ACC) response to emotional stimuli . Hyper-activation has also been reported in amygdala while GAD patients faced angry faces  and rated fearful faces , but with inconsistency too. For instance, studies have found activation in amygdala reduced  or no difference  while processing fearful faces in GAD patients, compared with HCs. Besides, there are inconclusive results too when it comes to the functional connectivity between amygdala and PFC, which was either decreased , increased  or no difference  in GAD patients, relative to HCs under emotional processing task. Similarly, mixed results have also displayed in functional connectivity between amygdala and cingulate area [10, 11]. In detail, one study has found that functional connectivity between posterior cingulate cortex (PCC) and amygdala was decreased , while another has found increased .
Emotional dysregulation model also predicts that GAD patients are incapable of appropriately regulating their emotions (e.g. maladaptive emotional management) [3, 4]. This has been mostly tested by emotion-regulation related task. Using reappraisal task (in which participants were required to maintain, downregulate or upregulate emotional response to emotional pictures), one study has found that GAD patients displayed reduced activation in PFC (e.g. dorsolateral PFC and dorsomedial PFC) during both reappraisal (reduce emotional response) and maintenance condition . Consistent with this, another study has also found GAD patients displayed reduced activation in PFC (e.g. dACC) in Downregulate trials relative to the View trials . Both studies have concluded that PFC’s hypo-activation found in GAD patients during emotional regulation task demonstrated their insufficient top-down control of limbic area [14, 15]. In line with this conclusion, another two studies have found similar results in GAD patients while required to make emotional adaptation. Etkin et al.  has found that GAD patients failed to activate ACC and a negative top-down ACC-amygdala connectivity under emotional regulation task, and this results was replicated in another study using the same paradigm by Etkin and Schatzberg , in which they have also found reduced ACC activation and decreased connectivity between ACC and amygdala.
Emotional dysregulation includes two separated, but mutually interplayed, abnormalities: abnormal emotional reactivity and regulation of reactivity . Studies (as indicated above) have demonstrated that GAD patients displayed inadequate top-down control while emotion-regulation is required. This is also supported by review studies [5, 19] and meta-analytic study . As to the origin of inadequate top-down control (in other words, PFC’s hypo-activation), it has been suggested that over-responsiveness of limbic circuitry could be the contributor that leads to the fatigue of top-down system (e.g. PFC), which makes the PFC unable to exert top-down control when needed [5, 14]. However, studies using emotional processing task appears to be in favor of overactive top-down control hypothesis, as indicated. What’s more, extant studies have fallen short of providing more direct evidence to this hypothesis. For one thing, results lack of consistency when it comes to PFC and limbic area’s reactivity to emotional stimuli (e.g. amygdala and its connectivity with other regions, as indicated above) ; secondly, very few studies aimed to identify how limbic area interacts with PFC. This makes it hard to interpret whether it’s inadequate/overactive top-down control in GAD. Therefore, for better understanding GAD’s pathology, more effort should be dedicated to clarifying GAD patients’ abnormal reactivity in PFC-limbic circuitry to emotional stimuli.
The current study
This study aimed to identify the abnormal reactivity in the PFC-limbic circuitry of GAD patients under emotional face (happy, neutral, fearful) identification task from the perspective of brain network response. Through brain network analysis, the connection and the orientation of connection between brain regions were determined, which put the relations between PFC and limbic area into straight testing, adding evidence to inadequate/overactive top-down control hypothesis. Noteworthy, identifying the abnormalities in brain network response of GAD patients is important, for 1) it can deepen our understanding of GAD’s pathology: very few studies performed brain network analysis in GAD patients under emotional task, even for those did, most of them used functional connectivity analysis (which is correlation-based) to define brain networks, which normally took two brain regions at a time. When it comes to reveal how brain regions interacts with each other, correlation-based functional connectivity analysis is still facing some major challenges and having a long way to go (for detail, please see reference ). However, by brain network analysis, it not only can unveil how multiple brain regions interact with each other at the same time, but also can make causal inference between brain regions. Which can’t be done by general linear model (GLM) (GLM only addresses the level of response by multiple brain regions ); 2) it can be a powerful tool to differentiate different anxiety disorders: several studies have failed using GLM to identify GAD-specific neurobiological basis in compared with other anxiety disorders (e.g. panic disorder and/or social anxiety disorder) [23, 24]. However, as indicated above, brain network analysis can unveil how brain regions interact with each other, it can unearth potential differential response pattern in brain network between two groups even when GLM results are not significant .
Brain network analysis on fMRI data has, in recent years, gained significant interest, but current methods used on fMRI data have not performed well . Smith et al.  has tested 38 extant methods used to identify networks on fMRI data and found none of them reached statistical significance with valid outcomes. A new algorithm called Multiple sample Greedy Equivalence Search (IMaGES) proposed by Ramsey et al.  has shown promise in identifying brain networks (finding and recalling connections) with excellent accuracy (over 95% in finding and recalling connections) . When IMaGES is combined with an orientation algorithm named Linear non-gaussian Orientation Fixed Structure (LOFS), which can determine the direction of connections, it reaches a precision of over 90% in determining the direction of connections . The combination of IMaGES and LOFS have been widely used in brain network analysis (please see reference [22, 28,29,30,31,32]).
One advantage of using IMaGES and LOFS is that it can include a dozen of brain regions when modeling their interactions without losing its performance. Besides, it needs no priori model to be set up like SEM or DCM . Therefore, with the main focus on PFC-limbic circuitry, we specifically selected PFC, ACC, amygdala, insula and hippocampus as regions of interest (ROIs), upon which we performed brain network analysis. PFC, ACC and amygdala were all reported to act abnormally in GAD, as indicated above. Insula and hippocampus were also included because both insula and hippocampus have been found to be implicated in GAD’s pathology as wellFootnote 1 [1, 33, 34]. What’s more, all of these brain regions are considered as part of emotion circuitry .
We focused on comparing the orientation and density of connections in brain network between GAD patients and HCs. As there are no canonical brain network response pattern that can be used as a reference, we compared brain network characteristics of patients with those of HCs under each condition. Based on emotional dysregulation model and previous findings, we predicted: 1) more bottom-up (connections that originate from limbic area to prefrontal area) but less top-down connections in GAD patients versus HCs; 2) the connection strength between PFC and amygdala would be differential in two groups.
Patients with GAD and HCs were enrolled in this study. Groups were well matched on age and gender. All participants with GAD were consecutively recruited at the Department of Medical Psychology and Department of Mood Disorders of Nanjing Brain Hospital, affiliated with Nanjing Medical University. They were required to meet the following inclusion criteria: (1) A primary diagnosis of GAD by an experienced psychiatrist based on the Diagnostic and Statistical Manual of Mental Disorders (fifth edition, DSM-5TM); (2) Confirmation of GAD diagnosis using the Mini-International Neuropsychiatric Interview (MINI); (3) Free of psychiatric medications for at least 6 months prior to study enrollment; (4) Age 20–60 years old; (5) Right-handed and able to complete all study activities. Exclusion criteria for GAD participants were: (1) Having any neurological disorders affecting the central and/or peripheral nervous systems; (2) Any comorbid psychiatric disorders including depression, panic disorder, bipolar disorder, obsessive-compulsive disorder, schizophrenia, alcohol abuse and/or dependence, social phobia or eating disorder; (3) Severe physical illness, pregnancy and/or breastfeeding; (4) Suicide attempts in the past year; (5) Inability to complete MRI, and (6) A major life change in the last year as defined by death of spouse, unemployment, severe illness, serious injury, legal disputes, property loss, traffic accident, natural disasters or divorce.
HCs were recruited through online advertising and offline posters. Inclusion criteria included: (1) Age 20–60 years old; (2) Right-handed and able to complete all study activities. Exclusion criteria were: (1) Comorbid neurological disorders; (2) History of any symptoms consistent with a psychiatric disorder; (3) Pregnancy and/or breastfeeding; (4) History of psychological consulting or psychotropic medication within 3 months prior to study enrollment; (5) Inability to complete MRI; (6) Major life change in the last year as defined by death of spouse, unemployment, severe illness, serious injury, legal disputes, property loss, accidents, natural disasters or divorce.
Participants were required to identify emotional faces by pressing specific buttons corresponding to certain kinds of emotional faces (button “1” for negative, button “2” for neutral, button “3” for positive) while undergoing fMRI scanning. Emotional faces were all derived from the Chinese Affective Picture System (CAPS) . Specifically, 20 images (10 were female, 10 were male) for each kind of emotion (happy, fearful, neutral) were randomly selected from its corresponding category (e.g. happy for positive, fearful for negative) for a total of 60 images. Of note, the keypad was used to ensure participants’ focus on the task. E-prime was used to program and present the emotional stimuli. A block design was employed for this study with each block consisting of 5 of the same kind of images. Each image was presented for 4 s. Before each block was presented, a 20-s cross was displayed on the screen. The entire task consisted of 12 blocks and lasted for 8 min in total. Images were not repeated. Balanced incomplete Latin square design and full Latin square design were used for within-participants and between-participants, respectively. Before the formal experiment, participants were well informed of the rules and practiced 3 times (images used in practice were not those that used in formal experiment).
fMRI data acquisition
fMRI data was obtained on a 3.0 Tesla Siemens Medical System scanner at Nanjing Brain Hospital. In order to reduce head movement and noise, each participant’s head was secured with foam pads and earplugs were placed in ear canals. During the 8-min (240 volume, 36 slices for each volume) fMRI scan, each participant was instructed to perform the task they had practiced before scanning. Imaging data were acquired using the echo planar imaging (EPI) sequence according to the following scan parameters: Acquisition matrix = 64 × 64, field of view (FOV) = 240 mm × 240 mm, slice thickness = 4 mm, spacing between slices = 4 mm, 36 slices, repetition time (TR) = 2000 ms, echo time (TE) = 30 ms, flip angle (FA) = 90°.
Imaging data processing
Imaging data were preprocessed using Data Processing Assistant Resting-State fMRI (DPARSF) software (http://rfmri.org/DPARSF) implemented in MATLAB. The basic procedures performed on images were slice timing, realignment, normalization to Montreal Neurological Institute (MNI) space with 3 mm × 3 mm × 3 mm resolution and smoothing with a 4 mm full width at half maximum (FWHM) Gaussian kernel. Participants were excluded if head motion exceeded 2.5 mm or 2.5 degrees. Individual and group level analysis were then performed using SPM8 software (http://www.fil.ion.ucl.ac.uk/spm).
Selection of ROIs
ROIs for brain network analysis were pre-hypothesized. We specifically focused on PFC-limbic circuitry, as we indicated previously. ROIs included PFC, ACC, amygdala, insula and hippocampus. All of them were found to implicate in emotional processing (e.g. fear)  and show abnormalities in GAD patients. ROIs were defined based on ALL (Automated Anatomical Labeling) atlas [37,38,39]. To identify clusters of ROIs that were significantly activated under different condition (positive, neutral, negative), one-sample t test was performed on all participants (GAD patients and HCs) under each condition to obtain condition-specific activation map. The significant activated clusters in pre-hypothesized ROIs were used as masks to extract timeseries for brain network analysis.
Brain network analysis
Brain network analysis was achieved with graphical causal modeling using IMaGES and the LOFS algorithm implemented in the TETRAD IV (version 5.3.0–1; http://www.phil.cmu.edu/projects/tetrad) software [27, 40]. Brain network analysis was conducted separately for each group under each condition (HC group under positive, neutral, negative condition; GAD group under positive, neutral, negative condition). Specifically, the time series of each ROI were extracted from each participant under each condition using REST  and were fed to IMaGES in order to identify the brain network shared by multiple participants under each condition. Unlike structural equation modeling (SEM) and dynamic causal modeling (DCM), IMaGES searched over all possible graphs (representatives of brain networks, in which nodes represent brain regions and connections between nodes represent connectivity between them) and through Bayesian Information Criterion (BIC) scoring, IMaGES picked the winning model with no need to posit a parametric model a priori . IMaGES started with an empty graph for a set of nodes and took two stages to find the graph that was shared by multiple participants. Instead of searching over each and every directed acyclic graph (DAG), IMaGES searched over Markov equivalence classes to guarantee efficiency . Since the data for brain network analysis was based upon an indirect measurement of neural activity and easily produced false triangles, IMaGES avoided spurious connections by increasing penalty discount . Penalty discount was used to obtain the first graph with no triangulations . At the forward searching stage, IMaGES added an additional connection (between nodes) at a time whose addition would most improve the BIC score. With each addition of a connection, the graph with the best BIC score would always be chosen . When the addition of a connection no longer improved the BIC score, it initiated a backward searching stage where it eliminated, one at a time, connection whose removal most improved the BIC score. Once no further improvements can be made, it stopped and output a pattern (representative of the brain network) . After IMaGES, the pattern was then submitted to LOFS to determine the connections’ direction for the graph by relying on the assumption that the residuals of the correct model with independent non-gaussian errors was doomed to be less Gaussian than the residuals of any incorrect model [25, 40]. The graph output by LOFS detailed the connections and the direction of the connections in the network which were described with arrows. In total, six graphs representing brain networks were obtained. In our study the Anderson–Darling score  was used to estimate the degree of non-Gaussianity while implementing LOFS. Afterwards, a structural equation modelling (SEM) estimator was employed to estimate the connection strength (represented by SEM coefficients) with a regression optimizer.
In total, 42 participants (23 GAD patients, 19 HCs) were recruited for this study. Of them, two participants (one in the GAD group, one in the HC group) were excluded before data analysis because of poor data quality. For the remaining 40 participants, six GAD patients and four HCs were excluded for not meeting head motion criteria (less than 2.5 mm or 2.5 degrees). For the excluded participants, no significant difference on any characteristic (age: t = 1.5, p = 0.2; gender: χ2 = 0.17, p = 0.68). The remaining 30 participants were well matched on age and gender. See Table 1.
fMRI between group analysis
Between group analysis found that abnormal activations in PFC, ACC, insula, hippocampus in GAD patients, relative to HCs, but no differences were found in amygdala’s activation between two groups (please see supplemental material Table 1).
There were significant activations in each pre-hypothesized ROI (Table 2 and Fig. 1) (corrected with voxel-wise FDR correction p < 0.05). the peak coordinates for clusters were: right superior frontal gyrus (RSFG) (MNI: − 24, 42, 51; t = 5.05), left superior frontal gyrus (LSFG) (MNI: 15, 57, 39; t = 3.78), left anterior cingulate gyrus (LACC) (MNI: − 3, 48, 12; t = 4.59), right insula (RInsula) (MNI: 39, − 15, 18; t = 4.53), left medial frontal gyrus (LMFG) (MNI: − 3, 60, 15; t = 3.45), right parahippocampal gyrus (RParah) (MNI: 33, − 42, − 6; t = 5.99), left parahippocampal gyrus (LParah) (MNI: − 33, − 39, − 9; t = 5.99) and left amygdala (LAmygdala) (MNI: − 18, − 6, − 18; t = 3.91) under negative condition; LACC (MNI: − 3, 54, 0; t = 4.57), LSFG (MNI: − 24, 33, 51; t = 3.93), RInsula (MNI: 39, − 18, 15; t = 5.40), LParah (MNI: − 33, − 42, − 9; t = 5.04), RParah (MNI: 33, − 39, − 12; t = 5.43) and LAmygdala (MNI: − 21, 0, − 18; t = 3.90) under neutral condition; RInsula (MNI: 42, − 15, 18; t = 4.83), LParah (MNI: − 27, − 42, − 9; t = 5.76), RParah (MNI: 18, − 42, − 9; t = 5.15) and LAmygdala (MNI: − 18, − 6, − 18; t = 4.99) under positive condition (Table 2; Fig. 1).
Brain network analysis
Six graphs were obtained in total. Each graph was chosen with the lowest BIC score (BIC = 112 for HCs under positive, BIC = 72 for HCs under neutral, BIC = 79 for HCs under negative and BIC = 115 for GAD under positive, BIC = 106 for GAD under neutral, BIC = 0 for GAD under negative). Each graph represented how brain regions interacted under corresponding condition. Orientation of the connections represented causal relations between ROIs. Numbers along with the connections represent the strength. The connections highlighted in blue were shared by two groups under corresponding condition. Connections common to two groups included: LSFG-RSFG, LACC-LMFG, RInsula-LAmygdala for positive condition; LSFG-RSFG, LACC-LMFG, RSFG-LACC, LACC-LAmygdala, RInsula-RParah for neutral condition; LSFG-RSFG, LACC-LMFG, RParah-LAmygdala for negative condition. None of them was significantly different in strength between two groups (two sample t test, all p values > 0.2). Connections common to all conditions in GAD group included: LSFG-RSFG, LMFG-LACC, LParah-RInsula, RInsula-RParah. Four separate Repeated measures ANOVA with condition as within-subjects variable showed a main effect of condition only for connection RInsula-RParah (F (2, 30) = 10.02, p < 0.0001, Partial Eta Squared = 0.40), but not for connection LSFG-RSFG (F (2,30) = 0.16, p > 0.5), connection LMFG-LACC (F (2,30) = 1.55, p > 0.1), connection LParah-RInsula (F (2,30) = 0.38, p > 0.5). for connection RInsula-RParah, subsequent comparisons revealed that under positive condition, it was stronger than that of under neutral condition and negative condition (positive vs neutral, t (15) = − 3.1, p < 0.01, d = 0.80; positive vs negative t (15) = − 4.48, p < 0.001, d = 1.15) (both passed a Bonferroni multiple comparisons correction). Connection common to both positive and negative condition in GAD group was LMFG-RSFG, no significant difference was found between two conditions (paired t test, t (15) = − 0.98, p > 0.05). Connections common to both neutral and negative condition in GAD group included: RInsula-LACC, RParah-LMFG. None of them was found significant different between two conditions (paired t test, RInsula-LACC, t (15) = 1.17, p > 0.1; RParah-LMFG, t (15) = 0.49, p > 0.1).
the number of connections is eight vs seven (HCs vs GAD) for positive condition, seven vs eight for neutral condition and seven vs ten for negative condition. The brain networks in GAD were denser than that of HCs under neutral and negative condition (densest under negative condition), but sparser than that of HCs under positive condition. There were always bottom-up connections in GAD group across condition (one of them under positive condition, two of them under neutral condition, three of them under negative condition), but only one top-down connection found under neutral condition. Under both positive and negative condition, top-down connections originating from LSFG to limbic area was found in HCs, while no top-down connection was found in GAD under any of two condition. Actually, the connection from LSFG was either disconnected with LParah under positive condition, or superseded by bottom-up connections when it reached to RSFG in GAD. There was one top-down connection in GAD under neutral condition (from ACC to LAmygdala), but this was not found in HCs. All the bottom-up connections in GAD were getting through the RInsula, then to PFC area across condition, which brought the Rinsula as a critical conjunction node in brain networks of GAD group. While in HCs, bottom-up/top-down connections were always getting through the LAmygdala. LACC was another critical conjunction node in brain networks of HCs across three conditions, with three or four connections connected (like the LAmygdala). While in GAD, it was either LMFG or LACC or both of them. LAmygdala in GAD was less connected across three condition (only one connection under neutral/positive condition, two under negative condition) relative to HCs. In PFC subnetwork, the orientation of connection between LACC and LMFG was reversed in GAD across conditions (Fig. 2).
Correlation analysis found that, under positive condition, the strength of LSFG-RSFG was positively correlated with HAMA score (Pearson correlation coefficient = 0.53, p < 0.05; see Fig. 3); No other correlation was found. (all p values were two tailed).
In this study, we used brain network analysis to identify the differential brain network response to emotional stimuli between GAD patients and HCs. The brain network comprised of several pre-hypothesized brain regions that were consistently reported having abnormal activation/functional connectivity in GAD patients during emotional tasks. Firstly, we replicated previous findings by GLM analysis that PFC area and limbic area were found abnormally activated in GAD patients during emotional tasks. Secondly, by brain network analysis, there were three main findings: 1) relative to HCs, GAD patients exhibited more bottom-up connections but less top-down connection; 2) regardless of condition, the insula was more connected, but the amygdala was less connected in GAD patients, relative to HCs; 3) GAD patients showed a different brain network response from HCs under neutral condition.
In keeping with our first prediction, we found more bottom-up but less top-down connection in GAD patients, relative to HCs. We concluded that this more bottom-up connection response pattern (less top-down but more bottom-up connection) in the brain network of GAD patients was likely to indicate inadequate top-down control. Specifically, under positive condition, what was common to GAD patients and HCs was that both groups had one bottom-up connection from limbic area to PFC. the bottom-up connections from limbic area to PFC during emotional task was highly likely to indicate an excitation process, since limbic area is the key part responsible for emotion generating (for example: the amygdala, insula and hippocampus [1, 43, 44]). However, besides the bottom-up connection, HCs also had one top-down connection from PFC (LSFG) to limbic area, which was not seen in GAD patients. The difference was even bigger under negative condition, in which GAD patients showed three bottom-up connections but no top-down connection, while HCs showed one top-down connection (LSFG to RSFG to LACC to LAmygdala). PFC is considered to be a key player in the cognitive modulation of emotion [35, 45] and its projecting to limbic area is considered critical to emotion modulation . The top-down connections from LSFG to LAmygdala (to RInsula) shown in HCs under both positive and negative condition may demonstrate the cognitive modulation of emotion from PFC, preventing overresponse during emotional stimuli confronting . Relatively, the lack of top-down connection and only having bottom-up connection in GAD patients under both negative and positive condition possibly indicated the imbalance between PFC’s cognitive modulation and limbic area’s emotional reactivity. More specifically, more bottom-up connection from limbic area but less top-down connection from PFC in GAD patients likely reflected the “over-responsiveness” of limbic area and PFC’s hypo-function, respectively. (note: over-responsiveness here meant denser connectivity from limbic area to PFC).
Across three conditions, the insula was more connected, but the amygdala was less connected in GAD patients, relative to that of HCs. The insula is associated with interoception generation from the body and is implicated in all kinds of subjective feelings (e.g. includes emotional awareness) . Abnormal activation in insula is implicated in multiple kinds of anxiety disorders (e.g. social phobia , social anxiety disorder [24, 49]). Actually, it is proposed that inability to discriminate typical fluctuations in interoception from potential aversive body signal is key to anxiety disorder . In GAD, it has reported abnormal connectivity between insula and amygdala under emotional task  or resting state . GAD patients are characterized by excessive, uncontrollable and sustained worry, which is largely independent of environment . Besides, they often have other clinical symptoms, like feeling on edge, muscle tension, etc. which are all integrated by insula to produce interoception and subject feelings. Therefore, it was reasonable to indicate that the enhanced-connectivity in insula was a reflection of the abnormal processing of interoception in GAD .
The amygdala was less connected regardless of condition in GAD patients, relative to HCs. The amygdala is implicated in emotional processing, emotional learning and cognitive evaluation of emotional stimuli and it is closely connected with cortical and subcortical regions [52, 53]. In HCs, the connections of amygdala with both subcortical and cortical regions was preserved. But in GAD patients, the amygdala was either disconnected with PFC (under positive condition), or disconnected with limbic area (under neutral condition), or less connected with limbic area (under negative condition), this may imply the dysfunction of amygdala in GAD patients like other studies [7, 13]), but from a perspective of brain network response. In detail, we utilized a more comprehensive and rigorous method to identify connections, a Bayesian network-based approach, in which determination of any connection between pairs of brain regions was assessed against the backdrop of the whole brain network . This differed from the traditional correlation-based functional connectivity analyses as most studies used, which account for only two brain regions at a time and therefore can’t infer causal effect . Therefore, correlation-based analysis may more easily find connections between two nodes, but they cannot be identified by using IMaGES.
Across condition, the insula directly/indirectly predicted the activity of amygdala in GAD patients, a reverse pattern relative to HCs. As discussed above, the over-connected of insula was implicated in GAD patients’ abnormal processing of interoception. The result that insula predicted amygdala’s activity may indicate that incapable of discriminately processing typical interoception and aversive body signal lead to increased emotional response .
Under neutral condition, the brain network response in GAD patients was different from that of HCs. Specifically, two bottom-up and one top-down connection were found in GAD patients. Patients with anxiety disorder is impaired in discriminating between threat and safety , tend to overgeneralize conditioned fear, respond to neutral stimuli as if the neutral stimuli are threat-related . More bottom-up connections in GAD patients may indicate that GAD patients overgeneralized fearful reaction to nonthreatening stimuli, mistakenly treated the neutral human face as it was fearful, which triggered over-responsiveness in limbic area (reflected by the number of bottom-up connections). One top-down connection from ACC to amygdala was found too. ACC is implicated in emotional regulation and threat-related appraisal [16, 56]. However, we suggested the top-down connection under neutral condition was unlikely to reflect ACC’s regulation over amygdala’s activity. Specifically, in this subnetwork, it was the insula predicted ACC’s activity (directly/indirectly), then the ACC got to predict amygdala’s activity. Both connections from the insula to ACC may be caused by fear overgeneralization (as discussed above), therefore the ACC’s activity predicted by insula may implicate threat-related appraisal (rather than emotional regulation), which then progressively triggered amygdala’s activity. Consistent with this “not regulation” explanation, study using fear generalization paradigm has also found that GAD patients showed insufficient PFC regulation over fear-generalized stimuli .
Inconsistent with the second prediction, among the connections that were common to both GAD patients and HCs, we didn’t find any pair of connections differ in their strength. The connection between LSFG and RSFG was exactly the same between two groups (same connection with the same orientation), and it didn’t vary across condition. No significant difference in the strength of this connection between two groups may indicate the communication between LSFG and RSFG was preserved in GAD patients. However, this connection was positively correlated with HAMA score in GAD patients under positive condition. To the best of our knowledge, this was the first time that connectivity between LSFG and RSFG was reported to positively correlate with anxiety level. This added evidence to the significance of lateral PFC in clinical symptom. Regardless of condition, the LMFG was more connected in GAD patients, compared with LACC, a reverse pattern in HCs. The medial PFC is implicated in introspective thinking and a dysregulation of this region may underlie the uncontrollable worry in GAD patients .
Together with GLM results. In this study, the amygdala’s activity during tasks did not differ between GAD patients and HCs. No difference in amygdala’s activity between two groups has been observed in previous studies [12, 59], it may need specific fearful stimuli to elicit abnormal amygdala response in GAD patients . We didn’t find increased activity in amygdala during presentation of fearful faces in either group, which is somewhat at odds with previous studies [12, 60, 61]. Between group analysis revealed that both limbic area and PFC were under-activated in GAD patients regardless of condition. Integrated with the brain network response pattern, it may suggest that GAD patients’ brain network response pattern is rigid, showing more bottom-up connections even when the limbic area is not as active as in HCs.
Several limitations are important to note. First, the accuracy and reaction time of the behavioral results were not recorded in this study, as other studies [59, 62]. Lack of behavioral results may make the result less convincing than it should be. Second, the relatively small sample size. Twenty-three GAD patients and 19 HCs were recruited. However, because of poor data quality and excessive head motion during scanning by couple of participants, the final sample size was 16 GAD patients and 14 HCs and this may be the reason no other significant correlation between clinical measurements and connectivity was found. Last, no mood rating was administered with participants after scanning, which may confine the explanation of the result.
Based on the current results, there are several implications for future research. First, compared with other anxiety disorders (e.g. social anxiety disorder), less studies have found aberrations in the insula of GAD patients. However, in this study, we found a highly distinct role of insula in GAD patients during emotional tasks. As indicated, the insula is implicated in interoception and has close link with autonomic system , future studies may benefit from combining fMRI technique with physiological measurements. Second, in this study, brain network analysis provided several new insights on GAD’s pathology (e.g. the aberrations in the insula). However, GAD and other anxiety disorders (e.g. social anxiety disorder and panic disorder) have been reported to share similar underlying neural mechanism [23, 24]. For example, one study has found that increased amygdala activation was found in patients with different diagnosis (GAD or panic disorder or social anxiety disorder) during fearful faces processing . This may imply over-response in limbic area is common to all kinds of anxiety disorder. However, in the same study they have also found patients with panic disorder showed a unique hyperactivation in the posterior insula and a shared dysregulation in temporal regions by patients with panic disorder or social anxiety disorder, an implication of disorder specific neural mechanism. In another study, researchers haven’t found GAD specific, nor panic-disorder specific neural mechanism during emotional regulation task. Both groups showed a hypoactivation in PFC, implying insufficient top-down control . These two studies may indicate that different anxiety disorders share the same hypoactivation in PFC and overresponse in limbic area. But in the meantime, they also have disorder specific neural mechanism (like the insula). Hence, it is important for future studies to use the same kind of brain network analysis to identify disorder-common and -specific circuitry among several anxiety disorders, for it not only can better identify different anxiety disorders based on neuroimaging techniques, but also can drive more specialized treatment to each disorder.
Third, by including neutral condition in the brain network analysis, we found evidence to fear overgeneralization hypothesis concerning GAD patients. However, in addition to fear overgeneralization, findings of more bottom-up connections in GAD patients under neutral condition also aligned well with the prediction of intolerance of uncertainty. Studies found GAD patients are intolerant of uncertainty , tend to have negative feelings while facing uncertainty  and reduction in the score of intolerance of uncertainty often predicted the improvements of GAD symptoms [66, 67]. Neutral facial expressions are the most ambiguous of all the stimuli and more likely to be perceived as negative even in health controls , study has found. Hence, undoubtedly, intolerance of uncertainty in GAD patients could also be a main contributor to the findings under neutral condition. Last, neutral facial expressions condition was often treated as “neutral” in previous studies of GAD. However, the current study provided direct evidence against that setting, at least from the perspective of neuroimaging.
Through brain network analysis, we found that relative to HCs, GAD patients exhibited more bottom-up connections but less top-down connection in their emotional reactivity brain network, which may indicate inadequate top-down control. In addition, compared with HCs, the insula was more connected with other brain regions in GAD patients, which may indicate abnormal interoception processing in GAD patients. Lastly, brain network response was different between GAD patients and HCs under neutral condition, possibly an indicative of fear overgeneralization and intolerance of uncertainty in GAD patients.
Availability of data and materials
The task fMRI datasets used during the current study are available from the corresponding author on reasonable request.
Firstly, hippocampus is considered part of the limbic system responsible for fear generation and was part of common neurocircuitry implicated in all anxiety disorders ; secondly, connectivity between insula and amygdala was impaired in GAD patients [9, 33]; thirdly, both insula and hippocampus showed increased activation under emotional task in GAD patients .
Generalized anxiety disorder
Functional magnetic resonance imaging
Anterior cingulate cortex
Multiple sample Greedy Equivalence Search
Linear non-gaussian Orientation Fixed Structure
Region of interest
Hamilton anxiety rating scale
Hamilton depression rating scale
Left/right Superior frontal gyrus
Left Medial frontal gyrus
Duval ER, Javanbakht A, Liberzon I. Neural circuits in anxiety and stress disorders: a focused review. Ther Clin Risk Manag. 2015;11:115–26.
Bell CC, Bell CC. Diagnostic and statistical manual of mental disorders, fourth edition, text revision: DSM-IV-TR quick reference to the diagnostic criteria from DSM-IV-TR. JAMA. 2001;285(6):811–2.
Behar E, DiMarco I, Hekler E, Mohlman J, Staples A. Current theoretical models of generalized anxiety disorder (GAD): conceptual review and treatment implications. J Anxiety Disord. 2009;23:1011–23.
Mennin DS, Heimberg RG, Turk CL, Fresco DM. Preliminary evidence for an emotion dysregulation model of generalized anxiety disorder. Behav Res Ther. 2005;43(10):1281–310.
Mochcovitch MD, da Rocha Freire RC, Garcia RF, Nardi AE. A systematic review of fMRI studies in generalized anxiety disorder: evaluating its neural and cognitive basis. J Affect Disord. 2014;167:336–42.
Hilbert K, Lueken U, Beesdo-Baum K. Neural structures, functioning and connectivity in generalized anxiety disorder and interaction with neuroendocrine systems: a systematic review. J Affect Disord. 2014;158:114–26.
Monk CS, Telzer EH, Mogg K, Bradley BP, Mai X, Louro HM, et al. Amygdala and ventrolateral prefrontal cortex activation to masked angry faces in children and adolescents with generalized anxiety disorder. Arch Gen Psychiatry. 2008;65(5):568–76.
Blair K, Shaywitz J, Smith BW, Rhodes R, Geraci M, Jones M, et al. Response to emotional expressions in generalized social phobia and generalized anxiety disorder: evidence for separate disorders. Am J Psychiatr. 2008;165(9):1193–202.
McClure EB, Monk CS, Nelson EE, Parrish JM, Adler A, Blair RJ, et al. Abnormal attention modulation of fear circuit function in pediatric generalized anxiety disorder. Arch Gen Psychiatry. 2007;64(1):97–106.
Strawn J, Buszek S, Weber W, Chu W-J, Whitsel R, Adler C, et al. Neurocircuitry of generalized anxiety disorder in adolescents: a pilot functional neuroimaging and functional connectivity study. Depress Anxiety. 2012;29(11):939–47.
Buff C, Brinkmann L, Neumeister P, Feldker K, Heitmann C, Gathmann B, et al. Specifically altered brain responses to threat in generalized anxiety disorder relative to social anxiety disorder and panic disorder. NeuroImage Clin. 2016;12:698–706.
Palm ME, Elliott R, McKie S, Deakin JF, Anderson IM. Attenuated responses to emotional expressions in women with generalized anxiety disorder. Psychol Med. 2011;41(5):1009–18.
Robinson OJ, Krimsky M, Lieberman L, Allen P, Vytal K, Grillon C. Towards a mechanistic understanding of pathological anxiety: the dorsal medial prefrontal-amygdala ‘aversive amplification’ circuit in unmedicated generalized and social anxiety disorders. Lancet Psychiatry. 2014;1(4):294–302.
Ball TM, Ramsawh HJ, Campbell-Sills L, Paulus MP, Stein MB. Prefrontal dysfunction during emotion regulation in generalized anxiety and panic disorders. Psychol Med. 2013;43(7):1475–86.
Blair KS, Geraci M, Smith BW, Hollon N, DeVido J, Otero M, et al. Reduced dorsal anterior cingulate cortical activity during emotional regulation and top-down Attentional control in generalized social phobia, generalized anxiety disorder, and comorbid generalized social phobia/generalized anxiety disorder. Biol Psychiatry. 2012;72(6):476–82.
Etkin A, Prater KE, Hoeft F, Menon V, Schatzberg AF. Failure of anterior cingulate activation and connectivity with the amygdala during implicit regulation of emotional processing in generalized anxiety disorder. Am J Psychiatr. 2010;167(5):545–54.
Etkin A, Schatzberg A. Common abnormalities and disorder-specific compensation during implicit regulation of emotional processing in generalized anxiety and major depressive disorders. Am J Psychiatry. 2011;168:968–78.
Lewis AR, Zinbarg RE, Durbin CE. Advances, problems, and challenges in the study of emotion regulation: a commentary. J Psychopathol Behav Assess. 2010;32(1):83–91.
Goossen B, van der Starre J, Van der Heiden C. A review of neuroimaging studies in generalized anxiety disorder: “so where do we stand?”. J Neural Transm. 2019;126:1203–16.
Wang HY, Zhang XX, Si CP, Xu Y, Liu Q, Bian HT, et al. Prefrontoparietal dysfunction during emotion regulation in anxiety disorder: a meta-analysis of functional magnetic resonance imaging studies. Neuropsychiatr Dis Treat. 2018;14:1183–98.
Reid AT, Headley DB, Mill RD, Sanchez-Romero R, Uddin LQ, Marinazzo D, et al. Advancing functional connectivity research from association to causation. Nat Neurosci. 2019;22(11):1751–60.
Mills-Finnerty C, Hanson C, Hanson SJ. Brain network response underlying decisions about abstract reinforcers. Neuroimage. 2014;103:48–54.
Fonzo G, Ramsawh H, Flagan T, Sullivan S, Letamendi A, Simmons A, et al. Common and disorder-specific neural responses to emotional faces in generalised anxiety, social anxiety and panic disorders. Br J Psychiatry. 2015;206:206–15.
Ball TM, Sullivan S, Flagan T, Hitchcock CA, Simmons A, Paulus MP, et al. Selective effects of social anxiety, anxiety sensitivity, and negative affectivity on the neural bases of emotional face processing. NeuroImage. 2012;59(2):1879–87.
Smith SM, Miller KL, Salimi-Khorshidi G, Webster M, Beckmann CF, Nichols TE, et al. Network modelling methods for FMRI. Neuroimage. 2011;54(2):875–91.
Ramsey JD, Hanson SJ, Hanson C, Halchenko YO, Poldrack RA, Glymour C. Six problems for causal inference from fMRI. Neuroimage. 2010;49(2):1545–58.
Ramsey JD, Hanson SJ, Glymour C. Multi-subject search correctly identifies causal connections and most causal directions in the DCM models of the Smith et al. simulation study. Neuroimage. 2011;58(3):838–48.
Manelis A, Almeida JRC, Stiffler R, Lockovich JC, Aslam HA, Phillips ML. Anticipation-related brain connectivity in bipolar and unipolar depression: a graph theory approach. Brain. 2016;139:2554–66.
Manelis A, Reder LM. Effective connectivity among the working memory regions during preparation for and during performance of the n-back task. Front Hum Neurosci. 2014;8:593.
Boukrina O, Hanson SJ, Hanson C. Modeling activation and effective connectivity of VWFA in same script bilinguals. Hum Brain Mapp. 2014;35(6):2543–60.
Boukrina O, Graves WW. Neural networks underlying contributions from semantics in reading aloud. Front Hum Neurosci. 2013;7:518.
Sun J, Hu X, Huang X, Liu Y, Li K, Li X, et al. Inferring consistent functional interaction patterns from natural stimulus FMRI data. Neuroimage. 2012;61(4):987–99.
Roy AK, Fudge JL, Kelly C, Perry JSA, Daniele T, Carlisi C, Benson B, Castellanos FX, Milham MP, Pine DS, Ernst M. Intrinsic functional connectivity of amygdala-based networks in adolescent generalized anxiety disorder. J Am Acad Child Adolesc Psychiatry. 2013;52(3):290–9.
Buff C, Schmidt C, Brinkmann L, Gathmann B, Tupak S, Straube T. Directed threat imagery in generalized anxiety disorder. Psychol Med. 2018;48(4):617–28.
Davidson RJ. Anxiety and affective style: role of prefrontal cortex and amygdala. Biol Psychiatry. 2002;51(1):68–80.
Bai L, Ma H, Huang YX, Luo YJ. The development of native Chinese affective picture system-a pretest in 46 college students. Chin Ment Health J. 2005;19:719–22.
Maldjian JA, Laurienti PJ, Burdette JH. Precentral gyrus discrepancy in electronic versions of the Talairach atlas. NeuroImage. 2004;21(1):450–5.
Maldjian JA, Laurienti PJ, Kraft RA, Burdette JH. An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets. NeuroImage. 2003;19(3):1233–9.
Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical Parcellation of the MNI MRI single-subject brain. NeuroImage. 2002;15(1):273–89.
Mumford JA, Ramsey JD. Bayesian networks for fMRI: a primer. Neuroimage. 2014;86:573–82.
Xiao-Wei S, Zhang-Ye D, Xiang-Yu L, Su-Fang L, Xi-Nian Z, Chao-Zhe Z, et al. REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS One. 2011;6(9):1–12.
Anderson TW, Darling DA. Asymptotic theory of certain “goodness of fit” criteria based on stochastic processes. Ann Math Stat. 1952;23.
Wager T, Davidson M, Hughes B, Lindquist M, Ochsner K. Prefrontal-subcortical pathways mediating successful emotion regulation. Neuron. 2008;59:1037–50.
Etkin A, Wager T. Functional neuroimaging of anxiety: a meta-analysis of emotional processing in PTSD, social anxiety disorder, and specific phobia. Am J Psychiatry. 2007;164:1476–88.
Ochsner KN, Bunge SA, Gross JJ, Gabrieli JDE. Rethinking feelings: an FMRI study of the cognitive regulation of emotion. J Cogn Neurosci. 2002;14(8):1215–29.
Hartley CA, Phelps EA. Changing fear: the neurocircuitry of emotion regulation. Neuropsychopharmacol. 2010;35(1):136–46.
Craig AD. How do you feel — now? The anterior insula and human awareness. Nat Rev Neurosci. 2009;10(1):59–70.
Lueken U, Kruschwitz JD, Muehlhan M, Siegert J, Hoyer J, Wittchen H-U. How specific is specific phobia? Different neural response patterns in two subtypes of specific phobia. NeuroImage. 2011;56(1):363–72.
Boehme S, Ritter V, Koranyi S, Stangier U, Strauss B, Miltner W, et al. Brain activation during anticipatory anxiety in social anxiety disorder. Soc Cogn Affect Neurosci. 2013;9:1413–8.
Paulus MP, Stein MB. Interoception in anxiety and depression. Brain Struct Funct. 2010;214(5–6):451–63.
Association AP. Diagnostic and Statistical Manual of Mental Disorders 4th edition, text revision (DSM-IV-TR). 2000.
Davis M, Whalen PJ. The amygdala: vigilance and emotion. Mol Psychiatry. 2001;6(1):13–34.
Gallagher M, Chiba AA. The amygdala and emotion. Curr Opin Neurobiol. 1996;6(2):221–7.
Cha J, Greenberg T, Carlson JM, Dedora DJ, Hajcak G, Mujica-Parodi LR. Circuit-wide structural and functional measures predict ventromedial prefrontal cortex fear generalization: implications for generalized anxiety disorder. J Neurosci. 2014;34(11):4043–53.
Lissek S. Toward an account of clinical anxiety predicated on basic, neutrally mapped mechanisms of pavlovian fear-learning: the case for conditioned overgeneralization. Depress Anxiety. 2012;29(4):257–63.
Kalisch R, Gerlicher AM. Making a mountain out of a molehill: on the role of the rostral dorsal anterior cingulate and dorsomedial prefrontal cortex in conscious threat appraisal, catastrophizing, and worrying. Neurosci Biobehav Rev. 2014;42:1–8.
Greenberg T, Carlson JM, Cha J, Hajcak G, Mujica-Parodi LR. Ventromedial prefrontal cortex reactivity is altered in generalized anxiety disorder during fear generalization. Depress Anxiety. 2013;30(3):242–50.
Paulesu E, Torti T, Danelli L, Ferri F, Scialfa G, Sberna M, et al. Neural correlates of worry in generalized anxiety disorder and in normal controls: a functional MRI study. Psychol Med. 2009;40:117–24.
Whalen PJ, Johnstone T, Somerville LH, Nitschke JB, Polis S, Alexander AL, et al. A functional magnetic resonance imaging predictor of treatment response to venlafaxine in generalized anxiety disorder. Biol Psychiatry. 2008;63(9):858–63.
Whalen PJ, Shin LM, McInerney SC, Fischer H, Wright CI, Rauch SL. A functional MRI study of human amygdala responses to facial expressions of fear versus anger. Emotion (Washington, DC). 2001;1(1):70–83.
Adolphs R, Tranel D, Damasio H, Damasio AR. Fear and the human amygdala. J Neurosci. 1995;15(9):5879–91.
Holzel BK, Hoge EA, Greve DN, Gard T, Creswell JD, Brown KW, et al. Neural mechanisms of symptom improvements in generalized anxiety disorder following mindfulness training. NeuroImage Clin. 2013;2:448–58.
de Morree H, Rutten G-J, Szabó B, Sitskoorn M, Kop W. Effects of insula resection on autonomic nervous system activity. J Neurosurg Anesthesiol. 2015;28:153–8.
Dugas MJ, Gagnon F, Ladouceur R, Freeston MH. Generalized anxiety disorder: a preliminary test of a conceptual model. Behav Res Ther. 1998;36(2):215–26.
Dugas MJ, Koerner N. Cognitive-behavioral treatment for generalized anxiety disorder: current status and future directions. J Cogn Psychother. 2005;19(1):61–81.
Hebert EA, Dugas MJ. Behavioral experiments for intolerance of uncertainty: challenging the unknown in the treatment of generalized anxiety disorder. Cogn Behav Pract. 2019;26(2):421–36.
Boswell JF, Thompson-Hollands J, Farchione TJ, Barlow DH. Intolerance of uncertainty: a common factor in the treatment of emotional disorders. J Clin Psychol. 2013;69(6):630–45.
Lee E, Kang JI, Park IH, Kim JJ, An SK. Is a neutral face really evaluated as being emotionally neutral? Psychiatry Res. 2008;157(1–3):77–85.
We deeply thank Dr. Christina Hines (MD) for her contribution in language editing.
This study was supported by National Natural Science Foundation of China (81971289, 81871344); Natural Science Foundation of Jiangsu Province (BK20161109); the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province, China (18KJB190003); key research and development program (Social Development) project of Jiangsu province (BE2019609). The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. And this research was approved by the Medical Research Ethics Committee of Nanjing Brain Hospital affiliated to Nanjing Medical University. Written informed consent was obtained from all individual participants included in the study.
Consent for publication
The authors declare that they have no competing interests.
About this article
Cite this article
Li, J., Zhong, Y., Ma, Z. et al. Emotion reactivity-related brain network analysis in generalized anxiety disorder: a task fMRI study. BMC Psychiatry 20, 429 (2020). https://doi.org/10.1186/s12888-020-02831-6
- Generalized anxiety disorder
- Emotion reactivity
- Brain network analysis
- Top-down control
- Task-based fMRI