Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Genome-wide Association Study of Cannabis Dependence Severity, Novel Risk Variants, and Shared Genetic Risks.

Author information

Affiliations

  • 1. Section of Biomedical Genetics, Department of Medicine, Boston University School of Medicine, Boston, Massachusetts.
      Authors
    • Sherva R 1
    • Koesterer R 1
    (2 authors)
  • 2. Program in Computational Biology and Bioinformatics, Yale University, New Haven, Connecticut.
      Authors
    • Wang Q 2
    (1 author)
  • 3. Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia 4Mental Illness Research Education and Clinical Center, Veterans Affairs (VA) Stars and Stripes Healthcare Network, Philadelphia VA Medical Center, Philadelphi.
      Authors
    • Kranzler H 3
    (1 author)
  • 4. Program in Computational Biology and Bioinformatics, Yale University, New Haven, Connecticut5Department of Genetics, Yale School of Medicine, West Haven, Connecticut6Department of Biostatistics, Yale School of Public Health, New Haven, Connecticut7VA Coop.
      Authors
    • Zhao H 4
    (1 author)
  • 5. Department of Psychiatry, Yale School of Medicine, New Haven, Connecticut.
      Authors
    • Herman A 5
    (1 author)
    • Show all (7)

ORCIDs linked to this article

JAMA Psychiatry, 01 May 2016, 73(5):472-480
https://doi.org/10.1001/jamapsychiatry.2016.0036 PMID: 27028160 PMCID: PMC4974817

Free full text in Europe PMC

A comment on this article appears in "Genome-wide Significant Associations for Cannabis Dependence Severity: Relevance to Psychiatric Disorders." JAMA Psychiatry. 2016 May 1;73(5):443-4. doi: 10.1001/jamapsychiatry.2016.0046.

Abstract 

Importance

Cannabis dependence (CAD) is a serious problem worldwide and is of growing importance in the United States because cannabis is increasingly available legally. Although genetic factors contribute substantially to CAD risk, at present no well-established specific genetic risk factors for CAD have been elucidated.

Objective

To report findings for DSM-IV CAD criteria from association analyses performed in large cohorts of African American and European American participants from 3 studies of substance use disorder genetics.

Design, setting, and participants

This genome-wide association study for DSM-IV CAD criterion count was performed in 3 independent substance dependence cohorts (the Yale-Penn Study, Study of Addiction: Genetics and Environment [SAGE], and International Consortium on the Genetics of Heroin Dependence [ICGHD]). A referral sample and volunteers recruited in the community and from substance abuse treatment centers included 6000 African American and 8754 European American participants, including some from small families. Participants from the Yale-Penn Study were recruited from 2000 to 2013. Data were collected for the SAGE trial from 1990 to 2007 and for the ICGHD from 2004 to 2009. Data were analyzed from January 2, 2013, to November 9, 2015.

Main outcomes and measures

Criterion count for DSM-IV CAD.

Results

Among the 14 754 participants, 7879 were male, 6875 were female, and the mean (SD) age was 39.2 (10.2) years. Three independent regions with genome-wide significant single-nucleotide polymorphism associations were identified, considering the largest possible sample. These included rs143244591 (β = 0.54, P = 4.32 × 10-10 for the meta-analysis) in novel antisense transcript RP11-206M11.7;rs146091982 (β = 0.54, P = 1.33 × 10-9 for the meta-analysis) in the solute carrier family 35 member G1 gene (SLC35G1); and rs77378271 (β = 0.29, P = 2.13 × 10-8 for the meta-analysis) in the CUB and Sushi multiple domains 1 gene (CSMD1). Also noted was evidence of genome-level pleiotropy between CAD and major depressive disorder and for an association with single-nucleotide polymorphisms in genes associated with schizophrenia risk. Several of the genes identified have functions related to neuronal calcium homeostasis or central nervous system development.

Conclusions and relevance

These results are the first, to our knowledge, to identify specific CAD risk alleles and potential genetic factors contributing to the comorbidity of CAD with major depression and schizophrenia.

Free full text 

Logo of nihpaLink to Publisher's site
JAMA Psychiatry. Author manuscript; available in PMC 2017 May 1.
Published in final edited form as:
PMCID: PMC4974817
NIHMSID: NIHMS804093
PMID: 27028160

Genome-wide Association Study of Cannabis Dependence Severity, Novel Risk Variants, and Shared Genetic Risks

Richard Sherva, PhD, Qian Wang, MS, Henry Kranzler, MD, Hongyu Zhao, PhD, Ryan Koesterer, MS, Aryeh Herman, PsyD, Lindsay A. Farrer, PhD, and Joel Gelernter, MD
Author information Copyright and License information Disclaimer

Abstract

IMPORTANCE

Cannabis dependence (CAD) is a serious problem worldwide and is of growing importance in the United States because cannabis is increasingly available legally. Although genetic factors contribute substantially to CAD risk, at present no well-established specific genetic risk factors for CAD have been elucidated.

OBJECTIVE

To report findings for DSM-IV CAD criteria from association analyses performed in large cohorts of African American and European American participants from 3 studies of substance use disorder genetics.

DESIGN, SETTING, AND PARTICIPANTS

This genome-wide association study for DSM-IV CAD criterion count was performed in 3 independent substance dependence cohorts (the Yale-Penn Study, Study of Addiction: Genetics and Environment [SAGE], and International Consortium on the Genetics of Heroin Dependence [ICGHD]). A referral sample and volunteers recruited in the community and from substance abuse treatment centers included 6000 African American and 8754 European American participants, including some from small families. Participants from the Yale-Penn Study were recruited from 2000 to 2013. Data were collected for the SAGE trial from 1990 to 2007 and for the ICGHD from 2004 to 2009. Data were analyzed from January 2, 2013, to November 9, 2015.

MAIN OUTCOMES AND MEASURES

Criterion count for DSM-IV CAD.

RESULTS

Among the 14 754 participants, 7879 were male, 6875 were female, and the mean (SD) age was 39.2 (10.2) years. Three independent regions with genome-wide significant single-nucleotide polymorphism associations were identified, considering the largest possible sample. These included rs143244591 (β = 0.54, P = 4.32 × 10−10 for the meta-analysis) in novel antisense transcript RP11-206M11.7; rs146091982 (β = 0.54, P = 1.33 × 10−9 for the meta-analysis) in the solute carrier family 35 member G1 gene (SLC35G1); and rs77378271 (β = 0.29, P = 2.13 × 10−8 for the meta-analysis) in the CUB and Sushi multiple domains 1 gene (CSMD1). Also noted was evidence of genome-level pleiotropy between CAD and major depressive disorder and for an association with single-nucleotide polymorphisms in genes associated with schizophrenia risk. Several of the genes identified have functions related to neuronal calcium homeostasis or central nervous system development.

CONCLUSIONS AND RELEVANCE

These results are the first, to our knowledge, to identify specific CAD risk alleles and potential genetic factors contributing to the comorbidity of CAD with major depression and schizophrenia.

After nicotine, cannabis is the most widely abused drug worldwide.1 In the United States, the accelerated decriminalization of cannabis is based on the erroneous perception that it is relatively harmless.2 In fact, cannabis use produces craving,3 dependence,4 and drug-seeking behavior,5 as with the use of other substances. Despite these risks, the prevalence of cannabis use and cannabis use disorders has dramatically increased since 2001,6 and the political momentum to increase availability has continued. Use of cannabis early in life is associated with an increased risk for schizophrenia (SCZ),7 and sets of SCZ-associated risk alleles predict cannabis use.8 Cannabis use is also a risk factor for depressive symptoms,9 and a twin study showed cannabis dependence (CAD) to be associated with an elevated risk for major depressive disorder (MDD).10 Substance use and other psychiatric illnesses may share common genetic risk factors; or reverse causation, self-medication, or confounding by other factors may explain their co-occurrence.

Despite knowledge of the neurobiology of the endocannabinoid system and its response to tetrahydrocannabinol, little is known about specific genetic factors influencing susceptibility to CAD or cannabis abuse. A twin study showed that several aspects of cannabis use are heritable, including an early opportunity to use (h2 = 72%), early onset of use (h2 = 80%), lifetime use of cannabis 11 or more times (h2 = 76%), and cannabis abuse or dependence (h2 = 21%–72%), where h2 is hertiability.11–13 Possible evidence of linkage of CAD on chromosome 1614 and linkage and association encompassing the neuregulin 1 gene (NRG1 [OMIM 142445]; known as a possible SCZ risk gene15) on chromosome 816 have been found. Despite several genome-wide association studies (GWAS) on cannabis-related traits, no genome-wide significant (GWS) associations were observed for initiation of use17 or for CAD.18 Herein we report on findings for DSM-IV CAD criteria from association analyses performed in large cohorts of African American and European American participants from 3 studies of substance use disorder genetics who underwent genotyping with genome-wide microarrays. The primary cohort has been used in previous studies to identify genes associated with opioid (OD),19 cocaine (CD),20 alcohol (AD),21 and nicotine (ND) dependence22 and posttraumatic stress disorder.23

Methods

Participants and Diagnostic Procedures

The samples included 6000 African American and 8754 European American participants (race was assigned based on genetic data; eMethods in the Supplement) from the following 3 studies: (1) the Yale-Penn Study cohort of small nuclear families and unrelated individuals (2020 individuals in 850 families and 6951 unrelated individuals), collected to study the genetics of substance dependence19–21; (2) the GWAS data set from the Study of Addiction: Genetics and Environment (SAGE),24–27 collected to study the genetics of AD, ND, and CD (183 individuals in 89 families and 3707 unrelated individuals); and (3) the GWAS International Consortium on the Genetics of Heroin Dependence (ICGHD),28,29 a collaboration formed to identify genes associated with heroin dependence risk (66 individuals in 33 families and 1827 unrelated individuals). The SAGE and ICGHD data sets are publicly available via application. The present study received institutional review board approval from all participating institutions, and written informed consent was obtained from all study participants.

Key Points

Question

What specific genetic variants contribute to cannabis dependence risk?

Findings

Three regions had genome-wide significant evidence of association with cannabis dependence and evidence of genetic overlap between cannabis dependence and schizophrenia and major depressive disorder.

Meaning

Cannabis dependence has a genetic risk component that may overlap with other psychiatric disorders.

Participants from the Yale-Penn Study were recruited from 2000 to 2013. These participants were administered the Semi-Structured Assessment for Drug Dependence and Alcoholism30 to derive DSM-IV diagnoses of lifetime CAD and other major psychiatric traits. Data were collected for the SAGE trial from 1990 to 2007, and participants underwent phenotyping with the Semi-Structured Assessment for the Genetics of Alcoholism.31 Data were collected for the ICGHD from 2004 to 2009, and participants completed a comprehensive psychiatric diagnostic interview based on the Semi-Structured Assessment of the Genetics of Alcoholism–Australia.31 The method of phenotyping was similar across the 3 samples. Additional information about recruitment, genotyping, imputation, and quality control for the study cohorts is provided in eMethods in the Supplement.

Statistical Analysis

Data were analyzed from January 2, 2013, to November 9, 2015. Association analyses were performed using a count of DSM-IV CAD criteria (0–7) as the outcome variable and the imputed minor allele dosage (adjusted for sex, age, and the first 3 ancestry principal components) as a predictor variable. This ordinal trait model has greater power to detect genetic associations than a univariate model based on disease status because of greater information content and improved specificity of the dependence measure. Association tests were performed using linear association models embedded in generalized estimating equations to correct for correlations among related individuals.32 Analyses were performed separately within each data set and population group, and the results were combined by meta-analysis using the inverse variance method implemented in the program METAL.33 Genomic inflation factors (λ) were calculated within each subpopulation, and P values were corrected accordingly. We performed a second correction for the λ factor calculated after the meta-analysis.

For the primary analysis, individuals were included regardless of cannabis exposure. As secondary analyses, individuals who reported never having used cannabis were excluded, and the primary model was repeated adjusting for the criterion counts for AD, CD, and OD. Participants from 2 genotyping batches in the Yale-Penn cohort (Yale-Penn 1 and Yale-Penn 2) were combined with the SAGE sample to form a discovery data set. A sample consisting of the ICGHD data and additional samples from the Yale-Penn cohort who did not undergo genotyping at the time of the discovery analyses (Yale-Penn 3) were used to replicate the top associations.

Cross-Disorder Analysis

We attempted to uncover shared genetic variation between CAD and 5 psychiatric disorders, including SCZ, MDD, bipolar affective disorder, attention-deficit/hyperactivity disorder, and autism spectrum disorder using the GWAS analysis reported herein and publicly available GWAS results from the Psychiatric Genomics Consortium (http://www.med.unc.edu/pgc/).34 To explore cross-disorder genetic relationships, we used stratified quintile-quintile (QQ) plots to evaluate the relative enrichment of single-nucleotide polymorphisms (SNPs) associated with both disorders. The QQ plots, which contrast the observed distribution of P values with the expected distribution under the null hypothesis (uniform in GWAS), were used to assess P value inflation in the GWAS results. Grouping associated SNPs for one disorder and comparing (across groups) the QQ plots of another disorder, however, could also reveal the enrichment of GWAS signals between disorders, which made them suitable for cross-disorder enrichment screening.

We also applied a statistical framework for pleiotropy analysis, Genetic Analysis Incorporating Pleiotropy and Annotation (GPA).35 The GPA was built as a mixture model with parameters estimated using an efficient expectation-maximization algorithm, where associated SNPs were modeled with a β [α, 1] distribution and unassociated SNPs with a uniform [0, 1] distribution. A likelihood ratio test assessed the significance of pleiotropy between disorders. The GPA also detected the SNPs that were pleiotropic by calculating the posterior probability of association with both disorders.

Results

Participant demographic characteristics and the correlation between the criterion counts for CAD and other substance use disorder traits are shown in Table 1. The DSM-IV CAD criterion counts were significantly (P < .05) correlated with the criteria counts for AD, CD, OD, and ND. The correlations varied by sample and population and ranged from r2 = 0.15 for OD to r2 = 0.61 for CD criteria. The CAD criterion counts were significantly heritable in European American (19%–25%; P = .006) but not African American (10%–11%; P = .08) participants. eFigure 1 in the Supplement shows a histogram of the CAD criterion count in African American and European American participants in each cohort; 3 or more criteria indicate a diagnosis of CAD. The criterion count distribution is very similar in African American and European American participants. In the Yale-Penn sample, where comorbid psychiatric diagnoses were available, CAD was significantly associated with MDD in African American participants (odds ratio, 1.07; P = .006) but not SCZ, bipolar affective disorder, attention-deficit/hyperactivity disorder, or autism spectrum disorder. Cannabis dependence was not associated with any of these disorders in European American participants.

Table 1

Demographic Characteristics of the Study Cohorts

SampleaAge, Mean (SD), yFemale Sex, No./Total No.DSM-IV CAD Criterion Count, Mean (SD)CAD Diagnosis, No./Total No. of ParticipantsCorrelation With DSM-IV CAD Criterion Count, r2
ADCDNDOD
Yale-Penn African American41.3 (9.7)2209/47501.7 (2.2)1296/47500.400.300.360.15
Yale-Penn European American38.6 (11.5)1712/42212.0 (2.3)1388/42210.270.320.310.23
SAGE African American39.9 (7.3)638/12501.4 (2.2)276/12500.430.430.390.17
SAGE European American38.4 (9.7)1478/26401.0 (1.9)434/26400.510.610.400.41
ICGHD36.2 (9.1)838/18933.2 (2.5)1062/18930.240.240.250.33

Abbreviations: AD, alcohol dependence; CAD, cannabis dependence; CD, cocaine dependence; ICGHD, International Consortium on the Genetics of Heroin Dependence; ND, nicotine dependence; OD, opioid dependence; SAGE, Study of Addiction: Genetics and Environment.

aSamples are described in the Participants and Diagnostic Procedure subsection of the Methods section.

GWAS Results

Manhattan and QQ plots for the meta-analysis discovery GWAS results for African American and European American Yale-Penn 1 and 2 and SAGE cohorts are displayed in eFigures 2 and 3 in the Supplement. We found little evidence of P value inflation. Table 2 shows associations in the discovery sample with P < 1.0 × 10−5 in African American or European American participants or the combined meta-analysis, trimmed for linkage disequilibrium. eTable 1 in the Supplement shows the same results, together with additional information about each SNP, including the results within each discovery sample subgroup, after excluding individuals with no cannabis exposure, and after adjusting for comorbid substance use disorders. We identified GWS associations with reliably imputed SNPs in 3 distinct regions (Table 2), 2 specific to African American participants and 1 in the combined sample. First, rs186825689 (P = 1.86 × 10−8 for the African American meta-analysis) is located 12.4 kb upstream from the gene encoding S100 calcium binding protein (S100B) with contributions from both informative African American samples. Second, rs143244591 (P = 2.18 × 10−8 for the African American meta-analysis) maps to a novel antisense transcript RP11-206M11.7 (Havana gene: OTTHUMG00000159583) located in the gene of the same name on chromosome 3 with at least nominally significant evidence in each of the 3 African American samples. Third, rs77378271(P = 2.76 × 10−8 for the European American meta-analysis) is an intronic SNP in the CUB and Sushi multiple domains 1 gene (CSMD1 [OMIM 608397]) with evidence of association in 3 of the 6 samples. We also identified consistent, non-GWS evidence of association in the combined sample of European American and African American participants with a large block of SNPs in and around the phosphatidylinositol 4-kinase type 2β gene (PI4K2B [OMIM 612101]), with consistent effect direction in every European American and African American population tested (minimum P = 1.74 × 10−7 for the meta-analysis). This signal was GWS when individuals without cannabis exposure were excluded (minimum P = 2.98 × 10−8 for the meta-analysis).

Table 2

SNPs Associated With DSM-IV CAD at P < 1.0 × 10−5 in the Discovery Meta-analysis Trimmed for SNPs in Linkage Disequilibrium

ChromosomeBase Pair
Positiona		Effect
Allele		Reference
Allele		SNPGeneEffect Allele
Frequency
P Value for Meta-analysis
African
American
Cohort		European
American
Cohort		African
American
Cohort		European
American
Cohort		All Participants
188729683CTrs74823926NA0.960.975.26 × 10−76.36 × 10−11.40 × 10−5
239166173TGrs114383460ARHGEF330.96NA1.09 × 10−6NANA
278028838TArs12621150NA0.850.907.67 × 10−11.05 × 10−62.91 × 10−4
2100451676TCrs7586604AFF30.390.164.15 × 10−56.19 × 10−31.06 × 10−6
2103764414GArs144605126NA0.96NA8.67 × 10−7NANA
2118490901GArs150064803NA0.95NA3.41 × 10−7NANA
2167214714GArs143020225SCN9A0.95NA7.19 × 10−7NANA
3149013935GArs143244591bRP11-206M11.70.96NA2.18 × 10−8bNANA
425201318ATrs73252553PI4K2B0.960.901.18 × 10−32.28 × 10−51.66 × 10−7
4119716950ACrs28595532SEC24D0.920.942.02 × 10−71.08E-011.13 × 10−6
511892384TCrs114311699CTNND20.96NA3.78 × 10−7NANA
5177746600GCrs10066744COL23A10.96NA4.82 × 10−7NANA
651221457AGrs17665889NA0.940.901.51 × 10−19.41 × 10−82.58 × 10−4
784952631AGrs12534830NA0.880.667.76 × 10−31.54 × 10−54.52 × 10−7
83073489AGrs77378271bCSMD10.960.942.13 × 10−12.76 × 10−8b4.60 × 10−8b
929364327GTrs10969106NANA0.97NA7.39 × 10−8NA
1031981385TCrs115553536NA0.94NA6.46 × 10−7NANA
1043592809GCrs74400468RET0.960.975.53 × 10−21.59 × 10−66.46 × 10−7
1070490106ATrs12218439CCAR10.960.961.01 × 10−62.01 × 10−11.13 × 10−4
1095659958AGrs146091982SLC35G10.95NA1.95 × 10−7NANA
1120561010CGrs73443003NA0.75NA1.31 × 10−6NANA
1181433204AAAAGrs200453611NA0.910.849.43 × 10−76.24 × 10−16.81 × 10−4
11108899423GTAGrs200391037NA0.960.881.02 × 10−13.72 × 10−61.32 × 10−6
1256274155TCrs193047854NA0.97NA7.06 × 10−7NANA
2021706604AATrs199783889NA0.94NA3.32 × 10−7NANA
2118019319TCrs78068107NA0.96NA1.02 × 10−6NANA
2148006053ACrs186825689bNA0.96NA1.86 × 10−8bNANA
aIndicates in human genome assembly build 37.
bIndicates genome-wide significant SNPs and P values.

Replication Results

The SNPs in Table 2 were tested for CAD association in the 2 replication samples (ICGHD and Yale-Penn 3). Table 3 shows the replication cohort-specific results for these SNPs, with the meta-analysis results from the discovery phase and the discovery + replication phase. The smallest P value in the ICGHD cohort among the 13 SNPs that could be reliably imputed and analyzed (this cohort was European Australian) was at rs74823926 (P = .064) in an intergenic region on chromosome 1. Several associations, however, were replicated in the Yale-Penn 3 sample (Table 3). The P values for 2 of the 3 GWS SNPs improved after meta-analysis with the replication cohorts (rs143244591 in RP11-206M11.7, from 1.38 × 10−8 to 4.32 × 10−10; rs77378271 in CSMD1, from 2.84 × 10−8 to 2.13 × 10−8), as did the P value for another SNP, rs146091982 in the solute carrier family 35 member G1 (SLC35G1 [Ensembl ENSG00000176273]) (from 1.31 × 10−7 to 1.33 × 10−9). The signal in PI4K2B also improved (P = 5.57 × 10−8 for the full meta-analysis). However, rs186825689 near S100B was no longer GWS (P = 8.27 × 10−8) in the full meta-analysis. The Figure shows Manhattan plots for the regions encompassing RP11-206M11.7 (Figure, A), SLC35G1 (Figure, B), CSMD1 (Figure, C), and PI4K2B (Figure, D) in the discovery sample and after meta-analysis with the replication samples.

An external file that holds a picture, illustration, etc. Object name is nihms804093f1.jpg
Regional Manhattan Plots of Association Results for DSM-IV Cannabis Dependence Criterion Count in 4 Genomic Regions

Association results from single-nucleotide polymorphisms (SNPs) in 4 regions. A, The 148.8- to 149.2-MB region encompassing RP11-206M11.7 on chromosome 3 in the Yale-Penn and Study of Addiction: Genetics and Environment (SAGE) African American participants. B, The 95.3- to 96-MB region encompassing SLC35G1 on chromosome 10 in the Yale-Penn and SAGE African American participants. C, The 2.8- to 4.8-MB region on chromosome 8 encompassing CSMD1 in the Yale-Penn, SAGE, and International Consortium on the Genetics of Heroin Dependence (ICGHD) African American and European American participants. D, The 25.07- to 25.43-MB region encompassing PI4K2B on chromosome 4 in the Yale-Penn, SAGE, and ICGHD African American and European American participants. In A and B, the SNPs are color coded according to the correlation coefficient (r2) in the 1000 Genomes African samples with the most significant SNP. In C and D, results from the African American and European American participants were combined, and no linkage disequilibrium information was displayed. The light purple circle represents the −log10 P value for the most significant regional SNP in the meta-analysis of the discovery samples; the purple diamond, the result for that SNP after meta-analysis with the replication sample(s). The light blue line and right y-axis show the observed recombination rate.

Table 3

Association Results in the Discovery and Replication Samples for the SNPs Shown in Table 2

SNPP Value for Replication Cohort
P Value for Meta-analysis
Directiond
Yale-Penn 3 Cohorta
ICGHD CohortDiscovery Cohortb
Discovery + Replication Cohortc
All Participants
African AmericanEuropean AmericanAfrican AmericanEuropean AmericanAfrican AmericanEuropean American
rs748239269.40 × 10−37.54 × 10−16.36 × 10−25.26 × 10−76.36 × 10−13.34 × 10−61.59 × 10−18.62 × 10−6+x++++−++
rs1143834607.22 × 10−1NANA1.09 × 10−6NA1.93 × 10−6NA1.93 × 10−6+x + xxx−xx
rs126211509.68 × 10−11.15 × 10−13.61 × 10−17.67 × 10−11.05 × 10−67.75 × 10−17.92 × 10−65.13 × 10−4++++−+−−+
rs75866042.80 × 10−12.78 × 10−16.04 × 10−14.15 × 10−56.19 × 10−31.76 × 10−51.69 × 10−21.81 × 10−6−−−−−−−+−
rs1446051261.18 × 10−1NANA8.67 × 10−7NA3.84 × 10−6NA3.84 × 10−6+x + x + x−xx
rs1500648031.74 × 10−1NANA3.41 × 10−7NA1.41 × 10−6NA1.41 × 10−6+x + x + x−xx
rs1430202251.53 × 10−1NANA7.19 × 10−7NA1.95 × 10−7NA1.95 × 10−7+x + xxx + xx
rs1432445913.24 × 10−3NANA2.18 × 10−8eNA4.32 × 10−10eNA4.32 × 10−10e+x + x + x+xx
rs732525536.40 × 10−18.32 × 10−24.10 × 10−11.18 × 10−32.28 × 10−52.25 × 10−35.45 × 10−65.57 × 10−8++++++−++
rs285955328.06 × 10−16.07 × 10−22.24 × 10−12.02 × 10−71.08 × 10−11.89 × 10−71.52 × 10−13.60 × 10−6++++++++−
rs1143116997.30 × 10−1NANA3.78 × 10−7NA2.75 × 10−7NA2.75 × 10−7+x + x + x+xx
rs100667443.89 × 10−1NANA4.82 × 10−7NA2.27 × 10−7NA2.27 × 10−7+x + xxx + xx
rs176658892.17 × 10−11.17 × 10−17.65 × 10−11.51 × 10−19.41 × 10−82.69 × 10−15.39 × 10−68.41 × 10−4−+++−++−+
rs125348305.85 × 10−18.67 × 10−13.99 × 10−17.76 × 10−31.54 × 10−51.05 × 10−23.40 × 10−51.12 × 10−6++++++−−+
rs773782719.25 × 10−24.19 × 10−2f7.95 × 10−12.13 × 10−12.76 × 10−8e1.07 × 10−15.16 × 10−82.13 × 10−8e+−++x+++−
rs10969106NA5.50 × 10−18.47 × 10−1NA7.39 × 10−8NA1.74 × 10−71.74 × 10−7x + x + x+xx+
rs1155535361.80 × 10−1NANA6.46 × 10−7NA2.14 × 10−6NA2.14 × 10−6+x + x + x−xx
rs744004683.49 × 10−26.61 × 10−15.66 × 10−15.53 × 10−21.59 × 10−61.22 × 10−22.46 × 10−59.82 × 10−7++−+++++−
rs122184392.29 × 10−11.01 × 10−13.21 × 10−11.01 × 10−62.01 × 10−14.91 × 10−62.13 × 10−14.20 × 10−4x++−++−−+
rs1460919828.84 × 10−4NANA1.95 × 10−7NA1.33 × 10−9eNA1.33 × 10−9e+x + x + x+xx
rs734430032.53 × 10−2NANA1.31 × 10−6NA1.20 × 10−7NA1.20 × 10−7−x−x−x−xx
rs2004536117.29 × 10−15.20 × 10−19.40 × 10−19.43 × 10−76.24 × 10−11.56 × 10−65.34 × 10−11.12 × 10−3++++++−++
rs2003910379.28 × 10−13.09 × 10−2f4.19 × 10−11.02 × 10−13.72 × 10−61.04 × 10−14.32 × 10−51.10 × 10−5++++++−−+
rs1930478546.74 × 10−1NANA7.06 × 10−7NA5.51 × 10−7NA5.51 × 10−7+x + xxx + xx
rs1997838894.38 × 10−1NANA3.32 × 10−7NA1.12 × 10−6NA1.12 × 10−6+x + xxx−xx
rs780681072.90 × 10−2NANA1.02 × 10−6NA1.31 × 10−5NA1.31 × 10−5+x + xxx−xx
rs1868256894.51 × 10−1NANA1.86 × 10−8eNA8.27 × 10−8NA8.27 × 10−8+x + xxx−xx

Abbreviations: ICGHD, International Consortium on the Genetics of Heroin Dependence; NA, not applicable; SAGE, Study of Addiction: Genetics and Environment; SNP, single-nucleotide polymorphism.

aIndicates participants in the Yale-Penn cohort who did not undergo genotyping at the time of the discovery analyses.
bIndicates participants in the Yale-Penn cohort (Yale-Penn 1 and 2) who underwent genotyping for the discovery analysis and in the SAGE cohort.
cIndicates all cohorts.
dFor the effect direction, the order of the symbols is the Yale-Penn 1 African American cohort, Yale-Penn 1 European American cohort, SAGE African American cohort, SAGE European American cohort, Yale-Penn 2 African American cohort, Yale-Penn 2 European American cohort, Yale-Penn 3 African American cohort, Yale-Penn 3 European American cohort, and the ICGHD European American cohort. + Indicates effect allele (listed in Table 2) is associated with an increase in cannabis dependence (CAD) criterion count; −, effect allele is associated with a decrease in CAD criterion count; and x, a valid effect estimate could not be obtained. See the Statistical Analysis subsection of the Methods section for an explanation of Yale-Penn 1, 2, and 3.
eIndicates genome-wide significant P values.
fIndicates P < .05 in the replication sample.

Cross-Disorder Analysis Results

The QQ plots of 5 Psychiatric Genomics Consortium traits (SCZ, bipolar affective disorder, autism spectrum disorder, attention-deficit/hyperactivity disorder, and MDD) were stratified based on our CAD GWAS results at significance levels of P < .05, P < .01, P < 1 × 10−3, and P < 1 × 10−4. We observed enrichment of the MDD GWAS signal in the CAD GWAS (eFigure 4 in the Supplement) in European American participants, but no clear enrichment for the other 4 psychiatric disorders in either population group (eFigure 5 in the Supplement).

We used GPA to test the significance of pleiotropy between CAD and the same 5 psychiatric disorders (eMethods in the Supplement). For each disease pair, we estimated the percentage of SNPs shared by 2 diseases and tested the significance of pleiotropy (eTable 2 in the Supplement). The European American population yielded significant evidence of CAD-MDD pleiotropy (P = 2.39 × 10−5); genome wide, 1.7% of all imputed SNPs were estimated to be associated with both CAD and MDD. Of these, rs10954732 in P450 oxidoreductase(POR[OMIM124015]) had the largest posterior probability (although not significant) of association with both traits (P = 2.59 × 10−6 for CAD; P = .02 for MDD; posterior probability, 0.70).

Discussion

We report herein the first GWS results for CAD to our knowledge. The sample includes a large proportion (18%–36%, depending on race and cohort) of individuals with CAD from 2 ancestral populations in 3 independent cohorts. We identified 3 regions with GWS SNPs imputed to the 1000 Genomes reference panel that implicate several biological processes and provide insight into the biology of CAD, including evidence of an inflammatory component in the disorder, which may also mediate risk for SCZ36 and MDD.37,38 The smallest P value observed (P = 4.32 × 10−10) was at rs143244591 in RP11-206M11.7. Little is known about this antisense transcript or which, if any, genes it regulates. Minor alleles were protective. The next most significant locus was SLC35G1 (rs146091982, P = 1.33 × 10−9), a potential member of the drug/metabolite transporter superfamily (EamA, previously DUF6). Ubiquitously expressed, SLC35G1 binds stromal interaction molecule 1, a calcium sensor that communicates the calcium load within the endoplasmic reticulum to store-operated channels in the plasma membrane39 when calcium stores in the endoplasmic reticulum are depleted.40 The SLC35G1–stromal interaction molecule 1 complex likely regulates the activity of the transporters that coordinate cytosolic calcium through modulation of pump activities.40 The third GWS locus, CSMD1 (rs77378271; P = 2.13 × 10−8), is highly expressed in the growth cones of developing central nervous system neurons, where it likely acts as a regulator of complement activation and inflammation.41 Different SNPs in CSMD1 have been associated with SCZ at the GWS level.42 Thus, CSMD1 is the second gene to be implicated in both disorders (after NRG116) and may explain at least part of their shared genetic susceptibility.

Two other established SCZ risk genes, RIMS1 (OMIM 606629) (minimum SNP, P = 1.59 × 10−5) and MEF2C (minimum SNP, P = 5.22 × 10−5), showed suggestive association with CAD. MEF2C is highly expressed in developing mammalian neurons and is thought to mediate calcium-dependent survival of neurons that have made the appropriate synaptic connections.43 From a biological perspective, RIMS1 is immediately relevant; RIMS1 acts as a scaffold protein that regulates synaptic vesicle exocytosis, affecting cannabinoid receptor 1 (CR1)–mediated long-term suppression of γ-aminobutyric acid release, ultimately mediating presynaptic forms of long-term plasticity.44 Minor alleles at rs142305709 in RIMS1 were associated with fewer CAD criteria in African American participants. We observed at least a nominally significant signal in both Yale-Penn African American analysis subsets and a non-significant trend in SAGE African American participants.

Limitations of the GWAS findings should be noted. One of the significant SNPs identified (rs143244591 on chromosome 3) has little supportive evidence for association from other SNPs in the region, possibly owing to low linkage disequilibrium. However, despite stringent imputation quality thresholds for including SNPs in the analysis (r2≥0.8) and evidence of an association in the replication sample, this signal may represent an imputation artifact. Second, although none of the GWS SNPs identified in the full GWAS analysis are rare, they could be described as infrequent, with minor allele frequencies in a range sometimes associated with false-positive results (4%–6%). Also, of the GWS regions, only CSMD1 showed evidence of associations in European American and African American participants. The region containing PI4K2B, which became GWS after excluding unexposed individuals (see below), was also at least nominally associated with CAD in both populations. The 2 African American–specific SNPs were rare or monomorphic in European American participants. The lack of association in European American participants could be owing to different linkage disequilibrium patterns or the absence of causal variants. The Yale-Penn samples who underwent genotyping on the HumanOmni1-Quad and Human Core Exome chips showed more consistent results than the corresponding SAGE population, which is not surprising insofar as SAGE participants were recruited from different areas and ascertained using different criteria (AD, CD, and OD in Yale-Penn and primarily AD and ND in SAGE). The difference in ascertainment criteria (use of licit vs illicit drugs) across studies likely explains the fact that the proportion of cannabis-exposed individuals varied significantly across cohorts (2293 in SAGE population [76.9%] and 7626 in the Yale-Penn population [85.0%]). The limitations of phenotypic distribution and population differences are more relevant to the Australian ICGHD replication cohort and may explain the lack of replication in this cohort. Despite this, we obtained statistically significant evidence for formal replication for the SNP in SLC35G1 and stronger evidence for association at many of the top SNPs after including the replication samples. Finally, these cohorts have higher rates of polysubstance dependence than the general population and may not be generalizable to individuals who only use cannabis.

Effect of Exposure Status and Comorbidity

Because the inclusion of genetically at-risk individuals who never initiated cannabis use might have influenced our results, we repeated the primary analyses in the discovery cohort after removing unexposed individuals. Two of the 3 regions identified remained GWS (eTable 1 in the Supplement). The P value for rs143244591 on chromosome 3 improved slightly (P = 1.13 × 10−8, meta-analysis exposed) and was associated at P ≤ .02 in each of the African American subgroups. The signal at rs77378271 in CSMD1 was almost identical (P = 2.95 × 10−8, meta-analysis exposed) and showed association at P < 5.07 × 10−3 in 2 of the 3 European American subgroup and at P = 4.46 × 10−4 in 1 of the African American subgroups. In addition, the block of SNPs in and around PI4K2B became GWS with a consistent effect direction (minor alleles being protective) in every European American and African American population tested and became GWS (minimum P = 2.98 × 10−8, meta-analysis exposed, at rs147170184). The evidence for pleiotropy between CAD and MDD was attenuated substantially (P = .60) after excluding unexposed participants. That the removal of unexposed individuals from the analysis had a relatively minor effect on the primary findings and actually improved the strength of some suggests that any loss in power owing to the smaller sample was offset by an increase in phenotypic precision. In the pleiotropy analysis, which relies on genome-level association results and is not limited to the most significantly associated SNPs, the power loss apparently outweighed any increase in precision. The significance of each of the top SNPs was modestly attenuated after adjusting for the DSM-IV criterion counts for AD, CD, and OD (eTable 1 in the Supplement).

Ion Homeostasis and Addiction

The previously published GWAS of OD19 and CD20 in a subset of this sample each identified risk genes and pathways involved in the regulation of neuronal calcium and potassium, and the pathway involving synaptic long-term potentiation was also identified for OD. Also, a cross-disorder analysis identified calcium signaling in neurons as a pathway mediating 5 psychiatric diseases, including SCZ and MDD.34 The GWS association in SLC35G1 and GWS (in the discovery sample only) associations in and around S100B suggest ion homeostasis may play a role in CAD risk.

Shared Risk for CAD and Other Psychiatric Disorders

Many previous studies7,8,45,46 have focused on the relationship between CAD and SCZ, whereas the correlation between CAD and MDD has received much less attention. Although depressive disorders are highly comorbid with CAD in clinical settings,47 to our knowledge no previous genomics study has explored CAD-MDD pleiotropy. We found some evidence for genetic correlation between the risks for CAD and MDD. The existence of shared genetic factors for CAD-MDD is supported by the overlap in SNPs nominally associated with both traits, although we found no significant evidence of pleiotropy at any single SNP. We also found limited support for the possibility that such a relationship exists for CAD and SCZ based on relatively strong signals for both traits with variants in CSMD1 (although not the same variants). Nongenetic explanations such as patients with SCZ or MDD mediating the symptoms of these disorders with cannabis use might also explain the comorbidity. These analyses are exploratory, and follow-up studies to validate and extend these findings are necessary.

Conclusions

This study provided the first GWS evidence to our knowledge for SNPs associated with CAD via GWAS in 3 distinct genomic locations. These findings will lead our understanding of genetic vulnerability to CAD in new directions that can inform our understanding of the biology of CAD. We obtained entirely novel evidence of genetic overlap between CAD and MDD and conclude that CSMD1 may be a candidate gene that affects the risk for CAD and SCZ, a topic of considerable research interest.7,48–50 These results also suggest that common pathways (nervous system development, inflammation, and ion homeostasis) mediate the risk for multiple psychiatric disorders and dependence on multiple substances, including cannabis.

Supplementary Material

Supplemental

Acknowledgments

Funding/Support: This study was supported by grants RC2 DA028909, R01 DA12690, R01 DA12849, R01 DA18432, R01 AA11330, and R01 AA017535 from the National Institutes of Health (NIH), the Veterans Affairs Connecticut Healthcare Center, and the Philadelphia Veterans Affairs Mental Illness Research, Education and Clinical Center. Funding support for the Study of Addiction: Genetics and Environment (SAGE) was provided by grant U01 HG004422 from the NIH Genes, Environment and Health Initiative. SAGE is one of the genome-wide association studies funded as part of the Gene Environment Association Studies under the Genes, Environment and Health Initiative. The portion of the genotyping that was performed at the Johns Hopkins University Center for Inherited Disease Research was supported by grant U01HG004438 from the NIH Genes, Environment and Health Initiative, and by contract HHSN268200782096C (High Throughput Genotyping for Studying the Genetic Contributions to Human Disease) from the National Institute on Alcohol Abuse and Alcoholism, the National Institute on Drug Abuse, NIH. Support for collection of datasets and samples was provided by grant U10 AA008401 (Collaborative Study on the Genetics of Alcoholism), grant P01 CA089392 (Collaborative Genetic Study of Nicotine Dependence), and grant R01 DA013423 (Family Study of Cocaine Dependence) from the NIH. The International Consortium on the Genetics of Heroin Dependence (principal investigator, Elliot Nelson, MD) was funded by HHSN268200782096C (NIH contract: High Throughput Genotyping for Studying the Genetic Contributions to Human Disease) and HHSN268201100011I (NIH contract: High Throughput Genotyping for Studying the Genetic Contributions to Human Disease).

Footnotes

Additional Information: The publicly available datasets used for analysis were obtained from dbGaP (http://www.ncbi.nlm.nih.gov/gap) through accession numbers phs000092.v1.p and phs000277.v1.p1.

Author Contributions: Dr Sherva had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Study concept and design: Kranzler, Koesterer, Farrer, Gelernter.

Acquisition, analysis, or interpretation of data: All authors.

Drafting of the manuscript: Sherva, Wang.

Critical revision of the manuscript for important intellectual content: Sherva, Kranzler, Zhao, Koesterer, Herman, Farrer, Gelernter.

Statistical analysis: Sherva, Wang, Zhao, Koesterer, Farrer.

Obtained funding: Kranzler, Gelernter.

Administrative, technical, or material support: Kranzler, Farrer, Gelernter.

Study supervision: Kranzler, Zhao, Farrer, Gelernter.

Conflict of Interest Disclosures: Dr Kranzler reports being a consultant or an advisory board member for Alkermes, Indivior, Lundbeck, and Otsuka (unrelated to the present study) and being a member of the American Society of Clinical Psychopharmacology’s Alcohol Clinical Trials Initiative, which is supported by AbbVie, Ethypharm, Lilly, Lundbeck, and Pfizer. No other disclosures were reported.

Role of the Funder/Sponsor: The funding sources had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Additional Contributions: Work in recruitment and assessment was performed at Yale University School of Medicine and the APT Foundation by James Poling, PhD; at McLean Hospital by Roger Weiss, MD; at the Medical University of South Carolina by Kathleen Brady, MD, PhD, and Raymond Anton, MD; and at the University of Pennsylvania by David Oslin, MD. Genotyping services for a part of our genome-wide association study were provided by the Center for Inherited Disease Research and Yale University (Center for Genome Analysis), which is fully funded by Federal contract N01-HG-65403 from the NIH to The Johns Hopkins University. Ann Marie Lacobelle, BS, Catherine Aldi, BS, and Christa Robinson, BS, provided technical assistance. The Semi-structured Assessment for Drug Dependence and Alcoholism interviewers, led by Yari Nuñez and Michelle Slivinsky, devoted substantial time and effort to phenotype the study sample. John Farrell, PhD, Section of Biomedical Genetics, Boston University School of Medicine, provided database management assistance. None of these individuals were compensated for this contribution to this study. Assistance with phenotype harmonization, genotype cleaning, and general study coordination was provided by the Gene Environment Association Studies Coordinating Center (supported by grant U01 HG004446 from the NIH). Assistance with data cleaning was provided by the National Center for Biotechnology Information.

References

1. United Nation Office on Drugs and Crime (UNODC) [Accessed February 18, 2016];Cannabis: A Short Review. https://www.unodc.org/documents/drug-prevention-and-treatment/cannabis_review.pdf. Published 2012.
2. Volkow ND, Compton WM, Weiss SR. Adverse health effects of marijuana use. N Engl J Med. 2014;371(9):879. [Abstract] [Google Scholar]
3. Filbey FM, Schacht JP, Myers US, Chavez RS, Hutchison KE. Marijuana craving in the brain. Proc Natl Acad Sci U S A. 2009;106(31):13016–13021. [Europe PMC free article] [Abstract] [Google Scholar]
4. Weinstein AM, Gorelick DA. Pharmacological treatment of cannabis dependence. Curr Pharm Des. 2011;17(14):1351–1358. [Europe PMC free article] [Abstract] [Google Scholar]
5. Wölfling K, Flor H, Grüsser SM. Psychophysiological responses to drug-associated stimuli in chronic heavy cannabis use. Eur J Neurosci. 2008;27(4):976–983. [Abstract] [Google Scholar]
6. Hasin DS, Saha TD, Kerridge BT, et al. Prevalence of marijuana use disorders in the United States between 2001–2002 and 2012–2013. JAMA Psychiatry. 2015;72(12):1235–1242. [Europe PMC free article] [Abstract] [Google Scholar]
7. Arseneault L, Cannon M, Witton J, Murray RM. Causal association between cannabis and psychosis: examination of the evidence. Br J Psychiatry. 2004;184:110–117. [Abstract] [Google Scholar]
8. Power RA, Verweij KJ, Zuhair M, et al. Genetic predisposition to schizophrenia associated with increased use of cannabis. Mol Psychiatry. 2014;19(11):1201–1204. [Europe PMC free article] [Abstract] [Google Scholar]
9. Bovasso GB. Cannabis abuse as a risk factor for depressive symptoms. Am J Psychiatry. 2001;158(12):2033–2037. [Abstract] [Google Scholar]
10. Lynskey MT, Glowinski AL, Todorov AA, et al. Major depressive disorder, suicidal ideation, and suicide attempt in twins discordant for cannabis dependence and early-onset cannabis use. Arch Gen Psychiatry. 2004;61(10):1026–1032. [Abstract] [Google Scholar]
11. Lynskey MT, Agrawal A, Henders A, Nelson EC, Madden PA, Martin NG. An Australian twin study of cannabis and other illicit drug use and misuse, and other psychopathology. Twin Res Hum Genet. 2012;15(5):631–641. [Europe PMC free article] [Abstract] [Google Scholar]
12. Kendler KS, Ohlsson H, Maes HH, Sundquist K, Lichtenstein P, Sundquist J. A population-based Swedish Twin and Sibling Study of cannabis, stimulant and sedative abuse in men. Drug Alcohol Depend. 2015;149:49–54. [Europe PMC free article] [Abstract] [Google Scholar]
13. Agrawal A, Lynskey MT, Bucholz KK, et al. DSM-5 cannabis use disorder: a phenotypic and genomic perspective. Drug Alcohol Depend. 2014;134:362–369. [Europe PMC free article] [Abstract] [Google Scholar]
14. Ehlers CL, Gilder DA, Gizer IR, Wilhelmsen KC. Heritability and a genome-wide linkage analysis of a type II/B cluster construct for cannabis dependence in an American Indian community. Addict Biol. 2009;14(3):338–348. [Europe PMC free article] [Abstract] [Google Scholar]
15. Stefansson H, Sigurdsson E, Steinthorsdottir V, et al. Neuregulin 1 and susceptibility to schizophrenia. Am J Hum Genet. 2002;71(4):877–892. [Europe PMC free article] [Abstract] [Google Scholar]
16. Han S, Yang BZ, Kranzler HR, et al. Linkage analysis followed by association show NRG1 associated with cannabis dependence in African Americans. Biol Psychiatry. 2012;72(8):637–644. [Europe PMC free article] [Abstract] [Google Scholar]
17. Verweij KJ, Vinkhuyzen AA, Benyamin B, et al. The genetic aetiology of cannabis use initiation: a meta-analysis of genome-wide association studies and a SNP-based heritability estimation. Addict Biol. 2013;18(5):846–850. [Europe PMC free article] [Abstract] [Google Scholar]
18. Agrawal A, Lynskey MT, Hinrichs A, et al. GENEVA Consortium. A genome-wide association study of DSM-IV cannabis dependence. Addict Biol. 2011;16(3):514–518. [Europe PMC free article] [Abstract] [Google Scholar]
19. Gelernter J, Kranzler HR, Sherva R, et al. Genome-wide association study of opioid dependence: multiple associations mapped to calcium and potassium pathways. Biol Psychiatry. 2014;76(1):66–74. [Europe PMC free article] [Abstract] [Google Scholar]
20. Gelernter J, Sherva R, Koesterer R, et al. Genome-wide association study of cocaine dependence and related traits: FAM53B identified as a risk gene. Mol Psychiatry. 2014;19(6):717–723. [Europe PMC free article] [Abstract] [Google Scholar]
21. Gelernter J, Kranzler HR, Sherva R, et al. Genome-wide association study of alcohol dependence: significant findings in African- and European-Americans including novel risk loci. Mol Psychiatry. 2014;19(1):41–49. [Europe PMC free article] [Abstract] [Google Scholar]
22. Gelernter J, Kranzler HR, Sherva R, et al. Genome-wide association study of nicotine dependence in American populations: identification of novel risk loci in both African-Americans and European-Americans. Biol Psychiatry. 2015;77(5):493–503. [Europe PMC free article] [Abstract] [Google Scholar]
23. Xie P, Kranzler HR, Yang C, Zhao H, Farrer LA, Gelernter J. Genome-wide association study identifies new susceptibility loci for posttraumatic stress disorder. Biol Psychiatry. 2013;74(9):656–663. [Europe PMC free article] [Abstract] [Google Scholar]
24. National Center for Biotechnology Information. [Accessed March 4, 2016];Study of Addiction: Genetics and Environment (SAGE) phs000092.v1.p1. http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000092.v1.p1.
25. Bierut LJ, Strickland JR, Thompson JR, Afful SE, Cottler LB. Drug use and dependence in cocaine dependent subjects, community-based individuals, and their siblings. Drug Alcohol Depend. 2008;95(1–2):14–22. [Europe PMC free article] [Abstract] [Google Scholar]
26. Edenberg HJ. The collaborative study on the genetics of alcoholism: an update. Alcohol Res Health. 2002;26(3):214–218. [Europe PMC free article] [Abstract] [Google Scholar]
27. Bierut LJ, Madden PA, Breslau N, et al. Novel genes identified in a high-density genome wide association study for nicotine dependence. Hum Mol Genet. 2007;16(1):24–35. [Europe PMC free article] [Abstract] [Google Scholar]
28. National Center for Biotechnology Information. [Accessed March 4, 2016];A Genome-wide Association Study of Heroin Dependence. phs000277.v1.p1. http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000277.v1.p1.
29. Maloney E, Degenhardt L, Darke S, Mattick RP, Nelson E. Suicidal behaviour and associated risk factors among opioid-dependent individuals: a case-control study. Addiction. 2007;102(12):1933–1941. [Abstract] [Google Scholar]
30. Pierucci-Lagha A, Gelernter J, Feinn R, et al. Diagnostic reliability of the Semi-Structured Assessment for Drug Dependence and Alcoholism (SSADDA) Drug Alcohol Depend. 2005;80(3):303–312. [Abstract] [Google Scholar]
31. Bucholz KK, Cadoret R, Cloninger CR, et al. A new, semi-structured psychiatric interview for use in genetic linkage studies: a report on the reliability of the SSAGA. J Stud Alcohol. 1994;55(2):149–158. [Abstract] [Google Scholar]
32. Zeger SL, Liang KY. Longitudinal data analysis for discrete and continuous outcomes. Biometrics. 1986;42(1):121–130. [Abstract] [Google Scholar]
33. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–2191. [Europe PMC free article] [Abstract] [Google Scholar]
34. Cross-Disorder Group of the Psychiatric Genomics Consortium. Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet. 2013;381(9875):1371–1379. [Europe PMC free article] [Abstract] [Google Scholar]
35. Chung D, Yang C, Li C, Gelernter J, Zhao H. GPA: a statistical approach to prioritizing GWAS results by integrating pleiotropy and annotation. PLoS Genet. 2014;10(11):e1004787. [Europe PMC free article] [Abstract] [Google Scholar]
36. Singh B, Chaudhuri TK. Role of C-reactive protein in schizophrenia: an overview. Psychiatry Res. 2014;216(2):277–285. [Abstract] [Google Scholar]
37. Rosenblat JD, Cha DS, Mansur RB, McIntyre RS. Inflamed moods: a review of the interactions between inflammation and mood disorders. Prog Neuropsychopharmacol Biol Psychiatry. 2014;53:23–34. [Abstract] [Google Scholar]
38. Iwata M, Ota KT, Duman RS. The inflammasome: pathways linking psychological stress, depression, and systemic illnesses. Brain Behav Immun. 2013;31:105–114. [Europe PMC free article] [Abstract] [Google Scholar]
39. Yuan JP, Zeng W, Huang GN, Worley PF, Muallem S. STIM1 heteromultimerizes TRPC channels to determine their function as store-operated channels. Nat Cell Biol. 2007;9(6):636–645. [Europe PMC free article] [Abstract] [Google Scholar]
40. Krapivinsky G, Krapivinsky L, Stotz SC, Manasian Y, Clapham DE. POST, partner of stromal interaction molecule 1 (STIM1), targets STIM1 to multiple transporters. Proc Natl Acad Sci U S A. 2011;108(48):19234–19239. [Europe PMC free article] [Abstract] [Google Scholar]
41. Kraus DM, Elliott GS, Chute H, et al. CSMD1 is a novel multiple domain complement-regulatory protein highly expressed in the central nervous system and epithelial tissues. J Immunol. 2006;176(7):4419–4430. [Abstract] [Google Scholar]
42. Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511(7510):421–427. [Europe PMC free article] [Abstract] [Google Scholar]
43. Mao Z, Bonni A, Xia F, Nadal-Vicens M, Greenberg ME. Neuronal activity-dependent cell survival mediated by transcription factor MEF2. Science. 1999;286(5440):785–790. [Abstract] [Google Scholar]
44. Chevaleyre V, Heifets BD, Kaeser PS, Südhof TC, Castillo PE. Endocannabinoid-mediated long-term plasticity requires cAMP/PKA signaling and RIM1alpha. Neuron. 2007;54(5):801–812. [published correction appears in Neuron. 2007 Jul 5;55(1):169] [Europe PMC free article] [Abstract] [Google Scholar]
45. Ho BC, Wassink TH, Ziebell S, Andreasen NC. Cannabinoid receptor 1 gene polymorphisms and marijuana misuse interactions on white matter and cognitive deficits in schizophrenia. Schizophr Res. 2011;128(1–3):66–75. [Europe PMC free article] [Abstract] [Google Scholar]
46. Eggan SM, Lazarus MS, Stoyak SR, et al. Cortical glutamic acid decarboxylase 67 deficiency results in lower cannabinoid 1 receptor messenger RNA expression: implications for schizophrenia. Biol Psychiatry. 2012;71(2):114–119. [Europe PMC free article] [Abstract] [Google Scholar]
47. Mathews RR, Hall WD, Gartner CE. Depression and psychological distress in tobacco smokers and people with cannabis dependence in the National Survey of Mental Health and Wellbeing. Med J Aust. 2011;195(3):S12–S15. [Abstract] [Google Scholar]
48. D’Souza DC, Sewell RA, Ranganathan M. Cannabis and psychosis/schizophrenia: human studies. Eur Arch Psychiatry Clin Neurosci. 2009;259(7):413–431. [Europe PMC free article] [Abstract] [Google Scholar]
49. Leweke FM, Koethe D. Cannabis and psychiatric disorders: it is not only addiction. Addict Biol. 2008;13(2):264–275. [Abstract] [Google Scholar]
50. Issa BA. Cannabis use and psychiatric disorders. Neurosciences (Riyadh) 2007;12(3):265–266. [Abstract] [Google Scholar]

Full text links 

Read article at publisher's site: https://doi.org/10.1001/jamapsychiatry.2016.0036

Read article for free, from open access legal sources, via Unpaywall: https://jamanetwork.com/journals/jamapsychiatry/articlepdf/2504223/yoi160002.pdfUnpaywall

Citations & impact 

Impact metrics

100
Citations
Jump to Citations
2
Data citations
Jump to Data

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/6391451
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/6391451

Article citations

Go to all (100) article citations

Data 

Data behind the article

This data has been text mined from the article, or deposited into data resources.

BioStudies: supplemental material and supporting data

dbGaP - The database of Genotypes and Phenotypes

Diseases (4)

Ensembl Genome Browser

SNPs (Showing 29 of 29)

Similar Articles 

To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.

Funding 

Funders who supported this work.

NCATS NIH HHS (1)

NCI NIH HHS (2)

NHGRI NIH HHS (4)

NHLBI NIH HHS (3)

NIAAA NIH HHS (3)

NIDA NIH HHS (5)