Abstract
Childhood socioeconomic status (SES) impacts cognitive development and mental health, but its association with human structural brain development is not yet well characterized. Here, we analyzed 1243 longitudinally acquired structural MRI scans from 623 youth (299 female/324 male) to investigate the relation between SES and cortical and subcortical morphology between ages 5 and 25 years. We found positive associations between SES and total volumes of the brain, cortical sheet, and four separate subcortical structures. These associations were stable between ages 5 and 25. Surface-based shape analysis revealed that higher SES is associated with areal expansion of lateral prefrontal, anterior cingulate, lateral temporal, and superior parietal cortices and ventrolateral thalamic, and medial amygdalo-hippocampal subregions. Meta-analyses of functional imaging data indicate that cortical correlates of SES are centered on brain systems subserving sensorimotor functions, language, memory, and emotional processing. We further show that anatomical variation within a subset of these cortical regions partially mediates the positive association between SES and IQ. Finally, we identify neuroanatomical correlates of SES that exist above and beyond accompanying variation in IQ. Although SES is clearly a complex construct that likely relates to development through diverse, nondeterministic processes, our findings elucidate potential neuroanatomical mediators of the association between SES and cognitive outcomes.
SIGNIFICANCE STATEMENT Childhood socioeconomic status (SES) has been associated with developmental disparities in mental health, cognitive ability, and academic achievement, but efforts to understand underlying SES–brain relationships are ongoing. Here, we leverage a unique developmental neuroimaging dataset to longitudinally map the associations between SES and regional brain anatomy at high spatiotemporal resolution. We find widespread associations between SES and global cortical and subcortical volumes and surface area and localize these correlations to a distributed set of cortical, thalamic, and amygdalo-hippocampal subregions. Anatomical variation within a subset of these regions partially mediates the positive relationship between SES and IQ. Our findings help to localize cortical and subcortical systems that represent candidate biological substrates for the known relationships between SES and cognition.
Introduction
Early brain development occurs within the context of each child's experiences and environment, which vary significantly as a function of socioeconomic status (SES). Childhood SES is typically measured by factors including parental income, education, and occupation (Hollingshead, 1957; McLoyd, 1998) and has been associated with disparate outcomes in mental health, cognitive development, and academic achievement (Brooks-Gunn and Duncan, 1997; Sirin, 2005; Noble et al., 2007; Reiss, 2013). These associations are thought to arise through diverse causal pathways, including: (1) direct SES-linked effects on cognitive and health outcomes (Ritsher et al., 2001; Kendler et al., 2015), (2) the capacity of mental health or cognitive challenges to negatively impact SES (Tiikkaja et al., 2016), and (3) factors that simultaneously increase risk for lowered SES and cognitive difficulties (Trzaskowski et al., 2014; Hill et al., 2016).
The robust epidemiological data connecting childhood SES to behavioral and cognitive development have motivated recent neuroimaging studies aimed at mapping SES effects on brain anatomy. SES has been positively associated with total gray matter volume and less consistently with white matter volume (Luby et al., 2013; Hair et al., 2015; Gianaros et al., 2017), as well as with volume in regions of a priori interest, including the prefrontal cortex and hippocampus (Jednoróg et al., 2012; Noble et al., 2012; Hanson et al., 2013; Luby et al., 2013; Hair et al., 2015; Holz et al., 2015). Furthermore, a recent landmark study (Noble et al., 2015) mapped cross-sectional associations between socioeconomic factors and cortical surface area. These studies suggest that observed associations between childhood SES and neurodevelopmental outcomes may be centered on specific brain systems and raise as yet unresolved questions regarding when SES–neuroanatomy associations are established, whether they are predominantly cortical or subcortical, and how correlations between SES and neuroanatomy relate to accompanying variation in cognitive ability (Noble et al., 2007; Brito and Noble, 2014).
Here, we sought to build on the current understanding of SES and brain development in five key directions. First, to determine whether SES associations with neuroanatomy are developmentally dynamic or stable in this age range (Giedd et al., 1999), we model SES–brain associations in a longitudinal sample of 1243 structural neuroimaging scans from 623 healthy individuals ages 5–25. We combine two complementary approaches to analysis of longitudinal data: mixed models (Pinheiro et al., 2016) and direct analysis of intra-individual estimates of anatomical change computed from repeat scans. Second, we separately model SES associations with regional cortical thickness (CT) and surface area (SA), two biologically distinct phenotypes that together determine cortical volume (Raznahan et al., 2011). Findings remain sparse and mixed regarding the relative strength and spatial distribution of SES relations with SA and CT (Lawson et al., 2013; Noble et al., 2015). Third, we extend our analyses to assess the relation between SES and the anatomy of five major noncortical structures (henceforth “subcortical” structures) using multi-atlas methods that provide both bulk volume estimates and spatially fine-grained measures of shape. Simultaneous examination of SES associations with cortical and subcortical anatomy is critical given evidence that these brain systems function (Redgrave et al., 2010), develop (Raznahan et al., 2014) and connect (Draganski et al., 2008) in a topographically organized manner. Fourth, we formally characterize functional associations of those cortical regions that correlate with SES using the Neurosynth platform for meta-analysis of neuroimaging data (Yarkoni et al., 2011). Finally, given previously noted positive associations between SES and cognitive performance and the potential of brain anatomy to vary as a function of cognitive ability (Walhovd et al., 2016), we probe the complex associations between SES, neuroanatomy, and cognition through two methods that have been used separately in previous literature: (1) including IQ as a covariate in models of SES effects on neuroanatomy (Noble et al., 2012; Lawson et al., 2013) and (2) assessing whether structural phenotypes correlated with SES mediate the relation between SES and IQ (Hair et al., 2015; Noble et al., 2015).
Materials and Methods
Participants.
This study includes a longitudinal sample of 1243 structural magnetic resonance imaging (sMRI) brain scans from 623 healthy children and adolescents between 5 and 25 years old (Table 1). Participants were recruited through local advertisement for a study of typical brain development conducted at the National Institute of Mental Health Intramural Research Program between 1990 and 2010. Participants were screened and excluded on the basis of a history of mental health treatment, use of psychiatric medication, enrollment in special services at school, or diagnosis of any medical condition known to affect the nervous system. The research protocol was approved by the institutional review board at the National Institute of Mental Health and written informed consent or assent was obtained from all children who participated in the study, as well as consent from their parents if the child was under the age of 18.
Socioeconomic status.
Childhood socioeconomic status (SES) was quantified using the Amherst modification of the Hollingshead two-factor index (Hollingshead, 1957; Watt, 1976). Parental education and occupation were each coded on a seven-point scale and these two values were used to derive and record a single SES score that was used for analyses. When education and occupation were reported for two parents, the highest SES was used. In the resulting index, children from the most advantaged families receive the lowest Hollingshead score of 20, whereas children from the most disadvantaged families receive the highest Hollingshead score of 134. For ease of interpretation, we refer to SES variation using conventional directionality such that “higher SES” refers to a lower Hollingshead score. Accordingly, reported positive associations with SES in the manuscript indicate variables that increase in value with greater SES and are thus negatively correlated with Hollingshead score. In graphical representation of such associations, x-axis scales for SES are reversed such that right-most values index the highest SES levels (i.e., lowest Hollingshead scores).
General cognitive ability.
Full-scale IQ was estimated for each child in the sample using an age-appropriate Wechsler scale. The majority (N = 562) of children received the Wechsler Abbreviated Scale of Intelligence (Wechsler, 1999) and other scales used include the WAIS-R, WISC-R, and WISC-III and the WPPSI and WPPSI-III. Of the 344 individuals with repeat assessments, only 60 had IQ measures within 1 year of at least 2 scans. We therefore used the most recent complete IQ measurement for each individual.
Image acquisition and processing.
All sMRI scans were T1-weighted images collected on the same 1.5 T General Electric SIGNA scanner with contiguous 1.5 mm axial slices using a 3D spoiled-gradient recalled-echo sequence (echo time = 5 ms; repetition time = 24 ms; flip angle = 45°; acquisition matrix = 256 × 192; number of excitations = 1; field of view = 24 cm). All scans passed visual assessment for motion artifact before processing.
Native sMRI scans were submitted to the CIVET 1.1.10 pipeline for automated morphometric analysis (Ad-Dab'bagh et al., 2006). The CIVET pipeline uses a validated neural net approach to voxel classification to calculate gray matter and white matter volume estimates (Zijdenbos et al., 2002; Cocosco et al., 2003) after initial correction of images for radiofrequency intensity nonuniformities (Collins et al., 1994; Sled et al., 1998). Following tissue classification, each image was fitted with two deformable mesh models to identify the inner and outer surface of cortical gray matter, and these surfaces were used to calculate cortical thickness and surface area at 40,962 vertices on each cortical hemisphere as described previously (MacDonald et al., 2000; Raznahan et al., 2012).
Subcortical segmentation and surface extraction were completed automatically using the MAGeT Brain algorithm (Chakravarty et al., 2013; Raznahan et al., 2014). Scans were first registered to the ICBM 152 template and corrected for radiofrequency intensity nonuniformities (Collins et al., 1994; Sled et al., 1998). For the striatum, thalamus, and pallidum, the segmentation atlas was created using a 3D reconstruction of serial histological data warped to an MRI-based template (Chakravarty et al., 2006). The MAGeT pipeline then customized this atlas to 21 randomly selected subjects within the sample. All 1243 scans were then warped to this set of templates, providing a set of 21 candidate subcortical segmentations for each scan. For the hippocampus and amygdala, five reference atlases were generated from high-resolution and high-contrast T1- and T2-weighted images from three males and two females using a 3 T scanner with final supersampled isotropic voxel dimensions of 0.3 mm (Wood, 2011). The MAGeT pipeline again created automated segmentation atlases for 21 randomly selected subjects, resulting in 105 possible segmentations (five atlases × 21 templates) for each of the 1243 scans in our dataset. Each scan was labeled using the 21 striatum, thalamus, and pallidum segmentations and the 105 hippocampus and amygdala segmentations and the final segmentation was decided upon using a label-voting procedure such that the label occurring most frequently at each voxel was retained. These procedures provided estimates of total bilateral volume for the hippocampus, amygdala, thalamus, striatum, and pallidum in each scan. All scans used for analysis passed visual quality control of these final subcortical segmentations to exclude visible segmentation errors of the five subcortical structures under study.
Surface-based representations of all five subcortical structures were then estimated on their respective atlases using a marching cubes algorithm (Lerch et al., 2008). Next, the nonlinear portions of the 21 transformations mapping each subject to the 21 input templates were concatenated and averaged across the template library to limit noise and increase precision and accuracy. These surface-based representations were warped to fit each template and each surface was warped to match each subject. This procedure yields 21 possible surface representations per subject for the striatum, thalamus, and pallidum and 105 possible surface representations for the hippocampus and amygdala, which were merged by estimating the median coordinate representation at each location. Next, a third of the surface area of each triangle forming the surface representation was assigned to each vertex within the triangle. The surface area at each vertex is the sum of all such assignments from all connected triangles. Finally, surface area values were blurred with a surface-based diffusion-smoothing kernel (5 mm for the amygdala, hippocampus, striatum, and thalamus and 3 mm for the pallidum). This processing stream generated surface area values for a total of 26,401 vertices across the five subcortical structures in each scan.
Experimental design and statistical analyses.
The effect of SES on each anatomical metric of interest was modeled using a linear mixed-effects model with sex and centered age as fixed effects covariates and each individual's ID (to account for multiple longitudinal scans per individual) and family ID (to account for the presence of dizygotic twin pairs and siblings in the sample) as nested random effects. All linear mixed-effects models discussed below were fitted using the R package nlme version 3.1–128 (Pinheiro et al., 2016). Sex was treated as a categorical variable coded as “M” or “F.” SES and age variables were both centered. Individual ID and family ID were specified as nested random effects and were fitted with random intercepts.
The decision to present core results with SES, sex, and age as main effects was made after first ruling out the presence of extensive interactions between these variables (i.e., SES × sex or SES × age interactions) in predicting structure or vertex-level anatomical variation effect on each anatomical variable. For most structures, there was no significant interaction between fixed-effects variables; the few cases in which interactions were found will be noted in the text. Otherwise, anatomy for the ith family's jth individual's kth time point was modeled as follows: Dependent volumetric variables of interest included total brain volume (TBV; the sum of intracranial gray and white matter volume), total intracranial gray matter volume (GMV) and total intracranial white matter volume (WMV), total bilateral volume of the cortex (CV), total cortical surface area (SA) and mean cortical thickness (CT), and total bilateral volume of each subcortical structure (i.e., “bulk” hippocampal, amygdalar, thalamic, striatal, and pallidal volumes). To aid comparison of the associations between SES and these diverse volumetric indices, we also estimated the standardized effect for each SES–volume association by re-running model [A] above after centering and scaling all variables so that the resulting β1 coefficient would index the SD shift in volume with one SD increase in SES (i.e., a decrease of ∼18 Hollingshead points).
Vertex-level anatomical variables of interest included cortical thickness and surface area at each of 80,962 cortical vertices and SA at each of 26,401 subcortical vertices (hippocampus, 1152 left/1215 right; amygdala, 1473 left/1405 right; thalamus, 3016 left/3108 right; striatum, 6450 left/6178 right; pallidum, 1266 left/1138 right). Vertex-specific β1 coefficients were visualized on the corresponding cortical or subcortical surface after applying a false discovery rate (FDR) correction for multiple comparisons. FDR corrections were calculated separately across the left and right cortical hemispheres and the left and right subcortical structures with q (the expected proportion of false rejections of the null hypothesis) set to 0.05.
Finally, we probed the relation among SES, anatomical metrics, and cognition in two ways. First, for all structures that showed a significant main effect of SES in model [A], we tested the robustness of this effect to the inclusion of full-scale IQ as a covariate in the linear mixed-effects models. Separate main effects of SES and IQ were included after first ruling out the presence of a significant SES by IQ interactive effect. For these analyses, anatomy for the ith family's jth individual's kth time point was modeled as follows: Linear age terms were used in models [A] and [B] after verifying that this simple parameterization of age yielded identical SES findings to models run with age parameterized as a nonlinear spline (which could allow for nonlinear age effects). Additionally, models [A] and [B] were re-run for each anatomical metric within the subset of participants (n = 534) who self-endorsed the federal race category of “white” and all main effects held.
To complement our mixed-effects models and to further probe the nature of SES–anatomy associations, we directly tested the relationship between SES and intra-individual change in brain anatomy for the 344 individuals with more than one MRI scan. For each anatomical variable, we calculated an annual rate of change for each individual. Anatomical change was then modeled using a linear model with main effects of SES, mean scan age, and sex and the interaction between SES and mean scan age. The association between annual rate of change and SES was investigated for each of the anatomical variables considered in the mixed models above: TBV, GMV, WMV, CV, total SA, mean CT, five subcortical volumes, and vertexwise CT, cortical SA, and subcortical SA.
Finally, we conducted mediation analyses to investigate the extent to which the global and vertexwise structural associations with SES mediate the association between SES and IQ. For the mediation tests with a single global mediator, we used the Mediation package in R (Tingley et al., 2014), with SES as the independent variable (treatment) and full-scale IQ as the dependent variable. The mediate() function estimates the average causal mediation effect (ACME) and the average direct effect (ADE), which together sum to the total effect of the treatment (i.e., SES) on the outcome (i.e., IQ). The proportion mediated, which we report, represents the size of the ACME relative to the total effect. Within the map of cortical vertices that showed a significant association with SES, we tested the mediating role of each vertex using the MultiMed package in R (Boca et al., 2014), which implements a permutation approach with joint correction to test multiple mediators simultaneously.
Interpretation of anatomical results.
To systematically investigate the functional implications of SES effects, we submitted the MNI coordinates of peak SES effects on surface area and cortical thickness to Neurosynth (Yarkoni et al., 2011), an online platform that extracts and synthesizes brain activation patterns and psychological terms across >11,000 functional neuroimaging publications. It can be used to generate correlations between the meta-analytic coactivation map for a given point of interest in the brain and each of the terms in the Neurosynth database. Using Neurosynth, we identified cognitive and psychological terms that frequently cooccur in the literature with functional activations similar to the observed pattern of SES effects on cortical morphology.
Results
Participant characteristics
Participant characteristics are detailed in Table 1. Hollingshead SES in our sample ranged from a high of 20 to a low of 95 after three scans from individuals with exceptionally low SES scores (Hollingshead SES = 115, |z-score| > 4) were removed. The complete Hollingshead scale ranges from 20 to 134 and approximately one-quarter of our 623 participants received the highest SES score of 20 (indicating that one or both of their parents received a graduate professional degree and employment in the highest occupation category). Therefore, our sample's distribution of Hollingshead scale scores does not capture the most severely socioeconomically disadvantaged children and it is relatively enriched for children with highly educated and highly employed parents. Throughout the text below, SES is referred to using conventional directionality such that “higher SES” refers to a lower Hollingshead score, and reported positive associations with SES are thus negatively correlated with Hollingshead score. Within the cross-sectional sample of 623 individuals, SES and IQ were significantly, positively correlated, r(621) = 0.31, p < 0.001.
SES and measures of brain volume
We found a strong positive association between SES and total brain volume (β1A = 1217.0, p < 0.001). A positive association with SES was also seen for total white and gray matter volumes, as well as for total volume of the cortical sheet (Table 2). Separate examination of the two determinants of cortical volume identified a strong positive association between SES and total cortical surface area (β1A = 144.1, p < 0.001) and a weaker positive association between SES and mean cortical thickness (β1A = 0.00086, p = 0.01). Greater SES was also significantly associated with greater bilateral volume of all subcortical volumes examined except the pallidum (Table 2). Analysis of standardized effect sizes for SES associations with regional brain volumes revealed that SES was most strongly related to total brain volume and total surface area such that an SD increase in SES was associated with a 0.17 SD increase in each of these metrics (Fig. 1A). Among subcortical structures, SES had the greatest effect on thalamic volume such that a SD increase in SES was associated with a 0.15 SD increase in thalamic volume (Fig. 1A).
These above associations between SES and brain anatomy were all stable between ages 5 and 25 and comparable between males and females, with two exceptions: (1) the association between SES and total cortical surface area was modified by sex such that there was a stronger SES effect on surface area for males than females (SES × sex interaction effect: βinteraction = 122.2, p = 0.019; SES main effects: β1A,MALES = 198.7, p < 0.001/β1A,FEMALES = 77.6, p = 0.066) and (2) the association between SES and hippocampal volume was modified by age such that the effect of SES on hippocampal volume grows with age (SES × age interaction effect: βinteraction = 0.19, p = 0.003).
SES and cortical morphology
After establishing that SES was associated with total surface area and mean cortical thickness (Table 2), we tested for regional specificity of these associations through vertex-level analysis of SES associations at 80,962 points (vertices) across the cortical sheet. Vertex-level analyses established that the robust association between greater SES and greater total cortical SA was underpinned by statistically significant positive associations between SES and regional SA within a distributed set of largely bilateral cortical areas including the lateral prefrontal, anterior cingulate, and lateral temporal and superior parietal lobule regions (Fig. 1B). Associations between SES and vertex-level SA were not significantly modified by age or sex after correction for multiple comparisons.
Associations between SES and regional cortical thickness were much more localized than observed for regional SA. Specifically, we identified a single locus of significant positive association between SES and CT in the right supramarginal gyrus.
SES and subcortical morphology
We investigated the spatial specificity of observed associations between SES and bilateral subcortical volumes by modeling surface area at each vertex across the surface of the hippocampus, amygdala, thalamus, and striatum (Fig. 1C). Hippocampal effects were concentrated bilaterally on the medial surface and neighboring medial surfaces of the right amygdala also showed a significant positive association with SES. The association between SES and thalamic surface area was localized primarily to the ventral posterior and ventral lateral thalamus. No other subcortical structures showed statistically significant shape associations with SES after correcting for multiple comparisons across vertices.
SES and intra-individual anatomical change
To complement our mixed-model analyses and directly model intra-individual change, we further investigated the relationships between SES and annual rate of anatomical change of each of the cortical and subcortical metrics discussed above. Among the 344 individuals with more than one MRI scan, there was no significant relationship between SES and change in any of the global anatomical metrics (TTV, GMV, WMV, CV, total SA, mean CT, or any subcortical volume). These results are fully consistent with the previously noted lack of significant SES × age interactive effect on any global anatomical metric, with the sole exception of findings for hippocampal volume. Results from analyses of intra-individual change and from mixed-model analysis were also largely consistent for vertex-level measures of cortical surface area and thickness: both analytic approaches suggested that age does not modulate SES effects on thickness or surface area across most of the cortical sheet or on area across subcortical surfaces. We observed two spatially limited exceptions to this convergence between analytic techniques in that analysis of intra-individual change uniquely identified a statistically significant negative relationships between SES and cortical thickness change in the left middle temporal gyrus and the left superior parietal lobule after FDR correction and a statistically significant, negative relationship between SES and surface area change in a single locus on the left dorsal thalamus. A figure showing results of the vertexwise change analyses is available on figshare (McDermott, 2018).
Separating main effects of SES and IQ on anatomy
Because of the noted strong association between SES and IQ (Fig. 2A), we next re-ran linear mixed-effects models with full-scale IQ as a covariate to parse main effects of SES and cognition on neuroanatomy. We observed positive associations between IQ and all global measures of cortical and subcortical anatomy (Table 2). The main effect of SES on each global metric was reduced in magnitude by the addition of IQ to the model; for mean cortical thickness, bulk amygdalar volume, and bulk striatal volume, the main effect of SES was no longer significant (Table 2), although SES continued to show a significant positive association with TBV, GMV, WMV, CV, total cortical SA, and bilateral hippocampal and thalamic volumes after controlling for IQ. There were no significant SES × IQ interactive effects on any global metrics.
We next probed the spatial patterning of IQ versus SES main effects on cortical surface area. After FDR correction for multiple comparisons, childhood SES and IQ were both positively associated with regional cortical SA (Fig. 2B). Specifically, when including IQ as a covariate, the main effect of SES was restricted to the bilateral superior parietal, right orbitofrontal, left inferior temporal, and bilateral medial prefrontal cortices. The main effect of IQ, in contrast, was localized to the left inferior and middle temporal, left inferior parietal, and left medial frontal regions. Strikingly, the regional maps of separate SES and IQ effects on cortical surface area together (Fig. 2B) resemble the regional effects of SES when IQ is not included as a covariate (Fig. 1B). In other words, the main effect of SES on cortical surface area appears to be separable into effects related to the strong IQ differences across SES localized to the left-lateralized perisylvian and medial prefrontal cortices and SES–SA associations that are independent from associated variation in general cognitive ability localized to bilateral parietal and frontal cortices.
Structural measures mediate the association between SES and IQ
In addition to testing the effects of IQ as a covariate in SES–brain associations, we also probed the relation among SES, cognition, and anatomy by investigating whether any of the structural brain measures found to be significantly associated with SES might mediate the association between SES and IQ in our cross-sectional dataset of 623 individuals. Using the Mediate package in R (Tingley et al., 2014), we estimated the proportion of the total relation between SES and IQ that was accounted for by each anatomical mediator. Among global cortical measures, total brain volume (proportion mediated = 0.148, p < 0.001), gray matter volume (proportion mediated = 0.145, p < 0.001), white matter volume (proportion mediated = 0.077, p < 0.001), cortical volume (proportion mediated = 0.146, p < 0.001), total surface area (proportion mediated = 0.114, p < 0.001), and mean cortical thickness (proportion mediated = 0.049, p = 0.01) each partially mediated the relation between SES and IQ. Subcortical volumes also partially mediated the association (hippocampus, proportion mediated = 0.052, p < 0.001; amygdala, proportion mediated = 0.034; p = 0.02; thalamus, proportion mediated = 0.087; p < 0.001; striatum, proportion mediated = 0.033, p = 0.01).
Finally, we extended the mediation analysis to investigate whether there was a regional specificity to the partial mediation of SES and IQ by cortical and subcortical surface area. To do so, we submitted all of the vertices in each hemisphere that showed a significant association with SES to a multiple mediation analysis using the MultiMed package in R (Boca et al., 2014). Three cortical regions, all in the left hemisphere, showed significant mediation of the SES association with IQ: the middle temporal gyrus, supramarginal gyrus, and anterior cingulate cortex (Fig. 2C). Notably, these largely overlapped with the cortical vertices that showed shared main effects of SES and IQ on surface area (regions in purple, Fig. 2B).
Functional implications of anatomical results
Coordinates of peak regional SES associations with cortical surface area and cortical thickness were submitted to Neurosynth (Yarkoni et al., 2011), a meta-analytic repository of structure-function relationships. The top five cognitive and psychological terms with the highest meta-analytic coactivation correlation coefficient r are displayed in Table 3 after redundant terms and anatomical terms were removed. This platform provided evidence that the spatial effects of SES are localized to regions preferentially associated with sensorimotor functions, as well as language, memory, and emotional processing, which have each been shown to exhibit SES differences (Farah et al., 2006; Noble et al., 2007; Kim P et al., 2013).
Discussion
Here, with a large, longitudinal, single-site neuroimaging sample, we both replicate and extend findings on the relation between childhood SES and structural brain anatomy in a number of key directions.
First, consistent with previously replicated findings, we find that SES is positively associated with global brain measures including gray matter volume (Hanson et al., 2013; Luby et al., 2013; Mackey et al., 2015; Gianaros et al., 2017), cortical surface area (Noble et al., 2015; Gianaros et al., 2017), cortical thickness (Mackey et al., 2015; Noble et al., 2015), and hippocampal volume (Noble et al., 2012, 2015; Luby et al., 2013; Hair et al., 2015). Although previous findings regarding white matter volume and amygdala volume have been mixed (Jednoróg et al., 2012; Luby et al., 2013; Noble et al., 2015; Gianaros et al., 2017), we identify significant positive SES associations with both WMV and amygdala volume. We extend analyses of SES associations with subcortical anatomy to provide the first evidence that greater childhood SES is associated with larger bilateral volumes of the thalamus and striatum. The thalamic finding is especially notable given reported associations between thalamic volume and cognitive performance (Van Der Werf et al., 2001; Xie et al., 2012). Comparative analyses of SES associations using standardized effect sizes ranks total cortical SA and bilateral thalamic volume as the morphometric properties of the cortex and subcortex, respectively, which show the largest effect-size relations with SES.
It is of particular note that, with the exception of hippocampal volume, we found the above noted SES associations with global cortical and subcortical anatomical metrics to be fixed between ages 5 and 25. This conclusion was supported by two complementary analytic approaches: mixed-model testing for significant SES × age interactive effects and direct analysis of SES effects on intra-individual rates of anatomical change. Other cross-sectional studies with age ranges similar to the population examined here have likewise reported limited or no SES × age interactions (Lawson et al., 2013; Noble et al., 2015); however, a longitudinal study of children under the age of 4 found that SES disparities in GMV grow with age (Hanson et al., 2013). These findings suggest that SES and age may interact differently across developmental time, with developmentally dynamic associations early in life, which subsequently stabilize. Future longitudinal studies bridging the gap between infancy and adolescence may help to clarify the temporal dynamics of SES associations with brain development.
Second, our study advances understanding of the neuroanatomical correlates of SES by using surface-based algorithms to parse associations between SES and anatomy of the cortical sheet. We observe SES associations with cortical surface area in a number of largely bilateral regions including lateral prefrontal, anterior cingulate, lateral temporal, and superior parietal lobule regions. These SES–surface area associations largely correspond to the map of associations between parental education and cortical surface area presented by Noble et al. (2015). Functional interpretation with Neurosynth (Yarkoni et al., 2011) suggests that the cortical regions that we found to show morphological associations with SES are preferentially involved in networks that underlie sensorimotor functions, as well as language, memory, and emotional regulation.
Third, we provide the first systematic examination of SES associations with subcortical shape, which indicates that variations in childhood SES are associated with focal differences in hippocampal, amygdalar, and thalamic anatomy (Fig. 1C). Histological and functional connectivity neuroimaging studies suggest that the ventral lateral and ventral posterior thalamic nuclei are preferentially connected to frontoparietal cortical systems that subserve primarily sensorimotor functions (Jones, 1985; Kim D-J et al., 2013). Strikingly, we also detect strong associations between SES and surface area within these cortical targets of ventral thalamic nuclei, suggesting that anatomical correlations of childhood SES variation may be organized by the topography of corticosubcortical connectivity.
Finally, we demonstrate convergent anatomical correlates of the strong relation between SES and IQ using two independent methods. Modeling the main effects of both IQ and SES on vertexwise surface area reveals that the map of statistical SES effects can be fractionated into distinct SES and IQ effects. Additionally, we demonstrate that many of the anatomical correlates of SES significantly mediate the relation between SES and IQ. We expand upon previous results that found significant SES–cognition mediation by whole-brain surface area (Noble et al., 2015) by conducting cross-sectional mediation tests at each cortical vertex that had a significant SES–surface area association. Specifically, we find that vertices in the left middle temporal gyrus, supramarginal gyrus, and anterior cingulate, the same three regions that showed a main effect of IQ when controlling for SES, mediate the relation between SES and IQ. These three regions are known to be involved in language (Vannest et al., 2009) and cognitive control (Shackman et al., 2011).
A mediation analysis allows us to test one possible pathway between SES, brain anatomy, and cognition and we show that SES may exert some of its effect on cognition by altering structural brain development, particularly in regions associated with language and learning. However, it is important to note that this pathway represents only one possible set of interactions between childhood environment, anatomy, and cognition. Farah (2017) provides a succinct review of the main processes that may operate to exacerbate neural and psychological SES disparities: the social causation hypothesis suggests that the environmental conditions associated with different levels of SES influence brain structure and function, whereas the social selection hypothesis suggests that genetic factors in parents that both affect their cognition and predispose them to a certain SES are then transmitted to their children. These processes likely both operate in concert and interact to some degree. Because of the inherent observational nature of studies of SES, it is difficult to determine to what extent the anatomical correlates of SES that we report here may reflect the shared genetic effects on SES and brain development versus direct effects of SES on brain development. Nevertheless, our findings help to pinpoint cortical and subcortical systems that represent candidate biological substrates for these diverse causal pathways.
Caveats and future directions
One important limitation of this sample is the composite nature of the Hollingshead two-factor index of SES. The field has lately recognized the importance of differentiating between the effects of separate components of socioeconomic status (Duncan and Magnuson, 2012). A few studies have already begun to identify differential effects of parental education and family income on structural brain development (Lawson et al., 2013; Noble et al., 2015). However, when data collection for this sample began in the 1990s, the Hollingshead was selected for its widespread use and ease of measurement and only the single derived SES score was recorded initially. As future studies are designed, it will be important to report on more nuanced factors that compose SES and may have unique effects on brain development and thus serve as specific targets for intervention.
It is also important to note that the subjects in this sample are not representative of the socioeconomic distribution in the United States. Sample composition is known to affect conclusions about the normative trajectory of brain development (LeWinn et al., 2017). Additionally, some studies have suggested that socioeconomic variables relate most strongly to brain development among the most disadvantaged children (Noble et al., 2015) and the fact that we did not detect such a gradient may be due to our unrepresentative SES distribution. The average IQ of our subjects was also higher than the expected mean of 100 (Table 1); however, IQ is known to increase across generations in the general population, so the high mean IQ may be partially attributable to the use of the same IQ assessment across the multiple decades of this study (Flynn, 1987). Finally, most investigations of SES and brain development to date, including this study, have sampled from Western nations in general and the United States in particular. Further research is necessary to explore how SES associations with brain development play out across other nations and cultures. Although we cannot generalize across the socioeconomic, cognitive, or cultural spectrums, it is notable that we found detectible SES-linked differences in cognitive ability and structural brain development within a typically developing cohort lacking frank socioeconomic deprivation.
Conclusion
Childhood socioeconomic status is a complex construct that influences the physical and psychosocial environment in which a child develops. Here, we demonstrate regionally specific associations between childhood SES and both cortical and subcortical morphology. Our findings inform ongoing efforts to clarify the spatiotemporal patterning of SES-related neuroanatomical variation and its relation to cognitive outcomes such as IQ. Definition of these neuroanatomical associations may ultimately provide candidate biological substrates against which to test potential mechanistic pathways between SES and cognitive and health outcomes. The results presented here do not establish a direct causal pathway between SES and brain development, nor do they indicate that childhood SES exerts a deterministic effect on development. Rather, by resolving neuroanatomical substrates that vary closely with SES, we contribute new biological information to a growing field of multidisciplinary research that ultimately aims to reduce SES variation in health and achievement.
Footnotes
This work was supported by the Intramural Research Program of the National Institutes of Mental Health–National Institutes of Health (1ZIAMH002949-01 under protocol NCT00001246). We thank the participants and families who took part in this study.
The authors declare no competing financial interests.
- Correspondence should be addressed to Armin Raznahan at raznahana{at}mail.nih.gov