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

New Distance-Based approach for Genome-Wide Association Studies

With the raise of genome-wide association studies (GWAS), the analysis of typical GWAS data sets with thousands of potentially predictive single nucleotide-polymorphisms (SNPs) has become crucial in Biomedicine research. Here, we propose a new method to identify SNPs related to disease in case-control studies. The method, based on genetic distances between individuals, takes into account the possible population substructure, and avoids the issues of multiple testing. The method provides two ordered lists of SNPs; one with SNPs which minor alleles can be considered risk alleles for the disease, and another one with SNPs which minor alleles can be considered as protective. These two lists provide a useful tool to help the researcher to decide where to focus attention in a first stage

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

Chiari malformation type I: a case-control association study of 58 developmental genes

Chiari malformation type I (CMI) is a disorder characterized by hindbrain overcrowding into an underdeveloped posterior cranial fossa (PCF), often causing progressive neurological symptoms. The etiology of CMI remains unclear and is most likely multifactorial. A putative genetic contribution to CMI is suggested by familial aggregation and twin studies. Experimental models and human morphometric studies have suggested an underlying paraxial mesoderm insufficiency. We performed a case-control association study of 303 tag single nucleotide polymorphisms (SNP) across 58 candidate genes involved in early paraxial mesoderm development in a sample of 415 CMI patients and 524 sex-matched controls. A subgroup of patients diagnosed with classical, small-PCF CMI by means of MRI-based PCF morphometry (n = 186), underwent additional analysis. The genes selected are involved in signalling gradients occurring during segmental patterning of the occipital somites (FGF8, Wnt, and retinoic acid pathways and from bone morphogenetic proteins or BMP, Notch, Cdx and Hox pathways) or in placental angiogenesis, sclerotome development or CMI-associated syndromes. Single-marker analysis identified nominal associations with 18 SNPs in 14 genes (CDX1, FLT1, RARG, NKD2, MSGN1, RBPJ1, FGFR1, RDH10, NOG, RARA, LFNG, KDR, ALDH1A2, BMPR1A) considering the whole CMI sample. None of these overcame corrections for multiple comparisons, in contrast with four SNPs in CDX1, FLT1 and ALDH1A2 in the classical CMI group. Multiple marker analysis identified a risk haplotype for classical CMI in ALDH1A2 and CDX1. Furthermore, we analyzed the possible contributions of the most significantly associated SNPs to different PCF morphometric traits. These findings suggest that common variants in genes involved in somitogenesis and fetal vascular development may confer susceptibility to CMI

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

A causal effects of gut microbiota in the development of migraine

Background: The causal association between the gut microbiome and the development of migraine and its subtypes remains unclear. Methods: The single nucleotide polymorphisms concerning gut microbiome were retrieved from the gene-wide association study (GWAS) of the MiBioGen consortium. The summary statistics datasets of migraine, migraine with aura (MA), and migraine without aura (MO) were obtained from the GWAS meta-analysis of the International Headache Genetics Consortium (IHGC) and FinnGen consortium. Inverse variance weighting (IVW) was used as the primary method, complemented by sensitivity analyses for pleiotropy and increasing robustness. Results: In IHGC datasets, ten, five, and nine bacterial taxa were found to have a causal association with migraine, MA, and MO, respectively, (IVW, all P < 0.05). Genus.Coprococcus3 and genus.Anaerotruncus were validated in FinnGen datasets. Nine, twelve, and seven bacterial entities were identified for migraine, MA, and MO, respectively. The causal association still exists in family.Bifidobacteriaceae and order.Bifidobacteriales for migraine and MO after FDR correction. The heterogeneity and pleiotropy analyses confirmed the robustness of IVW results. Conclusion: Our study demonstrates that gut microbiomes may exert causal effects on migraine, MA, and MO. We provide novel evidence for the dysfunction of the gut-brain axis on migraine. Future study is required to verify the relationship between gut microbiome and the risk of migraine and its subtypes and illustrate the underlying mechanism between them

Two-stage case-control association study of dopamine-related genes and migraine

Background We previously reported risk haplotypes for two genes related with serotonin and dopamine metabolism: MAOA in migraine without aura and DDC in migraine with aura. Herein we investigate the contribution to migraine susceptibility of eight additional genes involved in dopamine neurotransmission. Methods We performed a two-stage case-control association study of 50 tag single nucleotide polymorphisms (SNPs), selected according to genetic coverage parameters. The first analysis consisted of 263 patients and 274 controls and the replication study was composed by 259 cases and 287 controls. All cases were diagnosed according to ICHD-II criteria, were Spanish Caucasian, and were sex-matched with control subjects. Results Single-marker analysis of the first population identified nominal associations of five genes with migraine. After applying a false discovery rate correction of 10%, the differences remained significant only for DRD2 (rs2283265) and TH (rs2070762). Multiple-marker analysis identified a five-marker T-C-G-C-G (rs12363125-rs2283265-rs2242592-rs1554929-rs2234689) risk haplotype in DRD2 and a two-marker A-C (rs6356-rs2070762) risk haplotype in TH that remained significant after correction by permutations. These results, however, were not replicated in the second independent cohort. Conclusion The present study does not support the involvement of the DRD1, DRD2, DRD3, DRD5, DBH, COMT, SLC6A3 and TH genes in the genetic predisposition to migraine in the Spanish population