Genomewide association study of cocaine dependence and related traits
Description:
The discovery GWAS dataset (n=5,697 subjects) was genotyped using the Illumina OmniQuad microarray (890,000 analyzed SNPs). Additional genotypes were imputed based on the 1000 Genomes reference panel. Top-ranked findings were evaluated by incorporating information from publicly available GWAS data from 4,063 subjects. Then, the most significant GWAS SNPs were genotyped in 2,549 independent subjects. One genomewide-significant (GWS) result: rs7086629 at the FAM53B (“family with sequence similarity 53, member B”) locus. This was supported in both AAs and EAs; p-value (meta-analysis of all samples) =4.28×10−8. The gene maps to the same chromosomal region as the maximum peak we observed in a previous linkage study. NCOR2 (nuclear receptor corepressor 1) SNP rs150954431 was associated with p=1.19×10−9 in the EA discovery sample. SNP rs2456778, which maps to CDK1 (“cyclin-dependent kinase 1”), was associated with cocaine-induced paranoia in AAs in the discovery sample only (p=4.68×10−8).
Authors:
J Gelernter, R Sherva, R Koesterer, L Almasy, H Zhao, HR Kranzler, L Farrer
List of positional candidate genes after correcting for multiple testing and controlling the false discovery rate from genome wide association studies (GWAS) retrieved from the NHGRI-EBI Catalog of published genome-wide association studies (http://www.ebi.ac.uk/gwas/). The disease/trait examined in this study, as reported by the authors, was Cocaine dependence. The EFO term cocaine dependence was annotated to this set after curation by NHGRI-EBI. Intergenic SNPS were mapped to both the upstream and downstream gene. P-value uploaded. This gene set was generated using gwas2gs v. 0.1.8 and the GWAS Catalog v. 1.0.1.
Authors:
J Gelernter, R Sherva, R Koesterer, L Almasy, H Zhao, HR Kranzler, L Farrer
Genes exhibiting midbrain expression levels that map genome-wide significantly to cocaine dose-response curve in mice
Description:
Mice from 39 BXD strains acquired cocaine IVSA (0.56 mg/kg/infusion). Mice from 29 BXD strains completed a full dose-response curve (0.032-1.8 mg/kg/infusion). We identified independent genetic correlations between cocaine IVSA and measures of environmental exploration and cocaine sensitization. We identified genome-wide significant quantitative trait loci (QTL) on chromosomes 7 and 11 associated with shifts in the dose-response curve and on chromosome 16 associated with sessions to acquire cocaine IVSA. Using publicly available gene expression data from the nucleus accumbens, midbrain, and prefrontal cortex of drug-naïve mice, we identified Aplp1 and Cyfip2 as positional candidates underlying the behavioral QTL on chromosomes 7 and 11, respectively. A genome-wide significant trans-eQTL linking Fam53b (a GWAS candidate for human cocaine dependence) on chromosome 7 to the cocaine IVSA behavioral QTL on chromosome 11 was identified in the midbrain; Fam53b and Cyfip2 were co-expressed genome-wide significantly in the midbrain. Shown here are Genes exhibiting midbrain expression levels that map genome-wide significantly to the behavioral QTL associated with infusions on the ascending limb of the dose-response curve (IAL phenotype).
Authors:
Price E Dickson, Mellessa M Miller, Michele A Calton, Jason A Bubier, Melloni N Cook, Daniel Goldowitz, Elissa J Chesler, Guy Mittleman
The total transcriptome including genes that are differentially expressed in cocaine addicts compared to control subjects. Post-mortem brain samples were collected from the dorsolateral prefrontal cortex (dlPFC) of the cocaine addict group and the control group. To assess gene expression, RNA-seq was performed. Data taken from Supplementary Table 2. Values presented are k.diff values. Data available from GEO with accession number GSE99349."
Authors:
Efrain A Ribeiro, Joseph R Scarpa, Susanna P Garamszegi, Andrew Kasarskis, Deborah C Mash, Eric J Nestler
Data from GEO GSE194368 and analyzed using GEO2R, only top gene shown. Authors identified transcriptional adaptations of GR signaling in the amygdala of humans with OUD. Thus, GRs, their coregulators and downstream systems may represent viable therapeutic targets to treat the “stress side” of OUD.
Authors:
Stephanie A Carmack, Janaina C M Vendruscolo, M Adrienne McGinn, Jorge Miranda-Barrientos, Vez Repunte-Canonigo, Gabriel D Bosse, Daniele Mercatelli, Federico M Giorgi, Yu Fu, Anthony J Hinrich, Francine M Jodelka, Karen Ling, Robert O Messing, Randall T Peterson, Frank Rigo, Scott Edwards, Pietro P Sanna, Marisela Morales, Michelle L Hastings, George F Koob, Leandro F Vendruscolo
Transcriptional alterations in dorsolateral prefrontal cortex and nucleus accumbens implicate neuroinflammation and synaptic remodeling in opioid use disorder. Transcriptomic profile of 20 control subjects and 20 OUD subjects in brain region DLPFC and NAC. Analyzed using GEO2R (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE174409) separately for each brain region, comparing OUD and control samples.
Authors:
Xiangning Xue, Wei Zong, Jill R Glausier, Sam-Moon Kim, Micah A Shelton, BaDoi N Phan, Chaitanya Srinivasan, Andreas R Pfenning, George C Tseng, David A Lewis, Marianne L Seney, Ryan W Logan
Positional candidates for QTL associated with cocaine IVSA phenotypes_cis-eQTL p-value
Description:
Mice from 39 BXD strains acquired cocaine IVSA (0.56 mg/kg/infusion). Mice from 29 BXD strains completed a full dose-response curve (0.032-1.8 mg/kg/infusion). We identified independent genetic correlations between cocaine IVSA and measures of environmental exploration and cocaine sensitization. We identified genome-wide significant quantitative trait loci (QTL) on chromosomes 7 and 11 associated with shifts in the dose-response curve and on chromosome 16 associated with sessions to acquire cocaine IVSA. Using publicly available gene expression data from the nucleus accumbens, midbrain, and prefrontal cortex of drug-naïve mice, we identified Aplp1 and Cyfip2 as positional candidates underlying the behavioral QTL on chromosomes 7 and 11, respectively. A genome-wide significant trans-eQTL linking Fam53b (a GWAS candidate for human cocaine dependence) on chromosome 7 to the cocaine IVSA behavioral QTL on chromosome 11 was identified in the midbrain; Fam53b and Cyfip2 were co-expressed genome-wide significantly in the midbrain. Shown here are Priority positional candidates for behavioral QTL associated with cocaine IVSA phenotypes_cis-eQTL p-value
Authors:
Price E Dickson, Mellessa M Miller, Michele A Calton, Jason A Bubier, Melloni N Cook, Daniel Goldowitz, Elissa J Chesler, Guy Mittleman
Genome-wide significant covariation of IVSA phenotypes and gene expression in the mouse brain
Description:
Mice from 39 BXD strains acquired cocaine IVSA (0.56 mg/kg/infusion). Mice from 29 BXD strains completed a full dose-response curve (0.032-1.8 mg/kg/infusion). We identified independent genetic correlations between cocaine IVSA and measures of environmental exploration and cocaine sensitization. We identified genome-wide significant quantitative trait loci (QTL) on chromosomes 7 and 11 associated with shifts in the dose-response curve and on chromosome 16 associated with sessions to acquire cocaine IVSA. Using publicly available gene expression data from the nucleus accumbens, midbrain, and prefrontal cortex of drug-naïve mice, we identified Aplp1 and Cyfip2 as positional candidates underlying the behavioral QTL on chromosomes 7 and 11, respectively. A genome-wide significant trans-eQTL linking Fam53b (a GWAS candidate for human cocaine dependence) on chromosome 7 to the cocaine IVSA behavioral QTL on chromosome 11 was identified in the midbrain; Fam53b and Cyfip2 were co-expressed genome-wide significantly in the midbrain. Shown here are Genome-wide significant covariation of IVSA phenotypes and gene expression in the nucleus accumbens (NAc), midbrain, and prefrontal cortex (PFC).
Authors:
Price E Dickson, Mellessa M Miller, Michele A Calton, Jason A Bubier, Melloni N Cook, Daniel Goldowitz, Elissa J Chesler, Guy Mittleman
Cerebellum Gene Expression Correlates for TAILCLIP_LAT_SEC measured in BXD RI Females & Males obtained using SJUT Cerebellum mRNA M430 (Mar05) RMA. The TAILCLIP_LAT_SEC measures Mechanical Nociception - Tail Clip Test under the domain Pain. The correlates were thresholded at a p-value of less than 0.001.
Authors:
Philip VM, Duvvuru S, Gomero B, Ansah TA, Blaha CD, Cook MN, Hamre KM, Lariviere WR, Matthews DB, Mittleman G, Goldowitz D, Chesler EJ
QTL for nicotine sensitivity on Chr7 at D7Mit66 (116.91 Mbp , Build 37)
Description:
nicotine sensitivity spans 91.91 - 141.91 Mbp (NCBI Build 37) on Chr7. This interval was obtained by using an interval width of 25 Mbp around the peak marker (Build 37, MGI, http://informatics.jax.org).
QTL for alcohol preference locus on Chr7 at D7Mit105 (126.73 Mbp , Build 37)
Description:
alcohol preference locus spans 101.73 - 151.73 Mbp (NCBI Build 37) on Chr7. This interval was obtained by using an interval width of 25 Mbp around the peak marker (Build 37, MGI, http://informatics.jax.org).
alcohol preference locus 14, female specific at D7Mit105 with a LOD score of 1.84 (p < 0.004) spans and preference correlation of 0.591 101.73 - 151.73 Mbp (NCBI Build 37) on Chr7. This interval was obtained by using an interval width of 25 Mbp around the peak marker (Build 37, MGI, http://informatics.jax.org).
QTL for METH responses for home cage activity on Chr7 at D7Mit12 (129.57 Mbp , Build 37)
Description:
METH responses for home cage activity spans 104.57 - 154.57 Mbp (NCBI Build 37) on Chr7. This interval was obtained by using an interval width of 25 Mbp around the peak marker (Build 37, MGI, http://informatics.jax.org).
QTL for METH responses for body temperature on Chr7 at Xmv76 (137.02 Mbp , Build 37)
Description:
METH responses for body temperature spans 112.02 - 162.02 Mbp (NCBI Build 37) on Chr7. This interval was obtained by using an interval width of 25 Mbp around the peak marker (Build 37, MGI, http://informatics.jax.org).
Genes associated with Homo sapiens that interact with the MeSH term 'Arsenic' (D001151). Incorporates data from 87 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Genes associated with Oryctolagus cuniculus that interact with the MeSH term 'Ionomycin' (D015759). Incorporates data from 6 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Genes associated with Homo sapiens that interact with the MeSH term 'Vitamin E' (D014810). Incorporates data from 6 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Genes associated with Homo sapiens that interact with the MeSH term 'Cyclosporine' (D016572). Incorporates data from 1 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Genes associated with Homo sapiens that interact with the MeSH term 'Selenium' (D012643). Incorporates data from 99 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Genes associated with Homo sapiens that interact with the MeSH term 'Valproic Acid' (D014635). Incorporates data from 1238 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Genes associated with Oryctolagus cuniculus that interact with the MeSH term 'Tetradecanoylphorbol Acetate' (D013755). Incorporates data from 2 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Genes associated with Homo sapiens that interact with the MeSH term 'Benzo(a)pyrene' (D001564). Incorporates data from 3 publications curated by the Comparative Toxicogenomics Database (CTD). ODE Gene scores represent number of supporting publications per gene.
Ethanol Induced Hypothermia Chr# 7 rs13479153(25722935) with right flanking marker rs3700068(4187548) and left marker rs3716088(140189839). This was mapped in 300 + (b6x129)F2 mice.
Average rotarod training latency Chr# 7 mCV23423763(68111945) with right flanking marker rs3700068(4187548) and left marker rs3663988(146505067). This was mapped in 300 + (b6x129)F2 mice.
QTL associated with "alcohol preference locus 12, male specific". This interval was obtained by using a fixed interval width of 25 Mbp around the peak marker (135707912)
GeneWeaver