- Research article
- Open Access
Chronic nicotine differentially affects murine transcriptome profiling in isolated cortical interneurons and pyramidal neurons
BMC Genomics volume 18, Article number: 194 (2017)
Nicotine is known to differentially regulate cortical interneuron and pyramidal neuron activities in the neocortex, while the underlying molecular mechanisms have not been well studied. In this study, RNA-sequencing was performed in acutely isolated cortical somatostatin (Sst)- positive interneurons and pyramidal neurons (Thy1) from mice treated with systemic nicotine for 14 days. We assessed the differentially expressed genes (DEGs) by nicotine in Sst- or Thy1- neurons, respectively, and then compared DEGs between Sst- and Thy1- neurons in the absence and presence of nicotine.
In Sst-neurons, the DEGs by nicotine were associated with glycerophospholipid and nicotinate and nicotinamide metabolism; while in Thy1-neurons those related to immune response and purine and pyrimidine metabolisms were affected. Under basal condition, the DEGs between Sst- and Thy1- neurons were frequently associated with signal transduction, phosphorylation and potassium channel regulation. However, some new DEGs between Sst- and Thy1- neurons were found after nicotine, the majority of which belong to mitochondrial respiratory chain complex.
Nicotine differentially affected subset of genes in Sst- and Thy1- neurons, which might contribute to the distinct effect of nicotine on interneuron and pyramidal neuron activities. Meanwhile, the altered transcripts associated with mitochondrial activity were found between interneurons and pyramidal neurons after chronic nicotine.
The inter-play between interneurons and pyramidal neurons in the neocortex forms the basis of inhibition and excitation and neural network function [1, 2]. Interneurons powerfully regulate pyramidal neuron functions, albeit with low synaptic inputs in number . Recent studies demonstrate that GABAergic interneurons play a critical role in local circuit function and behavior [4–6]. Dysfunction of GABAergic signaling is associated with age-related cognitive decline, schizophrenia, ischemia and Alzheimer Disease [7–10].
Nicotine has been found to improve working memory and learning by activating nicotinic acetylcholine receptors (nAChRs) [11, 12]. Emerging evidence has shown that nicotine differentially regulates interneuron and pyramidal neuron activities. For instance, nicotine layer-specifically regulates neuronal activities where distinct interneurons and pyramidal neurons are located . While nAChRs enhance AMPA receptor mediated current and firing rate in interneurons [14, 15], they cause a sustained reduction of NMDA receptor mediated currents in pyramidal neurons , suggesting an involvement of different molecular mechanisms. Although gene profiling studies have made great progress in identifying individual neurons in the cortex [17, 18], and microarray studies have demonstrated that nicotine causes different gene expression in neuroblastoma cell line [19, 20] and in distinct brain regions including the cortex [21, 22], how nicotine may specifically regulate the gene expression in interneurons in relation to pyramidal neurons in cortical circuit remains unknown.
In this study, we acutely isolated cortical somatostatin (Sst) labeled interneurons and thymus cell antigen 1 (Thy1) labeled pyramidal neurons in mice treated with saline and systemic nicotine for 14 days, and the transcriptome profiling was compared in Sst- and Thy1-neuons with or without nicotine. We found that in Sst- neurons the most prominent genes affected by nicotine were associated with glycerophospholipid and nicotinate and nicotinamide metabolism, and in Thy1- neurons those associated with immune response and purine and pyrimidine metabolisms were influenced. In addition, the differentially expressed genes (DEGs) between Sst- and Thy1-neuons were associated with mitochondrial respiratory chain complex after nicotine treatment.
The following mice were used throughout the experiments: GIN mice (FVB-Tg [GadGFP] 45704Swn/J, The Jackson Laboratory, #003718) expressing enhanced green fluorescent protein (eGFP) in a group of somatostatin-positive interneurons (Sst), and thy1-YFP-H (YFPH) mice (B6.Cg-Tg [Thy1-YFP] HJrs/J, The Jackson Laboratory, #003782) expressing yellow fluorescent protein in layer V pyramidal neurons (Thy1) [23, 24]. A previous study has compared the genetic profiling between four lines of mice including GIN and YFPH used in our study. The authors found that the correlation coefficient between the samples from these two mice was not significantly different, indicating that GIN and YFPH mice do not show significant transcriptome differences in specific cell types, under basal condition . Although another major group of interneurons are parvalbumin (PV)-positive, they do not express nicotinic acetylcholine receptors in the cortex [25, 26] and thus were not chosen in this study.
Twelve male mice (20 g, 5–7weeks) were divided into four groups, with 3 GIN or YFPH mice in each group receiving nicotine or saline, respectively. Mice were implanted subcutaneously with an osmotic mini pump (model 2002; Alzet, Cupertino, CA) that delivered 48 mg/kg/d nicotine hydrogen tartrate (Sigma, St. Louis) or saline vehicle for 14 days [27, 28]. This dosage of nicotine is known to alter neuronal activity and gene expression in the brain of mice [29, 30]. Mice were housed in the 12/12 h light/dark cycle and were given food and water ad libitum. Genotypes were confirmed by PCR analysis before study. All protocols were approved by the Commission of Chongqing Medical University for ethics of experiments on animals and were in accordance with international standards.
Brain slice preparation and cell sorting
Mice were deeply anesthetized with sodium pentobarbital (100 mg/kg; i.p.) and decapitated . Brain slices (400 um) were prepared by a vibratome (Leica VP1200S, Germany) in ice-cold artificial cerebral spinal fluid (ACSF) containing (in mM):119 NaCl, 2.5 KCl, 1 MgCl2, 26 NaHCO3, 1 CaCl2, 1.25 NaH2PO4, 25 dextrose aerated with 95% O2/5% CO2 (pH7.4). The slices were then incubated in ACSF containing 50 μM APV, 20 μM DNQX and 100nM TTX for 30 min at 32 °C and then for at least 30 min at RT. The cortex were morphologically dissected and enzymatically digested by pronase E (1.5 mg/ml; Sigma Aldrich) in ACSF (containing 50 μM APV, 20 μM DNQX and 1 μM TTX) for 45 min, then triturated by three fire-polished Pasteur pipettes (400, 300 and 150 μm in inner diameter, respectively) in ACSF containing APV, DNQX, TTX and 1% BSA [31, 32]. The resultant pieces were allowed to settle for 2 min and the supernatant were collected for centrifugation at 800 g for 5 min, which removed the cellular debris and molecular contaminants . Fluorescent cells (Sst- or Thy1- positive) were carefully aspirated by a micropipette (30–50 μm) shaped by a pipette puller (P-97, Sutter, USA) . 10 cells were aspirated each time and 30 cells were collected from each mouse. Cells (10 in each) were then transferred to 3.5 μl lysis buffer containing 2% Triton X-100 and 5% RNaseOUT (20U/μl, Invitrogen) in nuclease free water and were immediately stored at − 80 °C. A total of 90 cells from 3 mice in each group were collected for final sequencing. This amount of cells has been proven to be enough to reliably measure the population transciptome .
RNA library construction and sequencing
SMARTer Ultra Low Input RNA for Illumina Sequencing kit (Clontech) was used to amplify cDNA from cell lysates, according to manufacturer’s instruction. This kit allows high-quality cDNA synthesis in a single cell containing as low as 10 pg RNA [30, 33]. The cDNA library was generated using Ovation® Ultralow DR Multiplex System (NuGEN). The amplified products were quantified by Qubit® 2.0 Fluorometer (Invitrogen) and validated by Agilent 2100 bioanalyzer (Agilent). TruSeq PE Cluster Kit v3 (Illumina) were used for cluster generation in an Illumina cBOT instrument. And the sequencing was performed on an Illumina HiSeq2500 instrument (Illumina) following the manufacturer’s protocol. Raw reads were performed with base-calling and quality-filtering process and then aligned to mouse genome (version: mmu10.p2) using the Tophat program . Differentially expressed genes (DEGs) were expressed as FPKM (Fragments Per Kilobase of exon model per Million mapped reads), which were filtered by EBSeq algorithm, after the significant analysis and false discovery rate (FDR) analysis under the following criteria: i) Fold Change > 2 or <0.5; ii) FDR < 0.05, in which the FDR was used to correct P values [35–37]. Gene ontology (GO) annotations from NCBI, UniProtand the Gene Ontology were used to elucidate the biological implications of DEGs . Pathway analysis was used to find out the significant pathway of DEGs according to KEGG (Kyoto encyclopedia of genes and genomes) database [39, 40]. Fisher’s exact and FDR tests were used to select the significant GO categories and pathways, and the threshold of significance was defined by P-value and FDR [41, 42]. The pathway and gene interaction networks were built based on KEGG database [43, 44].
Reaction mixtures contained of 5 × buffer (contained 2 mM MgSO4), 0.5 mM each of the dNTPs, 1 μM primers, 0.02U Phanta® Super Fidelity DNA Polymerase (Vazyme, Nanjing, China), and 0.5ul of the cDNA template was made from the SMARTer Kit. The thermal cycling program for the amplification was as follows: 95 °C for 3 min, 45 cycles of 95 °C for 30 s, 56 °C or 60 °C for 15 s, and 72 °C for 15 s followed by 72 °C for 7 min. The triplicate PCR products were mixed and visualized in a 1.5% agarose gels. The expression levels were visualized by a chemiluminescence system (Fusion Fx7, Fisher Biotech) and quantified by Quantity One software (Bio-Rad, CA, USA). Data were shown as means ± SEM (n = 3). The statistically significant differences between groups were accessed by paired Student’s t test using Graphpad Prism 5 software (GraphPad). Detailed sequences of the primers were shown in Additional file 1: Table S1.
Overview of RNA sequencing
All samples were subjected to massively paralleled paired-end cDNA sequencing. Before the read (sequencing fragment) mapping, clean reads were obtained from the raw reads (5GB) by removing the adaptor sequences from each library, reads with >5% ambiguous bases (noted as N) and low-quality reads containing more than 20% of bases with qualities of <20. Of all uniquely mapped reads, about 60% were aligned to the transcript exon, 10% at the intron, 25% at the UTR regions and the remaining at TES (transcription end site), TSS (transcription start site) and intergenic regions (Additional file 2: Figure S1A). Mapped reads (Additional file 1: Tables S2 & 3) were distributed consistently on the chromosomes (Additional file 2: Figure S1B-E).
To identify the purity of manual sorting, we measured the expression level of genes associated with non-neurons, such as glia, astrocyte, oligodendrocyte, microgila and red blood cells [17, 18, 45]. Generally, non-neuron marker genes such as glia marker Vim, astrocyte marker Gfap and red blood cell marker Hbb-b1 in control groups were expressed at very low level (Additional file 1: Table S4).
Nicotine induced DEGs related to different pathways in Sst- and Thy1-neurons
There were 789 and 711 DEGs (>2 fold change; FDR < 0.05) in Sst- and Thy1-neurons after nicotine treatment, respectively. Additional 20 common genes shared by Sst- and Thy1-neurons were both up- or down- regulated by nicotine (Additional file 3: Supplemental Excel S1). Pathway interaction network analysis revealed that 24 pathways were affected by nicotine in Sst-neurons (Fig. 1a, Additional file 4: Supplemental Excel S2). Not surprisingly, nicotine significantly activated NFκB signaling as previous reported . The genes associated with cancer and neuroactive ligand-receptor interaction were upregulated while those in metabolic pathways were downregulated by nicotine (Fig. 1a). As shown in Fig. 1b, an integrated pathway tree was depicted according to KEGG database, in which hyperactive pathways associated with glycerophospholipid metabolism and nicotinate and nicotinamide metabolism (red) and hypoactive metabolic pathways (blue) were shown. Ppap2b, Pld2, Pld1, Mboat1 and Lpl were grouped into glycerophospholipid metabolism, and Nmnat3, Nudt12 and Nmnat1 were categorized into nicotinate and nicotinamide metabolism, which were all activated by nicotine (Fig. 1c). Among these genes, Pld1, Pld2 and Paap2b owned the strongest degree of centrality, suggesting that glycerophospholipid metabolism has a powerful role in Sst- neurons after nicotine.
In Thy1-neurons, nicotine affected 12 pathways related to ABC transporters, calcium signaling, cytokine-cytokine receptor interaction, rheumatoid arthritis, natural killer cell mediated cytotoxicity and pyrimidine and purine metabolism (Fig. 2a, Additional file 5: Supplemental Excel S3). It was interesting to note that pathways in cytokine-cytokine receptor interaction, rheumatoid arthritis, and natural killer cell mediated cytotoxicity are all associated with immune response. According to KEGG database, gene networks related to immune response (Fig. 2b) and pyrimidine and purine metabolism (Fig. 2c) were built. Among these genes, Csf1r and Flt1, Nt5e and Nme4 were considered as the highly-connected genes in immune response (Csf1r and Flt1) and pyrimidine and purine metabolism (Nt5e and Nme4), respectively. These results indicated that in Thy1-neuons, nicotine significantly increased immune response genes while decreased pyrimidine and purine metabolism related genes.
DEGs between Sst- and Thy1-neurons under basal condition
There were total 3185 DEGs between Sst- and Thy1-neurons without nicotine treatment (Additional file 6: Supplemental Excel S4). GO analysis revealed that these DEGs were significantly enriched in intracellular signal transduction (P = 4.4 × 10−9), nervous system development (2.16 × −7) and potassium ion transport (2.66 × −6). Small GTPase mediated signal transduction (P = 7.8 × 10−7) and phosphorylation (P = 3.7 × 10−6) were also significant (Additional file 7: Supplemental Excel S5). The representative genes with the lowest FDR involved in these biological processes were listed in Fig. 3. As previously known that nervous system development genes such as neuronal differentiation 1/2/6 (Neurod1/2/6) were associated with glutamatergic neuron development in axon outgrowth and glutamatergic synaptogenesis [47, 48]. Moreover, potassium ion channel subtype genes Kcng1, Kcng2, Kcnq5, Kcnip3, Kcnh7, Kcnj4, Kcnj6, Kcns1, Kcnv1 and Kcnh3 were significantly higher in Thy1-neurons, while Kcnt2, Kcnj12, Kcnn1 and Kcnq4 were significantly higher Sst-neurons. Some of the DEGs might be cell-type specific. The genes with high expression pattern (FPKM >1000) relative to low expression pattern (Fold change >20, FDR = 0) were considered as specific genes . The top 80 Sst- and Thy1-neuron specific genes were shown in Additional file 8: Figure S2. Many of these were expressed in a manner consistent with previously published observations [17, 45], while some genes such as Thsd7a in Sst-neurons and Kcnh7 and Crym in Thy1-neurons were not previously reported, suggesting that these genes might represent new marker genes in these neurons.
Pathway analysis of DEGs between Sst- and Thy1- neurons in the absence and presence of nicotine
To further assess nicotine effect on DEGs between Sst- and Thy1-neurons, we performed pathway analysis according to KEGG data base in the absence and presence of nicotine (Additional file 9: Supplemental Excel S6). As shown in Fig. 4a & b, the DGEs between Sst- and Thy1-neurons associated with many common pathways remained exist regardless of nicotine application, which included retrograde endocannabinoid signaling, calcium signaling, axon guidance, glutamatergic synapse, cholinergic synapse, morphine addiction and GABAergic synapse, indicating that DEGs related to these pathways were not affected by nicotine. However, some new DEGs associated with Alzheimer Disease (AD, P = 3.57E-04), oxidative phosphorylation (P = 2.43E-03), Parkinson Disease (PD, P = 2.43E-03) and metabolic pathway (P = 3.37E-3) were present only in nicotine treated mice (marked by the rectangles in Fig. 4b). For most genes associated with AD and PD, significant changes were absent between Sst- and Thy1- neurons (Sst vs. Thy1, P > 0.05) in normal condition. After nicotine treatment however, these genes were significantly down-regulated in Sst-neurons compared with those in Thy1- neurons (Sst + Nic vs.Thy1 + Nic, Fig. 5a). Among these significantly changed genes, 62.5% in AD pathway and 74.1% in PD pathway were overlapped and were mitochondrial respiratory chain complex genes, which are closely associated with oxidative phosphorylation. In metabolic pathways, some top-ranked genes were shown in Fig. 5b. Different from AD/PD pathway genes which showed a clear trend of decrease (Sst vs. Thy1) in the presence of nicotine, metabolic pathway genes exhibited a complex profile. Both significantly decreased and increased genes in Sst neurons relative to Thy1 neurons were found after nicotine (Sst + Nic vs.Thy1 + Nic, Fig. 5b). Collectively, these results suggested that the relative alteration of mitochondrial respiratory chain complex genes between Sst- and Thy1- neurons may contribute to nicotinic regulation of neural network function, in which metabolic pathway may also play a role.
Validation of representative genes by RT-PCR
To validate RNA-seq findings, we examined the expression level of some hub genes described in Figs. 1 & 2 by RT-PCR. As expected, the expression of Dlx1 and Gdpd3 were very low in Thy1-neuron relative to the high level in Sst-neurons. And Dlx1 that is necessary for interneuron differentiation and migration [45, 49] was selectively expressed in Sst- neurons (Fig. 6a). Consistent with RNA-seq findings, Ppap2b, Pld1 and Crls1 were significantly up-regulated by nicotine in Sst neurons. In addition, Flt1 was up-regulated and Entpd3 and Nme4 were down-regulated by nicotine in Thy1-neurons (Fig. 6b & c).
Previous microarray studies have shown that in SH-SY5Y cells, brief exposure (1 h) of nicotine (1 mM) results in DEGs associated with 14 pathways, in which the Toll-like receptor and death receptor pathways involved in immune response were significantly affected . On the other hand, chronic nicotine causes regional DEGs, with prefrontal cortex (PFC) and nucleus accumbens (NAc) being most responsive to nicotine. The DEGs induced by nicotine in these areas are involved in phosotidylinositol signaling, calcium homeostasis and neuroprotection . In addition, a comparative study reveals that genes involved in protein modification are altered by systemic nicotine in a region specific manner . In the present study, we provide evidence that the effect of nicotine on DEGs is also dependent on distinct cell types. In addition, the relative DEGs between interneuron and pyramidal neurons might help understand how nicotine regulates brain functions at circuit level.
Effect of nicotine on lipid metabolism in interneurons
In our work, Pld1 (phospholipase D1), Pld2 and Ppap2b (phosphatidic acid phosphatase 2b, LPP3) in interneurons are significantly increased by nicotine. These results are consistent with the notion that glycerophospholipid signaling is affected by nicotine . Neuronal PLD can be activated by variety of stimulators such as neurotransmitters, hormones and growth factors . The basic functions of PLD are associated with vesicular trafficking, brain development and neuroprotection . However, different roles of PLD1 and PLD2 in Alzheimer Disease have been found. While PLD1 acts as a negative regulator of β-amyloid (Aβ) formation , PLD2 is required for Aβ-induced synaptotoxic action . PPAP2b is a lipid phosphohydrolase enzyme that catalyzes the conversion of phosphatidic acid (PA) to diacylglycerol (DAG). It also hydrolyzes lysophosphatidic acid (LPA), ceramide-1-phosphate (C1P) and sphingosine-1-phosphate (S1P), thus participating in a variety of cellular signaling . PPAP2b plays a key role in neuronal development through S1P signaling . However, the function of PPAP2b in interneurons is poorly understood.
Effect of nicotine on immune response and purine and pyrimidine metabolism in pyramidal neurons
An important pathway affected by nicotine is immune response. Two genes upregulated by nicotine are Csf1r (macrophage colony stimulating factor 1 receptor, M-CSFR) and Flt1 (fms-related tyrosine kinase 1, VEGFR1). CSF1R can be activated by CSF1 and interleukin-34 (IL-34), thus contributing to innate immunity by regulating the development of macrophage and microglia in the brain . FLT1 regulates VEGF mediated angiogenesis . It also binds to PlGF (placenta growth factor) that acts as a cytokine and is neuroprotective for cortical neurons [59, 60]. These results are consistent with the findings that nicotine regulates innate immune response in neuronal cells [19, 61]. The immune responses in the central nervous system are thought to modulate endocrine activity that controls cell migration, thermoregulation, drinking and feeding, among others . Furthermore, the chemokine mediated neuo-glial crosstalk plays important role in multiple sclerosis and AD .
Systemic application of nicotine leads to down-regulated genes associated with purine and pyrimidine metabolism, in which Nme4 (nucleoside diphosphate kinase D, NDPKD) and Nt5e (5′-nucleotidase, ecto, CD73) are considered to be highly-connected genes. NME4 binds to mitochondrial inner membrane through cardiolipin and is associated with short chain fatty acids metabolism, Kreb cycle and apoptosis . NT5E is a membrane-anchored protein that catalyzes the extracellular formation of adenosine from AMP. This protein mediates the inhibition of hippocampal synaptic plasticity and nociception [65, 66]. Although the functional consequences of these changes are not well understood, a recent study demonstrates that these genes are abnormally expressed in the brain of Parkinson disease .
Intrinsic transcriptome differences between Sst- and Thy1- neurons
Our study demonstrates that many intrinsic differences between interneurons and pyramidal neurons remain exist regardless of nicotine administration. Consistent with previous reports [17, 32], the majority of DEGs between Sst- and Thy1- neurons in our study are associated with intracellular signal transduction, nervous system development and ion channel expression and regulation (Fig. 3). Pathway analysis links these genes to neurotransmitters (glutamate and GABA), neuromodulators (e.g. dopamine, serotonin and acetylcholine) and drug addiction relative to cocaine, morphine and nicotine (Fig. 4a), suggesting that these DEGs might contribute to the integrated function of cortical interneurons and pyramidal neurons . In addition, we also find some new cell-specific genes, which include Thsd7a (thrombospondin, Type I, domain containing 7A) in Sst-neurons and Kcnh7 (potassium voltage-gated channel subfamily H member 7, HERG-3) and Crym (crystallin Mu, NADP-regulated thyroid-hormone-binding protein) in Thy1-neurons. The functional role of these genes is currently unknown.
Nicotine induces relative DEGs associated with mitochondrial respiratory chain between Sst- and Thy1- neurons
Major functions of mitochondria in neurons include the regulation of synaptic plasticity [68, 69]. Our study demonstrates that chronic nicotine causes significantly decreased genes of mitochondrial respiratory chain in interneurons relative to pyramidal neurons, suggesting that nicotine favors pyramidal activity. Consistently, ectrophysiological studies have demonstrated that chronic exposure to nicotine results in persistent depression of interneuron, while glutamatergic neurotransmission is always increased [70, 71]. At circuit level, nicotine induces gamma oscillations in hippocampal neurons [72, 73], which are associated with the integrated function of pyramidal neurons and interneurons . Thus, the relative DEGs between Sst- and Thy1- neurons might play a key role in nicotinic regulation of synaptic plasticity and network function. Unlike mitochondrial genes which show a clear trend of relative decrease in interneurons after nicotine, both the decreased and the increased genes are found in metabolic pathway, suggesting the complicated regulation. Example genes in metabolic pathway include nsdhl (NAD (P) dependent steroid dehydrogenase-like), mthfd1l (methylenetetrahydrofolate dehydrogenase (NADP+ dependent) 1-Like) and dhrs3 (dehydrogenase 3). However, the functional role of these genes in neurons is not well-understood.
Is this study, most DEGs by nicotine are not enriched in Sst- or Thy1- neurons (FPKM <1000, Additional file 4: Supplemental Excel S2 & Additional file 5: Supplemental Excel S3), suggesting that nicotine might not have significant effect on major transcripts . However, the relative changes of mitochondrial genes between Sst- and Thy1- neurons are highly enriched (FPKM >1000, Additional file 9: Supplemental Excel S6), which implies that nicotine may play a more prominent role in the regulation of functional balance between interneurons and pyramidal neurons. In addition, although many of the genes affected by nicotine can be categorized into pathways, still the majority of genes cannot be grouped, which does not rule out their functional importance. Moreover, the limitation of this study is that nicotinic effect was assessed from two genetically different mouse strains, thus the potential bias of DEG results between Sst- and Thy1- neurons may be existed (Fig. 3), as compared to the findings by Sugino K and colleagues . Nonetheless, our study highlights the following findings: (1) Interneurons instead of pyramidal neurons might play a dominant role in nicotinic regulation of glycerophospholipid signaling in specific brain regions . (2) Pyramidal neurons might be important in nicotinic regulation of immune response and calcium signaling [61, 75]. (3) The relative alterations of mitochondria related genes may contribute to nicotinic regulation of synaptic activity and neural network function [68, 73]. However, the role of the highly connected genes affected by nicotine in individual neurons is currently unknown, which remains to be studied in the future.
Artificial cerebral spinal fluid
Bovine Serum Albumin
Differentially expressed genes
Fragments Per Kilobase of exon model per Million mapped reads
Gamma aminobutyric acid
Kyoto encyclopedia of genes and genomes
Nicotinic acetylcholine receptors
N-methyl-D-aspartic acid receptor
Thymus cell antigen 1
Nieuwenhuys R. The neocortex. An overview of its evolutionary development, structural organization and synaptology. Anat Embryol (Berl). 1994;190(4):307–37.
Silberberg G. Polysynaptic subcircuits in the neocortex: spatial and temporal diversity. Curr Opin Neurobiol. 2008;18(3):332–7.
Trevelyan AJ, Watkinson O. Does inhibition balance excitation in neocortex? Prog Biophys Mol Biol. 2005;87(1):109–43.
Hangya B, Pi HJ, Kvitsiani D, Ranade SP, Kepecs A. From circuit motifs to computations: mapping the behavioral repertoire of cortical interneurons. Curr Opin Neurobiol. 2014;26:117–24.
Lovett-Barron M, Losonczy A. Behavioral consequences of GABAergic neuronal diversity. Curr Opin Neurobiol. 2014;26:27–33.
Wester JC, McBain CJ. Behavioral state-dependent modulation of distinct interneuron subtypes and consequences for circuit function. Curr Opin Neurobiol. 2014;29:118–25.
Verret L, Mann EO, Hang GB, Barth AM, Cobos I, Ho K, Devidze N, Masliah E, Kreitzer AC, Mody I, et al. Inhibitory interneuron deficit links altered network activity and cognitive dysfunction in Alzheimer model. Cell. 2012;149(3):708–21.
Barth AM, Mody I. Changes in hippocampal neuronal activity during and after unilateral selective hippocampal ischemia in vivo. J Neurosci. 2011;31(3):851–60.
McQuail JA, Frazier CJ, Bizon JL. Molecular aspects of age-related cognitive decline: the role of GABA signaling. Trends Mol Med. 2015;21(7):450–60.
Tse MT, Piantadosi PT, Floresco SB. Prefrontal cortical gamma-aminobutyric acid transmission and cognitive function: drawing links to schizophrenia from preclinical research. Biol Psychiatry. 2015;77(11):929–39.
Levin ED, McClernon FJ, Rezvani AH. Nicotinic effects on cognitive function: behavioral characterization, pharmacological specification, and anatomic localization. Psychopharmacology (Berl). 2006;184(3-4):523–39.
Mudo G, Belluardo N, Fuxe K. Nicotinic receptor agonists as neuroprotective/neurotrophic drugs. Progress in molecular mechanisms. J Neural Transm. 2007;114(1):135–47.
Poorthuis RB, Bloem B, Verhoog MB, Mansvelder HD. Layer-specific interference with cholinergic signaling in the prefrontal cortex by smoking concentrations of nicotine. J Neurosci. 2013;33(11):4843–53.
Chen GJ, Xiong Z, Yan Z. Abeta impairs nicotinic regulation of inhibitory synaptic transmission and interneuron excitability in prefrontal cortex. Mol Neurodegener. 2013;8:3.
Tang B, Luo D, Yang J, Xu XY, Zhu BL, Wang XF, Yan Z, Chen GJ. Modulation of AMPA receptor mediated current by nicotinic acetylcholine receptor in layer I neurons of rat prefrontal cortex. Sci Rep. 2015;5:14099.
Chen G, Chen P, Tan H, Ma D, Dou F, Feng J, Yan Z. Regulation of the NMDA receptor-mediated synaptic response by acetylcholinesterase inhibitors and its impairment in an animal model of Alzheimer’s disease. Neurobiol Aging. 2008;29(12):1795–804.
Sugino K, Hempel CM, Miller MN, Hattox AM, Shapiro P, Wu C, Huang ZJ, Nelson SB. Molecular taxonomy of major neuronal classes in the adult mouse forebrain. Nat Neurosci. 2006;9(1):99–107.
Cahoy JD, Emery B, Kaushal A, Foo LC, Zamanian JL, Christopherson KS, Xing Y, Lubischer JL, Krieg PA, Krupenko SA, et al. A transcriptome database for astrocytes, neurons, and oligodendrocytes: a new resource for understanding brain development and function. J Neurosci. 2008;28(1):264–78.
Cui W-Y, Wang J, Wei J, Cao J, Chang SL, Gu J, Li MD. Modulation of innate immune-related pathways in nicotine-treated SH-SY5Y cells. Amino Acids. 2012;43(3):1157–69.
Dunckley T, Lukas RJ. Nicotinic modulation of gene expression in SH-SY5Y neuroblastoma cells. Brain Res. 2006;1116(1):39–49.
Kane JK, Konu O, Ma JZ, Li MD. Nicotine coregulates multiple pathways involved in protein modification/degradation in rat brain. Brain Res Mol Brain Res. 2004;132(2):181–91.
Mansvelder HD, Mertz M, Role LW. Nicotinic modulation of synaptic transmission and plasticity in cortico-limbic circuits. Semin Cell Dev Biol. 2009;20(4):432–40.
Porrero C, Rubio-Garrido P, Avendano C, Clasca F. Mapping of fluorescent protein-expressing neurons and axon pathways in adult and developing Thy1-eYFP-H transgenic mice. Brain Res. 2010;1345:59–72.
Feng G, Mellor RH, Bernstein M, Keller-Peck C, Nguyen QT, Wallace M, Nerbonne JM, Lichtman JW, Sanes JR. Imaging neuronal subsets in transgenic mice expressing multiple spectral variants of GFP. Neuron. 2000;28(1):41–51.
Rudy B, Fishell G, Lee S, Hjerling-Leffler J. Three groups of interneurons account for nearly 100% of neocortical GABAergic neurons. Dev Neurobiol. 2011;71(1):45–61.
Couey JJ, Meredith RM, Spijker S, Poorthuis RB, Smit AB, Brussaard AB, Mansvelder HD. Distributed network actions by nicotine increase the threshold for spike-timing-dependent plasticity in prefrontal cortex. Neuron. 2007;54(1):73–87.
Hilario MR, Turner JR, Blendy JA. Reward sensitization: effects of repeated nicotine exposure and withdrawal in mice. Neuropsychopharmacology. 2012;37(12):2661–70.
Nashmi R, Xiao C, Deshpande P, McKinney S, Grady SR, Whiteaker P, Huang Q, McClure-Begley T, Lindstrom JM, Labarca C, et al. Chronic nicotine cell specifically upregulates functional alpha 4* nicotinic receptors: basis for both tolerance in midbrain and enhanced long-term potentiation in perforant path. J Neurosci. 2007;27(31):8202–18.
Xiao C, Nashmi R, McKinney S, Cai H, McIntosh JM, Lester HA. Chronic nicotine selectively enhances alpha4beta2* nicotinic acetylcholine receptors in the nigrostriatal dopamine pathway. J Neurosci. 2009;29(40):12428–39.
Henley BM, Williams BA, Srinivasan R, Cohen BN, Xiao C, Mackey ED, Wold BJ, Lester HA. Transcriptional regulation by nicotine in dopaminergic neurons. Biochem Pharmacol. 2013;86(8):1074–83.
Brewer GJ, Torricelli JR. Isolation and culture of adult neurons and neurospheres. Nat Protoc. 2007;2(6):1490–8.
Hempel CM, Sugino K, Nelson SB. A manual method for the purification of fluorescently labeled neurons from the mammalian brain. Nat Protoc. 2007;2(11):2924–9.
Wu AR, Neff NF, Kalisky T, Dalerba P, Treutlein B, Rothenberg ME, Mburu FM, Mantalas GL, Sim S, Clarke MF, et al. Quantitative assessment of single-cell RNA-sequencing methods. Nat Methods. 2014;11(1):41–6.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Leng N, Dawson JA, Thomson JA, Ruotti V, Rissman AI, Smits BM, Haag JD, Gould MN, Stewart RM, Kendziorski C. EBSeq: an empirical Bayes hierarchical model for inference in RNA-seq experiments. Bioinformatics. 2013;29(8):1035–43.
Miao X, Luo Q, Qin X, Guo Y, Zhao H. Genome-wide mRNA-seq profiling reveals predominant down-regulation of lipid metabolic processes in adipose tissues of Small Tail Han than Dorset sheep. Biochem Biophys Res Commun. 2015;467(2):413–20.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
Draghici S, Khatri P, Tarca AL, Amin K, Done A, Voichita C, Georgescu C, Romero R. A systems biology approach for pathway level analysis. Genome Res. 2007;17(10):1537–45.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32(Database issue):D277–80.
Wang YF, Li J, He YZ, Yu HQ, Li Y, Gu XD, Li W, Li HW. iPSCs are transcriptionally and post-transcriptionally indistinguishable from fESCs. Front Biosci (Landmark Ed). 2012;17:1659–68.
Hu S, Yao G, Wang Y, Xu H, Ji X, He Y, Zhu Q, Chen Z, Sun Y. Transcriptomic changes during the pre-receptive to receptive transition in human endometrium detected by RNA-Seq. J Clin Endocrinol Metab. 2014;99(12):E2744–53.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Li C, Li H. Network-constrained regularization and variable selection for analysis of genomic data. Bioinformatics. 2008;24(9):1175–82.
Zeisel A, Munoz-Manchado AB, Codeluppi S, Lonnerberg P, La Manno G, Jureus A, Marques S, Munguba H, He L, Betsholtz C, et al. Brain structure. Cell types in the mouse cortex and hippocampus revealed by single-cell RNA-seq. Science. 2015;347(6226):1138–42.
Crowley-Weber CL, Dvorakova K, Crowley C, Bernstein H, Bernstein C, Garewal H, Payne CM. Nicotine increases oxidative stress, activates NF-kappaB and GRP78, induces apoptosis and sensitizes cells to genotoxic/xenobiotic stresses by a multiple stress inducer, deoxycholate: relevance to colon carcinogenesis. Chem Biol Interact. 2003;145(1):53–66.
Roybon L, Mastracci TL, Li J, Stott SR, Leiter AB, Sussel L, Brundin P, Li JY. The Origin, Development and Molecular Diversity of Rodent Olfactory Bulb Glutamatergic Neurons Distinguished by Expression of Transcription Factor NeuroD1. PLoS ONE. 2015;10(6), e0128035.
Bormuth I, Yan K, Yonemasu T, Gummert M, Zhang M, Wichert S, Grishina O, Pieper A, Zhang W, Goebbels S, et al. Neuronal basic helix-loop-helix proteins Neurod2/6 regulate cortical commissure formation before midline interactions. J Neurosci. 2013;33(2):641–51.
Wonders C, Anderson S. Beyond migration: Dlx1 regulates interneuron differentiation. Nat Neurosci. 2005;8(8):979–81.
Li MD, Konu O, Kane JK, Becker KG. Microarray technology and its application on nicotine research. Mol Neurobiol. 2002;25(3):265–85.
Klein J. Functions and pathophysiological roles of phospholipase D in the brain. J Neurochem. 2005;94(6):1473–87.
Kanaho Y, Funakoshi Y, Hasegawa H. Phospholipase D signalling and its involvement in neurite outgrowth. Biochim Biophys Acta. 2009;1791(9):898–904.
Cai D, Netzer WJ, Zhong M, Lin Y, Du G, Frohman M, Foster DA, Sisodia SS, Xu H, Gorelick FS, et al. Presenilin-1 uses phospholipase D1 as a negative regulator of beta-amyloid formation. Proc Natl Acad Sci U S A. 2006;103(6):1941–6.
Oliveira TG, Chan RB, Tian H, Laredo M, Shui G, Staniszewski A, Zhang H, Wang L, Kim TW, Duff KE, et al. Phospholipase d2 ablation ameliorates Alzheimer’s disease-linked synaptic dysfunction and cognitive deficits. J Neurosci. 2010;30(49):16419–28.
Brindley DN, Pilquil C. Lipid phosphate phosphatases and signaling. J Lipid Res. 2009;50(Suppl):S225–30.
Lopez-Juarez A, Morales-Lazaro S, Sanchez-Sanchez R, Sunkara M, Lomeli H, Velasco I, Morris AJ, Escalante-Alcalde D. Expression of LPP3 in Bergmann glia is required for proper cerebellar sphingosine-1-phosphate metabolism/signaling and development. Glia. 2011;59(4):577–89.
Stanley ER, Chitu V. CSF-1 Receptor Signaling in Myeloid Cells. Cold Spring Harb Perspect Biol. 2014;6:a021857.
Shibuya M. VEGFR and type-V RTK activation and signaling. Cold Spring Harb Perspect Biol. 2013;5(10):a009092.
Dewerchin M, Carmeliet P. PlGF: a multitasking cytokine with disease-restricted activity. Cold Spring Harb Perspect Med. 2012;2:8.
Du H, Li P, Pan Y, Li W, Hou J, Chen H, Wang J, Tang H. Vascular endothelial growth factor signaling implicated in neuroprotective effects of placental growth factor in an in vitro ischemic model. Brain Res. 2010;1357:1–8.
Cui WY, Li MD. Nicotinic modulation of innate immune pathways via alpha7 nicotinic acetylcholine receptor. J Neuroimmune Pharmacol. 2010;5(4):479–88.
Rostene W, Guyon A, Kular L, Godefroy D, Barbieri F, Bajetto A, Banisadr G, Callewaere C, Conductier G, Rovere C, et al. Chemokines and chemokine receptors: new actors in neuroendocrine regulations. Front Neuroendocrinol. 2011;32(1):10–24.
Cartier L, Hartley O, Dubois-Dauphin M, Krause KH. Chemokine receptors in the central nervous system: role in brain inflammation and neurodegenerative diseases. Brain Res Brain Res Rev. 2005;48(1):16–42.
Lacombe M-L, Tokarska-Schlattner M, Epand RF, Boissan M, Epand RM, Schlattner U. Interaction of NDPK-D with cardiolipin-containing membranes: Structural basis and implications for mitochondrial physiology. Biochimie. 2009;91(6):779–83.
Zhang D, Xiong W, Chu S, Sun C, Albensi BC, Parkinson FE. Inhibition of hippocampal synaptic activity by ATP, hypoxia or oxygen-glucose deprivation does not require CD73. PLoS ONE. 2012;7(6), e39772.
Sowa NA, Taylor-Blake B, Zylka MJ. Ecto-5′-nucleotidase (CD73) inhibits nociception by hydrolyzing AMP to adenosine in nociceptive circuits. J Neurosci. 2010;30(6):2235–44.
Garcia-Esparcia P, Hernandez-Ortega K, Ansoleaga B, Carmona M, Ferrer I. Purine metabolism gene deregulation in Parkinson’s disease. Neuropathol Appl Neurobiol. 2015.
Mattson MP, Gleichmann M, Cheng A. Mitochondria in neuroplasticity and neurological disorders. Neuron. 2008;60(5):748–66.
Mattson MP. Mitochondrial regulation of neuronal plasticity. Neurochem Res. 2007;32(4-5):707–15.
Mansvelder HD, Keath JR, McGehee DS. Synaptic Mechanisms Underlie Nicotine-Induced Excitability of Brain Reward Areas. Neuron. 2002;33(6):905–19.
Pidoplichko VI, Noguchi J, Areola OO, Liang Y, Peterson J, Zhang T, Dani JA. Nicotinic cholinergic synaptic mechanisms in the ventral tegmental area contribute to nicotine addiction. Learn Mem. 2004;11(1):60–9.
Wang Y, Wang Z, Wang J, Henderson Z, Wang X, Zhang X, Song J, Lu C. The modulation of nicotinic acetylcholine receptors on the neuronal network oscillations in rat hippocampal CA3 area. Sci Rep. 2015;5:9493.
Zhang X, Ge XY, Wang JG, Wang YL, Wang Y, Yu Y, Li PP, Lu CB. Induction of long-term oscillations in the gamma frequency band by nAChR activation in rat hippocampal CA3 area. Neuroscience. 2015;301:49–60.
Kann O, Papageorgiou IE, Draguhn A. Highly energized inhibitory interneurons are a central element for information processing in cortical networks. J Cereb Blood Flow Metab. 2014;34(8):1270–82.
Quik M, Perez XA, Bordia T. Nicotine as a potential neuroprotective agent for Parkinson’s disease. Mov Disord. 2012;27(8):947–57.
This work was supported by NSFC grants (numbers 81171197 & 81220108010) to G-J. C.
Availability of data and materials
The NCBI-GEO accession number for the data is GSE89899. All other data supporting the findings is contained within the manuscript.
G-JC and ZY designed research; JY and A-YL performed research; BT and DL helped with brain slice preparation and neuronal isolation; Y-JL helped with PCR; B-LZ and X-FW helped with paper writing and graph making; JY analyzed data; JY and G-JC wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
In this study, the animals used were mice. No human materials were included. All protocols were approved by the Commission of Chongqing Medical University for ethics of experiments on animals and were in accordance with international standards.
Paired primer sequences used for RT-PCR. Table S2. Clean Summary. Table S3. Mapping result statistics. Table S4. Gene expression levels associated with non-neurons. (DOC 101 kb)
Overview of RNA-sequencing quality. A: Percentage of mapped reads onto the regions of exon, intron, 5′-untranslated region (5′-UTR), 3′-untranslated region (3′-UTR), transcription start site (TSS), transcriptionend site (TES) and intergenic region (InterGenic) in all samples. B-E: Distribution of reads on chromosomes in independent samples of Sst- and Thy1-neurons, in the absence (- Nic) and presence of nicotine (+ Nic). (TIF 8790 kb)
List of DEGs affected by nicotine in Sst- and Thy1-neurons. (XLSX 203 kb)
Pathways and related genes affected by nicotine in Sst neurons. (XLSX 33 kb)
Pathways and related genes affected by nicotine in Thy1 neurons. (XLSX 23 kb)
List of DEGs between Sst- and Thy1-neurons in the absence of nicotine. (XLSX 392 kb)
GO analysis of DEGs between Sst- and Thy1-neurons in the absence of nicotine. (XLSX 76 kb)
Heat map of top 80 cell-type-specific genes expressed in Sst- and Thy1-neurons. The genes with high expression pattern (FPKM > 1000, top and bottom) relative to low expression pattern (Fold change >20, FDR = 0) are selected. For low-expression genes, those with FPKM > 10 are shown on the top, while those with FPKM < 10 are shown on the bottom (n = 3). Gene symbols with expression pattern consistent with published observations are shown in red. (TIF 1891 kb)
Comparisons of pathways in the absence and presence of nicotine. (XLSX 172 kb)
About this article
Cite this article
Yang, J., Liu, AY., Tang, B. et al. Chronic nicotine differentially affects murine transcriptome profiling in isolated cortical interneurons and pyramidal neurons. BMC Genomics 18, 194 (2017). https://doi.org/10.1186/s12864-017-3593-x
- Pyramidal neuron