Trans-ancestry epigenome-wide association meta-analysis of DNA methylation with lifetime cannabis use

Cannabis is widely used worldwide, yet its links to health outcomes are not fully understood. DNA methylation can serve as a mediator to link environmental exposures to health outcomes. We conducted an epigenome-wide association study (EWAS) of peripheral blood-based DNA methylation and lifetime cannabis use (ever vs. never) in a meta-analysis including 9436 participants (7795 European and 1641 African ancestry) from seven cohorts. Accounting for effects of cigarette smoking, our trans-ancestry EWAS meta-analysis revealed four CpG sites significantly associated with lifetime cannabis use at a false discovery rate of 0.05 ( p < 5.85 × 10 - 7 ) : cg22572071 near gene ADGRF1, cg15280358 in ADAM12, cg00813162 in ACTN1, and cg01101459 near LINC01132. Additionally, our EWAS analysis in participants who never smoked cigarettes identified another epigenome-wide significant CpG site, cg14237301 annotated to APOBR. We used a leave-one-out approach to evaluate methylation scores constructed as a weighted sum of the significant CpGs. The best model can explain 3.79% of the variance in lifetime cannabis use. These findings unravel the DNA methylation changes associated with lifetime cannabis use that are independent of cigarette smoking and may serve as a starting point for further research on the mechanisms through which cannabis exposure impacts health outcomes.

Associations of polygenic risk scores for smoking heaviness and lifetime cannabis use with tobacco and cannabis co-use trajectories among African Americans

BACKGROUND

We aimed to identify distinct trajectories of tobacco, cannabis, and their co-use among African Americans, and to investigate whether these patterns were associated with polygenic risk scores (PRS) for tobacco and cannabis use.

METHOD

Participants (N=428 participants; 50.9% male) were initially recruited for an elementary school-based prevention in a Mid-Atlantic city when they were in first grade. From ages 14-26, participants reported on their frequency of tobacco and cannabis use in the past year during annual assessments. DNA was collected from participants at age 21. PRS for smoking heaviness (i.e., cigarettes per day) and lifetime cannabis use were created based on genome-wide association study results derived from Liu et al. (2019) and Pasman et al. (2018), respectively.

RESULTS

We identified five distinct trajectories of tobacco and cannabis co-use, including (1) Low Tobacco and Cannabis Use, (2) Adolescent Limited Tobacco and Cannabis Use, (3) Experimental Cannabis, Young Adult Increasing Tobacco, (4) Experimental Tobacco, Young Adult Increasing Cannabis, and (5) High, Chronic Tobacco and Cannabis Use. Compared to the Low Tobacco and Cannabis Use subgroup, individuals in the High, Chronic Tobacco and Cannabis Use subgroup had greater PRS for smoking heaviness, and individuals in the Experimental Cannabis, Young Adult Increasing Tobacco subgroup had higher PRS for lifetime cannabis use.

CONCLUSIONS

Polygenic risk for lifetime cannabis use and smoking heaviness is associated with the developmental progression of tobacco and cannabis co-use among African Americans, furthering knowledge on the etiology of co-use in this population.

Genome-wide DNA methylation association study of recent and cumulative marijuana use in middle aged adults

Marijuana is a widely used psychoactive substance in the US and medical and recreational legalization has risen over the past decade. Despite the growing number of individuals using marijuana, studies investigating the association between epigenetic factors and recent and cumulative marijuana use remain limited. We therefore investigated the association between recent and cumulative marijuana use and DNA methylation levels. Participants from the Coronary Artery Risk Development in Young Adults Study with whole blood collected at examination years (Y) 15 and Y20 were randomly selected to undergo DNA methylation profiling at both timepoints using the Illumina MethylationEPIC BeadChip. Recent use of marijuana was queried at each examination and used to estimate cumulative marijuana use from Y0 to Y15 and Y20. At Y15 (n = 1023), we observed 22 and 31 methylation markers associated (FDR P ≤ 0.05) with recent and cumulative marijuana use and 132 and 16 methylation markers at Y20 (n = 883), respectively. We replicated 8 previously reported methylation markers associated with marijuana use. We further identified 640 cis-meQTLs and 198 DMRs associated with recent and cumulative use at Y15 and Y20. Differentially methylated genes were statistically overrepresented in pathways relating to cellular proliferation, hormone signaling, and infections as well as schizophrenia, bipolar disorder, and substance-related disorders. We identified numerous methylation markers, pathways, and diseases associated with recent and cumulative marijuana use in middle-aged adults, providing additional insight into the association between marijuana use and the epigenome. These results provide novel insights into the role marijuana has on the epigenome and related health conditions.

Marijuana use and DNA methylation-based biological age in young adults

BACKGROUND

Marijuana is the third most commonly used drug in the USA and efforts to legalize it for medical and recreational use are growing. Despite the increase in use, marijuana's effect on aging remains understudied and understanding the effects of marijuana on molecular aging may provide novel insights into the role of marijuana in the aging process. We therefore sought to investigate the association between cumulative and recent use of marijuana with epigenetic age acceleration (EAA) as estimated from blood DNA methylation.

RESULTS

A random subset of participants from The Coronary Artery Risk Development in Young Adults (CARDIA) Study with available whole blood at examination years (Y) 15 and Y20 underwent epigenomic profiling. Four EAA estimates (intrinsic epigenetic age acceleration, extrinsic epigenetic age acceleration, PhenoAge acceleration, and GrimAge acceleration) were calculated from DNA methylation levels measured at Y15 and Y20. Ever use and cumulative marijuana-years were calculated from the baseline visit to Y15 and Y20, and recent marijuana use (both any and number of days of use in the last 30 days) were calculated at Y15 and Y20. Ever use of marijuana and each additional marijuana-year were associated with a 6-month (P < 0.001) and a 2.5-month (P < 0.001) higher average in GrimAge acceleration (GAA) using generalized estimating equations, respectively. Recent use and each additional day of recent use were associated with a 20-month (P < 0.001) and a 1-month (P < 0.001) higher GAA, respectively. A statistical interaction between marijuana-years and alcohol consumption on GAA was observed (P = 0.011), with nondrinkers exhibiting a higher GAA (β = 0.21 [95% CI 0.05, 0.36], P = 0.008) compared to heavy drinkers (β = 0.05 [95% CI - 0.09, 0.18], P = 0.500) per each additional marijuana-year. No associations were observed for the remaining EAA estimates.

CONCLUSIONS

These findings suggest cumulative and recent marijuana use are associated with age-related epigenetic changes that are related to lifespan. These observed associations may be modified by alcohol consumption. Given the increase in use and legalization, these findings provide novel insight on the effect of marijuana use on the aging process as captured through blood DNA methylation.

DNA Methylation-Based Biomarkers of Environmental Exposures for Human Population Studies

PURPOSE OF REVIEW

This manuscript orients the reader to the underlying motivations of environmental biomarker development for human population studies and provides the foundation for applying these novel biomarkers in future research. In this review, we focus our attention on the DNA methylation-based biomarkers of (i) smoking, among adults and pregnant women, (ii) lifetime cannabis use, (iii) alcohol consumption, and (iv) cumulative exposure to lead.

RECENT FINDINGS

Prior environmental exposures and lifestyle modulate DNA methylation levels. Exposure-related DNA methylation changes can either be persistent or reversible once the exposure is no longer present, and this combination of both persistent and reversible changes has essential value for biomarker development. Here, we present available biomarkers representing past and cumulative exposures using individual DNA methylation profiles. In the present work, we describe how the field of environmental epigenetics can leverage machine learning algorithms to develop exposure biomarkers and reduce problems of misreporting exposures or limited access technology. We emphasize the crucial role of the individual DNA methylation profiles in those predictions, providing a summary of each biomarker, and highlighting their advantages, and limitations. Future research can cautiously leverage these DNA methylation-based biomarkers to understand the onset and progression of diseases.

Investigating causality between liability to ADHD and substance use, and liability to substance use and ADHD risk, using Mendelian randomization

Attention-deficit hyperactivity disorder (ADHD) has consistently been associated with substance use, but the nature of this association is not fully understood. To inform intervention development and public health messages, a vital question is whether there are causal pathways from ADHD to substance use and/or vice versa. We applied bidirectional Mendelian randomization, using summary-level data from the largest available genome-wide association studies (GWAS) on ADHD, smoking (initiation, cigarettes per day, cessation, and a compound measure of lifetime smoking), alcohol use (drinks per week, alcohol problems, and alcohol dependence), cannabis use (initiation), and coffee consumption (cups per day). Genetic variants robustly associated with the "exposure" were selected as instruments and identified in the "outcome" GWAS. Effect estimates from individual genetic variants were combined with inverse-variance weighted regression and five sensitivity analyses (weighted median, weighted mode, MR-Egger, generalized summary data-based MR, and Steiger filtering). We found evidence that liability to ADHD increases likelihood of smoking initiation and heaviness of smoking among smokers, decreases likelihood of smoking cessation, and increases likelihood of cannabis initiation. There was weak evidence that liability to ADHD increases alcohol dependence risk but not drinks per week or alcohol problems. In the other direction, there was weak evidence that smoking initiation increases ADHD risk, but follow-up analyses suggested a high probability of horizontal pleiotropy. There was no clear evidence of causal pathways between ADHD and coffee consumption. Our findings corroborate epidemiological evidence, suggesting causal pathways from liability to ADHD to smoking, cannabis use, and, tentatively, alcohol dependence. Further work is needed to explore the exact mechanisms mediating these causal effects.

Investigating the causal effect of cannabis use on cognitive function with a quasi-experimental co-twin design

BACKGROUND

It is unclear whether cannabis use causes cognitive decline; several studies show an association between cannabis use and cognitive decline, but quasi-experimental twin studies have found little support for a causal effect. Here, we evaluate the association of cannabis use with general cognitive ability and executive functions (EFs) while controlling for genetic and shared environmental confounds in a longitudinal twin study.

METHODS

We first examined the phenotypic associations between cannabis initiation, frequency, and use disorder with cognitive abilities, while also controlling for pre-use general cognitive ability and other substance involvement. We tested the concurrent association between the cannabis use variables and cognitive abilities in late adolescence and young adulthood and the longitudinal association between cannabis use variables during adolescence and young adulthood cognitive abilities. Next, we used multilevel models to test whether these relations reflect between- and/or within-twin pair associations.

RESULTS

Phenotypically, cannabis use was related to poorer cognitive functioning, although most associations were negligible after accounting for other substance use. Nevertheless, there were few significant within-family twin-specific associations, except that age 17 cannabis frequency was associated with worse age 23 Common EF and general cognitive ability.

CONCLUSIONS

We found little support for a potential causal effect of cannabis use on cognition, consistent with previous twin studies. Results suggest that cannabis use may not cause decline in cognitive ability among a normative sample of cannabis users.

Psychosocial moderation of polygenic risk for cannabis involvement: the role of trauma exposure and frequency of religious service attendance

Cannabis use and disorders (CUD) are influenced by multiple genetic variants of small effect and by the psychosocial environment. However, this information has not been effectively incorporated into studies of gene-environment interaction (GxE). Polygenic risk scores (PRS) that aggregate the effects of genetic variants can aid in identifying the links between genetic risk and psychosocial factors. Using data from the Pasman et al. GWAS of cannabis use (meta-analysis of data from the International Cannabis Consortium and UK Biobank), we constructed PRS in the Collaborative Study on the Genetics of Alcoholism (COGA) participants of European (N: 7591) and African (N: 3359) ancestry. The primary analyses included only individuals of European ancestry, reflecting the ancestral composition of the discovery GWAS from which the PRS was derived. Secondary analyses included the African ancestry sample. Associations of PRS with cannabis use and DSM-5 CUD symptom count (CUDsx) and interactions with trauma exposure and frequency of religious service attendance were examined. Models were adjusted for sex, birth cohort, genotype array, and ancestry. Robustness models were adjusted for cross-term interactions. Higher PRS were associated with a greater likelihood of cannabis use and with CUDsx among participants of European ancestry (p < 0.05 and p < 0.1 thresholds, respectively). PRS only influenced cannabis use among those exposed to trauma (R2: 0.011 among the trauma exposed vs. R2: 0.002 in unexposed). PRS less consistently influenced cannabis use among those who attend religious services less frequently; PRS × religious service attendance effects were attenuated when cross-term interactions with ancestry and sex were included in the model. Polygenic liability to cannabis use was related to cannabis use and, less robustly, progression to symptoms of CUD. This study provides the first evidence of PRS × trauma for cannabis use and demonstrates that ignoring important aspects of the psychosocial environment may mask genetic influences on polygenic traits.

Systematic Review of Polygenic Gene-Environment Interaction in Tobacco, Alcohol, and Cannabis Use

Studies testing the effect of single genetic variants on substance use have had modest success. This paper reviewed 39 studies using polygenic measures to test interaction with any type of environmental exposure (G×E) in alcohol, tobacco, and cannabis use. Studies using haplotype combinations, sum scores of candidate-gene risk alleles, and polygenic scores (PS) were included. Overall study quality was moderate, with lower ratings for the polygenic methods in the haplotype and candidate-gene score studies. Heterogeneity in investigated environmental exposures, genetic factors, and outcomes was substantial. Most studies (N = 30) reported at least one significant G×E interaction, but overall evidence was weak. The majority (N = 26) found results in line with differential susceptibility and diathesis-stress frameworks. Future studies should pay more attention to methodological and statistical rigor, and focus on replication efforts. Additional work is needed before firm conclusions can be drawn about the importance of G×E in the etiology of substance use.

Exploring the relationship between polygenic risk for cannabis use, peer cannabis use and the longitudinal course of cannabis involvement

BACKGROUND AND AIMS

Few studies have explored how polygenic propensity to cannabis use unfolds across development, and no studies have yet examined this question in the context of environmental contributions such as peer cannabis use. Outlining the factors that contribute to progression from cannabis initiation to problem use over time may ultimately provide insights into mechanisms for targeted interventions. We sought to examine the relationships between polygenic liability for cannabis use, cannabis use trajectories from ages 12-30 years and perceived peer cannabis use at ages 12-17 years.

DESIGN

Mixed-effect logistic and linear regressions were used to examine associations between polygenic risk scores, cannabis use trajectory membership and perceived peer cannabis use.

SETTING

United States.

PARTICIPANTS

From the Collaborative Study on the Genetics of Alcoholism (COGA) study, a cohort of 1167 individuals aged 12-26 years at their baseline (i.e. first) interview.

MEASUREMENTS

Key measurements included life-time cannabis use (yes/no), frequency of past 12-month cannabis use, maximum life-time frequency of cannabis use, cannabis use disorder (using DSM-5 criteria) and perceived peer cannabis use. Polygenic risk scores (PRS) were created using summary statistics from a large (n = 162 082) genome-wide association study (GWAS) of cannabis use.

FINDINGS

Three trajectories reflecting no/low (n = 844), moderate (n = 137) and high (n = 186) use were identified. PRS were significantly associated with trajectory membership [P = 0.002-0.006, maximum conditional R2  = 1.4%, odds ratios (ORs) = 1.40-1.49]. Individuals who reported that most/all of their best friends used cannabis had significantly higher PRS than those who reported that none of their friends were users [OR = 1.35, 95% confidence interval (CI) = 1.04, 1.75, P = 0.023]. Perceived peer use itself explained up to 11.3% of the variance in trajectory class membership (OR = 1.50-4.65). When peer cannabis use and the cannabis use PRS were entered into the model simultaneously, both the PRS and peer use continued to be significantly associated with class membership (P < 0.01).

CONCLUSIONS

Genetic propensity to cannabis use derived from heterogeneous samples appears to correlate with longitudinal increases in cannabis use frequency in young adults.

Overlap of heritable influences between cannabis use disorder, frequency of use and opportunity to use cannabis: trivariate twin modelling and implications for genetic design

BACKGROUND

The genetic component of Cannabis Use Disorder may overlap with influences acting more generally on early stages of cannabis use. This paper aims to determine the extent to which genetic influences on the development of cannabis abuse/dependence are correlated with those acting on the opportunity to use cannabis and frequency of use.

METHODS

A cross-sectional study of 3303 Australian twins, measuring age of onset of cannabis use opportunity, lifetime frequency of cannabis use, and lifetime DSM-IV cannabis abuse/dependence. A trivariate Cholesky decomposition estimated additive genetic (A), shared environment (C) and unique environment (E) contributions to the opportunity to use cannabis, the frequency of cannabis use, cannabis abuse/dependence, and the extent of overlap between genetic and environmental factors associated with each phenotype.

RESULTS

Variance components estimates were A = 0.64 [95% confidence interval (CI) 0.58-0.70] and E = 0.36 (95% CI 0.29-0.42) for age of opportunity to use cannabis, A = 0.74 (95% CI 0.66-0.80) and E = 0.26 (95% CI 0.20-0.34) for cannabis use frequency, and A = 0.78 (95% CI 0.65-0.88) and E = 0.22 (95% CI 0.12-0.35) for cannabis abuse/dependence. Opportunity shares 45% of genetic influences with the frequency of use, and only 17% of additive genetic influences are unique to abuse/dependence from those acting on opportunity and frequency.

CONCLUSIONS

There are significant genetic contributions to lifetime cannabis abuse/dependence, but a large proportion of this overlaps with influences acting on opportunity and frequency of use. Individuals without drug use opportunity are uninformative, and studies of drug use disorders must incorporate individual exposure to accurately identify aetiology.

Does thinner right entorhinal cortex underlie genetic liability to cannabis use?

BACKGROUND

Although alterations in medial temporal lobe structures have been previously associated with use of cannabis, one of the most widely used illicit drugs, whether such alterations are a cause or effect of cannabis use has been unclear.

METHODS

In this cross-sectional observational study involving 404 twins/siblings, we have compared cortical thickness and surface area between groups of gender-matched sibling-pairs (concordant cannabis unexposed, concordant exposed and discordant for cannabis exposure) using permutation tests after controlling for potential confounds. Bi-variate polygenic model was used to assess the genetic and environmental contributions underlying cortical morphological phenotypes and frequency of cannabis use.

RESULTS

Cortical thickness of the right entorhinal cortex was significantly lower in the concordant exposed siblings compared to both discordant unexposed and discordant exposed groups [false discovery rate (FDR)-corrected, q < 0.05]. The association between the right entorhinal cortex thickness and frequency of cannabis use is due to the contribution of significant shared additive genetic (ρg = -0.19 ± 0.08; p = 0.02) factors but not unique environment (ρe = 0.05 ± 0.09; p = 0.53). Significantly lower surface area of the right entorhinal cortex in discordant exposed group compared with the discordant unexposed group furnishes preliminary evidence in support of causal effect of cannabis use (FDR-corrected, q < 0.05). However, bi-variate polygenic model-based analysis did not show any significant effect.

CONCLUSIONS

Shared genetic liability may underlie the association between cannabis exposure and thinner right entorhinal cortex. Prospective longitudinal studies are necessary to definitively disentangle the cause-effect relationships of cannabis use.

Genome-wide association meta-analysis of age at first cannabis use

BACKGROUND AND AIMS

Cannabis is one of the most commonly used substances among adolescents and young adults. Earlier age at cannabis initiation is linked to adverse life outcomes, including multi-substance use and dependence. This study estimated the heritability of age at first cannabis use and identified associations with genetic variants.

METHODS

A twin-based heritability analysis using 8055 twins from three cohorts was performed. We then carried out a genome-wide association meta-analysis of age at first cannabis use in a discovery sample of 24 953 individuals from nine European, North American and Australian cohorts, and a replication sample of 3735 individuals.

RESULTS

The twin-based heritability for age at first cannabis use was 38% [95% confidence interval (CI) = 19-60%]. Shared and unique environmental factors explained 39% (95% CI = 20-56%) and 22% (95% CI = 16-29%). The genome-wide association meta-analysis identified five single nucleotide polymorphisms (SNPs) on chromosome 16 within the calcium-transporting ATPase gene (ATP2C2) at P < 5E-08. All five SNPs are in high linkage disequilibrium (LD) (r2  > 0.8), with the strongest association at the intronic variant rs1574587 (P = 4.09E-09). Gene-based tests of association identified the ATP2C2 gene on 16q24.1 (P = 1.33e-06). Although the five SNPs and ATP2C2 did not replicate, ATP2C2 has been associated with cocaine dependence in a previous study. ATP2B2, which is a member of the same calcium signalling pathway, has been associated previously with opioid dependence. SNP-based heritability for age at first cannabis use was non-significant.

CONCLUSION

Age at cannabis initiation appears to be moderately heritable in western countries, and individual differences in onset can be explained by separate but correlated genetic liabilities. The significant association between age of initiation and ATP2C2 is consistent with the role of calcium signalling mechanisms in substance use disorders.

The interplay between externalizing disorders polygenic risk scores and contextual factors on the development of marijuana use disorders

Externalizing disorders have been extensively linked to substance use problems. However, less is known about whether genetic factors underpinning externalizing disorders and environmental features interact to predict substance use disorders (i.e., marijuana abuse and dependence) among urban African Americans. We examined whether polygenic risk scores (PRS) for conduct disorder (CD) and attention-deficit hyperactivity disorder (ADHD) interacted with contextual factors (i.e., parental monitoring, community disadvantage) to influence risk for marijuana use disorders in a sample of African American youth. Participants (N=1,050; 44.2% male) were initially recruited for an elementary school-based universal prevention trial in a Mid-Atlantic city and followed through age 20. Participants reported on their parental monitoring in sixth grade and whether they were diagnosed with marijuana abuse or dependence at age 20. Blood or saliva samples were genotyped using the Affymetrix 6.0 microarrays. The CD and ADHD PRS were created based on genome-wide association studies conducted by Dick et al. (2010) and Demontis et al. (2017), respectively. Community disadvantage was calculated based on census data when participants were in sixth grade. There was an interaction between the CD PRS and community disadvantage such that a higher CD PRS was associated with greater risk for a marijuana use disorder at higher levels of neighborhood disadvantage. This finding should be interpreted with caution owing to the number of significance tests performed. Implications for etiological models and future research directions are presented.

Associations between personality disorders and cannabis use and cannabis use disorder: a population-based twin study

BACKGROUND AND AIMS

Individual differences in DSM-IV personality disorders (PDs) are associated with increased prevalence of substance use disorders. Our aims were to determine which combination of PDs trait scores best predict cannabis use (CU) and cannabis use disorder (CUD), and to estimate the size and significance of genetic and environmental risks in PD traits shared with CU and CUD.

DESIGN

Linear mixed-effects models were used to identify PD traits for inclusion in twin analyses to explore the genetic and environmental associations between the traits and cannabis use.

SETTING

Cross-sectional data were obtained from Norwegian adult twins in a face-to-face interview in 1999-2004 as part of a population-based study of mental health.

PARTICIPANTS

Subjects were 1419 twins (μage  = 28.2 years, range = 19-36) from the Norwegian Institute of Public Health Twin Panel with complete PD and cannabis data.

MEASUREMENTS

PD traits were assessed using DSM-IV criteria. Life-time CU and CUD were based on DSM-IV abuse and dependence criteria, including withdrawal and craving.

FINDINGS

After adjusting for age and sex, antisocial [β = 0.23, 95% confidence interval (CI) = 0.19-0.28] and borderline PDs (β = 0.20, 95% CI = 0.14-0.26) were associated strongly with CU. Antisocial (β = 0.26, 95% CI = 0.21-0.31) and borderline PDs (β = 0.12, 95% CI = 0.06-0.18) were also linked strongly to CUD. Genetic risks in antisocial and borderline PD traits explained 32-60% of the total variance in CU and CUD. Dependent and avoidant PDs explained 11 and 16% of the total variance in CU and CUD, respectively.

CONCLUSIONS

Individual differences in the liability to cannabis use and cannabis use disorder appear to be linked to genetic risks correlated with antisocial and borderline personality disorder traits.

Marijuana use differentially affects cannabinoid receptor expression in early gestational human endometrium and placenta

Marijuana is one of the most abused drugs among pregnant women leading to maternal and fetal abnormalities. Cannabinoids are the active ingredients of marijuana, which interact with cannabinoid receptors such as CNR1 and CNR2 to activate cellular signaling pathways. Human endometrium and placenta are known to express CNR1 and CNR2 and can respond to cannabinoid signaling. In this study, we show that marijuana use significantly increases mRNA or protein expression of CNR1 and CNR2 in human endometrium from the first and early second trimester pregnancies, with minor effects on placental expression of CNRs.