Using both twin and molecular genetic strategies, we have learned a great deal in recent years about the relationship between the genetic liabilities to major psychiatric and substance use disorders (Consortium, 2013; Cross-Disorder Group of the Psychiatric Genomics Consortium, 2019; Kendler et al., Reference Kendler, Aggen, Knudsen, Roysamb, Neale and Reichborn-Kjennerud2011). We have also begun to clarify the associations between genetic liabilities to our key disorders and genetic propensity for educational attainment (Bulik-Sullivan et al., Reference Bulik-Sullivan, Finucane, Anttila, Gusev, Day, Loh and Neale2015; Lee et al., Reference Lee, Wedow, Okbay, Kong, Maghzian, Zacher and Cesarini2018; Rosoff, Kaminsky, McIntosh, Smith, & Lohoff, Reference Rosoff, Kaminsky, McIntosh, Smith and Lohoff2020). Considering the genetic aptitude for educational attainment (GAEA) and the genetic risks for alcohol use disorder (AUD) and drug use disorder (DUD), heritability has been shown to be substantial for all three phenotypes from twin models (Baker, Treloar, Reynolds, Heath, & Martin, Reference Baker, Treloar, Reynolds, Heath and Martin1996; Branigan, McCallum, & Freese, Reference Branigan, McCallum and Freese2013; Kendler et al., Reference Kendler, Ji, Edwards, Ohlsson, Sundquist and Sundquist2015a; Tambs, Sundet, Magnus, & Berg, Reference Tambs, Sundet, Magnus and Berg1989; Tsuang et al., Reference Tsuang, Lyons, Eisen, Goldberg, True, Lin and Eaves1996; Verhulst, Neale, & Kendler, Reference Verhulst, Neale and Kendler2015). High educational attainment, socio-economic status (SES), and intelligence are consistently associated with lower levels of alcohol consumption and/or problematic drinking in Sweden (Sjölund, Hemmingsson, & Allebeck, Reference Sjölund, Hemmingsson and Allebeck2015a; Sjölund, Hemmingsson, Gustafsson, & Allebeck, Reference Sjölund, Hemmingsson, Gustafsson and Allebeck2015b; Wennberg, Andersson, & Bohman, Reference Wennberg, Andersson and Bohman2002; Zettergren & Bergman, Reference Zettergren and Bergman2014), and similar trends are typically seen elsewhere (Barr, Silberg, Dick, & Maes, Reference Barr, Silberg, Dick and Maes2018; Müller et al., Reference Müller, Kowalewski, Metzler, Stettbacher, Rössler and Vetter2013; Rogne, Pedersen, & Von Soest, Reference Rogne, Pedersen and Von Soest2021). In two twin studies, the heritability of alcohol consumption was moderated by SES such that heritability was lower in subjects with higher SES (Davis & Slutske, Reference Davis and Slutske2018; Hamdi, Krueger, & South, Reference Hamdi, Krueger and South2015).
We seek to expand on these two twin studies in four important ways. First, we examined AUD rather than alcohol consumption. Second, instead of using an overt measure such as SES, we instead look at the genetic propensity to high educational attainment. This permits us to model directly the interaction at a genetic level, an approach rarely implemented in the past. Third, we added a second substance use disorder, DUD, to see if the findings generalize. Finally, instead of using twin models, we take advantage of a recently developed method to estimate genetic risks or aptitudes from extended pedigrees in a Swedish national sample – the family genetic risk score (FGRS) (Kendler et al., Reference Kendler, Ohlsson, Mościcki, Sundquist, Edwards and SundquistIn press; Kendler, Ohlsson, Sundquist, & Sundquist, Reference Kendler, Ohlsson, Sundquist and Sundquist2021a, Reference Kendler, Ohlsson, Sundquist and Sundquist2021b). So, our specific question is whether the genetic risks for AUD and DUD are moderated by the GAEA. In line with prior work, we predict that, if moderation occurs, high levels of GAEA will attenuate the impact of genetic liability on the risk for AUD or DUD.
Methods
We collected information on individuals from Swedish population-based registers with national coverage linking each person's unique personal identification number, which for confidentiality, was replaced with a serial number by Statistics Sweden. This study was approved by the Regional Ethical Review Board of Lund (No. 2008/409, 2012/795, and 2016/679).
Our database consisted of all individuals born in Sweden from 1960 to 1980 of parents themselves born in Sweden. In the database, we included date of first registration for AUD and DUD utilizing ICD-8, 9, 10 codes from Swedish national primary care, specialist, and hospital registries as well as criminal registers and prescribed drug registers. We also included individual educational attainment (see Appendix Table 1 for full definitions). In the database, individual FGRS for the two disorders as well as a GAEA were included. Similar to prior studies (Kendler et al., Reference Kendler, Ohlsson, Mościcki, Sundquist, Edwards and SundquistIn press; Kendler et al., Reference Kendler, Ohlsson, Sundquist and Sundquist2021a, Reference Kendler, Ohlsson, Sundquist and Sundquist2021b), the FGRS and GAEA were based on selected 1st, 2nd, 3rd, 4th, and 5th degree relatives to the probands with a mean of 40.1 relatives per proband. Briefly (see Appendix Table 2 for details), we first calculated the morbid risk for the phenotype in our sample of relatives based on age at first registration and then we transformed the binary trait into an underlying liability distribution, with the threshold that divides the population into those unaffected and affected for the disorder. Thereafter, we calculated the mean z-score for relatives with the disorder and the mean z-score for individuals without the disorder. For first-degree relatives, we also multiplied the z-score with a factor that sought to subtract the influence of cohabitation separately for siblings and parent–offspring pairs, thereby helping to remove from the FGRS any residual shared-environmental effect. For parent–offspring pairs, this was calculated by comparing the resemblance, by logistic regression, for father–offspring pairs where the biological father sired and raised his child (i.e. a father in an intact family) to the resemblance between children and their not-lived-with fathers, that is, those who sired their offspring but never lived with or near them when they were growing up. We have examined such not-lived-with fathers in several prior extended adoption studies (e.g. Kendler et al., Reference Kendler, Ji, Edwards, Ohlsson, Sundquist and Sundquist2015a Reference Kendler, Ohlsson, Maes, Sundquist, Lichtenstein and Sundquist2015b). For sibling pairs, we compared the resemblance in half-sibs who were v. were not reared together. As seen in Table 2 step 3 in the Appendix, the correction factors – the degree of the resemblance for our individual diagnoses that was retained after discounting the effect of shared environment equaled 0.99 (AUD) and 0.92 (DUD) for parent–offspring pairs and 0.69 (AUD) and 0.52 (DUD) for sibling pairs. Within each type of relative, we then had two components: the sum of the z-score and the total weighted number of relatives. These two components were weighted according to the genetic resemblance to the proband. For each proband, we summed the two components across all groups of relatives and used the quotient between the two components. Finally, to obtain the individual FGRS, we multiplied the quotient with a shrinkage factor based on the variance of the z-score across all relatives, the variance in the mean z-score across all probands, and the number of weighted number of relatives for each proband. So that the FGRSs would be more comparable across traits and to reduce the effect of register coverage, we standardized the FGRS by year of birth and county of residence into a z-score with mean = 0 and s.d. = 1. The calculation for GAEA followed the same principles but was simpler as we summed the weighted (by genetic resemblance) educational attainment across all relatives. The educational attainment was discounted for the effect of shared environment by a factor of 0.82 for parents and 0.73 for siblings. The weighted sum was then further weighted by number of relatives and standardized as described above.
To investigate whether GAEA moderates the impact of the familial genetic risk, we used Aalen's linear hazards model (Aalen, Reference Aalen1989). Follow-up time in months was measured from age 15 until time of first registration for AUD (DUD), death, emigration, or end of follow-up (31 December 2017). In the first model, aside from year of birth and sex, we include FGRSAUD (DUD) and the GAEA as well as their interaction. The results from this model are presented as the excess number of cases per 10 000 person-years, and the interaction is measured on the additive scale as we have defended previously (Kendler & Gardner, Reference Kendler and Gardner2010). To investigate the degree to which the interaction is a result of the effect of GAEA on educational attainment, we fitted three additional models, where the follow-up time started at age 25 so as to increase the chance that individuals had largely completed their education. Individuals who were registered for AUD or DUD prior to age 25 were excluded. Model A is the same model as described above, while model B also includes a main effect for educational attainment and model C further includes two additional two-way interaction terms [between educational attainment and FGRSAUD (FGRSDUD), and between educational attainment and GAEA] as well as a three-way interaction term between FGRSAUD (FGRSDUD), GAEA, and educational attainment. In all models, we investigated the proportionality assumption and there were no major violations. Statistical analyses were performed using SAS statistical software, version 9.4 (SAS Institute, 2012) and the R-package Timereg in R (Martinussen & Scheike, Reference Martinussen and Scheike2006; Scheike & Zhang, Reference Scheike and Zhang2011; Team, 2020).
Results
Descriptive statistics
Table 1 provides descriptive statistics for the AUD and DUD. AUD had a modestly higher population prevalence (5.4%) than DUD (3.4%), with the expected male excess in both disorders. In affected individuals, the genetic risk for each of the disorders was more strongly elevated for DUD (+0.72 s.d.) than for AUD (+0.58 s.d.). The average GAEA was substantially below the population mean for individuals affected with DUD (−0.40 s.d.) and AUD (−0.34 s.d.).
Model fitting: primary models
We first predicted risk for AUD controlling for sex and year of birth (Table 2). The model contained both the linear effects of the FGRS for AUD and the GAEA and their interaction. As expected, the FGRS for AUD had a strong and substantial effect on AUD risk while GAEA had a protective effect, approximately half as strong. Most importantly, we found a robust negative interaction between the FGRSAUD and the GAEA (Fig. 1a). For those with a high GAEA [light blue line (very light grey) at the bottom of the figure], the slope of increased risk with increasing FGRSAUD was relatively flat – that is, a large increase in genetic risk produced only a modest increase in affected cases. However, this slope got steeper as the GAEA got lower. At the lowest level of GAEA, the same change in FGRSAUD risk produced more than three times as many affected individuals as seen with those at the highest level of GAEA.
Fig. 1. (a) The interaction effects in the prediction of alcohol use disorder (AUD) between the family genetic risk score for alcohol use disorder and the genetic aptitude for educational attainment. The x-axis is the level of the family genetic risk score for AUD in standard deviation units. The y-axis is the excess number of AUD cases predicted per 10 000 person years. The colored (grey-scale) lines reflect the level of the genetic aptitude for educational attainment in standard deviation units. For example, the light blue (very light grey) line at the bottom reflects a quite high genetic aptitude for educational attainment (+2 s.d.) while the dark blue line (very dark grey) at the top of the figure reflects a quite low genetic aptitude for educational attainment (−2 s.d.). (b) The interaction effects in the prediction of drug use disorder (DUD) between the family genetic risk score for DUD and the genetic aptitude for educational attainment. The x-axis is the level of the family genetic risk score for DUD in standard deviation units. The y-axis is the excess number of DUD cases predicted per 10 000 person years. The colored lines reflect the level of the genetic aptitude for educational attainment in standard deviation units. For example, the light blue (very light grey) line at the bottom reflects a quite high genetic aptitude for educational attainment (+2 s.d.) while the dark blue line (very dark grey) at the top of the figure reflects a quite low genetic aptitude for educational attainment (−2 s.d.).
For DUD, both the main effect of the FGRSDUD and the GAEA were slightly smaller than seen with AUD as was the interaction between them (Table 2) but the pattern of their interaction, as illustrated in Fig. 1b, was very similar to that seen for AUD.
Model fitting with educational attainment
GAEA appeared to strongly moderate the impact of the genetic risk for AUD and DUD. However, this moderation by GAEA could be largely driven by the impact of GAEA on educational attainment itself. To investigate this, we fitted three models to our original sample from which we censored individuals with a registration of AUD/DUD prior to <age 25 to insure they had achieved their final level of educational attainment (Table 3). Model A is the same model fitted to our entire sample in Table 2. The results were similar although the strength of the interaction effects between the FGRS and GAEA declined for both AUD and DUD suggesting that these interactions are strongest in the early onset cases censored from this sample.
Model B included a main effect for educational attainment which robustly predicted risk for AUD and DUD. Importantly, the magnitudes of the interaction between GAEA and the FGRSAUD, and FGRSDUD observed in this model were nearly identical to those seen in model A. These results indicate that little of the observed interaction between GAEA and the genetic risks for AUD and DUD was mediated through actual educational attainment.
Model C contained all the parameters in model B to which we added a three-way interaction between the relevant FGRS, GAEA, and educational attainment. This interaction was significant for both AUD and DUD. We illustrate this interaction in Fig. 2a and b presenting the results obtained for individuals with high educational attainment (+1 s.d.), mean educational attainment, and very low educational attainment (−2 s.d.). The pattern was similar across both AUD and DUD. The largest interaction between GAEA and FGRS was seen in those with very low educational attainment. The interaction was present but attenuated for those who obtained a mean level of education. However, for those at high educational attainment, the interaction effects disappeared.
Fig. 2. (a) Results of a three-way interaction analysis for the prediction of alcohol use disorder (AUD) from the family genetic risk score for AUD, the genetic aptitude for educational attainment and the actual level of achieved educational attainment. Each of the three panels is configured like that of Fig. 1a. The top panel then depicts the interaction between the family genetic risk score for AUD and the genetic aptitude for educational attainment in the prediction of AUD in individuals with high levels of achieved educational attainment (+1 s.d.). The middle figure presents the same analyses for those with an average level of achieved educational attainment. The lowest figure presents the same analyses for those very low levels of achieved educational attainment (−2 s.d.). (b) Results of a three-way interaction analysis for the prediction of drug use disorder (DUD) from the family genetic risk score for DUD, the genetic aptitude for educational attainment and the actual level of achieved educational attainment. Each of the three panels is configured like that of Fig. 1b. The top panel then depicts the interaction between the family genetic risk score for DUD and the genetic aptitude for educational attainment in the prediction of DUD in individuals with high levels of achieved educational attainment (+1 s.d.). The middle figure presents the same analyses for those with an average level of achieved educational attainment. The lowest figure presents the same analyses for those very low levels of achieved educational attainment (−2 s.d.).
Discussion
The goal of this paper was to determine whether the GAEA moderates the impact of the genetic risks for AUD and DUD on rates of, respectively, AUD and DUD. We addressed this question with three analyses. We review them in turn.
First, in a national sample of individuals affected with AUD and DUD, we found, consistent with prior studies of the phenotype of educational attainment (Davis & Slutske, Reference Davis and Slutske2018; Hamdi et al., Reference Hamdi, Krueger and South2015), that high GAEA attenuated the impact of genetic liability on the risk of both AUD and DUD. This effect was clearly seen in the slope of the FGRS curves depicted in Fig. 1. For both AUD and DUD, the slope of the curve was shallowest for those with the highest GAEA and became progressively steeper as the GAEA became lower.
Second, while the Swedish registry does not contain the fine-grained phenotypic assessments that would permit us to test the specific mechanisms through which the GAEA protected against the genetic risk for AUD and DUD, we could evaluate one important hypothesis – the degree to which the protective effect of GAEA results from its impact on the actual level of educational attainment. Contrary to our expectation, our analyses demonstrated that almost none of the moderating effects of GAEA on genetic risk for AUD and DUD was mediated through the completed years of education. Instead, the resilience to genetic risk for AUD and DUD associated with a high GAEA likely results from the direct effects of GAEA.
We would speculate that the direct protective effect of GAEA on genetic risk for AUD and DUD likely involves its impact on cognition, personality, and development. A number of studies have used twin and molecular-genetic strategies to disentangle the components of the GAEA. Consistent with prior studies (Belsky et al., Reference Belsky, Moffitt, Corcoran, Domingue, Harrington, Hogan and Williams2016; Johnson, Deary, & Iacono, Reference Johnson, Deary and Iacono2009; Tambs et al., Reference Tambs, Sundet, Magnus and Berg1989), Krapohl et al. (Reference Krapohl, Rimfeld, Shakeshaft, Trzaskowski, McMillan, Pingault and Dale2014) found that the strongest component of the heritability of EA is intelligence. That high levels of general cognitive functioning could mediate the protective effect of GAEA in our Swedish sample is supported by a range of studies showing that high IQ in Sweden is related to lower levels of problematic drinking and/or the negative consequences thereof (Müller et al., Reference Müller, Kowalewski, Metzler, Stettbacher, Rössler and Vetter2013; Sjolund, Allebeck, & Hemmingsson, Reference Sjolund, Allebeck and Hemmingsson2012; Sjölund et al., Reference Sjölund, Hemmingsson and Allebeck2015a; Zettergren & Bergman, Reference Zettergren and Bergman2014). Other more specific cognitive processes could also be involved in protective effects such as working memory (Ellingson, Fleming, Vergés, Bartholow, & Sher, Reference Ellingson, Fleming, Vergés, Bartholow and Sher2014).
However, Krapohl et al. (Reference Krapohl, Rimfeld, Shakeshaft, Trzaskowski, McMillan, Pingault and Dale2014) also found that substantial proportions of the genetic propensity to high educational attainment resulted from the personality trait of self-efficacy. Other studies have shown that a polygenic risk score (PRS) for educational attainment predicts high levels of conscientiousness, agreeableness (Smith-Woolley, Selzam, & Plomin, Reference Smith-Woolley, Selzam and Plomin2019), and self-control (Belsky et al., Reference Belsky, Moffitt, Corcoran, Domingue, Harrington, Hogan and Williams2016), personality traits which, along with the correlated constructs of self-efficacy and constraint, predict low levels of impulsivity, risk taking, and substance misuse (Belcher, Volkow, Moeller, & Ferre, Reference Belcher, Volkow, Moeller and Ferre2014; Bogg & Roberts, Reference Bogg and Roberts2004; Kendler et al., Reference Kendler, Neale, Sullivan, Corey, Gardner and Prescott1999; Settles et al., Reference Settles, Fischer, Cyders, Combs, Gunn and Smith2012; Tang, Posner, Rothbart, & Volkow, Reference Tang, Posner, Rothbart and Volkow2015).
Finally, strong academic performance in adolescence, predicted by GAEA, is associated with reduced risks for psychoactive substance use and subsequent AUD and DUD (Fothergill et al., Reference Fothergill, Ensminger, Green, Crum, Robertson and Juon2008; Gauffin, Vinnerljung, Fridell, Hesse, & Hjern, Reference Gauffin, Vinnerljung, Fridell, Hesse and Hjern2013; Schulenberg, Bachman, O'Malley, & Johnston, Reference Schulenberg, Bachman, O'Malley and Johnston1994). Evidence that this relationship is likely causal comes from individual and school-based interventions for children and adolescents showing that improving academic performance results in lower psychoactive substance use (Eggert, Thompson, Herting, Nicholas, & Dicker, Reference Eggert, Thompson, Herting, Nicholas and Dicker1994; Fletcher, Bonell, Sorhaindo, & Strange, Reference Fletcher, Bonell, Sorhaindo and Strange2009; Lewis et al., Reference Lewis, Bavarian, Snyder, Acock, Day, DuBois and Flay2012). Furthermore, recent prospective cohort studies employing instrumental variable and co-relative designs have shown that high academic performance was likely causally related to a lowered risk for both DUD (Kendler et al., Reference Kendler, Ohlsson, Fagan, Lichtenstein, Sundquist and Sundquist2018b) and AUD (Kendler et al., Reference Kendler, Ohlsson, Fagan, Lichtenstein, Sundquist and Sundquist2020). These findings are consistent with the predictions of social control (Hirschi, Reference Hirschi1969) and social development theories (Kellam et al., Reference Kellam, Brown, Poduska, Ialongo, Wang, Toyinbo and Wilcox2008), which posit that students who succeed academically develop positive attachments to school. These attachments tend to facilitate a commitment to a prosocial lifestyle which in turn reduces the risk for AUD or DUD. By contrast, those who perform poorly in school are more likely to adopt antisocial lifestyles including an elevated risk for AUD and DUD. Thus, high GAEA would predict a developmental pathway marked by academic achievement which leads to the adoption of prosocial attitudes which could protect again the impact of high genetic risk for DUD or AUD.
Our third analysis examined whether the magnitude of the interactions between GAEA and FGRSAUD and FGRSDUD depended on the level of educational attainment. We found a robust three-way interaction which revealed a much stronger moderation of the genetic risk for AUD and DUD by GAEA in individuals with low v. high educational attainment. That is, much of the buffering by GAEA of the impact of elevated genetic risk for AUD and DUD in the Swedish population occurs in those with low levels of education. If targeted interventions could be developed on the basis of our findings with the goal of reducing future risk for substance use disorders, such interventions would have the highest impact if directed to those with high familial risk for AUD and/or DUD and low educational attainment.
Limitations
These findings should be viewed in the context of six potential methodological limitations. First, our results are dependent on the quality of the diagnosis of AUD and DUD in Swedish registries. While such registry data have important advantages (e.g. no refusals or reporting biases), it will not identify all the same cases as an interview-based assessment. Our affected subjects are likely, on average, more severely ill than those meeting DSM-5 criteria (American Psychiatric Association, 2013) at interview, although our lifetime prevalence of AUD and DUD is only moderately lower than those identified in nearby Norway (Kringlen, Torgersen, & Cramer, Reference Kringlen, Torgersen and Cramer2001). The validity of our definitions is supported by the high rates of concordance observed across our ascertainment methods (Kendler, Lönn, Salvatore, Sundquist, & Sundquist, Reference Kendler, Lönn, Salvatore, Sundquist and Sundquist2018a; Kendler et al., Reference Kendler, Ji, Edwards, Ohlsson, Sundquist and Sundquist2015a, Reference Kendler, Ohlsson, Maes, Sundquist, Lichtenstein and Sundquist2015b). Furthermore, the pattern of resemblance in relatives for AUD and DUD seen in Sweden is similar to those found in other samples based on personal interviews (Prescott & Kendler, Reference Prescott and Kendler1999; Tsuang et al., Reference Tsuang, Lyons, Eisen, Goldberg, True, Lin and Eaves1996).
Second, our FGRS, used previously in this journal (Kendler et al., Reference Kendler, Ohlsson, Mościcki, Sundquist, Edwards and SundquistIn press), assesses genetic risk on the basis of the aggregation of disorders in close and distant biological relatives and is not equivalent to a molecular PRS. Our corrections for shared environmental effects with parents and siblings are approximate. We also correct for year of birth to control for possible cohort effects and county of birth control for regional variations in diagnostic or policing practices. Appendix Table 3 shows that our final genetic risk scores are not highly sensitive to key steps in their calculation, as their deletion produces results that correlate highly with those from the full model and have similar predictive power. We show (Appendix Fig. 1) that the FGRS for DUD and AUD is relatively stable across regions within Sweden and have found, when controlling for pedigree structure, that the FGRS correlates very highly with a recently published quantitative family-history score using a different analytic approach (Hujoel, Gazal, Loh, Patterson, & Price, Reference Hujoel, Gazal, Loh, Patterson and Price2020).
We also explored, in several different analyses, the impact of the degree of relatives included in the FGRS calculation. In Appendix Table 4, we show the high correlations in the FGRS reported in the MS for 1st through 5th degree relatives with those obtained by eliminating more distant relatives and in Appendix Table 5, show the polychoric correlation between AUD/DUD and FGRSAUD, FGRSDUD, and FGRSEDU as a function of the breadth of relatives considered. The results are relatively stable until we trim down to only 1st and 2nd degree relatives. In Table 6, we recalculate our main results with various trimmings of the classes of relatives considered which shows that our findings are not sensitive to the particular set of relatives that contribute to the FGRS.
Third, as in other samples (Compton, Thomas, Stinson, & Grant, Reference Compton, Thomas, Stinson and Grant2007; Kessler et al., Reference Kessler, Crum, Warner, Nelson, Schulenberg and Anthony1997), we observe substantial comorbidity between AUD and DUD. We therefore re-ran our key results in four groups: AUD with DUD only cases censored, DUD with AUD cases only censored, AUD with comorbid DUD-AUD cases censored, and DUD with comorbid DUD-AUD cases censored. As seen in Appendix Table 7 and Fig. 2, the broad pattern of moderation seen in our main analyses was evident in all four of these groups.
Fourth, would our findings differ if analyzed with a more standard multiplicative Cox model, rather than the additive model employed? We present the results for AUD and DUD using such a model in Appendix Table 8 and Fig. 3a and b. The same pattern of results is observed.
Fifth, ADHD diagnoses are genetically correlated with risk for substance use disorders (Hicks, Iacono, & McGue, Reference Hicks, Iacono and McGue2012; Kendler & Myers, Reference Kendler and Myers2014) and might contribute to their association with GAEA, but were not considered here.
Finally, there is substantial current interest in the relationship of the genetic influences on AUD and alcohol consumption (e.g. Kranzler et al., Reference Kranzler, Zhou, Kember, Smith, Justice, Damrauer and Reid2019). While it would therefore be of interest to see if we obtained similar or different results for the interactions of GAEA on genetic control of alcohol consumption as we do for FGRSAUD, no registry data are available in Sweden for individual alcohol consumption levels.
Conclusions
The genetic propensity to high educational attainment substantially moderates genetic risk for both AUD and DUD. The impact of the genetic liability to AUD and DUD on the risk of illness is substantially attenuated in those with high v. low GAEA. This effect is not appreciably mediated by the level of educational attainment. Rather, it likely arises because individuals with high GAEA tend to have higher levels of intelligence, conscientiousness, self-control, and follow a developmental pathway that includes good school performance leading to prosocial attitudes. All these factors, in aggregate, seem to protect even those at high genetic risk, from developing AUD and DUD. The protective effects of a genetic propensity to high educational attainment on risk for substance use disorders are particularly strong in those with low educational attainment. These naturalistic findings could form the basis of prevention efforts in high-risk youth. Our findings are also of theoretical interest and encourage further investigations, including those using polygenic risks scores, that examine the degree to which genetic liabilities for psychiatric and substance use disorders moderate the impact on one another on disorder risk.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291721005134.
Financial support
Grant Support: This project was supported by grants R01AA023534 and R01DA030005 from the National Institutes of Health and from the Swedish Research Council (2016-01176), as well as Avtal om Läkarutbildning och Forskning (ALF) funding from Region Skåne. The funding source had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Conflict of interest
None.
Access to data and data analysis
Kristina Sundquist had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
