Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Associations Between Cannabis Use, Polygenic Liability for Schizophrenia, and Cannabis-related Experiences in a Sample of Cannabis Users

BACKGROUND AND HYPOTHESIS

Risk for cannabis use and schizophrenia is influenced in part by genetic factors, and there is evidence that genetic risk for schizophrenia is associated with subclinical psychotic-like experiences (PLEs). Few studies to date have examined whether genetic risk for schizophrenia is associated with cannabis-related PLEs.

STUDY DESIGN

We tested whether measures of cannabis involvement and polygenic risk scores (PRS) for schizophrenia were associated with self-reported cannabis-related experiences in a sample ascertained for alcohol use disorders (AUDs), the Collaborative Study on the Genetics of Alcoholism (COGA). We analyzed 4832 subjects (3128 of European ancestry and 1704 of African ancestry; 42% female; 74% meeting lifetime criteria for an AUD).

STUDY RESULTS

Cannabis use disorder (CUD) was prevalent in this analytic sample (70%), with 40% classified as mild, 25% as moderate, and 35% as severe. Polygenic risk for schizophrenia was positively associated with cannabis-related paranoia, feeling depressed or anhedonia, social withdrawal, and cognitive difficulties, even when controlling for duration of daily cannabis use, CUD, and age at first cannabis use. The schizophrenia PRS was most robustly associated with cannabis-related cognitive difficulties (β = 0.22, SE = 0.04, P = 5.2e-7). In an independent replication sample (N = 1446), associations between the schizophrenia PRS and cannabis-related experiences were in the expected direction and not statistically different in magnitude from those in the COGA sample.

CONCLUSIONS

Among individuals who regularly use cannabis, genetic liability for schizophrenia-even in those without clinical features-may increase the likelihood of reporting unusual experiences related to cannabis use.

Cnih3 Deletion Dysregulates Dorsal Hippocampal Transcription across the Estrous Cycle

In females, the hippocampus, a critical brain region for coordination of learning, memory, and behavior, displays altered physiology and behavioral output across the estrous or menstrual cycle. However, the molecular effectors and cell types underlying these observed cyclic changes have only been partially characterized to date. Recently, profiling of mice null for the AMPA receptor trafficking gene Cnih3 have demonstrated estrous-dependent phenotypes in dorsal hippocampal synaptic plasticity, composition, and learning/memory. We therefore profiled dorsal hippocampal transcriptomes from female mice in each estrous cycle stage, and contrasted it with that of males, across wild-type (WT) and Cnih3 mutants. In wild types, we identified only subtle differences in gene expression between the sexes, while comparing estrous stages to one another revealed up to >1000 differentially expressed genes (DEGs). These estrous-responsive genes are especially enriched in gene markers of oligodendrocytes and the dentate gyrus, and in functional gene sets relating to estrogen response, potassium channels, and synaptic gene splicing. Surprisingly, Cnih3 knock-outs (KOs) showed far broader transcriptomic differences between estrous cycle stages and males. Moreover, Cnih3 knock-out drove subtle but extensive expression changes accentuating sex differential expression at diestrus and estrus. Altogether, our profiling highlights cell types and molecular systems potentially impacted by estrous-specific gene expression patterns in the adult dorsal hippocampus, enabling mechanistic hypothesis generation for future studies of sex-differential neuropsychiatric function and dysfunction. Moreover, these findings suggest an unrecognized role of Cnih3 in buffering against transcriptional effects of estrous, providing a candidate molecular mechanism to explain estrous-dependent phenotypes observed with Cnih3 loss.

Multi-trait genome-wide association study of opioid addiction: OPRM1 and beyond

Opioid addiction (OA) is moderately heritable, yet only rs1799971, the A118G variant in OPRM1, has been identified as a genome-wide significant association with OA and independently replicated. We applied genomic structural equation modeling to conduct a GWAS of the new Genetics of Opioid Addiction Consortium (GENOA) data together with published studies (Psychiatric Genomics Consortium, Million Veteran Program, and Partners Health), comprising 23,367 cases and effective sample size of 88,114 individuals of European ancestry. Genetic correlations among the various OA phenotypes were uniformly high (rg > 0.9). We observed the strongest evidence to date for OPRM1: lead SNP rs9478500 (p = 2.56 × 10-9). Gene-based analyses identified novel genome-wide significant associations with PPP6C and FURIN. Variants within these loci appear to be pleiotropic for addiction and related traits.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

Exploring the genetic overlap of suicide-related behaviors and substance use disorders

Suicide-related behaviors are heterogeneous and transdiagnostic, and may demonstrate varying levels of genetic overlap with different substance use disorders (SUDs). We used linkage disequilibrium score regression, genomic structural equation models, and Mendelian randomization to examine the genetic relationships between several SUDs and suicide-related behaviors. Our analyses incorporated summary statistics from the largest genome-wide association studies (GWAS) of problematic alcohol use, the Fagerström test for nicotine dependence, cannabis use disorder, and opioid use disorder (Ns ranging from 46,213-435,563) and GWAS of ever self-harmed, suicide attempt, and suicide death (Ns ranging from 18,223-117,733). We also accounted for genetic liability to depression (N = 500,199) and risk tolerance (N = 315,894). Suicide-related behaviors were significantly genetically correlated with each other and each SUD, but there was little evidence of causal relationships between the traits. Simultaneously correlating a common SUD factor with each specific suicide indicator while controlling for depression and risk tolerance revealed significant, positive genetic correlations between the SUD factor and suicide-related behaviors (rg  = 0.26-0.45, SE = 0.08-0.09). In the model, depression's association with suicide death (β = 0.42, SE = 0.06) was weaker compared to ever-self harmed and suicide attempt (β = 0.58, SE = 0.05 and β = 0.50, SE = 0.06, respectively). We identify a general level of genetic overlap between SUDs and suicide-related behaviors, which is independent of depression and risk tolerance. Additionally, our findings suggest that genetic and behavioral contributions to suicide death may somewhat differ from nonlethal suicide-related behaviors.

Sex Differences in the Role of CNIH3 on Spatial Memory and Synaptic Plasticity

BACKGROUND

CNIH3 is an AMPA receptor (AMPAR) auxiliary protein prominently expressed in the dorsal hippocampus (dHPC), a region that plays a critical role in spatial memory and synaptic plasticity. However, the effects of CNIH3 on AMPAR-dependent synaptic function and behavior have not been investigated.

METHODS

We assessed a gain-of-function model of Cnih3 overexpression in the dHPC and generated and characterized a line of Cnih3-/- C57BL/6 mice. We assessed spatial memory through behavioral assays, protein levels of AMPAR subunits and synaptic proteins by immunoblotting, and long-term potentiation in electrophysiological recordings. We also utilized a super-resolution imaging workflow, SEQUIN (Synaptic Evaluation and Quantification by Imaging of Nanostructure), for analysis of nanoscale synaptic connectivity in the dHPC.

RESULTS

Overexpression of Cnih3 in the dHPC improved short-term spatial memory in female mice but not in male mice. Cnih3-/- female mice exhibited weakened short-term spatial memory, reduced dHPC synapse density, enhanced expression of calcium-impermeable AMPAR (GluA2-containing) subunits in synaptosomes, and attenuated long-term potentiation maintenance compared with Cnih3+/+ control mice; Cnih3-/- males were unaffected. Further investigation revealed that deficiencies in spatial memory and changes in AMPAR composition and synaptic plasticity were most pronounced during the metestrus phase of the estrous cycle in female Cnih3-/- mice.

CONCLUSIONS

This study identified a novel effect of sex and estrous on CNIH3's role in spatial memory and synaptic plasticity. Manipulation of CNIH3 unmasked sexually dimorphic effects on spatial memory, synaptic function, AMPAR composition, and hippocampal plasticity. These findings reinforce the importance of considering sex as a biological variable in studies of memory and hippocampal synaptic function.

Using Monozygotic Twins to Dissect Common Genes in Posttraumatic Stress Disorder and Migraine

Epigenetic mechanisms have been associated with genes involved in Posttraumatic stress disorder (PTSD). PTSD often co-occurs with other health conditions such as depression, cardiovascular disorder and respiratory illnesses. PTSD and migraine have previously been reported to be symptomatically positively correlated with each other, but little is known about the genes involved. The aim of this study was to understand the comorbidity between PTSD and migraine using a monozygotic twin disease discordant study design in six pairs of monozygotic twins discordant for PTSD and 15 pairs of monozygotic twins discordant for migraine. DNA from peripheral blood was run on Illumina EPIC arrays and analyzed. Multiple testing correction was performed using the Bonferroni method and 10% false discovery rate (FDR). We validated 11 candidate genes previously associated with PTSD including DOCK2, DICER1, and ADCYAP1. In the epigenome-wide scan, seven novel CpGs were significantly associated with PTSD within/near IL37, WNT3, ADNP2, HTT, SLFN11, and NQO2, with all CpGs except the IL37 CpG hypermethylated in PTSD. These results were significantly enriched for genes whose DNA methylation was previously associated with migraine (p-value = 0.036). At 10% FDR, 132 CpGs in 99 genes associated with PTSD were also associated with migraine in the migraine twin samples. Genes associated with PTSD were overrepresented in vascular smooth muscle, axon guidance and oxytocin signaling pathways, while genes associated with both PTSD and migraine were enriched for AMPK signaling and longevity regulating pathways. In conclusion, these results suggest that common genes and pathways are likely involved in PTSD and migraine, explaining at least in part the co-morbidity between the two disorders.

A large-scale genome-wide association study meta-analysis of cannabis use disorder

BACKGROUND

Variation in liability to cannabis use disorder has a strong genetic component (estimated twin and family heritability about 50-70%) and is associated with negative outcomes, including increased risk of psychopathology. The aim of the study was to conduct a large genome-wide association study (GWAS) to identify novel genetic variants associated with cannabis use disorder.

METHODS

To conduct this GWAS meta-analysis of cannabis use disorder and identify associations with genetic loci, we used samples from the Psychiatric Genomics Consortium Substance Use Disorders working group, iPSYCH, and deCODE (20 916 case samples, 363 116 control samples in total), contrasting cannabis use disorder cases with controls. To examine the genetic overlap between cannabis use disorder and 22 traits of interest (chosen because of previously published phenotypic correlations [eg, psychiatric disorders] or hypothesised associations [eg, chronotype] with cannabis use disorder), we used linkage disequilibrium score regression to calculate genetic correlations.

FINDINGS

We identified two genome-wide significant loci: a novel chromosome 7 locus (FOXP2, lead single-nucleotide polymorphism [SNP] rs7783012; odds ratio [OR] 1·11, 95% CI 1·07-1·15, p=1·84 × 10-9) and the previously identified chromosome 8 locus (near CHRNA2 and EPHX2, lead SNP rs4732724; OR 0·89, 95% CI 0·86-0·93, p=6·46 × 10-9). Cannabis use disorder and cannabis use were genetically correlated (rg 0·50, p=1·50 × 10-21), but they showed significantly different genetic correlations with 12 of the 22 traits we tested, suggesting at least partially different genetic underpinnings of cannabis use and cannabis use disorder. Cannabis use disorder was positively genetically correlated with other psychopathology, including ADHD, major depression, and schizophrenia.

INTERPRETATION

These findings support the theory that cannabis use disorder has shared genetic liability with other psychopathology, and there is a distinction between genetic liability to cannabis use and cannabis use disorder.

FUNDING

National Institute of Mental Health; National Institute on Alcohol Abuse and Alcoholism; National Institute on Drug Abuse; Center for Genomics and Personalized Medicine and the Centre for Integrative Sequencing; The European Commission, Horizon 2020; National Institute of Child Health and Human Development; Health Research Council of New Zealand; National Institute on Aging; Wellcome Trust Case Control Consortium; UK Research and Innovation Medical Research Council (UKRI MRC); The Brain & Behavior Research Foundation; National Institute on Deafness and Other Communication Disorders; Substance Abuse and Mental Health Services Administration (SAMHSA); National Institute of Biomedical Imaging and Bioengineering; National Health and Medical Research Council (NHMRC) Australia; Tobacco-Related Disease Research Program of the University of California; Families for Borderline Personality Disorder Research (Beth and Rob Elliott) 2018 NARSAD Young Investigator Grant; The National Child Health Research Foundation (Cure Kids); The Canterbury Medical Research Foundation; The New Zealand Lottery Grants Board; The University of Otago; The Carney Centre for Pharmacogenomics; The James Hume Bequest Fund; National Institutes of Health: Genes, Environment and Health Initiative; National Institutes of Health; National Cancer Institute; The William T Grant Foundation; Australian Research Council; The Virginia Tobacco Settlement Foundation; The VISN 1 and VISN 4 Mental Illness Research, Education, and Clinical Centers of the US Department of Veterans Affairs; The 5th Framework Programme (FP-5) GenomEUtwin Project; The Lundbeck Foundation; NIH-funded Shared Instrumentation Grant S10RR025141; Clinical Translational Sciences Award grants; National Institute of Neurological Disorders and Stroke; National Heart, Lung, and Blood Institute; National Institute of General Medical Sciences.

Leveraging genome-wide data to investigate differences between opioid use vs. opioid dependence in 41,176 individuals from the Psychiatric Genomics Consortium

To provide insights into the biology of opioid dependence (OD) and opioid use (i.e., exposure, OE), we completed a genome-wide analysis comparing 4503 OD cases, 4173 opioid-exposed controls, and 32,500 opioid-unexposed controls, including participants of European and African descent (EUR and AFR, respectively). Among the variants identified, rs9291211 was associated with OE (exposed vs. unexposed controls; EUR z = -5.39, p = 7.2 × 10-8). This variant regulates the transcriptomic profiles of SLC30A9 and BEND4 in multiple brain tissues and was previously associated with depression, alcohol consumption, and neuroticism. A phenome-wide scan of rs9291211 in the UK Biobank (N > 360,000) found association of this variant with propensity to use dietary supplements (p = 1.68 × 10-8). With respect to the same OE phenotype in the gene-based analysis, we identified SDCCAG8 (EUR + AFR z = 4.69, p = 10-6), which was previously associated with educational attainment, risk-taking behaviors, and schizophrenia. In addition, rs201123820 showed a genome-wide significant difference between OD cases and unexposed controls (AFR z = 5.55, p = 2.9 × 10-8) and a significant association with musculoskeletal disorders in the UK Biobank (p = 4.88 × 10-7). A polygenic risk score (PRS) based on a GWAS of risk-tolerance (n = 466,571) was positively associated with OD (OD vs. unexposed controls, p = 8.1 × 10-5; OD cases vs. exposed controls, p = 0.054) and OE (exposed vs. unexposed controls, p = 3.6 × 10-5). A PRS based on a GWAS of neuroticism (n = 390,278) was positively associated with OD (OD vs. unexposed controls, p = 3.2 × 10-5; OD vs. exposed controls, p = 0.002) but not with OE (p = 0.67). Our analyses highlight the difference between dependence and exposure and the importance of considering the definition of controls in studies of addiction.

Retrospectively assessed subjective effects of initial opioid use differ between opioid misusers with opioid use disorder (OUD) and those who never progressed to OUD: Data from a pilot and a replication sample

Attempts to identify opioid users with increased risk of escalating to opioid use disorder (OUD) have had limited success. Retrospectively assessed subjective effects of initial opioid misuse were compared in a pilot sample of opioid misusers (nonmedical use ≤60 times lifetime) who had never met criteria for OUD (N = 14) and heroin-addicted individuals in treatment for OUD (N = 15). Relative to opioid misusers without a lifetime OUD diagnosis, individuals with OUD reported greater euphoria and other positive emotions, activation, pruritus, and internalizing symptoms. Consistent with these findings, proxy Addiction Research Center Inventory (ARCI) Amphetamine Group, and Morphine Benzedrine Group scale mean item scores were significantly higher in those with OUD. Replication was attempted in opioid misusers with (N = 25) and without OUD (N = 25) who were assessed as part of an ongoing genetic study. We observed similar significant between-group differences in individual subjective effect items and ARCI scale mean item scores in the replication sample. We, thus confirm findings from prior reports that retrospectively assessed subjective responses to initial opioid exposure differ significantly between opioid users who do, and do not, progress to OUD. Our report extends these findings in comparisons limited to opioid misusers. Additional research will be necessary to examine prospectively whether the assessment of subjective effects after initial use has predictive utility in the identification of individuals more likely to progress to OUD.

Single cell transcriptomics reveals opioid usage evokes widespread suppression of antiviral gene program

Chronic opioid usage not only causes addiction behavior through the central nervous system, but also modulates the peripheral immune system. However, how opioid impacts the immune system is still barely characterized systematically. In order to understand the immune modulatory effect of opioids in an unbiased way, here we perform single-cell RNA sequencing (scRNA-seq) of peripheral blood mononuclear cells from opioid-dependent individuals and controls to show that chronic opioid usage evokes widespread suppression of antiviral gene program in naive monocytes, as well as in multiple immune cell types upon stimulation with the pathogen component lipopolysaccharide. Furthermore, scRNA-seq reveals the same phenomenon after a short in vitro morphine treatment. These findings indicate that both acute and chronic opioid exposure may be harmful to our immune system by suppressing the antiviral gene program. Our results suggest that further characterization of the immune modulatory effects of opioid is critical to ensure the safety of clinical opioids.

Over 100 million of opioid prescriptions are issued yearly in the USA alone, but the impact of opioid use on the immune system is barely characterized. Here the authors report antiviral immune response is blunted in several types of blood cells from opioid-dependent individuals, and when healthy donor cells are exposed to morphine in a dish.

International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci

The risk of posttraumatic stress disorder (PTSD) following trauma is heritable, but robust common variants have yet to be identified. In a multi-ethnic cohort including over 30,000 PTSD cases and 170,000 controls we conduct a genome-wide association study of PTSD. We demonstrate SNP-based heritability estimates of 5-20%, varying by sex. Three genome-wide significant loci are identified, 2 in European and 1 in African-ancestry analyses. Analyses stratified by sex implicate 3 additional loci in men. Along with other novel genes and non-coding RNAs, a Parkinson's disease gene involved in dopamine regulation, PARK2, is associated with PTSD. Finally, we demonstrate that polygenic risk for PTSD is significantly predictive of re-experiencing symptoms in the Million Veteran Program dataset, although specific loci did not replicate. These results demonstrate the role of genetic variation in the biology of risk for PTSD and highlight the necessity of conducting sex-stratified analyses and expanding GWAS beyond European ancestry populations.

Childhood Trauma and Two Stages of Alcohol Use in African American and European American Women: Findings from a Female Twin Sample

The current investigation assessed for moderating effects of childhood trauma on genetic and environmental contributions to timing of alcohol use initiation and alcohol use disorder in African American (AA) and European American (EA) women. Data were drawn from diagnostic telephone interviews conducted with 3786 participants (14.6% AA) in a longitudinal female twin study. Childhood trauma was defined alternately as child maltreatment and more broadly to include other events (e.g., witnessing violence). Phenotypic associations between childhood trauma and alcohol outcomes were estimated using logistic regression analyses. Twin modeling was conducted to test for moderating effects of childhood trauma on the contributions of genetic and environmental factors to timing of initiation and alcohol use disorder. Under both definitions, childhood trauma was associated with early initiation (relative risk ratios: 1.90, 1.72) and alcohol use disorder (odds ratios: 1.92, 1.76). Yet gene by environment effects were observed only for child maltreatment and timing of initiation in EA women, with heritable influences less prominent in those who had experienced child maltreatment (0.35, 95% CI: 0.05-0.66 vs. 0.52, 95% CI: 0.30-0.73). We found more similarities than differences in the association of childhood trauma with alcohol outcomes across racial/ethnic groups, trauma type, and stages of alcohol use. However, findings suggest that the relative contribution of genetic factors to alcohol outcomes differs by childhood maltreatment history in EA women specifically in the earliest stage of alcohol use.

Common genetic contributions to high-risk trauma exposure and self-injurious thoughts and behaviors

BACKGROUND

Prior research has documented shared heritable contributions to non-suicidal self-injury (NSSI) and suicidal ideation (SI) as well as NSSI and suicide attempt (SA). In addition, trauma exposure has been implicated in risk for NSSI and suicide. Genetically informative studies are needed to determine common sources of liability to all three self-injurious thoughts and behaviors, and to clarify the nature of their associations with traumatic experiences.

METHODS

Multivariate biometric modeling was conducted using data from 9526 twins [59% female, mean age = 31.7 years (range 24-42)] from two cohorts of the Australian Twin Registry, some of whom also participated in the Childhood Trauma Study and the Nicotine Addiction Genetics Project.

RESULTS

The prevalences of high-risk trauma exposure (HRT), NSSI, SI, and SA were 24.4, 5.6, 27.1, and 4.6%, respectively. All phenotypes were moderately to highly correlated. Genetic influences on self-injurious thoughts and behaviors and HRT were significant and highly correlated among men [rG = 0.59, 95% confidence interval (CI) (0.37-0.81)] and women [rG = 0.56 (0.49-0.63)]. Unique environmental influences were modestly correlated in women [rE = 0.23 (0.01-0.45)], suggesting that high-risk trauma may confer some direct risk for self-injurious thoughts and behaviors among females.

CONCLUSIONS

Individuals engaging in NSSI are at increased risk for suicide, and common heritable factors contribute to these associations. Preventing trauma exposure may help to mitigate risk for self-harm and suicide, either directly or indirectly via reductions in liability to psychopathology more broadly. In addition, targeting pre-existing vulnerability factors could significantly reduce risk for life-threatening behaviors among those who have experienced trauma.

Publisher Correction: Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways

The original version of this Article contained errors in the spelling of the authors Fan Liu and M. Arfan Ikram, which were incorrectly given as Fan Lui and Arfan M. Ikram. In addition, the original version of this Article also contained errors in the author affiliations which are detailed in the associated Publisher Correction.

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.

Transancestral GWAS of alcohol dependence reveals common genetic underpinnings with psychiatric disorders

Liability to alcohol dependence (AD) is heritable, but little is known about its complex polygenic architecture or its genetic relationship with other disorders. To discover loci associated with AD and characterize the relationship between AD and other psychiatric and behavioral outcomes, we carried out the largest genome-wide association study to date of DSM-IV-diagnosed AD. Genome-wide data on 14,904 individuals with AD and 37,944 controls from 28 case-control and family-based studies were meta-analyzed, stratified by genetic ancestry (European, n = 46,568; African, n = 6,280). Independent, genome-wide significant effects of different ADH1B variants were identified in European (rs1229984; P = 9.8 × 10-13) and African ancestries (rs2066702; P = 2.2 × 10-9). Significant genetic correlations were observed with 17 phenotypes, including schizophrenia, attention deficit-hyperactivity disorder, depression, and use of cigarettes and cannabis. The genetic underpinnings of AD only partially overlap with those for alcohol consumption, underscoring the genetic distinction between pathological and nonpathological drinking behaviors.

Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways

The total number of acquired melanocytic nevi on the skin is strongly correlated with melanoma risk. Here we report a meta-analysis of 11 nevus GWAS from Australia, Netherlands, UK, and USA comprising 52,506 individuals. We confirm known loci including MTAP, PLA2G6, and IRF4, and detect novel SNPs in KITLG and a region of 9q32. In a bivariate analysis combining the nevus results with a recent melanoma GWAS meta-analysis (12,874 cases, 23,203 controls), SNPs near GPRC5A, CYP1B1, PPARGC1B, HDAC4, FAM208B, DOCK8, and SYNE2 reached global significance, and other loci, including MIR146A and OBFC1, reached a suggestive level. Overall, we conclude that most nevus genes affect melanoma risk (KITLG an exception), while many melanoma risk loci do not alter nevus count. For example, variants in TERC and OBFC1 affect both traits, but other telomere length maintenance genes seem to affect melanoma risk only. Our findings implicate multiple pathways in nevogenesis.

Predictors of transitions across stages of heroin use and dependence prior to treatment-seeking among people in treatment for opioid dependence

BACKGROUND AND AIMS

Little is known about transition pathways among heroin users prior to treatment. This study examined the demographic and clinical predictors of transition speed from heroin use, to dependence, to first treatment episode.

METHODS

1149 heroin-dependent participants recruited from opioid agonist treatment clinics in Sydney, Australia, underwent a structured interview. Age of onset (AOO) was collected for heroin use, dependence and treatment-seeking, childhood maltreatment, psychiatric history and other substance dependence. Discrete-time survival analyses modelled years from onset of use to dependence, and from dependence to treatment-seeking, including demographic and clinical covariates.

FINDINGS

Median AOO for first heroin use, dependence and treatment-seeking was 18 years (inter-quartile range, or IQR = 6), 21 years (IQR = 7), and 24 years (IQR = 10) respectively. In adjusted models, younger birth cohorts (vs. born <1960), greater childhood maltreatment and later AAO of first heroin use were associated with more rapid transitions from heroin use to dependence. Living independently, parental violence, and alcohol dependence were associated with slower transitions. Earlier treatment-seeking was associated with younger birth cohorts, having dependent children and later AOO of dependence. Delayed treatment-seeking was associated with <10 years school education, living independently, depression and alcohol dependence.

CONCLUSIONS

In this treatment sample, onset of heroin use occurred during late adolescence, suggesting the need for targeted interventions in mid-adolescence. Transitions to heroin dependence, then treatment-seeking, occurred during early adulthood. Rapid transitions from use to dependence were associated with younger birth cohorts, greater exposure to childhood maltreatment, and later onset of use.

Genome-wide association study identifies a novel locus for cannabis dependence

Despite moderate heritability, only one study has identified genome-wide significant loci for cannabis-related phenotypes. We conducted meta-analyses of genome-wide association study data on 2080 cannabis-dependent cases and 6435 cannabis-exposed controls of European descent. A cluster of correlated single-nucleotide polymorphisms (SNPs) in a novel region on chromosome 10 was genome-wide significant (lowest P=1.3E-8). Among the SNPs, rs1409568 showed enrichment for H3K4me1 and H3K427ac marks, suggesting its role as an enhancer in addiction-relevant brain regions, such as the dorsolateral prefrontal cortex and the angular and cingulate gyri. This SNP is also predicted to modify binding scores for several transcription factors. We found modest evidence for replication for rs1409568 in an independent cohort of African American (896 cases and 1591 controls; P=0.03) but not European American (EA; 781 cases and 1905 controls) participants. The combined meta-analysis (3757 cases and 9931 controls) indicated trend-level significance for rs1409568 (P=2.85E-7). No genome-wide significant loci emerged for cannabis dependence criterion count (n=8050). There was also evidence that the minor allele of rs1409568 was associated with a 2.1% increase in right hippocampal volume in an independent sample of 430 EA college students (fwe-P=0.008). The identification and characterization of genome-wide significant loci for cannabis dependence is among the first steps toward understanding the biological contributions to the etiology of this psychiatric disorder, which appears to be rising in some developed nations.

Phenotypic and familial associations between childhood maltreatment and cannabis initiation and problems in young adult European-American and African-American women

BACKGROUND

Childhood maltreatment is a known risk factor for cannabis initiation and problem use, but the extent to which this association is attributable to shared familial influences is unknown. We estimate the magnitude of associations between childhood maltreatment, timing of cannabis initiation, and cannabis-related problems, in European-American (EA) and African-American (AA) women, and parse the relative influence of additive genetic (A), shared environmental (C), and individual-specific environmental (E) factors on these constructs and their covariation.

METHODS

Data were from diagnostic telephone interviews conducted with 3786 participants (14.6% AA) in a population-based study of female twins. Logistic regression analyses and twin modeling were used to test for associations, and estimate the relative contributions of genetic and environmental influences to childhood maltreatment and cannabis outcomes and their covariation.

RESULTS

Maltreatment was significantly associated with increased likelihood of cannabis initiation before age 15 among EAs (OR=6.33) and AAs (OR=3.93), but with increased likelihood of later initiation among EAs only (OR=1.68). Maltreatment was associated with cannabis problems among both groups (EA OR=2.32; AA OR=2.03). Among EA women, the covariation between maltreatment and cannabis outcomes was primarily attributable to familial environment (rC=0.67-0.70); among AAs, only individual-specific environment contributed (rE=0.37-0.40).

CONCLUSION

Childhood maltreatment is a major contributor to early initiation of cannabis as well as progression to cannabis problems in both AA and EA women. Distinctions by race/ethnicity are not in the relative contribution of genetic factors, but rather in the type of environmental influences that contribute to stages of cannabis involvement.

Major depressive disorder, suicidal thoughts and behaviours, and cannabis involvement in discordant twins: a retrospective cohort study

BACKGROUND

Early and frequent cannabis use are associated with an increased likelihood of major depressive disorder (MDD) as well as suicidal thoughts and behaviours. We identify associations between aspects of cannabis use, MDD, and suicidal thoughts and behaviours and examine whether such associations persist after accounting for those predisposing factors, including genetic liability and early family environment, that are shared by identical twins who are discordant for cannabis exposure. Any residual association in such identical pairs might be indicative of individual-specific pathways that might be of a causal nature.

METHODS

We did a logistic regression analysis of cannabis use from retrospective data on same-sex male and female twin pairs drawn from 3 studies that had recruited twins from the Australian Twin Registry, 1992-93 (sample 1), 1996-2000 (sample 2), and 2005-09 (sample 3). We studied associations between early use and frequent use of cannabis and MDD, suicidal ideation (ever and persistent), and suicide plan and attempt in the full sample as well as in pairs of monozygotic and dizygotic twins that were discordant for each measure of cannabis involvement at a single timepoint. Significant monozygotic associations were further adjusted for covariates, such as early alcohol or nicotine use, early dysphoric or anhedonic mood, conduct disorder, and childhood sexual abuse. Interactions between each cannabis measure and sex, sample or study effects, and birth year category were also examined as covariates.

FINDINGS

In 13 986 twins (6181 monozygotic and 7805 dizygotic), cannabis use ranged from 1345 (30·4%) of 4432 people in sample 1 to 2275 (69·0%) of 3299 in sample 3. Mean age of first cannabis use ranged from 17·9 years (SD 3·3) in sample 3 to 21·1 years (5·2) in sample 1, and frequent use (≥100 times) was reported by 214 (15·9%) of 1345 users in sample 1 and 499 (21·9%) of 2275 in sample 3. The prevalence of suicidal ideation ranged from 1102 (24·9%) of 4432 people in sample 1 to 1644 (26·3%) of 6255 people in sample 2 and 865 (26·2%) of 3299 people in sample 3. Prevalence of MDD ranged from 901 (20·3%) people in sample 1 to 1773 (28·3%) in sample 2. The monozygotic twin who used cannabis frequently was more likely to report MDD (odds ratio 1·98, 95% CI 1·11-3·53) and suicidal ideation (2·47, 1·19-5·10) compared with their identical twin who had used cannabis less frequently, even after adjustment for covariates. For early cannabis use, the monozygotic point estimate was not significant but could be equated to the significant dizygotic estimate, suggesting a possible association with suicidal ideation.

INTERPRETATION

The increased likelihood of MDD and suicidal ideation in frequent cannabis users cannot be solely attributed to common predisposing factors.

FUNDING

National Institute on Drug Abuse, National Institutes of Health, Australian National Health and Medical Research Council.