The Impact of Early Life Stress on the Genetic Influence on Brain and Cognitive Development in Children
=======================================================================================================

* Hee-Hwan Wang
* Seo-Yoon Moon
* Hyun-Jin Kim
* KaKyeong Kim
* Woo-Young Ahn
* Yoonjung Yoonie Joo
* Jiook Cha

## Abstract

Early life stress (ELS), such as *abuse*, neglect, and maltreatment, is a well-known risk factor for mental illness. However, it is unclear how ELS affects the brain and cognitive development. Identifying specific relationships of ELS with the genetic and brain-related underpinnings of cognitive development may reveal biological mechanisms responsible for the negative impact of ELS and those that lead to individual differences in sensitivity (or resilience) to ELS. In this study, to investigate the interlinked processes of cognitive development, we analyzed the multimodal data of DNA genotypes, brain imaging (MRI), and neuropsychological assessment (NIH Toolbox) outcomes of 4,276 children (ages 9 to 10 years, European ancestry) from the Adolescent Brain Cognitive Development (ABCD) study. We estimated the genetic influence on cognitive capacity using genome-wide polygenic scores (GPSs). Our regression and mediation analyses revealed significant causal relationships for the gene-brain-cognition pathway: Brain structural development significantly mediated the genetic influence on cognitive development (*partial mediation effect = 0*.*016, P**FWE**<0*.*001*). Interestingly, within the triangular relationship, we found a significant moderation effect of *abuse* only on the gene-to-brain pathway (*Index of Moderated Mediation = −0*.*007; 95% CI= −0*.*012 ∼ −0*.*002; P**FWE**<0*.*05*). These findings indicate the negative modulatory effects of ELS on the genetic influence on brain structural development that lead to disadvantageous neurocognitive development in prepubertal children.

Keywords
*   genes-brain-cognition
*   Early Life Stress (ELS)
*   moderated mediation
*   path modeling
*   genome-wide polygenic scores (GPS)

## Introduction

Early life stress (ELS), such as *abuse* and *neglect*, is a well-known major risk factor for maladaptive cognition, behaviors, and psychiatric disorders1–4 with long-term sequelae5. ELS modulates the expression of the genes responsible for stress physiology, emotion regulation, cognitive control, and learning and memory6–9,10. Animal literature shows a causal effect of ELS on cognition and behaviors11–20 and the epi-/genetic mechanisms within the brain21–24. In humans, similar mechanisms may be involved in psychopathology25,26, stress physiology7, or emotion regulation27–29.

Although extensive literature shows associations of ELS with poor cognitive outcomes and brain functional and structural representations1,3,8, exactly how ELS affects cognitive capacity in humans remains unclear. Moreover, the impact of ELS on the development of the brain and cognition in children remains unknown. Addressing these questions will provide much-needed insight into who, when, and where to intervene regarding the negative neurocognitive outcomes of ELS. Considering the well-known impact of ELS on cognition3,10, we hypothesize that ELS is likely to modulate genetic influences on brain development as well as cognitive development3,5.

Cognitive development involves multiple dynamic processes that are regulated by the complex interplay of genetic and environmental factors30–32. Despite its significant impact, however, less is known about the biological mechanisms of how early environmental factors such as ELS, interact with genetic influences and impact cognitive development in children. The recent development of the genome-wide polygenic score approach (GPS) allows individual genetic loadings for cognitive capacity (e.g., accounting for up to 13% of the variance in cognitive outcomes in adults33) to be estimated quantitatively. Using this rigorous computational approach, we aimed to test whether ELS modulates the gene-brain-cognition development process.

This study had two main questions. First, does a gene-brain-cognition pathway exist during developmental periods? If so, does ELS negatively affect the causal pathway? We investigated these questions using multimodal data, including genetics, brain imaging, psychological data, and cognitive testing in large samples of prepubertal children (ages 9 to 10 years) with a wide range of sociodemographic characteristics.

## Methods

### Study Participants

Participants were enrolled in the Adolescent Brain Cognitive Development (ABCD) study, an ongoing longitudinal multisite study of brain development and health of youth in the United States. We used the deidentified neuroimaging, demographic, genetic, and behavioral data from the ABCD 2.0 Data Release34, these data were collected across 21 sites from 2015 to 2018. The participants were 11,875 children from the U.S. who were aged 9-10 years35. After removing participants who had at least one missing value from genetic, neuroimaging, and ELS data, 4,276 participants were included in our analyses.

### Genotype data for Genome-wide polygenic scores (GPS)

Saliva samples were collected from study participants and were genotyped using the Affymetrix NIDA Smokescreen array (733,293 SNPs) at the Rutgers University Cell and DNA Repository. We removed any inferiorly genotyped SNPs with thresholds of (i) genotype call rate<95%, (ii) sample call rate<95%, and (iii) rare variants with minor allele frequency (MAF)<1%. We imputed the data using the Michigan Imputation Server36 based on the 1000 Genome phase 3 reference panel with Eagle ver2.4 phased output37. For the imputed 12,046,090 SNPs, we additionally removed data from any individuals with >5% missing genotypes; extreme heterozygosity (F coefficient > 3 standard deviation from the population mean); and SNPs with <0.4 imputation INFO score, >5% missingness rate, <1% MAF and Hardy-Weinberg equilibrium (p<10−6). Genetic ancestry was determined with the fastSTRUCTURE algorithm38, available from ABCD release 3.0. Our analysis was restricted to individuals with European ancestry to control for population stratification. We estimated both kinship coefficients (K.C.s) and principal components (P.C.s) to control familial relatedness and ancestry admixture using PC-Air39 and PC-Relate40. We selected unrelated samples that were inferred to be more distant than 4th-degree relatives (K.C.>0.022) and removed any genetic outliers that fell significantly outside (>6 S.D. limits) the center in P.C. space. In the rest of this paper, we used genotype data (11,301,999 variants) of 8,523 unrelated European samples after Q.C., and the first ten ancestrally informative P.C.s were used as covariates in all the reported analyses (**Supplementary Figure 1**).

To assess the subject-level genetic propensity for cognitive ability, we computed genome-wide polygenic scores (GPS) for cognitive performance (CP) and educational attainment (EA) using publicly available GWAS summary statistics based on the European population30. The GPSs were constructed using PRSice241, adjusting for the first ten principal components (P.C.s) of the genotype data to control for population stratification.

### Brain Imaging-Anatomical Imaging

T1-weighted (T1w) 3D structural MRI images were acquired in the ABCD study. We processed the images following established protocols42,43: The gradient nonlinearity distortion correction method was performed on structural MRI scans to improve geometric accuracy and image intensity reproducibility44. Based on tissue segmentation and sparse spatial smoothing, intensity nonuniformity was corrected. Then, the data were then resampled with 1 mm isotropic voxels into rigid alignment with an atlas-derived brain.

Cortical surface reconstruction was applied using the following procedures: structural MRI scans were processed using FreeSurfer v6.0 ([https://surfer.nmr.mgh.harvard.edu](https://surfer.nmr.mgh.harvard.edu)) for cortical surface reconstruction45, which includes skull-stripping46, white matter segmentation and initial mesh creation45, correction of topological defects, surface optimization47,48, and nonlinear registration to a spherical surface-based atlas49.

### Brain Imaging-Diffusion Spectrum Imaging

We used the diffusion spectrum images from the ABCD study that were preprocessed using the following protocol50 by the ABCD Data Analysis and Informatics Center (DAIC). Eddy current distortion correction was used with a nonlinear estimation using diffusion gradient orientations and amplitudes to predict the pattern of distortion51. Head motion was corrected by registering images synthesized from tensor fit52. Diffusion gradients were adjusted for head rotation52,53. To identify and replace dark slices due to abrupt head motion, we used robust diffusion tensor estimation54. B0 distortion was corrected with the reversing gradient method55. Gradient nonlinearity distortion correction was applied44. The data were resampled to a standard orientation with an isotropic resolution of 1.5 mm.

To estimate accurate brain imaging phenotypes, we used individual connectome data. We applied MRtrix356 for whole-brain white matter tract estimation and individualized connectome generation. For connectivity metrics, we used streamline counts associated with fiber connection strength57,58 associated with fiber integrity. We decreased noise59, and performed bias correction with the Advanced Normalization Tools (ANTs) pipeline’s N4 algorithm60. To obtain a connectivity index with a white matter pathway61, we performed probabilistic tractography by second-order integration over fiber orientation distributions62, with random seeding across the brain and target streamline counts of 20 million. These initial tractograms were filtered from preliminary tractograms using spherical-deconvolution informed filtering (2:1 ratio). With a final streamline count of 10 million, we generated an 84 × 84 whole-brain connectome matrix for each participant using the T1-based parcellation and segmentation in FreeSurfer. This pipeline ensured that individual participants’ connectomes were restricted to their own neuroanatomy. We carried out the computation on the supercomputers at Argonne Leadership Computing Facility Theta and Texas Advanced Computing Center Stampede2.

### Dimensionality reduction of brain data

To extract brain morphometric representations correlated with genetic influence on intelligence, we first performed generalized linear model analysis using cognitive capacity GPSs as independent variables and neuroimaging data as dependent variables adjusted to covariates (e.g., age, sex, maternal education, income, BMI, study site, and marital status). Among three different types of neuroimaging data (structural MRI, DTI-count, and DTI-FA), only brain morphometric features (i.e., structural MRI) showed significant associations with cognitive capacity GPSs (*P**FWE* *< 0*.*05*; **Figure 3, Supplementary Figure 2, Supplementary Table 1**)63. We then performed principal component analysis to extract representations from brain morphological features associated with CP GPS or EA GPS (R version 3.4.1). We used the first principal components as representations of the brain.

### NIH Toolbox Cognition Battery data

NIH Toolbox Cognitive Function Tests64 were used to assess various levels of general cognitive ability of the participants. For each child, summary scores from the NIH Toolbox Cognition Battery were provided, including the crystallized intelligence composite score (mainly measuring prior learning and past experiences about language), the fluid intelligence composite score (mainly measuring abstract reasoning and learning ability in novel situations), and total intelligence composite score, which is the combination of both crystallized and fluid intelligence composite scores65,66. The crystallized intelligence composite score is a composite of the *Picture Vocabulary Test* and the *Oral Reading Recognition Test* outcomes, and the fluid intelligence composite score is a composite of the *Dimensional Change Card Sort Test*, the *Flanker Inhibitory Contro*l and *Attention Test*, the *Picture Sequence Memory Test*, the *List Sorting Working Memory Test*, and the *Pattern Comparison Processing Speed Test* outcomes. Total intelligence composite score is a representation of general intelligence and is an aggregation of all the tests67.

### Early life stress

We derived ELS measures based on child exposure domains in the ABCD study68 (**Supplementary Table 2**). ELS measures were divided into three main categories: *household challenges, neglect*, and *abuse*. Subcategories of ELS data were as follows: *parental separation or divorce, criminal household member, household substance abuse, mental illness in household, mother treated violently* in *household challenges, emotional neglect, physical neglect* in *neglect*; and *physical abuse* and *sexual abuse* in *abuse*. We extracted the items of each subscale from the following various measurement tools (participant-or parent-reported): ABCD Family Environment Scale-Family Conflict Subscale Modified from PhenX, ABCD Diagnostic Interview for DSM-5 Traumatic Events, ABCD Family History Assessment, ABCD Parent Demographics Survey, ABCD Children’s Report of Parental Behavioral Inventory, and ABCD Parental Monitoring Survey. We averaged the measurements for each subcategory and transformed them into z scores. The higher the score, the more stressful experiences children had.

### Statistical Analysis

We performed a generalized linear model for three types of relationships — GPS-brain, GPS-intelligence composite scores, and brain-intelligence composite scores. We included the following covariates into those models: age, sex, maternal education, parental income, BMI, study site, and marital status. Family ID was not included because our analysis was performed in the dataset of unrelated individuals after excluding any relatives fourth-degree or closer during the genetic QC process. We used familywise error (i.e., Bonferroni correction) in each type of relationship for multiple comparison corrections.

For path modeling, we first tested an initial mediation model (**Figure 1-a**) to check whether the key relationship of the GPS-brain-cognition pathway was significant. Before examining the role of ELS in the GPS-brain-cognition pathway, we tested whether ELS directly affects cognition and whether ELS moderates the gene-cognition pathway (**Figure 1-b**). After we determined the potential moderation effect of ELS on the GPS-brain-cognition pathway, we evaluated the first-stage moderated mediation models and the second-stage moderated mediation models, following the framework of moderated mediation analysis69–71. These models assessed how ELS and its interaction with genetics or the brain, would affect the triangular pathway. We repeated the analyses with different subtypes of ELS and compared their effects.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/30/2021.12.27.21268445/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/F1)

Figure 1. Schematic overview of the tested path models.
(a) The baseline mediation model has the brain factor as a mediator, (b) The moderation model has ELS as a moderator, (c) The first-stage moderated mediation model for testing the moderating role of ELS on the relationship between the brain and genomic factors, (c) The second-stage moderated mediation model for testing the moderating role of ELS on the relationship between the brain and Intelligence composite scores.

Mediation and moderated mediation analyses were performed in R environment v3.4.1 using the *lavaan* v0.6-7 package with 1000 bootstrapping replicates. Mediation models included the same covariates that were used in the regression model. *Parental separation or divorce* had to be excluded in a moderated mediation model because of model convergence issues.

## Results

### Participants

After quality control, our study cohort comprised 4,276 unrelated children of European ancestry with a mean age of 9-10 years from the ABCD study35. **The** demographic characteristics of the patients, stratified by records of ELS experience, are summarized in **Table 1** (**Supplementary Table 3**) **and Figure 2**. All covariates, except for BMI, showed significant differences in experience with one or more type of ELS (*abuse: P < 0*.*001, neglect: P < 0*.*05, and household challenges: P < 0*.*05*) (**Supplementary Table 4**).

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/T1)

Table 1. Demographic Characteristics of participants With European ancestry (N =4,276).
The total participants are divided into five intervals regarding the Composite score of ELS.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/30/2021.12.27.21268445/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/F2)

Figure 2. Flow chart of data selection and research design.
From the ABCD 2.0 Data release, genetic, brain imaging, and cognitive assessment data for 11,875 participants were collected. For the 8,496 individuals after initial quality control and GPS calculalion, we additionally removed samples of non-European ancestry or with any missing values, for a total of 4,276 participants included in the analysis.

### Correlation among GPS, Brain, and Intelligence

#### Cognitive Capacity GPS-Brain

Out of 992 brain morphometric features, 169 and 44 brain features correlated significantly with CP GPS (*P**FWE**< 0*.*05*) and EA GPS, respectively (*P**FWE**< 0*.*05*), when adjusted for age, sex, maternal education, income, BMI, study site, and marital status (**Figure 3**). Of note, global brain features - including total gray matter volume, total cortex volume, subcortical gray matter volume, and total cerebral white matter volume - were shown to have significant associations with both GPSs (**Supplementary Table 1**). White matter fiber counts and fractional anisotropy of the structural connectomes showed no significant correlations with either GPS (**Supplementary Figure 2**). To obtain cognitive capacity GPS-related brain representations, we linearly reduced the dimensionality of brain morphometric features correlated with each GPS by performing principal component analysis. We verified that only the first components substantially explained the variance of each GPS-related brain morphometric feature (Brain CP GPS: 44.5%; Brain EA GPS: 54.9%) (**Supplementary Figurea 3, 4, 5**) and were significantly correlated with cognitive capacity GPS (Brain CP GPS: *β = 0.10, P**FWE**< 0.001;* Brain EA GPS: *β = 0.15, P**FWE**< 0.001*) (**Supplementary Table 5**).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/30/2021.12.27.21268445/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/F3)

Figure 3. Manhattan plots of GPS-brain linear regression.
The results of generalized linear models between brain morphometric features and (a) CP GPS and (b) between brain morphometric features and EA GPS are presented usi ng the Bonferroni-cor rected p values. The dotted line indicates the significance threshold(*P* *FWE* =0.05). TheX-axis is the category of brain features: global features (Global), brain cortical area (Area), cortical volume (Volume), subcortical volume (Sub. Vol.), average cortical curvature(Avg. Curv), and cortical thickness (Thickness).

#### Brain-Intelligence Composite Score

Both *Brain* *CP GPS* and *Brain* *EA GPS* were significantly correlated with all three types of intelligence composite scores (*P**FWE**< 0.001*). Among the different types of intelligence composite scores, both brain representations showed the highest effect size on crystallized intelligence composite scores (*Brain* *CP GPS*: *β = 0.18; Brain* *EA GPS*: *β = 0.18*), with the lowest effect size on fluid intelligence composite scores (*Brain* *CP GPS*: *β = 0.05; Brain* *EA GPS*: β*= 0.05*). Total intelligence composite score was moderately correlated with brain representations (*Brain* *CP GPS*: *β = 0.12; Brain* *EA GPS*: *β = 0.12*) (**Supplementary Table 6**).

#### GPS-Intelligence Composite Score

Both CP and EA GPS were significantly correlated with the outcomes of all three of the different types of intelligence measures (*P**FWE**< 0.05)*. Test statistics showed that composite scores of crystallized intelligence appeared to have the strongest associations with both GPSs (*CP GPS:* β*= 0.17; EA GPS: β = 0.15*), followed by total intelligence composite scores (*CP GPS: β = 0.14; EA GPS: β = 0.10*) and the fluid intelligence composite score (*CP GPS: β = 0.09; EA GPS: β = 0.05*) (**Supplementary Table 6**).

### Mediation Analysis

Our analysis of the GPS-brain-cognition pathway (**Figure 1-a**) showed that the brain morphometric representations significantly mediated the causal effect of both GPSs on all types of intelligence (*P**FWE**< 0.001)*. Both GPSs showed the strongest direct effect on crystallized intelligence composite scores (β ***direct effect of CP GPS*** *= 0.15;* β ***direct effect of EA GPS*** *= 0.13)* compared to that on fluid intelligence composite scores (β ***direct effect of CP GPS*** *= 0.09;* β ***direct effect of EA GPS*** *= 0.04)* and total intelligence composite scores (β ***direct effect of CP GPS*** *= 0.13;* β ***direct effect of EA GPS*** *= 0.09)*. Both GPSs also showed the strongest indirect effect size on crystallized intelligence composite scores, which were determined by brain morphometric representation (β ***indirect effect of CP GPS*** *= 0.016;* β ***indirect effect of EA GPS*** *= 0.016)*, compared to fluid intelligence composite scores (β ***indirect effect of CP GPS*** *= 0.004;* β ***indirect effect of EA GPS*** *= 0.005)* and total intelligence composite scores (β ***indirect effect of CP GPS*** *= 0.010;* β ***indirect effect of EA GPS*** *= 0.010)* (**Table 2, Supplementary Table 7, Figure 4**).

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/T2)

Table 2. Results of mediation analysis

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/30/2021.12.27.21268445/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/F4)

Figure 4. Path diagrams for crystallized intelligence mediation models.
Standardizes regression coefficients for each path. Red dotted lines indicate the total indirect effect.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/12/30/2021.12.27.21268445/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/F5)

Figure 5. Effect of abuse on the EA GPS-brain-cognition pathway.
The left panel shows standardized regression coefficients for each path. Red dotted lines indicate the total indirect effect. Gray dotted lines indicate the moderation effect of ELS on the GPS-brain relationship. The right panel shows differences in the indirect impact caused by composite ELS scores (transparent gray areas indicate confidence intervals).

#### Path Analysis: Moderated Mediation effect of ELS

Prior to moderated mediation analysis, we first assessed the correlation between ELS and intelligence composite scores and the evaluated the moderation effects of ELS on the relationship between cognitive capacity GPSs and intelligence composite scores (**Figure 1-b**). Among the different types of ELS, only *neglect* showed significant correlations with all types of intelligence composite scores (**Supplementary Table 8**). In moderation analysis, no significant moderation effects of the ELS variables were detected (*PFWE > 0.05*) (**Supplementary Table 9**).

The following first-stage moderated mediation analysis (**Figure 1-c**) showed the significant negative moderation effect of *abuse* on the EA GPS-brain-intelligence composite score relationship. Notably, *abuse* itself did not affect EA GPS-related brain morphometric representation; however, the interaction between EA GPS and *abuse* showed a significant negative association with brain morphometric representation (*β = −0.04*, ***95% CI*** *= −0.07 ∼ − 0.014, PFWE < 0.05)*. Among the different kinds of intelligence composite scores, crystallized intelligence composite scores were the most severely affected by the interplay of ELS and cognitive capacity GPSs through the brain (***Index of Moderated Mediation Crystallized Intelligence Score*** *= −0.007*, ***95% CI*** *= −0.012 ∼ −0.002; PFWE < 0.05*) compared to the effect on fluid intelligence composite scores (***Index of Moderated Mediation Fluid Intelligence Score*** *= − 0.002*, ***95% CI*** *= −0.004∼0.0; PFWE < 0.01*) and total intelligence composite scores (***Index of Moderated Mediation Total Intelligence Score*** *= −0.005*, ***95% CI*** *= −0.009 ∼ −0.002; PFWE < 0.01*) (**Table 3, Supplementary Table 10**). In the second-stage moderated mediation models (**Figure 1-d**), we found no significant interactions between the brain and ELS on intelligence composite scores (**Supplementary Table 11**). Taken together, the negative moderation effects of *abuse* on the causal pathway between genetic and phenotypic intelligence were only observed when the brain served as a mediator. The moderated mediation effect was only significant for *abuse*, but not for other ELS variables, even in subscales (**Supplementary Tables 12, 13**).

View this table:
[Table 3.](http://medrxiv.org/content/early/2021/12/30/2021.12.27.21268445/T3)

Table 3. Moderated Mediation effect of *Abuse* on Educational Attainment GPS-Brain-Cognition

## Discussion

By leveraging the multimodal genetic, neuroimaging, and cognitive assessment data of 4,267 children of European ancestry, we examined the complex interplay among brain structure and genomic influence on cognitive development and how ELS affects the gene-brain-cognition triad. Our analysis reveals likely causal neurodevelopment pathways under the genetic control of intelligence and ELS modulation. We found that brain structural development partially mediates the genomic contributions to cognitive development in young children. Of note, ELS (*abuse*) specifically modulates the influence of cognitive capacity GPSs on the brain, this modulatory effect then affects cognition. These results show the potential causal relationships between ELS (*abuse*) and the genetic pathway for neurocognitive development, accounting for biological how it may lead to cognitive deficits in children.

Extending the existing literature reporting the influence of genetics on the brain and cognition3,72–77, our study shows a more specific multifactorial model in which brain structural development mediates the genetic impact on cognitive development. Of note, this study shows that the association of the cognitive capacity GPS was only significant with gray matter phenotypes but not with white matter phenotypes. This result may be related to the developmental characteristics observed during preadolescence (9 to 10 years old), i.e., marked changes in gray matter, followed by relatively greater white matter changes or myelination78. Additionally, the brain mediation effect may be widespread involving several brain regions, considering that the brain mediator was the first principal component derived from hundreds of brain morphological features. Future research should examine how genetic influence is linked to longitudinal changes in brain structure over one’s lifetime.

The key contribution of this study to the literature is the discovery that ELS modulates the magnitude of the impact of genomic influence on neurocognitive development. This finding may account for a potential mechanism of the long-term effect of ELS on cognitive deficits. Note that our brain imaging results shows that ELS modulated the genomic influence on the cortical morphology of the various brain regions responsible for cognition and intelligence. In addition to the literature documenting the impact of ELS on affective processing, stress physiology and the relevant brain system3,72–77, this finding shows that ELS impacts cognitive brain circuitry.

Our moderated mediation results support genetic differential susceptibility to ELS. The genetic differential susceptibility model proposes that genetic variations of individuals regulate sensitivity or susceptibility to environmental influences79–81. In our analysis, as a result of the ELS moderation effects on the gene-to-brain pathway, children with the same level of ELS showed different cognitive outcomes depending on the magnitude of the EA GPS. That is, a child with a higher EA GPS showed a greater negative impact of ELS on the cognitive outcome, whereas one with a lower EA GPS showed a smaller negative impact.

Perhaps, in the case of individual with a high GPS, ELS may negatively regulate the expression of the genes responsible for the brain and cognitive development (sequentially as supported in our model). Consequently, otherwise, positive effects of high GPSs on cognitive outcomes may be eliminated under the effect of ELS, resulting in lower cognitive outcomes compared with the individuals with high GPSs and no ELS. A crucial outstanding question is whether this modulatory effect of ELS on the cognitive capacity GPS is linked to epigenetics. Animal research shows that ELS induces epigenetic modification in DNA methylation sites that are associated with a wide array of neural events, such as neurotransmitter biosynthesis, neurological system processes, glial cell proliferation, neurogenesis in the hippocampus, neural migration in the cerebral cortex, neuroplasticity, and neurodevelopmental delay82–84; many of these events are critically related to cognitive development. This line of research may allow precise risk stratification and the development of more individualized intervention strategies.

Another novel aspect of this study was the application of the GPS approach, rather than focusing on selected common variants of DNA, which permitted integrative assessment of the polygenic contribution to cognitive development and revealed its impact on the wide array of brain circuits. A recent study reported a similar finding of the association between cognitive capacity GPS and brain volumes in adult samples in the UK Biobank73. Our study showed that the gene-brain-cognition relationship manifests as early as preadolescence.

Our observations suggested that *physical or sexual abuse* significantly modulates the “gene-to-brain” path, which then leads to “cognition” but not the “gene-to-cognition” path. This finding is in line with previous ELS studies reporting that *physical neglect* and *abuse* negatively impact genetic regulation of brain systems, such as hippocampal glucocorticoid receptors, neurogenesis, and regulation of brain development6,7,9,82.

The direct effect of cognitive capacity GPSs on cognition was tenfold larger in magnitude than the indirect effects on the brain. Note that the gray matter representations derived from ROI-level morphometry may only partially account for the gene-to-brain effects. Indeed, a morphometric analysis may not fully detect subtle changes in the brain tissues under the control of epigenetic (or any other environmental) mechanisms, such as vascularization, neurogenesis, and synaptogenesis. Future research may elucidate sensitive brain representations related to genetic influence and environmental modulation.

Another possible explanation of the small indirect effect on the brain may be related to the limitations of linear models in testing complex nonlinear relationships. We used PCA to extract the brain representation (i.e., principal components) from thousands of brain morphometric variables. We used this method for the ease of statistical modeling (mediation analysis); however, since this method uses a linear and orthogonal transformation of data, we admit that this representation might be too simplistic to capture nonlinear relationships among the brain variables.

Our study showed the utility of our multitrait (cognitive performance and educational attainment) genomic approach for revealing the different patterns of genetic influences on the brain and cognition. Years of education (educational attainment, EA) is a widely used proxy for an intelligence phenotype. This is because of its high genetic correlation with intelligence and the ease of assessment compared to an evaluation of cognitive ability that would require behavioral tests (hence, it is unfavorable in large GWASs). For example, the literature shows that EA GPS correlates with cognitive functioning85 and the development of behaviors86.

Although both cognitive capacity GPSs showed significant effects on the gene-brain-cognition pathway, only EA GPS showed a significant modulatory effect of ELS on the pathway; CP GPS did not have a notable effect on the pathway. This observed ELS effect on EA GPS may be related to genetic loading not only for cognitive capacity but also for noncognitive skills and traits required for successful education attainment87. On the other hand, CP GPS correlated with a greater number of brain morphometric features than EA GPS did. This may reflect that CP GPS represents the genetic influence primarily on biological processes (i.e., brain development) directly linked to cognitive performance, but not so much on noncognitive processes. These points may be used to guide to future studies examining the gene-environment interaction in cognitive development.

Our results indicated greater magnitudes of the genetic influence on crystallized intelligence than on fluid intelligence during development. Crystallized intelligence is the ability to apply prior knowledge to problem solving; fluid intelligence is the ability to reason in novel situations without prior knowledge65,88,89. Since the genetic influence (both direct and indirect) was up to three times larger on crystallized intelligence compared to that on fluid intelligence, the impact of ELS via the gene-brain pathway was proportionally greater on crystallized intelligence than on fluid intelligence. This result adds granularity regarding the specificity of ELS impacts.

A few limitations of the study should be noted. First, it was not within the scope of this study to examine brain function, which is influenced by ELS through the triangular relationship. It is well known that ELS leads to changes in brain function3,8. Previous animal studies support that the function of adaptive gene-brain-cognition feedback systems may be affected by ELS8. Thus, the comprehensive impact of ELS on this triangular relationship needs to be further investigated in terms of brain function to reveal all the causal mechanisms. Second, we created composite scores of ELS with several proxy measures instead of using the established measure for childhood trauma, i.e., the Adverse Childhood Experiences (ACE) questionnaire90. This was mainly due to the unavailability of the ACE questionnaire data in the ABCD study. However, it is of note that we included questionnaires similar to those of the ACE questionnaire and included several ELS-related measurements68, thus aggregating various adverse environmental factors.

## Supporting information

Supplementary Tables [[supplements/268445_file10.xlsx]](pending:yes)

Supplementary Figures [[supplements/268445_file11.docx]](pending:yes)

## Data Availability

All data used are available online at NIMH Data Archive

[https://nda.nih.gov](https://nda.nih.gov) 

## Author contributions

*   Study concept and design: Y.J., J.C.

*   Acquisition, analysis, or interpretation of data: Y.J., H.W., S.M., K.K., H.K.

*   Drafting of the manuscript: Y.J., H.W., S.M., J.C.

*   Critical revision of the manuscript for important intellectual content: Y.J., J.K., H.W. S.M., K.K., W.A.

*   Statistical analysis: H.W., S.M.

*   Obtained funding: J.C.

*   Study supervision: Y.J. J.C.

## Data and Code Availability

Codes and data are freely available for reproducibility ([https://github.com/Transconnectome/ELS\_on_genes-brain-cognition](https://github.com/Transconnectome/ELS_on_genes-brain-cognition))

## Competing Interest Disclosures

None of the authors have significant competing financial, professional, or personal interests that might have influenced the performance or presentation of the work described in the manuscript.

## Acknowledgments

This work was supported by the New Faculty Startup Fund from Seoul National University and a Research grant from Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, South Korea (2021R1I1A1A01054995).

*   Received December 27, 2021.
*   Revision received December 27, 2021.
*   Accepted December 30, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Saleh, A. et al. Effects of early life stress on depression, cognitive performance and brain morphology. Psychol. Med. 47, 171–181 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291716002403&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27682320&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

2.  2.Pesonen, A.-K. et al. Cognitive ability and decline after early life stress exposure. Neurobiol. Aging 34, 1674–1679 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neurobiolaging.2012.12.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23337341&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317417100017&link_type=ISI) 

3.  3.Pechtel, P. & Pizzagalli, D. A. Effects of early life stress on cognitive and affective function: an integrated review of human literature. Psychopharmacology 214, 55–70 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00213-010-2009-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

4.  4.De Bellis, M. D. Developmental traumatology: the psychobiological development of maltreated children and its implications for research, treatment, and policy. Dev. Psychopathol. 13, 539–564 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0954579401003078&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11523847&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170540800007&link_type=ISI) 

5.  5.Peña, C. J. et al. Early life stress confers lifelong stress susceptibility in mice via ventral tegmental area OTX2. Science 356, 1185–1188 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNTYvNjM0My8xMTg1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMzAvMjAyMS4xMi4yNy4yMTI2ODQ0NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

6.  6.Niwa, M. et al. Adolescent stress-induced epigenetic control of dopaminergic neurons via glucocorticoids. Science 339, 335–339 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMzkvNjExNy8zMzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8zMC8yMDIxLjEyLjI3LjIxMjY4NDQ1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

7.  7.Fuge, P. et al. Interaction of Early Life Stress and Corticotropin-Releasing Hormone Receptor Gene: Effects on Working Memory. Biological Psychiatry vol. 76 888–894 (2014).
    
    

8.  8.Chen, Y. & Baram, T. Z. Toward Understanding How Early-Life Stress Reprograms Cognitive and Emotional Brain Networks. Neuropsychopharmacology 41, 197–206 (2016).
    
    

9.  9.McGowan, P. O. et al. Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse. Nat. Neurosci. 12, 342–348 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.2270&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19234457&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263577900020&link_type=ISI) 

10. 10.Koenen, K. C., Moffitt, T. E., Caspi, A., Taylor, A. & Purcell, S. Domestic violence is associated with environmental suppression of IQ in young children. Dev. Psychopathol. 15, 297–311 (2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0954579403000166&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12931829&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183838200004&link_type=ISI) 

11. 11.Adjimann, T. S., Argañaraz, C. V. & Soiza-Reilly, M. Serotonin-related rodent models of early-life exposure relevant for neurodevelopmental vulnerability to psychiatric disorders. Transl. Psychiatry 11, 280 (2021).
    
    

12. 12.Dong, E., Guidotti, A., Zhang, H. & Pandey, S. C. Prenatal stress leads to chromatin and synaptic remodeling and excessive alcohol intake comorbid with anxiety-like behaviors in adult offspring. Neuropharmacology 140, 76–85 (2018).
    
    

13. 13.Ivy, A. S. et al. Hippocampal dysfunction and cognitive impairments provoked by chronic early-life stress involve excessive activation of CRH receptors. J. Neurosci. 30, 13005–13015 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIzMC8zOS8xMzAwNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzMwLzIwMjEuMTIuMjcuMjEyNjg0NDUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

14. 14.Spinelli, S. et al. Early-life stress induces long-term morphologic changes in primate brain. Arch. Gen. Psychiatry 66, 658–665 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archgenpsychiatry.2009.52&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19487631&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266566800011&link_type=ISI) 

15. 15.Sun, X., Zhang, Y., Li, X., Liu, X. & Qin, C. Early-Life Neglect Alters Emotional and Cognitive Behavior in a Sex-Dependent Manner and Reduces Glutamatergic Neuronal Excitability in the Prefrontal Cortex. Frontiers in Psychiatry vol. 11 (2021).
    
    

16. 16.Bock, J., Riedel, A. & Braun, K. Differential changes of metabolic brain activity and interregional functional coupling in prefronto-limbic pathways during different stress conditions: functional imaging in freely behaving rodent pups. Front. Cell. Neurosci. 6, 19 (2012).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22590453&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

17. 17.Molet, J. et al. Fragmentation and high entropy of neonatal experience predict adolescent emotional outcome. Transl. Psychiatry 6, e702 (2016).
    
    

18. 18.Brunson, K. L. et al. Mechanisms of late-onset cognitive decline after early-life stress. J. Neurosci. 25, 9328–9338 (2005).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIyNS80MS85MzI4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMzAvMjAyMS4xMi4yNy4yMTI2ODQ0NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

19. 19.Radley, J. J. et al. Repeated stress alters dendritic spine morphology in the rat medial prefrontal cortex. J. Comp. Neurol. 507, 1141–1150 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cne.21588&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18157834&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252523000010&link_type=ISI) 

20. 20.Aisa, B., Tordera, R., Lasheras, B., Del Río, J. & Ramírez, M. J. Cognitive impairment associated to HPA axis hyperactivity after maternal separation in rats. Psychoneuroendocrinology vol. 32 256–266 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2006.12.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17307298&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245741300005&link_type=ISI) 

21. 21.Alviña, K., Jodeiri Farshbaf, M. & Mondal, A. K. Long term effects of stress on hippocampal function: Emphasis on early life stress paradigms and potential involvement of neuropeptide Y. J. Neurosci. Res. 99, 57–66 (2021).
    
    

22. 22.Kronman, H. et al. Long-term behavioral and cell-type-specific molecular effects of early life stress are mediated by H3K79me2 dynamics in medium spiny neurons. Nat. Neurosci. 24, 667–676 (2021).
    
    

23. 23.Arabadzisz, D. et al. Primate early life stress leads to long-term mild hippocampal decreases in corticosteroid receptor expression. Biol. Psychiatry 67, 1106–1109 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2009.12.016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20132928&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277918000014&link_type=ISI) 

24. 24.Wang, X.-D. et al. Forebrain CRF1 Modulates Early-Life Stress-Programmed Cognitive Deficits. Journal of Neuroscience vol. 31 13625–13634 (2011).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIzMS8zOC8xMzYyNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzMwLzIwMjEuMTIuMjcuMjEyNjg0NDUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

25. 25.De Bellis, M. D. et al. Brain structures in pediatric maltreatment-related posttraumatic stress disorder: a sociodemographically matched study. Biol. Psychiatry 52, 1066–1078 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-3223(02)01459-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12460690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180352700004&link_type=ISI) 

26. 26.Jackowski, A. P. et al. Corpus callosum in maltreated children with posttraumatic stress disorder: a diffusion tensor imaging study. Psychiatry Res. 162, 256–261 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pscychresns.2007.08.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18296031&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255122100009&link_type=ISI) 

27. 27.Gatt, J. M. et al. Interactions between BDNF Val66Met polymorphism and early life stress predict brain and arousal pathways to syndromal depression and anxiety. Mol. Psychiatry 14, 681–695 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2008.143&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19153574&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267284800005&link_type=ISI) 

28. 28.Caspi, A. et al. Influence of life stress on depression: moderation by a polymorphism in the 5-HTT gene. Science 301, 386–389 (2003).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMDEvNTYzMS8zODYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8zMC8yMDIxLjEyLjI3LjIxMjY4NDQ1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

29. 29.McLaughlin, K. A. et al. Maltreatment Exposure, Brain Structure, and Fear Conditioning in Children and Adolescents. Neuropsychopharmacology 41, 1956–1964 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2015.365&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26677946&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

30. 30.Lee, J. J. et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat. Genet. 50, 1112–1121 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0147-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30038396&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

31. 31.van der Sluis, S., Willemsen, G., de Geus, E. J. C., Boomsma, D. I. & Posthuma, D. Gene-environment interaction in adults’ IQ scores: measures of past and present environment. Behav. Genet. 38, 348–360 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10519-008-9212-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18535898&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257122300002&link_type=ISI) 

32. 32.Harden, K. P., Turkheimer, E. & Loehlin, J. C. Genotype by environment interaction in adolescents’ cognitive aptitude. Behav. Genet. 37, 273–283 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10519-006-9113-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16977503&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245407100002&link_type=ISI) 

33. 33.Plomin, R. & von Stumm, S. The new genetics of intelligence. Nat. Rev. Genet. 19, 148–159 (2018).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

34. 34.Jernigan, T. & Jernigan, T. Adolescent brain cognitive development study (ABCD) - annual release 2.0. (2019) doi:10.15154/1503209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15154/1503209&link_type=DOI) 

35. 35.Karcher, N. R. & Barch, D. M. The ABCD study: understanding the development of risk for mental and physical health outcomes. Neuropsychopharmacology 46, 131–142 (2021).
    
    

36. 36.Das, S. et al. Next-generation genotype imputation service and methods. Nat. Genet. 48, 1284–1287 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3656&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27571263&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

37. 37.Loh, P.-R. et al. Reference-based phasing using the Haplotype Reference Consortium panel. Nature Genetics vol. 48 1443–1448 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3679&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27694958&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

38. 38.Raj, A., Stephens, M. & Pritchard, J. K. fastSTRUCTURE: variational inference of population structure in large SNP data sets. Genetics 197, 573–589 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6OToiMTk3LzIvNTczIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMzAvMjAyMS4xMi4yNy4yMTI2ODQ0NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

39. 39.Conomos, M. P., Miller, M. B. & Thornton, T. A. Robust inference of population structure for ancestry prediction and correction of stratification in the presence of relatedness. Genet. Epidemiol. 39, 276–293 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21896&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25810074&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

40. 40.Conomos, M. P., Reiner, A. P., Weir, B. S. & Thornton, T. A. Model-free Estimation of Recent Genetic Relatedness. Am. J. Hum. Genet. 98, 127–148 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2015.11.022&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26748516&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

41. 41.Choi, S. W. & O’Reilly, P. F. PRSice-2: Polygenic Risk Score software for biobank-scale data. Gigascience 8, (2019).
    
    

42. 42.Garavan, H. et al. Recruiting the ABCD sample: Design considerations and procedures. Developmental Cognitive Neuroscience vol. 32 16–22 (2018).
    
    

43. 43.Casey, B. J. et al. The Adolescent Brain Cognitive Development (ABCD) study: Imaging acquisition across 21 sites. Developmental Cognitive Neuroscience vol. 32 43–54 (2018).
    
    

44. 44.Jovicich, J. et al. Reliability in multi-site structural MRI studies: effects of gradient non-linearity correction on phantom and human data. Neuroimage 30, 436–443 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2005.09.046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16300968&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236682200009&link_type=ISI) 

45. 45.Dale, A. M., Fischl, B. & Sereno, M. I. Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 9, 179–194 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/nimg.1998.0395&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9931268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078608900001&link_type=ISI) 

46. 46.Ségonne, F. et al. A hybrid approach to the skull stripping problem in MRI. Neuroimage 22, 1060–1075 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2004.03.032&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15219578&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222423200004&link_type=ISI) 

47. 47.Fischl, B., Liu, A. & Dale, A. M. Automated manifold surgery: constructing geometrically accurate and topologically correct models of the human cerebral cortex. IEEE Trans. Med. Imaging 20, 70–80 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/42.906426&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11293693&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167324900007&link_type=ISI) 

48. 48.Ségonne, F., Pacheco, J. & Fischl, B. Geometrically accurate topology-correction of cortical surfaces using nonseparating loops. IEEE Trans. Med. Imaging 26, 518–529 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TMI.2006.887364&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17427739&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245599200010&link_type=ISI) 

49. 49.Fischl, B., Sereno, M. I., Tootell, R. B. H. & Dale, A. M. High-resolution intersubject averaging and a coordinate system for the cortical surface. Human Brain Mapping vol. 8 272–284 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/(SICI)1097-0193(1999)8:4<272::AID-HBM10>3.0.CO;2-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10619420&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000084127200010&link_type=ISI) 

50. 50.Hagler, D. J., Jr. et al. Image processing and analysis methods for the Adolescent Brain Cognitive Development Study. Neuroimage 202, 116091 (2019).
    
    

51. 51.Zhuang, J. et al. Correction of eddy-current distortions in diffusion tensor images using the known directions and strengths of diffusion gradients. J. Magn. Reson. Imaging 24, 1188–1193 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmri.20727&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17024663&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

52. 52.Hagler, D. J. et al. Automated white-matter tractography using a probabilistic diffusion tensor atlas: Application to temporal lobe epilepsy. Human Brain Mapping vol. 30 1535–1547 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.20619&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18671230&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265727600012&link_type=ISI) 

53. 53.Leemans, A. & Jones, D. K. TheB-matrix must be rotated when correcting for subject motion in DTI data. Magnetic Resonance in Medicine vol. 61 1336–1349 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.21890&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19319973&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266429900009&link_type=ISI) 

54. 54.Chang, L.-C., Jones, D. K. & Pierpaoli, C. RESTORE: robust estimation of tensors by outlier rejection. Magn. Reson. Med. 53, 1088–1095 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.20426&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15844157&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228796900014&link_type=ISI) 

55. 55.Holland, D., Kuperman, J. M. & Dale, A. M. Efficient correction of inhomogeneous static magnetic field-induced distortion in Echo Planar Imaging. Neuroimage 50, 175–183 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2009.11.044&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19944768&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274810100016&link_type=ISI) 

56. 56.Tournier, J.-D. et al. MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. Neuroimage 202, 116137 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/J.NEUROIMAGE.2019.116137&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

57. 57.Cha, J. et al. Abnormal reward circuitry in anorexia nervosa: A longitudinal, multimodal MRI study. Hum. Brain Mapp. 37, 3835–3846 (2016).
    
    

58. 58.Cha, J. et al. Neural Correlates of Aggression in Medication-Naive Children with ADHD: Multivariate Analysis of Morphometry and Tractography. Neuropsychopharmacology 40, 1717–1725 (2015).
    
    

59. 59.Veraart, J. et al. Denoising of diffusion MRI using random matrix theory. NeuroImage vol. 142 394–406 (2016).
    
    

60. 60.Tustison, N. J. et al. N4ITK: improved N3 bias correction. IEEE Trans. Med. Imaging 29, 1310–1320 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/TMI.2010.2046908&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20378467&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278535800009&link_type=ISI) 

61. 61.Ciccarelli, O. et al. Probabilistic diffusion tractography: a potential tool to assess the rate of disease progression in amyotrophic lateral sclerosis. Brain 129, 1859–1871 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awl100&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16672290&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238761200026&link_type=ISI) 

62. 62.Calamante, F., Tournier, J.-D., Jackson, G. D. & Connelly, A. Track-density imaging (TDI): super-resolution white matter imaging using whole-brain track-density mapping. Neuroimage 53, 1233–1243 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2010.07.024&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20643215&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

63. 63.Bedre, R. reneshbedre/bioinfokit: Bioinformatics data analysis and visualization toolkit. (2021) doi:10.5281/zenodo.4422035.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.5281/zenodo.4422035&link_type=DOI) 

64. 64.Gershon, R. C. et al. NIH toolbox for assessment of neurological and behavioral function. Neurology 80, S2–6 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3182872e5f&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23479538&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

65. 65.Cattell, R. B. Theory of fluid and crystallized intelligence: A critical experiment. Journal of Educational Psychology vol. 54 1–22 (1963).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/h0046743&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1963CCC1100001&link_type=ISI) 

66. 66.Weintraub, S. et al. Cognition assessment using the NIH Toolbox. Neurology 80, S54–64 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1212/WNL.0b013e3182872ded&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23479546&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

67. 67.Akshoomoff, N. et al. VIII. NIH Toolbox Cognition Battery (CB): composite scores of crystallized, fluid, and overall cognition. Monogr. Soc. Res. Child Dev. 78, 119–132 (2013).
    
    

68. 68.Hoffman, E. A. et al. Stress exposures, neurodevelopment and health measures in the ABCD study. Neurobiology of Stress vol. 10 100157 (2019).
    
    

69. 69.Edwards, J. R. & Lambert, L. S. Methods for integrating moderation and mediation: a general analytical framework using moderated path analysis. Psychol. Methods 12, 1–22 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/1082-989X.12.1.1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17402809&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244988200001&link_type=ISI) 

70. 70.Hayes, A. F. An Index and Test of Linear Moderated Mediation. Multivariate Behav. Res. 50, 1–22 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00273171.2014.962683&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26609740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

71. 71.Preacher, K. J., Rucker, D. D. & Hayes, A. F. Addressing Moderated Mediation Hypotheses: Theory, Methods, and Prescriptions. Multivariate Behav. Res. 42, 185–227 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00273170701341316&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26821081&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247879400007&link_type=ISI) 

72. 72.Thompson, P. M. et al. Genetic influences on brain structure. Nat. Neurosci. 4, 1253–1258 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn758&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11694885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000172525000023&link_type=ISI) 

73. 73.Elliott, M. L. et al. A Polygenic Score for Higher Educational Attainment is Associated with Larger Brains. Cereb. Cortex 29, 3496–3504 (2019).
    
    

74. 74.Mitchell, B. L. et al. Educational attainment polygenic scores are associated with cortical total surface area and regions important for language and memory. Neuroimage 212, 116691 (2020).
    
    

75. 75.Posthuma, D. et al. The association between brain volume and intelligence is of genetic origin. Nature neuroscience vol. 5 83–84 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn0202-83&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11818967&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173558400002&link_type=ISI) 

76. 76.Pol, H. E. H. et al. Genetic Contributions to Human Brain Morphology and Intelligence. Journal of Neuroscience vol. 26 10235–10242 (2006).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIyNi80MC8xMDIzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEyLzMwLzIwMjEuMTIuMjcuMjEyNjg0NDUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

77. 77.van Leeuwen, M. et al. A genetic analysis of brain volumes and IQ in children. Intelligence vol. 37 181–191 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.intell.2008.10.005&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264911900008&link_type=ISI) 

78. 78.Gogtay, N. et al. Dynamic mapping of human cortical development during childhood through early adulthood. Proc. Natl. Acad. Sci. U. S. A. 101, 8174–8179 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAxLzIxLzgxNzQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8xMi8zMC8yMDIxLjEyLjI3LjIxMjY4NDQ1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

79. 79.Belsky, J. Variation in Susceptibility to Environmental Influence: An Evolutionary Argument. Psychological Inquiry vol. 8 182–186 (1997).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1207/s15327965pli0803_3&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XH52900003&link_type=ISI) 

80. 80.Boyce, W. T. & Ellis, B. J. Biological sensitivity to context: I. An evolutionary-developmental theory of the origins and functions of stress reactivity. Dev. Psychopathol. 17, 271–301 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0954579405050145&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16761546&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000229244700001&link_type=ISI) 

81. 81.Silveira, P. P. et al. Genetic Differential Susceptibility to Socioeconomic Status and Childhood Obesogenic Behavior: Why Targeted Prevention May Be the Best Societal Investment. JAMA Pediatr. 170, 359–364 (2016).
    
    

82. 82.Cecil, C. A. M. et al. Epigenetic signatures of childhood abuse and neglect: Implications for psychiatric vulnerability. J. Psychiatr. Res. 83, 184–194 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpsychires.2016.09.010&link_type=DOI) 

83. 83.Yang, B.-Z. et al. Child abuse and epigenetic mechanisms of disease risk. Am. J. Prev. Med. 44, 101–107 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amepre.2012.10.012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23332324&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

84. 84.Labonté, B. et al. Genome-wide epigenetic regulation by early-life trauma. Arch. Gen. Psychiatry 69, 722–731 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archgenpsychiatry.2011.2287&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22752237&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305907600012&link_type=ISI) 

85. 85.Rietveld, C. A. et al. GWAS of 126,559 individuals identifies genetic variants associated with educational attainment. Science 340, 1467–1471 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNDAvNjEzOS8xNDY3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMTIvMzAvMjAyMS4xMi4yNy4yMTI2ODQ0NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

86. 86.Belsky, D. W. et al. The Genetics of Success: How Single-Nucleotide Polymorphisms Associated With Educational Attainment Relate to Life-Course Development. Psychol. Sci. 27, 957–972 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0956797616643070&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27251486&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

87. 87.Demange, P. A. et al. Investigating the genetic architecture of noncognitive skills using GWAS-by-subtraction. Nat. Genet. 53, 35–44 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-00754-2&link_type=DOI) 

88. 88.Merrifield, P. R. & Cattell, R. B. Abilities: Their Structure, Growth, and Action. American Educational Research Journal vol. 12 516 (1975).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3102/00028312012004516&link_type=DOI) 

89. 89.Christoforou, A. et al. GWAS-based pathway analysis differentiates between fluid and crystallized intelligence. Genes Brain Behav. 13, 663–674 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/gbb.12152&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24975275&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 

90. 90.Felitti, V. J. et al. Relationship of childhood abuse and household dysfunction to many of the leading causes of death in adults. The Adverse Childhood Experiences (ACE) Study. Am. J. Prev. Med. 14, 245–258 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0749-3797(98)00017-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9635069&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F12%2F30%2F2021.12.27.21268445.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000074088700001&link_type=ISI)