The Genetically Informed Neurobiology of Addiction (GINA) model.

Addictions are heritable and unfold dynamically across the lifespan. One prominent neurobiological theory proposes that substance-induced changes in neural circuitry promote the progression of addiction. Genome-wide association studies have begun to characterize the polygenic architecture undergirding addiction liability and revealed that genetic loci associated with risk can be divided into those associated with a general broad-spectrum liability to addiction and those associated with drug-specific addiction risk. In this Perspective, we integrate these genomic findings with our current understanding of the neurobiology of addiction to propose a new Genetically Informed Neurobiology of Addiction (GINA) model.

Reproducible brain-wide association studies require thousands of individuals.

Magnetic resonance imaging (MRI) has transformed our understanding of the human brain through well-replicated mapping of abilities to specific structures (for example, lesion studies) and functions1-3 (for example, task functional MRI (fMRI)). Mental health research and care have yet to realize similar advances from MRI. A primary challenge has been replicating associations between inter-individual differences in brain structure or function and complex cognitive or mental health phenotypes (brain-wide association studies (BWAS)). Such BWAS have typically relied on sample sizes appropriate for classical brain mapping4 (the median neuroimaging study sample size is about 25), but potentially too small for capturing reproducible brain-behavioural phenotype associations5,6. Here we used three of the largest neuroimaging datasets currently available-with a total sample size of around 50,000 individuals-to quantify BWAS effect sizes and reproducibility as a function of sample size. BWAS associations were smaller than previously thought, resulting in statistically underpowered studies, inflated effect sizes and replication failures at typical sample sizes. As sample sizes grew into the thousands, replication rates began to improve and effect size inflation decreased. More robust BWAS effects were detected for functional MRI (versus structural), cognitive tests (versus mental health questionnaires) and multivariate methods (versus univariate). Smaller than expected brain-phenotype associations and variability across population subsamples can explain widespread BWAS replication failures. In contrast to non-BWAS approaches with larger effects (for example, lesions, interventions and within-person), BWAS reproducibility requires samples with thousands of individuals.

GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis.

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Pleiotropy and genetically inferred causality linking multisite chronic pain to substance use disorders.

Individuals suffering from chronic pain develop substance use disorders (SUDs) more often than others. Understanding the shared genetic influences underlying the comorbidity between chronic pain and SUDs will lead to a greater understanding of their biology. Genome-wide association statistics were obtained from the UK Biobank for multisite chronic pain (MCP, Neffective = 387,649) and from the Million Veteran Program and the Psychiatric Genomics Consortium meta-analyses for alcohol use disorder (AUD, Neffective = 296,974), cannabis use disorder (CanUD, Neffective = 161,053), opioid use disorder (OUD, Neffective = 57,120), and problematic tobacco use (PTU, Neffective = 270,120). SNP-based heritability was estimated for each of the traits and genetic correlation (rg) analyses were performed to assess MCP-SUD pleiotropy. Bidirectional Mendelian Randomization analyses evaluated possible causal relationships. Finally, to identify and characterize individual loci, we performed a genome-wide pleiotropy analysis and a brain-wide analysis using imaging phenotypes available from the UK Biobank. MCP was positively genetically correlated with AUD (rg = 0.26, p = 7.55 × 10-18), CanUD (rg = 0.37, p = 8.21 × 10-37), OUD (rg = 0.20, p = 1.50 × 10-3), and PTU (rg = 0.29, p = 8.53 × 10-12). Although the MR analyses supported bi-directional relationships, MCP had larger effects on AUD (pain-exposure: beta = 0.18, p = 8.21 × 10-4; pain-outcome: beta = 0.07, p = 0.018), CanUD (pain-exposure: beta = 0.58, p = 2.70 × 10-6; pain-outcome: beta = 0.05, p = 0.014) and PTU (pain-exposure: beta = 0.43, p = 4.16 × 10-8; pain-outcome: beta = 0.09, p = 3.05 × 10-6) than the reverse. The genome-wide analysis identified two SNPs pleiotropic between MCP and all SUD investigated: IHO1 rs7652746 (ppleiotropy = 2.69 × 10-8), and CADM2 rs1248857 (ppleiotropy = 1.98 × 10-5). In the brain-wide analysis, rs7652746 was associated with multiple cerebellum and amygdala imaging phenotypes. When analyzing MCP pleiotropy with each SUD separately, we found 25, 22, and 4 pleiotropic variants for AUD, CanUD, and OUD, respectively. To our knowledge, this is the first large-scale study to provide evidence of potential causal relationships and shared genetic mechanisms underlying MCP-SUD comorbidity.

Alcohol milestones and internalizing, externalizing, and executive function: longitudinal and polygenic score associations.

BACKGROUND: Although the link between alcohol involvement and behavioral phenotypes (e.g. impulsivity, negative affect, executive function [EF]) is well-established, the directionality of these associations, specificity to stages of alcohol involvement, and extent of shared genetic liability remain unclear. We estimate longitudinal associations between transitions among alcohol milestones, behavioral phenotypes, and indices of genetic risk. METHODS: Data came from the Collaborative Study on the Genetics of Alcoholism (n = 3681; ages 11-36). Alcohol transitions (first: drink, intoxication, alcohol use disorder [AUD] symptom, AUD diagnosis), internalizing, and externalizing phenotypes came from the Semi-Structured Assessment for the Genetics of Alcoholism. EF was measured with the Tower of London and Visual Span Tasks. Polygenic scores (PGS) were computed for alcohol-related and behavioral phenotypes. Cox models estimated associations among PGS, behavior, and alcohol milestones. RESULTS: Externalizing phenotypes (e.g. conduct disorder symptoms) were associated with future initiation and drinking problems (hazard ratio (HR)⩾1.16). Internalizing (e.g. social anxiety) was associated with hazards for progression from first drink to severe AUD (HR⩾1.55). Initiation and AUD were associated with increased hazards for later depressive symptoms and suicidal ideation (HR⩾1.38), and initiation was associated with increased hazards for future conduct symptoms (HR = 1.60). EF was not associated with alcohol transitions. Drinks per week PGS was linked with increased hazards for alcohol transitions (HR⩾1.06). Problematic alcohol use PGS increased hazards for suicidal ideation (HR = 1.20). CONCLUSIONS: Behavioral markers of addiction vulnerability precede and follow alcohol transitions, highlighting dynamic, bidirectional relationships between behavior and emerging addiction.

Investigation of convergent and divergent genetic influences underlying schizophrenia and alcohol use disorder.

BACKGROUND: Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders. METHODS: We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific. RESULTS: We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001). CONCLUSIONS: Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

Characterizing Alcohol Expectancies in the ABCD Study: Associations with Sociodemographic Factors, the Immediate Social Environment, and Genetic Propensities.

Alcohol expectancies (AEs) are associated with likelihood of alcohol initiation and subsequent alcohol use disorders. It is unclear whether genetic predisposition to alcohol use and/or related traits contributes to shaping how one expects to feel when drinking alcohol. We used the Adolescent Brain Cognitive Development study to examine associations between genetic propensities (i.e., polygenic risk for problematic alcohol use, depression, risk-taking), sociodemographic factors (i.e., parent income), and the immediate social environment (i.e., peer use and disapproval toward alcohol) and positive and negative AEs in alcohol-naïve children (max analytic N = 5,352). Mixed-effect regression models showed that age, parental education, importance of the child's religious beliefs, adverse childhood experiences, and peer disapproval of alcohol use were associated with positive and/or negative AEs, to varying degrees. Overall, our results suggest several familial and psychosocial predictors of AEs but little evidence of contributions from polygenic liability to problematic alcohol use or related phenotypes.

A Phenome-Wide Association Study (PheWAS) of Late Onset Alzheimer Disease Genetic Risk in Children of European Ancestry at Middle Childhood: Results from the ABCD Study.

Genetic risk for Late Onset Alzheimer Disease (AD) has been associated with lower cognition and smaller hippocampal volume in healthy young adults. However, whether these and other associations are present during childhood remains unclear. Using data from 5556 genomically-confirmed European ancestry youth who completed the baseline session of the ongoing the Adolescent Brain Cognitive DevelopmentSM Study (ABCD Study®), our phenome-wide association study estimating associations between four indices of genetic risk for late-onset AD (i.e., AD polygenic risk scores (PRS), APOE rs429358 genotype, AD PRS with the APOE region removed (ADPRS-APOE), and an interaction between ADPRS-APOE and APOE genotype) and 1687 psychosocial, behavioral, and neural phenotypes revealed no significant associations after correction for multiple testing (all ps > 0.0002; all pfdr > 0.07). These data suggest that AD genetic risk may not phenotypically manifest during middle-childhood or that effects are smaller than this sample is powered to detect.

Genome-wide association study in individuals of European and African ancestry and multi-trait analysis of opioid use disorder identifies 19 independent genome-wide significant risk loci.

Despite the large toll of opioid use disorder (OUD), genome-wide association studies (GWAS) of OUD to date have yielded few susceptibility loci. We performed a large-scale GWAS of OUD in individuals of European (EUR) and African (AFR) ancestry, optimizing genetic informativeness by performing MTAG (Multi-trait analysis of GWAS) with genetically correlated substance use disorders (SUDs). Meta-analysis included seven cohorts: the Million Veteran Program, Psychiatric Genomics Consortium, iPSYCH, FinnGen, Partners Biobank, BioVU, and Yale-Penn 3, resulting in a total N = 639,063 (Ncases = 20,686;Neffective = 77,026) across ancestries. OUD cases were defined as having a lifetime OUD diagnosis, and controls as anyone not known to meet OUD criteria. We estimated SNP-heritability (h2SNP) and genetic correlations (rg). Based on genetic correlation, we performed MTAG on OUD, alcohol use disorder (AUD), and cannabis use disorder (CanUD). A leave-one-out polygenic risk score (PRS) analysis was performed to compare OUD and OUD-MTAG PRS as predictors of OUD case status in Yale-Penn 3. The EUR meta-analysis identified three genome-wide significant (GWS; p ≤ 5 × 10-8) lead SNPs-one at FURIN (rs11372849; p = 9.54 × 10-10) and two OPRM1 variants (rs1799971, p = 4.92 × 10-09; rs79704991, p = 1.11 × 10-08; r2 = 0.02). Rs1799971 (p = 4.91 × 10-08) and another OPRM1 variant (rs9478500; p = 1.95 × 10-08; r2 = 0.03) were identified in the cross-ancestry meta-analysis. Estimated h2SNP was 12.75%, with strong rg with CanUD (rg = 0.82; p = 1.14 × 10-47) and AUD (rg = 0.77; p = 6.36 × 10-78). The OUD-MTAG resulted in a GWAS Nequivalent = 128,748 and 18 independent GWS loci, some mapping to genes or gene regions that have previously been associated with psychiatric or addiction phenotypes. The OUD-MTAG PRS accounted for 3.81% of OUD variance (beta = 0.61;s.e. = 0.066; p = 2.00 × 10-16) compared to 2.41% (beta = 0.45; s.e. = 0.058; p = 2.90 × 10-13) explained by the OUD PRS. The current study identified OUD variant associations at OPRM1, single variant associations with FURIN, and 18 GWS associations in the OUD-MTAG. The genetic architecture of OUD is likely influenced by both OUD-specific loci and loci shared across SUDs.

Alcohol use and grey matter structure: Disentangling predispositional and causal contributions in human studies.

Alcohol use is a growing global health concern and economic burden. Alcohol involvement (i.e., initiation, use, problematic use, alcohol use disorder) has been reliably associated with broad spectrum grey matter differences in cross-sectional studies. These findings have been largely interpreted as reflecting alcohol-induced atrophy. However, emerging data suggest that brain structure differences also represent pre-existing vulnerability factors for alcohol involvement. Here, we review evidence from human studies with designs (i.e., family-based, genomic, longitudinal) that allow them to assess the plausibility that these correlates reflect predispositional risk factors and/or causal consequences of alcohol involvement. These studies provide convergent evidence that grey matter correlates of alcohol involvement largely reflect predisposing risk factors, with some evidence for potential alcohol-induced atrophy. These conclusions highlight the importance of study designs that can provide causal clues to cross-sectional observations. An integrative model may best account for these data, in which predisposition to alcohol use affects brain development, effects which may then be compounded by the neurotoxic consequences of heavy alcohol use.

Genome-wide association study of problematic opioid prescription use in 132,113 23andMe research participants of European ancestry.

The growing prevalence of opioid use disorder (OUD) constitutes an urgent health crisis. Ample evidence indicates that risk for OUD is heritable. As a surrogate (or proxy) for OUD, we explored the genetic basis of using prescription opioids 'not as prescribed'. We hypothesized that misuse of opiates might be a heritable risk factor for OUD. To test this hypothesis, we performed a genome-wide association study (GWAS) of problematic opioid use (POU) in 23andMe research participants of European ancestry (N = 132,113; 21% cases). We identified two genome-wide significant loci (rs3791033, an intronic variant of KDM4A; rs640561, an intergenic variant near LRRIQ3). POU showed positive genetic correlations with the two largest available GWAS of OUD and opioid dependence (rg = 0.64, 0.80, respectively). We also identified numerous additional genetic correlations with POU, including alcohol dependence (rg = 0.74), smoking initiation (rg = 0.63), pain relief medication intake (rg = 0.49), major depressive disorder (rg = 0.44), chronic pain (rg = 0.42), insomnia (rg = 0.39), and loneliness (rg = 0.28). Although POU was positively genetically correlated with risk-taking (rg = 0.38), conditioning POU on risk-taking did not substantially alter the magnitude or direction of these genetic correlations, suggesting that POU does not simply reflect a genetic tendency towards risky behavior. Lastly, we performed phenome- and lab-wide association analyses, which uncovered additional phenotypes that were associated with POU, including respiratory failure, insomnia, ischemic heart disease, and metabolic and blood-related biomarkers. We conclude that opioid misuse can be measured in population-based cohorts and provides a cost-effective complementary strategy for understanding the genetic basis of OUD.

Characterisation of the genetic relationship between the domains of sleep and circadian-related behaviours with substance use phenotypes.

Sleep problems and substance use frequently co-occur. While substance use can result in specific sleep deficits, genetic pleiotropy could explain part of the relationship between sleep and substance use and use disorders. Here we use the largest publicly available genome-wide summary statistics of substance use behaviours (N = 79,729-632,802) and sleep/activity phenotypes to date (N = 85,502-449,734) to (1) assess the genetic overlap between substance use behaviours and both sleep and circadian-related activity measures, (2) estimate clusters from genetic correlations and (3) test processes of causality versus genetic pleiotropy. We found 31 genetic correlations between substance use and sleep/activity after Bonferroni correction. These patterns of overlap were represented by two genetic clusters: (1) tobacco use severity (age of first regular tobacco use and smoking cessation) and sleep health (sleep duration, sleep efficiency and chronotype) and (2) substance consumption/problematic use (drinks per day and cigarettes per day, cannabis use disorder, opioid use disorder and problematic alcohol use) and sleep problems (insomnia, self-reported short sleep duration, increased number of sleep episodes, increased sleep duration variability and diurnal inactivity) and measures of circadian-related activity (L5, M10 and sleep midpoint). Latent causal variable analyses determined that horizontal pleiotropy (rather than genetic causality) underlies a majority of the associations between substance use and sleep/circadian related measures, except one plausible genetically causal relationship for opioid use disorder on self-reported long sleep duration. Results show that shared genetics are likely a mechanism that is at least partly responsible for the overlap between sleep and substance use traits.

Etiological Overlap Between Sex Under the Influence and Number of Lifetime Sexual Partners.

Drug and alcohol use is associated with risky sexual behavior (RSB). It is unclear whether this association is due to correlated liabilities (e.g., third variables influencing both traits), or whether use of drugs and alcohol during sexual decision making increases RSB. This study addresses this question by fitting a series of biometrical models using over 800 twin pairs assessed in early adulthood (m = 25.21 years). Measures included an index of sex under the influence (e.g., frequency that drugs or alcohol affect sexual decision making), number of lifetime sexual partners, and a general measure of substance use. Analyses suggest the covariance among these measures is explained by both genetic and environmental correlated liabilities. The overlap was not specific to sex under the influence, but was shared with a measure of general substance use. Models testing necessary but not sufficient parameters for direction of causation suggest that sex under the influence is unlikely to cause an increase in RSB; more evidence for reverse causation was found.

Genes identified in rodent studies of alcohol intake are enriched for heritability of human substance use.

BACKGROUND: Rodent paradigms and human genome-wide association studies (GWAS) on drug use have the potential to provide biological insight into the pathophysiology of addiction. METHODS: Using GeneWeaver, we created rodent alcohol and nicotine gene-sets derived from 19 gene expression studies on alcohol and nicotine outcomes. We partitioned the SNP heritability of these gene-sets using four large human GWAS: (1) alcoholic drinks per week, (2) problematic alcohol use, (3) cigarettes per day, and (4) smoking cessation. We benchmarked our findings with curated human alcohol and nicotine addiction gene-sets and performed specificity analyses using other rodent gene-sets (e.g., locomotor behavior) and other human GWAS (e.g., height). RESULTS: The rodent alcohol gene-set was enriched for heritability of drinks per week, cigarettes per day, and smoking cessation, but not problematic alcohol use. However, the rodent nicotine gene-set was not significantly associated with any of these traits. Both rodent gene-sets showed enrichment for several non-substance-use GWAS, and the extent of this relationship tended to increase as a function of trait heritability. In general, larger gene-sets demonstrated more significant enrichment. Finally, when evaluating human traits with similar heritabilities, both rodent gene-sets showed greater enrichment for substance use traits. CONCLUSION: Our results suggest that rodent gene expression studies can help to identify genes that contribute to the heritability of some substance use traits in humans, yet there was less specificity than expected. We outline various limitations, interpretations, and considerations for future research.

Genetic and Environmental Influence on the Human Functional Connectome.

Detailed mapping of genetic and environmental influences on the functional connectome is a crucial step toward developing intermediate phenotypes between genes and clinical diagnoses or cognitive abilities. We analyzed resting-state functional magnetic resonance imaging data from two adult twin samples (Nos = 446 and 371) to quantify genetic and environmental influence on all pairwise functional connections between 264 brain regions (~35 000 functional connections). Nonshared environmental influence was high across the whole connectome. Approximately 14-22% of connections had nominally significant genetic influence in each sample, 4.6% were significant in both samples, and 1-2% had heritability estimates greater than 30%. Evidence of shared environmental influence was weak. Genetic influences on connections were distinct from genetic influences on a global summary measure of the connectome, network-based estimates of connectivity, and movement during the resting-state scan, as revealed by a novel connectome-wide bivariate genetic modeling procedure. The brain's genetic organization is diverse and not as one would expect based solely on structure evident in nongenetically informative data or lower resolution data. As follow-up, we make novel classifications of functional connections and examine highly localized connections with particularly strong genetic influence. This high-resolution genetic taxonomy of brain connectivity will be useful in understanding genetic influences on brain disorders.

Exploring the genetic overlap of suicide-related behaviors and substance use disorders.

Suicide-related behaviors are heterogeneous and transdiagnostic, and may demonstrate varying levels of genetic overlap with different substance use disorders (SUDs). We used linkage disequilibrium score regression, genomic structural equation models, and Mendelian randomization to examine the genetic relationships between several SUDs and suicide-related behaviors. Our analyses incorporated summary statistics from the largest genome-wide association studies (GWAS) of problematic alcohol use, the Fagerström test for nicotine dependence, cannabis use disorder, and opioid use disorder (Ns ranging from 46,213-435,563) and GWAS of ever self-harmed, suicide attempt, and suicide death (Ns ranging from 18,223-117,733). We also accounted for genetic liability to depression (N = 500,199) and risk tolerance (N = 315,894). Suicide-related behaviors were significantly genetically correlated with each other and each SUD, but there was little evidence of causal relationships between the traits. Simultaneously correlating a common SUD factor with each specific suicide indicator while controlling for depression and risk tolerance revealed significant, positive genetic correlations between the SUD factor and suicide-related behaviors (rg  = 0.26-0.45, SE = 0.08-0.09). In the model, depression's association with suicide death (ß = 0.42, SE = 0.06) was weaker compared to ever-self harmed and suicide attempt (ß = 0.58, SE = 0.05 and ß = 0.50, SE = 0.06, respectively). We identify a general level of genetic overlap between SUDs and suicide-related behaviors, which is independent of depression and risk tolerance. Additionally, our findings suggest that genetic and behavioral contributions to suicide death may somewhat differ from nonlethal suicide-related behaviors.

The Speed of Progression to Tobacco and Alcohol Dependence: A Twin Study.

We investigated the etiological role of genetic and environmental influences for two milestones of tobacco and alcohol use: age of initiation, and speed of progression to dependence (latency). Study participants included 1352 monozygotic and 1422 dizygotic twins (mean age at assessment = 24.31). Earlier ages of initiation significantly increased the likelihood of developing dependence, but were associated with longer dependence latencies for tobacco and alcohol. Latencies to dependence were heritable traits for tobacco (a2 = 0.63) and alcohol (a2 = 0.64). Genetic influences contributing to early age of initiation were associated with faster latencies to dependence but sometimes were counteracted by environmental factors, the extent to which depended on substance and, sometimes, sex. Our findings may have important implications for public policy and add to the literature by characterizing the genetic and environmental contributions to the speed of progression to tobacco and alcohol dependence.

Etiology of Stability and Growth of Internalizing and Externalizing Behavior Problems Across Childhood and Adolescence.

Internalizing and externalizing behaviors are heritable, and show genetic stability during childhood and adolescence. Less work has explored how genes influence individual differences in developmental trajectories. We estimated ACE biometrical latent growth curve models for the Teacher Report Form (TRF) and parent Child Behavior Checklist (CBCL) internalizing and externalizing scales from ages 7 to 16 years in 408 twin pairs from the Colorado Longitudinal Twin Study. We found that Intercept factors were highly heritable for both internalizing and externalizing behaviors (a2 = .61-.92), with small and nonsignificant environmental influences for teacher-rated data but significant nonshared environmental influences for parent-rated data. There was some evidence of heritability of decline in internalizing behavior (Slopes for teacher and parent ratings), but the Slope genetic variance was almost entirely shared with that for the Intercept when different than zero. These results suggest that genetic effects on these developmental trajectories operate primarily on initial levels and stability, with no significant unique genetic influences for change. Finally, cross-rater analyses of the growth factor scores revealed moderate to large genetic and environmental associations between growth factors derived from parents' and teachers' ratings, particularly the Intercepts.

Do executive functions explain the covariance between internalizing and externalizing behaviors?

This study examined whether executive functions (EFs) might be common features of internalizing and externalizing behavior problems across development. We examined relations between three EF latent variables (a common EF factor and factors specific to updating working memory and shifting sets), constructed from nine laboratory tasks administered at age 17, to latent growth intercept (capturing stability) and slope (capturing change) factors of teacher- and parent-reported internalizing and externalizing behaviors in 885 individual twins aged 7 to 16 years. We then estimated the proportion of intercept-intercept and slope-slope correlations predicted by EF as well as the association between EFs and a common psychopathology factor (P factor) estimated from all 9 years of internalizing and externalizing measures. Common EF was negatively associated with the intercepts of teacher-rated internalizing and externalizing behavior in males, and explained 32% of their covariance; in the P factor model, common EF was associated with the P factor in males. Shifting-specific was positively associated with the externalizing slope across sex. EFs did not explain covariation between parent-rated behaviors. These results suggest that EFs are associated with stable problem behavior variation, explain small proportions of covariance, and are a risk factor that that may depend on gender.