Hostname: page-component-cd9895bd7-dzt6s Total loading time: 0 Render date: 2024-12-21T11:50:57.224Z Has data issue: false hasContentIssue false

Genetic influences on externalizing psychopathology overlap with cognitive functioning and show developmental variation

Published online by Cambridge University Press:  31 March 2021

Josephine Mollon*
Affiliation:
Department of Psychiatry, Boston Children’s Hospital, Harvard Medical School, Boston, Massachusetts, USA
Emma E. M. Knowles
Affiliation:
Department of Psychiatry, Boston Children’s Hospital, Harvard Medical School, Boston, Massachusetts, USA
Samuel R. Mathias
Affiliation:
Department of Psychiatry, Boston Children’s Hospital, Harvard Medical School, Boston, Massachusetts, USA
Amanda Rodrigue
Affiliation:
Department of Psychiatry, Boston Children’s Hospital, Harvard Medical School, Boston, Massachusetts, USA
Tyler M. Moore
Affiliation:
Brain Behavior Laboratory, Department of Psychiatry, Perelman School of Medicine, Penn-CHOP Lifespan Brain Institute, University of Pennsylvania, Philadelphia, Pennsylvania, USA
Monica E. Calkins
Affiliation:
Brain Behavior Laboratory, Department of Psychiatry, Perelman School of Medicine, Penn-CHOP Lifespan Brain Institute, University of Pennsylvania, Philadelphia, Pennsylvania, USA
Ruben C. Gur
Affiliation:
Brain Behavior Laboratory, Department of Psychiatry, Perelman School of Medicine, Penn-CHOP Lifespan Brain Institute, University of Pennsylvania, Philadelphia, Pennsylvania, USA
Juan Manuel Peralta
Affiliation:
South Texas Diabetes and Obesity Institute, School of Medicine, University of Texas of the Rio Grande Valley, Brownsville, Texas, USA
Daniel J. Weiner
Affiliation:
Analytic and Translational Genetics Unit, Department of Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, USA Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA
Elise B. Robinson
Affiliation:
Analytic and Translational Genetics Unit, Department of Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, USA Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA
Raquel E. Gur
Affiliation:
Brain Behavior Laboratory, Department of Psychiatry, Perelman School of Medicine, Penn-CHOP Lifespan Brain Institute, University of Pennsylvania, Philadelphia, Pennsylvania, USA
John Blangero
Affiliation:
South Texas Diabetes and Obesity Institute, School of Medicine, University of Texas of the Rio Grande Valley, Brownsville, Texas, USA
Laura Almasy
Affiliation:
Department of Genetics, Perelman School of Medicine, Penn-CHOP Lifespan Brain Institute, University of Pennsylvania, Philadelphia, Pennsylvania, USA
David C. Glahn
Affiliation:
Department of Psychiatry, Boston Children’s Hospital, Harvard Medical School, Boston, Massachusetts, USA Olin Neuropsychiatry Research Center, Institute of Living, Hartford, Connecticut, USA
*
*Author for correspondence: Josephine Mollon, E-mail: josephine.mollon@childrens.harvard.edu

Abstract

Background

Questions remain regarding whether genetic influences on early life psychopathology overlap with cognition and show developmental variation.

Methods

Using data from 9,421 individuals aged 8–21 from the Philadelphia Neurodevelopmental Cohort, factors of psychopathology were generated using a bifactor model of item-level data from a psychiatric interview. Five orthogonal factors were generated: anxious-misery (mood and anxiety), externalizing (attention deficit hyperactivity and conduct disorder), fear (phobias), psychosis-spectrum, and a general factor. Genetic analyses were conducted on a subsample of 4,662 individuals of European American ancestry. A genetic relatedness matrix was used to estimate heritability of these factors, and genetic correlations with executive function, episodic memory, complex reasoning, social cognition, motor speed, and general cognitive ability. Gene × Age analyses determined whether genetic influences on these factors show developmental variation.

Results

Externalizing was heritable (h2 = 0.46, p = 1 × 10−6), but not anxious-misery (h2 = 0.09, p = 0.183), fear (h2 = 0.04, p = 0.337), psychosis-spectrum (h2 = 0.00, p = 0.494), or general psychopathology (h2 = 0.21, p = 0.040). Externalizing showed genetic overlap with face memory (ρg = −0.412, p = 0.004), verbal reasoning (ρg = −0.485, p = 0.001), spatial reasoning (ρg = −0.426, p = 0.010), motor speed (ρg = 0.659, p = 1x10−4), verbal knowledge (ρg = −0.314, p = 0.002), and general cognitive ability (g)(ρg = −0.394, p = 0.002). Gene × Age analyses revealed decreasing genetic variance (γg = −0.146, p = 0.004) and increasing environmental variance (γe = 0.059, p = 0.009) on externalizing.

Conclusions

Cognitive impairment may be a useful endophenotype of externalizing psychopathology and, therefore, help elucidate its pathophysiological underpinnings. Decreasing genetic variance suggests that gene discovery efforts may be more fruitful in children than adolescents or young adults.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2021. Published by Cambridge University Press on behalf of the European Psychiatric Association

Introduction

Psychiatric symptoms in early life are associated with poor cognition [Reference Calkins, Merikangas, Moore, Burstein, Behr and Satterthwaite1]. For example, psychotic symptoms in childhood and adolescence are associated with cognitive impairment [Reference Bora, Lin, Wood, Yung, McGorry and Pantelis2,Reference Brewer, Wood, Phillips, Francey, Pantelis and Yung3]. There is also evidence for IQ deficits in children with conduct problems [Reference Deater-Deckard, Mullineaux, Beekman, Petrill, Schatschneider and Thompson4], vocabulary deficits in children with aggression [Reference Dionne, Tremblay, Boivin, Laplante and Pérusse5], visuospatial deficits in children with hyperactivity [Reference Seguin, Parent, Tremblay and Zelazo6], and social deficits in adolescents with externalizing problems [Reference Gur, Moore, Calkins, Ruparel and Gur7]. Small, generalized deficits have also been reported in children with symptoms of anxiety and depression [Reference Lundy, Silva, Kaemingk, Goodwin and Quan8].

Early life psychopathology is also underpinned by genes with evidence of substantial genetic effects on childhood and adolescent psychopathology [Reference Wesseldijk, Fedko, Bartels, Nivard, van Beijsterveldt and Boomsma9,Reference Polderman, Benyamin, de Leeuw, Sullivan, van Bochoven and Visscher10], as well as specific psychotic [Reference Polanczyk, Moffitt, Arseneault, Cannon, Ambler and Keefe11,Reference Zavos, Freeman, Haworth, McGuire, Plomin and Cardno12], externalizing [Reference Hannigan, Walaker, Waszczuk, McAdams and Eley13Reference Huizink, Van Den Berg, van der Ende and Verhulst15], internalizing [Reference Ormel, Oldehinkel, Ferdinand, Hartman, De Winter and Veenstra16], and anxiety/depression symptoms [Reference Gregory and Eley17Reference Nivard, Dolan, Kendler, Kan, Willemsen and van Beijsterveldt21]. Since cognition in early life is also influenced by genes [Reference Haworth, Wright, Luciano, Martin, de Geus and van Beijsterveldt22,Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23], recent studies have examined whether there is a genetic component to cognitive correlates of early life psychopathology. Evidence from twin and family studies shows genetic overlap between psychopathology and cognitive ability [Reference Jacobs, Rijsdijk, Derom, Danckaerts, Thiery and Derom24], psychopathology and executive functions [Reference Harden, Engelhardt, Mann, Patterson, Grotzinger and Savicki25], attention deficit hyperactivity disorder (ADHD) and IQ [Reference Kuntsi, Eley, Taylor, Hughes, Asherson and Caspi26], ADHD and executive functions [Reference Bidwell, Willcutt, Defries and Pennington27,Reference Doyle, Biederman, Seidman, Reske-Nielsen and Faraone28], and inattention and attention regulation [Reference Wang, Deater-Deckard, Petrill and Thompson29]. GWAS evidence shows genetic overlap between ADHD and intelligence [Reference Demontis, Walters, Martin, Mattheisen, Als and Agerbo30-Reference Savage, Jansen, Stringer, Watanabe, Bryois and de Leeuw32], academic underperformance [Reference O’Connell, Shadrin, Smeland, Bahrami, Frei and Bettella31,Reference Liu, Li, Viding, Asherson and Pingault33,Reference Johnson, McGue and Iacono34], and executive functions [Reference Donati, Dumontheil and Meaburn35]. However, a comprehensive examination of genetic underpinnings of cognitive impairment in early life psychopathology, that is, comprising multiple dimensions of psychopathology and cognition, is lacking.

Gaps in knowledge also exist regarding whether genetic influences on early life psychopathology show developmental variation. A recent review reported static and dynamic effects on internalizing and externalizing symptoms between infancy and early adulthood [Reference Hannigan, Walaker, Waszczuk, McAdams and Eley13], with certain genetic effects remaining influential throughout development, and novel genetic factors also emerging. However, previous studies focused on developmental periods rather than examining age continuously. Another way to test for developmental variation in genetic influences is Gene × Age (G × A) interactions. G × A interactions can be tested using a cross-sectional design that models differences in psychopathology as a function of relatedness and similarity in age between individuals [Reference Glahn, Kent, Sprooten, Diego, Winkler and Curran36Reference Almasy, Towne, Peterson and Blangero38]. Thus, G × A analysis tests for fluctuations in action of genetic factors, as well as variation in genetic factors at different ages [Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23].

We used data from the Philadelphia Neurodevelopmental Cohort (PNC) [Reference Calkins, Merikangas, Moore, Burstein, Behr and Satterthwaite1], a large population-based sample aged 8–21, to derive factors of psychopathology using a bifactor model of item-level data from a psychiatric interview [Reference Shanmugan, Wolf, Calkins, Moore, Ruparel and Hopson39,Reference Moore, Calkins, Satterthwaite, Roalf, Rosen and Gur40]. In bifactor models, items load on two factors simultaneously: (a) a general factor that accounts for commonality of all items (here general psychopathology) and (b) specific factors that account for unique influence of specific domains over and above the general factor (specific dimensions of psychopathology) [Reference Knowles, Mathias, Pearlson, Barrett, Mollon and Denbow41,Reference Chen, West and Sousa42]. Thus, bifactor models estimate the contribution of items to the general factor after controlling for specific factors, and vice versa [Reference Shanmugan, Wolf, Calkins, Moore, Ruparel and Hopson39]. Therefore, the utility of bifactor models lies in their ability to parse out this shared variance between general and specific factors, such that there is no contamination between factors, meaning that the general factor cannot account for findings in the specific factors, or vice versa. In other words, bifactor models allow examination of the unique contribution of the general and specific factors to prediction of external factors, or of the unique contribution of external factors (here genetic factors) to the general and specific factors [Reference Reise, Moore and Haviland43,Reference Reise44]. Importantly, bifactor models accommodate orthogonal factor scores despite correlated latent factors [Reference Shanmugan, Wolf, Calkins, Moore, Ruparel and Hopson39,Reference Reise, Moore and Haviland43]. We then used a genetic relatedness matrix to establish whether these psychopathology factors (a) were heritable, (b) showed genetic overlap with cognition, and (b) showed G × A interactions. In line with previous evidence, we hypothesized that psychopathology factors, would be (a) heritable, (b) show negative genetic correlations with cognitive functioning, and (c) be influenced by developmentally dynamic genetic factors, that is, show G × A interactions.

Methods

Participants

PNC is a population-based sample from the greater Philadelphia area, comprising 9,421 individuals aged 8–21. The study has been described in detail [Reference Calkins, Merikangas, Moore, Burstein, Behr and Satterthwaite1]. Briefly, between 2006 and 2012, 50,293 adults were recruited by the Center for Applied Genomics at Children’s Hospital of Philadelphia and provided access to Electronic Medical Records (EMRs). EMRs were screened for eligibility for PNC participation, yielding 19,161 individuals, released to the recruitment team in weekly waves between 2009 and 2011. Potential probands (ages 18–21) or caregivers/legal guardians (ages 8–17) were sent letters introducing the study, and then contacted by phone to explain the study, verify eligibility, and schedule appointments. Participants provided written consent for genomic studies upon providing blood samples during the clinical visit. Inclusion criteria were: (a) ability to provide signed informed consent (parental consent for participants <18), (b) English language proficiency, and (c) physical and cognitive ability to participate in cognitive testing. Data are in dbGaP (https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000607.v3.p2).

Genetic analyses were limited to participants who identified as white non-Hispanic (European American), leaving 4,662 subjects with genetic, cognitive, and psychiatric data. Mean age was 13.8 (standard deviation [SD] = 3.6), 50.3% were male (n = 2,346).

Cognitive assessment

Participants completed the Penn computerized neurocognitive battery [Reference Gur, Richard, Hughett, Calkins, Macy and Bilker45,Reference Moore, Reise, Gur, Hakonarson and Gur46], which consists of 14 tests that capture functioning in five domains: (a) executive function (abstraction and mental flexibility, attention, working memory), (b) episodic memory (verbal, facial, spatial), (c) complex cognition (verbal reasoning, nonverbal reasoning, spatial processing), (d) social cognition (emotion identification, emotion differentiation, age differentiation), and (e) speed (motor, sensorimotor). Accuracy and reaction times are recorded for each test. All tests show moderate to very high reliability [Reference Gur, Ragland, Moberg, Turner, Bilker and Kohler47]. The battery also included the reading subtest of the Wide Range Achievement Test (WRAT), a measure of general verbal knowledge. As in our prior work [Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23], we derived a general composite score (g) as the first component of principal component analysis (PCA) of accuracy scores. We also derived a general composite score for speed (gs) as the first component of PCA of reaction times. To minimize the impact of missing data, multivariate imputation by chained equation (MICE) [Reference White, Royston and Wood48] was used to impute missing values using the mice package in R [Reference van Buuren and Groothuis-Oudshoorn49]. Imputation was based on age, sex and available cognitive data (participants missing <50% cognitive data) [Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23]. Subsequent analyses were conducted on imputed data.

Psychopathology assessment

Psychiatric symptoms were ascertained using a computerized, structured interview (GOASSESS) [Reference Calkins, Merikangas, Moore, Burstein, Behr and Satterthwaite1,Reference Shanmugan, Wolf, Calkins, Moore, Ruparel and Hopson39,Reference Satterthwaite, Connolly, Ruparel, Calkins, Jackson and Elliott50], a modified version of the Kiddie-Schedule for Affective Disorders and Schizophrenia [Reference Kaufman, Birmaher, Brent, Rao, Flynn and Moreci51]. GOASSESS was administered to caregivers/legal guardians (ages 8–10), participants and caregivers/legal guardians (ages 11–17), and participants (ages 18–21). Bachelor- and Master-level assessors underwent a 25-h training protocol comprising didactic sessions, assigned readings, and supervised pairwise practice. Assessors were certified through standardized procedures requiring observation by a certified clinical observer who rated proficiency on a 60-item checklist of interview procedures. Responses coded by the assessor were required to correspond to responses coded by a certified clinical observer. Assessors underwent repeat observation until meeting passing criteria [Reference Calkins, Merikangas, Moore, Burstein, Behr and Satterthwaite1].

Factor analysis to create factors of psychopathology

We applied a confirmatory bifactor model [Reference Reise, Moore and Haviland43] in Mplus [Reference Muthén and Muthén52] to 112 items from the GOASSESS [Reference Calkins, Merikangas, Moore, Burstein, Behr and Satterthwaite1,Reference Shanmugan, Wolf, Calkins, Moore, Ruparel and Hopson39,Reference Moore, Calkins, Satterthwaite, Roalf, Rosen and Gur40](Figure 1a) using mean- and variance-adjusted weighted least squared estimator. Five orthogonal factors of psychopathology were generated for 9,421 individuals with GOASSESS data: (a) anxious-misery (mood and anxiety), (b) externalizing (ADHD and conduct disorder), (c) fear (phobias), (d) psychosis-spectrum, and (e) a general factor of overall psychopathology. Since bifactor models estimate the contribution of items to an overall dimension (general psychopathology) after controlling for specific factors, and vice versa, all factors (both general and specific) are orthogonal. Thus, bifactor models parse out the shared variance between general and specific factors, such that there is no contamination between factors. Therefore, bifactor models allow examination of the unique contribution of external factors (here genetic factors) to the general and specific factors [Reference Reise, Moore and Haviland43,Reference Reise44]. Table S1 shows factor loadings, Table S2 shows correlations between factors (and with cognition), Figure S1 shows test information plots.

Figure 1. Schematic of (a) bifactor and (b) hierarchical models of 112 items from the GOASSESS structured interview.

Genotyping

Samples were genotyped on one of four Illumina arrays: HumanHap550, HumanHap610, OmniExpress, or Human1M. Genotyped data were imputed at the Broad Institute [Reference Germine, Robinson, Smoller, Calkins, Moore and Hakonarson53] using IMPUTE2 and reference haplotypes in Phase I of the 1000 genomes data (June 2011 release) that included ~37,138,905 variants from 1,094 individuals from Africa, Asia, Europe, and the United States. Imputed genotype data were used in subsequent analyses.

Estimation of genetic relatedness matrix

Empirical relatedness quantifies the proportion of alleles that are identical by descent between individuals. Empirical relatedness was calculated for all pairs of individuals using genotype data [Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23]. Briefly, 50k common autosomal single nucleotide proteins (SNPs) in approximate linkage equilibrium were selected from all available SNP variants after linkage disequilibrium (LD) pruning (r 2 > 0.1) using PLINK [Reference Chang, Chow, Tellier, Vattikuti, Purcell and Lee54]. Relatedness was estimated from these SNPs using IBDLD [Reference Han and Abney55](up to 50 SNPs within a 2-cm span). The matrix was inspected to ensure correct properties (trace equal to number of genotyped subjects, symmetry, positive semi-definiteness, range of diagonal, and off-diagonal elements). Distribution of estimated relatedness values has been presented previously [Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23]. Empirical genetic relatedness matrices are advantageous because quantitative genetic analyses (described below) that are traditionally applied to family data using genetic relatedness matrices calculated from pedigree information can, in turn, be applied to cohorts of related and/or unrelated individuals (such as the PNC) using empirical genetic relatedness instead.

Statistical and quantitative genetic analyses

R [56] was used to generate descriptive statistics and graphics. Genetic analyses were conducted in Sequential Oligogenetic Linkage Analysis Routines (SOLAR) [Reference Blangero, Goring, Kent, Williams, Peterson and Almasy57] on 4,662 participants who identified as white non-Hispanic (European American), regardless of relatedness. While analyses in SOLAR performed on family data are robust to population stratification, the PNC sample comprises both related and unrelated individuals. Therefore, we only included individuals of European American ancestry (the most common ancestry group in the PNC sample) in our analyses and excluded individuals of non-European ancestry. Relatedly, since rare variants that may explain a substantial proportion of phenotypic variance are not well captured by common SNPs, using related individuals is more powerful than using unrelated individuals when estimating heritability, but the combination of related and unrelated individuals, as in the PNC, is optimal. Moreover, related individuals (even distantly) are critical for detecting G × A interactions (described below). When using only unrelated individuals it is not possible to detect changes in genetic correlation over time since related pairs serve as a pseudo-longitudinal design where the same polygenotypes are observed at different ages.

Univariate and bivariate polygenic models

SOLAR implements linear mixed-effects models, which decompose the overall variance of a quantitative trait [Reference Almasy and Blangero58,Reference Almasy and Blangero59]. Traditionally, these analyses are performed on family data using matrices calculated from pedigree information, but can be applied to cohorts of related and unrelated individuals using relatedness estimated from genotype data [Reference Speed, Hemani, Johnson and Balding60]. Under a univariate polygenic model, the phenotypic variance (σ2p) of a trait is decomposed into genetic (σ2g) and environmental (σ2e) components. Environmental variance incapsulates all variance that is not genetic, including error. Narrow-sense heritability (h2) is the proportion of phenotypic variance accounted for by additive genetic variance (h2 = σ2g/σ2p). To determine whether heritability (h2) was significantly greater than 0, likelihood of the polygenic model was compared to that of a model with h2 constrained to 0. Under a bivariate polygenic model, phenotypic covariance between two traits is decomposed into genetic and environmental components to determine the extent to which traits are influenced by shared genetic effects. Since genetic correlations between traits are only meaningful if traits are heritable, we estimated heritability of all traits. Bivariate polygenic analyses were then applied to significantly heritable pairs of traits to estimate genetic (ρg), environmental (ρe), and phenotypic (ρp) correlations. The genetic correlation (ρg) denotes the correlation between latent additive genetic factors influencing both traits. The environmental correlation (ρe) denotes the correlation between nongenetic factors influencing both traits. To determine whether genetic (ρg) and environmental (ρe) correlations were significantly different from 0, likelihood of the bivariate polygenic model was compared to that of a model where the parameter of interest was constrained to 0.

Gene × Age interaction models

A polygenic model can be extended to examine Gene × Environment (G × E) interactions [Reference Glahn, Kent, Sprooten, Diego, Winkler and Curran36Reference Almasy, Towne, Peterson and Blangero38]. One consequence of G × E is that additive genetic variance is greater under certain environments than others. To test for this effect with a continuous environmental variable (age), the polygenic model is modified to include a linear function on the logarithm of σ2g. This linear function contains a free parameter, γ, reflecting change in σ2g unit of the environmental variable (age in years). A nonzero value of γ implies a heritable response to the environment, and therefore, a G × E interaction. This G × E interaction tests for fluctuations (with age) in action of genetic factors and a significant G × E interaction suggests a change in magnitude of effect of specific genetic factors (with age). A second consequence of G × E is that the trait exhibits imperfect pleiotropy with itself, that is, the genetic correlation between the trait measured under one environment and the trait measured under another environment is less than 1. This phenomenon can be examined in cross-sectional studies where individuals are tested under a single environment (timepoint), provided relatedness between individuals is known [Reference Blangero37]. To test for this effect, the polygenic model is modified to include another free parameter, λ, reflecting the rate of decay in genetic correlation (ρg) as difference in the environmental variable increases. A nonzero value of λ implies imperfect pleiotropy, and therefore, a G × E interaction. This G × E interaction tests for variation in genetic factors influencing the trait (at different ages) and a significant G × E interaction suggests a change (with age) in the genetic factors themselves. G × E interaction models were fitted to heritable traits, with age in years as the continuous environmental variable that is, Gene × Age interactions. See Table S3 for more information.

All models included age, age2, sex, and their interactions as covariates. To adjust for multiple testing, false discovery rate (FDR) was set at 5% [Reference Benjamini and Yekutieli61]. Rank-based inverse normal transformations were applied to all traits to ensure normality.

Sensitivity analyses

In addition to the bifactor model described above, we applied a confirmatory hierarchical model (Figure 1b) in Mplus [Reference Muthén and Muthén52] to generate four correlated factors of psychopathology and a general factor. We generated the same factors as the bifactor model: specific factors of (a) anxious-misery, (b) externalizing, (c) fear, and (d) psychosis-spectrum, and (e) a general factor. In bifactor models, general and specific factors are orthogonal, whereas in hierarchical models, the general factor is defined by the specific factors and thus general and specific factors are explicitly correlated. Univariate, bivariate, and G × A analyses, as described above, were repeated on these factors to examine the contaminating effect, that is shared variance of general and specific factors. Table S4 shows factor loadings, Table S5 shows correlations between factors (and with cognition), Figure S2 shows test information plots.

Results

Externalizing psychopathology is heritable

Significant heritability estimates were observed for general (h2 = 0.21, p = 0.040) and externalizing psychopathology (h2 = 0.46, p = 2 × 10−6), but only externalizing survived FDR correction (Figure 2 and Table 1). As reported previously [Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23], most accuracy (h2 range = 0.21–0.72) and reaction time (h2 range = 0.23–0.38) measures were also heritable. Subsequent bivariate and G × A analyses were run on significantly heritable traits (after FDR correction).

Figure 2. Heritability estimates for all neurocognitive measures and psychopathology factors. *Error bars represent standard errors (SEs).

Table 1. Heritability estimates for all traits, genetic, and phenotypic correlations between externalizing and cognition.

Bolded estimates significant after correction for multiple testing (FDR = 0.05).

Abbreviations: AWRAT, Wide Range Achievement Test; SE, standard error.

Overlapping genetic factors on externalizing and cognition suggest pleiotropic effects

Externalizing showed significant negative phenotypic correlations with most accuracy measures (range ρp = −0.042 to −0.127) and some reaction time measures (range ρp = −0.036 to −0.66; Table 1) after FDR correction, such that greater psychopathology was associated with poorer accuracy and slower reaction times. Significant negative genetic correlations were observed after FDR correction between externalizing and accuracy measures of face memory (ρg = −0.412, p = 0.004), verbal reasoning (ρg = −0.485, p = 0.001), spatial processing (ρg = −0.426, p = 0.010), general verbal knowledge (ρg = −0.314, p = 0.002), g (ρg = −0.44, p = 0.002), and sensorimotor speed (ρg = −0.659, p = 1 × 10−4) suggesting that genetic factors underlying externalizing overlap with lower accuracy and slower reaction times.

Genetic variance on externalizing decreases with age

Significant decrease in genetic variance was observed on externalizing (γg = −0.146, p = 0.004; Figure 3), suggesting that specific genetic factors influence externalizing psychopathology between childhood and early adulthood, but also that the magnitude of effect of these genetic factors decreases with age. Significant increase in environmental variance was also observed (γe = 0.059, p = 0.009). Decay in genetic correlation did not reach statistical significance (λ = 0.027, p = 0.423), such that we did not find evidence for change in genetic factors, that is, novel genetic influences. Data presented in Figure 3 were generated using Formula 5 in Table S3.

Figure 3. Estimated genetic variance, environmental variance and heritability by age for externalizing.

Results of genetic analyses are robust to factor analytic approach

Univariate, bivariate, and G × A analyses conducted on psychopathology factors derived from the hierarchical model generated similar results. Externalizing (h2 = 0.58, p = 1x10−7) remained heritable, but general psychopathology (h2 = 0.37, p = 0.001) anxious-misery (h2 = 0.35, p = 0.002) and fear (h2 = 0.25, p = 0.011) were also significantly heritable, although fear did not survive FDR correction (Table S6 and Figure S3). General psychopathology and anxious-misery showed significant phenotypic correlations with cognition (Table S6), but only externalizing showed significant genetic correlations with verbal reasoning (ρg = −0.43, p = 0.001), general verbal knowledge (ρg = −0.25, p = 0.005), and g (ρg = −0.37, p = 3x10−6), after FDR correction. Again, significant decrease in genetic variance (γg = −0.112, p = 0.029) and increase in environmental variance (γe = 0.089, p = 0.008) were observed for externalizing (Figure S4) Change in genetic variance on general psychopathology (γg = 0.017, p = 0.279) and anxious-misery (γg = 0.023, p = 0.269) was not significant. Data presented in Figure S4 were generated using Formula 5 in Table S3.

Discussion

Using a large, population-based cohort of individuals aged 8–21, we showed that externalizing psychopathology in the first two decades of life is under considerable genetic influence. Externalizing showed genetic overlap with lower performance on face memory, verbal reasoning, spatial processing, motor speed, verbal knowledge, and general cognitive ability. We did not find evidence for novel genetic factors on externalizing throughout this developmental period, rather we found a decrease in genetic variance, and increase in environmental variance. These findings have several implications for our understanding of early life psychopathology.

First, to the best of our knowledge, this is the first study to use a large, population-based cohort and genetic relatedness matrix to estimate heritability of data-driven factors of psychopathology (both orthogonal and correlated) throughout childhood and early adulthood. Our finding of heritable general psychopathology, anxious-misery, externalizing, and fear are in line with previous evidence [Reference Polanczyk, Moffitt, Arseneault, Cannon, Ambler and Keefe11Reference Huizink, Van Den Berg, van der Ende and Verhulst15]. We did not find evidence for genetic influences on the psychosis-spectrum factor, in contrast to prior evidence [Reference Polanczyk, Moffitt, Arseneault, Cannon, Ambler and Keefe11,Reference Zavos, Freeman, Haworth, McGuire, Plomin and Cardno12]. However, SNP heritability estimates of psychotic symptoms are more modest, with a report of nonsignificant and zero estimates for paranoia and hallucinations, respectively [Reference Sieradzka, Power, Freeman, Cardno, Dudbridge and Ronald62]. Moreover, while there has been progress in delineating molecular genetic underpinnings of adult schizophrenia [Reference Ripke, O’Dushlaine, Chambert, Moran, Kähler and Akterin63], the same cannot be said of childhood psychotic symptoms [Reference Zammit, Hamshere, Dwyer, Georgiva, Timpson and Moskvina64]. Notably, subjects were less likely to endorse items that loaded on the psychosis-spectrum factor (average 11% endorsing), than on externalizing (21%), anxious-misery (15%), and fear (16%). Similarly, clinical diagnostic rates of attention deficit, oppositional defiant, and conduct disorders in PNC are 16, 33, and 7%, respectively [Reference Moore, Calkins, Satterthwaite, Roalf, Rosen and Gur40], while rate of threshold psychotic symptoms is 4% [Reference Calkins, Moore, Merikangas, Burstein, Satterthwaite and Bilker65]. Psychotic symptoms also emerge later during development than symptoms relating to externalizing, anxious-misery, and fear [Reference Häfner, Maurer, Löffler and Riecher-Rössler66,Reference Kessler, Berglund, Demler, Jin, Merikangas and Walters67].

In contrast, we found externalizing psychopathology to be significantly and moderately heritable, with genetic factors explaining 46–58% of the variance. Equally substantial genetic influences on externalizing between ages 8 and 26 have been reported in twin and adoption studies [Reference Hannigan, Walaker, Waszczuk, McAdams and Eley13Reference Huizink, Van Den Berg, van der Ende and Verhulst15]. A very similar SNP heritability of 0.44 has also been reported for externalizing [Reference Pappa, Fedko, Mileva-Seitz, Hottenga, Bakermans-Kranenburg and Bartels68], but another study found SNP heritability estimates of zero for a range of externalizing problems [Reference Trzaskowski, Dale and Plomin69]. Differences between SNP and pedigree heritability estimates are likely due to rare variants not well captured by common SNPs, with SNP heritability representing the lower bound. Thus, using related individuals is more powerful than using unrelated individuals when estimating heritability, and the combination of all possible relationships, as in our sample, results in something of a hybrid between SNP and pedigree estimates. However, shared environment may contribute to heritability inflation when using related individuals. Nevertheless, rare variants account for a significant proportion of total heritability [Reference Zaitlen, Kraft, Patterson, Pasaniuc, Bhatia and Pollack70], and our heritability estimates are generally in line with previous studies.

Importantly, we found significant heritability estimates for general psychopathology and anxious-misery factors generated from a hierarchical, but not bifactor, model. This finding highlights the utility of bifactor models, in which all factors (both general and specific) are orthogonal [Reference Moore, Calkins, Satterthwaite, Roalf, Rosen and Gur40,Reference Knowles, Carless, de Almeida, Curran, McKay and Sprooten71], thereby allowing us to examine the unique contribution of genetic factors to general psychopathology and specific factors of anxious-misery, externalizing, fear, and psychosis-spectrum [Reference Reise, Moore and Haviland43,Reference Reise44]. This finding, as well as the finding of stronger genetic influences on externalizing than general psychopathology, is also in line with previous evidence of genetic signal on specific cognitive factors emerging only when variance associated with general cognitive ability is parsed out, that is when applying a bifactor model [Reference Knowles, Carless, de Almeida, Curran, McKay and Sprooten71]. In hierarchical models, on the other hand, variance associated with the general factor may account for significant findings in the specific factors, and vice versa. Thus, our findings suggest that genetic effects underlying early-life anxiety and depression may underlie more general psychopathology. Estimates of genetic influence on childhood anxiety [Reference Gregory and Eley17,Reference Thapar and McGuffin18], depression [Reference Rice, Harold and Thapar19,Reference Birmaher, Ryan, Williamson, Brent, Kaufman and Dahl20], and phobias [Reference Lichtenstein and Annas72,Reference Kendler, Gardner, Annas, Neale, Eaves and Lichtenstein73] are heterogeneous, ranging from small to large (h 2 = 0.15–0.77) for anxiety and depression [Reference Gregory and Eley17,Reference Rice, Harold and Thapar19], and for phobias (h 2 = 0.23–0.72) [Reference Lichtenstein and Annas72,Reference Kendler, Gardner, Annas, Neale, Eaves and Lichtenstein73].

Second, we found negative genetic correlations between externalizing and face memory, verbal reasoning, spatial processing, sensorimotor speed, verbal knowledge, and general cognitive ability (g). These findings are in line with previous evidence of genetic overlap between ADHD and general and executive cognitive functions [Reference Bidwell, Willcutt, Defries and Pennington27Reference Savage, Jansen, Stringer, Watanabe, Bryois and de Leeuw32,Reference Donati, Dumontheil and Meaburn35], but expand knowledge by showing that genetic overlap extends to externalizing psychopathology more generally, as well as to memory, complex reasoning, and speed functions. Similarly, in the same sample, genetic overlap was reported between inattention and memory, social cognition, executive function, and complex reasoning [Reference Micalizzi, Brick, Marraccini, Benca-Bachman, Palmer and Knopik74]. Genetic overlap between externalizing and cognitive functioning may be due to the same genetic factors influencing both traits. A recent genome-wide meta-analyses identified shared risk loci for ADHD and intelligence [Reference O’Connell, Shadrin, Smeland, Bahrami, Frei and Bettella31]. Rare genetic variants have also been identified for cognition [Reference Huguet, Schramm, Douard, Tamer, Main and Monin75,Reference Kendall, Bracher-Smith, Fitzpatrick, Lynham, Rees and Escott-Price76] and ADHD [Reference Harich, van der Voet, Klein, Cizek, Fenckova and Schenck77], with overlapping genes implicated in ADHD studies of common and rare variants [Reference Martin, O’Donovan, Thapar, Langley and Williams78]. However, while there is evidence of common and rare variants disrupting similar biological pathways in ADHD [Reference Stergiakouli, Hamshere, Holmans, Langley, Zaharieva and Hawi79], the neurobiological mechanisms underlying shared genetic influences on ADHD and intelligence remain unclear [Reference O’Connell, Shadrin, Smeland, Bahrami, Frei and Bettella31]. Future studies that utilize animal models are needed to elucidate the causal and biological pathways underlying these shared genetic influences. Alternatively, cognitive impairments may lead to externalizing, and/or vice versa. Importantly, different cognitive measures show different psychometric properties and associations with specific functions require replication. Future longitudinal studies incorporating behavioral, neuroimaging, and genetic data can further disentangle these associations. Nevertheless, our findings suggest that cognitive impairment may be a useful endophenotype [Reference Glahn, Knowles, McKay, Sprooten, Raventós and Blangero80,Reference Gottesman and Gould81] of externalizing psychopathology. Interestingly, the WRAT, a measure of verbal knowledge showed the strongest phenotypic correlation with externalizing, but the weakest (statistically significant) genetic correlation. This finding highlights the importance of elucidating genetic underpinnings of phenotypic associations to delineate biological etiology.

Finally, we found static and dynamic genetic influences on externalizing psychopathology between childhood and early adulthood, in line with previous evidence [Reference Hannigan, Walaker, Waszczuk, McAdams and Eley13Reference Huizink, Van Den Berg, van der Ende and Verhulst15,Reference Kan, Dolan, Nivard, Middeldorp, van Beijsterveldt and Willemsen82]. While we did not find evidence for novel genetic influences, we found a decrease in genetic variance, and increase in environmental variance. Our findings are in line with those of Huizink [Reference Huizink, Van Den Berg, van der Ende and Verhulst15], who found a decrease in genetic influences on externalizing from 43 to 29% between age 12 and 26, as well as an increase in environmental influences from 39 to 52% between these ages. Similarly, Wichers [Reference Wichers, Gardner, Maes, Lichtenstein, Larsson and Kendler14] found a decrease in genetic effects from 78 to 73% and an increase in environmental effects from 20 to 26% between age 8 and 20, but also reported novel genetic influences throughout adolescence [Reference Wichers, Gardner, Maes, Lichtenstein, Larsson and Kendler14]. Huizink, on the other hand, reported novel environmental, but not genetic, influences [Reference Huizink, Van Den Berg, van der Ende and Verhulst15]. Several phenomena may underlie these G × A interactions. Genes may become less expressed due to maturational processes involving hormonal and physical development [Reference Hannigan, Walaker, Waszczuk, McAdams and Eley13]. Increasing environmental influences likely reflect growing peer influences and substance use [Reference Miettunen, Murray, Jones, Maki, Ebeling and Taanila83]. Of note, environmental variance in this study incapsulates all variance that is not genetic. However, measurement error is unlikely to account for our findings since reliability of our factors is high (Figures S1 and S2). Moreover, method of symptom reporting differed by age, and genetic effects may differ by reporting method. For example, Scourfield et al. [Reference Scourfield, Van den Bree, Martin and McGuffin84] found heritability estimates of 54 and 35% for parent- and self-reported conduct problems, respectively. Nevertheless, G × A interaction analyses adjusting additionally for reporting method generated similar findings, with significant decrease in genetic variance (γg = −0.151, p = 0.002; γg = −0.112, p = 0.026), and increase in environmental variance (γe = 0.058, p = 0.006; γe = 0.092, p = 0.007) on externalizing factors derived from bifactor and hierarchical models, respectively. Future studies that are able to combine longitudinal, self- and parent-report symptom data will help elucidate these age-associated effects further. Conversely, we previously reported increasing genetic and environmental variance on general cognitive ability in this sample [Reference Mollon, Knowles, Mathias, Gur, Peralta and Weiner23]. Thus, while a proportion of genetic factors underlying externalizing and cognition overlap, other, nonoverlapping genetic influences may show diverging developmental trajectories. Advanced quantitative genetic methods may shed light on trajectories of shared genetic influences.

This study has limitations. First, our analyses were restricted to European American individuals and future studies should include other populations. Second, our data were cross-sectional and longitudinal studies with repeated assessments of identical measures and individuals are needed to fully establish age-related changes in genetic factors. Third, lower heritability estimates for general psychopathology and anxious-misery meant less power to detect genetic correlations with cognition and G × A interactions. Finally, although we used a large sample, and comprehensive assessments of cognition and psychopathology, our findings require replication.

Acknowledgment

The authors thank all study participants.

Conflicts of Interest

The authors declare no conflicts of interest.

Financial Support

This research was supported by National Institute of Mental Health grants R01 MH107248 and MH107235.

Supplementary Materials

To view supplementary material for this article, please visit http://dx.doi.org/10.1192/j.eurpsy.2021.21.

Data Availability Statement

The data that support the findings of this study are openly available in the database of Genotypes and Phenotypes (dbGaP): https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000607.v3.p2.

References

Calkins, ME, Merikangas, KR, Moore, TM, Burstein, M, Behr, MA, Satterthwaite, TD, et al. The Philadelphia Neurodevelopmental Cohort: constructing a deep phenotyping collaborative. J Child Psychol Psychiatry. 2015;56:1356–69.CrossRefGoogle ScholarPubMed
Bora, E, Lin, A, Wood, SJ, Yung, AR, McGorry, PD, Pantelis, C. Cognitive deficits in youth with familial and clinical high risk to psychosis: a systematic review and meta-analysis. Acta Psychiatr Scand. 2014;130:115.CrossRefGoogle ScholarPubMed
Brewer, WJ, Wood, SJ, Phillips, LJ, Francey, SM, Pantelis, C, Yung, AR, et al. Generalized and specific cognitive performance in clinical high-risk cohorts: a review highlighting potential vulnerability markers for psychosis. Schizophr Bull. 2006;32:538–55.CrossRefGoogle Scholar
Deater-Deckard, K, Mullineaux, PY, Beekman, C, Petrill, SA, Schatschneider, C, Thompson, LA. Conduct problems, IQ, and household chaos: a longitudinal multi-informant study. J Child Psychol Psychiatry. 2009;50:1301–8.Google ScholarPubMed
Dionne, G, Tremblay, R, Boivin, M, Laplante, D, Pérusse, D. Physical aggression and expressive vocabulary in 19-month-old twins. Dev Psychol. 2003;39:261.Google ScholarPubMed
Seguin, JR, Parent, S, Tremblay, RE, Zelazo, PD. Different neurocognitive functions regulating physical aggression and hyperactivity in early childhood. J Child Psychol Psychiatry. 2009;50:679–87.CrossRefGoogle ScholarPubMed
Gur, RE, Moore, TM, Calkins, ME, Ruparel, K, Gur, RC. Face processing measures of social cognition: a dimensional approach to developmental psychopathology. Biol Psychiatry. 2017;2:502–9.Google ScholarPubMed
Lundy, SM, Silva, GE, Kaemingk, KL, Goodwin, JL, Quan, SF. Cognitive functioning and academic performance in elementary school children with anxious/depressed and withdrawn symptoms. Open Pediatr Med J. 2010;4:19.CrossRefGoogle ScholarPubMed
Wesseldijk, LW, Fedko, IO, Bartels, M, Nivard, MG, van Beijsterveldt, CE, Boomsma, DI, et al. Psychopathology in 7-year-old children: Differences in maternal and paternal ratings and the genetic epidemiology. Am J Med Genet B Neuropsychiatr Genet. 2017;174:251–60.Google ScholarPubMed
Polderman, TJ, Benyamin, B, de Leeuw, CA, Sullivan, PF, van Bochoven, A, Visscher, PM, et al. Meta-analysis of the heritability of human traits based on fifty years of twin studies. Nat Genet. 2015;47:702–9.CrossRefGoogle ScholarPubMed
Polanczyk, G, Moffitt, TE, Arseneault, L, Cannon, M, Ambler, A, Keefe, RS, et al. Etiological and clinical features of childhood psychotic symptoms: results from a birth cohort. Arch Gen Psychiatry. 2010;67:328–38.Google ScholarPubMed
Zavos, HM, Freeman, D, Haworth, CM, McGuire, P, Plomin, R, Cardno, AG, et al. Consistent etiology of severe, frequent psychotic experiences and milder, less frequent manifestations: a twin study of specific psychotic experiences in adolescence. JAMA Psychiatry. 2014;71:1049–57.Google ScholarPubMed
Hannigan, L, Walaker, N, Waszczuk, M, McAdams, T, Eley, T. Aetiological influences on stability and change in emotional and behavioural problems across development: a systematic review. Psychopathol Rev. 2017;4:pr. 038315.Google ScholarPubMed
Wichers, M, Gardner, C, Maes, H, Lichtenstein, P, Larsson, H, Kendler, K. Genetic innovation and stability in externalizing problem behavior across development: a multi-informant twin study. Behav Genet. 2013;43:191201.CrossRefGoogle ScholarPubMed
Huizink, AC, Van Den Berg, MP, van der Ende, J, Verhulst, FC. Longitudinal genetic analysis of internalizing and externalizing problem behavior in adopted biologically related and unrelated sibling pairs. Twin Res Human Genet. 2007;10:5565.Google ScholarPubMed
Ormel, J, Oldehinkel, AJ, Ferdinand, RF, Hartman, CA, De Winter, AF, Veenstra, R, et al. Internalizing and externalizing problems in adolescence: general and dimension-specific effects of familial loadings and preadolescent temperament traits. Psychol Med. 2005;35:1825–35.CrossRefGoogle ScholarPubMed
Gregory, AM, Eley, TC. The genetic basis of child and adolescent anxiety. Anxiety Disord Child Adolesc. 2011:161–78.CrossRefGoogle Scholar
Thapar, A, McGuffin, P. Are anxiety symptoms in childhood heritable? J Child Psychol Psychiatry. 1995;36:439–47.Google ScholarPubMed
Rice, F, Harold, G, Thapar, A. The genetic aetiology of childhood depression: a review. J Child Psychol Psychiatry. 2002;43:6579.CrossRefGoogle ScholarPubMed
Birmaher, B, Ryan, ND, Williamson, DE, Brent, DA, Kaufman, J, Dahl, RE, et al. Childhood and adolescent depression: a review of the past 10 years. Part I. J Am Acad Child Adolesc Psychiatry. 1996;35:1427–39.Google ScholarPubMed
Nivard, MG, Dolan, CV, Kendler, KS, Kan, KJ, Willemsen, G, van Beijsterveldt, CE, et al. Stability in symptoms of anxiety and depression as a function of genotype and environment: a longitudinal twin study from ages 3 to 63 years. Psychol Med. 2015;45:1039–49.CrossRefGoogle ScholarPubMed
Haworth, CM, Wright, MJ, Luciano, M, Martin, NG, de Geus, EJ, van Beijsterveldt, CE, et al. The heritability of general cognitive ability increases linearly from childhood to young adulthood. Mol Psychiatry. 2010;15:1112–20.CrossRefGoogle ScholarPubMed
Mollon, J, Knowles, EEM, Mathias, SR, Gur, R, Peralta, JM, Weiner, DJ, et al. Genetic influence on cognitive development between childhood and adulthood. Mol Psychiatry. 2021;26(2):656–65.Google ScholarPubMed
Jacobs, N, Rijsdijk, F, Derom, C, Danckaerts, M, Thiery, E, Derom, R, et al. Child psychopathology and lower cognitive ability: a general population twin study of the causes of association. Mol Psychiatry. 2002;7:368–74.CrossRefGoogle ScholarPubMed
Harden, KP, Engelhardt, LE, Mann, FD, Patterson, MW, Grotzinger, AD, Savicki, SL, et al. Genetic associations between executive functions and a general factor of psychopathology. J Am Acad Child Adolesc Psychiatry. 2020;59:749–58.CrossRefGoogle Scholar
Kuntsi, J, Eley, TC, Taylor, A, Hughes, C, Asherson, P, Caspi, A, et al. Co-occurrence of ADHD and low IQ has genetic origins. Am J Med Genet B Neuropsychiatr Genet. 2004;124B:41–7.Google ScholarPubMed
Bidwell, LC, Willcutt, EG, Defries, JC, Pennington, BF. Testing for neuropsychological endophenotypes in siblings discordant for attention-deficit/hyperactivity disorder. Biol Psychiatry. 2007;62:991–8.Google ScholarPubMed
Doyle, AE, Biederman, J, Seidman, LJ, Reske-Nielsen, JJ, Faraone, SV. Neuropsychological functioning in relatives of girls with and without ADHD. Psychol Med. 2005;35:1121–32.CrossRefGoogle ScholarPubMed
Wang, Z, Deater-Deckard, K, Petrill, SA, Thompson, LA. Externalizing problems, attention regulation, and household chaos: a longitudinal behavioral genetic study. Dev Psychopathol. 2012;24:755–69.CrossRefGoogle ScholarPubMed
Demontis, D, Walters, RK, Martin, J, Mattheisen, M, Als, TD, Agerbo, E, et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet. 2019;51:6375.Google ScholarPubMed
O’Connell, KS, Shadrin, A, Smeland, OB, Bahrami, S, Frei, O, Bettella, F, et al. Identification of genetic loci shared between attention-deficit/hyperactivity disorder, intelligence, and educational attainment. Biol Psychiatry. 2020;87(12):1052–62.CrossRefGoogle ScholarPubMed
Savage, JE, Jansen, PR, Stringer, S, Watanabe, K, Bryois, J, de Leeuw, CA, et al. Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence. Nat Genet. 2018;50:912–9.CrossRefGoogle ScholarPubMed
Liu, C-Y, Li, Y, Viding, E, Asherson, P, Pingault, J-B. The developmental course of inattention symptoms predicts academic achievement due to shared genetic aetiology: a longitudinal twin study. Eur Child Adolesc Psychiatry. 2019;28(3):367–75.CrossRefGoogle ScholarPubMed
Johnson, W, McGue, M, Iacono, WG. Disruptive behavior and school grades: genetic and environmental relations in 11-year-olds. J Edu Psychol. 2005;97(3):391405.Google Scholar
Donati, G, Dumontheil, I, Meaburn, EL. Genome-wide association study of latent cognitive measures in adolescence: genetic overlap with intelligence and education. Mind Brain Edu. 2019;13:224–33.Google ScholarPubMed
Glahn, DC, Kent, JW Jr, Sprooten, E, Diego, VP, Winkler, AM, Curran, JE, et al. Genetic basis of neurocognitive decline and reduced white-matter integrity in normal human brain aging. PNAS. 2013;110:19006–11.Google ScholarPubMed
Blangero, J. Statistical genetic approaches to human adaptability. Hum Biol. 1993;65:5.Google ScholarPubMed
Almasy, L, Towne, B, Peterson, C, Blangero, J. Detecting genotype× age interaction. Genet Epidemiol. 200;21:S81924.Google Scholar
Shanmugan, S, Wolf, DH, Calkins, ME, Moore, TM, Ruparel, K, Hopson, RD, et al. Common and dissociable mechanisms of executive system dysfunction across psychiatric disorders in youth. Am J Psychiatry. 2016;173:517–26.Google ScholarPubMed
Moore, TM, Calkins, ME, Satterthwaite, TD, Roalf, DR, Rosen, AFG, Gur, RC, et al. Development of a computerized adaptive screening tool for overall psychopathology (“p”). J Psychiatr Res. 2019;116:2633.Google Scholar
Knowles, EEM, Mathias, SR, Pearlson, GD, Barrett, J, Mollon, J, Denbow, D, et al. Clinical correlates of subsyndromal depression in African American individuals with psychosis: the relationship with positive symptoms and comorbid substance dependence. Schizophr Res. 2019;206:333–46.CrossRefGoogle ScholarPubMed
Chen, FF, West, SG, Sousa, KH. A comparison of bifactor and second-order models of quality of life. Multivariate Behav Res. 2006;41:189225.CrossRefGoogle ScholarPubMed
Reise, SP, Moore, TM, Haviland, MG. Bifactor models and rotations: exploring the extent to which multidimensional data yield univocal scale scores. J Personality Assess. 2010;92:544–59.CrossRefGoogle ScholarPubMed
Reise, SP. Invited paper: the rediscovery of bifactor measurement models. Multivariate Behav Res. 2012;47:667–96.CrossRefGoogle ScholarPubMed
Gur, RC, Richard, J, Hughett, P, Calkins, ME, Macy, L, Bilker, WB, et al. A cognitive neuroscience-based computerized battery for efficient measurement of individual differences: standardization and initial construct validation. J Neurosci Methods. 2010;187:254–62.CrossRefGoogle ScholarPubMed
Moore, TM, Reise, SP, Gur, RE, Hakonarson, H, Gur, RC. Psychometric properties of the Penn computerized neurocognitive battery. Neuropsychology. 2015;29:235.CrossRefGoogle ScholarPubMed
Gur, RC, Ragland, JD, Moberg, PJ, Turner, TH, Bilker, WB, Kohler, C, et al. Computerized neurocognitive scanning: I. Methodology and validation in healthy people. Neuropsychopharmacology. 2001;25:766–76.Google ScholarPubMed
White, IR, Royston, P, Wood, AM. Multiple imputation using chained equations: issues and guidance for practice. Stat Med. 2011;30:377–99.CrossRefGoogle ScholarPubMed
van Buuren, Sv, Groothuis-Oudshoorn, K. mice: multivariate imputation by chained equations in R. J Stat Softw. 2010;45(3):168.Google Scholar
Satterthwaite, TD, Connolly, JJ, Ruparel, K, Calkins, ME, Jackson, C, Elliott, MA, et al. The Philadelphia Neurodevelopmental Cohort: a publicly available resource for the study of normal and abnormal brain development in youth. Neuroimage. 2016;124:1115–9.CrossRefGoogle Scholar
Kaufman, J, Birmaher, B, Brent, D, Rao, U, Flynn, C, Moreci, P, et al. Schedule for affective disorders and schizophrenia for school-age children-present and lifetime version (K-SADS-PL): initial reliability and validity data. J Am Acad Child Adolesc Psychiatry. 1997;36:980–8.Google ScholarPubMed
Muthén, LK, & Muthén, BO. Mplus. Statistical analysis with latent variables. User’s guide. 2009;7.Google Scholar
Germine, L, Robinson, E, Smoller, J, Calkins, M, Moore, T, Hakonarson, H, et al. Association between polygenic risk for schizophrenia, neurocognition and social cognition across development. Transl Psychiatry. 2016;6:e924.CrossRefGoogle ScholarPubMed
Chang, CC, Chow, CC, Tellier, LC, Vattikuti, S, Purcell, SM, Lee, JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.Google ScholarPubMed
Han, L, Abney, M. Using identity by descent estimation with dense genotype data to detect positive selection. Eur J Hum Genet. 2013;21:205–11.CrossRefGoogle ScholarPubMed
R Development Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2008.Google Scholar
Blangero, J, Goring, HH, Kent, JW Jr, Williams, JT, Peterson, CP, Almasy, L, et al. Quantitative trait nucleotide analysis using Bayesian model selection. Human Biol. 2005;77:541–59.CrossRefGoogle ScholarPubMed
Almasy, L, Blangero, J. Variance component methods for analysis of complex phenotypes. Cold Spring Harb Protoc. 2010;2010:pdb.top77.Google ScholarPubMed
Almasy, L, Blangero, J. Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998;62:1198–211.CrossRefGoogle ScholarPubMed
Speed, D, Hemani, G, Johnson, MR, Balding, DJ. Improved heritability estimation from genome-wide SNPs. Am J Hum Genet. 2012;91:1011–21.Google ScholarPubMed
Benjamini, Y, Yekutieli, D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001:1165–88.Google Scholar
Sieradzka, D, Power, RA, Freeman, D, Cardno, AG, Dudbridge, F, Ronald, A. Heritability of individual psychotic experiences captured by common genetic variants in a community sample of adolescents. Behav Genet. 2015;45:493502.CrossRefGoogle Scholar
Ripke, S, O’Dushlaine, C, Chambert, K, Moran, JL, Kähler, AK, Akterin, S, et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet. 2013;45:1150.CrossRefGoogle Scholar
Zammit, S, Hamshere, M, Dwyer, S, Georgiva, L, Timpson, N, Moskvina, V, et al. A population-based study of genetic variation and psychotic experiences in adolescents. Schizophr Bull. 2014;40:1254–62.CrossRefGoogle ScholarPubMed
Calkins, ME, Moore, TM, Merikangas, KR, Burstein, M, Satterthwaite, TD, Bilker, WB, et al. The psychosis spectrum in a young US community sample: findings from the Philadelphia Neurodevelopmental Cohort. World Psychiatry. 2014;13:296305.CrossRefGoogle Scholar
Häfner, H, Maurer, K, Löffler, W, Riecher-Rössler, A. The influence of age and sex on the onset and early course of schizophrenia. Br J Psychiatry. 1993;162:80–6.CrossRefGoogle ScholarPubMed
Kessler, RC, Berglund, P, Demler, O, Jin, R, Merikangas, KR, Walters, EE. Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the National Comorbidity Survey Replication. Arch Gen Psychiatry. 2005;62:593602.CrossRefGoogle ScholarPubMed
Pappa, I, Fedko, IO, Mileva-Seitz, VR, Hottenga, JJ, Bakermans-Kranenburg, MJ, Bartels, M, et al. Single nucleotide polymorphism heritability of behavior problems in childhood: genome-wide complex trait analysis. J Am Acad Child Adolesc Psychiatry. 2015;54:737–44.Google ScholarPubMed
Trzaskowski, M, Dale, PS, Plomin, R. No genetic influence for childhood behavior problems from DNA analysis. J Am Acad Child Adolesc Psychiatry. 2013;52:1048–56 e3.CrossRefGoogle ScholarPubMed
Zaitlen, N, Kraft, P, Patterson, N, Pasaniuc, B, Bhatia, G, Pollack, S, et al. Using extended genealogy to estimate components of heritability for 23 quantitative and dichotomous traits. PLoS Genet. 2013;9:e1003520.Google ScholarPubMed
Knowles, EE, Carless, MA, de Almeida, MA, Curran, JE, McKay, DR, Sprooten, E, et al. Genome-wide significant localization for working and spatial memory: Identifying genes for psychosis using models of cognition. Am J Med Genet B Neuropsychiatr Genet. 2014;165B:8495.Google ScholarPubMed
Lichtenstein, P, Annas, P. Heritability and prevalence of specific fears and phobias in childhood. J Child Psychol Psychiatry. 2000;41:927–37.CrossRefGoogle ScholarPubMed
Kendler, KS, Gardner, CO, Annas, P, Neale, MC, Eaves, LJ, Lichtenstein, P. A longitudinal twin study of fears from middle childhood to early adulthood: evidence for a developmentally dynamic genome. Arch Gen Psychiatry. 2008;65:421–9.CrossRefGoogle ScholarPubMed
Micalizzi, L, Brick, LA, Marraccini, ME, Benca-Bachman, CE, Palmer, RHC, Knopik, VS. Single nucleotide polymorphism heritability and differential patterns of genetic overlap between inattention and four neurocognitive factors in youth. Dev Psychopathol. 2020:111.Google Scholar
Huguet, G, Schramm, C, Douard, E, Tamer, P, Main, A, Monin, P, et al. Genome-wide analysis of gene dosage in 24,092 individuals estimates that 10,000 genes modulate cognitive ability. Mol. Psychiatry. 2021.Google ScholarPubMed
Kendall, KM, Bracher-Smith, M, Fitzpatrick, H, Lynham, A, Rees, E, Escott-Price, V, et al. Cognitive performance and functional outcomes of carriers of pathogenic copy number variants: analysis of the UK Biobank. Br J Psychiatry. 2019;214:297304.CrossRefGoogle ScholarPubMed
Harich, B, van der Voet, M, Klein, M, Cizek, P, Fenckova, M, Schenck, A, et al. From rare copy number variants to biological processes in ADHD. Am J Psychiatry. 2020;177:855–66.Google ScholarPubMed
Martin, J, O’Donovan, MC, Thapar, A, Langley, K, Williams, N. The relative contribution of common and rare genetic variants to ADHD. Transl Psychiatry. 2015;5:e506.CrossRefGoogle ScholarPubMed
Stergiakouli, E, Hamshere, M, Holmans, P, Langley, K, Zaharieva, I, Hawi, Z, et al. Investigating the contribution of common genetic variants to the risk and pathogenesis of ADHD. Am J Psychiatry. 2012;169:186–94.Google ScholarPubMed
Glahn, DC, Knowles, EE, McKay, DR, Sprooten, E, Raventós, H, Blangero, J, et al. Arguments for the sake of endophenotypes: examining common misconceptions about the use of endophenotypes in psychiatric genetics. Am J Med Genet B Neuropsychiatr Genet. 2014;165:122–30.CrossRefGoogle Scholar
Gottesman, II, Gould, TD. The endophenotype concept in psychiatry: etymology and strategic intentions. Am J Psychiatry. 2003;160:636–45.CrossRefGoogle ScholarPubMed
Kan, KJ, Dolan, CV, Nivard, MG, Middeldorp, CM, van Beijsterveldt, CE, Willemsen, G, et al. Genetic and environmental stability in attention problems across the lifespan: evidence from the Netherlands twin register. J Am Acad Child Adolesc Psychiatry. 2013;52:1225.CrossRefGoogle ScholarPubMed
Miettunen, J, Murray, GK, Jones, PB, Maki, P, Ebeling, H, Taanila, A, et al. Longitudinal associations between childhood and adulthood externalizing and internalizing psychopathology and adolescent substance use. Psychol Med. 2014;44:1727–38.CrossRefGoogle ScholarPubMed
Scourfield, J, Van den Bree, M, Martin, N, McGuffin, P. Conduct problems in children and adolescents: a twin study. Arch Gen Psychiatry. 2004;61:489–96.CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Schematic of (a) bifactor and (b) hierarchical models of 112 items from the GOASSESS structured interview.

Figure 1

Figure 2. Heritability estimates for all neurocognitive measures and psychopathology factors. *Error bars represent standard errors (SEs).

Figure 2

Table 1. Heritability estimates for all traits, genetic, and phenotypic correlations between externalizing and cognition.

Figure 3

Figure 3. Estimated genetic variance, environmental variance and heritability by age for externalizing.

Supplementary material: File

Mollon et al. supplementary material

Mollon et al. supplementary material

Download Mollon et al. supplementary material(File)
File 3.7 MB
Submit a response

Comments

No Comments have been published for this article.