Genome-wide association meta-analysis of cocaine dependence: shared genetics with comorbid conditions

Cocaine dependence is a complex psychiatric disorder that is highly comorbid with other psychiatric traits. Twin and adoption studies suggest that genetic variants contribute substantially to cocaine dependence susceptibility, which has an estimated heritability of 65-79%. Here we performed a meta-analysis of genome-wide association studies of cocaine dependence using four datasets from the dbGaP repository (2085 cases and 4293 controls, all of them selected by their European ancestry). Although no genome-wide significant hits were found in the SNP-based analysis, the gene-based analysis identified HIST1H2BD as associated with cocaine-dependence (10% FDR). This gene is located in a region on chromosome 6 enriched in histone-related genes, previously associated with schizophrenia (SCZ). Furthermore, we performed LD Score regression analysis with comorbid conditions and found significant genetic correlations between cocaine dependence and SCZ, ADHD, major depressive disorder (MDD) and risk taking. We also found, through polygenic risk score analysis, that all tested phenotypes are significantly associated with cocaine dependence status: SCZ (R2 = 2.28%; P = 1.21e-26), ADHD (R2 = 1.39%; P = 4.5e-17), risk taking (R2 = 0.60%; P = 2.7e-08), MDD (R2 = 1.21%; P = 4.35e-15), children's aggressive behavior (R2 = 0.3%; P = 8.8e-05) and antisocial behavior (R2 = 1.33%; P = 2.2e-16). To our knowledge, this is the largest reported cocaine dependence GWAS meta-analysis in European-ancestry individuals. We identified suggestive associations in regions that may be related to cocaine dependence and found evidence for shared genetic risk factors between cocaine dependence and several comorbid psychiatric traits. However, the sample size is limited and further studies are needed to confirm these results

MIR-9, miR-153 and miR-124 are down-regulated by acute exposure to cocaine in adopaminergic cell model and may contribute to cocaine dependence

Cocaine is one of the most used psychostimulant drugs worldwide. MicroRNAs are post-transcriptional regulators of gene expression that are highly expressed in brain, and several studies have shown that cocaine can alter their expression. In a previous study, we identified several protein-coding genes that are differentially expressed in a dopaminergic neuron-like model after an acute exposure to cocaine. Now, we used the prediction tool WebGestalt to identify miRNA molecules potentially involved in the regulation of these genes. Using the same cellular model, we found that seven of these miRNAs are down-regulated by cocaine: miR-124-3p, miR-124-5p, miR-137, miR-101-3p, miR-9-5p, miR-369-3p and miR-153-3p, the last three not previously related to cocaine. Furthermore, we found that three of the miRNA genes that are differentially expressed in our model (hsa-miR-9-1, hsa-miR-153-1 and hsa-miR-124-3) are nominally associated with cocaine dependence in a case-control study (2,085 cases and 4,293 controls). In summary, we highlighted novel miRNAs that may be involved in those cocaine-induced changes of gene expression that underlie addiction. Moreover, we identified genetic variants that contribute to cocaine dependence in three of these miRNA genes, supporting the idea that genes differentially expressed under cocaine may play an important role in the susceptibility to cocaine dependence

The pleiotropic contribution of genes in dopaminergic and serotonergic pathways to addiction and related behavioral traits

Introduction: Co-occurrence of substance use disorders (SUD) and other behavioral conditions, such as stress-related, aggressive or risk-taking behaviors, in the same individual has been frequently described. As dopamine (DA) and serotonin (5-HT) have been previously identified as key neurotransmitters for some of these phenotypes, we explored the genetic contribution of these pathways to SUD and these comorbid phenotypes in order to better understand the genetic relationship between them. Methods: We tested the association of 275 dopaminergic genes and 176 serotonergic genes with these phenotypes by performing gene-based, gene-set and transcriptome-wide association studies in 11 genome-wide association studies (GWAS) datasets on SUD and related behaviors. Results: At the gene-wide level, 68 DA and 27 5-HT genes were found to be associated with at least one GWAS on SUD or related behavior. Among them, six genes had a pleiotropic effect, being associated with at least three phenotypes: ADH1C, ARNTL, CHRNA3, HPRT1, HTR1B and DRD2. Additionally, we found nominal associations between the DA gene sets and SUD, opioid use disorder, antisocial behavior, irritability and neuroticism, and between the 5-HT-core gene set and neuroticism. Predicted gene expression correlates in brain were also found for 19 DA or 5-HT genes. Discussion: Our study shows a pleiotropic contribution of dopaminergic and serotonergic genes to addiction and related behaviors such as anxiety, irritability, neuroticism and risk-taking behavior, highlighting a role for DA genes, which could explain, in part, the co-occurrence of these phenotype

Differential expression of miR-1249-3p and miR-34b-5p between vulnerable and resilient phenotypes of cocaine addiction

Cocaine addiction is a complex brain disorder involving long-term alterations that leadto loss of control over drug seeking. The transition from recreational use to pathologi-cal consumption is different in each individual, depending on the interaction betweenenvironmental and genetic factors. Epigenetic mechanisms are ideal candidates tostudy psychiatric disorders triggered by these interactions, maintaining persistentmalfunctions in specific brain regions. Here we aim to study brain-region-specific epi-genetic signatures following exposure to cocaine in a mouse model of addiction tothis drug. Extreme subpopulations of vulnerable and resilient phenotypes wereselected to identify miRNA signatures for differential vulnerability to cocaine addic-tion. We used an operant model of intravenous cocaine self-administration to evalu-ate addictive-like behaviour in rodents based on the Diagnostic and StatisticalManual of Mental Disorders Fifth Edition criteria to diagnose substance use disor-ders. After cocaine self-administration, we performed miRNA profiling to comparetwo extreme subpopulations of mice classified as resilient and vulnerable to cocaineaddiction. We found that mmu-miR-34b-5p was downregulated in the nucleusaccumbens of vulnerable mice with high motivation for cocaine. On the other hand,mmu-miR-1249-3p was downregulated on vulnerable mice with high levels of motordisinhibition. The elucidation of the epigenetic profile related to vulnerability to cocaine addiction is expected to help find novel biomarkers that could facilitate theinterventions to battle this devastating disorder

An integrated analysis of genes and functional pathways for aggression in human and rodent models

Human genome-wide association studies (GWAS), transcriptome analyses of animal models, and candidate gene studies have advanced our understanding of the genetic architecture of aggressive behaviors. However, each of these methods presents unique limitations. To generate a more confident and comprehensive view of the complex genetics underlying aggression, we undertook an integrated, cross-species approach. We focused on human and rodent models to derive eight gene lists from three main categories of genetic evidence: two sets of genes identified in GWAS studies, four sets implicated by transcriptome-wide studies of rodent models, and two sets of genes with causal evidence from online Mendelian inheritance in man (OMIM) and knockout (KO) mice reports. These gene sets were evaluated for overlap and pathway enrichment to extract their similarities and differences. We identified enriched common pathways such as the G-protein coupled receptor (GPCR) signaling pathway, axon guidance, reelin signaling in neurons, and ERK/MAPK signaling. Also, individual genes were ranked based on their cumulative weights to quantify their importance as risk factors for aggressive behavior, which resulted in 40 top-ranked and highly interconnected genes. The results of our cross-species and integrated approach provide insights into the genetic etiology of aggression

Cross-disorder genetic analyses implicate dopaminergic signaling as a biological link between Attention-Deficit/Hyperactivity Disorder and obesity measures

Attention-Deficit/Hyperactivity Disorder (ADHD) and obesity are frequently comorbid, genetically correlated, and share brain substrates. The biological mechanisms driving this association are unclear, but candidate systems, like dopaminergic neurotransmission and circadian rhythm, have been suggested. Our aim was to identify the biological mechanisms underpinning the genetic link between ADHD and obesity measures and investigate associations of overlapping genes with brain volumes. We tested the association of dopaminergic and circadian rhythm gene sets with ADHD, body mass index (BMI), and obesity (using GWAS data of N = 53,293, N = 681,275, and N = 98,697, respectively). We then conducted genome-wide ADHD-BMI and ADHD-obesity gene-based meta-analyses, followed by pathway enrichment analyses. Finally, we tested the association of ADHD-BMI overlapping genes with brain volumes (primary GWAS data N = 10,720-10,928; replication data N = 9428). The dopaminergic gene set was associated with both ADHD (P = 5.81 × 10−3) and BMI (P = 1.63 × 10−5); the circadian rhythm was associated with BMI (P = 1.28 × 10−3). The genome-wide approach also implicated the dopaminergic system, as the Dopamine-DARPP32 Feedback in cAMP Signaling pathway was enriched in both ADHD-BMI and ADHD-obesity results. The ADHD-BMI overlapping genes were associated with putamen volume (P = 7.7 × 10−3; replication data P = 3.9 × 10−2) a brain region with volumetric reductions in ADHD and BMI and linked to inhibitory control. Our findings suggest that dopaminergic neurotransmission, partially through DARPP-32-dependent signaling and involving the putamen, is a key player underlying the genetic overlap between ADHD and obesity measures. Uncovering shared etiological factors underlying the frequently observed ADHD-obesity comorbidity may have important implications in terms of prevention and/or efficient treatment of these conditions

Exploring the contribution to ADHD of genes involved in Mendelian Disorders presenting with hyperactivity and/or inattention

Attention-deficit hyperactivity disorder (ADHD) is a complex neurodevelopmental disorder characterized by hyperactivity, impulsivity, and/or inattention, which are symptoms also observed in many rare genetic disorders. We searched for genes involved in Mendelian disorders presenting with ADHD symptoms in the Online Mendelian Inheritance in Man (OMIM) database, to curate a list of new candidate risk genes for ADHD.We explored the enrichment of functions and pathways in this gene list, and tested whether rare or common variants in these genes are associated with ADHD or with its comorbidities. We identified 139 genes, causal for 137 rare disorders, mainly related to neurodevelopmental and brain function. Most of these Mendelian disorders also present with other psychiatric traits that are often comorbid with ADHD. Using whole exome sequencing (WES) data from 668 ADHD cases, we found rare variants associated with the dimension of the severity of inattention symptoms in three genes: KIF11, WAC, and CRBN. Then, we focused on common variants and identified six genes associated with ADHD (in 19,099 cases and 34,194 controls): MANBA, UQCC2, HIVEP2, FOPX1, KANSL1, and AUH. Furthermore, HIVEP2, FOXP1, and KANSL1 were nominally associated with autism spectrum disorder (ASD) (18,382 cases and 27,969 controls), as well as HIVEP2 with anxiety (7016 cases and 14,475 controls), and FOXP1 with aggression (18,988 individuals), which is in line with the symptomatology of the rare disorders they are responsible for. In conclusion, inspecting Mendelian disorders and the genes responsible for them constitutes a valuable approach for identifying new risk genes and the mechanisms of complex disorders

Variants of the aggression-related RBFOX1 gene in a population representative birth cohort study: aggressiveness, personality and alcohol use disorder

Background: Recently, RBFOX1, a gene encoding an RNA binding protein, has consistently been associated with aggressive and antisocial behavior. Several loci in the gene have been nominally associated with aggression in genome-wide association studies, the risk alleles being more frequent in the general population. We have hence examined the association of four RBFOX1 single nucleotide polymorphisms, previously found related to aggressive traits, with aggressiveness, personality, and alcohol use disorder in birth cohort representative samples. Methods: We used both birth cohorts of the Estonian Children Personality Behavior and Health Study (ECPBHS; original n = 1,238). Aggressiveness was assessed using the Buss-Perry Aggression Questionnaire and the Lifetime History of Aggressiveness structured interview at age 25 (younger cohort) or 33 (older cohort). Big Five personality at age 25 was measured with self-reports and the lifetime occurrence of alcohol use disorder assessed with the MINI interview. RBFOX1 polymorphisms rs809682, rs8062784, rs12921846, and rs6500744 were genotyped in all participants. Given the restricted size of the sample, correction for multiple comparisons was not applied. Results: Aggressiveness was not significantly associated with the RBFOX1 genotype. RBFOX1 rs8062784 was associated with neuroticism and rs809682 with extraversion. Two out of four analyzed RBFOX1 variants, rs8062784 and rs12921846, were associated with the occurrence of alcohol use disorder. Conclusions: In the birth cohort representative sample of the ECPBHS, no association of RBFOX1 with aggressiveness was found, but RBFOX1 variants affected basic personality traits and the prevalence of alcohol use disorder. Future studies on RBFOX1 should consider the moderating role of personality and alcohol use patterns in aggressiveness

RBFOX1, encoding a splicing regulator, is a candidate gene for aggressive behavior

The RBFOX1 gene (or A2BP1) encodes a splicing factor important for neuronal development that has been related to autism spectrum disorder and other neurodevelopmental phenotypes. Evidence from complementary sources suggests that this gene contributes to aggressive behavior. Suggestive associations with RBFOX1 have been identified in genome-wide association studies (GWAS) of anger, conduct disorder, and aggressive behavior. Nominal association signals in RBFOX1 were also found in an epigenome-wide association study (EWAS) of aggressive behavior. Also, variants in this gene affect temporal lobe volume, a brain area that is altered in several aggression-related phenotypes. In animals, this gene has been shown to modulate aggressive behavior in Drosophila. RBFOX1 has also been associated with canine aggression and is upregulated in mice that show increased aggression after frustration of an expected reward. Associated common genetic variants as well as rare duplications and deletions affecting RBFOX1 have been identified in several psychiatric and neurodevelopmental disorders that are often comorbid with aggressive behaviors. In this paper, we comprehensively review the cumulative evidence linking RBFOX1 to aggression behavior and provide new results implicating RBFOX1 in this phenotype. Most of these studies (genetic and epigenetic analyses in humans, neuroimaging genetics, gene expression and animal models) are hypothesis-free, which strengthens the validity of the findings, although all the evidence is nominal and should therefore be taken with caution. Further studies are required to clarify in detail the role of this gene in this complex phenotype

Exploring allele specific methylation in drug dependence susceptibility

Drug dependence is a neuropsychiatric condition that involves genetic, epigenetic and environmental factors. Allele-specific methylation (ASM) is a common and stable epigenetic mechanism that involves genetic variants correlating with differential levels of methylation at CpG sites. We selected 182 single-nucleotide polymorphisms (SNPs) described to influence cis ASM in human brain regions to evaluate their possible contribution to drug dependence susceptibility. We performed a case-control association study in a discovery sample of 578 drug-dependent patients (including 428 cocaine-dependent subjects) and 656 controls from Spain, and then, we followed-up the significant associations in an independent sample of 1,119 cases (including 589 cocaine-dependent subjects) and 1,092 controls. In the discovery sample, we identified five nominal associations, one of them replicated in the follow-up sample (rs6020251). The pooled analysis revealed an association between drug dependence and rs6020251 but also rs11585570, both overcoming the Bonferroni correction for multiple testing. We performed the same analysis considering only cocaine-dependent patients and obtained similar results. The rs6020251 variant correlates with differential methylation levels of cg17974185 and lies in the first intron of the CTNNBL1 gene, in a genomic region with multiple histone marks related to enhancer and promoter regions in brain. Rs11585570 is an eQTL in brain and blood for the SCP2 and ECHDC2 genes and correlates with differential methylation of cg27535305 and cg13461509, located in the promoter regions of both genes. To conclude, using an approach that combines genetic and epigenetic data, we highlighted the CTNNBL1, SCP2 and ECHDC2 genes as potential contributors to drug dependence susceptibility