Strain- and region-specific gene expression profiles in mouse brain in response to chronic nicotine treatment

Multi-omic and multi-species meta-analyses of nicotine consumption

Cross-species translational approaches to human genomic analyses are lacking. The present study uses an integrative framework to investigate how genes associated with nicotine use in model organisms contribute to the genetic architecture of human tobacco consumption. First, we created a model organism geneset by collecting results from five animal models of nicotine exposure (RNA expression changes in brain) and then tested the relevance of these genes and flanking genetic variation using genetic data from human cigarettes per day (UK BioBank N = 123,844; all European Ancestry). We tested three hypotheses: (1) DNA variation in, or around, the 'model organism geneset' will contribute to the heritability to human tobacco consumption, (2) that the model organism genes will be enriched for genes associated with human tobacco consumption, and (3) that a polygenic score based off our model organism geneset will predict tobacco consumption in the AddHealth sample (N = 1667; all European Ancestry). Our results suggested that: (1) model organism genes accounted for ~5-36% of the observed SNP-heritability in human tobacco consumption (enrichment: 1.60-31.45), (2) model organism genes, but not negative control genes, were enriched for the gene-based associations (MAGMA, H-MAGMA, SMultiXcan) for human cigarettes per day, and (3) polygenic scores based on our model organism geneset predicted cigarettes per day in an independent sample. Altogether, these findings highlight the advantages of using multiple species evidence to isolate genetic factors to better understand the etiological complexity of tobacco and other nicotine consumption.

Discovery of a Role for Rab3b in Habituation and Cocaine Induced Locomotor Activation in Mice Using Heterogeneous Functional Genomic Analysis

Substance use disorders are prevalent and present a tremendous societal cost but the mechanisms underlying addiction behavior are poorly understood and few biological treatments exist. One strategy to identify novel molecular mechanisms of addiction is through functional genomic experimentation. However, results from individual experiments are often noisy. To address this problem, the convergent analysis of multiple genomic experiments can discern signal from these studies. In the present study, we examine genetic loci that modulate the locomotor response to cocaine identified in the recombinant inbred (BXD RI) genetic reference population. We then applied the GeneWeaver software system for heterogeneous functional genomic analysis to integrate and aggregate multiple studies of addiction genomics, resulting in the identification of Rab3b as a functional correlate of the locomotor response to cocaine in rodents. This gene encodes a member of the RAB family of Ras-like GTPases known to be involved in trafficking of secretory and endocytic vesicles in eukaryotic cells. The convergent evidence for a role of Rab3b includes co-occurrence in previously published genetic mapping studies of cocaine related behaviors; methamphetamine response and cocaine- and amphetamine-regulated transcript prepropeptide (Cartpt) transcript abundance; evidence related to other addictive substances; density of polymorphisms; and its expression pattern in reward pathways. To evaluate this finding, we examined the effect of RAB3 complex perturbation in cocaine response. B6;129-Rab3btm1Sud Rab3ctm1sud Rab3dtm1sud triple null mice (Rab3bcd -/-) exhibited significant deficits in habituation, and increased acute and repeated cocaine responses. This previously unidentified mechanism of the behavioral predisposition and response to cocaine is an example of many that can be identified and validated using aggregate genomic studies.

Neonatal Nicotine Exposure Primes Midbrain Neurons to a Dopaminergic Phenotype and Increases Adult Drug Consumption

BACKGROUND

Nicotine intake induces addiction through neuroplasticity of the reward circuitry, altering the activity of dopaminergic neurons of the ventral tegmental area. Prior work demonstrated that altered circuit activity can change neurotransmitter expression in the developing and adult brain. Here we investigated the effects of neonatal nicotine exposure on the dopaminergic system and nicotine consumption in adulthood.

METHODS

Male and female mice were used for two-bottle-choice test, progressive ratio breakpoint test, immunohistochemistry, RNAscope, quantitative polymerase chain reaction, calcium imaging, and DREADD (designer receptor exclusively activated by designer drugs)-mediated chemogenic activation/inhibition experiments.

RESULTS

Neonatal nicotine exposure potentiates drug preference in adult mice, induces alterations in calcium spike activity of midbrain neurons, and increases the number of dopamine-expressing neurons in the ventral tegmental area. Specifically, glutamatergic neurons are first primed to express transcription factor Nurr1, then acquire the dopaminergic phenotype following nicotine re-exposure in adulthood. Enhanced neuronal activity combined with Nurr1 expression is both necessary and sufficient for the nicotine-mediated neurotransmitter plasticity to occur.

CONCLUSIONS

Our findings illuminate a new mechanism of neuroplasticity by which early nicotine exposure primes the reward system to display increased susceptibility to drug consumption in adulthood.

The circadian gene Nr1d1 in the mouse nucleus accumbens modulates sociability and anxiety-related behaviour

Nuclear receptor subfamily 1, group D, member 1 (Nr1d1) (also known as Rev-erb alpha) has been linked to circadian rhythm regulation, mood-related behaviour and disorders associated with social deficits. Recent work from our laboratory found striking decreases in Nr1d1 in the nucleus accumbens (NAc) in the maternal condition and indirect evidence that Nr1d1 was interacting with numerous addiction and reward-related genes to modulate social reward. In this study, we applied our insights from the maternal state to nonparental adult mice to determine whether decreases in Nr1d1 expression in the NAc via adeno-associated viral (AAV) vectors and short hairpin RNA (shRNA)-mediated gene knockdown were sufficient to modulate social behaviours and mood-related behaviours. Knockdown of Nr1d1 in the NAc enhanced sociability and reduced anxiety, but did not affect depressive-like traits in female mice. In male mice, Nr1d1 knockdown had no significant behavioural effects. Microarray analysis of Nr1d1 knockdown in females identified changes in circadian rhythm and histone deacetylase genes and suggested possible drugs, including histone deacetylase inhibitors, that could mimic actions of Nr1d1 knockdown. Quantitative real-time PCR (qPCR) analysis confirmed expression upregulation of gene period circadian clock 1 (Per1) and period circadian clock 2 (Per2) with Nr1d1 knockdown. The evidence for roles for opioid-related genes opioid receptor, delta 1 (Oprd1) and preproenkephalin (Penk) was also found. Together, these results suggest that Nr1d1 in the NAc modulates sociability and anxiety-related behaviour in a sex-specific manner, and circadian, histone deacetylase and opioid-related genes may be involved in the expression of these behavioural phenotypes.

Genetic and neuroendocrine regulation of the postpartum brain

Changes in expression of hundreds of genes occur during the production and function of the maternal brain that support a wide range of processes. In this review, we synthesize findings from four microarray studies of different maternal brain regions and identify a core group of 700 maternal genes that show significant expression changes across multiple regions. With those maternal genes, we provide new insights into reward-related pathways (maternal bonding), postpartum depression, social behaviors, mental health disorders, and nervous system plasticity/developmental events. We also integrate the new genes into well-studied maternal signaling pathways, including those for prolactin, oxytocin/vasopressin, endogenous opioids, and steroid receptors (estradiol, progesterone, cortisol). A newer transcriptional regulation model for the maternal brain is provided that incorporates recent work on maternal microRNAs. We also compare the top 700 genes with other maternal gene expression studies. Together, we highlight new genes and new directions for studies on the postpartum brain.

Identification of Differentially Expressed Genes through Integrated Study of Alzheimer's Disease Affected Brain Regions

Background

Alzheimer’s disease (AD) is the most common form of dementia in older adults that damages the brain and results in impaired memory, thinking and behaviour. The identification of differentially expressed genes and related pathways among affected brain regions can provide more information on the mechanisms of AD. In the past decade, several studies have reported many genes that are associated with AD. This wealth of information has become difficult to follow and interpret as most of the results are conflicting. In that case, it is worth doing an integrated study of multiple datasets that helps to increase the total number of samples and the statistical power in detecting biomarkers. In this study, we present an integrated analysis of five different brain region datasets and introduce new genes that warrant further investigation.

Methods

The aim of our study is to apply a novel combinatorial optimisation based meta-analysis approach to identify differentially expressed genes that are associated to AD across brain regions. In this study, microarray gene expression data from 161 samples (74 non-demented controls, 87 AD) from the Entorhinal Cortex (EC), Hippocampus (HIP), Middle temporal gyrus (MTG), Posterior cingulate cortex (PC), Superior frontal gyrus (SFG) and visual cortex (VCX) brain regions were integrated and analysed using our method. The results are then compared to two popular meta-analysis methods, RankProd and GeneMeta, and to what can be obtained by analysing the individual datasets.

Results

We find genes related with AD that are consistent with existing studies, and new candidate genes not previously related with AD. Our study confirms the up-regualtion of INFAR2 and PTMA along with the down regulation of GPHN, RAB2A, PSMD14 and FGF. Novel genes PSMB2, WNK1, RPL15, SEMA4C, RWDD2A and LARGE are found to be differentially expressed across all brain regions. Further investigation on these genes may provide new insights into the development of AD. In addition, we identified the presence of 23 non-coding features, including four miRNA precursors (miR-7, miR570, miR-1229 and miR-6821), dysregulated across the brain regions. Furthermore, we compared our results with two popular meta-analysis methods RankProd and GeneMeta to validate our findings and performed a sensitivity analysis by removing one dataset at a time to assess the robustness of our results. These new findings may provide new insights into the disease mechanisms and thus make a significant contribution in the near future towards understanding, prevention and cure of AD.

KCNN Genes that Encode Small-Conductance Ca2+-Activated K+ Channels Influence Alcohol and Drug Addiction

Small-conductance Ca(2+)-activated K(+) (KCa2) channels control neuronal excitability and synaptic plasticity, and have been implicated in substance abuse. However, it is unknown if genes that encode KCa2 channels (KCNN1-3) influence alcohol and drug addiction. In the present study, an integrative functional genomics approach shows that genetic datasets for alcohol, nicotine, and illicit drugs contain the family of KCNN genes. Alcohol preference and dependence QTLs contain KCNN2 and KCNN3, and Kcnn3 transcript levels in the nucleus accumbens (NAc) of genetically diverse BXD strains of mice predicted voluntary alcohol consumption. Transcript levels of Kcnn3 in the NAc negatively correlated with alcohol intake levels in BXD strains, and alcohol dependence enhanced the strength of this association. Microinjections of the KCa2 channel inhibitor apamin into the NAc increased alcohol intake in control C57BL/6J mice, while spontaneous seizures developed in alcohol-dependent mice following apamin injection. Consistent with this finding, alcohol dependence enhanced the intrinsic excitability of medium spiny neurons in the NAc core and reduced the function and protein expression of KCa2 channels in the NAc. Altogether, these data implicate the family of KCNN genes in alcohol, nicotine, and drug addiction, and identify KCNN3 as a mediator of voluntary and excessive alcohol consumption. KCa2.3 channels represent a promising novel target in the pharmacogenetic treatment of alcohol and drug addiction.

Addiction and reward-related genes show altered expression in the postpartum nucleus accumbens

Motherhood involves a switch in natural rewards, whereby offspring become highly rewarding. Nucleus accumbens (NAC) is a key CNS region for natural rewards and addictions, but to date no study has evaluated on a large scale the events in NAC that underlie the maternal change in natural rewards. In this study we utilized microarray and bioinformatics approaches to evaluate postpartum NAC gene expression changes in mice. Modular Single-set Enrichment Test (MSET) indicated that postpartum (relative to virgin) NAC gene expression profile was significantly enriched for genes related to addiction and reward in five of five independently curated databases (e.g., Malacards, Phenopedia). Over 100 addiction/reward related genes were identified and these included: Per1, Per2, Arc, Homer2, Creb1, Grm3, Fosb, Gabrb3, Adra2a, Ntrk2, Cry1, Penk, Cartpt, Adcy1, Npy1r, Htr1a, Drd1a, Gria1, and Pdyn. ToppCluster analysis found maternal NAC expression profile to be significantly enriched for genes related to the drug action of nicotine, ketamine, and dronabinol. Pathway analysis indicated postpartum NAC as enriched for RNA processing, CNS development/differentiation, and transcriptional regulation. Weighted Gene Coexpression Network Analysis (WGCNA) identified possible networks for transcription factors, including Nr1d1, Per2, Fosb, Egr1, and Nr4a1. The postpartum state involves increased risk for mental health disorders and MSET analysis indicated postpartum NAC to be enriched for genes related to depression, bipolar disorder (BPD), and schizophrenia. Mental health related genes included: Fabp7, Grm3, Penk, and Nr1d1. We confirmed via quantitative PCR Nr1d1, Per2, Grm3, Penk, Drd1a, and Pdyn. This study indicates for the first time that postpartum NAC involves large scale gene expression alterations linked to addiction and reward. Because the postpartum state also involves decreased response to drugs, the findings could provide insights into how to mitigate addictions.

Prioritizing Genes Related to Nicotine Addiction Via a Multi-source-Based Approach

Nicotine has a broad impact on both the central and peripheral nervous systems. Over the past decades, an increasing number of genes potentially involved in nicotine addiction have been identified by different technical approaches. However, the molecular mechanisms underlying nicotine addiction remain largely unknown. Under such situation, prioritizing the candidate genes for further investigation is becoming increasingly important. In this study, we presented a multi-source-based gene prioritization approach for nicotine addiction by utilizing the vast amounts of information generated from for nicotine addiction study during the past years. In this approach, we first collected and curated genes from studies in four categories, i.e., genetic association analysis, genetic linkage analysis, high-throughput gene/protein expression analysis, and literature search of single gene/protein-based studies. Based on these resources, the genes were scored and a weight value was determined for each category. Finally, the genes were ranked by their combined scores, and 220 genes were selected as the prioritized nicotine addiction-related genes. Evaluation suggested the prioritized genes were promising targets for further analysis and replication study.

Potential translational targets revealed by linking mouse grooming behavioral phenotypes to gene expression using public databases

Rodent self-grooming is an important, evolutionarily conserved behavior, highly sensitive to pharmacological and genetic manipulations. Mice with aberrant grooming phenotypes are currently used to model various human disorders. Therefore, it is critical to understand the biology of grooming behavior, and to assess its translational validity to humans. The present in-silico study used publicly available gene expression and behavioral data obtained from several inbred mouse strains in the open-field, light-dark box, elevated plus- and elevated zero-maze tests. As grooming duration differed between strains, our analysis revealed several candidate genes with significant correlations between gene expression in the brain and grooming duration. The Allen Brain Atlas, STRING, GoMiner and Mouse Genome Informatics databases were used to functionally map and analyze these candidate mouse genes against their human orthologs, assessing the strain ranking of their expression and the regional distribution of expression in the mouse brain. This allowed us to identify an interconnected network of candidate genes (which have expression levels that correlate with grooming behavior), display altered patterns of expression in key brain areas related to grooming, and underlie important functions in the brain. Collectively, our results demonstrate the utility of large-scale, high-throughput data-mining and in-silico modeling for linking genomic and behavioral data, as well as their potential to identify novel neural targets for complex neurobehavioral phenotypes, including grooming.

Ryanodine receptor-2 upregulation and nicotine-mediated plasticity

Nicotine, the major psychoactive component of cigarette smoke, modulates neuronal activity to produce Ca2+-dependent changes in gene transcription. However, the downstream targets that underlie the long-term effects of nicotine on neuronal function, and hence behaviour, remain to be elucidated. Here, we demonstrate that nicotine administration to mice upregulates levels of the type 2 ryanodine receptor (RyR2), a Ca2+-release channel present on the endoplasmic reticulum, in a number of brain areas associated with cognition and addiction, notably the cortex and ventral midbrain. Nicotine-mediated RyR2 upregulation was driven by CREB, and caused a long-lasting reinforcement of Ca2+ signalling via the process of Ca2+-induced Ca2+ release. RyR2 upregulation was itself required for long-term phosphorylation of CREB in a positive-feedback signalling loop. We further demonstrate that inhibition of RyR-activation in vivo abolishes sensitization to nicotine-induced habituated locomotion, a well-characterised model for onset of drug dependence. Our findings, therefore, indicate that gene-dependent reprogramming of Ca2+ signalling is involved in nicotine-induced behavioural changes.

Detection of genetic association and a functional polymorphism of dynamin 1 gene with nicotine dependence in European and African Americans

Although it has been documented that dynamin 1 gene (DNM1) is significantly modulated by nicotine in animal models, its association with nicotine dependence (ND) in human population remained to be unexplored. To determine whether DNM1 is associated with ND, in this study, we genotyped seven single-nucleotide polymorphisms (SNPs) within this gene in 602 nuclear families of either African-American (AA) or European-American (EA) origin. Individual SNP-based association analysis revealed a significant association of SNP rs3003609 with smoking quantity (SQ; P=0.0031) and Heaviness of Smoking Index (HSI; P=0.0042) in the EA sample. Furthermore, our haplotype-based association analyses indicated that haplotypes T-G-T, formed by rs2502731-rs2229917-rs3003609 (at a frequency of 54%), G-T-A, formed by rs2229917-rs3003609-rs16930313 (at a frequency of 52%), and T-A-G, formed by rs3003609-rs16930313-rs7022174 (at a frequency of 52%) are significantly associated with SQ (Z=-2.44 to -2.92; P=0.015-0.0055) and HSI (Z=-2.52 to -2.67; P=0.012-0.0076) in the EA sample. In the AA sample, another haplotype, G-T-A, formed by rs7875406-rs2502731-rs2229917, at a frequency of 12% was significantly associated with SQ (Z=-2.58; P=0.0098). Finally, by using in vitro gene expression assays, we demonstrated that the T allele of rs3003609 in the exon 9 of DNM1 significantly decreases the expression of DNM1, by 27.1% at the mRNA and 22.0% at the protein level, suggesting that rs3003609 represents a functional polymorphism affecting DNM1 expression and may partly contributed to the observed association of the gene with ND in our samples. Taken together, our findings indicate that DNM1 is likely involved in the etiology of ND and represents a plausible candidate for further investigation in independent samples.