Decline in Genetic Influence on the Co-Occurrence of Alcohol, Marijuana, and Nicotine Dependence Symptoms From Age 14 to 29
Publication: American Journal of Psychiatry
Abstract
Objective
Cross-sectional studies have demonstrated high rates of comorbidity among substance use disorders. However, few studies have examined the developmental course of incident comorbidity and how it changes from adolescence to adulthood. The authors examine patterns of comorbidity among substance use disorders to gain insight into the effect of shared versus specific etiological influences on measures of substance abuse and dependence.
Method
The authors evaluated the pattern of correlations among nicotine, alcohol, and marijuana abuse and dependence symptom counts as well as their underlying genetic and environmental influences in a community-representative twin sample (N=3,762). Symptoms were assessed at ages 11, 14, 17, 20, 24, and 29 years. A single common factor was used to model the correlations among symptom counts at each age. The authors examined age-related changes in the influence of this general factor by testing for differences in the mean factor loading across time.
Results
Mean levels of abuse or dependence symptoms increased throughout adolescence, peaked around age 20, and declined from age 24 to age 29. The influence of the general factor was highest at ages 14 and 17, but decreased from age 17 to age 24. Genetic influences of the general factor declined considerably with age alongside an increase in nonshared environmental influences.
Conclusions
Adolescent substance abuse or dependence is largely a function of shared etiology. As young people age, their symptoms are increasingly influenced by substance-specific etiological factors. Heritability analyses revealed that the generalized risk is primarily influenced by genetic factors in adolescence, but nonshared environmental influences increase in importance as substance dependence becomes more specialized in adulthood.
Epidemiological studies have documented high rates of comorbidity among alcohol, nicotine, and illicit drug dependence disorders (1–5). Comorbidity suggests that these nosologically distinct disorders are caused, in part, by common etiological processes (6–8). In substance use disorders, the etiology shared among disorders indexed by comorbidity has been referred to as “externalizing” or “disinhibitory psychopathology” (9).
While comorbidity among substance use disorders has been extensively studied cross-sectionally—typically for lifetime prevalence rates—few studies have examined how comorbidity changes over time, an important topic given the large changes in incidence over the lifespan. Specifically, substance use tends to emerge in middle adolescence, increases substantially throughout adolescence, peaks in the early twenties, plateaus, and then decreases in the late twenties (10). Whether the rates of comorbidity among substance use disorders remain consistently high, however, is unknown. For example, rates of comorbidity might decline, suggesting that individuals begin to specialize in their substance use over time.
To investigate comorbidity across time, we used a large longitudinal twin study to examine changes in correlations among substance use disorder symptom counts over time. Assessments coincided with key developmental transitions in substance use, including before initiation (age 11), initiation (age 14), regular use (age 17), heavy use and dependence (ages 20 and 24), and the period when individuals decrease their use or exhibit patterns of persistent substance use problems (age 29). First, we examined patterns of mean-level change in nicotine, alcohol, and marijuana dependence symptoms. Second, we fitted a single common factor model to account for the correlations among incident symptom counts of nicotine, alcohol, and marijuana dependence at each age. This common factor model then allowed us to test for changes in the contributions of common and specific etiological influences on these symptom counts across time. Third, we used standard twin models to estimate genetic and environmental influences on the common factor over time.
In addition to analyzing the full sample, we also fitted separate models to a subsample of early-onset users. We did so because estimates of correlations at earlier ages may be affected by a minority of high-risk individuals who tend to exhibit dependence symptoms for multiple substances. Within this high-risk sample, correlations among dependence symptoms may remain high into adulthood. This could be obscured, however, in analyses of the full sample as substance use becomes more normative in adulthood. That is, if the correlations among substance use disorders are lower for the larger group of later-onset users, what is actually an artifact of early- instead of later-onset use would appear to be an overall decline in the correlations among substance use disorders. We thus performed our analyses both for the full sample and in a subsample of participants who had at least one symptom by their assessment at age 17, earlier than the age at which use of any of these substances becomes legal in the United States.
Method
Sample
Participants (N=3,762; 52% female) were drawn from the Minnesota Twin Family Study, a community-representative longitudinal study of Minnesota families (11). The younger twin cohort (N=2,510; 51% female) was first assessed at age 11 during the years 1988–2005. The older cohort (N=1,252; 54% female) was first assessed at age 17 during the years 1989–1996. Members of the age 11 cohort were invited to participate in follow-up assessments at ages 14 and 17, and all twins were invited to participate in follow-up assessments at ages 20, 24, and 29. Cohorts were combined for all analyses. Participants received modest payments for their assessments. Written assent or consent was obtained from all participants, including the parents of minor children, and all study protocols were approved by the University of Minnesota institutional review board.
Additional analyses were conducted with a subsample of participants who had at least one symptom by their assessment at age 17 for nicotine, alcohol, or marijuana dependence. This resulted in an early-use subsample of 580 male and 486 female participants.
Pooling across cohorts in the full sample, the actual mean ages at assessment were 11.8 (SD=0.4), 14.9 (SD=0.6), 17.8 (SD=0.7), 21.1 (SD=0.8), 25.0 (SD=0.9), and 29.5 (SD=0.7) years. Participation rates ranged from 87.3% to 93.6% for the follow-up assessments. To examine attrition, we compared the results from 17-year-olds who completed the adult assessments at ages 20, 24, and 29 against those who did not. For male participants, the Cohen’s d values for mean differences in dependence symptoms at age 17 between those who completed the later assessments (N=1,570) and those who did not (N=238) were 0.00, −0.08, and 0.09 for nicotine, alcohol, and marijuana, respectively. For female participants, the Cohen’s d values for similar comparisons were −0.19, −0.01, and 0.13 (p>0.05 in all cases).
Measures
Diagnostic symptom counts were obtained during in-person interviews with trained interviewers. In a consensus process, graduate students and staff with advanced training in clinical assessment reviewed cases to verify the presence of symptoms.
Assessments at ages 11 and 14 used the Diagnostic Interview for Children and Adolescents (12) to examine DSM-III-R nicotine dependence, alcohol dependence or abuse, and marijuana dependence or abuse. Assessments at later ages used a modified version of the Substance Abuse Module (13) of the Composite International Diagnostic Interview (CIDI [14]) to assess DSM-III-R symptoms of substance use disorders. Abuse and dependence symptoms were collapsed for alcohol and marijuana. We also obtained mothers’ reports of their children’s symptoms at ages 11, 14, and 17. The follow-up assessments at each age covered the interval that had elapsed since the last assessment. We used a best-estimate approach (15) whereby a symptom was considered present if reported by either the child or the mother. Diagnostic interrater reliability of substance use disorders was greater than 0.91 (16). To rule out possible informant effects, all analyses were repeated using only the child as the informant, and the pattern of results was identical.
Analysis of Change in Comorbidity
Bivariate correlations were computed using Spearman’s rank-order correlation statistic, which is robust to departures from bivariate normality (17). We used confirmatory factor analysis (18) to model the pattern of correlations among substance use disorders over time. For each assessment age, a single factor was fitted to account for correlations among symptom counts of alcohol, nicotine, and marijuana dependence. Because of prohibitively low variance, symptoms at age 11 were not included in the model. This resulted in a model with five general factors, one for each age of assessment and each representing the covariance among substance use disorders at each age. All factors were allowed to correlate, and all same-drug residuals were allowed to correlate across ages to account for the within-person correlated nature of the longitudinal data. Loadings were standardized. Thus, the variance of the symptom count variables can be modeled as a function of the general and residual (substance-specific) factors:
where “var” denotes variance. Since the variance of the general factor was set to 1, the variance in a standardized symptom count accounted for by the general factor is simply the square of the factor loading. To obtain a single estimate of comorbidity for each age of assessment (i.e., a single estimate of the variance in the three symptom counts accounted for by the general factor), we calculated the mean squared factor loadings at each age. The mean squared loading provides a reasonable metric to test for changes in correlations because it is directly proportional to the magnitude of correlations among the symptom count variables (i.e., higher correlations among the symptom counts result in greater mean variance in the symptom counts accounted for by the general factor). To test for significant differences, we constrained the mean squared loading to be the same across assessment waves, and we conducted likelihood ratio tests to evaluate the change in model fit according to standard practice (19).
In addition, because the sample was composed of twins, we used standard twin modeling to decompose the general factor variances into three components: additive genetic (A), shared environmental (C), and nonshared environmental (E) components. All analyses were conducted in the same model rather than separate phenotypic and biometric models. Because we set the variance of the general factor to 1, the sum of the A, C, and E variance components of the general factor is also 1. To obtain the average variance in the symptom counts accounted for by genetic influences on the factor, one simply substitutes the term “var(general factor)” in the above equation with the corresponding genetic, shared environmental, or nonshared environmental component of the general factor variance. For example, a mean squared loading of 0.4 at age 17 would indicate that 40% of the variance in the symptom counts was due to the general factor. If the additive genetic component A of the age 17 factor was 0.7, then the genetic influence of the general factor onto the symptom counts would be 0.4×0.7=0.28. That is, 28% of the symptom count variance would be due to the genetic influence of the factor.
Investigating the biometric decomposition of symptom count correlations is only applicable to the full sample because the subsample analysis of individuals symptomatic by age 17 was, by definition, a within-individual analysis and disregards unaffected or later-affected (after age 17) co-twins. Resulting biometric decompositions would be difficult to interpret because only those pairs of twins who were symptomatic by age 17 were included in the cross-twin correlations.
Analyses were conducted using R, version 2.10.1 (20), and the factor analysis was conducted using OpenMx, version 1.0.6 (21). Missing data was handled using full information maximum likelihood. We evaluated model fit using chi-square tests and the difference in the Akaike information criterion (22) between the saturated and alternative models (positive values indicate that the alternative provides better fit) separately for the four samples. We also report the root mean square error of approximation, where values of 0.06 or less indicate a very good fit (23). The chi-square test provides an exact test of model fit, but it is sensitive to sample size and the magnitude of correlations among measures and is always significant (p<0.05) when statistical power is high (e.g., in large samples). It is often significant even when the model provides an accurate and useful representation of the data (24). Therefore, other fit indices have been developed and are primarily used to evaluate model fit. The root mean square error of approximation attempts to correct deficiencies of the chi-square by adjusting for the degrees of freedom and sample size. Conventional root mean square error of approximation cutoffs are 0.08 for good fit and 0.05 for very good fit. The Akaike information criterion has strong theoretical properties in that the selected model is expected to fit best upon cross-validation (22).
Results
Longitudinal trends in symptom counts were somewhat different for male and female participants (Figure 1 and Table 1). For male participants, symptom counts significantly increased from age 11 to late adolescence, peaked during the early twenties, and declined thereafter. For female participants, a more prolonged plateau during the late teens and early twenties was observed, with a marked drop only by age 29. Mean symptom counts increased more rapidly for male than for female participants, with teenage boys maintaining higher mean-level symptoms after age 14. Variances also increased during adolescence and decreased during the late twenties.
| Symptom Count | |||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| N | % Ever Used Substance | Prevalence of Dependence (%) | Prevalence of Subclinical Dependence (% With One or Two Symptoms) | Mean | SD | Sex Differences in Mean Symptom Countb | Sex Differences in Variance of Symptom Countc | ||||||||
| Assessment Age (Years) | Substance | M | F | M | F | M | F | M | F | M | F | M | F | Effect Size | M/F Ratio |
| 11 | Nicotine | 1,233 | 1,282 | 9 | 4 | 0.2 | 0.1 | 0.3 | 0.0 | 0.02 | 0.00 | 0.25 | 0.08 | N/A | N/A |
| Alcohol | 1,233 | 1,281 | 8 | 2 | 0.0 | 0.0 | 0.1 | 0.0 | 0.00 | 0.00 | 0.03 | 0.00 | N/A | N/A | |
| Marijuana | 1,233 | 1,282 | 1 | 0.4 | 0.0 | 0.0 | 0.0 | 0.0 | 0.00 | 0.00 | 0.03 | 0.00 | N/A | N/A | |
| 14 | Nicotine | 1,148 | 1,191 | 35 | 27 | 4.5 | 4.6 | 3.5 | 2.8 | 0.22 | 0.31 | 0.99 | 1.15 | −0.08 | 0.74*** |
| Alcohol | 1,132 | 1,191 | 53 | 41 | 1.7 | 0.6 | 3.1 | 3.2 | 0.08 | 0.08 | 0.54 | 0.55 | 0.00 | 0.96 | |
| Marijuana | 1,132 | 1,191 | 13 | 10 | 1.7 | 1.4 | 2.4 | 1.5 | 0.09 | 0.08 | 0.67 | 0.62 | 0.02 | 1.17* | |
| 17 | Nicotine | 1,129 | 1,644 | 69 | 54 | 17.1 | 12.5 | 11.2 | 8.6 | 0.96 | 0.72 | 1.82 | 1.62 | 0.14** | 1.26*** |
| Alcohol | 1,420 | 1,644 | 79 | 77 | 11.1 | 4.9 | 23.6 | 13.3 | 0.74 | 0.37 | 1.55 | 1.19 | 0.27*** | 1.70*** | |
| Marijuana | 1,423 | 1,644 | 34 | 30 | 9.2 | 3.2 | 9.5 | 5.7 | 0.59 | 0.26 | 1.67 | 1.02 | 0.24*** | 2.68*** | |
| 20 | Nicotine | 1,393 | 1,331 | 87 | 72 | 27.4 | 17.0 | 15.1 | 12.0 | 1.50 | 0.90 | 1.95 | 1.63 | 0.33*** | 1.43*** |
| Alcohol | 1,101 | 1,331 | 96 | 96 | 24.6 | 6.2 | 30.0 | 15.6 | 1.56 | 0.44 | 1.93 | 1.12 | 0.73*** | 2.97*** | |
| Marijuana | 1,111 | 1,330 | 55 | 45 | 11.8 | 3.8 | 13.4 | 6.8 | 0.84 | 0.27 | 1.69 | 0.97 | 0.42*** | 3.04*** | |
| 24 | Nicotine | 1,009 | 1,213 | 90 | 77 | 26.4 | 17.3 | 14.1 | 12.5 | 1.34 | 0.90 | 1.83 | 1.60 | 0.26*** | 1.31*** |
| Alcohol | 1,148 | 1,213 | 98 | 98 | 23.0 | 6.9 | 32.0 | 17.8 | 1.39 | 0.52 | 1.76 | 1.23 | 0.58*** | 2.05*** | |
| Marijuana | 1,143 | 1,211 | 60 | 50 | 7.4 | 3.1 | 13.2 | 5.4 | 0.56 | 0.20 | 1.39 | 0.85 | 0.31*** | 2.67*** | |
| 29 | Nicotine | 957 | 636 | 90 | 79 | 23.1 | 15.7 | 16.2 | 10.9 | 1.22 | 0.79 | 1.77 | 1.55 | 0.26*** | 1.30** |
| Alcohol | 957 | 636 | 98 | 97 | 12.9 | 4.5 | 27.4 | 11.2 | 0.89 | 0.31 | 1.52 | 0.95 | 0.44*** | 2.56*** | |
| Marijuana | 955 | 635 | 60 | 46 | 5.5 | 0.6 | 8.0 | 3.7 | 0.42 | 0.09 | 1.25 | 0.52 | 0.32*** | 5.78*** | |
a M=male participants; F=female participants.
b
Sex differences in mean symptom count were evaluated with a likelihood ratio test correcting for within-family correlations. Values reported are Cohen’s d. Prevalence and symptom counts were measured in a since-last-assessment format.
c
The equality of variance across sex is given as a simple ratio of male variance to female variance, but again tested by correcting for within-family correlations.
*p<0.05. **p<0.01. ***p<0.001.
The prevalence rates for each disorder and the proportion of individuals who had ever used each substance for each age are reported in Table 1. The mean age at initiation was 14.4 years (SD=3.4) for nicotine, 15.5 years (SD=2.6) for alcohol use without parental permission, and 16.5 years (SD=2.5) for marijuana. The mean age at symptom onset was 17.5 years (SD=2.7) for nicotine, 18.2 years (SD=4.4) for alcohol, and 17.2 years (SD=2.3) for marijuana.
Cross-drug correlations decreased with age, as can be seen in the full bivariate correlation matrix in Figure 2. For example, for male participants, the correlation between nicotine and alcohol declined from 0.61 at age 14 to 0.26 at age 29.
FIGURE 2. Within and Across Age Correlations Between Substance Use Symptom Count Measuresa
a Results from male participants are reported in the lower triangle and from female participants in the upper triangle. Correlations (without decimals) are displayed within each colored box. To aid visualization, the matrix is a heat map, with hotter colors signifying higher correlations. The matrix is organized into blocks by age. Note the trend in the bolded diagonal blocks; the colors generally become cooler as one moves from the upper left to lower right, indicating a steady decrease in correlations among the substances over time. The off-diagonals have purposefully been partially obscured to focus the reader’s attention on the block diagonal without omitting relevant information about the cross-age correlations. Female participants in the younger cohort had just begun their age 29 assessment, and thus the age 14/age 29 block is empty.
Factor models were fitted to provide a formal test of changes in the cross-drug correlations. The model fit in the full male sample was very good (χ2=216.66, df=170, p=0.009; ΔAkaike information criterion=123.33; root mean square error of approximation=0.02). Model fit in the full female sample was also good (χ2=432.32, df=170, p<0.001; ΔAkaike information criterion=−92.32; root mean square error of approximation=0.05). The model fit was good in both the male (χ2=137.45, df=170, p=0.97; ΔAkaike information criterion=202.55; root mean square error of approximation<0.01) and female (χ2=133.42, df=106, p=0.04; ΔAkaike information criterion=78.6; root mean square error of approximation=0.04) early-use subsamples.
For each model, the size of the standardized factor loadings generally decreased with age, declining from 0.6–0.8 at age 14 to 0.3–0.6 at age 29, as reported in Table 2. The table also shows (as do the gray lines in Figure 3) that the mean squared loading, which served as our metric of change, decreased with age. Recall that the mean squared factor loading provides the average variance in the symptoms accounted for by the general factor. To illustrate how these estimates are calculated, we entered the loadings from Table 2 into equation 1 for the full female sample. The mean squared loading was (0.652+0.712+0.792)/3=0.52 at age 14 and (0.692+0.532+0.332)/3=0.29 at age 29. To test whether the decrease from age 14 to age 29 was statistically significant, we constrained the mean squared loadings to be same at each age and examined decrement in model fit. The likelihood ratio test indicated a significant decrease in the mean squared loading from age 14 to age 29 for all samples: full male sample (χ2=158.4, df=4, p=3.2×10−33), full female sample (χ2=146.97, df=4, p=9.1×10−31), male subsample of early users (χ2=54.66, df=4, p=3.8×10−11), and from age 14 to age 24 for the female subsample of early users (χ2=21.36, df=3, p=8.9×10−5). We tested for sex differences by comparing the fit of an unconstrained model to that of a model that constrained the mean squared loading to be equal across male and female participants. Relative to young men, young women showed an earlier decrease in comorbidity at age 20 (χ2=24.18, df=1, p=8.8×10−7), but not at any other ages.
| Factor Loading | Factor Variance Componentb | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Sample and Assessment Age | Nicotine | Alcohol | Marijuana | Mean Squared Loading | A | C | E | |||||||
| MLE | 95% CI | MLE | 95% CI | MLE | 95% CI | MLE | 95% CI | MLE | 95% CI | MLE | 95% CI | MLE | 95% CI | |
| Full female sample (N=1,907) | ||||||||||||||
| Age 14 | 0.65 | 0.59–0.70 | 0.71 | 0.66–0.83 | 0.79 | 0.74–0.84 | 0.52 | 0.48–0.56 | ||||||
| Age 17 | 0.60 | 0.55–0.64 | 0.78 | 0.74–0.81 | 0.73 | 0.68–0.83 | 0.51 | 0.47–0.54 | 0.72 | 0.56–0.88 | 0.14 | 0–0.29 | 0.14 | 0.10–0.19 |
| Age 20 | 0.51 | 0.45–0.58 | 0.57 | 0.50–0.64 | 0.60 | 0.54–0.65 | 0.32 | 0.28–0.36 | 0.42 | 0.14–0.74 | 0.21 | 0–0.47 | 0.37 | 0.25–0.49 |
| Age 24 | 0.60 | 0.53–0.60 | 0.54 | 0.46–0.61 | 0.39 | 0.31–0.45 | 0.27 | 0.24–0.31 | 0.61 | 0.31–0.61 | 0.09 | 0–0.34 | 0.30 | 0.19–0.42 |
| Age 29 | 0.69 | 0.60–0.79 | 0.53 | 0.43–0.62 | 0.33 | 0.25–0.41 | 0.29 | 0.24–0.33 | 0.36 | 0.09–0.61 | 0.12 | 0–0.39 | 0.53 | 0.36–0.68 |
| Full male sample (N=1,775) | ||||||||||||||
| Age 14 | 0.74 | 0.69–0.77 | 0.76 | 0.71–0.80 | 0.77 | 0.73–0.81 | 0.57 | 0.54–0.61 | 0.35 | 0.16–0.57 | 0.46 | 0.26–0.63 | 0.19 | 0.14–0.24 |
| Age 17 | 0.70 | 0.65–0.74 | 0.79 | 0.75–0.82 | 0.72 | 0.67–0.76 | 0.54 | 0.51–0.57 | 0.69 | 0.50–0.88 | 0.18 | 0.03–0.36 | 0.14 | 0.10–0.18 |
| Age 20 | 0.65 | 0.60–0.70 | 0.67 | 0.62–0.72 | 0.70 | 0.64–0.75 | 0.45 | 0.42–0.49 | 0.58 | 0.37–0.73 | 0.13 | 0.01–0.32 | 0.29 | 0.22–0.38 |
| Age 24 | 0.55 | 0.48–0.61 | 0.56 | 0.50–0.62 | 0.56 | 0.49–0.62 | 0.31 | 0.27–0.35 | 0.54 | 0.30–0.73 | 0.16 | 0.01–0.39 | 0.30 | 0.19–0.42 |
| Age 29 | 0.47 | 0.41–0.54 | 0.62 | 0.55–0.69 | 0.59 | 0.52–0.65 | 0.32 | 0.28–0.36 | 0.47 | 0.47–0.63 | 0.05 | 0–0.26 | 0.48 | 0.34–0.61 |
| Female subsamplec (N=486) | ||||||||||||||
| Age 14 | 0.49 | 0.35–0.61 | 0.64 | 0.52–0.74 | 0.82 | 0.69–0.93 | 0.44 | 0.36–0.51 | ||||||
| Age 17 | 0.32 | 0.21–0.42 | 0.58 | 0.46–0.69 | 0.77 | 0.65–0.89 | 0.34 | 0.29–0.40 | ||||||
| Age 20 | 0.23 | 0.12–0.35 | 0.52 | 0.32–0.73 | 0.68 | 0.49–0.91 | 0.26 | 0.20–0.34 | ||||||
| Age 24 | 0.31 | 0.17–0.44 | 0.63 | 0.48–0.76 | 0.43 | 0.28–0.59 | 0.22 | 0.16–0.29 | ||||||
| Male subsamplec (N=580) | ||||||||||||||
| Age 14 | 0.69 | 0.61–0.76 | 0.73 | 0.64–0.80 | 0.76 | 0.68–0.83 | 0.53 | 0.47–0.59 | ||||||
| Age 17 | 0.47 | 0.37–0.56 | 0.62 | 0.52–0.72 | 0.65 | 0.55–0.76 | 0.35 | 0.29–0.41 | ||||||
| Age 20 | 0.49 | 0.39–0.59 | 0.57 | 0.46–0.67 | 0.68 | 0.55–0.81 | 0.35 | 0.28–0.41 | ||||||
| Age 24 | 0.40 | 0.28–0.52 | 0.49 | 0.36–0.62 | 0.48 | 0.34–0.62 | 0.21 | 0.15–0.28 | ||||||
| Age 29 | 0.43 | 0.21–0.44 | 0.61 | 0.48–0.75 | 0.64 | 0.50–0.79 | 0.30 | 0.24–0.36 | ||||||
a MLE=the maximum likelihood estimate. The mean squared loading is the mean of the squared loadings listed for that age. Age 29 data are not available for the female subsample because the assessment is ongoing.
b
A=additive genetic variance component of the factors; C=environmental component; E=nonshared environmental component. A, C, and E estimates for the 14-year-old samples were poorly estimated because of lack of cross-twin covariance in the symptom counts at that age and are not provided for the female participants. A, C, and E estimates are not provided for the subsamples of early-onset users because these samples are, by definition, within individual and exclude co-twins who did not exhibit symptoms at the age 17 assessment. Resulting ACE estimates are therefore difficult to interpret.
c
Subsamples are participants with at least one symptom by age 17.
FIGURE 3. Percent of Symptom Count Variance Accounted for by the General Factor at Each Agea
a The gray lines in the figure show the decrease in correlations over time as expressed by decrease in the average percentage of variance accounted for by a general factor at each time point (i.e., the “mean squared loading” column in Table 2). Error bars indicate 95% confidence intervals. The full sample is shown on the left while the subsample of individuals who had at least one nicotine, alcohol, or marijuana dependence symptom by their age 17 assessment is shown on the right. The decline for both sexes in both samples was statistically significant (see Results). In the full male and female samples the blue, green, and red lines represent the proportions of phenotypic variance (gray lines) that are due to additive genetic (blue; A), shared environmental (red; C), and nonshared environmental (green; E) influence of the general factor at each age. Male values are always represented by the darker hue. These estimates can be computed directly from values provided in Table 2 by multiplying the mean squared loading by the corresponding ACE value (e.g., by the A value to obtain the additive genetic variance, which is plotted in blue). The majority of phenotypic decline in the full sample is because of a statistically significant decline in heritability as well as a statistically nonsignificant decline in shared environmental variance. In contrast, nonshared environmental variance significantly increased with age (see text). Note that age 14 estimates are not given because girls at this age had few symptoms. Genetic and environmental components are not given for the subsample as it was composed of selected individuals and not selected twin pairs.
Table 2 lists the additive genetic, shared environmental, and nonshared environmental components of each age’s general factor for each sample. While these estimates are useful, they must be multiplied by their respective factor loadings to determine the genetic and environmental impact on the individual symptom counts. To do this, we replaced the var(factor) term in equation 1 with the A, C, and E variance components. The resulting values are displayed for each age in Figure 3. Note that the scaling is such that adding the A, C, and E components reproduces the mean squared loading (the mean squared loading is also listed in Table 2). For example, for men at age 29, the phenotypic variance accounted for by the general factor was 32%, which was composed of 15% (0.32×0.47×100%) additive genetic variance, 15% (0.32×0.48×100%) nonshared environmental variance, and 2% (0.32×0.05×100%) shared environmental variance. Unfortunately, symptoms were expressed relatively infrequently by both members of a female twin pair at age 14, rendering the age 14 additive genetic, shared environmental, and nonshared environmental components poorly estimated and hence excluded from Figure 3 (along with the male estimates for consistency).
We tested for change in heritability, shared environmental, and nonshared environmental estimates (see Figure 3) by fixing relevant parameters to be equal at ages 17, 20, 24, and 29, testing for decrements in model fit. For example, to test for a change in heritability, we fixed the A component multiplied by the mean squared loading of each general factor to be equal, and we tested the fit of this model against the original model where all components are freely estimated. The decrease in heritability was significant for male (χ2=23.86, df=3, p=2.7×10−5) and female (χ2=26.6, df=3, p=7.0×10−6) participants. Thus, it appears that the vast majority of the phenotypic decline is the result of a decline in genetic variance. The increase in nonshared environment was also significant for male (χ2=15.73, df=3, p=0.001) and female (χ2=11.48, df=3, p=0.009) participants. Changes in shared environment did not reach significance for either sex.
Discussion
We examined changes in the correlations among symptoms of nicotine, alcohol, and marijuana abuse and dependence from age 11 to age 29. As shown in Figure 3, male and female participants experienced significant declines in these correlations from adolescence to adulthood. This was also true for an early-onset subsample of individuals with at least one symptom of nicotine, alcohol, or marijuana dependence by age 17. Correlation among disorders or their symptoms is evidence that those disorders share etiology. Our results suggest that the shared etiology contributing to nicotine, alcohol, and marijuana dependence symptoms diminish over time. That is, younger individuals tended to use these three substances indiscriminately, whereas older individuals began to show a preference for one substance over the others. Despite decreases in the correlations, the rates of use (Figure 1) continued to climb throughout late adolescence and early adulthood. Finally, after age 17 the correlations among symptoms became less attributable to pleiotropic genetic effects and increasingly a result of nonshared environmental influences (Figure 3), indicating that the types of etiological processes contributing to the variation in use of multiple drugs gradually changes during the transition to adulthood.
Several processes might contribute to the transition from general to specific influences. For one, adolescents are more impulsive and risk-taking than adults (25). Personality traits such as disinhibition and sensation-seeking are not predispositions to the use of any particular substance but rather the use of whatever substances might be available (9, 26). Additionally, neurodevelopmental changes relevant to behavioral disinhibition continue throughout adolescence. For example, by adolescence, the nucleus accumbens—important in reward sensitivity—is well developed but poorly regulated by a still maturing prefrontal cortex (27, 28), resulting in deficits of top-down control over the reward system that slowly improve in early adulthood. This developmental window is the same time period in which we observed decreases in the comorbidity among different substance use disorders, suggesting that disinhibitory mechanisms with known neurological substrates may be in play. Further supporting this hypothesis is that the earlier decline in comorbidity for female relative to male participants from age 17 to age 20 (Figure 3) is consistent with the earlier pubertal (29), cortical (30), and personality (31) maturation in girls compared with boys.
Our results are consistent with individual differences in drug reinforcement. That is, initial drug use that begins in adolescence is characterized by relatively indiscriminate experimentation. Because of individual differences in the reinforcing effects of different drugs, however, people may eventually restrict their use to those drugs that provide the greatest reinforcement. There are myriad etiological mechanisms relevant to individual differences in drug reinforcement, including drug metabolism effects, drug availability, social rewards and punishments, and differences in drug-specific neurological sensitivity that may be further moderated by drug and alcohol neurotoxic mechanisms (26, 32). Future research might consider these covariates in characterizing the transition to substance specialization.
While our results are consistent with the theory that disinhibition accounts for comorbidity among substance use disorders, they are not entirely inconsistent with the gateway hypothesis (33). The gateway hypothesis holds that using one drug leads to using other drugs, perhaps as a result of the experienced high and a greater desire to obtain bigger and better highs. If correct, the theory would predict the opposite pattern of comorbidity we observed; correlations should be low at younger ages and increase over time. At younger ages, most people would not yet have used their first or second gateway drug, and they would not have had time to explore other drugs. Over time, the correlations among drugs should increase as the initial gateway drug use would cause them to use other drugs. In fact, we found that as drug dependence symptoms increased (Figure 1), the magnitude of the associations among different drugs decreased (Figure 3); this is opposite to the pattern predicted by the gateway theory. This conclusion is consistent with a growing body of literature showing results that are inconsistent with parts of the gateway theory (1, 9, 34–36). That said, the increase in nonshared environmental effects observed in our sample could contain etiology analogous to a gateway process (i.e., environments that contribute to multiple drug dependence in one person but not in that person’s twin). While the range of possible environments is vast and the effect is small, it could very well include drug use (e.g., impaired cognition caused by neurotoxicity) or other risk factors for nonspecific drug use, such as occupational, social, and legal problems.
Our study has some limitations. The use of a community-representative sample, while itself a strength, may not apply to clinical populations. The fact that our results were consistent for the subsample of individuals who were symptomatic at age 17 suggests this is not a major concern. We used symptom counts because symptoms are relevant to the DSM clinical literature and provide measurement consistency from ages 11 to 29. In supplementary analyses we evaluated other measures of quantity and frequency of nicotine, alcohol, and marijuana use (not shown). These measures had higher means and variances at younger ages, and analyses resulted in the same trends as presented here for symptoms.
In the United States, adolescent development is confounded with substance-use-relevant environmental changes. Adolescents gradually experience greater autonomy and financial freedom from caregivers. The purchase of tobacco and alcohol becomes legal at ages 18 and 21, respectively, while marijuana use is always illegal. Undoubtedly, these environmental influences affected symptom means and correlations. However, one cannot disentangle these influences from other behavioral and neurological changes, at least within a single culture. Cross-cultural and cross-generational studies (e.g., comparisons between societies that differ in drug laws) are required to unravel maturational and environmental changes during development.
References
1.
Kendler KS, Jacobson KC, Prescott CA, Neale MC: Specificity of genetic and environmental risk factors for use and abuse/dependence of cannabis, cocaine, hallucinogens, sedatives, stimulants, and opiates in male twins. Am J Psychiatry 2003; 160:687–695
2.
Kendler KS, Prescott CA, Myers J, Neale MC: The structure of genetic and environmental risk factors for common psychiatric and substance use disorders in men and women. Arch Gen Psychiatry 2003; 60:929–937
3.
Kessler RC, Chiu WT, Demler O, Merikangas KR, Walters EE: Prevalence, severity, and comorbidity of 12-month DSM-IV disorders in the National Comorbidity Survey Replication. Arch Gen Psychiatry 2005; 62:617–627
4.
Moffitt TE, Caspi A, Taylor A, Kokaua J, Milne BJ, Polanczyk G, Poulton R: How common are common mental disorders? evidence that lifetime prevalence rates are doubled by prospective versus retrospective ascertainment. Psychol Med 2010; 40:899–909
5.
Krueger RF, Hicks BM, Patrick CJ, Carlson SR, Iacono WG, McGue M: Etiologic connections among substance dependence, antisocial behavior, and personality: modeling the externalizing spectrum. J Abnorm Psychol 2002; 111:411–424
6.
Krueger RF, Markon KE: Reinterpreting comorbidity: a model-based approach to understanding and classifying psychopathology. Annu Rev Clin Psychol 2006; 2:111–133
7.
Meehl PE: Comorbidity and taxometrics. Clin Psychol-Sci Pr 2001; 8:507–519
8.
Meehl PE: Factors and taxa, traits and types, differences of degree and differences in kind. J Pers 1992; 60:117–174
9.
Iacono WG, Malone SM, McGue M: Behavioral disinhibition and the development of early-onset addiction: common and specific influences. Annu Rev Clin Psychol 2008; 4:325–348
10.
Young SE, Corley RP, Stallings MC, Rhee SH, Crowley TJ, Hewitt JK: Substance use, abuse and dependence in adolescence: prevalence, symptom profiles and correlates. Drug Alcohol Depend 2002; 68:309–322
11.
Iacono WG, McGue M: Minnesota Twin Family Study. Twin Res 2002; 5:482–487
12.
Welner Z, Reich W, Herjanic B, Jung KG, Amado H: Reliability, validity, and parent-child agreement studies of the Diagnostic Interview for Children and Adolescents (DICA). J Am Acad Child Adolesc Psychiatry 1987; 26:649–653
13.
Robins LN, Babor TF, Cottler LB: Composite International Diagnostic Interview: Expanded Substance Abuse Module. St Louis, Washington University, Department of Psychiatry, 1987
14.
Robins LN, Wing J, Wittchen HU, Helzer JE, Babor TF, Burke J, Farmer A, Jablenski A, Pickens R, Regier DA, Sartorius N, Towle LH: The Composite International Diagnostic Interview: an epidemiologic instrument suitable for use in conjunction with different diagnostic systems and in different cultures. Arch Gen Psychiatry 1988; 45:1069–1077
15.
Leckman JF, Sholomskas D, Thompson WD, Belanger A, Weissman MM: Best estimate of lifetime psychiatric diagnosis: a methodological study. Arch Gen Psychiatry 1982; 39:879–883
16.
Iacono WG, Carlson SR, Taylor J, Elkins IJ, McGue M: Behavioral disinhibition and the development of substance-use disorders: findings from the Minnesota Twin Family Study. Dev Psychopathol 1999; 11:869–900
17.
Caruso JC, Cliff N: Empirical size, coverage, and power of confidence intervals for Spearman’s rho. Educ Psychol Meas 1997; 57:637–654
18.
Brown T: Confirmatory Factor Analysis for Applied Research. New York, Guilford, 2006
19.
Miller MB, Neale MC: Using nonlinear constraints in Mx to automate estimation of confidence intervals for parameters in genetic models. Behav Genet 1995; 25:279–280
20.
R Development Core Team: R: A Language and Environment for Statistical Computing. Vienna, R Foundation for Statistical Computing, 2011
21.
Boker S, Neale M, Maes H, Wilde M, Spiegel M, Brick T, Spies J, Estabrook R, Kenny S, Bates T, Mehta P, Fox J: OpenMx: an open source extended structural equation modeling framework. Psychometrika 2011; 76:306–317
22.
Vrieze SI: Model selection and psychological theory: a discussion of the differences between the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). Psychol Methods 2012; 17:228–243
23.
Hu LT, Bentler PM: Cutoff criteria for fit indexes in covariance structure analysis: conventional criteria versus new alternatives. Struct Equ Modeling 1999; 6:1–55
24.
Bentler PM: Comparative fit indexes in structural models. Psychol Bull 1990; 107:238–246
25.
Steinberg L, Cauffman E, Woolard J, Graham S, Banich M: Are adolescents less mature than adults? minors’ access to abortion, the juvenile death penalty, and the alleged APA “flip-flop.” Am Psychol 2009; 64:583–594
26.
Bechara A: Decision making, impulse control and loss of willpower to resist drugs: a neurocognitive perspective. Nat Neurosci 2005; 8:1458–1463
27.
Casey BJ, Getz S, Galvan A: The adolescent brain. Dev Rev 2008; 28:62–77
28.
Casey BJ, Tottenham N, Liston C, Durston S: Imaging the developing brain: what have we learned about cognitive development? Trends Cogn Sci 2005; 9:104–110
29.
Tanner JM, Whitehouse RH, Takaishi M: Standards from birth to maturity for height, weight, height velocity, and weight velocity: British children, 1965. II. Arch Dis Child 1966; 41:613–635
30.
Lenroot RK, Gogtay N, Greenstein DK, Wells EM, Wallace GL, Clasen LS, Blumenthal JD, Lerch J, Zijdenbos AP, Evans AC, Thompson PM, Giedd JN: Sexual dimorphism of brain developmental trajectories during childhood and adolescence. Neuroimage 2007; 36:1065–1073
31.
Klimstra TA, Hale WW, Raaijmakers QA, Branje SJT, Meeus WHJ: Maturation of personality in adolescence. J Pers Soc Psychol 2009; 96:898–912
32.
Goldstein RZ, Volkow ND: Drug addiction and its underlying neurobiological basis: neuroimaging evidence for the involvement of the frontal cortex. Am J Psychiatry 2002; 159:1642–1652
33.
Kandel DB, Jessor R: The gateway hypothesis revisited, in Stages and Pathways of Drug Involvement. Edited by, Kandel DB. New York, Cambridge University Press, 2002
34.
Tarter RE, Vanyukov M, Kirisci L, Reynolds M, Clark DB: Predictors of marijuana use in adolescents before and after licit drug use: examination of the gateway hypothesis. Am J Psychiatry 2006; 163:2134–2140
35.
Vanyukov MM, Tarter RE, Kirisci L, Kirillova GP, Maher BS, Clark DB: Liability to substance use disorders: 1. Common mechanisms and manifestations. Neurosci Biobehav Rev 2003; 27:507–515
36.
Irons DE, McGue M, Iacono WG, Oetting WS: Mendelian randomization: a novel test of the gateway hypothesis and models of gene-environment interplay. Dev Psychopathol 2007; 19:1181–1195
Information & Authors
Information
Published In
History
Received: 22 August 2011
Revision received: 1 February 2012
Revision received: 27 April 2012
Accepted: 24 May 2012
Published online: 1 October 2012
Published in print: October 2012
Authors
Funding Information
The authors report no financial relationships with commercial interests.Supported by grants DA 05147, DA 13240, DA 024417, and DA 025868 from the National Institute on Drug Abuse; grant AA 09367 from the National Institute on Alcohol Abuse and Alcoholism; and grant MH 017069 from NIMH.
Metrics & Citations
Metrics
Citations
Export Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
For more information or tips please see 'Downloading to a citation manager' in the Help menu.
View Options
View options
PDF/EPUB
View PDF/EPUBLogin options
Already a subscriber? Access your subscription through your login credentials or your institution for full access to this article.
Personal login Institutional Login Open Athens loginNot a subscriber?
PsychiatryOnline subscription options offer access to the DSM-5-TR® library, books, journals, CME, and patient resources. This all-in-one virtual library provides psychiatrists and mental health professionals with key resources for diagnosis, treatment, research, and professional development.
Need more help? PsychiatryOnline Customer Service may be reached by emailing [email protected] or by calling 800-368-5777 (in the U.S.) or 703-907-7322 (outside the U.S.).