-
PDF
- Split View
-
Views
-
Cite
Cite
Evan A Winiger, Spencer B Huggett, Alexander S Hatoum, Naomi P Friedman, Christopher L Drake, Kenneth P Wright, John K Hewitt, Onset of regular cannabis use and young adult insomnia: an analysis of shared genetic liability, Sleep, Volume 43, Issue 5, May 2020, zsz293, https://doi.org/10.1093/sleep/zsz293
- Share Icon Share
Abstract
Estimate the genetic and environmental influences on the relationship between onset of regular cannabis use and young adult insomnia.
In a population-based twin cohort of 1882 twins (56% female, mean age = 22.99, SD = 2.97) we explored the genetic/environmental etiology of the relationship between onset of regular cannabis use and insomnia-related outcomes via multivariate twin models.
Controlling for sex, current depression symptoms, and prior diagnosis of an anxiety or depression disorder, adult twins who reported early onset for regular cannabis use (age 17 or younger) were more likely to have insomnia (β = 0.07, p = 0.024) and insomnia with short sleep on weekdays (β = 0.08, p = 0.003) as young adults. We found significant genetic contributions for the onset of regular cannabis use (a2 = 76%, p < 0.001), insomnia (a2 = 44%, p < 0.001), and insomnia with short sleep on weekdays (a2 = 37%, p < 0.001). We found significant genetic correlations between onset of regular use and both insomnia (rA = 0.20, p = 0.047) and insomnia with short sleep on weekdays (rA = 0.25, p = 0.008) but no significant environmental associations between these traits.
We found common genetic liabilities for early onset of regular cannabis use and insomnia, implying pleiotropic influences of genes on both traits.
This study strengthens the small collection of research that shows an association of early cannabis use and sleep deficits in young adulthood, extending the effects to clinical sleep outcomes. Furthermore, it provides novel insight into shared genetics between early regular cannabis use and adult insomnia, implying a pleotropic genetic influence on both traits. Insomnia with short sleep is considered the most severe phenotype of insomnia and understanding developmental risk factors may prove to be useful for prevention efforts.
Introduction
Insomnia is one of the most common sleep disorders [1] with symptoms and diagnosis prevalence rates as high as 30% and 5%–10%, respectively, in the general population [2]. It is a burdensome diagnosis associated with increased health care costs [3], poor quality of life [4], and occupational impairments [5, 6]. The development, consistency, and severity of insomnia are often attributed to psychological factors: stress and cognitive-emotional arousal are frequently premorbid with diagnoses [7, 8], and insomnia has high comorbidity with psychiatric disorders like depression and anxiety [9, 10]. Cannabis use is cross-sectionally associated with increased rates of insomnia and insomnia symptoms [11–16]. Given evidence that early cannabis use can predict later sleep problems [17, 18], research is needed to understand the influence of early cannabis use on specific adult clinical sleep outcomes such as insomnia.
In addition to insomnia, cannabis use has been linked with poor subjective sleep quality [11, 15, 19–23], reduced time in the rapid eye movement phase of sleep [24, 25], eveningness (a preference for later sleep–wake timing) [26], later bed times [27], prolonged latency to sleep onset [28], and shorter sleep duration [27–33]. Both premorbid insomnia [12] and more generalized sleep problems [17, 34–38] significantly predict later cannabis use, but only a small number of studies provide evidence for early cannabis use predicting later sleep components such as tiredness, trouble sleeping [17], and sleep duration [18, 35], and no studies have focused on specific sleep disorders. With longitudinal evidence in both directions, it is possible that sleep problems could influence cannabis use, cannabis use could influence sleep problems, or an underlying shared liability such as common genetics could be responsible for their association.
Recent evidence suggests that genetics may play a role in the etiology of the relationship between early cannabis use and adult sleep problems [18], and that there could be a common genetic liability for cannabis use and sleep problems that explain their relationship. Several lines of evidence are consistent with this common genetic model: the presence of genes believed to be involved with circadian rhythm/sleep in genome-wide associations studies of lifetime cannabis use [39–41], clock gene genetic variants that are associated with cannabis addiction [42], and the possible role of the endocannabinoid system in the circadian rhythm/sleep–wake cycle [43–46]. Further research is needed on the possible shared genetics between cannabis use and clinical sleep outcomes such as insomnia.
Insomnia with short sleep duration is a subclassification of insomnia described as insomnia with an average sleep time of less than 6 h per night and elevated morbidity [47]. It is frequently conceptualized as the most biologically severe phenotype of the disorder [48] and is associated with distinctive consequences including biomarkers that are focused on physiologic hyperarousal [49–51], cardiometabolic morbidities and health related problems [52–55], and increased neurocognitive impairments [56, 57]. Furthermore, individuals with insomnia with short sleep may not respond as well to cognitive behavioral therapy for insomnia and might be at greater risk for chronic sleep disturbance [58], demonstrating the challenges of treating insomnia with short sleep as well as the need to identify developmental factors that could heighten vulnerability for this diagnosis. No prior study has looked at either the predictive phenotypic or potential genetic relationship between early cannabis use and insomnia disorder and insomnia with short sleep.
Twin studies can be used to estimate the genetic and environmental correlations between traits [59]. We used a young adult twin sample to dissect the relationship between the onset of regular cannabis use and both insomnia and insomnia with short sleep on weekdays. First, we tested if onset of regular cannabis use was associated with insomnia and insomnia with short sleep on weekdays, controlling for known correlates including sex [60], current depression symptoms [61], prior diagnoses of a depression or anxiety disorder [9, 10], and shift-work [62]. Next, we used univariate twin analyses to estimate the genetic and environmental etiological components of each trait individually. We then tested a sequence of nested bivariate twin models that decomposed the genetic and environmental overlap between onset of regular cannabis and both insomnia and insomnia with short sleep on weekdays.
Methods
Participants
Participants were 1882 individual twins (472 monozygotic [MZ] pairs, 304 dizygotic [DZ] same-sex pairs, 165 opposite-sex pairs, 56% female) from the Colorado community twin sample and longitudinal twin study [63] who completed both an online sleep survey and a substance use questionnaire. The average age was 22.99 years (range = 18–33, SD = 2.97) at the time of the online sleep survey, which was administered on average 2.31 years (SD = 1.75) before the substance use assessment (mean age = 25.29, range = 21–33, SD = 2.71). We excluded 39 participants due to reporting “small children keep awake” or being “pregnant” as reasons for sleep issues. Seven individuals were missing an insomnia classification because not enough information was provided to determine their diagnoses. Twenty-three subjects had missing responses on the sleep duration questionnaire. We assigned subjects NA for sleep duration if they responded as receiving more than 12 h (17 subjects) or 0 h of average weekday sleep (one subject). One subject was excluded due to an implausible response of age of first regular cannabis use being 1. All research protocols were reviewed and approved by the University of Colorado’s Investigational Review Board.
Procedure
Twins completed an online survey about sleep problems that took between 30 and 45 min (it also included the depression and anxiety questions used in our analysis). At a later time, they completed the online substance use questionnaire as part of a large battery of assessments that took between 2 and 3 h for the third wave of a longitudinal three-wave study.
Measures
Short sleep duration on weekdays.
Weekday self-reported sleep duration was assessed via a question asking, “During the past month, thinking about your average WEEKDAY, how long did you ACTUALLY sleep, EACH night (or your longest sleep period if you work a night shift or rotating shift)?” Participants entered responses in open text boxes for both hours and minutes. We averaged responses that were in the form of ranges (e.g. 7–9 h would become 8 h) (n = 41). We summed hours and minutes to create a total average weekday sleep measure (M = 7.07, SD = 1.31, Median = 7). We coded short sleep duration if the summed response was less than 6 h and compared outcomes to those with sleep duration responses that summed 6 or more. Weekday sleep was utilized because it captures the most normative insomnia sleep disturbances, free of the flexibility and variation of weekend schedules as well as the potential confounds from compensatory sleep of weekend sleep [54].
Insomnia.
We based insomnia diagnosis on DSM-IV-TR criteria [64]. The questionnaire asked how often in the past month participants experienced difficulty falling asleep, difficulty staying asleep, and had nonrefreshing sleep (never, sometimes, usually, or always). Participants who did not report “never” for any of these items were asked follow-up questions, including how long they had the sleep problem (years and months), and to what extent it interfered with their “daily functioning (daytime fatigue, ability to function at work/daily chores, concentration, memory, mood, etc.),” with answers “not at all interfering” = 0, “a little” =1, “somewhat” = 2, “much” = 3, and “very much interfering” = 4. Participants met criteria for insomnia if they responded to at least one of the three problems (falling asleep, staying asleep, or non-refreshing sleep) as “usually” or “always” for the duration of at least a month, and with it interfering with daily functioning at least “somewhat.” Insomnia (n = 391) was coded as a binary variable.
Insomnia with short sleep duration on weekdays.
A four-category ordinal variable was created to characterize insomnia with short sleep on weekdays for our twin models. We coded subjects who met criteria for both insomnia and short sleep duration on weekdays of less than 6 h as a 3 (n = 120), subjects with insomnia and with sleep duration 6 h or more as a 2 (n = 265), subjects who did not have insomnia but had short sleep duration on weekdays of less than 6 h as a 1 (n = 118), and subjects who did not have insomnia or sleep duration 6 h or more as a 0 (n = 1342). This ordinal categorization allowed us to score nearly all participants and resulted in similar group sizes for each category of insomnia.
Covariates of sleep.
Our phenotypic analyses controlled for known correlates of sleep: sex [60], current depression symptoms [61], prior diagnoses of a depression or anxiety disorder [9, 10], and shift-work [62]. Depression symptoms were assessed at the same time as the sleep questions with the Center for Epidemiological Studies-Depression (CES-D) scale [65]. We used a log-transformed sum of 19 items; the sleep disturbance question from the CES-D composite score was removed due to its direct overlap with our outcome measure. Prior anxiety and depressive disorder diagnoses were assessed at the time of the sleep survey via questions asking, “Have you ever been diagnosed with an anxiety disorder?” and “Have you ever been diagnosed with depression?” We coded anxiety and depression disorder as binary variables. Of those who endorsed having insomnia with short sleep on weekdays, 29% endorsed having a shift-work job such as a regular evening shift (n = 6), regular night shift (n = 8), or rotating shift (n = 21). As a form of quality control, we coded those who had a shift-work job as a 1 (n = 526) and those who did not as a 0 (n = 1348) and conducted a series of models with our other covariates that both included and excluded the shift-work variable as a control.
Onset of regular cannabis use.
Age of onset of regular cannabis use (M = 17.53 years, SD = 2.99, Median = 17) was assessed in a self-report addendum to the Composite International Diagnostic Interview Substance Abuse Module (CIDI-SAM) [66]. Participants who endorsed any lifetime cannabis use were asked “How old were you when you began using marijuana on a regular basis, that is at least once per month?” We coded responses as an ordinal variable: 2 = age of first regular use at age 17 or earlier as (early onset group) (n = 300); 1 = age of first regular use after age 17 (late onset group) (n = 246); 0 = never used cannabis regularly (n = 1336). This categorization was selected based on the mean and median age of 17 and because it divided the distribution of age of onset into approximately equal groups of early- and late-onset of regular cannabis use. Utilizing a three-category ordinal variable, rather than a binary categorization, improves the performance of multivariate genetic and environmental modeling [67].
Statistical analyses
Phenotypic models.
R version 3.5.1 [68] was used for all descriptive statistics. Mplus version 8.1 [69] was utilized for our phenotypic analysis, structural equation modeling, and comparison tests. To correct for the nonindependence of the twin pairs in the phenotypic analyses, we used the TYPE = COMPLEX command to cluster data by family. This method uses a weighted likelihood function to obtain scaled χ 2 and standard errors corrected for nonindependence using a sandwich estimator. This technique effectively corrects for the nonindependence of twin data [70]. Probit regression with a means and variances adjusted weighted least squares (WLSMV) estimator was used to examine predictors of insomnia and insomnia with short sleep on weekdays.
Univariate twin models.
We used biometrical twin ACE modeling [59] in Mplus to estimate the genetic and environmental contributions to the onset of regular cannabis use, insomnia, and insomnia with short sleep on weekdays. These models assume three latent factors are responsible for the variance of an individual trait: additive genetic effects (A), shared environmental (C), and nonshared or unique environmental factors (E). MZ twins correlate perfectly for A because they share all their all their additive genetic influences (alleles identical by descent). DZ twins correlate 0.5 for A because they share on average half of their alleles identical by descent. C is fixed to correlate at 1 for MZ and DZ twin pairs, by definition for the environmental influences they share. E is fixed to a zero correlation, by definition, for MZ and DZ twin pairs. Classical twin modeling leverages the differences between MZ and DZ pairs to estimate the contribution of each latent variable (A, C, or E) to a trait. Twins were randomly assigned to twin1 and twin2, but males from opposite sex twin pairs were assigned to twin1.
Because the data were ordinal, we assumed an underlying normal distribution and estimated threshold models using the WLSMV estimator [71]. We assessed fit with the omnibus χ 2 statistic, supplemented with root-mean-square error of approximation (RMSEA) and confirmatory fit index (CFI). RMSEA < 0.06 and CFI > 0.95 indicate good model fit [72]. We used χ 2 difference tests to check each individual trait for sex differences in the distributions of ordinal variable and the magnitude of A, C, and E estimates in males and females (scalar sex differences). We used χ 2 difference tests to determine best fitting models by checking if individual parameters (A, C, or E) could be dropped with no significant decrement in fit.
Bivariate twin models.
We used Cholesky decompositions to decompose the environmental/genetic variances and covariances between the onset of regular cannabis use and both insomnia and insomnia with short sleep on weekdays. These analyses allow us to estimate both the unique and overlapping contributions of genetic/environmental pathways between traits. Once the genetic overlap is estimated using the Cholesky decomposition, a simple conversion of the Cholesky path coefficients can provide the proportion of genetic variation shared between two traits, i.e. the genetic correlation:.
Results
Phenotypic analysis
Table 1 displays descriptive statistics and the frequencies of depressive symptoms, prior anxiety diagnoses, prior depression diagnoses, no insomnia and sleep 6 h or more on weekdays, no insomnia with short sleep on weekdays (<6 h), insomnia, and insomnia with short-sleep duration on weekdays for the full sample and by sex. Table 2 displays the descriptive statistics and frequencies for each group of regular cannabis onset (early, late, none).
Descriptive statistics and frequencies for all variables at the time of the sleep questionnaire (with the exception of age of onset of regular cannabis use which was taken from the substance use assessment)
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Full sample (n = 1882) . | Female (n = 1061) . | Male (n = 821) . |
Age (years) | 22.99 (2.97) | 22.79 (2.87) | 23.24 (3.07) |
Age of onset of regular cannabis use (years) | 17.53 (2.99) | 17.61 (3.16) | 17.48 (2.85) |
Sleep duration (h) | 7.08 (1.31) | 7.02 (1.31) | 7.14 (1.32) |
Current depression symptoms (CES-D) | 11.02 (8.62) | 12.05 (9.09) | 9.69 (7.77) |
Frequencies | |||
Full sample | Female | Male | |
Prior anxiety diagnosis | 9.15% (n = 171) | 12.26% (n = 129) | 5.14% (n = 42) |
Prior depression diagnosis | 15.19% (n = 286) | 19.22% (n = 202) | 10.27% (n = 84) |
No insomnia and sleep 6 h or more on weekdays | 72.74% (n = 1342) | 69.21% (n = 717) | 77.26% (n = 625) |
No insomnia with short sleep on weekdays (<6 h) | 6.39% (n = 118) | 6.37% (n = 66) | 6.43% (n = 52) |
Insomnia | 20.83% (n = 391) | 24.29% (n = 257) | 16.36% (n = 134) |
Insomnia with short sleep on weekdays (<6 h) | 6.50% (n = 120) | 7.14% (n = 74) | 5.69% (n = 46) |
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Full sample (n = 1882) . | Female (n = 1061) . | Male (n = 821) . |
Age (years) | 22.99 (2.97) | 22.79 (2.87) | 23.24 (3.07) |
Age of onset of regular cannabis use (years) | 17.53 (2.99) | 17.61 (3.16) | 17.48 (2.85) |
Sleep duration (h) | 7.08 (1.31) | 7.02 (1.31) | 7.14 (1.32) |
Current depression symptoms (CES-D) | 11.02 (8.62) | 12.05 (9.09) | 9.69 (7.77) |
Frequencies | |||
Full sample | Female | Male | |
Prior anxiety diagnosis | 9.15% (n = 171) | 12.26% (n = 129) | 5.14% (n = 42) |
Prior depression diagnosis | 15.19% (n = 286) | 19.22% (n = 202) | 10.27% (n = 84) |
No insomnia and sleep 6 h or more on weekdays | 72.74% (n = 1342) | 69.21% (n = 717) | 77.26% (n = 625) |
No insomnia with short sleep on weekdays (<6 h) | 6.39% (n = 118) | 6.37% (n = 66) | 6.43% (n = 52) |
Insomnia | 20.83% (n = 391) | 24.29% (n = 257) | 16.36% (n = 134) |
Insomnia with short sleep on weekdays (<6 h) | 6.50% (n = 120) | 7.14% (n = 74) | 5.69% (n = 46) |
CES-D = Center for Epidemiological Studies-Depression scale.
Descriptive statistics and frequencies for all variables at the time of the sleep questionnaire (with the exception of age of onset of regular cannabis use which was taken from the substance use assessment)
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Full sample (n = 1882) . | Female (n = 1061) . | Male (n = 821) . |
Age (years) | 22.99 (2.97) | 22.79 (2.87) | 23.24 (3.07) |
Age of onset of regular cannabis use (years) | 17.53 (2.99) | 17.61 (3.16) | 17.48 (2.85) |
Sleep duration (h) | 7.08 (1.31) | 7.02 (1.31) | 7.14 (1.32) |
Current depression symptoms (CES-D) | 11.02 (8.62) | 12.05 (9.09) | 9.69 (7.77) |
Frequencies | |||
Full sample | Female | Male | |
Prior anxiety diagnosis | 9.15% (n = 171) | 12.26% (n = 129) | 5.14% (n = 42) |
Prior depression diagnosis | 15.19% (n = 286) | 19.22% (n = 202) | 10.27% (n = 84) |
No insomnia and sleep 6 h or more on weekdays | 72.74% (n = 1342) | 69.21% (n = 717) | 77.26% (n = 625) |
No insomnia with short sleep on weekdays (<6 h) | 6.39% (n = 118) | 6.37% (n = 66) | 6.43% (n = 52) |
Insomnia | 20.83% (n = 391) | 24.29% (n = 257) | 16.36% (n = 134) |
Insomnia with short sleep on weekdays (<6 h) | 6.50% (n = 120) | 7.14% (n = 74) | 5.69% (n = 46) |
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Full sample (n = 1882) . | Female (n = 1061) . | Male (n = 821) . |
Age (years) | 22.99 (2.97) | 22.79 (2.87) | 23.24 (3.07) |
Age of onset of regular cannabis use (years) | 17.53 (2.99) | 17.61 (3.16) | 17.48 (2.85) |
Sleep duration (h) | 7.08 (1.31) | 7.02 (1.31) | 7.14 (1.32) |
Current depression symptoms (CES-D) | 11.02 (8.62) | 12.05 (9.09) | 9.69 (7.77) |
Frequencies | |||
Full sample | Female | Male | |
Prior anxiety diagnosis | 9.15% (n = 171) | 12.26% (n = 129) | 5.14% (n = 42) |
Prior depression diagnosis | 15.19% (n = 286) | 19.22% (n = 202) | 10.27% (n = 84) |
No insomnia and sleep 6 h or more on weekdays | 72.74% (n = 1342) | 69.21% (n = 717) | 77.26% (n = 625) |
No insomnia with short sleep on weekdays (<6 h) | 6.39% (n = 118) | 6.37% (n = 66) | 6.43% (n = 52) |
Insomnia | 20.83% (n = 391) | 24.29% (n = 257) | 16.36% (n = 134) |
Insomnia with short sleep on weekdays (<6 h) | 6.50% (n = 120) | 7.14% (n = 74) | 5.69% (n = 46) |
CES-D = Center for Epidemiological Studies-Depression scale.
Descriptive statistics and frequencies for all variables for each onset of regular cannabis use group
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Onset of regular cannabis use . | . | . |
. | 17 and younger (n = 300) . | After 17 (n = 246) . | No onset (n = 1336) . |
Sleep duration (h) | 6.99 (1.45) | 7.21 (1.41) | 7.07 (1.26) |
Current depression symptoms (CES-D) | 12.07 (8.87) | 12.17 (9.49) | 10.58 (8.35) |
Frequencies | |||
Onset of regular cannabis use | |||
17 and younger | After 17 | No onset | |
Prior anxiety diagnosis | 13.67% (n = 41) | 12.50% (n = 30) | 7.52% (n = 100) |
Prior depression diagnosis | 21.67% (n = 65) | 21.07% (n = 51) | 12.81% (n = 170) |
No insomnia and sleep ≥6 h on weekdays | 62.46% (n = 183) | 74.17% (n = 178) | 74.77% (n = 981) |
No insomnia with short sleep on weekdays (<6 h) | 8.53% (n = 25) | 4.58% (n = 11) | 6.25% (n = 82) |
Insomnia | 28.96% (n = 86) | 20.82% (n = 51) | 19.03% (n = 254) |
Insomnia with short sleep on weekdays (<6 h) | 9.89% (n = 29) | 7.92% (n = 19) | 5.49% (n = 72) |
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Onset of regular cannabis use . | . | . |
. | 17 and younger (n = 300) . | After 17 (n = 246) . | No onset (n = 1336) . |
Sleep duration (h) | 6.99 (1.45) | 7.21 (1.41) | 7.07 (1.26) |
Current depression symptoms (CES-D) | 12.07 (8.87) | 12.17 (9.49) | 10.58 (8.35) |
Frequencies | |||
Onset of regular cannabis use | |||
17 and younger | After 17 | No onset | |
Prior anxiety diagnosis | 13.67% (n = 41) | 12.50% (n = 30) | 7.52% (n = 100) |
Prior depression diagnosis | 21.67% (n = 65) | 21.07% (n = 51) | 12.81% (n = 170) |
No insomnia and sleep ≥6 h on weekdays | 62.46% (n = 183) | 74.17% (n = 178) | 74.77% (n = 981) |
No insomnia with short sleep on weekdays (<6 h) | 8.53% (n = 25) | 4.58% (n = 11) | 6.25% (n = 82) |
Insomnia | 28.96% (n = 86) | 20.82% (n = 51) | 19.03% (n = 254) |
Insomnia with short sleep on weekdays (<6 h) | 9.89% (n = 29) | 7.92% (n = 19) | 5.49% (n = 72) |
CES-D = Center for Epidemiological Studies-Depression scale.
Descriptive statistics and frequencies for all variables for each onset of regular cannabis use group
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Onset of regular cannabis use . | . | . |
. | 17 and younger (n = 300) . | After 17 (n = 246) . | No onset (n = 1336) . |
Sleep duration (h) | 6.99 (1.45) | 7.21 (1.41) | 7.07 (1.26) |
Current depression symptoms (CES-D) | 12.07 (8.87) | 12.17 (9.49) | 10.58 (8.35) |
Frequencies | |||
Onset of regular cannabis use | |||
17 and younger | After 17 | No onset | |
Prior anxiety diagnosis | 13.67% (n = 41) | 12.50% (n = 30) | 7.52% (n = 100) |
Prior depression diagnosis | 21.67% (n = 65) | 21.07% (n = 51) | 12.81% (n = 170) |
No insomnia and sleep ≥6 h on weekdays | 62.46% (n = 183) | 74.17% (n = 178) | 74.77% (n = 981) |
No insomnia with short sleep on weekdays (<6 h) | 8.53% (n = 25) | 4.58% (n = 11) | 6.25% (n = 82) |
Insomnia | 28.96% (n = 86) | 20.82% (n = 51) | 19.03% (n = 254) |
Insomnia with short sleep on weekdays (<6 h) | 9.89% (n = 29) | 7.92% (n = 19) | 5.49% (n = 72) |
Mean (SD) sample characteristics . | . | . | . |
---|---|---|---|
. | Onset of regular cannabis use . | . | . |
. | 17 and younger (n = 300) . | After 17 (n = 246) . | No onset (n = 1336) . |
Sleep duration (h) | 6.99 (1.45) | 7.21 (1.41) | 7.07 (1.26) |
Current depression symptoms (CES-D) | 12.07 (8.87) | 12.17 (9.49) | 10.58 (8.35) |
Frequencies | |||
Onset of regular cannabis use | |||
17 and younger | After 17 | No onset | |
Prior anxiety diagnosis | 13.67% (n = 41) | 12.50% (n = 30) | 7.52% (n = 100) |
Prior depression diagnosis | 21.67% (n = 65) | 21.07% (n = 51) | 12.81% (n = 170) |
No insomnia and sleep ≥6 h on weekdays | 62.46% (n = 183) | 74.17% (n = 178) | 74.77% (n = 981) |
No insomnia with short sleep on weekdays (<6 h) | 8.53% (n = 25) | 4.58% (n = 11) | 6.25% (n = 82) |
Insomnia | 28.96% (n = 86) | 20.82% (n = 51) | 19.03% (n = 254) |
Insomnia with short sleep on weekdays (<6 h) | 9.89% (n = 29) | 7.92% (n = 19) | 5.49% (n = 72) |
CES-D = Center for Epidemiological Studies-Depression scale.
Coding onset as an ordinal variable and controlling for sex, onset of regular cannabis use significantly predicted insomnia (standardized β = 0.14, p = <0.001) and insomnia with short sleep on weekdays (standardized β = 0.14, p = <0.001). Controlling for sex, depressive symptoms, and prior diagnoses of an anxiety or depression disorder, onset of regular cannabis use significantly predicted insomnia (standardized β = 0.07, p = 0.024) and insomnia with short sleep on weekdays (standardized β = 0.08, p = 0.003). Lastly, in a model including all covariates and shift-work, onset of regular cannabis use still significantly predicted insomnia (standardized β = 0.07, p = 0.025) and insomnia with short sleep on weekdays (standardized β = 0.08, p = 0.003), suggesting that shift-work did not explain the association of onset of early regular cannabis use with insomnia and insomnia with short sleep on weekdays.
Genetic analyses
Cross-twin and cross-twin cross-trait correlations.
We found higher cross-twin correlations among MZ twin pairs than DZ twin pairs for all traits (see Table 3), indicating genetic influences on onset of regular cannabis use, insomnia, and insomnia with short sleep on weekdays. Table 3 also displays cross-twin cross-trait correlations among MZ twin pairs and DZ twin pairs, suggesting common genetic influences between the onset of regular cannabis use and both insomnia and insomnia with short sleep on weekdays.
Cross-twin and cross-twin cross-trait correlations for all main study variables
Zygosity . | Onset of regular cannabis use . | Insomnia . | Insomnia with short sleep on weekdays (<6 h) . | Onset and insomnia . | Onset and insomnia with short sleep on weekdays (<6 h) . |
---|---|---|---|---|---|
MZ | 0.74* (p < 0.001) | 0.46* (p < 0.001) | 0.39* (p < 0.001) | 0.11 (p = 0.102) | 0.13* (p = 0.019) |
DZ | 0.49* (p < 0.001) | 0.16 (p = 0.086) | 0.15 (p = 0.134) | 0.07 (p = 0.227) | 0.06 (p = 0.220) |
Zygosity . | Onset of regular cannabis use . | Insomnia . | Insomnia with short sleep on weekdays (<6 h) . | Onset and insomnia . | Onset and insomnia with short sleep on weekdays (<6 h) . |
---|---|---|---|---|---|
MZ | 0.74* (p < 0.001) | 0.46* (p < 0.001) | 0.39* (p < 0.001) | 0.11 (p = 0.102) | 0.13* (p = 0.019) |
DZ | 0.49* (p < 0.001) | 0.16 (p = 0.086) | 0.15 (p = 0.134) | 0.07 (p = 0.227) | 0.06 (p = 0.220) |
MZ = monozygotic; DZ = dizygotic.
* p < 0.05.
Cross-twin and cross-twin cross-trait correlations for all main study variables
Zygosity . | Onset of regular cannabis use . | Insomnia . | Insomnia with short sleep on weekdays (<6 h) . | Onset and insomnia . | Onset and insomnia with short sleep on weekdays (<6 h) . |
---|---|---|---|---|---|
MZ | 0.74* (p < 0.001) | 0.46* (p < 0.001) | 0.39* (p < 0.001) | 0.11 (p = 0.102) | 0.13* (p = 0.019) |
DZ | 0.49* (p < 0.001) | 0.16 (p = 0.086) | 0.15 (p = 0.134) | 0.07 (p = 0.227) | 0.06 (p = 0.220) |
Zygosity . | Onset of regular cannabis use . | Insomnia . | Insomnia with short sleep on weekdays (<6 h) . | Onset and insomnia . | Onset and insomnia with short sleep on weekdays (<6 h) . |
---|---|---|---|---|---|
MZ | 0.74* (p < 0.001) | 0.46* (p < 0.001) | 0.39* (p < 0.001) | 0.11 (p = 0.102) | 0.13* (p = 0.019) |
DZ | 0.49* (p < 0.001) | 0.16 (p = 0.086) | 0.15 (p = 0.134) | 0.07 (p = 0.227) | 0.06 (p = 0.220) |
MZ = monozygotic; DZ = dizygotic.
* p < 0.05.
Univariate analyses.
To determine the genetic/environmental contributions to the variance of each individual trait we conducted univariate twin models. Table 4 reports the model fit statistics and model comparisons for each individual trait. We found significant sex differences in thresholds for onset of regular cannabis use (χ 2diff(2) = 52.227, p < 0.001), insomnia (χ 2diff(1) = 15.216, p < 0.001), and insomnia with short sleep on weekdays (χ 2diff(3) = 23.295, p < 0.001). We did not find scalar sex differences in the variance components for any traits (all χ 2diff(3) < 0.933, p > 0.817) suggesting no significant differences in the genetic and environmental contributions for these traits. Therefore, our models allowed separate thresholds for each sex, but equated A, C, and E parameters.
Fit and comparison tests of univariate twin models for onset of regular cannabis use, insomnia, and insomnia with short sleep on weekdays (<6 h)
Measure . | Model . | Model fit . | . | . | . | . | Standardized paths . | . | . | χ 2 test with full ACE model . | . | . |
---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | A . | C . | E . | df . | χ 2 . | P . |
Onset of regular cannabis use | ACE | 11.332 | 22 | 0.9697 | 0.000 | 1 | 0.708 | 0.486 | 0.512 | – | – | – |
AE | 14.557 | 23 | 0.9098 | 0.000 | 1 | 0.872 | – | 0.489 | 1 | 3.747 | 0.0529 | |
CE | 22.304 | 23 | 0.5019 | 0.000 | 1 | – | 0.814 | 0.580 | 1 | 13.060 | <0.001 | |
AC | 50.879 | 23 | 0.0007 | 0.080 | 0.937 | 1 | 0 | – | 1 | 47.438 | <0.001 | |
Insomnia | ACE | 3.553 | 12 | 0.9902 | 0.000 | 1 | 0.664 | 0 | 0.748 | – | – | – |
AE | 3.554 | 13 | 0.9951 | 0.000 | 1 | 0.664 | – | 0.748 | 1 | 0.00 | 0.9986 | |
CE | 9.706 | 13 | 0.7178 | 0.000 | 1 | – | 0.583 | 0.813 | 1 | 6.145 | 0.0132 | |
AC | 31.932 | 13 | 0.0025 | 0.088 | 0.497 | 1 | 0 | – | 1 | 28.351 | <0.001 | |
Insomnia with short sleep on weekdays (<6 h) | ACE | 17.461 | 32 | 0.9827 | 0.000 | 1 | 0.610 | 0 | 0.793 | – | – | – |
AE | 17.587 | 33 | 0.9871 | 0.000 | 1 | 0.610 | – | 0.793 | 1 | 0.00 | 1.00 | |
CE | 22.182 | 33 | 0.9235 | 0.000 | 1 | – | 0.532 | 0.847 | 1 | 5.802 | 0.0160 | |
AC | 109.127 | 33 | 0.00 | 0.111 | 0 | 1 | 0 | – | 1 | 115.650 | <0.001 |
Measure . | Model . | Model fit . | . | . | . | . | Standardized paths . | . | . | χ 2 test with full ACE model . | . | . |
---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | A . | C . | E . | df . | χ 2 . | P . |
Onset of regular cannabis use | ACE | 11.332 | 22 | 0.9697 | 0.000 | 1 | 0.708 | 0.486 | 0.512 | – | – | – |
AE | 14.557 | 23 | 0.9098 | 0.000 | 1 | 0.872 | – | 0.489 | 1 | 3.747 | 0.0529 | |
CE | 22.304 | 23 | 0.5019 | 0.000 | 1 | – | 0.814 | 0.580 | 1 | 13.060 | <0.001 | |
AC | 50.879 | 23 | 0.0007 | 0.080 | 0.937 | 1 | 0 | – | 1 | 47.438 | <0.001 | |
Insomnia | ACE | 3.553 | 12 | 0.9902 | 0.000 | 1 | 0.664 | 0 | 0.748 | – | – | – |
AE | 3.554 | 13 | 0.9951 | 0.000 | 1 | 0.664 | – | 0.748 | 1 | 0.00 | 0.9986 | |
CE | 9.706 | 13 | 0.7178 | 0.000 | 1 | – | 0.583 | 0.813 | 1 | 6.145 | 0.0132 | |
AC | 31.932 | 13 | 0.0025 | 0.088 | 0.497 | 1 | 0 | – | 1 | 28.351 | <0.001 | |
Insomnia with short sleep on weekdays (<6 h) | ACE | 17.461 | 32 | 0.9827 | 0.000 | 1 | 0.610 | 0 | 0.793 | – | – | – |
AE | 17.587 | 33 | 0.9871 | 0.000 | 1 | 0.610 | – | 0.793 | 1 | 0.00 | 1.00 | |
CE | 22.182 | 33 | 0.9235 | 0.000 | 1 | – | 0.532 | 0.847 | 1 | 5.802 | 0.0160 | |
AC | 109.127 | 33 | 0.00 | 0.111 | 0 | 1 | 0 | – | 1 | 115.650 | <0.001 |
Models in bold indicate the best fitting model based on χ 2 difference tests.
Fit and comparison tests of univariate twin models for onset of regular cannabis use, insomnia, and insomnia with short sleep on weekdays (<6 h)
Measure . | Model . | Model fit . | . | . | . | . | Standardized paths . | . | . | χ 2 test with full ACE model . | . | . |
---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | A . | C . | E . | df . | χ 2 . | P . |
Onset of regular cannabis use | ACE | 11.332 | 22 | 0.9697 | 0.000 | 1 | 0.708 | 0.486 | 0.512 | – | – | – |
AE | 14.557 | 23 | 0.9098 | 0.000 | 1 | 0.872 | – | 0.489 | 1 | 3.747 | 0.0529 | |
CE | 22.304 | 23 | 0.5019 | 0.000 | 1 | – | 0.814 | 0.580 | 1 | 13.060 | <0.001 | |
AC | 50.879 | 23 | 0.0007 | 0.080 | 0.937 | 1 | 0 | – | 1 | 47.438 | <0.001 | |
Insomnia | ACE | 3.553 | 12 | 0.9902 | 0.000 | 1 | 0.664 | 0 | 0.748 | – | – | – |
AE | 3.554 | 13 | 0.9951 | 0.000 | 1 | 0.664 | – | 0.748 | 1 | 0.00 | 0.9986 | |
CE | 9.706 | 13 | 0.7178 | 0.000 | 1 | – | 0.583 | 0.813 | 1 | 6.145 | 0.0132 | |
AC | 31.932 | 13 | 0.0025 | 0.088 | 0.497 | 1 | 0 | – | 1 | 28.351 | <0.001 | |
Insomnia with short sleep on weekdays (<6 h) | ACE | 17.461 | 32 | 0.9827 | 0.000 | 1 | 0.610 | 0 | 0.793 | – | – | – |
AE | 17.587 | 33 | 0.9871 | 0.000 | 1 | 0.610 | – | 0.793 | 1 | 0.00 | 1.00 | |
CE | 22.182 | 33 | 0.9235 | 0.000 | 1 | – | 0.532 | 0.847 | 1 | 5.802 | 0.0160 | |
AC | 109.127 | 33 | 0.00 | 0.111 | 0 | 1 | 0 | – | 1 | 115.650 | <0.001 |
Measure . | Model . | Model fit . | . | . | . | . | Standardized paths . | . | . | χ 2 test with full ACE model . | . | . |
---|---|---|---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | A . | C . | E . | df . | χ 2 . | P . |
Onset of regular cannabis use | ACE | 11.332 | 22 | 0.9697 | 0.000 | 1 | 0.708 | 0.486 | 0.512 | – | – | – |
AE | 14.557 | 23 | 0.9098 | 0.000 | 1 | 0.872 | – | 0.489 | 1 | 3.747 | 0.0529 | |
CE | 22.304 | 23 | 0.5019 | 0.000 | 1 | – | 0.814 | 0.580 | 1 | 13.060 | <0.001 | |
AC | 50.879 | 23 | 0.0007 | 0.080 | 0.937 | 1 | 0 | – | 1 | 47.438 | <0.001 | |
Insomnia | ACE | 3.553 | 12 | 0.9902 | 0.000 | 1 | 0.664 | 0 | 0.748 | – | – | – |
AE | 3.554 | 13 | 0.9951 | 0.000 | 1 | 0.664 | – | 0.748 | 1 | 0.00 | 0.9986 | |
CE | 9.706 | 13 | 0.7178 | 0.000 | 1 | – | 0.583 | 0.813 | 1 | 6.145 | 0.0132 | |
AC | 31.932 | 13 | 0.0025 | 0.088 | 0.497 | 1 | 0 | – | 1 | 28.351 | <0.001 | |
Insomnia with short sleep on weekdays (<6 h) | ACE | 17.461 | 32 | 0.9827 | 0.000 | 1 | 0.610 | 0 | 0.793 | – | – | – |
AE | 17.587 | 33 | 0.9871 | 0.000 | 1 | 0.610 | – | 0.793 | 1 | 0.00 | 1.00 | |
CE | 22.182 | 33 | 0.9235 | 0.000 | 1 | – | 0.532 | 0.847 | 1 | 5.802 | 0.0160 | |
AC | 109.127 | 33 | 0.00 | 0.111 | 0 | 1 | 0 | – | 1 | 115.650 | <0.001 |
Models in bold indicate the best fitting model based on χ 2 difference tests.
We then tested each etiological pathway with a χ 2 difference test to determine the most parsimonious model. Shared environmental influences were not significantly different from zero for all univariate models (all χ 2diff(1) < 3.747, p > 0.0529), but model comparisons determined A and E could not be dropped (all χ 2diff(1) > 5.825, p < 0.0160). In the most parsimonious A and E models, additive genetic and nonshared environmental factors contributed to the onset of regular cannabis use (a2 = 76% and e2 = 24%, all χ 2diff(1) > 43.691, p ≤ 0.001), insomnia (a2 = 44% and e2 = 56%, all χ 2diff(1) > 40.920, p = <0.001), and insomnia with short sleep on weekdays (a2 = 37% and e2 = 63%, all χ 2diff(1) > 40.479, p ≤ 0.001).
Bivariate twin models.
We conducted separate bivariate Cholesky decompositions for onset of regular cannabis with insomnia and onset of regular cannabis use with insomnia with short sleep on weekdays. We found no significant scalar sex differences for either of our bivariate models (all χ 2diff(9) < 4.583, p > 0.869), suggesting no genetic and/or environmental sex differences. Table 5 includes model fits for bivariate twin models without sex differences. We found that C paths could be dropped from both models without a significant decrement in fit (all χ 2diff(3) < 3.773, p > 0.287), similar to our univariate results.
Fit and comparison tests of bivariate twin models for onset of regular cannabis use with both insomnia and insomnia with short sleep on weekdays (<6 h)
Measure . | Model . | Model fit . | . | . | . | . | χ 2 test with full Cholesky ACE decomposition . | . | . |
---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | df . | χ 2 . | P . |
Onset of regular cannabis use and insomnia | Full Cholesky ACE | 40.923 | 51 | 0.8426 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 58.191 | 54 | 0.3238 | 0.020 | 0.991 | 3 | 18.375 | <0.001 | |
Dropping all C paths | 44.599 | 54 | 0.8154 | 0.000 | 1.000 | 3 | 3.736 | 0.2914 | |
Dropping all E paths | 142.629 | 54 | 0.0000 | 0.093 | 0.809 | 3 | 125.396 | <0.001 | |
Onset of regular cannabis use and insomnia with short sleep on weekdays (<6 h) | Full Cholesky ACE | 54.866 | 71 | 0.9215 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 71.443 | 74 | 0.5626 | 0.000 | 1.000 | 3 | 19.147 | <0.001 | |
Dropping all C paths | 58.465 | 74 | 0.9071 | 0.000 | 1.000 | 3 | 3.772 | 0.2872 | |
Dropping all E paths | 189.233 | 74 | 0.0000 | 0.091 | 0.758 | 3 | 182.649 | <0.001 |
Measure . | Model . | Model fit . | . | . | . | . | χ 2 test with full Cholesky ACE decomposition . | . | . |
---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | df . | χ 2 . | P . |
Onset of regular cannabis use and insomnia | Full Cholesky ACE | 40.923 | 51 | 0.8426 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 58.191 | 54 | 0.3238 | 0.020 | 0.991 | 3 | 18.375 | <0.001 | |
Dropping all C paths | 44.599 | 54 | 0.8154 | 0.000 | 1.000 | 3 | 3.736 | 0.2914 | |
Dropping all E paths | 142.629 | 54 | 0.0000 | 0.093 | 0.809 | 3 | 125.396 | <0.001 | |
Onset of regular cannabis use and insomnia with short sleep on weekdays (<6 h) | Full Cholesky ACE | 54.866 | 71 | 0.9215 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 71.443 | 74 | 0.5626 | 0.000 | 1.000 | 3 | 19.147 | <0.001 | |
Dropping all C paths | 58.465 | 74 | 0.9071 | 0.000 | 1.000 | 3 | 3.772 | 0.2872 | |
Dropping all E paths | 189.233 | 74 | 0.0000 | 0.091 | 0.758 | 3 | 182.649 | <0.001 |
Models in bold indicate the best fitting model via χ 2 difference tests.
Fit and comparison tests of bivariate twin models for onset of regular cannabis use with both insomnia and insomnia with short sleep on weekdays (<6 h)
Measure . | Model . | Model fit . | . | . | . | . | χ 2 test with full Cholesky ACE decomposition . | . | . |
---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | df . | χ 2 . | P . |
Onset of regular cannabis use and insomnia | Full Cholesky ACE | 40.923 | 51 | 0.8426 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 58.191 | 54 | 0.3238 | 0.020 | 0.991 | 3 | 18.375 | <0.001 | |
Dropping all C paths | 44.599 | 54 | 0.8154 | 0.000 | 1.000 | 3 | 3.736 | 0.2914 | |
Dropping all E paths | 142.629 | 54 | 0.0000 | 0.093 | 0.809 | 3 | 125.396 | <0.001 | |
Onset of regular cannabis use and insomnia with short sleep on weekdays (<6 h) | Full Cholesky ACE | 54.866 | 71 | 0.9215 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 71.443 | 74 | 0.5626 | 0.000 | 1.000 | 3 | 19.147 | <0.001 | |
Dropping all C paths | 58.465 | 74 | 0.9071 | 0.000 | 1.000 | 3 | 3.772 | 0.2872 | |
Dropping all E paths | 189.233 | 74 | 0.0000 | 0.091 | 0.758 | 3 | 182.649 | <0.001 |
Measure . | Model . | Model fit . | . | . | . | . | χ 2 test with full Cholesky ACE decomposition . | . | . |
---|---|---|---|---|---|---|---|---|---|
. | . | χ 2 . | df . | P . | RMSEA . | CFI . | df . | χ 2 . | P . |
Onset of regular cannabis use and insomnia | Full Cholesky ACE | 40.923 | 51 | 0.8426 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 58.191 | 54 | 0.3238 | 0.020 | 0.991 | 3 | 18.375 | <0.001 | |
Dropping all C paths | 44.599 | 54 | 0.8154 | 0.000 | 1.000 | 3 | 3.736 | 0.2914 | |
Dropping all E paths | 142.629 | 54 | 0.0000 | 0.093 | 0.809 | 3 | 125.396 | <0.001 | |
Onset of regular cannabis use and insomnia with short sleep on weekdays (<6 h) | Full Cholesky ACE | 54.866 | 71 | 0.9215 | 0.000 | 1.000 | – | – | – |
Dropping all A paths | 71.443 | 74 | 0.5626 | 0.000 | 1.000 | 3 | 19.147 | <0.001 | |
Dropping all C paths | 58.465 | 74 | 0.9071 | 0.000 | 1.000 | 3 | 3.772 | 0.2872 | |
Dropping all E paths | 189.233 | 74 | 0.0000 | 0.091 | 0.758 | 3 | 182.649 | <0.001 |
Models in bold indicate the best fitting model via χ 2 difference tests.
Figure 1 illustrates the best fitting A and E Cholesky decompositions and their genetic/environmental cross paths. The A and E Cholesky decomposition with the best fit between onset of regular cannabis use and insomnia indicated a significant genetic overlap (a cross-path = 0.13, rA = 0.20, χ 2diff(1) = 3.934, p = 0.047) between the traits. The best fitting A and E Cholesky decomposition between onset of regular cannabis use and insomnia with short sleep on weekdays also suggested a significant genetic overlap (a cross-path = 0.15, rA = 0.25, χ 2diff(1) = 6.942, p = 0.008). Neither of these Cholesky decompositions suggested significant environmental overlap (all χ 2diff(1) < 2.300, p > 0.129).

Best fitting bivariate AE Cholesky decompositions between onset of regular cannabis use and both insomnia (Panel A) and insomnia with short sleep on weekdays (<6 h) (Panel B). Solid lines with arrows represent non-standardized additive genetic (A) and unique environmental (E) variance paths for each univariate trait. Dashed lines with arrows represent additive genetic (A) and unique environmental (E) cross paths between onset of regular cannabis use and insomnia/insomnia with short sleep on weekdays (<6 h). Parameters that are significant (p < 0.05) are indicated via asterisks (all parameters significant besides the cross E paths in both decompositions).
Discussion
We tested the hypothesis that early onset of regular cannabis use predicted higher rates of young adult insomnia. Controlling for sex, current depression symptoms, and prior anxiety and depression diagnoses, we found that early onset of regular cannabis use was associated with increased rates of young adult insomnia and insomnia with short sleep on weekdays. Using a genetically informative twin design, we found both genetic and nonshared environmental contributions to the etiology of age of regular cannabis use onset, insomnia, and insomnia with short sleep on weekdays individually. These findings of significant genetic contribution are consistent with prior research on cannabis initiation [73] and insomnia [74] but this is the first report of genetic contributions to insomnia with short sleep on weekdays. Lastly, we found evidence of significant overlapping additive genetic influences on onset of regular cannabis use and both insomnia and insomnia with short sleep on weekdays, implying common genetic liabilities, but we did not find significant environmental overlap.
Our genetic and environmental variance estimates for insomnia resemble a prior twin study that utilized the same cohort [75], but our additional results on insomnia with short sleep on weekdays are novel. We found that additive genetics contributed more to insomnia alone, and that unique environment contributed more to insomnia with short sleep on weekdays. The higher environmental contribution to insomnia with short sleep on weekdays may reflect the fact that the unique environment largely contributes to the majority of the variance in adult sleep duration (e2 = 67–80% and a2 = 20–32%) [18, 76, 77]. Thus, inclusion of the short sleep with insomnia might lead the trait to have a larger contribution from unique environment. Our findings of genetic contribution to the onset of regular cannabis use are consistent with studies of lifetime cannabis use and problem cannabis use [73], and align with a prior twin study focused on the onset of regular cannabis use and adult sleep duration that used a subsample of the current study’s sample [18]. Although our best fitting and most parsimonious model for onset of regular cannabis use consisted of A and E parameter, the presence of some C (though only marginally significant) is consistent with previous studies of substance use [78].
These findings are consistent with prior studies [17, 18] and build on the current body of research regarding the relationship between early cannabis use and adult sleep problems. We find early regular cannabis use predicts higher rates of insomnia and insomnia with short sleep on weekdays in young adulthood. Thus, our findings provide evidence in support of the theory that early regular cannabis use is associated with subsequent sleep problems, extending this effect to clinical sleep outcomes. The differing time frames of the measures strengthen the support for possible long-lasting effects of early cannabis on sleep problems in young adulthood and our study reinforces the evidence of possible pleotropic influence of genes on early cannabis use and sleep [18].
Several possible behavioral, developmental, or neurological mechanisms have previously been proposed that could explain how early regular cannabis use might influence later sleep and insomnia in young adulthood. For example, early cannabis users could develop maladaptive behaviors or long-lasting sleep problems that could progress into adulthood affecting sleep [35]. Additionally, circadian pathway disturbances resulting from substance use [79–81] and early cannabis use could lead to developmental disruptions in circadian rhythms and late sleep timing that could advance into adulthood. Lastly, early cannabis use could result in altered brain development that affects adult insomnia [18], as there is evidence of alterations of the prefrontal cortex (PFC) for both traits. Functional and structural imaging studies have shown early cannabis use is associated with alterations in the PFC and early users suffer impairments in neurocognitive performance associated with those area [82]. Correspondingly, insomnia is associated with deficits in neurocognitive performance in executive function tasks, suggesting PFC impairments [56, 57, 83] and there is evidence of insomnia being related to structural alterations in regions of the PFC [84–86], changes in functional PFC connectivity as well as in networks involving the PFC [87–92], and variations in EEG sleep patterns in the PFC [93, 94].
However, our results suggest that the relationship of onset of regular cannabis use with insomnia and insomnia with short sleep on weekdays may also be due to pleiotropic influence of genes on these traits. A recent study using similar twin models found that shared genetics may play a role in the etiology of the relationship between early cannabis use and adult sleep duration [18], and there are various additional lines of evidence consistent with this common genetic liability theory. Endocannabinoids may play a large part in this relationship, as the endocannabinoid system may be involved in the circadian sleep–wake cycle [43–46]. Endocannabinoids influence sleep, and their levels vary with time of day [45]. Additionally, several large genome-wide association studies (GWASs) of sleep-related variables, including chronotype [95], sleep duration [96, 97], and insomnia [98] have found significant genes and genetic pathways linked with cannabis use [40, 99] or cannabinoid activity [100–106].
Genes related to circadian rhythm and sleep could also play an important role in this relationship. Several GWASs of lifetime cannabis use [39–41] have found significant genes associated that are believed to be involved in circadian rhythm and sleep behaviors [107–111], and a recent study found that several clock gene polymorphisms were significant risk factors for cannabis addiction [42]. Despite our best fitting model implying that shared genetics explain the association between these traits, it is entirely plausible there is a causal link between cannabis use and sleep deficits; a previous analysis found that a causal model between early regular cannabis use and sleep duration provided the best fit (although there are limitations to consider with that design) [18]. GWAS methods such as Mendelian randomization [112] can provide further information on causal mechanisms between the genetics of these traits. Further research should consider the role of genetics in the relationship between early regular cannabis use and later sleep (specifically endocannabinoid and clock/circadian/sleep genes), with a focus on clinically important outcomes such as insomnia.
Limitations
There are several limitations worth addressing. First, our sleep duration measure was self-reported and only a small collection of insomnia with short sleep studies has used self-reported measures [54, 113]. The reliability of self-reported reports, and their consistency with more commonly used objective measures like polysomnography (PSG) and actigraphy [114] are not perfect. But there are also critiques of objective laboratory sleep assessments, as they are commonly limited to one night or a few nights of recording and could lack accuracy with respect to at-home sleep variability [115]. Criticisms of actigraphy include its inability to distinguish sleep disorders [116] and its tendency to overestimate sleep times [117–119]. A recent comprehensive review of insomnia with short sleep called for research with self-reported as well as objective measures [114], and research on self-reported sleep duration measures is clearly needed as an alternative when objective measures are too costly and time intensive. Development of accurate, valid, and easy to use wearable technology that can better assess sleep in the home environment is needed. Our definition of short sleep (<6 h) with insomnia is consistent with the existing literature, but it should be noted that more than 7 h sleep is recommended for promoting health in adults [120], thus future research should explore other cutoffs of short sleep with insomnia.
Second, the lower bound for the age range at the time of the substance use assessment was 18, which is very close to the cut-off age for early regular cannabis use and perhaps does not allow enough time to pass to develop insomnia symptomology. Third, our measures were collected via self-report procedures and our onset of regular cannabis use measure was a retrospective question which could be prone to memory discrepancies or report bias [121]. Lastly, due to time interval between the substance use interview and sleep questionnaire, we could not separately assess the contributions of current cannabis or other substance use to insomnia as we did not have frequency of use measures collected at the time of administration of the sleep questionnaire. Lack of concurrent substance use measures is a limitation of the current study, as use of other substances including alcohol [122, 123] and tobacco [124–126], as well as cannabis [15], have been linked to sleep deficits and insomnia. Future studies should consider and include appropriate substance use covariates as well as objective measures of sleep to address the mentioned concerns.
Summary
Our findings are consistent with the theory that early cannabis use is associated with increased rates of young adult insomnia and insomnia with short sleep on weekdays. These results extend the current body of research regarding the relationship of early cannabis use with adult sleep to include clinical sleep outcomes as well as provide novel insight into shared genetics between these traits. This is also the first twin study to estimate, and find, genetic influences on insomnia with short sleep on weekdays, which exhibited similar etiology to other sleep traits. Further research is needed to understand the potential developmental impact of early regular cannabis use on adult sleep, specifically on clinically important outcomes such as insomnia. Future studies should focus on genetic factors and would benefit from using both self-reported and object measurers of sleep duration.
Ethical Approval
Documented consent was obtained from all study participants and all procedures of this study followed the ethical standards of the University of Colorado Institutional Review Board.
Funding
Supported by grants P60 DA011015, R01 DA042755, T32 DA017637, and funding from the Henry Ford Health System.
Conflict of interest statement: This study was not financially supported by industry. Christopher L. Drake has received grant support from Cephalon, Zeo Inc., Merck Inc., and Takeda Inc. and participated in speaking engagements for Cephalon and Asante Communications all unrelated to this study. Kenneth P. Wright Jr. has received research support from the National Institutes of Health, Office of Naval Research, Pac-12, and SomaLogic, Inc., consulting fees from or served as a paid member of scientific advisory boards for the Sleep Disorders Research Advisory Board—National Heart, Lung and Blood Institute, CurAegis Technologies, Circadian Therapeutics, Ltd and has received speaker/educational/travel consultant honorarium fees from the American Academy of Sleep Medicine, American College of Chest Physicians, American College of Sports Medicine, American Diabetes Association, Associated Professional Sleep Societies, Kellogg Company, Obesity Medicine Association, and The European Association for the Study of Obesity. All other authors have no potential conflicts of interest to disclose.
Acknowledgments
We would like to thank all authors for their insightful comments and suggestions regarding the design, analysis, and wording of the manuscript.
References
Comments