- Research article
- Open Access
Behavioral sensitization induced by methamphetamine causes differential alterations in gene expression and histone acetylation of the prefrontal cortex in rats
BMC Neuroscience volume 22, Article number: 24 (2021)
Methamphetamine (METH) is one of the most widely abused illicit substances worldwide; unfortunately, its addiction mechanism remains unclear. Based on accumulating evidence, changes in gene expression and chromatin modifications might be related to the persistent effects of METH on the brain. In the present study, we took advantage of METH-induced behavioral sensitization as an animal model that reflects some aspects of drug addiction and examined the changes in gene expression and histone acetylation in the prefrontal cortex (PFC) of adult rats.
We conducted mRNA microarray and chromatin immunoprecipitation (ChIP) coupled to DNA microarray (ChIP-chip) analyses to screen and identify changes in transcript levels and histone acetylation patterns. Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, were performed to analyze the differentially expressed genes. We then further identified alterations in ANP32A (acidic leucine-rich nuclear phosphoprotein-32A) and POU3F2 (POU domain, class 3, transcription factor 2) using qPCR and ChIP-PCR assays.
In the rat model of METH-induced behavioral sensitization, METH challenge caused 275 differentially expressed genes and a number of hyperacetylated genes (821 genes with H3 acetylation and 10 genes with H4 acetylation). Based on mRNA microarray and GO and KEGG enrichment analyses, 24 genes may be involved in METH-induced behavioral sensitization, and 7 genes were confirmed using qPCR. We further examined the alterations in the levels of the ANP32A and POU3F2 transcripts and histone acetylation at different periods of METH-induced behavioral sensitization. H4 hyperacetylation contributed to the increased levels of ANP32A mRNA and H3/H4 hyperacetylation contributed to the increased levels of POU3F2 mRNA induced by METH challenge-induced behavioral sensitization, but not by acute METH exposure.
The present results revealed alterations in transcription and histone acetylation in the rat PFC by METH exposure and provided evidence that modifications of histone acetylation contributed to the alterations in gene expression caused by METH-induced behavioral sensitization.
Methamphetamine (METH) has been one of the most widely abused illicit substances worldwide for many years. Abusers exhibit deficits in cognition, behavior and emotions with long-term use and have a high METH relapse potential. The neurobiological mechanism of METH addiction is attributed to adaptations in signaling cascades and neural plasticity in the central nervous system. The enduring effects of these drugs, even after long periods of abstinence, have suggested the presence of persistent molecular events caused by transcriptional, epigenetic and translational changes  within addiction-related brain regions , such as the nucleus accumbens (NAc), prefrontal cortex (PFC), and hippocampus. Notably, accumulating evidence has revealed an important role of epigenetic regulation in the effects of psychostimulants  (e.g., cocaine [4, 5]) and identified the involvement of modifications of histones present in chromatin [6, 7], DNA methylation [8, 9], and DNA hydroxymethylation in the development of addiction. With respect to METH, however, only a few studies have investigated the epigenetic effects of these drugs, mainly focusing on the NAc and striatum [10,11,12,13,14].
In animals, repeated exposure to METH produces behavioral sensitization, which is manifested as a progressive increase in locomotor responses that is elicited by a challenge dose of the drug after long-term withdrawal. Behavioral sensitization has been suggested to be a useful animal model of drug abuse that may reflect the properties of craving . Alterations to the associated brain regions, such the PFC, are suggested to mediate the etiology and maintenance of METH-induced behavioral sensitization . Repeated exposure to drugs such as cocaine or METH was also reported to produce long-lasting changes in synaptic plasticity in the ventral tegmental area (VTA), NAc, and PFC [17,18,19]. A study in nonhuman primates has shown that changes in the PFC following METH administration may persist for as long as 3.5 years .
In the present study, we focused on whether the changes in gene expression occurring in response to METH-induced behavioral sensitization are related to histone acetylation in the PFC. Given the importance of the PFC in the development of METH abuse, we screened the alterations in gene expression and histone acetylation in the PFC following METH-induced behavioral sensitization and further investigated the contribution of histone acetylation to the expression of the candidate genes in the different phases of the behavior model. Alterations in histone acetylation might, in part, affect the METH-induced changes in gene expression in the PFC.
Materials and methods
Adult male Sprague-Dawley rats (Grade SPF, Beijing Animal Center, Beijing, China) weighing 160–180 g were used in this study. Animals were maintained on a 12/12-hour light/dark cycle (lights on at 6:30 and off at 18:30), with an ambient temperature set at 20 °C to 22 °C. Animals were housed in 470 × 300 × 200-mm cages and provided access to rat maintenance food and water ad libitum. During housing, animals were monitored twice daily, and no abnormal events regarding the health status were detected. All experiments using animals were conducted in accordance with the NIH Guide for the Care and Use of Laboratory Animals (8th edition) and the guidelines of the National Research Council ( 398, China). All sections of this report adhere to the ARRIVE Guidelines for reporting animal research. A completed ARRIVE guidelines checklist is included.
METH-HCl was provided by the Beijing Public Security Bureau Forensic Medical Examination Center (China). The drug was dissolved in saline (0.9% NaCl).
The locomotor activity was measured in chambers (40 cm×40 cm×65 cm) using Digbehv spontaneous activity monitors  (Shanghai Jiliang Software Technology Co. Ltd., Shanghai, China). The total distance of horizontal locomotor activity was recorded with a video camera placed above the chamber and analyzed with the Digbehv software SuperState v3.0 system (AniLab Software & Instruments Co., Ningbo, China).
The behavioral sensitization procedures were conducted using the method described in our previous study . Behavioral sensitization consists of three phases (development, withdrawal and challenge phases) (Fig. 1a). The rats were administered METH (5 mg/kg, s.c.) or normal saline (NS) (1 ml/kg, s.c.) once daily in the acute phase or for 7 days in the development phase. On days 8–14, the rats were housed in home cages and did not receive injections in the withdrawal phase. On day 15, METH (1 mg/kg, s.c.) was administered to challenge the rats that were administered NS or METH in the development phase, and then the locomotor activity was recorded for 60 min in the challenge phase. The experiment was conducted using 40 rats (n =8 per group, 5 groups). A randomized METH and saline administration schedule was created using the standard = RAND() function in Microsoft Excel to reduce the number of animals used and the experimental error.
RNA isolation and cDNA synthesis
Twenty-four hours after the behavioral experiment, rats were anesthetized with pentobarbital (70 mg/kg, i.p.) and then euthanized by decapitation. PFC tissues were dissected (+ 5.16 mm to 3.00 mm from the bregma) according to the stereotactic atlas published by Paxinos and Watson. Frozen PFC tissues were thawed in TRIzol (Invitrogen, Cat. No. 15,596,018) and then processed according to the manufacturer’s protocol. Total RNA was reverse-transcribed using the ReverTra Ace qPCR RT Kit (Qiagen, Cat. No. 205,313).
Real-time PCR was performed using a SYBR Green PCR kit (Qiagen, Cat. No. 204,145) and quantified using the 2−ΔΔCt method. Real-time PCR was then performed using an ABI 7300 Real-Time PCR System (Applied Biosystems, Foster City, CA), amplified using SYBR Green (TOYOPO, Cat. No. QPS-201) and quantified using the 2−ΔΔCt method. The primers for the analyzed genes are shown in Additional file 1: Table S1. Each sample was assayed in duplicate and normalized to the GAPDH levels.
PFC tissues were crosslinked with 1% formaldehyde at 37 °C for 10 min, and then glycine was added to a final concentration of 0.125 M and incubated for 5 min to stop the crosslinking. Samples were centrifuged at 800g at 4 °C for 10 min to pellet tissues. The PFC tissues were sonicated on wet ice in 300 µl of nuclear lysis buffer using a sonicator (Sonics Vibra-cell™, Sonics & Materials Inc., Newtown, Connecticut, USA). The sonication conditions were used: amplitude 70%, pulse 10 sec, and pause 59 sec for 9 cycles. The supernatant was collected as described in the instructions of a chromatin immunoprecipitation kit (Millipore, Cat. No. 17408). The lysate was precipitated with antibodies against acetyl-histone H3 (Millipore, Cat. No. 06599) and acetyl-histone H4 (Millipore, Cat. No. 06866). The following primer sequences were used for ChIP-PCR: POU3F2: forward primer 5′-CAGTACCGCTGAATAACGC-3′ and reverse primer 5′-CTGAGCCTGCTTACGATGA-3′; ANP32A: forward primer 5′-CCATCGGAGCATGTTTCGT-3′ and reverse primer 5′-AGGGTGAGATCGCCGTCTT-3′. The PFC samples used for qPCR and ChIP assay were obtained from the rats in another behavioral experiment, which were different from the samples used in microarray experiments. The samples were processed separately to perform the qPCR and ChIP experiments.
Twenty-four hours after the behavioral experiment, the samples extracted from the PFCs of 8 rats per group were mixed and processed for the mRNA microarray or ChIP on chip DNA microarray. The mRNA microarray analysis was conducted using the Rat Whole Genome One Array that examines 19,803 unique genes using 24,358 probes (Phalanx Biotech Group, Hsinchu, Taiwan). The fluorescently amplified RNA targets were prepared from 2.5 µg of total RNA samples using the OneArray® Amino Allyl aRNA Amplification Kit (Phalanx Biotech Group), labeled with Cy5 dye (Amersham Pharmacia, Piscataway, NJ, USA), and then hybridized to the Rat Whole Genome OneArray® (ROA 1.1) (Phalanx Biotech Group) according to the Phalanx hybridization protocol. The slides were scanned using an Axon 4000 scanner (Molecular Devices, Sunnyvale, CA, USA). The Cy5 fluorescence intensity of each spot was analyzed using GenePix pro software (Molecular Devices). The signal intensity of each spot was loaded into the Rosetta Resolver System® (Rosetta Biosoftware) for data analysis.
For ChIP on chip, the DNA microarray analysis was performed using the Rat 3 × 720K RefSeq Promoter Array that examines 15,287 unique gene promoters overlapping 1,070 − 4,280 bp (Nimblegen). The input DNA was amplified using the WGA2 kit (Sigma, Cat. No. WGA2-50R) and purified with the QIAquick PCR purification kit (Qiagen). The product was labeled with Cy3 and Cy5 dyes and hybridized using a Roche NimbleGen hybridization kit. After drying, the slides were scanned with a Roche NimbleGen MS200, and the data were extracted using NimbleGen2.6 software. For each probe on the array, the log2 ratios of the Cy5-labeled ChIP DNA/Cy3-labeled input DNA were calculated using Signal MAP software. Peak detection was performed using NimbleScan software.
Differential gene expression analysis and functional enrichment
Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, were performed to identify which differentially expressed genes were significantly enriched in GO terms and metabolic pathways at the thresholds of a corrected P-value [also named false discovery rate (FDR)] ≤ 0.05 and log2 fold change (log2 FC). GO functional enrichment and KEGG pathway analyses were conducted with the KOBAS program (http://kobas.cbi.pku.edu.cn/home.do) .
All data are presented as means ± SEM. The challenge data were assessed using one-way analysis of variance (ANOVA) followed by a post hoc test. For other data, a t-test was used. P values < 0.05 were considered statistically significant.
METH-induced behavioral sensitization
The model of behavioral sensitization was established by 7 days of METH (5.0 mg/kg, s.c.) administration followed by 7 days of withdrawal and challenged with a single METH injection (1.0 mg/kg, s.c.) on day 15. One-way ANOVA showed a difference between the groups (F(2, 21) = 14.58, P < 0.001). Rats in the METH sensitization group (METH-METH group) showed significantly greater locomotor activity than either the saline control (NS-NS) or METH acute treatment group (NS-METH) (Tukey’s post hoc test, P < 0.001 compared with the NS-NS group; P < 0.01 compared with the NS-METH group, Fig. 1b), indicating the increased locomotor response of sensitization. Compared with the NS-NS group, rats in the NS-METH group showed longer distances (P < 0.05, Fig. 1b).
Changes in mRNA expression in the PFC of METH-induced behavioral sensitization rat
An mRNA microarray was performed to determine a genome-wide profile of METH-induced behavioral sensitization. The cut-off value was set at a fold change of ≤ 0.5 or ≥ 2.0 (P < 0.01). Figure 2 shows the clustering diagram and a heat map of the changes in mRNA expression in the NS-NS, NS-METH, and METH-METH groups. The expression of a large number of genes was altered by METH exposure (P < 0.01) (Fig. 2b). Among these, 275 genes were differentially expressed in the METH-METH group compared with the NS-METH group (97 genes were upregulated and 178 genes were downregulated), and 232 genes displayed differential expression in the NS-METH group compared with the NS-NS group (122 genes were upregulated and 110 genes were downregulated) (Additional file 2: Table S2).
Ontology and pathway analyses of differentially expressed mRNAs in the PFC of METH-induced sensitization rats
We performed a GO analysis and categorized the differentially expressed genes to investigate whether the clustering of differentially expressed genes induced by METH exposure correlates with functional categories. Figure 3a and Additional file 3: Table S3 show the top 10 ranked GO terms and the number of genes in each category. Compared with the acute METH treatment (NS-METH), the biological processes that were markedly increased in the METH-induced behavioral sensitization (METH-METH) group included biological regulation, cellular process and organism process. Regarding the molecular function, binding and protein binding were enriched in the METH-METH group.
We further performed a KEGG pathway enrichment analysis to assess the functional features of METH-induced behavioral sensitization-mediated alterations in gene sets. The top 10 significantly altered pathways are shown in Fig. 3b and Additional file 4: Table S4. Many signal transduction pathways were enriched in the METH-induced behavioral sensitization group, including cell adhesion, PI3K-AKT signaling pathway, olfactory transduction, and cell cycle (METH-METH group vs. NS-METH group, P < 0.0001). For the METH acute treatment (NS-METH) group, the genes categorized by the KEGG analysis were enriched in the cell adhesion, the pathway in cancer, and dopaminergic pathway.
Changes in histone acetylation in the PFC of METH-induced sensitization rats
ChIP coupled with a DNA microarray analysis revealed that METH increased histone acetylation (H3 or H4) on a large number of gene promoters (Fig. 4), in accordance with mRNA activation. Compared with the NS-METH group, the METH-METH group induced many more acetylation modifications of H3 than H4 on the promoters; specifically, 821 genes presented H3 hyperacetylation and 10 genes showed H4 hyperacetylation (Additional file 5: Table S5). The NS-METH treatment also induced hyperacetylation on the gene promoters, including 947 genes with H3 hyperacetylation and 4902 genes with H4 hyperacetylation. Very little H3/H4 hyperacetylation overlapped between the METH-METH and NS-METH groups.
Real‐time PCR confirmation of the candidate genes
Based on the mRNA microarray data, GO/KEGG enrichment analyses, and the analysis of the functional properties, the genes listed in Additional file 6: Table S6 were likely to be associated with METH-induced behavioral sensitization. Of these genes, AVP, BCL2L1, EGR1, E2F3, LNX2, SHOC2, STX2 and ZFP36 clustered into the top 10 GO and KEGG categories, and ANP32A, EML2, EXOG, HIRA, METRN, POU3F2, STK32, SYT8, and TRIM17 were included in the GO category. According to the analysis of the GO categories, the expression of LNX2, SHOC2 and STX2 was altered by a single METH injection, and the expression of the other genes listed above was affected by METH-induced behavioral sensitization (METH-METH) but not the acute METH treatment (NS-METH). Furthermore, all of the genes displayed H3/H4 hyperacetylation. We then used qPCR to confirm the expression of some genes of interest. METH challenge in the behavioral sensitization model caused a significant increase in the mRNA levels of acidic nuclear phosphoprotein 32 family member A (ANP32A), calcium/calmodulin-dependent protein kinase II inhibitor 1 (CAMK2N1), echinoderm microtubule-associated protein-like 2 (EML2) and POU class 3 homeobox 2 (POU3F2), and a decrease in syntaxin 2 (STX2), tripartite motif-containing 17 (TRIM17), and zinc finger protein 36 (ZFP36) (P < 0.05, one-way ANOVA followed by Turkey’s post hoc test, Figs. 5 and 6).
Alterations in ANP32A and POU3F2 expression in the development, withdrawal and challenge periods of METH-induced sensitization
Then, we selected ANP32A and POU3F2 and further measured the mRNA expression and histone modification of these two genes. In addition to the challenge of behavioral sensitization, the development and withdrawal phases were also involved to investigate the changes in ANP32A and POU3F2 expression throughout the whole process of METH-induced behavioral sensitization. After a chronic METH (5 mg/kg, s.c.) treatment for 7 days (development phase), the levels of the ANP32A mRNA and H4 acetylation were markedly increased compared with the NS group (P < 0.05, t-test, Fig. 6a, b). However, the expression of the ANP32A mRNA returned to the normal level, while H4 acetylation remained at a high level after 7 days of withdrawal (P < 0.05, t-test). Then, an injection of METH (1 mg/kg, s.c.) was administered as a challenge on day 15. One-way ANOVA showed a significant difference in mRNA expression (F(2, 18) = 8.824, P < 0.01) and H4 acetylation (F(2, 26) = 6.072, P < 0.01) between groups, and Tukey’s post hoc test showed that the challenge induced the expression of ANP32A mRNA (P < 0.01, METH-METH vs. NS-NS) and H4 acetylation (P < 0.01, METH-METH vs. NS-METH, P < 0.05, METH-METH vs. NS-NS) in the METH pretreatment group (METH-METH). The injection of METH alone also increased the level of the ANP32A mRNA (P < 0.05, NS-METH vs. NS-NS), but did not affect histone acetylation.
We next measured the change in levels of the POU3F2 mRNA and histone modification. The expression of the POU3F2 mRNA was significantly increased by the chronic METH treatment for 7 days and decreased after 7 days of withdrawal (P < 0.05, t-test), while H3 and H4 acetylation were not affected by the treatment (Fig. 6c, d). In the challenge phase, one-way ANOVA showed that challenge with a low dose of METH caused a significant difference in mRNA expression (F(2, 45) = 4.01, P < 0.05) and H3 and H4 hyperacetylation (H3: F(2, 36) = 4.17, P < 0.05; H4: F(2, 38) = 3.47, P < 0.05) between groups. Tukey’s post hoc test revealed a significant increase in POU3F2 mRNA expression (P < 0.05, METH-METH vs. NS-NS) and H3 and H4 hyperacetylation in the behavioral sensitization group (H3: P < 0.05, METH-METH vs. NS-NS, METH-METH vs. NS-METH; H4: P < 0.05, METH-METH vs. NS-NS). Notably, the single METH injection (NS-METH) also increased the expression of the POU3F2 mRNA (P < 0.05, vs. NS-NS), but not the levels of H3/H4 acetylation.
Behavioral sensitization is an animal model that reflects some aspects of drug addiction, such as drug seeking and relapse behaviors, but the neurobiological mechanism has not yet been completely elucidated. Drug-induced changes in gene expression and sustained changes in genes regulated by epigenetic mechanisms in addiction-associated regions are thought to contribute to long-lasting behavioral abnormalities and neuroadaptation. The present study examined the transcriptional changes in gene expression and alterations in histone acetylation in the PFC of adult rats that experienced the different phases of METH-induced behavioral sensitization to further analyze the neurological mechanism underlying the development of METH addiction. The mRNA microarray and ChIP microarray analyses revealed 275 differentially expressed genes and significant hyperacetylation, including 821 genes with H3 hyperacetylation and 10 genes with H4 hyperacetylation, caused by METH-induced behavioral sensitization in the PFC, suggesting that histone hyperacetylation is associated with the increased expression of a set of genes caused by METH-induced behavioral sensitization.
Based on accumulating evidence, drug exposure causes modifications of global histone acetylation or specific changes in the gene-promoter regions in a complex manner. As shown in the studies by Wang et al. chronic self- or experimenter-administered cocaine induced long-term increases in global H3 and H4 acetylation in the NAc shell , and cocaine abstinence caused the acetylation of histone H3K9 and histone H4K8 in the PFC . ChIP sequencing provides access to precise patterns of histone acetylation. According to Renthal et al. , chronic cocaine exposure produces a larger number of genes with H3 acetylation than H4 acetylation in the NAc. A set of genes, such as CBP in the NAc shell  or BDNF in the PFC , is activated by cocaine, and these changes are associated with the hyperacetylation of H3 and/or H4. Regarding METH exposure, Tracey et al.. reported time-dependent increases in the levels of acetylated H4K5 and H4K8 and decreases in the levels of acetylated H3K9, H3K18, and H4K16 in the NAc following a single METH injection . González et al.. observed increased levels of total H3 acetylation and decreased H4 acetylation after a single METH injection, while repeated METH injections decreased levels of total H3 and H4 acetylation in the mPFC [28, 29]. As shown in the present study, METH-induced behavioral sensitization caused a large number of genes to exhibit H3 hyperacetylation and few genes exhibited H4 hyperacetylation, while more H4 hyperacetylation than H3 hyperacetylation was induced by acute METH exposure in the PFC.
METH-induced behavioral sensitization led to a change in the transcriptional profile of the PFC in the present study, as manifested by the large number of differentially expressed genes categorized by GO and KEGG pathway enrichment, including cell adhesion, the PI3K-AKT signaling pathway, olfactory transduction, and the cell cycle, compared with the acute METH treatment control. In fact, exposure to drugs, such as cocaine and morphine, has consistently been shown to induce cell adhesion  and activate intracellular signaling pathways, such as PI3K-AKT  and the prolactin signaling pathway , to mediate the effects of addiction. Some potential effects of the genes identified and confirmed in this study using qPCR on behavior have been identified in the previous reports. The expression of some genes encoding proteins involved in behavioral sensitization decreased (e.g., STX2, ZFP36 and TRIM17) and the expression of some genes increased (e.g., ANP32A, CAMK2N1, EML2 and POU3F2) at the mRNA level induced by METH challenge. These genes encode not only transcription factors but also proteins involved in neuron differentiation, synaptic plasticity and intracellular signaling. STX2 was reported to play a role in CNS damage [33,34,35]. STX2 interacts with neurons and promotes structural changes in the cell membrane that may indicate the synaptic vesicle cycle . STX2 induces the PI3K-involved signaling pathway to release glutamate, which contributes to the CNS impairment . EML2 encodes echinoderm microtubule-associated protein-like 2 protein, which plays a role in the formation of the mitotic spindle and interphase microtubule network . The EML2 methylation status differs among individual cell types in the brain tissue , suggesting that EML2 is an important determinant of cell type-specific differentiation in the brain .
We further identified the changes in ANP32A and POU3F2 by conducting qPCR and H3/H4 acetylation measurements. ANP32A is mostly localized in the nucleus but is also present on the endoplasmic reticulum and in the cytoplasm . To date, ANP32A has been reported to play roles in a variety of cellular functions, including RNA transport, transcription, cell-mediated cytotoxicity, and tumor suppression [42,43,44]. In the CNS, ANP32A contributes to the survival and differentiation of neurons [45,46,47], as well as synaptic plasticity and memory in aged mice or tau transgenic mice used as a model of Alzheimer’s disease. In the present study, levels of the ANP32A mRNA were increased by acute or chronic treatment with METH, as well as METH challenge-induced behavioral sensitization. The ChIP assay showed a robust induction of the acetylation of H4 at the promoter of ANP32A caused by chronic METH use or challenge-induced behavioral sensitization, but not by acute METH treatment. The lack of change in H4 hyperacetylation in response to acute METH exposure but its potential increase in the development and challenge stages of METH-induced behavioral sensitization suggested that H4 hyperacetylation at the ANP32A gene may have an important role in modulating the motivation for METH reinforcement.
POU3F2 (also known as BRN2), a class III POU protein , is prominently expressed in the neocortex and upregulated in the progenitor cells of the subventricular zone, intermediate zone and outer layer of the neocortex during the early stages of embryonic brain development [49, 50]. POU3F2 contributes to neural formation and cell fate determination and regulates cortical neural migration, neurogenesis and the positioning of cortical neurons [49, 51]. Furthermore, POU3F2 has been suggested to regulate tyrosine hydroxylase and tryptophan hydroxylase 2, the rate-limiting enzymes involved in the synthesis of dopamine and 5-HT , both of which are important for the regulation of psychological and physiological functions. POU3F2Δ/Δ mice in which all three homopolymeric amino acid repeats were deleted from the POU3F2 transactivation domain displayed cognitive impairments in object recognition and object location tests , suggesting that POU3F2 is involved in cognitive function. In the present study, the expression of the POU3F2 mRNA was upregulated by METH exposure, including acute or chronic METH treatment, and METH challenge-induced behavioral sensitization. However, when METH was withdrawn after chronic treatment, the POU3F2 mRNA was downregulated instead. Notably, histone acetylation at the promoter of POU3F2 showed a specific change during the course of METH treatment, which manifested as both H4 and H3 hyperacetylation only in response to METH-induced behavioral sensitization but not acute or chronic METH treatment or withdrawal. The change in the patterns suggests that histone acetylation of POU3F2 is involved in the mechanism of neuroplasticity underlying the transition from recreational METH use to compulsive use. Moreover, POU3F2 is involved in schizophrenia [54, 55], the syndrome and pathology of which are very similar to METH addiction. Therefore, we speculated that POU3F2 potentially represents one of the common genetic bases underlying schizophrenia and METH use.
We identified the changes in gene expression and the diverse alterations in histone acetylation in the PFC in the different phases of METH-induced behavioral sensitization. By further analyzing the levels of acetylated histones bound to the promoters of some genes (e.g., ANP32A and POU3F2), we revealed that histone modifications might participate in the transcriptional changes in gene expression caused by the specific manners of METH treatment, which further confirmed the epigenetic mechanisms of gene regulation underlying the long-lasting changes induced by drugs. Meanwhile, METH-induced histone acetylation might be necessary, but not sufficient, to induce changes in gene expression. In the present study, the functions of the genes of interest in the brain were not explored. The transcript levels and histone acetylation of the genes associated with functions in specific brain regions may provide new knowledge that will enable researchers to elucidate the neurobiological mechanism underlying METH addiction and develop alternative pharmacological methods for treating psychiatric disorders, which should be extensively studied in the near future.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files. The sequencing data are available from Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE166096).
Kyoto Encyclopedia of Genes and Genomes
Acidic leucine-rich nuclear phosphoprotein-32A
POU domain, class 3, transcription factor 2
Robison AJ, Nestler EJ. Transcriptional and epigenetic mechanisms of addiction. Nat Rev Neurosci. 2011;12(11):623–37.
Godino A, Jayanthi S, Cadet JL. Epigenetic landscape of amphetamine and methamphetamine addiction in rodents. Epigenetics. 2015;10(7):574–80.
Nestler EJ. Epigenetic mechanisms of drug addiction. Neuropharmacology. 2014;76 Pt B:259–68.
Pierce RC, Fant B, Swinford-Jackson SE, Heller EA, Berrettini WH, Wimmer ME. Environmental, genetic and epigenetic contributions to cocaine addiction. Neuropsychopharmacology. 2018;43(7):1471–80.
Heller EA, Hamilton PJ, Burek DD, Lombroso SI, Pena CJ, Neve RL, Nestler EJ. Targeted epigenetic remodeling of the cdk5 gene in nucleus accumbens regulates cocaine- and stress-evoked behavior. The Journal of neuroscience: the official journal of the Society for Neuroscience. 2016;36(17):4690–7.
Rogge GA, Wood MA. The role of histone acetylation in cocaine-induced neural plasticity and behavior. Neuropsychopharmacology. 2013;38(1):94–110.
Kennedy PJ, Harvey E. Histone deacetylases as potential targets for cocaine addiction. CNS Neurol Disord Drug Target. 2015;14(6):764–72.
Urb M, Niinep K, Matsalu T, Kipper K, Herodes K, Zharkovsky A, Timmusk T, Anier K, Kalda A. The role of DNA methyltransferase activity in cocaine treatment and withdrawal in the nucleus accumbens of mice. Addict Biol. 2020;25(1):e12720.
Le Q, Yan B, Yu X, Li Y, Song H, Zhu H, Hou W, Ma D, Wu F, Zhou Y, et al. Drug-seeking motivation level in male rats determines offspring susceptibility or resistance to cocaine-seeking behaviour. Nat Commun. 2017;8:15527.
Cadet JL, Brannock C, Krasnova IN, Jayanthi S, Ladenheim B, McCoy MT, Walther D, Godino A, Pirooznia M, Lee RS. Genome-wide DNA hydroxymethylation identifies potassium channels in the nucleus accumbens as discriminators of methamphetamine addiction and abstinence. Mol Psychiatry. 2017;22(8):1196–204.
Shibasaki M, Mizuno K, Kurokawa K, Ohkuma S. L-type voltage-dependent calcium channels facilitate acetylation of histone H3 through PKCgamma phosphorylation in mice with methamphetamine-induced place preference. J Neurochem. 2011;118(6):1056–66.
Cadet JL, Jayanthi S, McCoy MT, Ladenheim B, Saint-Preux F, Lehrmann E, De S, Becker KG, Brannock C. Genome-wide profiling identifies a subset of methamphetamine (METH)-induced genes associated with METH-induced increased H4K5Ac binding in the rat striatum. BMC Genomics. 2013;14:545.
Martin TA, Jayanthi S, McCoy MT, Brannock C, Ladenheim B, Garrett T, Lehrmann E, Becker KG, Cadet JL. Methamphetamine causes differential alterations in gene expression and patterns of histone acetylation/hypoacetylation in the rat nucleus accumbens. PLoS One. 2012;7(3):e34236.
Li H, Li F, Wu N, Su RB, Li J. Methamphetamine induces dynamic changes of histone deacetylases in different phases of behavioral sensitization. CNS Neurosci Ther. 2014;20(9):874–6.
Steketee JD, Kalivas PW. Drug wanting: behavioral sensitization and relapse to drug-seeking behavior. Pharmacol Rev. 2011;63(2):348–65.
Gulley JM, Stanis JJ. Adaptations in medial prefrontal cortex function associated with amphetamine-induced behavioral sensitization. Neuroscience. 2010;166(2):615–24.
Kourrich S, Rothwell PE, Klug JR, Thomas MJ. Cocaine experience controls bidirectional synaptic plasticity in the nucleus accumbens. J Neurosci. 2007;27(30):7921–8.
Dong Y, Nasif FJ, Tsui JJ, Ju WY, Cooper DC, Hu XT, Malenka RC, White FJ. Cocaine-induced plasticity of intrinsic membrane properties in prefrontal cortex pyramidal neurons: adaptations in potassium currents. J Neurosci. 2005;25(4):936–40.
Onn SP, Grace AA. Amphetamine withdrawal alters bistable states and cellular coupling in rat prefrontal cortex and nucleus accumbens neurons recorded in vivo. J Neurosci. 2000;20(6):2332–45.
Selemon LD, Begovic A, Goldman-Rakic PS, Castner SA. Amphetamine sensitization alters dendritic morphology in prefrontal cortical pyramidal neurons in the non-human primate. Neuropsychopharmacology. 2007;32(4):919–31.
Xu CM, Wang J, Wu P, Xue YX, Zhu WL, Li QQ, Zhai HF, Shi J, Lu L. Glycogen synthase kinase 3beta in the nucleus accumbens core is critical for methamphetamine-induced behavioral sensitization. J Neurochem. 2011;118(1):126–39.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316-322.
Wang L, Lv Z, Hu Z, Sheng J, Hui B, Sun J, Ma L. Chronic cocaine-induced H3 acetylation and transcriptional activation of CaMKIIalpha in the nucleus accumbens is critical for motivation for drug reinforcement. Neuropsychopharmacology. 2010;35(4):913–28.
Sadakierska-Chudy A, Frankowska M, Jastrzebska J, Wydra K, Miszkiel J, Sanak M, Filip M. Cocaine administration and its withdrawal enhance the expression of genes encoding histone-modifying enzymes and histone acetylation in the rat prefrontal cortex. Neurotox Res. 2017;32(1):141–50.
Renthal W, Kumar A, Xiao G, Wilkinson M, Covington HE 3rd, Maze I, Sikder D, Robison AJ, LaPlant Q, Dietz DM, et al. Genome-wide analysis of chromatin regulation by cocaine reveals a role for sirtuins. Neuron. 2009;62(3):335–48.
Malvaez M, Mhillaj E, Matheos DP, Palmery M, Wood MA. CBP in the nucleus accumbens regulates cocaine-induced histone acetylation and is critical for cocaine-associated behaviors. The Journal of neuroscience: the official journal of the Society for Neuroscience. 2011;31(47):16941–8.
Sadri-Vakili G, Kumaresan V, Schmidt HD, Famous KR, Chawla P, Vassoler FM, Overland RP, Xia E, Bass CE, Terwilliger EF, et al. Cocaine-induced chromatin remodeling increases brain-derived neurotrophic factor transcription in the rat medial prefrontal cortex, which alters the reinforcing efficacy of cocaine. The Journal of neuroscience: the official journal of the Society for Neuroscience. 2010;30(35):11735–44.
Gonzalez B, Torres OV, Jayanthi S, Gomez N, Sosa MH, Bernardi A, Urbano FJ, Garcia-Rill E, Cadet JL, Bisagno V. The effects of single-dose injections of modafinil and methamphetamine on epigenetic and functional markers in the mouse medial prefrontal cortex: potential role of dopamine receptors. Prog Neuro-psychopharmacol Biol Psychiatry. 2019;88:222–34.
Gonzalez B, Jayanthi S, Gomez N, Torres OV, Sosa MH, Bernardi A, Urbano FJ, Garcia-Rill E, Cadet JL, Bisagno V. Repeated methamphetamine and modafinil induce differential cognitive effects and specific histone acetylation and DNA methylation profiles in the mouse medial prefrontal cortex. Prog Neuro-psychopharmacol Biol Psychiatry. 2018;82:1–11.
Zhong X, Drgonova J, Li CY, Uhl GR. Human cell adhesion molecules: annotated functional subtypes and overrepresentation of addiction-associated genes. Ann N Y Acad Sci. 2015;1349:83–95.
Cui Y, Zhang XQ, Xin WJ, Jing J, Liu XG. Activation of phosphatidylinositol 3-kinase/Akt-mammalian target of Rapamycin signaling pathway in the hippocampus is essential for the acquisition of morphine-induced place preference in rats. Neuroscience. 2010;171(1):134–43.
Laurent P, Becker JA, Valverde O, Ledent C, de Kerchove d’Exaerde A, Schiffmann SN, Maldonado R, Vassart G, Parmentier M. The prolactin-releasing peptide antagonizes the opioid system through its receptor GPR10. Nat Neurosci. 2005;8(12):1735–41.
Kawalia SB, Raschka T, Naz M, de Matos Simoes R, Senger P, Hofmann-Apitius M. Analytical strategy to prioritize Alzheimer’s disease candidate genes in gene regulatory networks using public expression data. JAD. 2017;59(4):1237–54.
Kramer J, Deppe M, Gobel K, Tabelow K, Wiendl H, Meuth SG. Recovery of thalamic microstructural damage after Shiga toxin 2-associated hemolytic-uremic syndrome. J Neurol Sci. 2015;356(1–2):175–83.
Obata F, Tohyama K, Bonev AD, Kolling GL, Keepers TR, Gross LK, Nelson MT, Sato S, Obrig TG. Shiga toxin 2 affects the central nervous system through receptor globotriaosylceramide localized to neurons. J Infect Dis. 2008;198(9):1398–406.
Sudhof TC, Petrenko AG, Whittaker VP, Jahn R. Molecular approaches to synaptic vesicle exocytosis. Prog Brain Res. 1993;98:235–40.
Obata F, Hippler LM, Saha P, Jandhyala DM, Latinovic OS. Shiga toxin type-2 (STX2) induces glutamate release via phosphoinositide 3-kinase (PI3K) pathway in murine neurons. Front Mol Neurosci. 2015;8:30.
Richards MW, O’Regan L, Roth D, Montgomery JM, Straube A, Fry AM, Bayliss R. Microtubule association of EML proteins and the EML4-ALK variant 3 oncoprotein require an N-terminal trimerization domain. Biochem J. 2015;467(3):529–36.
Duong CV, Emes RD, Wessely F, Yacqub-Usman K, Clayton RN, Farrell WE. Quantitative, genome-wide analysis of the DNA methylome in sporadic pituitary adenomas. Endocrine-related Cancer. 2012;19(6):805–16.
Samsuwan J, Muangsub T, Yanatatsaneejit P, Mutirangura A, Kitkumthorn N. Combined bisulfite restriction analysis for brain tissue identification. Forensic Sci Int. 2018;286:42–5.
Wang S, Wang Y, Lu Q, Liu X, Wang F, Ma X, Cui C, Shi C, Li J, Zhang D. The expression and distributions of ANP32A in the developing brain. BioMed Res Int. 2015;2015:207347.
Bai J, Brody JR, Kadkol SS, Pasternack GR. Tumor suppression and potentiation by manipulation of pp32 expression. Oncogene. 2001, 20(17):2153–60.
Adegbola O, Pasternack GR. Phosphorylated retinoblastoma protein complexes with pp32 and inhibits pp32-mediated apoptosis. J Biol Chem. 2005;280(16):15497–502.
Schafer ZT, Parrish AB, Wright KM, Margolis SS, Marks JR, Deshmukh M, Kornbluth S. Enhanced sensitivity to cytochrome c-induced apoptosis mediated by PHAPI in breast cancer cells. Cancer research. 2006;66(4):2210–8.
Khan MZ, Vaidya A, Meucci O. CXCL12-mediated regulation of ANP32A/Lanp, a component of the inhibitor of histone acetyl transferase (INHAT) complex, in cortical neurons. J Neuroimmune Pharmacol. 2011;6(1):163–70.
Kular RK, Cvetanovic M, Siferd S, Kini AR, Opal P. Neuronal differentiation is regulated by leucine-rich acidic nuclear protein (LANP), a member of the inhibitor of histone acetyltransferase complex. J Biol Chem. 2009;284(12):7783–92.
Opal P, Garcia JJ, McCall AE, Xu B, Weeber EJ, Sweatt JD, Orr HT, Zoghbi HY. Generation and characterization of LANP/pp32 null mice. Mol Cell Biol. 2004;24(8):3140–9.
Verrijzer CP, Van der Vliet PC. POU domain transcription factors. Biochim Biophys Acta. 1993;1173(1):1–21.
Fujii H, Hamada H. A CNS-specific POU transcription factor, Brn-2, is required for establishing mammalian neural cell lineages. Neuron. 1993;11(6):1197–206.
Hauptmann G, Gerster T. Combinatorial expression of zebrafish Brn-1- and Brn-2-related POU genes in the embryonic brain, pronephric primordium, and pharyngeal arches. Dev Dyn. 2000;218(2):345–58.
Wegner M, Drolet DW, Rosenfeld MG. POU-domain proteins: structure and function of developmental regulators. Curr Opin Cell Biol. 1993;5(3):488–98.
Scheuch K, Lautenschlager M, Grohmann M, Stahlberg S, Kirchheiner J, Zill P, Heinz A, Walther DJ, Priller J. Characterization of a functional promoter polymorphism of the human tryptophan hydroxylase 2 gene in serotonergic raphe neurons. Biol Psychiatry. 2007;62(11):1288–94.
Hashizume K, Yamanaka M, Ueda S. POU3F2 participates in cognitive function and adult hippocampal neurogenesis via mammalian-characteristic amino acid repeats. Genes Brain Behav. 2018;17(2):118–25.
Prata DP, Costa-Neves B, Cosme G, Vassos E. Unravelling the genetic basis of schizophrenia and bipolar disorder with GWAS: A systematic review. J Psychiatr Res. 2019;114:178–207.
Pearl JR, Colantuoni C, Bergey DE, Funk CC, Shannon P, Basu B, Casella AM, Oshone RT, Hood L, Price ND, et al. Genome-scale transcriptional regulatory network models of psychiatric and neurodegenerative disorders. Cell Syst. 2019;8(2):122–35 e127.
We thank Dr. Su-Wei Zhang for providing technical assistance in analyzing the data.
This study was supported by the National Key R&D Program of China (No. 2017YFC1310404) and the National Natural Science Foundation of China (Nos. 81673409 and U1502225). The grants supported the necessary funds for acquiring the materials and animals that used in the research, the salaries for laboratory and veterinary personnel, and support staff for data collection and analysis. The sponsors were not involved in the design, methods, data collection, or preparation of this paper.
Ethics approval and consent to participate
All animal experimental procedures and care of animals were performed in accordance with the NIH Guide for the Care and Use of Laboratory Animals (8th edition) and the guidelines of the National Research Council ( 398, China). Ethical approval was obtained from Beijing Institute of Pharmacology and Toxicology.
Consent for publication
The authors have no competing interests to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The primers for the genes in Real-time PCR.
Raw data of mRNA microarray in METH-induced behavioral sensitization.
The detailed information of top 10 GO terms and the categorized differential expression genes in the group of METH-METH vs. NS-METH, METH-METH vs. NS-NS, and NS-METH vs. NS-NS.
The detailed information of top 10 KEGG pathway enrichment and the categorized differential expression genes in the group of METH-METH vs. NS-METH, METH-METH vs. NS-NS, and NS-METH vs. NS-NS.
Raw data of ChIP-All peaks in METH-induced behavioral sensitization.
The information of genes of interest. The genes were selected based on the analysis of mRNA microarray, GO and KEGG enrichment.
About this article
Cite this article
Li, H., Chen, JA., Ding, QZ. et al. Behavioral sensitization induced by methamphetamine causes differential alterations in gene expression and histone acetylation of the prefrontal cortex in rats. BMC Neurosci 22, 24 (2021). https://doi.org/10.1186/s12868-021-00616-5
- Behavioral sensitization
- Differentially expressed genes
- Histone acetylation