Research article | Open | Published:
Stress susceptibility-specific phenotype associated with different hippocampal transcriptomic responses to chronic tricyclic antidepressant treatment in mice
BMC Neurosciencevolume 14, Article number: 144 (2013)
The effects of chronic treatment with tricyclic antidepressant (desipramine, DMI) on the hippocampal transcriptome in mice displaying high and low swim stress-induced analgesia (HA and LA lines) were studied. These mice displayed different depression-like behavioral responses to DMI: stress-sensitive HA animals responded to DMI, while LA animals did not.
To investigate the effects of DMI treatment on gene expression profiling, whole-genome Illumina Expression BeadChip arrays and qPCR were used. Total RNA isolated from hippocampi was used. Expression profiling was then performed and data were analyzed bioinformatically to assess the influence of stress susceptibility-specific phenotypes on hippocampal transcriptomic responses to chronic DMI. DMI treatment affected the expression of 71 genes in HA mice and 41 genes in LA mice. We observed the upregulation of Igf2 and the genes involved in neurogenesis (HA: Sema3f, Ntng1, Gbx2, Efna5, and Rora; LA: Otx2, Rarb, and Drd1a) in both mouse lines. In HA mice, we observed the upregulation of genes involved in neurotransmitter transport, the termination of GABA and glycine activity (Slc6a11, Slc6a9), glutamate uptake (Slc17a6), and the downregulation of neuropeptide Y (Npy) and corticotropin releasing hormone-binding protein (Crhbp). In LA mice, we also observed the upregulation of other genes involved in neuroprotection (Ttr, Igfbp2, Prlr) and the downregulation of genes involved in calcium signaling and ion binding (Adcy1, Cckbr, Myl4, Slu7, Scrp1, Zfp330).
Several antidepressant treatment responses are similar in individuals with different sensitivities to stress, including the upregulation of Igf2 and the genes involved in neurogenesis. However, the findings also reveal that many responses to antidepressant treatments, involving the action of individual genes engaged in neurogenesis, neurotransmitter transport and neuroprotection, depend on constitutive hippocampal transcriptomic profiles and might be genotype dependent. The results suggest that, when and if this becomes feasible, antidepressant treatment should take into consideration individual sensitivity to stress.
The hippocampus is an integral component of the limbic system and mechanisms that control mood [1, 2]. The hippocampus also plays an important role in stress, pain and nociceptive responses and might mediate depressive behaviors . Because the propensity to develop depression and the responses to antidepressant treatments might depend on the susceptibility of an individual to stress, we studied mice selected for high (HA line) or low (LA line) stress reactivity, as measured by the magnitude of swim stress-induced analgesia (SIA) . There were differences among the mice in terms of responses to a variety of antidepressants (desipramine, venlafaxine, and aminosenktide) and profound differences in tail suspension and forced swim tests, reflecting depression-like behaviors [5–7]. Under control conditions, clear differences in the constitutive, hippocampal transcriptomic profiles related to the upregulation of genes involved in calcium signaling in LA mice and the robust upregulation of genes encoding receptors for neurotensin (Ntsr2) and GABA (Gabard) in HA mice were observed . These findings suggested that the selective breeding for swim SIA affected many aspects of hippocampal neuron physiology, including metabolism, structural changes and cellular signaling. Furthermore, unpredictable chronic mild stress (CMS), used to trigger desipramine-reversible depressive-like behavior , produced mouse line-specific responses in hippocampal gene expression . The results suggest that hippocampal responses to stress depend on genotype, and the potential drug targets against the detrimental effects of stress might include glutamate transport, neurogenesis, axonogenesis, neurite development, or chromatin modifications, to name a few.
The present study was conducted to identify transcriptomic response involved in the action of the tricyclic antidepressant, desipramine, in the hippocampus. We used expression microarray analyses to characterize changes in the transcription profiles and identify genes that are affected after chronic treatment with DMI in the hippocampus of the two distinct mouse lines selected for stress reactivity. The results provide insight into the genetic background that might modulate the molecular response to antidepressants and mediate resistance to antidepressants.
Swiss-Webster male mice (weighing 25–30 g, 12 weeks old) exhibiting high (HA) and low (LA) swim stress-induced analgesia (SIA) were used . Briefly, the outbred male and female mice were scored for the latency of the hind paw nociceptive reflex on a 56°C hot plate at 2 min after the completion of a 3-min swimming test in 20°C water. The animals displaying the longest (50 – 60 s) and the shortest (<10 s) post-swim latencies of the nociceptive response were selected as progenitors of the HA and LA lines. A similar procedure was repeated in subsequent offspring, and the animals displaying the longest and the shortest post-swim hot plate latencies were mated to maintain the lines or selected for use in subsequent experiments. The HA and LA male mice used in the present experiment belonged to the 68th generation of albino Swiss–Webster mice that have been selectively bred in our laboratory for high (HA line) and low (LA line) magnitudes of swim-stress induced analgesia (SSIA) as described earlier .
Animals were males, 5 months old and weighed 41.5 ± 0.8 g and 39.5 ± 0.5 g (HA and LA lines, respectively). Mice from each line and treatment group were housed four to six per cage on a 12-h/12-h light dark cycle with unlimited access to standard murine chow food and water. All procedures were carried out between 9 am and 4 pm. Mice from HA and LA lines were injected and decapitated in a pseudorandom order.
HA and LA animals were injected intraperitoneally with desipramine (DMI, 7.5 mg/kg/day, Sigma-Aldrich, USA) or saline once a day for 21 consecutive days. of the experiment. Each experimental group contained 12 mice. DMI was selected based on the results obtained from our previous studies [5, 6]. These studies revealed that DMI shortened the immobility period of HA mice in the forced swim and tail suspension tests but was ineffective in LA animals. All procedures involving animals were performed in accordance with the Guiding Principles for the Care and Use of Research Animals, and were approved by the the Institutional Ethics Committee of the Institute of Genetics and Animal Breeding and then by the National Third Local Ethics Commission (permission No. 3/2005). Permission had been issued in accordance with the guidelines of the pertinent Polish Parliamentary Bill. All observations were performed by a trained observer who was blinded to the treatment conditions and was not informed of the treatment conditions in advance.
Two days after the injections, the animals were sacrificed through decapitation (HA DMI, n = 12; HA saline, n = 12; LA DMI, n = 12; LA saline, n = 12), and the brains were rapidly removed. The hippocampi were quickly dissected on ice-cold glass, inserted into freezing vials, frozen in liquid nitrogen and stored at -80°C until further analysis. Total RNA was isolated separately from each individual hippocampus using the NucleoSpin RNA II kit (Macherey-Nagel, Germany). A Nanodrop (Nanodrop, USA) and Bioanalyzer (Agilent, USA) were used to estimate the quantity and quality of each RNA sample, and the RIN (RNA Integrity Number) values ranged from 9.2 to 9.8 for all samples.
To minimize the influence of the individual differences between the animals and the variation introduced through dissection and tissue preparation, the individual RNA samples from three mice were pooled. The quantity and quality of the RNA in the pools were estimated using a Nanodrop and Bioanalyzer. Each pool was subsequently converted to biotinylated cRNA using the Illumina RNA Amplification Kit (Ambion Inc., USA) and 100 ng of total RNA, purified using the RNeasy kit (Qiagen, Germany) and hybridized to a single microarray. Four independent biological replicates of the microarray analysis were prepared for each group of mice.
Microarray, hybridization and fluorescent detection
The gene expression profiling was performed using MouseRef-8_V2 Expression BeadChip microarrays (Illumina, USA) according to MIAME guidelines . Each array contains >22,000 well-annotated RefSeq 50-mer oligonucleotide probes selected from the NCBI database (Release 16). The cRNA samples were applied to the arrays and assembled onto hybridization cartridges. Overall, 16 hybridizations (4 for each experimental group) were performed. The hybridization and washing of the arrays was performed according to the manufacturer’s hybridization protocol. The microarrays were scanned on a BeadArray Reader (Illumina, USA).
Data normalization, selection of the most differently expressed genes and hierarchical clustering
Raw microarray data were processed using the BeadArray and LIMMA package of the R statistical language (Bioconductor project; ). The data preprocessing involved the normalization of the expression levels using a quantile method preceded by log2 transformation. Linear model fitting was performed for the pre-processed dataset. The coefficients of the obtained fitted models described the differences between the RNA sources hybridized to the arrays. The empirical Bayesian analysis was performed to determine whether the contrast coefficients from the linear models were equal zero and identify differentially expressed genes. The genes showing log fold-changes greater than 0.5 (logFC > 0.5) and p < 0.05 were considered significantly differentially expressed. The Benjamini and Hochberg method , to control the false discovery rate (FDR), was used to correct the p values.
Gene lists enrichment analysis
The gene lists (GenBank accession numbers) generated from the microarray results were submitted to the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7. . The genes were classified into functional categories using biological process (BP), cellular components (CC), and molecular function (MF) ontologies. Analysis of the association of the genes with physiological or biochemical pathways was performed using the Kyoto Encyclopedia of Genes and Genomes database (KEGG; ). To identify significantly overrepresented biological categories and KEGG pathways within the lists of differentially expressed genes, the threshold for the enrichment analysis was set at p ≤ 0.05. For details see Lisowski et al., 2011 [10, 16].
Quantitative real-time PCR (qPCR) for microarray validation
To validate the results of the microarrays, quantitative real-time PCR (qPCR) with SYBR Green was performed as previously described , according to MIQE guidelines . Eight genes, belonging to different functional categories and significantly differing in the expression between the saline and desipramine experimental groups, were selected in each line. The qPCR assays were performed in triplicate on the same, non-pooled individual RNA samples, which were used in the microarray experiment. All amplified PCR fragments were sequenced to verify the resulting PCR product.
For reference, using a previously described methodology [18, 19], 10 common reference genes belonging to different functional classes (Gapdh, Tbp, Hprt1, Actb, B2m, Rpl13a, Hsp90ab1, Gusb, Tfrc, and Ppia) were assessed using qPCR, and two housekeeping genes stably expressed in hippocampi of DMI-treated mice were used: TATA box binding protein (Tbp) and hypoxanthine guanine phosphoribosyl transferase 1 (Hprt1). The normalization factor (NF) was obtained from the geometric mean of the raw expression data of Tbp and Hprt1.
The primers were designed using ExonPrimer software (Institute of Human Genetics, TUM/Helmholtz Center Munich, Germany; ) and Mus musculus GenBank sequences. All primers produced amplicons spanning two exons in highly conserved coding regions and included all known alternatively spliced mRNA variants. Additional file 1 summarizes the primer information, including gene names, primer sequences, amplicon lengths, annealing temperatures, and GenBank accession numbers.
The data from the triplicate qPCR reactions were normalized using the average cycle threshold (Ct) value and a mathematical model for the relative quantification of the real-time RT-PCR results . The qPCR data obtained from the saline and DMI-treated groups were analyzed using t-tests. The differences between the groups were considered significant at p < 0.05. The degree and significance of the correlation between fold-changes, as determined in the microarray and qPCR analyses, were evaluated using the Pearson moment correlation.
Genome-wide gene expression
Comparisons of the hippocampal transcriptomic profiles of the control and desipramine-treated animals revealed that DMI affected the expression of 71 genes in HA and 41 genes in LA mice, meeting the criteria of log fold-changes greater than 0.5 and adjusted p values <0.05. DMI treatment in HA mice, compared with saline injections, resulted in the upregulation of 57 transcripts and the downregulation of 14 transcripts. DMI treatment in LA animals resulted in upregulation of 28 transcripts and downregulation of 13 transcripts. The entire set of probes, with the corresponding gene names and the FC and p values that differed between the saline- and DMI-treated animals are presented in Additional files 2 and 3. All genes annotated to these probe sets were expressed in mouse hippocampus according to the Novartis Gene Expression Atlas . The cell-type classification analysis of differentially expressed genes, according to Cahoy , the GeneCards® database  screening and Ingenuity Pathway Analysis®  of the major canonical pathways, revealed that most of these genes are characteristic of neurons and oligodendrocytes.
Effect of desipramine on hippocampal transcriptome in HA line
DMI treatment in HA mice resulted in the upregulation of transcripts involved in 12 significant (p < 0.05) biological processes, 5 cellular components, 7 molecular function terms and 2 biochemical pathways. The biological process terms included neurotransmitter transport (p = 0.001) (Slc6a9, Slc17a6, Slc6a11, Rims3), cell motion (p = 0.002) (Vav3, Sema3f, Gbx2, Efna5, Elmo1), cell projection organization (p = 0.003) (Vav3, Ntng1, Gbx2, Efna5), angiogenesis (p = 0.005) (Gbx2, Arhgap24, Anxa2) axonogenesis/neuron differentiation/neuron projection (p = 0.004) (Ntng1, Gbx2, Efna5, Rora), behavior (p = 0.039) (Ptgds, Sema3f, Chrna4, Zic1) and neurogenesis (p = 0.048) (Sema3f, Ntng1, Gbx2, Efna5, Rora). The cellular component terms included elements of the cell junction (p = 0.0002) (Cbln1, Slc17a6, Syt9, Chrna4, Arhgap24, Calb2, Anxa2, Rims3), plasma membrane (p = 0.001) (Slc6a11, Ntng1, Syt9, Arhgap24, Calb2, Rims3, Anxa2, Ly6a, Gng8, Slc6a9, Cbln1, Slc17a6, Chrna4), and the synapse (p = 0.011) (Cbln1, Slc17a6, Syt9, Chrna4, Rims3). The molecular function terms included amino acid transmembrane transport (p = 0.007) (Slc6a9, Slc17a6, Slc6a11) and chemorepellent activity (p = 0.01) (Sema3f, Efna5). To identify significantly over-represented pathways in the list of DMI-upregulated genes, we searched the KEGG database. In the hippocampus of HA mice, the significantly overrepresented pathways included chemokine signaling (p = 0.003) (Gng8, Vav3, Prkcd, Elmo1) and axon guidance pathways (p = 0.021) (Sema3f, Ntng1, Efna5).
The downregulated transcripts were involved in 6 significant (p < 0.05) biological processes and 1 cellular component term, including mast cell (p = 0.009) and myeloid leukocyte activation (p = 0,02), the regulation of response to external stimulus (p = 0.023) and vesicle-mediated transport (p = 0.029) (Cplx2, Fcer1g, Npy, Arc) and coding elements of synapses (p = 0.026) (Phactr1, Arc, Cplx2).
Moreover, as shown in Additional file 1: Table S1, DMI treatment resulted in the upregulation of neuronatin (Nnat), a factor involved in the modulation of dendritic calcium levels during homeostatic plasticity, and the downregulation of neuropeptide Y (Npy), which is implicated in a diverse range of biological actions, including anxiolysis and the control of gonadotrophin-releasing hormone or pituitary hormone secretion. DMI also caused the downregulation of corticotropin-releasing hormone-binding protein (Crhbp), which inactivates CRH. For the list of GO functional category classifications of up- and down-regulated genes affected after DMI treatment in the HA line see Additional file 4.
Effect of desipramine on hippocampal transcriptome in the LA line
The genes upregulated in LA mice after DMI treatment are involved in 13 significant (p < 0.05) biological processes and 2 cellular component terms. The biological process terms included behavior (p = 0.021), including behavioral interactions (p = 0.045), mating (p = 0.016), reproduction (p = 0.028) and locomotory behavior (p = 0.048) (Pde1b, Ppp1r1b, Enpp2, Drd1a), brain structure development, including striatum (p = 0.015), subpallium (0.019), eye (p = 0.022) and forebrain (p = 0.025), and neurogenesis (p = 0.03) (Otx2, Rarb, Drd1a). The cellular component terms included elements of the extracellular space (p = 0.006) and extracellular region (p = 0.035) (Ttr, Kl, Enpp2, Sostdc1, 1500015o10rik, Tac1, Igf2, Igfbp2). The molecular function terms and biochemical pathways were not included in the list of upregulated genes in LA mice and were not statistically significant. The downregulated genes in the LA mice transcriptomes were not significantly enriched.
Notably, the transcripts upregulated after DMI treatment in LA mice include transthyretin (Ttr), the prolactin receptor (Prlr), insulin-like growth factor 2 (Igf2), and insulin-like growth factor binding protein 2 (Igfbp2), which are important genes involved in neuroprotection. For the list of GO functional category classifications of up- and down-regulated genes affected after DMI treatment in the LA line see Additional file 5.
Quantitative real-time PCR (qPCR) for microarray validation
To validate the microarray results, 16 (8 from HA and 8 from LA) differentially expressed genes (DEGs) were selected for quantification using qPCR. The DEGs were selected from among the significant functional groups associated with neurotransmitter transport, immunity, cell differentiation, neurogenesis, axonogenesis, and calcium signaling and axon guidance pathways. The following genes were selected for the HA line: Chrna4, Crhbp, Gbx2, Igf2, Nnat, Npy, Slc6a11, and Slc17a6. The following genes were selected for the LA line: Cckbr, Csrp1, Drd1a, Igf2, Igfbp2, Otx2, Prlr, and Ttr. Differences in the expression revealed using microarray analyses were confirmed. The expression levels are shown in Figure 1. The DEGs identified through microarray tests were consistent with the qPCR results (R = 0.94, p < 0.0001).
Early comparisons of the constitutive gene expression profiles between the HA and LA lines identified 205 transcripts with different levels of mRNA in the hippocampi of non-stressed animals . The results suggest that individual differences in depression-like behavioral patterns may be associated with the characteristics of the constitutive hippocampal transcriptome.
The present study indicates that two different basal hippocampal transcriptomic profiles influenced the response to DMI, although common changes induced through drug treatments were also observed. The putative genes associated with the action of the ADs in the brain based on the literature search indicated that chronic AD treatments increase the levels of the cAMP response element-binding protein (CREB) mRNA and its receptor trkB in the hippocampus [26–28]. Furthermore, other transcription factors, such as Fos family members, are also increased in the brain after the chronic administration of ADs [29, 30]. Antidepressant drugs produce antidepressant effects in the prefrontal cortex and hippocampus . Experiments have demonstrated that antidepressants significantly affect the cAMP cascade, and a homologous pathway exists for all classes of antidepressants. However, other observations, using different biochemical and molecular assays, suggest limited changes. Thus, genes encoding cytochrome b , glyceraldehyde-3-phosphate dehydrogenase (GAPDH) , the RING-H2 finger motif gene , the cysteine string protein , the tetraspan protein , the GRB2 homeobox protein and the ribosomal protein L35a , isoforms of protein kinase C and cAMP-dependent protein kinase , the Mss4 gene , VAMP2/synaptobrevin-2 , and GAP-43  are regulated through treatment with antidepressants in different parts of the rodent brain. Long-lasting neuroadaptations likely include complex changes in the gene expression in the limbic system of the brain. For example, changes in the expression of Crhbp and Npy (downregulated in HA) after AD treatment have been reported . Evidence of the importance of ion transport in the rat brain cortex has been shown in the context of antidepressant treatment [43, 44]. Specifically, changes in gene expression coding elements of calcium and sodium channels were observed after 28 days of treatment with amitryptiline . In the present study, downregulation of adenylate cyclase 1 (Adcy1) and the cholecystokinin B receptor (Cckbr) involved in calcium signaling occurred in the LA line. In the HA line, we observed the upregulation of neuronatin (Nnat). A recent study showed that Nnat adjusts dendritic calcium levels through the regulation of intracellular calcium storage, thus neuronatin might impact synaptic plasticity through the modulation of dendritic calcium levels during homeostatic plasticity, thereby potentially regulating neuronal excitability, receptor trafficking, and calcium dependent signaling in hippocampal neurons . The results obtained in the present study, however, also implicate a number of novel “AD genes”. Among these, alterations in the expression of genes involved in neurogenesis (Sema3f, Ntng1, Gbx2, Efna5, Rora, Otx2, Rarb, Drd1a) and neuroprotection (Ttr, Prl, Igf2) are of particular interest.
Effects of desipramine on genes related to neurogenesis
The chronic administration of antidepressants increases the proliferation and survival of neural stem cells and new neurons in the hippocampus [47–49]. Increased stem cell proliferation is required for the behavioral effects of antidepressants, suggesting that this effect might underlie the therapeutic action of antidepressant drugs on mood . It takes weeks for newborn neurons to differentiate and mature into fully functional hippocampal granule neurons , suggesting that the gradual maturation of newly born dentate granule neurons might explain the delayed action of antidepressant drugs . Specific gene expression profiles of neurogenesis in AD treatment have not been clearly identified, although according to Gene Ontology  classification, the upregulation of five genes in the HA line (Sema3F, Ntng1, Gbx2, Efna5, Rora) and three genes in the LA line (Otx2, Rarb, Drd1a) involved in neurogenesis, neurite outgrowth, and the proliferation and differentiation of neuronal progenitor cells were detected. Igf2 should be added to the group of genes related to neurogenesis, as this gene was upregulated in both lines after DMI treatment.
Effects of desipramine on genes related to the chronic stress response
The analyses performed in the present work indicated a robust upregulation of Ttr gene-coding transthyretin (TTR) and the Prlr gene-coding prolactin receptor in the LA line, and in both lines, we observed the upregulation of the Igf2 gene encoding somatomedin. TTR is a carrier protein and a major transporter of thyroid hormones and retinol in plasma and cerebrospinal fluid . The clinical features of defects in TTR include seizures, stroke-like episodes, dementia and psychomotor deterioration. However, the antidepressant-like effect of the histone deacetylase inhibitor sodium butyrate (SB) upregulates transthyretin (Ttr), with simultaneous increases in histone H4 acetylation at the promoter of the Ttr gene . Moreover, the genes, including Ttr, Prl and Igf2, whose hippocampal expression patterns were downregulated in four rat models of endogenous depression and chronic stress, represent a generalized molecular response to chronic stress . In addition, stress has previously been shown to increase the expression of neuropeptide Y [57–59], which in the present study was downregulated in the HA line after DMI treatment. Prolactin is a neuromodulator of behavioral and neuroendocrine stress in the rat [60, 61]. The downregulation of brain prolactin receptors increased anxiety-related behavior, demonstrating an anxiolytic effect of PRL in the brain . Furthermore the stress-induced increase of corticotropin secretion was reduced after the chronic intracerebroventricular infusion of PRL . Fujikawa  showed that PRL levels increase in response to stress acting on the central nervous system to protect against acute stress-induced hypocalcemia. Moreover, a recent study on gene expression profiling in neural stem cells (NSCs) showed the upregulation of insulin-like growth factor 2 (Igf2) in the dentate gyrus (DG) compared with that in immature neurons. IGF2 selectively controls the proliferation of DG NSCs in vitro and in vivo through AKT-dependent signaling. Thus, the gene expression profiling of NSCs and progeny cells revealed that IGF2 is a novel regulator of adult neurogenesis . Altogether, these findings suggest the possibility that the upregulation of transcripts for Ttr, Prlr in HA line, and Igf2 in both lines might be involved in AD actions.
Effects of desipramine on the genes involved in neurotransmitter transport
Among the genes involved in neurotransmission, microarray analysis revealed the upregulation of Slc6a, Slc17a6, Slc6a11 (also Glyt1, Vglut2 Gat3 respectively) and Rab3. Glyt1 terminates the action of glycine through high affinity sodium-dependent reuptake into presynaptic terminals and plays a role in the regulation of glycine levels in NMDA receptor-mediated neurotransmission. Vglut2 mediates the uptake of glutamate into synaptic vesicles at presynaptic nerve terminals of excitatory neural cells. Gat3 terminates the action of GABA through high affinity sodium-dependent reuptake into presynaptic terminals and Rab3 regulates synaptic membrane exocytosis. Because antidepressants inhibit GABA uptake, the upregulation of Gat3 is an important factor in AD therapy. Tordera  showed that repeated treatment with fluoxetine, paroxetine or desipramine increased the abundance of Vglut1 mRNA in the hippocampus, but did not increase the expression of Vglut2 and Vglut3 mRNA. The data obtained in the present study suggest that a course of AD drug treatment also increases the expression of Vglut2, another gene involved in the regulation of glutamate secretion, but only in vulnerable individuals, such as the HA mouse line.
The impact of desipramine treatment on HA and LA mouse lines
The objective of the present study was to characterize the transcriptional response to chronic desipramine treatment in the hippocampus of mouse lines displaying robust differences in the behavioral responses to ADs. In a previous study, we showed that desipramine, venlafaxine and aminosenktide shortened the immobility period of HA mice in the forced swim and tail suspension tests, but were ineffective in the LA line [5–7]. Moreover, at the molecular level, the lines were characterized by differences in the basal transcriptomic profile in the hippocampus, which indicates that selective breeding for swim SIA affected many aspects of hippocampal neuron physiology, including metabolism, structural changes, and cellular signaling . The genes involved in calcium signaling pathways are upregulated in LA mice, while HA mice are characterized by the robust upregulation of Klf5 and genes encoding receptors for neurotensin and GABA . Klf5 encodes the KLF5 zinc finger protein, which act as transcription factors that bind to GC box promoter elements and activate the transcription of these genes, depending on the DNA methylation status in the promoter regions. In turn, we observed that chronic mild stress (CMS) generated a self-contained response within each mouse line . In the HA line, CMS affected several genes involved in chromatin modifications, suggesting a role for the epigenetic regulation of the stress response in the HA line. KLF5 is primarily localized in the hippocampus, which is the principal source of glutamatergic neurotransmission. These findings suggest that KLF5 gene expression is involved in glutamatergic neurotransmission . Altogether, the previous and current results on the hippocampal transcription patterns of the genes involved in calcium signaling, GABA neurotransmission, ion transport, and glutamate transporters suggest a functional endpoint for structurally unrelated antidepressants and explain the differences observed in the lines in response to antidepressants.
The gene expression patterns during ADs treatment could be very complex. Nevertheless, our study revealed a number of novel potential targets for antidepressant therapy, such as the Npy, Crhbp, Ttr, Igfbp2, Prlr, and genes involved in neurogenesis that were not previously linked to antidepressant action. Potential drug targets against the detrimental effects of stress and depression include also glutamate transporters, GABA neurotransmission, and neuro- and axonogenesis, and suggest that pharmacotherapy might be genotype dependent. Single nucleotide polymorphisms (SNPs), epigenetic statuses or epimutations in the regulatory elements of the identified genes should be considered in the pharmacogenomics of antidepressants. Genome wide gene expression profiling using two lines of mice bred for stress reactivity facilitates the elucidation of broader mechanisms underlying antidepressants effects.
Femenia T, Gomez-Galan M, Lindskog M, Magara S: Dysfunctional hippocampal activity affects emotion and cognition in mood disorders. Brain Res. 2012, 1476: 58-70.
Price JL, Drevets WC: Neural circuits underlying the pathophysiology of mood disorders. Trends Cogn Sci. 2012, 16 (1): 61-71. 10.1016/j.tics.2011.12.011.
Robinson MJ, Edwards SE, Iyengar S, Bymaster F, Clark M, Katon W: Depression and pain. Front Biosci. 2009, 14: 5031-5051. 10.2741/3585.
Panocka I, Marek P, Sadowski B: Inheritance of stress-induced analgesia in mice - selective breeding study. Brain Res. 1986, 397 (1): 152-155. 10.1016/0006-8993(86)91379-X.
Juszczak GR, Sliwa AT, Wolak P, Tymosiak-Zielinska A, Lisowski P, Swiergiel AH: The usage of video analysis system for detection of immobility in the tail suspension test in mice. Pharmacol Biochem Behav. 2006, 85 (2): 332-338. 10.1016/j.pbb.2006.08.016.
Juszczak GR, Lisowski P, Sliwa AT, Swiergiel AH: Computer assisted video analysis of swimming performance in a forced swim test: Simultaneous assessment of duration of immobility and swimming style in mice selected for high and low swim-stress induced analgesia. Physiol Behav. 2008, 95 (3): 400-407. 10.1016/j.physbeh.2008.07.003.
Panocka I, Massi M, Lapo I, Swiderski T, Kowalczyk M, Sadowski B: Antidepressant-type effect of the NK3 tachykinin receptor agonist aminosenktide in mouse lines differing in endogenous opioid system activity. Peptides. 2001, 22 (7): 1037-1042. 10.1016/S0196-9781(01)00438-7.
Lisowski P, Stankiewicz AM, Goscik J, Wieczorek M, Zwierzchowski L, Swiergiel AH: Selection for stress-induced analgesia affects the mouse hippocampal transcriptome. J Mol Neurosci. 2012, 47 (1): 101-112. 10.1007/s12031-011-9692-2.
Willner P, Towell A, Sampson D, Sophokleous S, Muscat R: Reduction of sucrose preference by chronic unpredictable mild stress, and its restoration by a tricyclic antidepressant. Psychopharmacology (Berl). 1987, 93 (3): 358-364.
Lisowski P, Juszczak GR, Goscik J, Wieczorek M, Zwierzchowski L, Swiergiel AH: Effect of chronic mild stress on hippocampal transcriptome in mice selected for high and low stress-induced analgesia and displaying different emotional behaviors. Eur Neuropsychopharmacol. 2011, 21 (1): 45-62. 10.1016/j.euroneuro.2010.08.004.
Edgar R, Barrett T: NCBI GEO standards and services for microarray data. Nat Biotechnol. 2006, 24 (12): 1471-1472. 10.1038/nbt1206-1471.
Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc B Methodol. 1995, 57 (1): 289-300.
Database for Annotation, Visualization and Integrated Discovery. http://david.abcc.ncifcrf.gov.
Kyoto Encyclopedia of Genes and Genomes. http://www.genome.jp/kegg/pathway.html.
Lisowski P, Wieczorek M, Goscik J, Juszczak GR, Stankiewicz AM, Zwierzchowski L, Swiergiel AH: Effects of chronic stress on prefrontal cortex transcriptome in mice displaying different genetic backgrounds. J Mol Neurosci. 2012, 50() (1): 33-57. doi: 10.1007/s12031-012-9850-1. Epub 2012 Jul 27
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, et al: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55 (4): 611-622. 10.1373/clinchem.2008.112797.
Lisowski P, Pierzchala M, Goscik J, Pareek CS, Zwierzchowski L: Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid. J Appl Genet. 2008, 49 (4): 367-372. 10.1007/BF03195635.
Lisowski P, Robakowska-Hyzorek D, Blitek A, Kaczmarek MM, Gajewska A, Kochman H, Zwierzchowski L, Ziecik AJ, Kochman K: Development of real-time PCR assays in the study of gonadotropin subunits, follistatin and prolactin genes expression in the porcine anterior pituitary during the preovulatory period. Neuroendocrinol Lett. 2008, 29 (6): 958-964.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29 (9): e45-10.1093/nar/29.9.e45.
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-278. 10.1523/JNEUROSCI.4178-07.2008.
Ingenuity Pathway Analysis. http://www.ingenuity.com/.
Butterweck V, Winterhoff H, Herkenham M: St John’s wort, hypericin, and imipramine: a comparative analysis of mRNA levels in brain areas involved in HPA axis control following short-term and long-term administration in normal and stressed rats. Mol Psychiatry. 2001, 6 (5): 547-564. 10.1038/sj.mp.4000937.
Dowlatshahi D, MacQueen GM, Wang JF, Young LT: Increased temporal cortex CREB concentrations and antidepressant treatment in major depression. Lancet. 1998, 352 (9142): 1754-1755. 10.1016/S0140-6736(05)79827-5.
Saarelainen T, Hendolin P, Lucas G, Koponen E, Sairanen M, MacDonald E, Agerman K, Haapasalo A, Nawa H, Aloyz R, et al: Activation of the TrkB neurotrophin receptor is induced by antidepressant drugs and is required for antidepressant-induced behavioral effects. J Neurosci. 2003, 23 (1): 349-357.
Nestler EJ, Kelz MB, Chen JS: Delta FosB: a molecular mediator of long-term neural and behavioral plasticity. Brain Res. 1999, 835 (1): 10-17. 10.1016/S0006-8993(98)01191-3.
Hope BT, Nye HE, Kelz MB, Self DW, Iadarola MJ, Nakabeppu Y, Duman RS, Nestler EJ: Induction of a long-lasting AP-1 complex composed of altered Fos-like proteins in brain by chronic cocaine and other chronic treatments. Neuron. 1994, 13 (5): 1235-1244. 10.1016/0896-6273(94)90061-2.
Manji HK, Drevets WC, Charney DS: The cellular neurobiology of depression. Nature medicine. 2001, 7 (5): 541-547. 10.1038/87865.
Huang NY, Strakhova M, Layer RT, Skolnick P: Chronic antidepressant treatments increase cytochrome b mRNA levels in mouse cerebral cortex. J Mol Neurosci. 1997, 9 (3): 167-176. 10.1007/BF02800499.
Tohda M, Qi Z, Watanabe H: Influence of chronic treatment with imipramine on mRNA levels in rat brain: elevation of glyceraldehyde-3-phosphate dehydrogenase levels. Jpn J Pharmacol. 1999, 81 (4): 393-396. 10.1254/jjp.81.393.
Yamada M, Yamada M, Yamazaki S, Takahashi K, Nishioka G, Kudo K, Ozawa H, Yamada S, Kiuchi Y, Kamijima K, et al: Identification of a novel gene with RING-H2 finger motif induced after chronic antidepressant treatment in rat brain. Biochem Biophys Res Commun. 2000, 278 (1): 150-157. 10.1006/bbrc.2000.3773.
Yamada M, Yamada M, Yamazaki S, Takahashi K, Nara K, Ozawa H, Yamada S, Kiuchi Y, Oguchi K, Kamijima K, et al: Induction of cysteine string protein after chronic antidepressant treatment in rat frontal cortex. Neurosci Lett. 2001, 301 (3): 183-186. 10.1016/S0304-3940(01)01638-X.
Hua LV, Green M, Wong A, Warsh JJ, Li PP: Tetraspan protein CD151: a common target of mood stabilizing drugs?. Neuropsychopharmacology. 2001, 25 (5): 729-736. 10.1016/S0893-133X(01)00269-X.
Landgrebe J, Welzl G, Metz T, van Gaalen MM, Ropers H, Wurst W, Holsboer F: Molecular characterisation of antidepressant effects in the mouse brain using gene expression profiling. J Psychiatr Res. 2002, 36 (3): 119-129. 10.1016/S0022-3956(01)00061-9.
Rausch JL, Gillespie CF, Fei Y, Hobby HM, Stoming T, Ganapathy V, Leibach FH: Antidepressant effects on kinase gene expression patterns in rat brain. Neurosci Lett. 2002, 334 (2): 91-94. 10.1016/S0304-3940(02)01106-0.
Andriamampandry C, Muller C, Schmidt-Mutter C, Gobaille S, Spedding M, Aunis D, Maitre M: Mss4 gene is up-regulated in rat brain after chronic treatment with antidepressant and down-regulated when rats are anhedonic. Mol Pharmacol. 2002, 62 (6): 1332-1338. 10.1124/mol.62.6.1332.
Yamada M, Takahashi K, Tsunoda M, Nishioka G, Kudo K, Ohata H, Kamijima K, Higuchi T, Momose K, Yamada M: Differential expression of VAMP2/synaptobrevin-2 after antidepressant and electroconvulsive treatment in rat frontal cortex. Pharmacogenomics J. 2002, 2 (6): 377-382. 10.1038/sj.tpj.6500135.
Chen B, Wang JF, Sun XJ, Young LT: Regulation of GAP-43 expression by chronic desipramine treatment in rat cultured hippocampal cells. Biol Psychiatry. 2003, 53 (6): 530-537. 10.1016/S0006-3223(02)01551-2.
Conti B, Maier R, Barr AM, Morale MC, Lu X, Sanna PP, Bilbe G, Hoyer D, Bartfai T: Region-specific transcriptional changes following the three antidepressant treatments electro convulsive therapy, sleep deprivation and fluoxetine. Mol Psychiatry. 2007, 12 (2): 167-189. 10.1038/sj.mp.4001897.
Knuuttila JEA, Toronen P, Castren E: Effects of antidepressant drug imipramine on gene expression in rat prefrontal cortex. Neurochem Res. 2004, 29 (6): 1235-1244.
Palotas M, Palotas A, Puskas LG, Kitajka K, Pakaski M, Janka Z, Molnar J, Penke B, Kalman J: Gene expression profile analysis of the rat cortex following treatment with imipramine and citalopram. Int J Neuropsychopharmacol. 2004, 7 (4): 401-413. 10.1017/S1461145704004493.
Boehm C, Newrzella D, Herberger S, Schramm N, Eisenhardt G, Schenk V, Sonntag-Buck V, Sorgenfrei O: Effects of antidepressant treatment on gene expression profile in mouse brain: cell type-specific transcription profiling using laser microdissection and microarray analysis. J Neurochem. 2006, 97: 44-49.
Oyang EL, Davidson BC, Lee W, Poon MM: Functional characterization of the dendritically localized mRNA neuronatin in Hippocampal neurons. Plos One. 2011, 6 (9): e24879-10.1371/journal.pone.0024879. doi: 10.1371/journal.pone.0024879. Epub 2011 Sep 14
Jacobs BL, van Praag H, Gage FH: Adult brain neurogenesis and psychiatry: a novel theory of depression. Mol Psychiatry. 2000, 5 (3): 262-269. 10.1038/sj.mp.4000712.
Madsen TM, Treschow A, Bengzon J, Bolwig TG, Lindvall O, Tingstrom A: Increased neurogenesis in a model of electroconvulsive therapy. Biol Psychiatry. 2000, 47 (12): 1043-1049. 10.1016/S0006-3223(00)00228-6.
Malberg JE, Eisch AJ, Nestler EJ, Duman RS: Chronic antidepressant treatment increases neurogenesis in adult rat hippocampus. J Neurosci. 2000, 20 (24): 9104-9110.
Santarelli L, Saxe M, Gross C, Surget A, Battaglia F, Dulawa S, Weisstaub N, Lee J, Duman R, Arancio O, et al: Requirement of hippocampal neurogenesis for the behavioral effects of antidepressants. Science. 2003, 301 (5634): 805-809. 10.1126/science.1083328.
van Praag H, Schinder AF, Christie BR, Toni N, Palmer TD, Gage FH: Functional neurogenesis in the adult hippocampus. Nature. 2002, 415 (6875): 1030-1034. 10.1038/4151030a.
Wang JW, David DJ, Monckton JE, Battaglia F, Hen R: Chronic fluoxetine stimulates maturation and synaptic plasticity of adult-born hippocampal granule cells. J Neurosci. 2008, 28 (6): 1374-1384. 10.1523/JNEUROSCI.3632-07.2008.
Gene Ontology. http://www.geneontology.org/.
Landers KA, McKinnon BD, Li HK, Subramaniam VN, Mortimer RH, Richard K: Carrier-mediated thyroid hormone transport into placenta by placental transthyretin. J Clin Endocrinol Metab. 2009, 94 (7): 2610-2616. 10.1210/jc.2009-0048.
Yamawaki Y, Fuchikami M, Morinobu S, Segawa M, Matsumoto T, Yamawaki S: Antidepressant-like effect of sodium butyrate (HDAC inhibitor) and its molecular mechanism of action in the rat hippocampus. World J Biol Psychiatry. 2012, 13 (6): 458-467. 10.3109/15622975.2011.585663.
Andrus BM, Blizinsky K, Vedell PT, Dennis K, Shukla PK, Schaffer DJ, Radulovic J, Churchill GA, Redei EE: Gene expression patterns in the hippocampus and amygdala of endogenous depression and chronic stress models. Mol Psychiatry. 2012, 17 (1): 49-61. 10.1038/mp.2010.119.
Zambello E, Fuchs E, Abumaria N, Rygula R, Domenici E, Caberlotto L: Chronic psychosocial stress alters NPY system: different effects in rat and tree shrew. Prog Neuropsychopharmacol Biol Psychiatry. 2010, 34 (1): 122-130. 10.1016/j.pnpbp.2009.10.011.
Thorsell A, Slawecki CJ, El Khoury A, Mathe AA, Ehlers CL: The effects of social isolation on neuropeptide Y levels, exploratory and anxiety-related behaviors in rats. Pharmacol Biochem Behav. 2006, 83 (1): 28-34. 10.1016/j.pbb.2005.12.005.
Heilig M, Widerlov E: Neurobiology and clinical aspects of Neuropeptide-Y. Crit Rev Neurobiol. 1995, 9 (2–3): 115-136.
Blume A, Torner L, Liu Y, Subburaju S, Aguilera G, Neumann ID: Prolactin activates mitogen-activated protein kinase signaling and corticotropin releasing hormone transcription in Rat hypothalamic neurons. Endocrinology. 2009, 150 (4): 1841-1849. 10.1210/en.2008-1023.
Torner L, Toschi N, Pohlinger A, Landgraf R, Neumann ID: Anxiolytic and anti-stress effects of brain prolactin: Improved efficacy of antisense targeting of the prolactin receptor by molecular modeling. J Neurosci. 2001, 21 (9): 3207-3214.
Fujikawa T, Soya H, Tamashiro KLK, Sakai RR, McEwen BS, Nakai N, Ogata M, Suzuki I, Nakashima K: Prolactin prevents acute stress-induced hypocalcemia and ulcerogenesis by acting in the brain of rat. Endocrinology. 2004, 145 (4): 2006-2013.
Bracko O, Singer T, Aigner S, Knobloch M, Winner B, Ray J, Clemenson GD, Suh H, Couillard-Despres S, Aigner L, et al: Gene expression profiling of neural stem cells and their neuronal progeny reveals IGF2 as a regulator of adult Hippocampal neurogenesis. J Neurosci. 2012, 32 (10): 3376-3387. 10.1523/JNEUROSCI.4248-11.2012.
Tordera RM, Pei Q, Sharp T: Evidence for increased expression of the vesicular glutamate transporter, VGLUT1, by a course of antidepressant treatment. J Neurochem. 2005, 94 (4): 875-883. 10.1111/j.1471-4159.2005.03192.x.
Yanagi M, Hashimoto T, Kitamura N, Fukutake M, Komure O, Nishiguchi N, Kawamata T, Maeda K, Shirakawa O: Expression of Kruppel-like factor 5 gene in human brain and association of the gene with the susceptibility to schizophrenia. Schizophr Res. 2008, 100 (1–3): 291-301.
This research was supported through funding from the following grants: European Commission Framework 6 Integrated Project NEWMOOD (LSHMCT-2004-503474), and Polish Ministry of Science and Higher Education Grants: Iuventus Plus (IP2011 030371), and Mobility Plus (DPN/MOB109/II/2012).
The authors declare no competing interests.
PL carried out the RNA extractions, RNA quality assessment, microarray analysis, gene expression data extraction, and normalization, bioinformatics, drafted the manuscript and was responsible for the study concept and design. GRJ carried out the DMI treatment. JG performed the statistical analysis. AMS carried out the qPCR assays. MW dissected the brain structures. LZ participated in the design of the study. AHS supervised the project and contributed to the final draft of the paper. All authors read and approved the final manuscript.