Hostname: page-component-7b9c58cd5d-f9bf7 Total loading time: 0 Render date: 2025-03-15T14:43:49.642Z Has data issue: false hasContentIssue false

Assessing the interplay between multigenic and environmental influences on adolescent to adult pathways of antisocial behaviors

Published online by Cambridge University Press:  22 November 2017

James J. Li*
Affiliation:
University of Wisconsin–Madison
*
Address correspondence and reprint requests to: James J. Li, Department of Psychology, University of Wisconsin–Madison, 1202 West Johnson Street, Madison, WI 53706; E-mail: james.li@wisc.edu.
Rights & Permissions [Opens in a new window]

Abstract

The current investigation utilized a developmental psychopathology approach to test the hypothesis that multigenic (i.e., dopaminergic and serotonergic genes) and multienvironmental factors interactively contribute to developmental pathways of antisocial behavior (ASB). A sample of 8,834 Caucasian individuals from the National Longitudinal Study of Adolescent to Adult Health (Add Health) were used to (a) examine the developmental pathways of ASB from age 13 to 32 using growth mixture modeling, (b) compute weighted multigenic risk scores (Add Health MRS) for ASB from six well-characterized polymorphisms in dopamine and serotonin genes, and (c) test the interaction between the Add Health MRS and a measures of support (incorporating indicators of both positive and negative support from parents and schools). Four pathways of adolescent to adult ASB emerged from the growth mixture models: low, adolescence-peaked, high decline, and persistent. Add Health MRS predicted the persistent ASB pathway, but not other ASB pathways. Males with high Add Health MRS, but not low MRS, had significantly greater odds of being in the adolescence-peaked pathway relative to the low pathway at low levels of school connectedness. Nonfamilial environmental influences during adolescence may have a cumulative impact on the development of ASB, particularly among males with greater underlying genetic risks.

Type
Special Issue Articles
Copyright
Copyright © Cambridge University Press 2017 

Antisocial behaviors (ASB) consist of behaviors and attitudes that violate societal norms, including both aggressive (e.g., physical fighting and assault with a weapon) and nonaggressive rule breaking (e.g., theft, property destruction, and truancy) dimensions (Burt, Reference Burt2009; Frick et al., Reference Frick, Lahey, Loeber, Tannenbaum, Van Horn, Christ and Hanson1993; Loeber & Schmaling, Reference Loeber and Schmaling1985). Although ASB are relatively common during adolescence (Moffitt, Reference Moffitt1993), they are highly impairing as adolescents engaged in ASB are at increased risk for academic failure (Farrington, Reference Farrington1995), developing psychopathology (Castellani et al., Reference Castellani, Pastorelli, Eisenberg, Gerbino, Di Giunta, Ceravolo and Milioni2014), and are more likely to engage in criminal activity in adulthood (Mercer, Keijers, Crocetti, Branje, & Meeus, Reference Mercer, Keijsers, Crocetti, Branje and Meeus2016). Furthermore, ASB are a tremendous cost to society, including the increased economic burden they convey on our health care system, schools and communities, and the juvenile justice system (Foster & Jones, Reference Foster and Jones2005). The indirect costs of juvenile incarceration alone have been estimated to be between $8 and $21 billion annually at state and local government levels (Justice Policy Institute, 2014). Given the pervasive and substantial impact of ASB, a better understanding of the multitude of developmental pathways that contribute to ASB are paramount.

Adolescent to Adult Pathways of ASB

Adolescence is marked by drastic changes in social, psychological, and brain development (Casey, Jones, & Hare, Reference Casey, Jones and Hare2008). In contrast to their childhood years, adolescents are more peer directed in their social interactions and simultaneously seek greater separation from their caregivers, leading to increased risky behaviors and a higher likelihood of engagement in ASB (Kelley, Schochet, & Landry, Reference Kelley, Schochet and Landry2004). Not all adolescents engage in ASB, and most adolescents who do tend to desist as they age into adulthood (Moffitt, Reference Moffitt1993, Reference Moffitt2001) given that prefrontal and subcortical regions of the brain (which pertain to executive control and emotion reactivity) become more fully developed by adulthood (Casey et al., Reference Casey, Jones and Hare2008). However, a small subset of youth also persist in their ASB well into adulthood (Moffit, Reference Moffitt1993, Reference Moffitt2001). Sampson and Laub (Reference Sampson and Laub2003) examined the rates of criminal offending (i.e., violent crime, property offenses, and alcohol/drug crimes) in a sample of 500 delinquent males who were followed from childhood (youngest was age 7) to late adulthood (oldest was age 70) and found evidence of two distinct pathways where peak criminal offending was exhibited during the mid-30s. Massoglia (Reference Massoglia2006) used latent class models to characterize the structure of ASB during adolescence (ages 11–17) and examined whether individuals transition out of ASB by adulthood (ages 21–27) using data from the National Survey of Youth. Four distinct adolescent ASB classes emerged from the analysis: normative, predatory (i.e., aggressive and violent), drug use, and pervasive (i.e., chronic offenders), and a significant proportion of adolescents who engaged in predatory, drug use, or pervasive ASB transitioned between these classes as they aged into adulthood, rather than having completely desisted in their ASB. Several studies have similarly identified unique groups of individuals who peak in ASB during adolescence but persist in modest levels of ASB well into their mid-20s (DiLalla & Gottesman, Reference DiLalla and Gottesman1989; Moffitt, Reference Moffitt, Lahey, Moffitt and Caspi2003; Nagin & Tremblay, Reference Nagin and Tremblay2005; Odgers et al., Reference Odgers, Moffitt, Broadbent, Dickson, Hancox, Harrington and Caspi2008), including studies using data from the National Longitudinal Study of Adolescent to Adult Health (Add Health; Barnes, Beaver, & Boutwell, Reference Barnes, Beaver and Boutwell2011; Tung & Lee, Reference Tung and Lee2016; Zheng & Cleveland, Reference Zheng and Cleveland2015). Taken together, these studies suggest that there are likely to be unique pathways of ASB from adolescence into adulthood. It has been postulated that individual differences in biological and environmental antecedents may be associated with distinct developmental pathways of ASB (Burt, Reference Burt2009; Farrington, Reference Farrington1987; Ingoldsby & Shaw, Reference Ingoldsby and Shaw2002; Taylor, Iacono, & McGue, Reference Taylor, Iacono and McGue2000).

Genetic Predispositions for ASB

Behavioral genetic studies suggest that the developmental trajectories of ASB differ in their heritability. For example, based on a twin subsample (N = 2,284 twin pairs) from the Add Health, genetic influences accounted for 56% to 70% of the variance within the “life-course persistent offender” pathway (i.e., individuals who positively endorsed delinquency items across all three waves of the Add Health study) and 35% of the variance in the “adolescence-limited offender” pathway (i.e., individuals who positively endorsed delinquency items at Wave 3, but not at Waves 1 or 2; Barnes et al., Reference Barnes, Beaver and Boutwell2011). More recently, Zheng and Cleveland (Reference Zheng and Cleveland2015) used latent class growth analysis to characterize homogenous pathways of ASB within a twin sample, also from Add Health (N = 356 twin pairs) and found that additive genetic influences had a more significant influence for males in the “chronic” class than “desister” or “decliner” classes of ASB. Nonshared environmental factors were significantly more relevant for desisters and decliners than either shared environmental or genetic factors (Zheng & Cleveland, Reference Zheng and Cleveland2015). These studies suggest that genetic influences may play a more important role for certain (i.e., chronic or persistent) developmental pathways of ASB, although the precise genetic architecture underlying these pathways has yet to be conclusively determined.

Although identifying the genes for ASB remains an ongoing effort, genes regulating dopamine and serotonin transmission have been a central focus because of the prominent role they play in human and animal aggression (see recent reviews by Ficks & Waldman, Reference Ficks and Waldman2014; and Veroude et al., Reference Veroude, Zhang-James, Fernàndez-Castillo, Bakker, Cormand and Faraone2016). The dopamine system is principally involved in behavioral activation and reward processing (Everitt & Robbins, Reference Everitt and Robbins2000). This system has long been implicated in modulating several ASB-related traits such as impulsivity and physical aggression in mice (Tidey & Miczek, Reference Tidey and Miczek1996) and in humans (Buckholtz et al., Reference Buckholtz, Treadway, Cowan, Woodward, Li, Ansari and Kessler2010; Schlüter et al., Reference Schlüter, Winz, Henkel, Prinz, Rademacher, Schmaljohann and Mottaghy2013). Genetic variation in polymorphisms relevant to dopamine transmission have been linked to various subtypes of ASB as well (Burt & Mikolajewski, Reference Burt and Mikolajewski2008; Windhorst et al., Reference Windhorst, Mileva-Seitz, Rippe, Tiemeier, Jaddoe, Verhulst and Bakermans-Kranenburg2016). Separate studies have also reported on the association of the serotonin system, which is involved in emotion regulation and inhibition, and ASB-related phenotypes in humans (Davidson, Putnam, & Larson, Reference Davidson, Putnam and Larson2000; Moore, Scarpa, & Raine, Reference Moore, Scarpa and Raine2002) and in animal models (e.g., Gibbons, Barr, Bridger, & Liebowitz, Reference Gibbons, Barr, Bridger and Liebowitz1979; Holmes, Murphy, & Crawley, Reference Holmes, Murphy and Crawley2002). Genetic influences on serotonergic systems have also been implicated in human ASB (Haberstick, Smolen, & Hewitt, Reference Haberstick, Smolen and Hewitt2006).

However, most studies have studied dopamine and serotonin as separate systems in relation to ASB, despite increasing evidence that they are interconnected (D'Souza & Craig, Reference D'Souza and Craig2008; Seo, Patrick, & Kennealy, Reference Seo, Patrick and Kennealy2008). For example, immediately following a threat, mice were found to have increased prefrontal dopaminergic activity and a corresponding decrease in serotonergic functioning (Van Erp & Miczek, Reference Van Erp and Miczek2000). Human studies have similarly shown that high levels of dopamine activity corresponded with lower levels of serotonin activity in association with human aggression and psychopathic traits (Soderstrom, Blennow, Manhem, & Forsman, Reference Soderstrom, Blennow, Manhem and Forsman2001). These studies suggest that genetic variation in these systems should be considered from a multigenic perspective, such that dopaminergic and serotonergic gene variants may collectively account for a greater portion of the phenotypic variance than they would independently (Plomin, Haworth, & Davis, Reference Plomin, Haworth and Davis2009). Few studies have examined the multigenic influence of dopamine and serotonin genes on ASB. Hence, the current study focuses on polymorphisms across six well-characterized dopaminergic and serotonergic genes (i.e., dopamine active transporter 1 [DAT1], dopamine receptor D4 [DRD4], dopamine receptor D2 [DRD2], catechol-O-methyltransferase [COMT], monoamine oxidase A [MAOA], and serotonin transporter linked polymorphic regon [5-HTTLPR]) in relation to ASB.

Methods to quantify multigenic influences from known biological systems have been unsatisfactory, however. One common approach is to sum the number of putative risk alleles across several polymorphisms, such that total scores represent an individual's cumulative genetic risk or susceptibility to an outcome (e.g., Beaver & Belsky, Reference Beaver and Belsky2012; Belsky & Beaver, Reference Belsky and Beaver2011; Nikolova, Ferrell, Manuck, & Hariri, Reference Nikolova, Ferrell, Manuck and Hariri2011; Stice, Yokum, Burger, Epstein, & Smolen, Reference Stice, Yokum, Burger, Epstein and Smolen2012; Thibodeau, Cicchetti, & Rogosch, Reference Thibodeau, Cicchetti and Rogosch2015). This method is used in studies where the genetic variants are selected based on prior knowledge of their biological function, rather than through an empirical investigation of the most highly associated genetic variants across the genome (Derringer et al., Reference Derringer, Krueger, Dick, Saccone, Grucza, Agrawal and Nurnberger2010; Marceau et al., Reference Marceau, Palmer, Neiderhiser, Smith, McGeary and Knopik2016; Stice et al., Reference Stice, Yokum, Burger, Epstein and Smolen2012). The primary advantages of this approach are that it reduces the burden of multiple testing and consequently improves power (Windhorst et al., Reference Windhorst, Mileva-Seitz, Rippe, Tiemeier, Jaddoe, Verhulst and Bakermans-Kranenburg2016). However, this approach may not be biologically plausible; first, the effect of each risk allele on the outcome cannot be assumed to be equivalent at each locus, given that genome-wide association studies (GWAS) of ASB have demonstrated that some loci confer larger effects on ASB than others (Pappa et al., Reference Pappa, St. Pourcain, Benke, Cavadino, Hakulinen, Nivard and Evans2016; Salvatore & Dick, Reference Salvatore and Dick2016; Tielbeek et al., Reference Tielbeek, Medland, Benyamin, Byrne, Heath, Madden and Verweij2012). Second, assumptions regarding the risk variant at each locus need to be made to compute this score, which will not likely account for the possibility of inverse or reciprocal effects of dopamine and serotonergic genes (Seo et al., Reference Seo, Patrick and Kennealy2008). Third, findings have been mixed with regard to whether genetic risk scores computed in this way are predictive of ASB and ASB-related traits (Gizer & Waldman, Reference Gizer and Waldman2012; Thibodeau et al., Reference Thibodeau, Cicchetti and Rogosch2015).

An alternative strategy is to derive robust prediction estimates for each locus so that the genetic risk score also reflects the relative effects of each locus on the outcome. Obtaining prediction estimates for each locus usually requires a comparable but independent data set (i.e., a testing set) so that the model is not overfitted to the data (i.e., the extent to which a model for one data set might not be generalizable to a different data set). Alternatively, one can split the data in half, such that one half of the data would be used as the training set and the other half would be used as the validation set to avoid issues with overfitting. However, this approach has the significant disadvantage of reducing the sample size of the resultant models by half, a sacrifice that genetically informed studies can ill afford. Furthermore, the results would depend entirely on the random sets that were used for the training and validation sample, which may lead to unstable estimates if the process were repeated. A variation of this approach that does not require an independent data set per se is a technique called k-fold cross-validation (Kohavi, Reference Kohavi1995). In this approach, a single data set is partitioned k number of times to create k unique and equally sized validation samples. The remaining k – 1 partitions are used to measure or validate how well the model performs in terms of its predictive ability. To attain stability of the model, this process is repeated k times (i.e., folds) on the different partitions, and the prediction estimates are then averaged over each fold (Cordell, Reference Cordell2009). This method has been previously used in polygenic risk score predictions for cardiovascular disease and adult psychopathology (e.g., International Schizophrenia Consortium, 2009; Simonson, Wills, Keller, & McQueen, Reference Simonson, Wills, Keller and McQueen2011). To address the conceptual limitations of a total risk alleles approach, the current study used k-fold cross-validation to derive robust prediction estimates for ASB from multiple dopaminergic and serotonergic gene variants, which were then be used to derive weighted multigenic risk scores (MRS).

Gene × Environment (G × E) Interaction Supportive Parenting and School Connectedness

The expression of ASB may depend not only on an individual's underlying genetic susceptibility for ASB but also on whether the individual is exposed to certain types of environments (G × E; Caspi et al., Reference Caspi, McClay, Moffitt, Mill, Martin, Craig and Poulton2002; Jaffee et al., Reference Jaffee, Caspi, Moffitt, Dodge, Rutter, Taylor and Tully2005). Home environments characterized by high instability and low parental involvement and warmth are well-established risk factors for ASB (Catalano & Hawkins, Reference Catalano, Hawkins and Hawkins1996; Doyle & Markiewicz, Reference Doyle and Markiewicz2005), and many G × E studies of ASB have focused on the moderating role of family-related adversity. Meta-analyses suggest that polymorphic variation in dopaminergic and serotonergic genes increase risk for ASB in the presence, but not necessarily in the absence, of severe family adversity (e.g., maltreatment; Byrd & Manuck, Reference Byrd and Manuck2014; Kim-Cohen et al., Reference Kim-Cohen, Caspi, Taylor, Williams, Newcombe, Craig and Moffitt2006; Tielbeek et al., Reference Tielbeek, Karlsson Linnér, Beers, Posthuma, Popma and Polderman2016). However, a myopic focus on the extreme (i.e., negative) end of the environment spectrum betrays the significant amount of variation that exists within an individual's environment, including his or her exposure to more supportive influences from family. Emerging evidence suggests higher frequency warm and supportive maternal behavior at age 12 leads to promotive changes in the amygdala and prefrontal cortex in the offspring 4 years later, suggesting that positive parenting may impact brain development even during the period of adolescence (Whittle et al., Reference Whittle, Simmons, Dennison, Vijayakumar, Schwartz, Yap and Allen2014).

One promising theory posits that genetic influences may confer heightened sensitivity to both adverse and enriched environments, a phenomenon known as differential susceptibility (Belsky & Pluess, Reference Belsky and Pluess2009). Although tests of differential susceptibility are still emerging, there is promising evidence from studies employing rigorous, developmentally sensitive methods to study these effects in the context of parenting (Bakermans-Kranenburg & van IJzendoorn, Reference Bakermans-Kranenburg and van Ijzendoorn2011). For instance, Tung and Lee (Reference Tung and Lee2016) found that perceived parental support (i.e., closeness and warmth) moderated the association between adolescent 5-HTTLPR genotype and their trajectories of ASB measured into early adulthood (mid-20s). Specifically, individuals with the long–long (LL) genotype were least likely to be in the “adolescent peak” ASB trajectory class if they reported high-quality parental support, but were the most likely to be in this trajectory of ASB if they reported low-quality parental support compared to individuals who carried at least one short allele (Tung & Lee, Reference Tung and Lee2016). Belsky and Beaver (Reference Belsky and Beaver2011) found that adolescent youths carrying a greater number of sensitivity alleles (using a simple composite of alleles) in DAT1, DRD2, DRD4, MAOA, and 5-HTTLPR had simultaneously higher and lower scores on self-regulation (e.g., impulsivity) when they self-reported low- and high-quality supportive parenting, respectively. Thus, genes may potentially contribute to an individual's sensitivity to variations in the family environment, leading to developmental outcomes that are “for better and for worse” (Belsky & Pluess, Reference Belsky and Pluess2009).

However, meta-analytic evidence also show that familial influences tend to decrease while nonfamilial influences increase with age as it relates to the development of ASB (Ferguson, Reference Ferguson2010; Rhee & Waldman, Reference Rhee and Waldman2002), suggesting that nonfamilial factors may also play a significant role in moderating genetic susceptibility to ASB during the period of adolescence and early adulthood. The school environment represents a potentially important ecological context that may plausibly interact with genetic factors underlying adolescent to adult trajectories of ASB. The extent to which adolescents felt accepted and supported by others in their school environment predicted lower levels of delinquency and alcohol use relative to those who did not feel as connected with their school environment (Crosnoe, Erickson, & Dornbusch, Reference Crosnoe, Erickson and Dornbusch2002). Several studies have also found that adolescents who felt disengaged from their schools or had poor relationships with their teachers were more likely to engage in delinquency during adolescence and had more substance use and mental health problems in their later years (Bond et al., Reference Bond, Butler, Thomas, Carlin, Glover, Bowes and Patton2007; Jacobson & Rowe, Reference Jacobson and Rowe1999; Resnick et al., Reference Resnick, Bearman, Blum, Bauman, Harris, Jones and Ireland1997; Schochet, Dadds, Ham, & Montague, 2006). Research from affective neuroscience suggests that social connectedness more broadly may impact patterns of brain development. For instance, social isolation (i.e., feeling disconnected or rejected) predicted greater activity in the anterior cingulate cortex among adolescents, a region of the brain that is involved in emotion regulation (Masten et al., Reference Masten, Eisenberger, Borofsky, Pfeifer, McNealy, Mazziotta and Dapretto2009). In a study of adults, individuals who reported being more socially connected with others at baseline had increases in vagal tone over time, an index of autonomic flexibility that is linked to well-being (Kok & Fredrickson, Reference Kok and Fredrickson2010). Social connectedness (and specifically, school connectedness) may also impact gene expression. Twin studies have shown that monozygotic twins may have different perceptions of their school experience (such as involvement in class, relationship with staff and teachers, and peer-related stressors), leading to individual differences in their behavior and academic achievement (Asbury, Almeida, Hibel, Harlaar, & Plomin, Reference Asbury, Almeida, Hibel, Harlaar and Plomin2008; Walker & Plomin, Reference Walker and Plomin2006). Finally, animal models have demonstrated the promotive effects of environmental enrichment more broadly (i.e., a combination of social relationships, cognitive, sensory, and physical resources) on behavioral and neuronal development, even after early exposure to stress (Hackman, Farah, & Meaney, Reference Hackman, Farah and Meaney2010). Despite its biological plausibility, studies of G × E have largely ignored the role that the school environment may play in the development of ASB.

It is also noteworthy that previous G × E studies focusing on nonfamilial environmental influences on ASB have largely examined the role of deviant peers (Cleveland, Wiebe, & Rowe, Reference Cleveland, Wiebe and Rowe2005; Kendler et al., Reference Kendler, Jacobson, Gardner, Gillespie, Aggen and Prescott2007; Lee, Reference Lee2011). Few studies have investigated the potential moderating role of the school environment independent from deviant peer effects. This is a crucial gap in the literature because adolescents are also exposed to varying degrees of structure and support from their teachers, safety from staff and administrators, and access to school-based programs and activities that may buffer the individual from negative peer or familial influences as they relate to the development of ASB (Bond et al., Reference Bond, Butler, Thomas, Carlin, Glover, Bowes and Patton2007; Mrug & Windle, Reference Mrug and Windle2009). One exception is the study by Crosnoe et al. (Reference Crosnoe, Erickson and Dornbusch2002), which showed that high levels of school connectedness attenuated the negative effects of deviant peer affiliation on adolescent ASB, suggesting that support from teachers and other positive characteristics of the school may buffer the negative effects of deviant peers on ASB. Hence, the degree to which adolescents feel supported, accepted, and connected to their schools, independent from the influence of their peers, may serve as not only a potential risk factor for future negative outcomes but also a source of support and enrichment that leads to positive behavioral adjustment, making this a strong candidate for assessing differential susceptibility.

In addition, some studies have shown that adolescents with a genetic propensity to be antisocial may select (or attract) peers who are also antisocial (Kendler, Jacobson, Myers, & Eaves, Reference Kendler, Jacobson, Myers and Eaves2008), which is known as a gene–environment correlation (rGE). Hence, there is some possibility that G × E studies focusing on peer closeness (or any other plausible environmental exposure) may have been driven by an rGE effect rather than a true G × E (Jaffee & Price, Reference Jaffee and Price2007). As relatively few G × E studies have systematically tested for concurrent rGE effects (Knafo & Jaffee, Reference Knafo and Jaffee2013), the possibility that environmental influences, including school connectedness and supportive parenting, and ASB may be partially explained by an individual's genetic constitution is also an important inquiry in studies of G × E.

Current Study

The current study uses a rigorous developmental psychopathology approach to examine the interplay between multigenic and multi-environmental influences on the pathways of ASB. The aims and hypotheses of the study were as follows:

  1. 1. To examine pathways of ASB from adolescence to adulthood. Using data from a large, nationally representative prospective sample of adolescents (i.e., Add Health, Waves 1–4) who have been followed into their early 30s, it is hypothesized that later onset developmental pathways of ASB will be identified in line with previous studies that have similarly examined ASB into the later years of adulthood.

  2. 2. To assess the relative influence of each dopamine and serotonin gene variant (i.e., DAT1, DRD4, DRD2, MAOA, COMT, and 5-HTTLPR) on ASB and compute weighted MRS on the basis of each polymorphism's relative contribution to ASB. A multigenic approach is hypothesized to more reliably predict the developmental pathways of ASB than each single polymorphism would alone, and possibly provide a stronger signal to detect G × E.

  3. 3. To test developmental G × E effects from the differential susceptibility perspective, whereby positive and negative dimensions of the adolescents’ family and school environment were examined as potential moderators of genetic susceptibility in relation to the developmental pathways of ASB. Given the increased importance of nonfamilial influences on behavioral outcomes during adolescence (Ferguson, Reference Ferguson2010; Rhee & Waldman, Reference Rhee and Waldman2002), it is hypothesized that variation in how connected an adolescent feels in his/her school will moderate the association between genetic susceptibility for ASB and developmental pathways of ASB, such that individuals with greater genetic susceptibility for ASB will be more likely to engage in ASB under conditions of low school connectedness, but the same individuals will also be protected from ASB under conditions of high school connectedness.

Method

Participants

Add Health began in 1994 and ascertained a stratified random sample of adolescents in Grades 7–12 from high schools across the United States. Data were collected from adolescents, parents, fellow students, school administrators, siblings, friends, and romantic partners across four waves of data collection: Wave 1 (1994–1995, Grades 7–12, N = 20,745), Wave 2 (1995–1996, Grades 8–12, N = 14,738), Wave 3 (2001–2002, ages 18–26, N = 15,197), and Wave 4 (2007–2008, ages 24–32, N = 15,701). Data collection for Wave 5 is currently ongoing. At Wave 1, the average age was 15.22 (SD = 1.65, range = 12–20), the racial/ethnic composition was 62.1% Caucasian (including Hispanic or Latino), 23.0% Black or African American, 7.1% Asian or Pacific Islander, 1.2% Native American, and 6.6% other, and 49.5% were males. Due to the nonrandom association of genotypes and race/ethnicityFootnote 1 (i.e., population stratification; Hutchison, Stallings, McGeary, & Bryan, Reference Hutchison, Stallings, McGeary and Bryan2004), analyses for the current study were performed for the Caucasian subsample of individuals in which genetic data were available, resulting in a sample of 8,834 with full genotypic and phenotypic information germane to the current investigation (individuals with mixed ancestry were excluded from the sample). Genotypes in the current study were nonrandomly associated with race/ethnicity. Full details of the study design and available measures can be obtained at http://www.cpc.unc.edu/projects/addhealth/design.

Measures

ASB

ASB was assessed during in-home interviews (“delinquency scale” and “fighting and violence”) conducted at Waves 1–4. Ten identical (or highly similar) items were assessed at each wave, reflecting the presence of nonaggressive rule-breaking behaviors (e.g., property damage, stealing, and selling drugs) and aggressive behaviors (e.g., pulling a knife or gun on someone and engaging in a fight). Items were summed to create a composite score (range = 0–10) at each wave.Footnote 2 The scale demonstrated good to adequate internal consistency across waves (Cronbach αs = 0.75, 0.75, 0.61, and 0.65 for Waves 1, 2, 3, and 4, respectively).

Supportive parenting

Supportive parenting was assessed during Wave 1 using 12 items from the “protective factors,” “relations with parents,” and the “personality and family” in-home interviews. The items assessed maternal warmth (e.g., “Your mother is warm and loving toward you”), care (e.g., “How much do you feel that your parents care about you?”), closeness (e.g., “How close do you feel to your mother?”), communication quality (e.g., “You are satisfied with the way your mother and you communicate with each other”), understanding (e.g., “How much do you feel that people in your family understand you?”), and the overall quality of the parental relationship (e.g., “Overall, you are satisfied with your relationship with you mother”). This construct has been used by several previous Add Health studies, demonstrating good internal consistency (Cronbach α = 0.83) and predictive validity (Borowsky, Ireland, & Resnik, Reference Borowsky, Ireland and Resnick2001; Li, Berk, & Lee, Reference Li, Berk and Lee2013; Resnick et al., Reference Resnick, Bearman, Blum, Bauman, Harris, Jones and Ireland1997). All 12 items were rated on a 5-point scale, where 0 = not at all and 4 = very much. Scores for each item were summed, and the resultant score was then standardized.

School connectedness

Six items from Wave 1 assessed the degree to which participants felt connected to their school. These items came from the “relations with parents” and “academics and education” in-home interviews. The items assessed feelings of belongingness (e.g., “You feel close to people at your school,” “You feel like you are part of your school,” and “You feel close to people at school”), teacher support (e.g., “Your teachers care a lot about you” and “Teachers treat students fairly”), and safety (e.g., “You feel safe in your school”). Given that the focus of the investigation was on the risk/protective role of the broader school context independent from peer effects, items pertaining to peer closeness were not included. All items were rated on a 5-point scale, where 0 = strongly disagree and 4 = strongly agree. Some items were reverse coded to maintain the consistency of the scale so that higher scores represent higher levels of school connectedness. The school connectedness scale has demonstrated good predictive validity in previous Add Health studies (Batanova & Loukas, Reference Batanova and Loukas2014; McNeely, Nonnemaker, & Blum, Reference McNeely, Nonnemaker and Blum2002). The scale also demonstrated adequate internal consistency (Cronbach α = 0.73). Scores for each item were summed, and the resultant score was then standardized.

Covariates

Peer closeness, substance use behaviors, and household income were included as covariates. Peer closeness and substance use behaviors are known to covary with the variables of interest (i.e., supportive parenting, school connectedness, and ASB) and were thus added into the multinomial logistic regression models to avoid the possibility of confounding due to peer and substance use effects. Peer closeness was assessed during the Wave 1 in-home interview. Items came from the “friends and family” section of the interview. Adolescents were asked to nominate at least one male and one female friend and then answer five yes or no questions pertaining to how frequently they saw their friends (e.g., “Did you go to your friend's house during the past 7 days?” and “Did you meet with your friend after school to hang out or go somewhere during the past 7 days?”) and how close they felt to the friends they nominated (e.g., “Did you talk to your friend about a problem during the past 7 days?” and “Did you talk to your friend the past 7 days?”). The final variable was the composite score of the items. Substance use behaviors included cigarette smoking (“Have you ever smoked cigarettes regularly, that is, at least 1 cigarette every day for 30 days?”), heavy alcohol consumption (“Over the past 12 months, on how many days did you drink 5 or more drinks in a row”), and marijuana use (“During the past 30 days, how many times did you use marijuana?”) assessed during in-home interviews conducted at each wave. Each substance use item was dichotomized and summed to create a composite score, which was then averaged across all four waves. Substance use items were dichotomized in the following ways: cigarette smoking (i.e., having smoked one cigarette a day, every day, for past 30 days), heavy alcohol consumption (i.e., at least one instance in which five or more drinks were consumed in a row over the past year), and marijuana use (i.e., at least one instance in which marijuana was consumed in the past 30 days). Household income was coded on an ordinal scale, where 1 = $0–$20,000 and 6 ≥ $100,000.

Genotyping

DNA samples were collected from participant saliva at Wave 4. Full details of the DNA extraction and genotyping procedures are available in Add Health documentation (see Smolen et al., Reference Smolen, Whitsel, Tabor, Killeya-Jones, Cuthbertson, Hussey and Harris2013). The current investigation focused on six dopaminergic and serotonergic genes associated with ASB, including three variable number tandem repeats (VNTRs; DRD4, DAT1, and MAOA) and three single nucleotide polymorphisms (SNPs; rs1800497 [DRD2 TaqIA], rs4680 [val158met COMT], and rs25531 [triallelic 5-HTTLPR]).

All VNTRs were genotyped using primer sequences that were written 5′ to 3′. DRD4 is located on 11p15.5 and contains a 48 base pair (bp) polymorphism in the third exon, resulting in 2–11 repeats (R). The most common variants are those that consist of 2R, 4R, and 7R; larger variants (i.e., 7R) are associated with a reduction in dopamine binding in vitro compared to smaller variants (Van Tol et al., Reference Van Tol, Wu, Guan, Ohara, Bunzow, Civelli and Jovanovic1992). The 7R allele has most frequently been linked to human aggressive behaviors (Buchmann et al., Reference Buchmann, Zohsel, Blomeyer, Hohm, Hohmann, Jennen-Steinmetz and Esser2014; Hohmann et al., Reference Hohmann, Becker, Fellinger, Banaschewski, Schmidt, Esser and Laucht2009; Nobile et al., Reference Nobile, Giorda, Marino, Carlet, Pastore, Vanzin and Battaglia2007). DAT1 (SLC6A3), located on 5p15.3, contains a 40-bp VNTR in the 3′-untranslated region of the gene. The most common variants consist of 9R and 10R, and the 10R variant is associated with increased expression of the dopamine transporter protein (resulting in decreased striatal dopamine levels; VanNess, Owens, & Kilts, Reference VanNess, Owens and Kilts2005) relative to the 9R variant. The 10R allele and has been linked with aggression and ASB in both adolescent and adult samples (Beaver, Wright, & Walsh, Reference Beaver, Wright and Walsh2008; Burt & Mikolajewski, Reference Burt and Mikolajewski2008; Guo, Roettger, & Shih, Reference Guo, Roettger and Shih2007). MAOA is on Xp11.3-11.4 and contains a 30-bp VNTR in the 5′-regulatory region of the gene. Shorter repeats (2R and 3R) are commonly grouped as “low” activity variants and longer repeats (3.5R, 4R, and 5R) are grouped as “high” activity variants (Sabol, Hu, & Hamer, Reference Sabol, Hu and Hamer1998). The low-activity variant is associated with reduced transcription of an enzyme that metabolizes monoamines relative to the high-activity variant and has been linked to increased risk for aggression and ASB (Buckholtz & Meyer-Lindenberg, Reference Buckholtz and Meyer-Lindenberg2008; Ficks & Waldman, Reference Ficks and Waldman2014). The low-activity variant of MAOA is also an important marker of risk for ASB among children exposed to maltreatment (Kim-Cohen et al., Reference Kim-Cohen, Caspi, Taylor, Williams, Newcombe, Craig and Moffitt2006; Widom & Brzustowicz, Reference Widom and Brzustowicz2006).

DRD2 is located on 11q2.2–2.3 and contains the TaqIA polymorphism (rs1800497). A1 allele carriers have been found to have lower density of the D2 receptor in the striatum and nucleus accumbens compared to individuals homozygous for the A2 allele (Jönsoon et al., Reference Jönsson, Nöthen, Grünhage, Farde, Nakashima, Propping and Sedvall1999; Noble, Gottschalk, Fallon, Ritchie, & Wu, Reference Noble, Gottschalk, Fallon, Ritchie and Wu1997; Pohjalainen et al., Reference Pohjalainen, Rinne, Någren, Lehikoinen, Anttila, Syvälahti and Hietala1998), resulting in lower DA signaling. Individuals with the A1 allele have been found to have a higher risk of developing ASB compared to individuals without this allele (Guo, Roettger, & Shih, Reference Guo, Roettger and Shih2007). 5-HTTLPR (solute carrier family C6, member 4 [SLC6A4]) maps to 17q11.1–17q12 and contains a 44-bp polymorphism in the 5′-regulatory region (Heils et al., Reference Heils, Teufel, Petri, Stöber, Riederer, Bengel and Lesch1996). The most common alleles contain 14 (short) and 16 (long) repeat units. The rs25531 SNP was also assayed, allowing for the determination of long A (LA) and long G (LG) alleles. The LG allele transcriptionally similar to the short allele of SLC6A4 and is frequently combined with SLC6A4 (Hu et al., Reference Hu, Lipsky, Zhu, Akhtar, Taubman, Greenberg and Murphy2006). The short allele is associated with decreased transcriptional efficiency of the serotonin transporter proteins (relative to the long allele) and is associated with ASB (Haberstick et al., Reference Haberstick, Smolen and Hewitt2006; Liao, Hong, Shih, & Tsai, Reference Liao, Hong, Shih and Tsai2004), although null associations have also been reported (Sakai et al., Reference Sakai, Nishikawa, Leyton, Benkelfat, Young and Diksic2006). Finally, the val158met polymorphism in COMT is located on 22q11.2 and transcribes an enzyme that metabolizes catecholamine, which serves to modulate the transmission of serotonin, norepinephrine, and dopamine (Cumming, Brown, Damsma, & Fibiger, Reference Cumming, Brown, Damsma and Fibiger1992). The valine allele leads to a three- to fourfold increase in activity of the catecholamine enzyme (compared to the methionine allele; Zubieta et al., Reference Zubieta, Heitzeg, Smith, Bueller, Xu, Xu and Goldman2003) and has been most commonly found to be associated with ASB (e.g., DeYoung et al., Reference DeYoung, Getchell, Koposov, Yrigollen, Haeffel, Klinteberg and Grigorenko2010; Qian et al., Reference Qian, Liu, Wang, Yang, Guan and Faraone2009). Based on recent empirical precedent, each genetic variant was coded according to the number of putative risk alleles in relation to ASB (e.g., Thibodeau et al., Reference Thibodeau, Cicchetti and Rogosch2015).

Statistical analysis

Step 1. Growth mixture models (GMM) of ASB

ASB sum scores were modeled longitudinally using GMM in Mplus 7.3 (Muthen & Muthen, Reference Muthén and Muthén2015). GMM is a group-based analytic method that identifies subpopulations characterized by their observed trajectories. In contrast to latent growth curve analysis, GMM allows for within-class variation of the growth parameters (Muthen & Muthen, Reference Muthén and Muthén2015). To model growth trajectories from adolescence and into adulthood, data were restructured such that time was represented by age rather than by wave, resulting in a large amount of “missing data by design” (Little, Reference Little2013). Missingness was accounted for by full information maximum likelihood estimation in Mplus. In order to account for the nonnormality of ASB sum scores in the GMM (Walters & Ruscio, Reference Walters and Ruscio2013), a zero-inflated Poisson model was specified such that two growth models were estimated: the first growth model described the count part of the outcome for all individuals who are able to assume values of zero and above, and the second growth model described the inflation part of the outcome (i.e., the probability of being able to assume any value except zero; Muthen & Muthen, Reference Muthén and Muthén2015). To account for the nonlinear growth of ASB across this age range (i.e., Bongers, Koot, van der Ende, & Verhulst, Reference Bongers, Koot, van der Ende and Verhulst2003; Odgers et al., Reference Odgers, Moffitt, Broadbent, Dickson, Hancox, Harrington and Caspi2008; Thompson & Tabone, Reference Thompson and Tabone2010), a quadratic function was specified in the model. The optimal number of classes was determined based the Bayesian information criterion (BIC) and the Lo–Mendell–Rubin likelihood ratio test (LRT; Lo, Mendell, & Rubin, Reference Lo, Mendell and Rubin2001). Simulation studies have shown that BIC and LRT indices are more reliable indicators of the “true” number of classes in GMM (Nylund, Asparouhov, & Muthen, Reference Nylund, Asparouhov and Muthén2007). Models that have lower BIC values indicate better fit to the data, whereas LRT provides a hypothesis test where a significant p value (p < .05) indicates rejection of the null hypothesis (k classes – 1) in favor of a model with at least k classes (Lo et al., Reference Lo, Mendell and Rubin2001). In addition to examining traditional measures of fit, models were also evaluated on the basis of their interpretability (i.e., meaningfulness of classes and consistency with existing theories). Regarding potential sex differences, although prevalence rates are known to differ between males and females, developmental pathways of ASB have not been found to differ by sex (Côté, Zoccolillo, Tremblay, Nagin, & Vitaro, Reference Côté, Zoccolillo, Tremblay, Nagin and Vitaro2001; Mazerolle, Brame, Paternoster, Piquero, & Dean, Reference Mazerolle, Brame, Paternoster, Piquero and Dean2000; Moffitt, Caspi, Harrington, & Milne, Reference Moffitt, Caspi, Harrington and Milne2002; Odgers et al., Reference Odgers, Moffitt, Broadbent, Dickson, Hancox, Harrington and Caspi2008). Thus, to reduce concerns due to multiple testing, the GMM was conducted for the full genotypic sample rather than separately for each sex.

Step 2. Calculating Add Health MRS for ASB

Stringent tests of rGE were conducted in order to rule out a confounding effect on the G × E models. In line with empirical precedent (i.e., Thibodeau et al., Reference Thibodeau, Cicchetti and Rogosch2015), each genotype was regressed on supportive parenting, school connectedness, and covariates (i.e., household income and peer closeness). Polymorphisms that were significantly associated with either environment variable were excluded from the Add Health MRS (Thibodeau et al., Reference Thibodeau, Cicchetti and Rogosch2015).

After rGE confounds were accounted for, k-fold cross-validation (Kohavi, Reference Kohavi1995) was used to derive robust prediction estimates of ASB for each polymorphism. The data were divided randomly into k equal-sized partitions. By convention, k = 5 was used for the current analysis (Pereira, Mithcell, & Botvinick, Reference Pereira, Mitchell and Botvinick2009). For each k-fold, a model was fit to the data (the model is described below) using the other k – 1 partitions (i.e., the “testing” set), with the fold left out used as the validation set. More explicitly, each testing set was composed of a different 80% of the sample (N = 7,067) from the total sample (N = 8,834), while each validation set consisted of the remaining 20%. In total, each of the five validation sets was composed of a different 20% of the total sample so that there was no overlap between any of the validation sets. For each of the five testing sets, beta regression models were conducted (and subsequently cross-validated) to examine associations between each polymorphism and the most severe GMM class that emerged (i.e., persistent ASB). This outcome was selected because the persistent ASB pathway has a high heritability relative to other ASB pathways (Burt, Reference Burt2009; Zheng & Cleveland, Reference Zheng and Cleveland2015), and genetic influences for aggressive and persistent forms of ASB may be partially shared with other forms of ASB (Burt, Reference Burt2013) as well as externalizing problems more generally (Lahey, Van Hulle, Singh, Waldman, & Rathouz, Reference Lahey, Van Hulle, Singh, Waldman and Rathouz2011). Separate cross-validations were conducted for males and females because the beta regression for females did not include MAOA genotype, given that MAOA is located on the X chromosome and heterozygosity can only be present in females but not males (Ficks & Waldman, Reference Ficks and Waldman2014). Furthermore, inferences about MAOA expression among heterozygous female carriers has yet to be determined, as few studies of ASB have included this group in their investigations (Carrel & Willard, Reference Carrel and Willard2005; Ficks & Waldman, Reference Ficks and Waldman2014; Kim-Cohen et al., Reference Kim-Cohen, Caspi, Taylor, Williams, Newcombe, Craig and Moffitt2006; Li et al., Reference Li, Berk and Lee2013). Prediction estimates for each polymorphism were calculated by averaging the standardized beta coefficients from each of the five folds (Cordell, Reference Cordell2009), which was then used to calculate Add Health MRS from the following formula (Derringer et al., Reference Derringer, Krueger, Dick, Saccone, Grucza, Agrawal and Nurnberger2010),

$$\eqalign{{\rm Add\ Health \,MRS} = \; &({\rm Genotyp}{\rm e}_{\rm 1}\,{\rm \times}\,{\rm \beta }_{\rm 1}) + ({\rm Genotyp}{\rm e}_{\rm 2}\,{\rm \times}\,{\rm \beta }_{\rm 2}) \cr \noalign {\vskip-18pt} \cr &+ \cdots + ({\rm Genotyp}{\rm e}_{\rm i}\,{\rm \times}\,{\rm \beta }_i),}$$

where genotypei was dummy-coded according to each putative risk allele (see Table 1 for the risk allele for each genotype) and β i was the average standardized beta coefficient for each genotype from the k-fold cross-validation analysis. As specified by others (Marcaeu et al., 2016), the p value threshold was set to 1 so that all polymorphisms were included into the score regardless of its significance. Thus, higher Add Health MRS indicates greater risk and lower scores suggests lower risk for ASB.

Table 1. Allele frequencies and Hardy–Weinberg equilibrium (HWE) tests

Step 3. Testing G × E Interactions

After the best fitting model from the GMM was determined and the Add Health MRS was derived, multinomial logistic regressions were conducted separately for males and females to examine the association between Add Health MRS, environmental factors (i.e., supportive parenting and community support), and their interactions on the log-likelihood of ASB class membership relative to the most prevalent (i.e., low) class. All models included a stratification variable (region), sample weights (GSWGT1), and the sampling unit variable (PSUSCID), while also covarying for the effects of household income and peer closeness. Per expert recommendation, all covariate × genotype and covariate × environment cross products were included simultaneously to account for the effects of the covariate on the primary interactions of interest (Keller, Reference Keller2014).

Results

Demographic comparison, principal component analysis, and rGE

Table 2 displays the mean comparisons between sexes for each of the primary variables in the study. Females were younger in age and reported fewer counts of substance use and ASB across each wave, relative to males. Females also scored lower on supportive parenting and reported slightly higher school connectedness relative to males. There were no sex differences in terms of household income. Genotype distributions did not deviate significantly from Hardy–Weinberg equilibrium (Table 1).

Table 2. Demographic information of the Caucasian-only genotypic sample by sex

Note: Substance use behaviors included cigarette smoking, heavy alcohol consumption, and marijuana use.

A principal component analysis was performed on all supportive parenting and school connectedness items from the Wave 1 in-home survey data of Add Health (Table 3). Three components had eigenvalues greater than 1, although the first two components accounted for 44% of the total variance and the remaining component only accounted for 5% of the variance. As expected, the first two components differentiated between the items that pertained to school connectedness and supportive parenting, suggesting that both environmental variables constitute separable components for the present analysis. Items with high loadings (i.e., loadings >0.30) on Components 1 and 3 were the supportive parenting items (e.g., “close to mom” and “good communication”). These two components were also strongly correlated (r = .284, p < .001), indicating that they may be reasonably tapping into a similar factor. Component 2 consisted of high loadings on items exclusively related to school connectedness (e.g., “feel close to people at school” and “teachers care about you”).

Table 3. Principal component analysis of supportive parenting and school connectedness variables in Add Health (oblimin rotation)

Tests of rGE were conducted using multiple linear regression, where each candidate gene was regressed on supportive parenting and school connectedness (as separate models). For males, no significant associations were detected between supportive parenting and DAT1 (β = 0.006, p = .75), DRD4 (β = 0.007, p = .72), DRD2 (β = 0.006, p = .73), COMT (β = 0.001, p = .98), 5-HTTLPR (β = 0.019, p = .30), or MAOA (β = –0.019, p = .29). There were also no significant associations detected for school connectedness and DAT1 (β = 0.02, p = .39), DRD4 (β = 0.015, p = .39), DRD2 (β = –0.010, p = .59), 5-HTTLPR (β = 0.018, p = .30), or MAOA (β = –0.002, p = .91), although a significant association was found for COMT (β = –0.039, p = .03). For females, no significant associations were detected between supportive parenting and DAT1 (β = 0.005, p = .75), DRD4 (β = –0.011, p = .53), COMT (β = 0.007, p = .67), or 5-HTTLPR (β = 0.003, p = .88), but there was an association for DRD2 (β = –0.044, p = .01). There were no significant associations detected for school connectedness and DAT1 (β = –0.025, p = .15), DRD4 (β = –0.003, p = .86), DRD2 (β = –0.035, p = .66), 5-HTTLPR (β = –0.031, p = .07), or COMT (β = 0.020, p = .25). Because of the possibility of an rGE confound, COMT and DRD2 were excluded from the Add Health MRS computation for males and females, respectively.Footnote 3

GMM of ASB and sex differences

Growth curves of ASB were modeled from age 13 to 32 using GMM (Figure 1). A four-class solution was the optimal model (Table 4). According to the model, most of the sample (53.2% total; 62.6% of females, 42.56% of males) belonged in the low pathway, where rates of mean counts of ASB were uniformly low across the entire age range. An adolescence-peaked ASB pathway also emerged, consisting of 23.1% of the sample (23.9% of females, 22.3% of males). Individuals belonging to this class were characterized by relatively high rates of ASB during middle adolescence (between the ages of 13 and 16) and a steep decline beginning at age 17. By age 22, individuals in the adolescence-peaked pathway were indistinguishable from their low-ASB peers in terms of mean counts of ASB. The persistent ASB pathway represented 4.0% of the sample (2.0% of females, 6.3% of males) and was characterized by a relatively moderate but persistent trend of ASB throughout adolescence and adulthood. The final class that emerged was a high decline ASB pathway, representing 19.7% of the sample (11.5% of females, 28.8% of males). Similar to the adolescent-peaked pathway, high decline ASB individuals showed the highest initial rate of ASB during adolescence, but unlike those in the adolescence-peaked pathway, they exhibited a more gradual decline in mean counts of ASB that extended well into their adult years. Consistent with prior studies (e.g., Ferguson, Reference Ferguson2010; Moffitt et al., Reference Moffitt, Caspi, Harrington and Milne2002), there was a significant overall sex difference in ASB pathway membership (χ2 = 605.42, df = 3, p < .01), such that males were more likely to belong in the persistent and high decline ASB pathways than females, females were more likely to be in the low ASB pathway than males, and males and females did not differ on prevalence of adolescent-peaked pathway membership.

Figure 1. Estimated means of the latent trajectories of antisocial behavior.

Table 4. Fit statistics for growth mixture model of ASB sum scores from age 13 to 32

Note: AIC, Akaike information criterion; BIC, Bayesian information criterion; LMR-LRT, Lo–Mendell–Rubin likelihood ratio test for k classes (null) versus k + 1.

Add Health multilocus genetic risk scores

Each genotype was beta regressed on the probability of persistent ASB, controlling for household income and peer closeness, in a fivefold cross-validation analysis. In males, DRD4 genotype was significantly associated with the probability of persistent ASB in four of the five folds, such that carrying a 7R allele was negatively associated with ASB. In females, none of the genotypes conferred a significant main effect on the probability of the persistent ASB pathway. All results, including tables for the male and female cross-validation models, are reported in the online-only supplementary materials.

Prediction estimates used to calculate Add Health MRS were derived from the average of the standardized beta coefficients across the fivefold cross-validation analysis. Tables S.1 and S.2 in the online-only supplementary materials include the histogram of weighted Add Health MRS for males and females. Table 5 presents the results of discrete probabilities of ASB class membership, mean supportive parenting, and mean school connectedness when regressed on weighted Add Health MRS. As expected, Add Health MRS were positively associated with the probability of persistent ASB class membership, but not probabilities for membership in low, adolescence-peaked ASB, or high decline ASB classes. As an additional test of rGE, Add Health MRS were not associated with either supportive parenting or school connectedness, suggesting that the individual's genetic liability for ASB was not simultaneously predictive of the quality of their own environment.Footnote 4

Table 5. Regressions results from each k-fold cross-validation: Add Health multigenic risk scores predicting discrete probabilities of latent antisocial behavior class membership, supportive parenting, and school connectedness by sex

Note: Estimates as shown are standardized beta coefficients.

*p < .05. **p < .01.

Add Health MRS × Supportive Parenting interaction

Multinomial logistic regressions were conducted separately for males and females, predicting the log-odds of ASB pathway membership (relative to low ASB) from Add Health MRS, supportive parenting, and their interaction and controlling for household income, average substance use behavior across waves, and peer closeness (see Table 6). For males, no significant interaction was detected between supportive parenting and Add Health MRS in predicting membership in the adolescence-peaked, persistent, and high decline ASB pathways relative to the low ASB pathway. Similarly, in females no significant interactions were detected between supportive parenting and Add Health MRS in predicting membership in the adolescence-peaked, persistent, and high decline ASB pathways compared to the low ASB pathway.

Table 6. Multinomial logistic regressions: Supportive parenting and multigenic risk scores (MRS)

Note: The low class served as the reference class for the multinomial regressions.

MRS × School Connectedness interaction

Multinomial logistic regressions were also conducted to predict the log-odds of ASB class membership from Add Health MRS, school connectedness, and their interaction (see Figure 2), controlling for household income, substance use behavior, and peer closeness (see Table 7). In males, a significant interaction was detected between school connectedness and Add Health MRS in predicting membership in the adolescence-peaked ASB pathway (logits = –6.78, SE = 3.30, p = .05) relative to the low class.Footnote 5 In females, no significant interactions emerged between school connectedness and Add Health MRS for females in predicting membership in any of the pathways of the ASB.

Figure 2. School Connectedness × Add Health Multigenic Risk Score interaction predicting odds of adolescence-peaked pathway antisocial behavior membership (relative to low antisocial behavior) in males.

Table 7. Multinomial logistic regressions: School connectedness and multigenic risk scores (MRS)

Note: The low class served as the reference class for the multinomial regressions.

Using a median split of Add Health MRS (median = 0), post hoc tests of the interaction between school connectedness and Add Health MRS revealed that the association between school connectedness and the log odds of membership in the adolescence-peaked pathway was robust for high (logits = –0.33, SE = 0.08, p < .001, n = 248) and low Add Health MRS (logits = –0.32, SE = 0.07, p < .001, n = 236) males. To assess for differential susceptibility, odds ratios and confidence intervals were assessed “at the margins” of school connectedness (lowest and highest values of the z score) comparing high and low Add Health MRS groups. Nonoverlapping confidence intervals would indicate a significant difference in odds ratios between groups. No evidence of differential susceptibility emerged. High Add Health MRS males who reported the lowest levels of school connectedness (–4 SD) had a 41% increased odds of being in the adolescence-peaked pathway compared to low ASB pathway, odds ratio (OR) = 1.41, 95% confidence interval (CI) [1.13, 1.49], and low Add Health MRS males at the lowest level of school connectedness had no significant increase in odds of belonging to this pathway relative to the low ASB pathway, OR = 1.00, 95% CI [0.89, 1.21]. These odds did not significantly differ from one another, however, given their nonoverlapping CIs. At high levels of school connectedness (+2 SD), both high and low Add Health MRS males were equally protected from being in the adolescence-peaked pathway relative to the low pathway: high Add Health MRS, OR = 0.15, 95% CI [0.04, 0.35]; low Add Health MRS, OR = 0.28, 95% CI [0.14, 0.73].

Discussion

The current study used rigorous quantitative methods to examine the interplay between multigenic and multienvironmental influences on ASB. First, four unique pathways of ASB emerged from the growth mixture analysis: low, adolescence-peaked, persistent, and a high decline ASB classes. The high-decline class accounted for nearly 20% of the total sample. Second, the k-fold cross-validation technique was used to derive robust prediction estimates for ASB across six polymorphisms in dopamine and serotonin genes, such that weighted MRS were predictive of the persistent pathway of ASB in both males and females. Third, a significant G × E effect was detected for male pathways of ASB. Specifically, compared to males with low Add Health MRS, males with high Add Health MRS had somewhat greater odds of adolescence-peaked pathway membership at low levels of school connectedness. No evidence of differential susceptibility was detected, however, such that males with high and low MRS were not differentially protected from membership into the adolescence-peaked pathway at high levels of school connectedness. Furthermore, no interactive effects between MRS and supportive parenting was detected in predicting the developmental pathways of ASB for either males or females.

The adolescent to adult pathways for ASB were highly consistent with those described in the broader developmental literature (D'Unger, Land, & McCall, Reference D'Unger, Land and McCall2002; Moffitt, Reference Moffitt1993; Odgers et al., Reference Odgers, Moffitt, Broadbent, Dickson, Hancox, Harrington and Caspi2008) as well as to those described using a limited subsample from Add Health (e.g., Tung & Lee, Reference Tung and Lee2016; Zheng & Cleveland, Reference Zheng and Cleveland2013). As hypothesized, a high decline ASB pathway was also identified, representing nearly 20% of the total sample. Similar to the adolescence-peaked pathway of ASB, individuals in the high decline pathway exhibited peak rates of ASB during adolescence, but they also had the highest initial rate of ASB (even compared to the persistent class) and a much slower decline in their ASB throughout adulthood relative to the adolescence-peaked pathway. The existence of the high decline ASB class had been previously found for both males and females in a small Add Health subsample and other samples as well, despite focusing only on the first three waves of data (Nagin & Tremblay, Reference Nagin and Tremblay2005; Zheng & Cleveland, Reference Zheng and Cleveland2013, Reference Zheng and Cleveland2015). Prevalence rates were found to be roughly equivalent between sexes in a previous study (12% vs. 11% for males and females, respectively; Zheng & Cleveland, Reference Zheng and Cleveland2013), but the current results differed from the previous finding in that nearly 29% of males belonged in the high decline pathway versus only 12% of females. One reason for the divergence in these results is that the pathways combined aggressive and nonaggressive rule-breaking forms of ASB, which were undifferentiated in the current study.Footnote 6 For instance, Zheng and Cleveland (Reference Zheng and Cleveland2013) reported that high decliner males had consistently high rates of nonaggressive rule breaking and aggressive delinquency over time, whereas female high decliners were primarily characterized by their high nonaggressive rule-breaking behaviors that slowly declined over time. The higher proportion of males in the high decliner pathway may potentially be accounted for by some males who are exclusively aggressive in their ASB, which is very rare in females (Moffitt, Reference Moffitt1993). Although there was no factor analytic evidence of subtypes across time points (i.e., by age group; results available upon request), subsequent studies that differentiate between aggressive and nonaggressive forms of ASB may yet identify unique, potentially sex-specific developmental pathways of ASB.

Weighted MRS for ASB, which were derived from robust prediction estimates of each dopamine and serotonin polymorphism, significantly predicted the persistent pathway of ASB but not the other developmental pathways of ASB. First, it is important to note that the choice of phenotype in the genetic association cross-validation models likely had a significant impact on predictive ability of the resultant Add Health MRS in the full sample (Wray et al., Reference Wray, Yang, Hayes, Price, Goddard and Visscher2013). The current study used the probability of persistent ASB as the outcome choice given that (a) the genetic influences for chronic and aggressive forms of ASB are highest relative to other developmental pathways (Burt, Reference Burt2009; Zhen & Cleveland, 2015) and (b) genetic influences for the persistent ASB pathway may overlap (partially) with other (i.e., nonaggressive) ASB pathways (Burt, Reference Burt2013) and with externalizing problems more generally (Lahey et al., Reference Lahey, Van Hulle, Singh, Waldman and Rathouz2011). Although Add Health MRS significantly predicted the persistent ASB pathway, it did not predict adolescence-peaked or high decline ASB pathways, providing some evidence of a genetic distinction between the developmental pathways, but only with respect to the dopaminergic and serotonergic polymorphisms that were assessed. Thus, it cannot be ruled out that the current results would differ if a different phenotype were selected (e.g., antisocial personality disorder diagnosis and criminal arrest record). As related phenotypes should theoretically tap into similar genetic signals that underlie a single latent trait, future studies should compare the predictive ability of Add Health MRS and assess for the accuracy of the prediction across these related (but different) phenotypes. Second, the current study referred to the extensive literature regarding the risk variant for each polymorphism as it pertained to ASB, yet the results showed that some risk variants were negatively associated with ASB. For example, a negative association between the DRD4 7R allele and the probability of persistent ASB was detected for males, such that carrying more copies of the 7R allele predicted a lower probability of the persistent ASB. Although this contrasts with a large portion of prior literature regarding the risk status of the 7R allele on ASB (e.g., Buchmann et al., Reference Buchmann, Zohsel, Blomeyer, Hohm, Hohmann, Jennen-Steinmetz and Esser2014; Hohmann et al., Reference Hohmann, Becker, Fellinger, Banaschewski, Schmidt, Esser and Laucht2009; Nobile et al., Reference Nobile, Giorda, Marino, Carlet, Pastore, Vanzin and Battaglia2007), some studies found that the 7R allele was associated with less ASB in adolescent samples as well (Beaver et al., Reference Beaver, Wright and Walsh2008; Janssens et al., Reference Janssens, Goossens, Van Den Noortgate, Colpin, Verschueren and Van Leeuwen2015). Studies making assumptions about the risk effect of a certain allele may obfuscate the possibility of a protective effect, particularly given the evidence that serotonergic and dopaminergic systems may have reciprocal effects in relation to ASB (Seo et al., Reference Seo, Patrick and Kennealy2008). Hence, a quantitative approach (e.g., cross-validation analysis) to genetic scoring may potentially provide a better understanding about the interplay between serotonin and dopaminergic function underlying ASB, since the direction and effect size for each SNP can be estimated from a model that simultaneously (rather than independently) incorporates all of the relevant markers in the prediction of ASB. Third, the overall genetic effects on ASB might have been “washed out” given the contrasting direction of allelic effects on ASB in the genetic association cross-validation analysis. Furthermore, others have noted that the combination of positive and negative genetic effects is more biologically realistic (Marcaeu et al., 2016). That this composite score still predicted a certain form of ASB suggests that aggregate genetic effects on ASB are robust.

Add Health MRS interacted with school connectedness in predicting membership into the adolescence-peaked pathway relative to the low ASB pathway. This interaction was evident even after accounting for the possibility of rGE confounds, and controlling for household income, substance use behavior (averaged across adolescence and adulthood), and peer closeness. Furthermore, secondary analysis indicated that the effect of DRD4 genotype on ASB may have been a significant driver of this result, but the weighted Add Health MRS yielded an even stronger signal to detect the G × E effect. Overall, the finding reinforces the pivotal role of nonfamilial influences on ASB during the transitory period of development (i.e., from early adolescence into adulthood) and demonstrate how school-related experiences (independent from peer closeness) during adolescence may have cumulative interactive effects among those with a vulnerable genetic constitution for not only the adolescence-peaked form of ASB, as might be expected from prior literature (Moffitt, Reference Moffitt1993, 2002). Adolescents who feel unaccepted and unsupported at school may be less likely to follow social norms and expectations, resulting in a higher likelihood of engaging in delinquency and aggression into their later adult years (Catalano & Hawkins, Reference Catalano, Hawkins and Hawkins1996). Genetic influences for ASB are likely strengthened under these conditions, as prior behavioral genetic studies have found that genetic factors for externalizing problems (including substance use and ASB) are increased under conditions where the individual has greater access to and affiliation with deviant peers (Harden, Hill, Turkheimer, & Emery, Reference Harden, Hill, Turkheimer and Emery2008; Samek, Hicks, Keyes, Iacono, & McGue, Reference Samek, Hicks, Keyes, Iacono and McGue2016) and is living in a more disadvantaged neighborhood (Burt, Klump, Gorman-Smith, & Neiderhiser, Reference Burt, Klump, Gorman-Smith and Neiderhiser2016; Tuvblad, Grann, & Lichtenstein, Reference Tuvblad, Grann and Lichtenstein2006).

In addition, the interaction between Add Health MRS and school connectedness on adolescence-peaked ASB was only observed for males. There is evidence that genetic and environmental influences differentially contribute to developmental pathways of ASB across sexes (Van Hulle, Rodgers, D'Onofrio, Waldman, & Lahey, Reference Van Hulle, Rodgers, D'Onofrio, Waldman and Lahey2007; Zheng & Cleveland, Reference Zheng and Cleveland2015). Males and females may have different patterns of socialization that affect their engagement in ASB, and others have suggested that socialization processes and familial factors may be more relevant for females than for males (Fontaine, Carbonneau, Vitaro, Barker, & Tremblay, Reference Fontaine, Carbonneau, Vitaro, Barker and Tremblay2009; Zheng & Cleveland, Reference Zheng and Cleveland2015). Hence, the male-specificity of the G × E effect was unexpected, as prior studies have shown that adolescent girls may be more sensitive to social evaluation, acceptance, and connections with teachers and peers as it relates to general mental health outcomes than adolescent boys (Calvete & Cardenoso, Reference Calvete and Cardeñoso2005; Shochet et al., Reference Shochet, Dadds, Ham and Montague2006). Findings have been inconsistent regarding the differential effects of school connectedness for males and females as it pertains to ASB as well (see Batanova & Loukas, Reference Batanova and Loukas2014; Crosnoe et al., Reference Crosnoe, Erickson and Dornbusch2002; Kuperminc & Allen, Reference Kuperminc and Allen2001). Overall, the current findings provide preliminary support for the broader literature regarding the differential contribution of genetic and environmental influences on developmental pathways of ASB across sex, but results should be replicated before making firm conclusions regarding the male-specific G × E effects that were detected.

Contrary to hypotheses, there was no evidence of differential susceptibility, as low levels of perceived school connectedness exacerbated the risk for the adolescence-peaked pathway of ASB among males with high Add Health MRS, but these individuals were not differentially protected by high levels of school connectedness compared to males with low Add Health MRS. First, it should be noted that relatively few studies of differential susceptibility have focused on environments assessed during the adolescent period (for exceptions, see Belsky & Beaver, Reference Belsky and Beaver2011; Latendresse et al., Reference Latendresse, Bates, Goodnight, Lansford, Budde, Goate and Dick2011; Li, Berk, & Lee, Reference Li, Berk and Lee2013; Tung & Lee, Reference Tung and Lee2016), and whether these effects might influence developmental trajectories of ASB measured into adulthood. A large portion of the genetically informed studies of differential susceptibility have focused on early rearing environments, assessed during infancy and the school-age years (see Bakermans-Kranenburg & van IJzendoorn, Reference Bakermans-Kranenburg and van Ijzendoorn2011), but one meta-analysis that focused on 5-HTTLPR did not find any evidence that age moderated the association between the sensitivity genotype (i.e., short–short and short–long carriers) and various developmental outcomes as a function of either adverse or enriched environments (van IJzendoorn, Belsky, & Bakermans-Kranenburg, Reference van IJzendoorn, Belsky and Bakermans-Kranenburg2012). In the case of ASB specifically, other genes may come online at the later stages of development (i.e., genetic innovation; Pingault et al., Reference Pingault, Viding, Galéra, Greven, Zheng, Plomin and Rijsdijk2015, Wichers, Wigman, & Myin-Germeys, Reference Wichers, Wigman and Myin-Germeys2015), but what these genes are, and how these genes operate in the context of G × E is currently unknown. Second, despite having assessed parenting and school connectedness across broad dimensions (i.e., negative to positive support), differential susceptibility in G × E may be difficult to detect nonetheless due to measurement artifacts and sample size. In the current study, it is possible that high levels of school connectedness or supportive parenting may not have been sufficiently rewarding or salient for dopamine activation (Tomer et al., Reference Tomer, Slagter, Christian, Fox, King, Murali and Davidson2014; Tripp & Wickens, Reference Tripp and Wickens2008). This may also explain why studies examining highly aversive environmental stimuli (e.g., stressful life events and maltreatment) have found relatively more consistent G × E effects (see meta-analyses by Byrd & Manuck, Reference Byrd and Manuck2014; Kim-Cohen et al., Reference Kim-Cohen, Caspi, Taylor, Williams, Newcombe, Craig and Moffitt2006; Tielbeek et al., Reference Tielbeek, Karlsson Linnér, Beers, Posthuma, Popma and Polderman2016). A consideration for future studies is to rigorously measure the impact of the environmental criterion in G × E, insofar as the salience of effect at either end of the environment spectrum has a large and inverse effect on the phenotype. Experimental studies of G × E, for example, are powerful tests for differential susceptibility because the salience of the environment (i.e., intervention effects) is manipulated and directly measurable by the researchers (Brody, Yu, Beach, & Philibert, Reference Brody, Yu, Beach and Philibert2015). Another possible explanation for the lack of differential susceptibility effect is that deriving weighted genetic risk scores for ASB may have precluded identification of genetic effects that are contingent upon exposure to environmental conditions (Caspi, Hariri, Holmes, Uher, & Moffitt, Reference Caspi, Hariri, Holmes, Uher and Moffitt2010). Robust genotype-to-phenotype associations detected in studies may not readily guarantee robust G × E findings. One future direction is to characterize genetic variation in phenotypes that are plausibly related to how an individual reacts to environmental factors (e.g., stress sensitivity and social vulnerability; Caspi et al., Reference Caspi, Hariri, Holmes, Uher and Moffitt2010).

Supportive parenting did not moderate genetic risk for developing ASB, which is consistent with behavior genetic findings of reduced familial environmental influences, relative to nonfamilial environmental influences on ASB as a function of age (Ferguson, Reference Ferguson2010; Rhee & Waldman, Reference Rhee and Waldman2002). The findings contrast with at least two previous reports demonstrating the moderating role of supportive parenting on genetic susceptibility for ASB-related phenotypes (i.e., Belsky & Beaver, Reference Belsky and Beaver2011; Tung & Lee, Reference Tung and Lee2016). Both Tung and Lee (Reference Tung and Lee2016) and Belsky and Beaver (Reference Belsky and Beaver2011) used the Add Health sample in their studies as well, but they also employed fewer waves of data, a different method of measuring genetic susceptibility, and a different outcome measure compared to the current study. These are important study differences that may explain nonreplication, but it should also be noted that current study undertook very rigorous steps toward measuring each component of the G × E equation to ensure that the findings (or lack thereof) were robust, including a systematic approach for combining the effects of functionally relevant polymorphisms, assessment of positive and negative dimensions of familial and nonfamilial environmental influences on development, and the use of GMM to assess for age 13–32 trajectories of ASB. The current study focused on a limited set of polymorphisms as well, and it is possible that familial factors may interact with other genetic variants as they pertain to ASB that were not assessed here. Understanding interplay between familial and genetic influences remains a high priority for scientific inquiry because of its potential implications for intervention for adolescents. For instance, Brody, Chen, and Beach (Reference Brody, Chen and Beach2013) conducted a randomized control trial of adolescents involved in a family-centered alcohol prevention program and found that youths who carried the A1 and 7R alleles in the DRD2 and DRD4 genes, respectively, had the greatest increase in alcohol use during a 2-year period if they were in the control condition but were also the most likely to benefit from the prevention program (by virtue of having the smallest increase in alcohol use behaviors) compared to those with fewer risk alleles.

Several study limitations should be noted, in addition to those that were previously mentioned. First, only Caucasians were analyzed in the current analysis, which limits the generalizability of the current findings. However, the advantage of a homogenous racial/ethnic sample is that it minimizes concerns related to population stratification (genotypes were nonrandomly associated with race/ethnicity in the present sampleFootnote 7 ). Future studies should consider the strong possibility that genetic and environmental influences may be specific to certain populations. For example, in the substance use literature recent evidence from GWAS suggest that genetic variants for alcohol dependence differ between African Americans and Caucasians (Gelernter et al., Reference Gelernter, Kranzler, Sherva, Almasy, Koesterer, Smith and Wodarz2014). Although G × E investigations focusing on non-Caucasian samples are still quite rare in the literature (although see Thibodeau et al., Reference Thibodeau, Cicchetti and Rogosch2015, for a recent exception), it is possible that G × E effects may be differ depending on the population being assessed. Second, the Add Health MRS was not cross-validated outside of the Add Health sample (which is why the term “Add Health MRS” rather than simply MRS was used throughout). As such, the Add Health MRS may not predict persistent ASB in a different sample. While the k-fold cross-validation does provide safeguards against overfitting, concerns about the generalizability of the Add Health MRS are warranted in the absence of an independent validation sample. Third, substance use behaviors are also known to covary with dimensions of ASB, such that they may have shared etiologies (Kendler, Prescott, Myers, & Neale, Reference Kendler, Prescott, Myers and Neale2003; Krueger, Markon, Patrick, Benning, & Kramer, Reference Krueger, Markon, Patrick, Benning and Kramer2007). The current investigation included substance use behavior as a covariate in the multinomial logistic models, which partially accounts for the covariation of substance use and ASB pathways in the interpretation of G × E. However, there may be unique genetic and environmental influences on pathways unique to substance use, as well for a “general” factor of externalizing psychopathology that may account for the covariation between the different externalizing dimensions (Caspi et al., Reference Caspi, Houts, Belsky, Goldman-Mellor, Harrington, Israel and Moffitt2014). Fourth, childhood data on ASB were not collected on the sample, since Add Health began collecting data on participants when they were in high school. The absence of this data makes it difficult to compare the pathways derived from the current study to those theorized by Moffitt (Reference Moffitt1993, Reference Moffitt, Lahey, Moffitt and Caspi2003). Fifth, the reliance on self-report for the variables in the study may have confounded the interpretations of the results due to shared method variance. Because ASB was self-reported at each wave, there may have also been underreporting or a biased sample at each subsequent wave since individuals with severe ASB tend to be more difficult to be retained in longitudinal studies (Farrington, Reference Farrington1998; Huizinga & Elliott, Reference Huizinga and Elliott1986). Hence, the use of additional informants or using other sources of ASB, such as official criminal records, may have improved measurement accuracy (Jaffee, Caspi, Moffitt, & Taylor, Reference Jaffee, Caspi, Moffitt and Taylor2004). Sixth, the polymorphisms examined in the current study were limited by what was publicly available in Add Health. A larger set of genotypes with relevance to dopamine and serotonin neurotransmission would have been preferred (Plomin et al., Reference Plomin, Haworth and Davis2009). Similarly, genome-wide information would help to identify several more unique genetic variants (perhaps outside of the dopamine or serotonin system) that may, in aggregate, account for a significant portion of the phenotypic variance in the ASB pathways. Although genome-wide information on the Add Health data set is not yet publicly available (i.e., via dbGap), future studies should consider the interactive role of polygenic risk (derived from GWAS) and environmental influences on the development of ASB and related phenotypes.

The current study provides several important contributions to the literature and directions for future study. First, the findings illustrate the utility of integrating a strong developmental framework toward understanding the role of multiple genetic and environmental influences on pathways of ASB. In psychiatric genetics, a developmental approach is crucial for the study of many psychiatric disorders because genetic/environmental associations for a particular age group may not be generalizable to all age groups (Dick, Reference Dick2011). Employing longitudinal strategies, such as GMM, allows researchers to examine how genetic and environmental influences predict different patterns of development, which is a critical area of focus in the field of developmental psychopathology. Second, this study builds upon the advances in the field of psychiatric genetics by computing weighted genetic risk scores for ASB from a set of previously well-characterized polymorphisms in dopamine and serotonin genes. This approach provided a strong signal to detect a G × E effect and differs from the increasingly common strategy of summing the number of putative risk alleles for each marker by accounting for the differential effects that each marker has on ASB. GWAS-based investigations can now be used to test for enrichment across biological systems in relation to an outcome, potentially providing more targeted investigations for G × E (see Edwards, Bacanu, Bigdeli, Moscati, & Kendler, Reference Edwards, Bacanu, Bigdeli, Moscati and Kendler2016, for a recent example of dopamine-gene pathway analysis for schizophrenia in a GWAS). Biologically based polygenic approaches will likely become more widely implemented in G × E studies as they potentially provide a better understanding about how genes in a targeted biological system and environments interact to influence behavior. Moving forward, more studies are concurrently needed to validate the multitude of genetic scoring strategies that are developed (Derringer et al., Reference Derringer, Krueger, Dick, Saccone, Grucza, Agrawal and Nurnberger2010; Marcaeu et al., 2016).

Supplementary Material

To view the supplementary material for this article, please visit https://doi.org/10.1017/S0954579417001511.

Footnotes

The author was supported in part by a core grant to the Waisman Center from the National Institute of Child Health and Human Development (P30-HD03352). The author also acknowledges contributions from Ben Domingue and Garvesh Raskutti for their assistance in the preparation of this manuscript, including their analytic expertise and critical feedback. This research uses data from the National Longitudinal Study of Adolescent to Adult Health (Add Health), a program project directed by Kathleen Mullan Harris and designed by J. Richard Udry, Peter S. Bearman, and Kathleen Mullan Harris at the University of North Carolina at Chapel Hill, and funded by Grant P01-HD31921 from the Eunice Kennedy Shriver National Institute of Child Health and Human Development, with cooperative funding from 23 other federal agencies and foundations. Special acknowledgment is due Ronald R. Rindfuss and Barbara Entwisle for assistance in the original design. Information on how to obtain the Add Health data files is available on the Add Health website (http://www.cpc.unc.edu/addhealth). No direct support was received from Grant P01-HD31921 for this analysis.

1. Allelic frequencies of each racial/ethnic subgroup are presented online-only supplementary Table S.1.

2. Although studies have traditionally separated the dimensions of ASB into aggression (overt) and rule-breaking (covert) domains, confirmatory factor models comparing unidimensional, two-factor models, and bifactor models across each wave of data consistently showed that unidimensional and bifactor solutions were superior fits to the data, over the two-factor solution. This is also consistent with findings by Tackett, Daoud, De Bolle, and Burt (Reference Tackett, Daoud, De Bolle and Burt2013). Results (including factor loadings and fit indices) are available upon request from the author.

3. Genotypes were not associated with household income or peer influence. Results available upon request.

4. MRS was not significantly associated with household income or peer closeness.

5. Given that a significant effect of DRD4 genotype on the probability of persistent pathway was detected among males in three of the five folds of the genotypic cross-validation analysis, the interaction between DRD4 genotype and school connectedness on ASB pathways was also tested in males using multinomial logistic regression. The interaction between MRS and school connectedness was significant in predicting the adolescence-peaked pathway (logits = –0.28, SE = 0.11, p = .01), but not the persistent (logits = 0.34, SE = 0.22, p = .13) or high decline (logits = 0.13, SE = 0.12, p = .28) pathways relative to low ASB. Hence, DRD4 genotype likely contributed a larger effect on the significant MRS × School Connectedness interaction for male pathways of ASB relative to other genotypes. However, a cumulative weighted risk score that included other genetic variants provided an even stronger signal to detect G × E effects overall.

6. Please see Footnote 1 regarding the theoretical and empirical basis for examining a single dimension of ASB in the current study.

7. See online-only supplementary Table S.2.

References

Asbury, K., Almeida, D., Hibel, J., Harlaar, N., & Plomin, R. (2008). Clones in the classroom: A daily diary study of the nonshared environmental relationship between monozygotic twin differences in school experience and achievement. Twin Research and Human Genetics, 11, 586595.CrossRefGoogle ScholarPubMed
Bakermans-Kranenburg, M. J., & van Ijzendoorn, M. H. (2011). Differential susceptibility to rearing environment depending on dopamine-related genes: New evidence and a meta-analysis. Development and Psychopathology, 23, 3952.CrossRefGoogle ScholarPubMed
Barnes, J. C., Beaver, K. M., & Boutwell, B. B. (2011). Examining the genetic underpinnings to Moffitt's developmental taxonomy: A behavioral genetic analysis. Criminology, 49, 923954.CrossRefGoogle Scholar
Batanova, M., & Loukas, A. (2014). Unique and interactive effects of empathy, family, and school factors on early adolescents’ aggression. Journal of Youth and Adolescence, 43, 18901902.CrossRefGoogle ScholarPubMed
Beaver, K. M., & Belsky, J. (2012). Gene-environment interaction and the intergenerational transmission of parenting: Testing the differential-susceptibility hypothesis. Psychiatric Quarterly, 83, 2940.CrossRefGoogle ScholarPubMed
Beaver, K. M., Wright, J. P., & Walsh, A. (2008). A gene-based evolutionary explanation for the association between criminal involvement and number of sex partners. Biodemography and Social Biology, 54, 4755.CrossRefGoogle ScholarPubMed
Belsky, J., & Beaver, K. M. (2011). Cumulative-genetic plasticity, parenting and adolescent self-regulation. Journal of Child Psychology and Psychiatry, 52, 619626.CrossRefGoogle ScholarPubMed
Belsky, J., & Pluess, M. (2009). Beyond diathesis stress: Differential susceptibility to environmental influences. Psychological Bulletin, 135, 885.CrossRefGoogle ScholarPubMed
Bond, L., Butler, H., Thomas, L., Carlin, J., Glover, S., Bowes, G., & Patton, G. (2007). Social and school connectedness in early secondary school as predictors of late teenage substance use, mental health, and academic outcomes. Journal of Adolescent Health, 40, 357.CrossRefGoogle ScholarPubMed
Bongers, I. L., Koot, H. M., van der Ende, J., & Verhulst, F. C. (2003). The normative development of child and adolescent problem behavior. Journal of Abnormal Psychology, 112, 179.CrossRefGoogle ScholarPubMed
Borowsky, I. W., Ireland, M., & Resnick, M. D. (2001). Adolescent suicide attempts: Risks and protectors. Pediatrics, 107, 485493.CrossRefGoogle ScholarPubMed
Brody, G. H., Chen, Y. F., & Beach, S. R. (2013). Differential susceptibility to prevention: GABAergic, dopaminergic, and multilocus effects. Journal of Child Psychology and Psychiatry, 54, 863871.CrossRefGoogle ScholarPubMed
Brody, G. H., Yu, T., Beach, S. R., & Philibert, R. A. (2015). Prevention effects ameliorate the prospective association between nonsupportive parenting and diminished telomere length. Prevention Science, 16, 171180.CrossRefGoogle ScholarPubMed
Buchmann, A. F., Zohsel, K., Blomeyer, D., Hohm, E., Hohmann, S., Jennen-Steinmetz, C., … Esser, G. (2014). Interaction between prenatal stress and dopamine D4 receptor genotype in predicting aggression and cortisol levels in young adults. Psychopharmacology, 231, 30893097.CrossRefGoogle ScholarPubMed
Buckholtz, J. W., & Meyer-Lindenberg, A. (2008). MAOA and the neurogenetic architecture of human aggression. Trends in Neurosciences, 31, 120129.CrossRefGoogle ScholarPubMed
Buckholtz, J. W., Treadway, M. T., Cowan, R. L., Woodward, N. D., Li, R., Ansari, M. S., … Kessler, R. M. (2010). Dopaminergic network differences in human impulsivity. Science, 329, 532.CrossRefGoogle ScholarPubMed
Burt, S. A. (2009). Are there meaningful etiological differences within antisocial behavior? Results of a meta-analysis. Clinical Psychology Review, 29, 163178.CrossRefGoogle ScholarPubMed
Burt, S. A. (2013). Do etiological influences on aggression overlap with those on rule breaking? A meta-analysis. Psychological Medicine, 43, 18011812.CrossRefGoogle ScholarPubMed
Burt, S. A., Klump, K. L., Gorman-Smith, D., & Neiderhiser, J. M. (2016). Neighborhood disadvantage alters the origins of children's nonaggressive conduct problems. Clinical Psychological Science, 4, 511526.CrossRefGoogle ScholarPubMed
Burt, S. A., & Mikolajewski, A. J. (2008). Preliminary evidence that specific candidate genes are associated with adolescent-onset antisocial behavior. Aggressive Behavior, 34, 437445.CrossRefGoogle ScholarPubMed
Byrd, A. L., & Manuck, S. B. (2014). MAOA, childhood maltreatment, and antisocial behavior: Meta-analysis of a gene-environment interaction. Biological Psychiatry, 75, 917.CrossRefGoogle ScholarPubMed
Calvete, E., & Cardeñoso, O. (2005). Gender differences in cognitive vulnerability to depression and behavior problems in adolescents. Journal of Abnormal Child Psychology, 33, 179192.CrossRefGoogle ScholarPubMed
Carrel, L., & Willard, H. F. (2005). X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature, 434, 400404.CrossRefGoogle ScholarPubMed
Casey, B. J., Jones, R. M., & Hare, T. A. (2008). The adolescent brain. Annals of the New York Academy of Sciences, 1124, 111126.CrossRefGoogle ScholarPubMed
Caspi, A., Hariri, A. R., Holmes, A., Uher, R., & Moffitt, T. E. (2010). Genetic sensitivity to the environment: The case of the serotonin transporter gene and its implications for studying complex diseases and traits. Focus, 8, 398416.CrossRefGoogle Scholar
Caspi, A., Houts, R. M., Belsky, D. W., Goldman-Mellor, S. J., Harrington, H., Israel, S., … Moffitt, T. E. (2014). The p factor: One general psychopathology factor in the structure of psychiatric disorders? Clinical Psychological Science, 2, 119137.CrossRefGoogle Scholar
Caspi, A., McClay, J., Moffitt, T. E., Mill, J., Martin, J., Craig, I. W., … Poulton, R. (2002). Role of genotype in the cycle of violence in maltreated children. Science, 297, 851854.CrossRefGoogle ScholarPubMed
Castellani, V., Pastorelli, C., Eisenberg, N., Gerbino, M., Di Giunta, L., Ceravolo, R., & Milioni, M. (2014). Hostile, aggressive family conflict trajectories during the transition to adulthood: Associations with adolescent Big Five and emerging adulthood adjustment problems. Journal of Adolescence, 37, 647658.CrossRefGoogle ScholarPubMed
Catalano, R. F., & Hawkins, J. D. (1996). A theory of antisocial behavior. In Hawkins, J. D. (Ed.), Delinquency and crime: Current theories. Cambridge: Cambridge University Press.Google Scholar
Cleveland, H. H., Wiebe, R. P., & Rowe, D. C. (2005). Sources of exposure to smoking and drinking friends among adolescents: A behavioral genetic evaluation. Journal of Genetic Psychology, 166, 153169.Google ScholarPubMed
Cordell, H. J. (2009). Detecting gene–gene interactions that underlie human diseases. Nature Reviews Genetics, 10, 392404.CrossRefGoogle ScholarPubMed
Côté, S., Zoccolillo, M., Tremblay, R. E., Nagin, D., & Vitaro, F. (2001). Predicting girls’ conduct disorder in adolescence from childhood trajectories of disruptive behaviors. Journal of the American Academy of Child & Adolescent Psychiatry, 40, 678684.CrossRefGoogle ScholarPubMed
Crosnoe, R., Erickson, K. G., & Dornbusch, S. M. (2002). Protective functions of family relationships and school factors on the deviant behavior of adolescent boys and girls reducing the impact of risky friendships. Youth & Society, 33, 515544.CrossRefGoogle Scholar
Cumming, P., Brown, E., Damsma, G., & Fibiger, H. (1992). Formation and clearance of interstitial metabolites of dopamine and serotonin in the rat striatum: An in vivo microdialysis study. Journal of Neurochemistry, 59, 19051914.CrossRefGoogle ScholarPubMed
Davidson, R. J., Putnam, K. M., & Larson, C. L. (2000). Dysfunction in the neural circuitry of emotion regulation—A possible prelude to violence. Science, 289, 591594.CrossRefGoogle ScholarPubMed
Derringer, J., Krueger, R. F., Dick, D. M., Saccone, S., Grucza, R. A., Agrawal, A., … Nurnberger, J. I. (2010). Predicting sensation seeking from dopamine genes a candidate-system approach. Psychological Science, 21, 12821290.CrossRefGoogle Scholar
DeYoung, C. G., Getchell, M., Koposov, R. A., Yrigollen, C. M., Haeffel, G. J., Klinteberg, B., … Grigorenko, E. L. (2010). Variation in the catechol-O-methyltransferase Val158Met polymorphism associated with conduct disorder and ADHD symptoms among adolescent male delinquents. Psychiatric Genetics, 20, 20.CrossRefGoogle Scholar
Dick, D. M. (2011). Gene-environment interaction in psychological traits and disorders. Annual Review of Clinical Psychology, 7, 383.CrossRefGoogle ScholarPubMed
DiLalla, L. F., & Gottesman, I. I. (1989). Heterogeneity of causes for delinquency and criminality: Lifespan perspectives. Development and Psychopathology, 1, 339349.CrossRefGoogle Scholar
Doyle, A. B., & Markiewicz, D. (2005). Parenting, marital conflict and adjustment from early- to mid-adolescence: Mediated by adolescent attachment style? Journal of Youth and Adolescence, 34, 97110.CrossRefGoogle Scholar
D'Souza, U. M., & Craig, I. W. (2008). Functional genetic polymorphisms in serotonin and dopamine gene systems and their significance in behavioural disorders. Progress in Brain Research, 172, 7398.CrossRefGoogle ScholarPubMed
D'Unger, A. V., Land, K. C., & McCall, P. L. (2002). Sex differences in age patterns of delinquent/criminal careers: Results from Poisson latent class analyses of the Philadelphia cohort study. Journal of Quantitative Criminology, 18, 349375.CrossRefGoogle Scholar
Edwards, A. C., Bacanu, S. A., Bigdeli, T. B., Moscati, A., & Kendler, K. S. (2016). Evaluating the dopamine hypothesis of schizophrenia in a large-scale genome-wide association study. Schizophrenia Research. Advance online publication.CrossRefGoogle Scholar
Everitt, B. J., & Robbins, T. W. (2000). Second-order schedules of drug reinforcement in rats and monkeys: Measurement of reinforcing efficacy and drug-seeking behaviour. Psychopharmacology, 153, 1730.CrossRefGoogle ScholarPubMed
Farrington, D. P. (1987). Early precursors of frequent offending. In Families, schools, and delinquency prevention (pp. 2750). New York: Springer.CrossRefGoogle Scholar
Farrington, D. P. (1995). The development of offending and antisocial behaviour from childhood: Key findings from the Cambridge Study in Delinquent Development. Journal of Child Psychology and Psychiatry, 6, 929964.CrossRefGoogle Scholar
Farrington, D. P. (1998). Predictors, causes, and correlates of male youth violence. Crime and Justice, 24, 421475.CrossRefGoogle Scholar
Ferguson, C. J. (2010). Genetic contributions to antisocial personality and behavior: A meta-analytic review from an evolutionary perspective. Journal of Social Psychology, 150, 160180.CrossRefGoogle ScholarPubMed
Ficks, C. A., & Waldman, I. D. (2014). Candidate genes for aggression and antisocial behavior: A meta-analysis of association studies of the 5HTTLPR and MAOA-uVNTR. Behavior Genetics, 44, 427444.CrossRefGoogle ScholarPubMed
Fontaine, N., Carbonneau, R., Vitaro, F., Barker, E. D., & Tremblay, R. E. (2009). Research review: A critical review of studies on the developmental trajectories of antisocial behavior in females. Journal of Child Psychology and Psychiatry, 50, 363385.CrossRefGoogle ScholarPubMed
Foster, E. M., & Jones, D. E. (2005). The high costs of aggression: Public expenditures resulting from conduct disorder. American Journal of Public Health, 95, 17671772.CrossRefGoogle ScholarPubMed
Frick, P. J., Lahey, B. B., Loeber, R., Tannenbaum, L., Van Horn, Y., Christ, M. A. G., … Hanson, K. (1993). Oppositional defiant disorder and conduct disorder: A meta-analytic review of factor analyses and cross-validation in a clinic sample. Clinical Psychology Review, 13, 319340.CrossRefGoogle Scholar
Gelernter, J., Kranzler, H. R., Sherva, R., Almasy, L., Koesterer, R., Smith, A. H., … Wodarz, N. (2014). Genome-wide association study of alcohol dependence: Significant findings in African- and European-Americans including novel risk loci. Molecular Psychiatry, 19, 4149.CrossRefGoogle ScholarPubMed
Gibbons, J. L., Barr, G. A., Bridger, W. H., & Liebowitz, S. F. (1979). Manipulations of dietary tryptophan: Effects on mouse killing and brain serotonin in the rat. Brain Research, 169, 139153.CrossRefGoogle ScholarPubMed
Gizer, I. R., & Waldman, I. D. (2012). Double dissociation between lab measures of inattention and impulsivity and the dopamine transporter gene (DAT1) and dopamine D4 receptor gene (DRD4). Journal of Abnormal Psychology, 121, 1011.CrossRefGoogle ScholarPubMed
Guo, G., Roettger, M. E., & Shih, J. C. (2007). Contributions of the DAT1 and DRD2 genes to serious and violent delinquency among adolescents and young adults. Human Genetics, 121, 125136.CrossRefGoogle ScholarPubMed
Haberstick, B. C., Smolen, A., & Hewitt, J. K. (2006). Family-based association test of the 5HTTLPR and aggressive behavior in a general population sample of children. Biological Psychiatry, 59, 836843.CrossRefGoogle Scholar
Hackman, D. A., Farah, M. J., & Meaney, M. J. (2010). Socioeconomic status and the brain: Mechanistic insights from human and animal research. Nature Reviews Neuroscience, 11, 651659.CrossRefGoogle ScholarPubMed
Harden, K. P., Hill, J. E., Turkheimer, E., & Emery, R. E. (2008). Gene-environment correlation and interaction in peer effects on adolescent alcohol and tobacco use. Behavior Genetics, 38, 339347.CrossRefGoogle ScholarPubMed
Heils, A., Teufel, A., Petri, S., Stöber, G., Riederer, P., Bengel, D., & Lesch, K. P. (1996). Allelic variation of human serotonin transporter gene expression. Journal of Neurochemistry, 66, 26212624.CrossRefGoogle ScholarPubMed
Hohmann, S., Becker, K., Fellinger, J., Banaschewski, T., Schmidt, M. H., Esser, G., & Laucht, M. (2009). Evidence for epistasis between the 5-HTTLPR and the dopamine D4 receptor polymorphisms in externalizing behavior among 15-year-olds. Journal of Neural Transmission, 116, 16211629.CrossRefGoogle ScholarPubMed
Holmes, A., Murphy, D. L., & Crawley, J. N. (2002). Reduced aggression in mice lacking the serotonin transporter. Psychopharmacology, 161, 160167.CrossRefGoogle ScholarPubMed
Hu, X. Z., Lipsky, R. H., Zhu, G., Akhtar, L. A., Taubman, J., Greenberg, B. D., … Murphy, D. L. (2006). Serotonin transporter promoter gain-of-function genotypes are linked to obsessive-compulsive disorder. American Journal of Human Genetics, 78, 815826.CrossRefGoogle ScholarPubMed
Huizinga, D., & Elliott, D. S. (1986). Reassessing the reliability and validity of self-report delinquency measures. Journal of Quantitative Criminology, 2, 293327.CrossRefGoogle Scholar
Hutchison, K. E., Stallings, M., McGeary, J., & Bryan, A. (2004). Population stratification in the candidate gene study: Fatal threat or red herring? Psychological Bulletin, 130, 66.CrossRefGoogle ScholarPubMed
Ingoldsby, E. M., & Shaw, D. S. (2002). Neighborhood contextual factors and early-starting antisocial pathways. Clinical Child and Family Psychology Review, 5, 2155.CrossRefGoogle ScholarPubMed
International Schizophrenia Consortium. (2009). Common polygenic variation contributes to risk of schizophrenia that overlaps with bipolar disorder. Nature, 460, 748752.CrossRefGoogle Scholar
Jacobson, K. C., & Rowe, D. C. (1999). Genetic and environmental influences on the relationships between family connectedness, school connectedness, and adolescent depressed mood: Sex differences. Developmental Psychology, 35, 926.CrossRefGoogle ScholarPubMed
Jaffee, S. R., Caspi, A., Moffitt, T. E., Dodge, K. A., Rutter, M., Taylor, A., & Tully, L. A. (2005). Nature × Nurture: Genetic vulnerabilities interact with physical maltreatment to promote conduct problems. Development and Psychopathology, 17, 6784.CrossRefGoogle ScholarPubMed
Jaffee, S. R., Caspi, A., Moffitt, T. E., & Taylor, A. (2004). Physical maltreatment victim to antisocial child: Evidence of an environmentally mediated process. Journal of Abnormal Psychology, 113, 44.CrossRefGoogle ScholarPubMed
Jaffee, S. R., & Price, T. S. (2007). Gene–environment correlations: A review of the evidence and implications for prevention of mental illness. Molecular Psychiatry, 12, 432442.CrossRefGoogle ScholarPubMed
Janssens, A., Goossens, L., Van Den Noortgate, W., Colpin, H., Verschueren, K., & Van Leeuwen, K. (2015). Parents’ and adolescents’ perspectives on parenting evaluating conceptual structure, measurement invariance, and criterion validity. Assessment, 22, 473489.CrossRefGoogle ScholarPubMed
Jönsson, E. G., Nöthen, M. M., Grünhage, F., Farde, L., Nakashima, Y., Propping, P., & Sedvall, G. C. (1999). Polymorphisms in the dopamine D2 receptor gene and their relationships to striatal dopamine receptor density of healthy volunteers. Molecular Psychiatry, 4, 290296.CrossRefGoogle ScholarPubMed
Justice Policy Institute. (2014). Sticker shock: Calculating the full price tag for youth incarceration. Washington, DC: Author.Google Scholar
Keller, M. C. (2014). Gene × Environment interaction studies have not properly controlled for potential confounders: The problem and the (simple) solution. Biological Psychiatry, 75, 1824.CrossRefGoogle Scholar
Kelley, A. E., Schochet, T., & Landry, C. F. (2004). Risk taking and novelty seeking in adolescence: Introduction to part I. Annals of the New York Academy of Sciences, 1021, 2732.CrossRefGoogle ScholarPubMed
Kendler, K. S., Jacobson, K. C., Gardner, C. O., Gillespie, N., Aggen, S. A., & Prescott, C. A. (2007). Creating a social world: A developmental twin study of peer-group deviance. Archives of General Psychiatry, 64, 958965.CrossRefGoogle ScholarPubMed
Kendler, K. S., Jacobson, K., Myers, J. M., & Eaves, L. J. (2008). A genetically informative developmental study of the relationship between conduct disorder and peer deviance in males. Psychological Medicine, 38, 10011011.CrossRefGoogle ScholarPubMed
Kendler, K. S., Prescott, C. A., Myers, J., & Neale, M. C. (2003). The structure of genetic and environmental risk factors for common psychiatric and substance use disorders in men and women. Archives of General Psychiatry, 60, 929937.CrossRefGoogle ScholarPubMed
Kim-Cohen, J., Caspi, A., Taylor, A., Williams, B., Newcombe, R., Craig, I. W., & Moffitt, T. E. (2006). MAOA, maltreatment, and gene–environment interaction predicting children's mental health: New evidence and a meta-analysis. Molecular Psychiatry, 11, 903913.CrossRefGoogle ScholarPubMed
Knafo, A., & Jaffee, S. R. (2013). Gene–environment correlation in developmental psychopathology. Development and Psychopathology, 25, 16.CrossRefGoogle ScholarPubMed
Kohavi, R. (1995). A study of cross-validation and bootstrap for accuracy estimation and model selection. In Proceedings of the 14th International Joint Conference on Artificial Engineering (pp. 11371145). San Francisco, CA: Morgan Kaufmann.Google Scholar
Kok, B. E., & Fredrickson, B. L. (2010). Upward spirals of the heart: Autonomic flexibility, as indexed by vagal tone, reciprocally and prospectively predicts positive emotions and social connectedness. Biological Psychology, 85, 432436.CrossRefGoogle ScholarPubMed
Krueger, R. F., Markon, K. E., Patrick, C. J., Benning, S. D., & Kramer, M. D. (2007). Linking antisocial behavior, substance use, and personality: An integrative quantitative model of the adult externalizing spectrum. Journal of Abnormal Psychology, 116, 645.CrossRefGoogle ScholarPubMed
Kuperminc, G. P., & Allen, J. P. (2001). Social orientation: Problem behavior and motivations toward interpersonal problem solving among high risk adolescents. Journal of Youth and Adolescence, 30, 597622.CrossRefGoogle ScholarPubMed
Lahey, B. B., Van Hulle, C. A., Singh, A. L., Waldman, I. D., & Rathouz, P. J. (2011). Higher-order genetic and environmental structure of prevalent forms of child and adolescent psychopathology. Archives of General Psychiatry, 68, 181189.CrossRefGoogle ScholarPubMed
Latendresse, S. J., Bates, J. E., Goodnight, J. A., Lansford, J. E., Budde, J. P., Goate, A., … Dick, D. M. (2011). Differential susceptibility to adolescent externalizing trajectories: Examining the interplay between CHRM2 and peer group antisocial behavior. Child Development, 82, 17971814.CrossRefGoogle ScholarPubMed
Lee, S. S. (2011). Deviant peer affiliation and antisocial behavior: Interaction with monoamine oxidase A (MAOA) genotype. Journal of Abnormal Child Psychology, 39, 321332.CrossRefGoogle ScholarPubMed
Li, J. J., Berk, M. S., & Lee, S. S. (2013). Differential susceptibility in longitudinal models of gene–environment interaction for adolescent depression. Development and Psychopathology, 25, 9911003.CrossRefGoogle ScholarPubMed
Liao, D. L., Hong, C. J., Shih, H. L., & Tsai, S. J. (2004). Possible association between serotonin transporter promoter region polymorphism and extremely violent crime in Chinese males. Neuropsychobiology, 50, 284287.CrossRefGoogle ScholarPubMed
Little, T. D. (2013). Longitudinal structural equation modeling. New York: Guilford Press.Google Scholar
Lo, Y., Mendell, N. R., & Rubin, D. B. (2001). Testing the number of components in a normal mixture. Biometrika, 88, 767778.CrossRefGoogle Scholar
Loeber, R., & Schmaling, K. B. (1985). Empirical evidence for overt and covert patterns of antisocial conduct problems: A meta-analysis. Journal of Abnormal Child Psychology, 13, 337353.CrossRefGoogle Scholar
Marceau, K., Palmer, R. H., Neiderhiser, J. M., Smith, T. F., McGeary, J. E., & Knopik, V. S. (2016). Passive rGE or developmental gene-environment cascade? An investigation of the role of xenobiotic metabolism genes in the association between smoke exposure during pregnancy and child birth weight. Behavior Genetics, 46, 365377.CrossRefGoogle ScholarPubMed
Massoglia, M. (2006). Desistance or displacement? The changing patterns of offending from adolescence to young adulthood. Journal of Quantitative Criminology, 22, 215239.CrossRefGoogle Scholar
Masten, C. L., Eisenberger, N. I., Borofsky, L. A., Pfeifer, J. H., McNealy, K., Mazziotta, J. C., & Dapretto, M. (2009). Neural correlates of social exclusion during adolescence: Understanding the distress of peer rejection. Social Cognitive and Affective Neuroscience, 4, 143157.CrossRefGoogle ScholarPubMed
Mazerolle, P., Brame, R., Paternoster, R., Piquero, A., & Dean, C. (2000). Onset age, persistence, and offending versatility: Comparisons across gender. Criminology, 38, 11431172.CrossRefGoogle Scholar
McNeely, C. A., Nonnemaker, J. M., & Blum, R. W. (2002). Promoting school connectedness: Evidence from the national longitudinal study of adolescent health. Journal of School Health, 72, 138146.CrossRefGoogle ScholarPubMed
Mercer, N., Keijsers, L., Crocetti, E., Branje, S., & Meeus, W. (2016). Adolescent abstention from delinquency: Examining the mediating role of time spent with (delinquent) peers. Journal of Research on Adolescence. Advance online publication.CrossRefGoogle ScholarPubMed
Moffitt, T. E. (1993). Adolescence-limited and life-course-persistent antisocial behavior: A developmental taxonomy. Psychological Review, 100, 674.CrossRefGoogle ScholarPubMed
Moffitt, T. E. (2001). Sex differences in antisocial behaviour: Conduct disorder, delinquency, and violence in the Dunedin Longitudinal Study. New York: Cambridge University Press.CrossRefGoogle Scholar
Moffitt, T. E. (2003). Life course persistent and adolescence-limited antisocial behavior: A 10-year research review and a research agenda. In Lahey, B. B., Moffitt, T. E., & Caspi, A. (Eds.), Causes of conduct disorder and juvenile delinquency (pp. 4975). New York: Guilford Press.Google Scholar
Moffitt, T. E., Caspi, A., Harrington, H., & Milne, B. J. (2002). Males on the life-course-persistent and adolescence-limited antisocial pathways: Follow-up at age 26 years. Development and Psychopathology, 14, 179207.CrossRefGoogle ScholarPubMed
Moore, T. M., Scarpa, A., & Raine, A. (2002). A meta-analysis of serotonin metabolite 5-HIAA and antisocial behavior. Aggressive Behavior, 28, 299316.CrossRefGoogle Scholar
Mrug, S., & Windle, M. (2009). Bidirectional influences of violence exposure and adjustment in early adolescence: Externalizing behaviors and school connectedness. Journal of Abnormal Child Psychology, 37, 611623.CrossRefGoogle ScholarPubMed
Muthén, L. K., & Muthén, B. O. (2015). Mplus user's guide (7th ed.). Los Angeles, CA: Author.Google Scholar
Nagin, D. S., & Tremblay, R. E. (2005). Developmental trajectory groups: Fact or a useful statistical fiction? Criminology, 43, 873904.CrossRefGoogle Scholar
Nikolova, Y. S., Ferrell, R. E., Manuck, S. B., & Hariri, A. R. (2011). Multilocus genetic profile for dopamine signaling predicts ventral striatum reactivity. Neuropsychopharmacology, 36, 19401947.CrossRefGoogle ScholarPubMed
Nobile, M., Giorda, R., Marino, C., Carlet, O., Pastore, V., Vanzin, L., … Battaglia, M. (2007). Socioeconomic status mediates the genetic contribution of the dopamine receptor D4 and serotonin transporter linked promoter region repeat polymorphisms to externalization in preadolescence. Development and Psychopathology, 19, 11471160.CrossRefGoogle ScholarPubMed
Noble, E. P., Gottschalk, L. A., Fallon, J. H., Ritchie, T. L., & Wu, J. C. (1997). D 2 dopamine receptor polymorphism and brain regional glucose metabolism. American Journal of Medical Genetics, 74, 162166.3.0.CO;2-W>CrossRefGoogle Scholar
Nylund, K. L., Asparouhov, T., & Muthén, B. O. (2007). Deciding on the number of classes in latent class analysis and growth mixture modeling: A Monte Carlo simulation study. Structural Equation Modeling, 14, 535569.CrossRefGoogle Scholar
Odgers, C. L., Moffitt, T. E., Broadbent, J. M., Dickson, N., Hancox, R. J., Harrington, H., … Caspi, A. (2008). Female and male antisocial trajectories: From childhood origins to adult outcomes. Development and Psychopathology, 20, 673716.CrossRefGoogle ScholarPubMed
Pappa, I., St. Pourcain, B., Benke, K., Cavadino, A., Hakulinen, C., Nivard, M. G., … Evans, D. M. (2016). A genome-wide approach to children's aggressive behavior: The EAGLE consortium. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics, 171, 562572.CrossRefGoogle ScholarPubMed
Pereira, F., Mitchell, T., & Botvinick, M. (2009). Machine learning classifiers and fMRI: A tutorial overview. NeuroImage, 45, S199S209.CrossRefGoogle ScholarPubMed
Pingault, J. B., Viding, E., Galéra, C., Greven, C. U., Zheng, Y., Plomin, R., & Rijsdijk, F. (2015). Genetic and environmental influences on the developmental course of attention-deficit/hyperactivity disorder symptoms from childhood to adolescence. JAMA Psychiatry, 72, 651658.CrossRefGoogle ScholarPubMed
Plomin, R., Haworth, C. M., & Davis, O. S. (2009). Common disorders are quantitative traits. Nature Reviews Genetics, 10, 872878.CrossRefGoogle ScholarPubMed
Pohjalainen, T., Rinne, J. O., Någren, K., Lehikoinen, P., Anttila, K., Syvälahti, E. K. G., & Hietala, J. (1998). The A1 allele of the human D~ 2 dopamine receptor gene predicts low D~ 2 receptor availability in healthy volunteers. Molecular Psychiatry, 3, 256260.CrossRefGoogle Scholar
Qian, Q. J., Liu, J., Wang, Y. F., Yang, L., Guan, L. L., & Faraone, S. V. (2009). Attention deficit hyperactivity disorder comorbid oppositional defiant disorder and its predominately inattentive type: Evidence for an association with COMT but not MAOA in a Chinese sample. Behavioral and Brain Functions, 5, 1.CrossRefGoogle ScholarPubMed
Resnick, M. D., Bearman, P. S., Blum, R. W., Bauman, K. E., Harris, K. M., Jones, J., … Ireland, M. (1997). Protecting adolescents from harm: Findings from the National Longitudinal Study on Adolescent Health. Journal of the American Medical Association, 278, 823832.CrossRefGoogle ScholarPubMed
Rhee, S. H., & Waldman, I. D. (2002). Genetic and environmental influences on antisocial behavior: A meta-analysis of twin and adoption studies. Psychological Bulletin, 128, 490.CrossRefGoogle ScholarPubMed
Sabol, S. Z., Hu, S., & Hamer, D. (1998). A functional polymorphism in the monoamine oxidase A gene promoter. Human Genetics, 103, 273279.CrossRefGoogle ScholarPubMed
Sakai, Y., Nishikawa, M., Leyton, M., Benkelfat, C., Young, S. N., & Diksic, M. (2006). Cortical trapping of α-[11 C] methyl-L-tryptophan, an index of serotonin synthesis, is lower in females than males. NeuroImage, 33, 815824.CrossRefGoogle ScholarPubMed
Salvatore, J. E., & Dick, D. M. (2016). Genetic influences on conduct disorder. Neuroscience & Biobehavioral Reviews. Advance online publication.Google Scholar
Samek, D. R., Hicks, B. M., Keyes, M. A., Iacono, W. G., & McGue, M. (2016). Antisocial peer affiliation and externalizing disorders: Evidence for Gene × Environment × Development interaction. Development and Psychopathology. Advance online publication.Google Scholar
Sampson, R. J., & Laub, J. H. (2003). Life-course desisters? Trajectories of crime among delinquent boys followed to age 70. Criminology, 41, 555592.CrossRefGoogle Scholar
Schlüter, T., Winz, O., Henkel, K., Prinz, S., Rademacher, L., Schmaljohann, J., … Mottaghy, F. M. (2013). The impact of dopamine on aggression: An [18F]-FDOPA PET Study in healthy males. Journal of Neuroscience, 33, 1688916896.CrossRefGoogle ScholarPubMed
Seo, D., Patrick, C. J., & Kennealy, P. J. (2008). Role of serotonin and dopamine system interactions in the neurobiology of impulsive aggression and its comorbidity with other clinical disorders. Aggression and Violent Behavior, 13, 383395.CrossRefGoogle ScholarPubMed
Shochet, I. M., Dadds, M. R., Ham, D., & Montague, R. (2006). School connectedness is an underemphasized parameter in adolescent mental health: Results of a community prediction study. Journal of Clinical Child and Adolescent Psychology, 35, 170179.CrossRefGoogle ScholarPubMed
Simonson, M. A., Wills, A. G., Keller, M. C., & McQueen, M. B. (2011). Recent methods for polygenic analysis of genome-wide data implicate an important effect of common variants on cardiovascular disease risk. BMC Medical Genetics, 12, 146.CrossRefGoogle ScholarPubMed
Smolen, A., Whitsel, E. A., Tabor, J., Killeya-Jones, L. A., Cuthbertson, C. C., Hussey, J. M., … Harris, K. M. (2013). Add Health documentation: Candidate genes. Retrieved from http://www.cpc.unc.edu/projects/addhealth/documentation/guides/DNA_documentation.pdf Google Scholar
Soderstrom, H., Blennow, K., Manhem, A., & Forsman, A. (2001). CSF studies in violent offenders: 5-HIAA as a negative and HVA as a positive predictor of psychopathy. Journal of Neural Transmission, 108, 869878.CrossRefGoogle ScholarPubMed
Stice, E., Yokum, S., Burger, K., Epstein, L., & Smolen, A. (2012). Multilocus genetic composite reflecting dopamine signaling capacity predicts reward circuitry responsivity. Journal of Neuroscience, 32, 1009310100.CrossRefGoogle ScholarPubMed
Tackett, J. L., Daoud, S. L., De Bolle, M., & Burt, S. A. (2013). Is relational aggression part of the externalizing spectrum? A bifactor model of youth antisocial behavior. Aggressive Behavior, 39, 149159.CrossRefGoogle ScholarPubMed
Taylor, J., Iacono, W. G., & McGue, M. (2000). Evidence for a genetic etiology of early-onset delinquency. Journal of Abnormal Psychology, 109, 634643.CrossRefGoogle ScholarPubMed
Thibodeau, E. L., Cicchetti, D., & Rogosch, F. A. (2015). Child maltreatment, impulsivity, and antisocial behavior in African American children: Moderation effects from a cumulative dopaminergic gene index. Development and Psychopathology, 27, 16211636.CrossRefGoogle ScholarPubMed
Thompson, R., & Tabone, J. K. (2010). The impact of early alleged maltreatment on behavioral trajectories. Child Abuse & Neglect, 34, 907916.CrossRefGoogle ScholarPubMed
Tidey, J. W., & Miczek, K. A. (1996). Social defeat stress selectively alters mesocorticolimbic dopamine release: An in vivo microdialysis study. Brain Research, 721, 140149.CrossRefGoogle ScholarPubMed
Tielbeek, J. J., Karlsson Linnér, R., Beers, K., Posthuma, D., Popma, A., & Polderman, T. J. (2016). Meta-analysis of the serotonin transporter promoter variant (5-HTTLPR) in relation to adverse environment and antisocial behavior. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics, 171, 748760.CrossRefGoogle ScholarPubMed
Tielbeek, J. J., Medland, S. E., Benyamin, B., Byrne, E. M., Heath, A. C., Madden, P. A., … Verweij, K. J. (2012). Unraveling the genetic etiology of adult antisocial behavior: A genome-wide association study. PLOS ONE, 7, e45086.CrossRefGoogle ScholarPubMed
Tomer, R., Slagter, H. A., Christian, B. T., Fox, A. S., King, C. R., Murali, D., … Davidson, R. J. (2014). Love to win or hate to lose? Asymmetry of dopamine D2 receptor binding predicts sensitivity to reward versus punishment. Journal of Cognitive Neuroscience, 26, 10391048.CrossRefGoogle ScholarPubMed
Tripp, G., & Wickens, J. R. (2008). Research review: Dopamine transfer deficit: A neurobiological theory of altered reinforcement mechanisms in ADHD. Journal of Child Psychology and Psychiatry, 49, 691704.CrossRefGoogle ScholarPubMed
Tung, I., & Lee, S. S. (2016). Latent trajectories of adolescent antisocial behavior: Serotonin transporter linked polymorphic region (5-HTTLPR) genotype influences sensitivity to perceived parental support. Development and Psychopathology. Advance online publication.Google Scholar
Tuvblad, C., Grann, M., & Lichtenstein, P. (2006). Heritability for adolescent antisocial behavior differs with socioeconomic status: Gene–environment interaction. Journal of Child Psychology and Psychiatry, 47, 734743.CrossRefGoogle ScholarPubMed
Van Erp, A. M., & Miczek, K. A. (2000). Aggressive behavior, increased accumbal dopamine, and decreased cortical serotonin in rats. Journal of Neuroscience, 20, 93209325.CrossRefGoogle ScholarPubMed
Van Hulle, C. A., Rodgers, J. L., D'Onofrio, B. M., Waldman, I. D., & Lahey, B. B. (2007). Sex differences in the causes of self-reported adolescent delinquency. Journal of Abnormal Psychology, 116, 236.CrossRefGoogle ScholarPubMed
van IJzendoorn, M. H., Belsky, J., & Bakermans-Kranenburg, M. J. (2012). Serotonin transporter genotype 5HTTLPR as a marker of differential susceptibility? A meta-analysis of child and adolescent gene-by-environment studies. Translational Psychiatry, 2, e147.CrossRefGoogle ScholarPubMed
VanNess, S. H., Owens, M. J., & Kilts, C. D. (2005). The variable number of tandem repeats element in DAT1 regulates in vitro dopamine transporter density. BMC Genetics, 6, 1.CrossRefGoogle ScholarPubMed
Van Tol, H. H., Wu, C. M., Guan, H. C., Ohara, K., Bunzow, J. R., Civelli, O., … Jovanovic, V. (1992). Multiple dopamine D4 receptor variants in the human population. Nature, 358, 149152.Google ScholarPubMed
Veroude, K., Zhang-James, Y., Fernàndez-Castillo, N., Bakker, M. J., Cormand, B., & Faraone, S. V. (2016). Genetics of aggressive behavior: An overview. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics, 171, 343.CrossRefGoogle Scholar
Walker, S. O., & Plomin, R. (2006). Nature, nurture, and perceptions of the classroom environment as they relate to teacher-assessed academic achievement: A twin study of nine-year-olds. Educational Psychology, 26, 541561.CrossRefGoogle Scholar
Walters, G. D., & Ruscio, J. (2013). Trajectories of youthful antisocial behavior: Categories or continua? Journal of Abnormal Child Psychology, 41, 653666.CrossRefGoogle ScholarPubMed
Whittle, S., Simmons, J. G., Dennison, M., Vijayakumar, N., Schwartz, O., Yap, M. B., … Allen, N. B. (2014). Positive parenting predicts the development of adolescent brain structure: A longitudinal study. Developmental Cognitive Neuroscience, 8, 717.CrossRefGoogle ScholarPubMed
Wichers, M., Wigman, J. T. W., & Myin-Germeys, I. (2015). Micro-level affect dynamics in psychopathology viewed from complex dynamical system theory. Emotion Review, 7, 362367.CrossRefGoogle Scholar
Widom, C. S., & Brzustowicz, L. M. (2006). MAOA and the “cycle of violence”: Childhood abuse and neglect, MAOA genotype, and risk for violent and antisocial behavior. Biological Psychiatry, 60, 684689.CrossRefGoogle ScholarPubMed
Windhorst, D. A., Mileva-Seitz, V. R., Rippe, R. C., Tiemeier, H., Jaddoe, V. W., Verhulst, F. C., … Bakermans-Kranenburg, M. J. (2016). Beyond main effects of gene-sets: Harsh parenting moderates the association between a dopamine gene-set and child externalizing behavior. Brain and Behavior. Advance online publication.CrossRefGoogle Scholar
Wray, N. R., Yang, J., Hayes, B. J., Price, A. L., Goddard, M. E., & Visscher, P. M. (2013). Pitfalls of predicting complex traits from SNPs. Nature Reviews Genetics, 14, 507515.CrossRefGoogle ScholarPubMed
Zheng, Y., & Cleveland, H. H. (2013). Identifying gender-specific developmental trajectories of nonviolent and violent delinquency from adolescence to young adulthood. Journal of Adolescence, 36, 371381.CrossRefGoogle ScholarPubMed
Zheng, Y., & Cleveland, H. H. (2015). Differential genetic and environmental influences on developmental trajectories of antisocial behavior from adolescence to young adulthood. Journal of Adolescence, 45, 204213.CrossRefGoogle ScholarPubMed
Zubieta, J. K., Heitzeg, M. M., Smith, Y. R., Bueller, J. A., Xu, K., Xu, Y., … Goldman, D. (2003). COMT val158met genotype affects μ-opioid neurotransmitter responses to a pain stressor. Science, 299, 12401243.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Allele frequencies and Hardy–Weinberg equilibrium (HWE) tests

Figure 1

Table 2. Demographic information of the Caucasian-only genotypic sample by sex

Figure 2

Table 3. Principal component analysis of supportive parenting and school connectedness variables in Add Health (oblimin rotation)

Figure 3

Figure 1. Estimated means of the latent trajectories of antisocial behavior.

Figure 4

Table 4. Fit statistics for growth mixture model of ASB sum scores from age 13 to 32

Figure 5

Table 5. Regressions results from each k-fold cross-validation: Add Health multigenic risk scores predicting discrete probabilities of latent antisocial behavior class membership, supportive parenting, and school connectedness by sex

Figure 6

Table 6. Multinomial logistic regressions: Supportive parenting and multigenic risk scores (MRS)

Figure 7

Figure 2. School Connectedness × Add Health Multigenic Risk Score interaction predicting odds of adolescence-peaked pathway antisocial behavior membership (relative to low antisocial behavior) in males.

Figure 8

Table 7. Multinomial logistic regressions: School connectedness and multigenic risk scores (MRS)

Supplementary material: File

Li supplementary material 1

Li supplementary material

Download Li supplementary material 1(File)
File 37.5 KB