Genetic Architecture and Molecular Neuropathology of Human Cocaine Addiction

Validation of Enhancer Regions in Primary Human Neural Progenitor Cells using Capture STARR-seq

Genome-wide association studies (GWAS) and expression analyses implicate noncoding regulatory regions as harboring risk factors for psychiatric disease, but functional characterization of these regions remains limited. We performed capture STARR-sequencing of over 78,000 candidate regions to identify active enhancers in primary human neural progenitor cells (phNPCs). We selected candidate regions by integrating data from NPCs, prefrontal cortex, developmental timepoints, and GWAS. Over 8,000 regions demonstrated enhancer activity in the phNPCs, and we linked these regions to over 2,200 predicted target genes. These genes are involved in neuronal and psychiatric disease-associated pathways, including dopaminergic synapse, axon guidance, and schizophrenia. We functionally validated a subset of these enhancers using mutation STARR-sequencing and CRISPR deletions, demonstrating the effects of genetic variation on enhancer activity and enhancer deletion on gene expression. Overall, we identified thousands of highly active enhancers and functionally validated a subset of these enhancers, improving our understanding of regulatory networks underlying brain function and disease.

Regulation of hippocampal mossy fiber-CA3 synapse function by a Bcl11b/C1ql2/Nrxn3(25b+) pathway

The transcription factor Bcl11b has been linked to neurodevelopmental and neuropsychiatric disorders associated with synaptic dysfunction. Bcl11b is highly expressed in dentate gyrus granule neurons and is required for the structural and functional integrity of mossy fiber-CA3 synapses. The underlying molecular mechanisms, however, remained unclear. We show in mice that the synaptic organizer molecule C1ql2 is a direct functional target of Bcl11b that regulates synaptic vesicle recruitment and long-term potentiation at mossy fiber-CA3 synapses in vivo and in vitro. Furthermore, we demonstrate C1ql2 to exert its functions through direct interaction with a specific splice variant of neurexin-3, Nrxn3(25b+). Interruption of C1ql2-Nrxn3(25b+) interaction by expression of a non-binding C1ql2 mutant or by deletion of Nrxn3 in the dentate gyrus granule neurons recapitulates major parts of the Bcl11b as well as C1ql2 mutant phenotype. Together, this study identifies a novel C1ql2-Nrxn3(25b+)-dependent signaling pathway through which Bcl11b controls mossy fiber-CA3 synapse function. Thus, our findings contribute to the mechanistic understanding of neurodevelopmental disorders accompanied by synaptic dysfunction.

USP7/Maged1-mediated H2A monoubiquitination in the paraventricular thalamus: an epigenetic mechanism involved in cocaine use disorder

The risk of developing drug addiction is strongly influenced by the epigenetic landscape and chromatin remodeling. While histone modifications such as methylation and acetylation have been studied in the ventral tegmental area and nucleus accumbens (NAc), the role of H2A monoubiquitination remains unknown. Our investigations, initially focused on the scaffold protein melanoma-associated antigen D1 (Maged1), reveal that H2A monoubiquitination in the paraventricular thalamus (PVT) significantly contributes to cocaine-adaptive behaviors and transcriptional repression induced by cocaine. Chronic cocaine use increases H2A monoubiquitination, regulated by Maged1 and its partner USP7. Accordingly, Maged1 specific inactivation in thalamic Vglut2 neurons, or USP7 inhibition, blocks cocaine-evoked H2A monoubiquitination and cocaine locomotor sensitization. Additionally, genetic variations in MAGED1 and USP7 are linked to altered susceptibility to cocaine addiction and cocaine-associated symptoms in humans. These findings unveil an epigenetic modification in a non-canonical reward pathway of the brain and a potent marker of epigenetic risk factors for drug addiction in humans.

Loss of CTRP10 results in female obesity with preserved metabolic health

Obesity is a major risk factor for type 2 diabetes, dyslipidemia, cardiovascular disease, and hypertension. Intriguingly, there is a subset of metabolically healthy obese (MHO) individuals who are seemingly able to maintain a healthy metabolic profile free of metabolic syndrome. The molecular underpinnings of MHO, however, are not well understood. Here, we report that CTRP10/C1QL2-deficient mice represent a unique female model of MHO. CTRP10 modulates weight gain in a striking and sexually dimorphic manner. Female, but not male, mice lacking CTRP10 develop obesity with age on a low-fat diet while maintaining an otherwise healthy metabolic profile. When fed an obesogenic diet, female Ctrp10 knockout (KO) mice show rapid weight gain. Despite pronounced obesity, Ctrp10 KO female mice do not develop steatosis, dyslipidemia, glucose intolerance, insulin resistance, oxidative stress, or low-grade inflammation. Obesity is largely uncoupled from metabolic dysregulation in female KO mice. Multi-tissue transcriptomic analyses highlighted gene expression changes and pathways associated with insulin-sensitive obesity. Transcriptional correlation of the differentially expressed gene (DEG) orthologous in humans also show sex differences in gene connectivity within and across metabolic tissues, underscoring the conserved sex-dependent function of CTRP10. Collectively, our findings suggest that CTRP10 negatively regulates body weight in females, and that loss of CTRP10 results in benign obesity with largely preserved insulin sensitivity and metabolic health. This female MHO mouse model is valuable for understanding sex-biased mechanisms that uncouple obesity from metabolic dysfunction.

Cis-meQTL for cocaine use-associated DNA methylation in an HIV-positive cohort show pleiotropic effects on multiple traits

BACKGROUND

Cocaine use (CU) is associated with psychiatric and medical diseases. Little is known about the mechanisms of CU-related comorbidities. Findings from preclinical and clinical studies have suggested that CU is associated with aberrant DNA methylation (DNAm) that may be influenced by genetic variants [i.e., methylation quantitative trait loci (meQTLs)]. In this study, we mapped cis-meQTLs for CU-associated DNAm sites (CpGs) in an HIV-positive cohort (Ntotal = 811) and extended the meQTLs to multiple traits.

RESULTS

We conducted cis-meQTL analysis for 224 candidate CpGs selected for their association with CU in blood. We identified 7,101 significant meQTLs [false discovery rate (FDR) < 0.05], which mostly mapped to genes involved in immunological functions and were enriched in immune pathways. We followed up the meQTLs using phenome-wide association study and trait enrichment analyses, which revealed 9 significant traits. We tested for causal effects of CU on these 9 traits using Mendelian Randomization and found evidence that CU plays a causal role in increasing hypertension (p-value = 2.35E-08) and decreasing heel bone mineral density (p-value = 1.92E-19).

CONCLUSIONS

These findings suggest that genetic variants for CU-associated DNAm have pleiotropic effects on other relevant traits and provide new insights into the causal relationships between cocaine use and these complex traits.

Genetic pathways regulating the longitudinal acquisition of cocaine self-administration in a panel of inbred and recombinant inbred mice

To identify addiction genes, we evaluate intravenous self-administration of cocaine or saline in 84 inbred and recombinant inbred mouse strains over 10 days. We integrate the behavior data with brain RNA-seq data from 41 strains. The self-administration of cocaine and that of saline are genetically distinct. We maximize power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci are identified in the genome-wide association study. A transcriptome-wide association study highlights the Trpv2 ion channel as a key locus for cocaine self-administration as well as identifying 17 additional genes, including Arhgef26, Slc18b1, and Slco5a1. We find numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance of Trpv2, an ionotropic cannabinoid receptor, for the response to cocaine.

Gene expression genetics of the striatum of Diversity Outbred mice

Brain transcriptional variation is a heritable trait that mediates complex behaviors, including addiction. Expression quantitative trait locus (eQTL) mapping reveals genomic regions harboring genetic variants that influence transcript abundance. In this study, we profiled transcript abundance in the striatum of 386 Diversity Outbred (J:DO) mice of both sexes using RNA-Seq. All mice were characterized using a behavioral battery of widely-used exploratory and risk-taking assays prior to transcriptional profiling. We performed eQTL mapping, incorporated the results into a browser-based eQTL viewer, and deposited co-expression network members in GeneWeaver. The eQTL viewer allows researchers to query specific genes to obtain allelic effect plots, analyze SNP associations, assess gene expression correlations, and apply mediation analysis to evaluate whether the regulatory variant is acting through the expression of another gene. GeneWeaver allows multi-species comparison of gene sets using statistical and combinatorial tools. This data resource allows users to find genetic variants that regulate differentially expressed transcripts and place them in the context of other studies of striatal gene expression and function in addiction-related behavior.

Gene expression genetics of the striatum of Diversity Outbred mice

Brain transcriptional variation is a heritable trait that mediates complex behaviors, including addiction. Expression quantitative trait locus (eQTL) mapping reveals genomic regions harboring genetic variants that influence transcript abundance. In this study, we profiled transcript abundance in the striatum of 386 Diversity Outbred (J:DO) mice of both sexes using RNA-Seq. All mice were characterized using a behavioral battery of widely-used exploratory and risk-taking assays prior to transcriptional profiling. We performed eQTL mapping, incorporated the results into a browser-based eQTL viewer, and deposited co-expression network members in GeneWeaver. The eQTL viewer allows researchers to query specific genes to obtain allelic effect plots, analyze SNP associations, assess gene expression correlations, and apply mediation analysis to evaluate whether the regulatory variant is acting through the expression of another gene. GeneWeaver allows multi-species comparison of gene sets using statistical and combinatorial tools. This data resource allows users to find genetic variants that regulate differentially expressed transcripts and place them in the context of other studies of striatal gene expression and function in addiction-related behavior.

Cardiolipin and OPA1 Team up for Methamphetamine-Induced Locomotor Activity by Promoting Neuronal Mitochondrial Fusion in the Nucleus Accumbens of Mice

Mitochondria are highly dynamic organelles with coordinated cycles of fission and fusion occurring continuously to satisfy the energy demands in the complex architecture of neurons. How mitochondria contribute to addicted drug-induced adaptable mitochondrial networks and neuroplasticity remains largely unknown. Through liquid chromatography-mass spectrometry-based lipidomics, we first analyzed the alteration of the mitochondrial lipidome of three mouse brain areas in methamphetamine (METH)-induced locomotor activity and conditioned place preference. The results showed that METH remodeled the mitochondrial lipidome of the hippocampus, nucleus accumbens (NAc), and striatum in both models. Notably, mitochondrial hallmark lipid cardiolipin (CL) was specifically increased in the NAc in METH-induced hyperlocomotor activity, which was accompanied by an elongated giant mitochondrial morphology. Moreover, METH significantly boosted mitochondrial respiration and ATP generation as well as the copy number of mitochondrial genome DNA in the NAc. By screening the expressions of mitochondrial dynamin-related proteins, we found that repeated METH significantly upregulated the expression of long-form optic atrophy type 1 (L-OPA1) and enhanced the interaction of L-OPA1 with CL, which may promote mitochondrial fusion in the NAc. On the contrary, neuronal OPA1 depletion in the NAc not only recovered the dysregulated mitochondrial morphology and synaptic vesicle distribution induced by METH but also attenuated the psychomotor effect of METH. Collectively, upregulated CL and OPA1 cooperate to mediate METH-induced adaptation of neuronal mitochondrial dynamics in the NAc, which correlates with the psychomotor effect of METH. These findings propose a potential therapeutic approach for METH addiction by inhibiting neuronal mitochondrial fusion.

Genome- and transcriptome-wide splicing associations with alcohol use disorder

Genetic mechanisms of alternative mRNA splicing have been shown in the brain for a variety of neuropsychiatric traits, but not substance use disorders. Our study utilized RNA-sequencing data on alcohol use disorder (AUD) in four brain regions (n = 56; ages 40-73; 100% 'Caucasian'; PFC, NAc, BLA and CEA) and genome-wide association data on AUD (n = 435,563, ages 22-90; 100% European-American). Polygenic scores of AUD were associated with AUD-related alternative mRNA splicing in the brain. We identified 714 differentially spliced genes between AUD vs controls, which included both putative addiction genes and novel gene targets. We found 6463 splicing quantitative trait loci (sQTLs) that linked to the AUD differentially spliced genes. sQTLs were enriched in loose chromatin genomic regions and downstream gene targets. Additionally, the heritability of AUD was enriched for DNA variants in and around differentially spliced genes associated with AUD. Our study also performed splicing transcriptome-wide association studies (TWASs) of AUD and other drug use traits that unveiled specific genes for follow-up and splicing correlations across SUDs. Finally, we showed that differentially spliced genes between AUD vs control were also associated with primate models of chronic alcohol consumption in similar brain regions. Our study found substantial genetic contributions of alternative mRNA splicing in AUD.

Interaction between maternal immune activation and peripubertal stress in rats: impact on cocaine addiction-like behaviour, morphofunctional brain parameters and striatal transcriptome

Substance use disorders are more prevalent in schizophrenia, but the causal links between both conditions remain unclear. Maternal immune activation (MIA) is associated with schizophrenia which may be triggered by stressful experiences during adolescence. Therefore, we used a double-hit rat model, combining MIA and peripubertal stress (PUS), to study cocaine addiction and the underlying neurobehavioural alterations. We injected lipopolysaccharide or saline on gestational days 15 and 16 to Sprague-Dawley dams. Their male offspring underwent five episodes of unpredictable stress every other day from postnatal day 28 to 38. When animals reached adulthood, we studied cocaine addiction-like behaviour, impulsivity, Pavlovian and instrumental conditioning, and several aspects of brain structure and function by MRI, PET and RNAseq. MIA facilitated the acquisition of cocaine self-administration and increased the motivation for the drug; however, PUS reduced cocaine intake, an effect that was reversed in MIA + PUS rats. We found concomitant brain alterations: MIA + PUS altered the structure and function of the dorsal striatum, increasing its volume and interfering with glutamatergic dynamics (PUS decreased the levels of NAA + NAAG but only in LPS animals) and modulated specific genes that could account for the restoration of cocaine intake such as the pentraxin family. On its own, PUS reduced hippocampal volume and hyperactivated the dorsal subiculum, also having a profound effect on the dorsal striatal transcriptome. However, these effects were obliterated when PUS occurred in animals with MIA experience. Our results describe an unprecedented interplay between MIA and stress on neurodevelopment and the susceptibility to cocaine addiction.

The Genetically Informed Neurobiology of Addiction (GINA) model

Addictions are heritable and unfold dynamically across the lifespan. One prominent neurobiological theory proposes that substance-induced changes in neural circuitry promote the progression of addiction. Genome-wide association studies have begun to characterize the polygenic architecture undergirding addiction liability and revealed that genetic loci associated with risk can be divided into those associated with a general broad-spectrum liability to addiction and those associated with drug-specific addiction risk. In this Perspective, we integrate these genomic findings with our current understanding of the neurobiology of addiction to propose a new Genetically Informed Neurobiology of Addiction (GINA) model.

Sleep-mediated regulation of reward circuits: implications in substance use disorders

Our modern society suffers from both pervasive sleep loss and substance abuse-what may be the indications for sleep on substance use disorders (SUDs), and could sleep contribute to the individual variations in SUDs? Decades of research in sleep as well as in motivated behaviors have laid the foundation for us to begin to answer these questions. This review is intended to critically summarize the circuit, cellular, and molecular mechanisms by which sleep influences reward function, and to reveal critical challenges for future studies. The review also suggests that improving sleep quality may serve as complementary therapeutics for treating SUDs, and that formulating sleep metrics may be useful for predicting individual susceptibility to SUDs and other reward-associated psychiatric diseases.

A Systematic Review of Genetic Polymorphisms Associated with Bipolar Disorder Comorbid to Substance Abuse

It is currently unknown which genetic polymorphisms are involved in substance use disorder (SUD) comorbid with bipolar disorder (BD). The research on polymorphisms in BD comorbid with SUD (BD + SUD) is summarized in this systematic review. We looked for case-control studies that genetically compared adults and adolescents with BD and SUD, healthy controls, and BD without SUD. PRISMA was used to create our protocol, which is PROSPERO-registered (identification: CRD4221270818). The following bibliographic databases were searched indefinitely until December 2021 to identify potentially relevant articles: PubMed, PsycINFO, Scopus, and Web of Science. This systematic review, after the qualitative analysis of the study selection, included 17 eligible articles. In the selected studies, 66 polymorphisms in 29 genes were investigated. The present work delivers a group of potentially valuable genetic polymorphisms associated with BD + SUD: rs11600996 (ARNTL), rs228642/rs228682/rs2640909 (PER3), PONQ192R (PON1), rs945032 (BDKRB2), rs1131339 (NR4A3), and rs6971 (TSPO). It is important to note that none of those findings have been confirmed by two or more studies; thus, we believe that all the polymorphisms identified in this review require additional evidence to be confirmed.

ΔFOSB: A Potentially Druggable Master Orchestrator of Activity-Dependent Gene Expression

ΔFOSB is a uniquely stable member of the FOS family of immediate early gene AP1 transcription factors. Its accumulation in specific cell types and tissues in response to a range of chronic stimuli is associated with biological phenomena as diverse as memory formation, drug addiction, stress resilience, and immune cell activity. Causal connections between ΔFOSB expression and the physiological and behavioral sequelae of chronic stimuli have been established in rodent and, in some cases, primate models for numerous healthy and pathological states with such preclinical observations often supported by human data demonstrating tissue-specific ΔFOSB expression associated with several specific syndromes. However, the viability of ΔFOSB as a target for therapeutic intervention might be questioned over presumptive concerns of side effects given its expression in such a wide range of cell types and circumstances. Here, we summarize numerous insights from the past three decades of research into ΔFOSB structure, function, mechanisms of induction, and regulation of target genes that support its potential as a druggable target. We pay particular attention to the potential for targeting distinct ΔFOSB isoforms or distinct ΔFOSB-containing multiprotein complexes to achieve cell type or tissue specificity to overcome off-target concerns. We also cover critical gaps in knowledge that currently limit the exploitation of ΔFOSB's therapeutic possibilities and how they may be addressed. Finally, we summarize both current and potential future strategies for generating small molecules or genetic tools for the manipulation of ΔFOSB in the clinic.

Analyses of polymorphisms of intron 2 of OPRK1 (kappa-opioid receptor gene) in association with opioid and cocaine dependence diagnoses in an African-American population

RATIONALE

The dynorphin/kappa-opioid receptor (KOR) system (encoded by PDYN and OPRK1 genes respectively) is highly regulated by repeated exposure to drugs of abuse, including mu-opioid agonists and cocaine. These changes in the dynorphin/KOR system can then influence the rewarding effects of these drugs of abuse. Activation of the dynorphin/KOR system is also thought to have a role in the pro-addictive effects of stress. Recent in vitro assays showed that the OPRK1 intron 2 may function as a genomic enhancer in the regulation KOR expression, and contains a glucocorticiod-responsive sequence site. We hypothesize that SNPs in intron 2 of OPRK1 are associated with categorical opioid or cocaine dependence diagnoses, as well as with dimensional aspects of drug use (i.e., magnitude of drug exposure).

METHODS

This study includes 577 subjects ≥ 18 years old, with African ancestry (AA) from the USA. They were divided into three groups: 152 control subjects, 142 persons with lifetime opioid dependence diagnosis (OD), and 283 subjects with lifetime cocaine dependence diagnosis (CD). Five SNPs (rs16918909, rs7016778, rs997917, rs6473797, rs10111937) that span 10 Kb nucleotides in intron 2 of OPRK1 were used for the association analyses. Genotyping was performed with the Smokescreen® array or sequencing of PCR-amplified DNA fragments. Association analyses for OD and CD diagnoses and the OPRK1 intron 2 alleles were carried out with Fisher's exact test. The Kreek-McHugh-Schluger-Kellogg (KMSK) scales were used for dimensional measure of maximum exposure to specific drugs, using Mann-Whitney tests.

RESULTS

Two SNPs, rs997917 and rs10111937 showed point-wise significant allelic association (p < 0.05) with CD diagnosis, and rs10111937 showed a point-wise significance in association with OD. None of these single SNP associations with categorical diagnoses were significant after correction for multiple testing (p_corr > 0.05). However, significant associations of several genotype patterns (diplotypes) were found with cocaine dependence, but none for opioid dependence. The most significant genotype pattern with cocaine dependence diagnosis occurred for rs6473797 and rs10111937 (p_corr = 0.036, odds ratio = 1.92, FDR < 0.05), and survived correction for multiple testing. Dimensional analyses with KMSK scores show that persons with either rs997917 or rs10111937 variants had greater exposure to cocaine, compared to those with prototype allele (Mann-Whitney tests, point-wise).

CONCLUSIONS

This study provides additional support of potential importance of regulatory regions of intron 2 of the OPRK1 gene in development of cocaine and opioid dependence diagnoses, in a population with African-American ancestry.

RNA Solutions: Synthesizing Information to Support Transcriptomics (RNASSIST)

MOTIVATION

Transcriptomics is a common approach to identify changes in gene expression induced by a disease state. Standard transcriptomic analyses consider differentially expressed genes (DEGs) as indicative of disease states so only a few genes would be treated as signals when the effect size is small, such as in brain tissue. For tissue with small effect sizes, if the DEGs do not belong to a pathway known to be involved in the disease, there would be little left in the transcriptome for researchers to follow up with.

RESULTS

We developed RNA Solutions: Synthesizing Information to Support Transcriptomics (RNASSIST), a new approach to identify hidden signals in transcriptomic data by linking differential expression and co-expression networks using machine learning. We applied our approach to RNA-seq data of post-mortem brains that compared the Alcohol Use Disorder (AUD) group with the control group. Many of the candidate genes are not differentially expressed so would likely be ignored by standard transcriptomic analysis pipelines. Through multiple validation strategies, we concluded that these RNASSIST-identified genes likely play a significant role in AUD.

AVAILABILITY AND IMPLEMENTATION

The RNASSIST algorithm is available at https://github.com/netrias/rnassist and both the software and the data used in RNASSIST are available at https://figshare.com/articles/software/RNAssist_Software_and_Data/16617250.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

GABAA receptor subtypes and benzodiazepine use, misuse, and abuse

Benzodiazepines have been in use for over half a century. While they remain highly prescribed, their unfavorable side-effect profile and abuse liability motivated a search for alternatives. Most of these efforts focused on the development of benzodiazepine-like drugs that are selective for specific GABA_A receptor subtypes. While there is ample evidence that subtype-selective GABA_A receptor ligands have great potential for providing symptom relief without typical benzodiazepine side-effects, it is less clear whether subtype-selective targeting strategies can also reduce misuse and abuse potential. This review focuses on the three benzodiazepine properties that are relevant to the DSM-5-TR criteria for Sedative, Hypnotic, or Anxiolytic Use Disorder, namely, reinforcing properties of benzodiazepines, maladaptive behaviors related to benzodiazepine use, and benzodiazepine tolerance and dependence. We review existing evidence regarding the involvement of different GABA_A receptor subtypes in each of these areas. The reviewed studies suggest that α1-containing GABA_A receptors play an integral role in benzodiazepine-induced plasticity in reward-related brain areas and might be involved in the development of tolerance and dependence to benzodiazepines. However, a systematic comparison of the contributions of all benzodiazepine-sensitive GABA_A receptors to these processes, a mechanistic understanding of how the positive modulation of each receptor subtype might contribute to the brain mechanisms underlying each of these processes, and a definitive answer to the question of whether specific chronic modulation of any given subtype would result in some or all of the benzodiazepine effects are currently lacking from the literature. Moreover, how non-selective benzodiazepines might lead to the maladaptive behaviors listed in DSM and how different GABA_A receptor subtypes might be involved in the development of these behaviors remains unexplored. Considering the increasing burden of benzodiazepine abuse, the common practice of benzodiazepine misuse that leads to severe dependence, and the current efforts to generate side-effect free benzodiazepine alternatives, there is an urgent need for systematic, mechanistic research that provides a better understanding of the brain mechanisms of benzodiazepine misuse and abuse, including the involvement of specific GABA_A receptor subtypes in these processes, to establish an informed foundation for preclinical and clinical efforts.

Genome-Wide Association Study on Three Behaviors Tested in an Open Field in Heterogeneous Stock Rats Identifies Multiple Loci Implicated in Psychiatric Disorders

Many personality traits are influenced by genetic factors. Rodents models provide an efficient system for analyzing genetic contribution to these traits. Using 1,246 adolescent heterogeneous stock (HS) male and female rats, we conducted a genome-wide association study (GWAS) of behaviors measured in an open field, including locomotion, novel object interaction, and social interaction. We identified 30 genome-wide significant quantitative trait loci (QTL). Using multiple criteria, including the presence of high impact genomic variants and co-localization of cis-eQTL, we identified 17 candidate genes (Adarb2, Ankrd26, Cacna1c, Cacng4, Clock, Ctu2, Cyp26b1, Dnah9, Gda, Grxcr1, Eva1a, Fam114a1, Kcnj9, Mlf2, Rab27b, Sec11a, and Ube2h) for these traits. Many of these genes have been implicated by human GWAS of various psychiatric or drug abuse related traits. In addition, there are other candidate genes that likely represent novel findings that can be the catalyst for future molecular and genetic insights into human psychiatric diseases. Together, these findings provide strong support for the use of the HS population to study psychiatric disorders.

Ibrutinib as a potential therapeutic for cocaine use disorder

Cocaine use presents a worldwide public health problem with high socioeconomic cost. No current pharmacologic treatments are available for cocaine use disorder (CUD) or cocaine toxicity. To explore pharmaceutical treatments for tthis disorder and its sequelae we analyzed gene expression data from post-mortem brain tissue of individuals with CUD who died from cocaine-related causes with matched cocaine-free controls (n = 71, Mage = 39.9, 100% male, 49% with CUD, 3 samples/brain regions). To match molecular signatures from brain pathology with potential therapeutics, we leveraged the L1000 database honing in on neuronal mRNA profiles of 825 repurposable compounds (e.g., FDA approved). We identified 16 compounds that were negatively associated with CUD gene expression patterns across all brain regions (padj < 0.05), all of which outperformed current targets undergoing clinical trials for CUD (all padj > 0.05). An additional 43 compounds were positively associated with CUD expression. We performed an in silico follow-up potential therapeutics using independent transcriptome-wide in vitro (neuronal cocaine exposure; n = 18) and in vivo (mouse cocaine self-administration; n = 12-15) datasets to prioritize candidates for experimental validation. Among these medications, ibrutinib was consistently linked with the molecular profiles of both neuronal cocaine exposure and mouse cocaine self-administration. We assessed the therapeutic efficacy of ibrutinib using the Drosophila melanogaster model. Ibrutinib reduced cocaine-induced startle response and cocaine-induced seizures (n = 61-142 per group; sex: 51% female), despite increasing cocaine consumption. Our results suggest that ibrutinib could be used for the treatment of cocaine use disorder.

Genetics of substance use disorders: a review

Substance use disorders (SUDs) are prevalent and result in an array of negative consequences. They are influenced by genetic factors (h2 = ~50%). Recent years have brought substantial progress in our understanding of the genetic etiology of SUDs and related traits. The present review covers the current state of the field for SUD genetics, including the epidemiology and genetic epidemiology of SUDs, findings from the first-generation of SUD genome-wide association studies (GWAS), cautions about translating GWAS findings to clinical settings, and suggested prioritizations for the next wave of SUD genetics efforts. Recent advances in SUD genetics have been facilitated by the assembly of large GWAS samples, and the development of state-of-the-art methods modeling the aggregate effect of genome-wide variation. These advances have confirmed that SUDs are highly polygenic with many variants across the genome conferring risk, the vast majority of which are of small effect. Downstream analyses have enabled finer resolution of the genetic architecture of SUDs and revealed insights into their genetic relationship with other psychiatric disorders. Recent efforts have also prioritized a closer examination of GWAS findings that have suggested non-uniform genetic influences across measures of substance use (e.g. consumption) and problematic use (e.g. SUD). Additional highlights from recent SUD GWAS include the robust confirmation of loci in alcohol metabolizing genes (e.g. ADH1B and ALDH2) affecting alcohol-related traits, and loci within the CHRNA5-CHRNA3-CHRNB4 gene cluster influencing nicotine-related traits. Similar successes are expected for cannabis, opioid, and cocaine use disorders as sample sizes approach those assembled for alcohol and nicotine.

The Association between a MAOB Variable Number Tandem Repeat Polymorphism and Cocaine and Opiate Addictions in Polyconsumers

Genetic analysis of the association between alcohol, cocaine, and opiate addiction and variable number tandem repeat (VNTR) polymorphisms in monoamine oxidase B (MAOB) and serotonergic 5-hydroxytryptamine (serotonin) receptor 1B and 2C (HTR1B 21 and HTR2C) pathway genes was performed in a sample of 302 polyconsumers. Our genetic association analysis revealed a significant association between a 184 base pair (bp) VNTR polymorphism in the MAOB gene and addiction to cocaine and opiates. This work highlights new genetic marker associations in cocaine and opiate polyconsumer addictions. These data help to clarify and quantify the complex role of genetics in addictive disorders, as well as their future contribution to the prevention (genetic counselling), diagnosis (genetic diagnosis of vulnerability), and treatment (pharmacogenomics) of these disorders.

Cocaine induces paradigm-specific changes to the transcriptome within the ventral tegmental area

During the initial stages of drug use, cocaine-induced neuroadaptations within the ventral tegmental area (VTA) are critical for drug-associated cue learning and drug reinforcement processes. These neuroadaptations occur, in part, from alterations to the transcriptome. Although cocaine-induced transcriptional mechanisms within the VTA have been examined, various regimens and paradigms have been employed to examine candidate target genes. In order to identify key genes and biological processes regulating cocaine-induced processes, we employed genome-wide RNA-sequencing to analyze transcriptional profiles within the VTA from male mice that underwent one of four commonly used paradigms: acute home cage injections of cocaine, chronic home cage injections of cocaine, cocaine-conditioning, or intravenous-self administration of cocaine. We found that cocaine alters distinct sets of VTA genes within each exposure paradigm. Using behavioral measures from cocaine self-administering mice, we also found several genes whose expression patterns corelate with cocaine intake. In addition to overall gene expression levels, we identified several predicted upstream regulators of cocaine-induced transcription shared across all paradigms. Although distinct gene sets were altered across cocaine exposure paradigms, we found, from Gene Ontology (GO) term analysis, that biological processes important for energy regulation and synaptic plasticity were affected across all cocaine paradigms. Coexpression analysis also identified gene networks that are altered by cocaine. These data indicate that cocaine alters networks enriched with glial cell markers of the VTA that are involved in gene regulation and synaptic processes. Our analyses demonstrate that transcriptional changes within the VTA depend on the route, dose and context of cocaine exposure, and highlight several biological processes affected by cocaine. Overall, these findings provide a unique resource of gene expression data for future studies examining novel cocaine gene targets that regulate drug-associated behaviors.

Do gene expression findings from mouse models of cocaine use recapitulate human cocaine use disorder in reward circuitry?

Animal models of drug use have investigated possible mechanisms governing human substance use traits for over 100 years. Most cross-species research on drug use/addiction examines behavioral overlap, but studies assessing neuromolecular (e.g. RNA) correspondence are lacking. Our study utilized transcriptome-wide data from the hippocampus and ventral tegmental area (VTA)/midbrain from a total of 35 human males with cocaine use disorder/controls and 49 male C57BL/6J cocaine/saline administering/exposed mice. We hypothesized differential expressed genes and systems of co-expressed genes (gene networks) would show appreciable overlap across mouse cocaine self-administration and human cocaine use disorder. We found modest, but significant relationships between differentially expressed genes associated with cocaine self-administration (short access) and cocaine use disorder within reward circuitry. Differentially expressed genes underlying models of acute cocaine exposure (cocaine), context re-exposure and cocaine + context re-exposure were not consistently associated with human CUD across brain regions. Investigating systems of co-expressed genes, we found several validated gene networks with weak to moderate conservation between cocaine/saline self-administering mice and disordered cocaine users/controls. The most conserved hippocampal and VTA gene networks demonstrated substantial overlap (2029 common genes) and included both novel and previously implicated targets for cocaine use/addiction. Lastly, we conducted (expression-based) phenome-wide association studies of the nine common hub genes across conserved gene networks. Common hub genes were associated with dopamine/serotonin function, cocaine self-administration and other relevant mouse traits. Overall, our study pinpointed and characterized conserved brain-related RNA patterns across mouse cocaine self-administration and human cocaine use disorder. We offer recommendations for future research and add to the dialogue surrounding pre-clinical animal research for human disease.