Mood disorders (major depressive disorder and bipolar disorder), posttraumatic stress disorder, and anxiety disorders (generalized anxiety disorder, panic disorder, agoraphobia, and specific and social phobia) are highly comorbid (
1) and collectively account for more than 65% of nonfatal disease burden attributable to psychiatric disorders (
2). Up to 90% of patients with an anxiety disorder meet criteria for a concurrent mood disorder (
2,
3), and as many as 70% of individuals with mood disorders meet criteria for an anxiety disorder during their lifetime (
4,
5). Negative affective states are shared and central clinical features of these disorders (
6), including bipolar disorder, where depressive symptoms are the dominant psychopathology (
7).
Meta-analyses of brain imaging studies on mood, posttraumatic stress, and anxiety disorders have shown that each of these disorders is associated with abnormalities in task-related brain engagement (summarized in eTable 1 in the Supplement). The findings of these diagnosis-specific meta-analyses show conspicuous divergence (eTable 1 in the Supplement) that has been attributed to low numbers of contributing studies, reporting bias from region-of-interest (ROI) analyses, and inadequate correction for multiple comparisons (
8,
9). Of note, methodological improvements over time have led to a progressive reduction in the number clusters of case-control differences reported in diagnosis-specific meta-analyses (eFigure 1 in the Supplement). Using data from task-related functional magnetic resonance imaging (fMRI) studies published in the last 15 years, we demonstrated that diagnostic differences in the brain regions implicated in mood and anxiety disorders largely reflected the association with ROI analyses (
9). By contrast, when only whole-brain analyses were considered, there were large pairwise correlations between the diagnosis-specific profiles (ρ range, 0.79-0.82; all
P < .001) (
9).
Here, we extend this line of research in 2 distinct ways. First, we sought to identify brain regions where aberrant task-related activation was most likely to show transdiagnostic convergence across major depressive disorder, bipolar disorder, and anxiety and posttraumatic stress disorders. To achieve this, we capitalized on activation likelihood estimation (ALE) meta-analytic tools (
10-
13) to synthesize coordinates of case-control differences in what is, to our knowledge, the largest sample of fMRI articles comprising the body of the relevant literature over the last 15 years. Second, we anchored the analysis plan to the Research Domain Criteria (RDoC) framework (
14) proposed by the US National Institute of Mental Health. The RDoC framework is the best approximation to a criterion approach to the classification of the array of activation tasks used in the primary studies and enables a principled interpretation of results in terms of dysfunction in clearly defined cognitive processes. Based on current neurobiological (
15,
16), we predicted that transdiagnostic clusters of aberrant brain activity would converge in regions within the prefrontal, insular, and anterior cingulate cortex and in subcortical regions (particularly the amygdala/hippocampus and striatum) that support the adaptive regulation of cognition and affect.
Method
Literature Search and Article Eligibility
We applied the Preferred Reporting Items for Systematic Reviews and Meta-analyses criteria (
http://www.prismastatement.org/) to identify articles that used whole-brain analyses of task-related fMRI to compare healthy adults with adult patients who received a diagnosis of major depressive disorder, bipolar disorder, generalized anxiety disorder, panic disorder, agoraphobia, specific and social phobias, and posttraumatic stress disorder (details of the search and article eligibility criteria in the eMethods and eFigure 2 in the Supplement). Because we used data from published studies, no institutional review board approval was sought and patient consent was not obtained.
Database Construction
We use the term article to denote the published manuscript and the term experiment to denote the coordinates of case-control differences reported in each article. Accordingly, from each article, we extracted coordinates of case-control differences derived from whole-brain analyses only. These were then coded according to the strength of the magnetic field of the scanner, the diagnostic classification system, symptom severity, the direction of change in brain activity in patients compared with healthy individuals (hypoactivation or hyperactivation), and the corresponding RDoC domain and construct. The coding of tasks according to their corresponding RDoC domain and construct is described in the eMethods and shown in eTable 2 in the Supplement. For example, tasks such as the n-back and the Sternberg were assigned to the construct of working memory and the domain of cognitive systems where as various facial affect processing tasks were assigned to the construct of social communication and the domain of social processes. This allowed us to create 3 groupings of tasks labeled by their type (eg, facial affect processing) and their RDoC construct and domain. For each article, the symptom severity of the clinical sample was coded based on the mean psychopathology rating reported. To accommodate the use of different instruments across studies and clinical populations, symptom levels were labeled as minimal/mild, moderate, or severe (details in eMethods in the Supplement). Furthermore, for each experiment, but separately for patient and control groups, we coded their diagnostic status, sample size, age, and sex (percentage of men). In patients, medication status was coded as the percentage of patients receiving any psychotropic medication in each study sample. Further details of the database construction are provided in the eMethods in the Supplement.
Activation Likelihood Estimation
We used ALE, implemented in MATLAB (MathWorks), to test whether the whole-brain coordinates of case-control differences across experiments and disorders converged into discrete clusters with a nonrandom spatial distribution (
10-
13). The fundamental assumption of the ALE is that each voxel has the same a priori chance of differentiating patients from control individuals (null hypothesis). Consequently, ROI analyses were excluded because they violate this assumption and their inclusion would artificially bias results in favor of voxels within these regions. The main outcome of the ALE analysis are the clusters (ie, grouping of brain regions) in which the coordinates of the experiments converge. Per best-practice standards (
10-
13) statistically significant clusters were identified using a cluster-level family wise error–corrected threshold of
P less than .05 (cluster-forming threshold at voxel-level
P < .001). Additionally, for each suprathreshold cluster, we extracted the per voxel probability of functional change from the modeled activation maps. These values represent the probability of identifying a functional change for a mean voxel within the clusters derived from the modeled activation maps. Details of the procedures involved are described in the eMethods in the Supplement.
We analyzed coordinates of hypoactivation or hyperactivation in patients compared with healthy individuals separately to enhance interpretability. First, we identified suprathreshold clusters of hypoactivation and hyperactivation by pooling coordinates from all diagnoses and tasks and then conducted follow-up analyses to identify the effect of moderators. For the follow-up analyses, we extracted per-voxel probabilities of functional change for each cluster and conducted nonparametric Kruskal-Wallis tests and Spearman correlations to calculate the contribution of age, sex, RDoC domain/construct, diagnosis, symptom severity, and medication.
In generating the modeled activation maps, we pooled coordinates across diagnoses for 2 reasons. First, the disorders considered here are highly comorbid and hence the pooled analyses accommodate uncertainty about their symptomatic and syndromal boundaries. Moreover, comorbidity is not always reported in primary studies and therefore it is difficult to estimate its prevalence in the samples examined and its potential contribution to the neuro imaging results. Second, pooling results across diagnoses balances power, specificity, and sensitivity and allows for a data-driven quantification of the diagnosis-specific contribution to each suprathreshold cluster. We conducted supplemental diagnosis-specific analyses, which are presented in the eMethods and eResults in the Supplement.
In generating the modeled activation maps, we pooled coordinates from all the tasks used in the primary experiments based on 2 considerations. First, traditional neuropsychologic formulations tend to consider cognitive tasks as relatively specific to a particular process. Advances in cognitive and affective neuroscience have led to the recognition that the association between brain structure and function is pluripotent (one-to-many) and degenerate (many-to-one) (
17,
18). Therefore, any given task engages brain regions outside those predicted by the cognitive mechanisms attributed to that particular task, while a single brain area may be activated by disparate tasks that may not share cognitive components (
17,
18). Our approach accommodates pluripotency and offers a more realistic representation of the relevance of cognitive domains to case-control differences. Following the identification of suprathreshold clusters from the pooled analyses, we estimated the contribution of tasks to each cluster. For these follow-up analyses, tasks were grouped according to their assigned RDoC domain/construct; the use of the RDoC frame-work provided an organizing principle for the multitude of tasks used in the primary studies. Compared with other classifications that are primarily driven by convention, the RDoC framework has a clearly defined origin and rationale (
14).
Finally, we used an alternate meta-analytic algorithm to confirm the reproducibility of the results of the main analyses and conducted several ancillary meta-analyses focusing on each diagnosis separately and using alternate classification of tasks (described in the eMethods, eResults, and eTable 8 in the Supplement).
Discussion
Meta-analyses of 367 task-related fMRI experiments in mood disorders, posttraumatic stress disorder, and anxiety disorders, comprising data from 4507 patients and 4755 control individuals, detected statistically robust transdiagnostic clusters of hypoactivation in the inferior prefrontal cortex/ insula, the inferior parietal lobule, and the putamen. These regions are part of a right-dominant brain system that supports contextual shifting and stopping of mental operations and behavioral responses (
19-
26). Specifically, the right inferior prefrontal cortex is critically involved in the inhibition of contextually inappropriate cognitive, affective, and motor responses (
19-
21); similarly, the putamen, particularly on the right, is essential for terminating contextually inappropriate motor and cognitive processes (
22). The anterior insula (
23,
24) and the right inferior parietal lobule (
25,
26) participate in the generation of salience-related signals that either initiate or terminate the engagement of attentional and working memory networks in response to changing demands. The insula, and particularly the anterior portion on the right, is thought to have a major role in integrating interoceptive information with information from other brain regions, thus supporting the formation of the conscious experience of an embodied self (
27,
28). This awareness of negative emotional states may act as a salient trigger for the insula and the adjacent inferior frontal regions to engage mechanisms of cognitive control. Notably, experiments involving domains of nonaffective cognition, affective processing, and social cognition showed a similar range of contributions to these clusters of hypoactivation (respective range: 28%-50% and 29%-56%) (
Figure 1B; eTable 7 in the Supplement). We therefore infer that the dominant abnormality in mood disorders, post traumatic stress disorder, and anxiety disorders involves a diagnosis-general disruption in salience processing (including interoceptive processing) and inhibitory control. These results contradict early hypotheses, which stipulated that affective morbidity results from right-sided fronto-parietal hyperactivity in response to negative/withdrawal stimuli (
29), but are in line with evidence that emphasizes the role of deficient cognitive control (
30-
32).
Further support derives from studies showing that deficits in the ability to stop and shift ongoing affective states and thoughts are the most significant predictors of affective symptoms and syndromes (
33-
35). Neurocognitive studies in mood and anxiety disorders also indicate a general disruption in cognitive control because they consistently report deficits of large effect size in stopping and shifting responses in a range of tasks (
32,
36,
37). Thus, impaired engagement of brain regions that subserve salience processing and inhibitory control present a plausible explanatory mechanism for the affective and nonaffective abnormalities observed in patients. In a separate meta-analysis (
38) of functional neuroimaging studies that was limited only to tasks of cognitive control, hypoactivation in the right inferior prefrontal/insular cortex was also reported as a transdiagnostic feature of schizophrenia, bipolar disorder, major depressive disorder, anxiety disorders, and substance use (
38). When considered together, these findings point to the possibility that abnormalities in brain regions involved in switching and stopping may underpin the vulnerability to develop any and all forms of psychopathology. Interestingly, similar arguments have been put forward for a single dimension of psychopathology, termed factor
p, as a main predictor of individuals’ liability for all mental disorders (
39). The relationship between the
p factor and disrupted engagement in salience/inhibitory control regions presents an intriguing avenue for future research.
We also identified 3 transdiagnostic clusters of hyperactivation in patients compared with healthy individuals centered in the left amygdala/parahippocampal gyrus, the left thalamus, and the perigenual/dorsal anterior cingulate cortex that were attributable mainly to experiments mapping to RDoC domains relating to affective and social processing (
Figure 2B). The clusters identified appear plausible because they comprised regions consistently associated with affective morbidity (
15,
16). The perigenual/dorsal anterior cingulate cortex is known to exert a regulatory influence on emotional experience and appraisal (
40) while the amygdala and parahippocampal gyrus, particularly on the left, are involved in emotional memory formation and retrieval (
41). The dorsal anterior cingulate cortex is also closely involved in the generation of internal autonomic and their associated expressive emotional responses (
42). Its relative hyperactivation in patients is consistent with the notion of increased arousal in response to stress that may be a trait feature of mood and anxiety disorders but may also reflect increased stress response to the fMRI tasks (
43).
Notwithstanding, these clusters were only detectable at a liberal statistical threshold, indicating greater inconsistency across primary studies that may indicate that hyperactivation in patients compared with healthy individuals may be more sensitive to variations in fMRI task design (eg, type or duration of stimuli or task instructions) and neuroimaging acquisition and analysis parameters. Detailed investigations that could directly address these issues would require more data than are currently available in the entire literature corpus.
Limitations
We placed substantial emphasis on the rigor and reproducibility of our methods to address ongoing concerns about the disparity in the number and localization of clusters in previous meta-analyses (eTable 1 in the Supplement). To further enhance reproducibility, we classified experiments based on the RDoC framework, which offers a structured approach to classification for fMRI tasks in future studies. We only included studies in adults, and therefore these findings may not generalize to pediatric or geriatric groups. We did not consider studies that failed to find case-control differences because such practice could only be justified if negative studies were sufficiently powered. We did not find an effect of symptom severity on the transdiagnostic clusters. This observation should be viewed with caution because of the variable instruments used to rate psychopathology and the reliance on group means from each study sample. Medication status did not have a statistically significant moderator association with the results reported. Medication has been shown to have mostly normalizing effects (
44) and may have attenuated case-control differences in the primary studies. We examined disorders with significant symptomatic and syndromal overlap for which we had comparable amount of data across diagnoses. We decided to exclude task-related fMRI articles on schizophrenia because the disproportionately larger number of studies (>250)
9 would have skewed the results. Given the observed power in this study, the results are statistically robust, but as the literature expands it is possible that additional transdiagnostic or disease specific clusters may emerge.