Evidence for Dissociable Linkage of Dimensions of Psychopathology to Brain Structure in Youths

A total of 1,601 participants completed multimodal neuroimaging as part of the PNC (15, 16). The institutional review boards of the University of Pennsylvania and the Children’s Hospital of Philadelphia approved the study procedures. Participants over age 18 provided written informed consent after receiving a complete description of the study, and we obtained assent and parental consent for minors. The final sample consisted of 1,394 youths; the demographic characteristics of the sample are summarized in Table 1 (see the online supplement for details). In this final sample, 155 participants (11%) were taking psychotropic medications at the time of imaging and were evaluated in sensitivity analyses, as described below.

As described previously (15, 16), psychiatric symptoms were assessed using a structured screening interview (GOASSESS) based on a modified version of the Kiddie Schedule for Affective Disorders and Schizophrenia (see the online supplement). The lifetime prevalence of each disorder is summarized in Table 1. An exploratory factor analysis of 112 item-level symptoms (see the online supplement) identified four correlated dimensions of psychopathology: anxious-misery, psychosis, behavioral, and fear (Figure 1A), which show a high degree of overlap across dimensions and diagnostic screening categories. To increase specificity, we then used a confirmatory item bifactor analysis, which yielded five orthogonal dimensions of psychopathology: anxious-misery, psychosis, behavioral/externalizing, fear, and overall psychopathology (symptom burden across all psychiatric disorders; see Figure 1B and 1C). Additionally, three cognitive factors (executive function/complex reasoning, social cognition, and episodic memory) derived from factor analysis (see the online supplement) were included as predictors in the multivariate analyses described below.

Image acquisition and processing have been reported in detail elsewhere (15). Briefly, imaging data were acquired on the same MRI scanner using the same imaging sequences for all participants (see the online supplement). Three highly trained image analysts independently assessed structural image quality control using manually derived ratings (18) (see the online supplement), and the average quality rating across the three raters was included as a covariate in all models to control for the confounding influence of variation in image quality.

Structural image processing utilized Advanced Normalization Tools (ANTs; see the online supplement). Previous large-scale evaluation studies have shown that this procedure is highly accurate and more sensitive to individual differences over the lifespan than comparable techniques. To avoid registration bias and maximize sensitivity, a custom adolescent template and tissue priors were used. This procedure yielded two maps for each subject: a cortical thickness image and a volume image (log transformed determinant of the Jacobian of the deformation field).

We used nonnegative matrix factorization (NMF) to identify networks where brain structure covaried consistently across participants (17). Details regarding the implementation of NMF have been presented elsewhere (17, 19) and in the online supplement (see Figure S1 for a schematic of the NMF procedure). NMF yields networks that are highly interpretable and have improved statistical power compared with standard mass univariate analyses (see the online supplement for details).

NMF networks were calculated from cortical thickness maps; to allow for correspondence across image types, the loadings for these networks were then applied to the volume maps, resulting in only cortical regions being analyzed. Results were similar when networks were calculated directly from volume maps. NMF networks were visualized on the inflated Population-Average, Landmark-, and Surface-Based (PALS) cortical surfaces using the Caret software package.

After identifying networks, we conducted analyses to identify associations between brain structure and dimensions of psychopathology. Given that structural maturation is a nonlinear process, we modeled both linear and nonlinear age effects using penalized splines within generalized additive models using the R package voxel, which relies on the mgcv package in R (20). Generalized additive models assess a penalty on nonlinearity to avoid overfitting and to capture both linear and nonlinear effects in a data-driven fashion. All models used the restricted maximum likelihood framework, which produces estimates of variance and covariance parameters. Based on previous work documenting sex differences in cortical thickness and volume (21), we included sex in the model. In addition, we added mean image quality ratings (described above) as an additional model covariate to ensure that image quality did not drive the observed associations (18). For each network, we examined associations between each dimension of psychopathology and cortical thickness or volume, using the following formula: NMF structural network = spline(age) + sex + image quality rating + anxious-misery + psychosis + behavioral + fear + overall psychopathology.

Interactions between fear and age, fear and sex, and age and sex were evaluated and found to be nonsignificant for cortical thickness. The same interactions for overall psychopathology and anxious-misery were nonsignificant for volume; significant age-by-sex interactions are reported below. To control for multiple testing across networks, we controlled the false discovery rate (q<0.05).

We conducted sensitivity analyses to ensure that our results were robust to methodological choices and were not influenced by confounding variables. First, we repeated the analyses described above using standard anatomical regions. Anatomical regions were delineated using a top-performing, highly accurate multi-atlas labeling tool with joint label fusion implemented in ANTs (see the online supplement). Second, we included maternal education level as an additional covariate and excluded the participants who were taking psychotropic medications at the time of imaging (N=155, 11%). Third, we evaluated whether global differences in total gray matter volume or average cortical thickness were driving the observed effects. Fourth, we used a traditional case-control approach to illustrate the usefulness of the bifactor model in increasing the specificity of the results by taking into account the overlapping variance between disorders.

In addition to univariate associations between each structural covariance network and each dimension of psychopathology, we also conducted a multivariate analysis to examine whether structural networks predicted psychopathology above and beyond demographic characteristics and cognitive performance. To do this, we compared a reduced model including age, sex, and three cognitive performance factors to a full model predicting fear with age, sex, cognitive factors, and all 18 cortical thickness networks using an F test. Similarly, we also performed the same analyses with the volume networks as predictors of overall psychopathology and anxious-misery. Adjusted R² is reported for the models.

For the code and a companion wiki detailing the analytic procedures used in this study, see https://github.com/PennBBL/KaczkurkinPark_BifactorStructure/wiki/ KaczkurkinPark_BifactorStructure. Data from the PNC can be accessed at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000607.v3.p2. The NMF code can be found at https://github.com/asotiras/brainlets.

Structural covariance networks were delineated using NMF at multiple resolutions. The final 18-network solution was chosen on the basis of two considerations. First, we evaluated the gradient of reconstruction error (see Figure S2 in the online supplement), which shows only nominal decrements in error beyond 14 networks. Second, we checked the split-half reliability at this resolution, which revealed an adjusted Rand index of 0.93 for the 18-network solution, suggesting that this solution is highly reproducible. This resolution is also consistent with previous reports (19). Accordingly, the 18-network solution was used for all subsequent analyses. As in previous work using NMF (17, 19), the structural covariance networks identified were highly symmetric bilaterally (Figure 2).

Having identified 18 cortical thickness covariance networks using NMF, we next examined associations with the dimensions of psychopathology summarized by the bifactor analysis. Results revealed that increased fear was associated with reduced cortical thickness in 13 networks after false discovery rate correction (Figure 3; see also Table S1 in the online supplement), with the relationship between fear and cortical thickness showing small effect sizes (partial r values ≤−0.12). Regions affected included the posterior cingulate and temporal-parietal junction (network 16, see Figure 3A), the anterior and subgenual cingulate cortex, and the anterior insula (network 13, see Figure 3C). Furthermore, there were relatively widespread associations across the temporal, orbitofrontal, and occipital cortex. Notably, the association between fear and cortical thickness was specific, and cortical thickness was only weakly associated with other dimensions of psychopathology. The anxious-misery and behavioral dimensions were not associated with cortical thickness in any network, whereas psychosis and overall psychopathology were associated with diminished thickness in only a single network (networks 5 and 16, respectively). However, the results for psychosis and overall psychopathology did not remain significant during sensitivity analyses (described below). Interactions between fear and age, fear and sex, and age and sex were nonsignificant for cortical thickness.

Next, we repeated these analyses, but instead quantified volume (rather than cortical thickness) within each structural covariance network. These analyses revealed that overall psychopathology was associated with reduced volume in all networks, suggesting a global association with gray matter volume (partial r values ≤−0.14; see Table S2 in the online supplement). As illustrated in Figure 4A and 4B, total gray matter volume was negatively associated with overall psychopathology symptoms. In contrast, anxious-misery symptoms were associated with increased volume in most networks (partial r values ≤−0.12; see Table S2). Specifically, total gray matter volume (Figure 4C) was positively associated with anxious-misery symptoms, as were all networks except one (network 7; Figure 4D). Furthermore, fear symptoms were associated with reduced volume in multiple regions (see Table S2), and the behavioral dimension showed reduced volume in the superior parietal cortex (network 6) and the fusiform cortex (network 10). However, associations with the fear and behavioral dimensions did not remain significant during sensitivity analyses (see below). Psychopathology-by-age and psychopathology-by-sex interactions were nonsignificant for volume. Significant age-by-sex interactions were apparent in all volume networks (false discovery rate–corrected p values ≥0.027) except networks 4 and 8.

To ensure that our results were not specific to NMF networks, we evaluated associations using a highly accurate anatomical brain parcellation. Parcellation results aligned with analyses using NMF networks for all dimensions including fear, overall psychopathology, and anxious-misery (see Figure S3 in the online supplement). In addition, we conducted sensitivity analyses to evaluate potentially confounding variables. Nearly all associations between cortical thickness and fear and between volume and overall psychopathology or anxious-misery remained significant after excluding participants taking psychotropic medications (11%) and including maternal education as an additional covariate (see Table S3 in the online supplement). However, the associations between network thickness and fear were no longer significant when average cortical thickness was included as an additional covariate (false discovery rate–corrected p values ≥0.485), suggesting a distributed effect. Similarly, the associations between network volume and either overall psychopathology or anxious-misery were no longer significant when total gray matter volume was added as a covariate, consistent with a global effect (false discovery rate–corrected p values ≥0.229). Finally, using a traditional case-control approach, we found that cortical thickness and volume were reduced in most diagnostic categories (see Table S4 in the online supplement), illustrating the lack of specificity when using this approach.

Next, we tested whether structural covariance networks provided information about psychopathology above and beyond demographic characteristics and cognitive performance. We found a significant difference between a reduced model with only age, sex, and the three cognitive factors and a full model in which fear was predicted by age, sex, cognitive factors, and the 18 cortical thickness networks (F=2.32, df=1367, 1385, p=0.001). The correlation (r) between the actual fear scores and the predicted fear scores in the full model (N=1,389) was 0.28 (p<0.001) (see Figure S4A in the online supplement). However, while the proportion of variance in fear explained by the predictors improved in the full model, it was still relatively modest (adjusted R²=0.06). For volume networks, there was a significant difference between the reduced and full models for overall psychopathology (F=1.78, df=1367, 1385, p=0.023) but not anxious-misery (p=0.218). The correlation between the actual and predicted overall psychopathology scores in the full model (N=1,389) was 0.35 (p<0.001) (see Figure S4B in the online supplement). Compared with the reduced model, the proportion of variance in overall psychopathology explained by the predictors also showed a small improvement in the full model (adjusted R²=0.11).

We found that cortical thickness was reduced in association with fear symptoms in most networks. This result is consistent with case-control studies in adults showing cortical thinning in posttraumatic stress disorder, specific phobia, and social anxiety disorder (4, 22, 23). Our results are also broadly convergent with fear and anxiety networks identified using task-based functional MRI including the salience and ventral attention networks, which are critical for processing of emotionally salient information and attention bias to threat (24). These networks include regions that were affected in this study, including the anterior insula, anterior cingulate cortex, temporo-parietal junction, and ventrolateral prefrontal cortex. The anterior cingulate has been associated with self-regulation, including emotional processing (25), salience processing (26), and attention (27), and has demonstrated functional connectivity with the insula (25, 26), which is implicated in fear. However, our results also extend beyond these salience and emotion regulation networks, suggesting a more widespread effect. While our results are generally consistent with previous case-control research, they build on this literature by showing that global reduced cortical thickness is associated with the spectrum of fear symptoms across disorders in a community-based developmental sample with substantial comorbidity.

It is increasingly recognized that psychopathology exists on a continuum. In a dimensional framework, individuals may be characterized by a profile of symptoms that span categorical boundaries (1, 3). Critically, the bifactor model captures comorbidity and individual variation in the overall level of psychopathology through the general p factor, which may contribute to the nonspecificity of biomarkers found across disorders (2). Recent research has demonstrated that p is linked to individual variation in cognition (28), executive function (29), cerebral blood flow (30), and genomics (31). In this study, we add to this growing literature by demonstrating that higher levels of p were associated with global reductions in gray matter volume in youths. This finding is consistent with work by Goodkind et al. reporting common gray matter volume loss across multiple psychiatric disorders in adults (8) and with previous work showing reduced volume associated with p (13, 14). In contrast to the reduction of gray matter volume seen with higher levels of overall psychopathology, higher levels of anxious-misery symptoms were associated with significant increases in global gray matter volume. Although the magnitudes of these effects were small, they are consistent with effect sizes reported for other variables associated with brain structure, such as IQ (32), and for other forms of psychopathology (5).

Associations with specific model dimensions, such as fear or anxious-misery, must be understood within the context of the bifactor model. Notably, scores from the bifactor model are uncorrelated and represent the burden of a specific dimension while accounting for the overall burden of general psychopathology. In contrast, dimensions derived from traditional factor-analytic models can be highly correlated with each other, which reduces the specificity of observed associations (1, 3). As in any factor model, an individual receives a score from each dimension of the bifactor model, and groups can be compared on these scores. For example, while patients with depression, on average, have high scores on the anxious-misery dimension, because of prominent comorbidity, they also tend to have high levels of overall psychopathology. Accordingly, even though the specific anxious-misery subfactor is associated with greater gray matter volume, this effect may be overwhelmed by the countervailing impact of the overall psychopathology dimension. Indeed, such an example illustrates the advantages of a hierarchical dimensional framework for parsing heterogeneous categorical clinical diagnoses. In contrast, case-control analyses in each diagnostic category yielded nonspecific results, further illustrating the advantage of this approach. Given that traditional diagnostic categories do not take into account the substantial heterogeneity within and comorbidity among disorders, reliance on categorical diagnoses may impede the development of clinically useful neurobiological markers (33–35).

Whereas the majority of translational psychiatric imaging studies have considered adults, in this study we document associations with dimensions of psychopathology in a community-based sample of youths. The PNC is a quasi-epidemiological sample; however, the clinical screening diagnoses were globally consistent with population rates (36). It is important to interpret our results in the context of normative cortical maturation, where reductions in cortical thickness and volume due to both myelination and pruning follow distinct trajectories throughout development (21, 37). Reduced volume or thickness associated with psychopathology is consistent with at least three potential developmental aberrations. First, our findings could represent a structural abnormality that is present from early life (or even in utero) and fixed throughout development. Such a deficit could potentially be linked to maternal infections during pregnancy or obstetric complications at birth (38). Second, these results could alternatively stem from a flattened trajectory of cortical expansion in early childhood, with reduced peak cortical volume and thickness. Finally, these results may be consistent with an accelerated course (or earlier onset) of the normative process of cortical thinning and gray matter volume loss. Intriguingly, recent evidence from multiple lines of research, including epigenetics and translational neuroimaging, suggests that childhood adversity may accelerate the process of cortical development (39). To disambiguate these possibilities, it will be necessary to follow large samples of youths longitudinally from early in life and to acquire detailed data regarding pregnancy and the childhood environment.