Skip to main content

Advertisement

The effects of olanzapine on genome-wide DNA methylation in the hippocampus and cerebellum

Clinical Epigeneticsvolume 6, Article number: 1 (2014) | Download Citation

Abstract

Background

The mechanism of action of olanzapine in treating schizophrenia is not clear. This research reports the effects of a therapeutic equivalent treatment of olanzapine on DNA methylation in a rat model in vivo.

Genome-wide DNA methylation was assessed using a MeDIP-chip analysis. All methylated DNA immunoprecipitation (MeDIP), sample labelling, hybridization and processing were performed by Arraystar Inc (Rockville, MD, USA). The identified gene promoters showing significant alterations to DNA methylation were then subjected to Ingenuity Pathway Analysis (Ingenuity System Inc, CA, USA).

Results

The results show that olanzapine causes an increase in methylation in 1,140, 1,294 and 1,313 genes and a decrease in methylation in 633, 565 and 532 genes in the hippocampus, cerebellum and liver, respectively. Most genes affected are tissue specific. Only 41 affected genes (approximately 3%) showed an increase and no gene showed a decrease in methylation in all three tissues. Further, the two brain regions shared 123 affected genes (approximately 10%). The affected genes are enriched in pathways affecting dopamine signalling, molecular transport, nervous system development and functions in the hippocampus; ephrin receptor signalling and synaptic long-term potentiation in the cerebellum; and tissue morphology, cellular assembly and organization in the liver. Also, the affected genes included those previously implicated in psychosis.

Conclusions

The known functions of affected genes suggest that the observed epigenetic changes may underlie the amelioration of symptoms as well as accounting for certain adverse effects including the metabolic syndrome. The results give insights into the mechanism of action of olanzapine, therapeutic effects and the side effects of antipsychotics.

Background

Schizophrenia is one of the most devastating of psychiatric disorders [1]. The treatment of schizophrenia requires the suppression of hallucinations, delusions, agitation and the behavioural problems that accompany these symptoms [2]. Psychotherapy and rehabilitation can be undertaken when the acute symptoms start to subside through antipsychotic drug treatment.

The first antipsychotic drug, chlorpromazine, introduced in the early 1950s, was a major breakthrough because, unlike previously used sedative drugs, it could ameliorate hallucinations and delusions without overly sedating the patient [3]. Many other antipsychotic drugs were subsequently introduced [4], but these have not significantly advanced the treatment of schizophrenia. The early promise of the second-generation antipsychotics (atypical antipsychotics), such as clozapine and olanzapine, has been replaced by an acceptance that they are no more effective than the first-generation drugs [5]. However, second-generation antipsychotics have recently shown positive effects on verbal cognition [6]. Second-generation drugs have fewer neurological side effects but, unfortunately, many induce weight gain and the metabolic syndrome [7–10].

Our current understanding of the cause of schizophrenia is based on the pharmacological effects of the antipsychotic drugs used to treat the illness: they all bind to post-synaptic dopamine receptors especially D2 and the affinity at D2 receptors is both necessary and sufficient for the antipsychotic effects [11]. This, coupled with the observation that drugs that release dopamine into the synaptic cleft can induce the positive symptoms of schizophrenia (behaviours and feelings that are not real but imaginary), led to the hypothesis that excessive dopamine transmission in certain brain regions may cause the symptoms of schizophrenia [12]. A post-synaptic blockade occurs rapidly after a person ingests an antipsychotic drug. In contrast, the therapeutic effects of antipsychotics take days or weeks to accrue [13]. This suggests that downstream effects are important. One possibility is that the post-synaptic dopamine blockade causes a downstream cascade that has a therapeutic effect through altered gene transcription [14, 15]. A downstream effect, such as altered transcription, would explain the delay in the onset of therapeutic action. Other clinical observations also demonstrate the need for a more complex model than a post-synaptic dopamine blockade. Patients frequently fail to respond to an antipsychotic but subsequently show a robust response to a different drug despite the fact that both block the D2 receptor [16]. Moreover, many patients with schizophrenia show only a partial response to antipsychotics or fail to respond at all [17]. A refinement of the dopamine hypothesis proposes that an increase of D2 levels in the striatum may cause hallucinations and delusions and reduced D1 levels in the frontal lobes may cause cognitive deficits [18, 19]. This model is compatible with the delayed treatment effect but cannot explain the individual responsiveness to antipsychotics.

Epigenetic changes are another mechanism used to explain these clinical observations. They also offer an alternative therapeutic target for this serious disease: after 50 years of frustration we need to move our focus beyond post-synaptic dopamine receptors. Epigenetic changes associated with a drug can alter the expression of a single or variable number of genes without altering the gene sequence(s) [20]. Specifically, DNA methylation is a core epigenetic mechanism that involves the covalent binding of a methyl group to the 5-carbon position of cytosine leading to altered gene expression [21]. It is influenced by stochastic events including exposure to a variety of environments such as drug treatment [22, 23]. If DNA methylation plays a role in drug response, the drug or its metabolite must modify the methylation profile of the genome [20].

Limited research demonstrates that antipsychotic drugs can alter DNA methylation and gene expression [24]. However, most of this research has been conducted using variable post-mortem human brain tissues [25] and inappropriate non-brain cell types [26] that are not always ideal. Therefore, we have assessed the effects of a therapeutic equivalent dose of olanzapine (2.5 mg/kg per day for 21 days), a commonly used antipsychotic, on DNA methylation in rat brains using rat methylation arrays. The results demonstrate for the first time that the effect of olanzapine on DNA methylation is widespread and tissue specific, which may account for its efficacy and adverse effects.

Results

First, we assessed the locomotor activity of rats split into olanzapine-treated and vehicle-treated groups. Activity was significantly decreased (P = 0.001) in the olanzapine-treated group compared to the vehicle-treated group (Additional file 1: Figure S1a). Further, olanzapine-treated rats significantly increased in weight (P = 0.004) compared to the control group (Additional file 1: Figure S1b). Second, we assessed gene-specific DNA methylation across (almost) all genes in response to olanzapine treatment on two brain regions (hippocampus and cerebellum) and the liver, as detailed below.

Olanzapine causes widespread and tissue-specific changes in genome-wide methylation in rats in vivo

Widespread changes in gene-specific DNA methylation were apparent in all three tissues studied (hippocampus, cerebellum and liver), as shown by the heat map for the hippocampus (Figure 1A). The results identified genes where there was an increase or a decrease in methylation in drug-treated rats compared to controls. Specifically, almost twice as many genes showed an increase compared to the genes that showed a decrease in methylation, in response to olanzapine in each of the three tissues (Additional file 2: Table S1A, Additional file 3: Table S1B, Additional file 4: Table S2A, Additional file 5: Table S2B, Additional file 6: Table S3A and Additional file 7: Table S3B). Also, the set of genes affected differs across the three tissues. Approximately 75% of genes with an increase (Figure 1B) and over 90% of the genes with a decrease (Figure 1C) in methylation, following olanzapine treatment, were specific to a given tissue. Further, there was a small number of genes with a similar pattern of increase (total 164) or decrease (total 24) in the two brain regions and a smaller number with an increase (41) in all three tissues. The tissue-specific results are novel and were further assessed as follows.

Figure 1
figure1

Heat map of differential DNA methylation and the number of genes showing a) increased and b) decreased methylation. (A) Heat map of differential DNA methylation enrichment peaks of genes in the hippocampus, following olanzapine treatment of rats. The average normalized log2-ratio scan values were used to calculate the M' value (M' = Average(log2MeDIPE/InputE) – Average(log2MeDIPC/InputC)) for each probe. Venn diagrams depicting the number of genes that show an increase (B) or decrease (C) in methylation in brain tissues and the liver, following olanzapine treatment of rats.

Full size image

Pathways and associated network functions of genes that had changes in methylation in the hippocampus following olanzapine treatment

The genes (Additional file 2: Table S1A and Additional file 3: Table S1B) that had an increase (Table 1a) or a decrease (Table 1b) in methylation in the rat hippocampus following olanzapine treatment were assessed by pathway analysis. Genes with increased methylation were predominantly enriched in the dopamine-DARPP32 feedback in cAMP signalling canonical pathway (P = 1.6 × 10–3). The associated network functions that were affected included metabolic diseases and neurological disorders (Table 1, Additional file 1: Figure S2). In addition to the changes caused by olanzapine in methylation in psychosis-related canonical pathways, the results showed decreased methylation of genes involved in CDC42 and calcium signalling (P = 2.5 × 10–3) in the hippocampus. These genes affect nervous system development and function. Also, the cellular effects of sildenafil (Viagra) were revealed as an interesting canonical pathway (P = 9.8 × 10–3) for the 123 genes that had increased methylation as a result of olanzapine treatment in the two brain regions (Additional file 1: Figure S3).

Table 1 Top pathways and associated networks identified by pathway analysis for the hippocampus following olanzapine treatment
Full size table

Pathways and associated network functions of genes that had changes in methylation in the cerebellum following olanzapine treatment

We also analysed genes that had increased methylation in the cerebellum (Additional file 4: Table S2A and Additional file 5: Table S2B). The most significant pathway identified for the cerebellum was for ephrin receptor signalling (P = 5.23 × 10–4) (Table 2a; Figure 2). Synaptic long-term potentiation (P = 2.94 × 10–3), which is implicated in learning and plasticity, was among the most significant pathways identified (Additional file 1: Figure S4). Moreover, pathways involved in signalling (Erk/Mapk, circadian rhythm and protein kinase A) were identified (Table 2a). Interestingly, genes with reduced methylation were also involved in pathways such as ephrin B signalling (Table 2b).

Table 2 Top networks identified by pathway analysis for the cerebellum following olanzapine treatment
Full size table
Figure 2
figure2

Ephrin receptor signalling is the most significant canonical pathway in an olanzapine-treated cerebellum. An asterisk indicates a gene previously implicated in psychosis (from Ingenuity Pathway Analysis).

Full size image

Pathways and associated network functions of genes that had changes in methylation in the liver following olanzapine treatment

The results from a non-brain tissue sample also showed that the effects of olanzapine are not restricted to brain regions; it may also affect liver the (Additional file 6: Table S3A and Additional file 7: Table S3B). Genes that had an increase (Additional file 8: Table S4a) or a decrease (Additional file 8: Table S4b) in methylation in the liver, were in several pathways including lipid metabolism, cell death and organ morphology. Such genes include DRD1/DRD2, NMDAR and PTEN (Additional file 1: Figure S5).

Discussion

There are a number of methods to test the effects of antipsychotic drugs, including the locomotor activity test [27] and the prepulse inhibition test [28]. In this study, we used the locomotor activity test. Significantly reduced locomotor activity in olanzapine-treated rats in this experiment (Figure 1A,B) has suggested that the drug administration paradigm employed was sufficient to cause therapeutically relevant effects in rats. Comparable therapeutic doses in rats were effective in previous studies and resulted in locomotor-suppressive effects [27]. Moreover, the significant increase in weight of olanzapine-treated rats in this and previous studies [29] indicated that the paradigm adapted might also be capable of causing metabolic disturbances, as seen in patients taking olanzapine for a long time [30]. Interestingly, the molecular results showed that olanzapine treatment caused genome-wide DNA methylation changes (Figure 1). Further, the results showed that most genes affected were tissue specific (hippocampus, cerebellum or liver). Also, the gene-specific methylation changes affected a number of networks that were tissue-specific, as expected. More importantly, the identified networks support two known effects of olanzapine, discussed in the following sections. The first is the recovery from psychosis [31] and the second is the adverse effects [32] of olanzapine. Olanzapine-induced DNA methylation changes in genes involved in canonical pathways may alter the associated network functions. However, further study is required to analyse the effects (on a protein level) of, specifically, the gene-specific methylation changes on each identified network.

We argue that the two manifestations could be attributed to tissue-specific alterations that disturb the coordinated expression of genes critical in the identified networks (Tables 1 and 2). This model is backed by a number of observations. First, the phenotypic effect of olanzapine is not immediate; rather it takes days or weeks after the initiation of treatment [13]. This may be the time that is needed for gene-specific methylation to alter the expression of the specific genes [33, 34]. Also, patients may not respond to this drug, depending on their CYP 1A2 genotype, which can metabolize this drug [35] or acquire resistance. Patients may need to take a different type of antipsychotic drug [16]. During this time, a patient may be affected by metabolic disorders, weight gain and related adverse effects [30]. We will discuss the specific mechanisms of the effects of olanzapine in the following section.

Olanzapine-based psychosis recovery may involve changes in gene methylation

We argue that an increase or a decrease in methylation of specific gene promoters, following olanzapine treatment, may decrease or increase their transcriptional efficiency [36, 37], specifically for the hippocampus, which is one of the primary sites responsible for psychotic symptoms [15, 38, 39]. Further, the pattern of transcriptional efficiency may also be modulated by other factors such as chromatin structure and elongation efficiency [38, 40, 41]. We acknowledge that the prefrontal cortex and nucleus accumbens, which are also implicated in psychosis [40, 42, 43], may need to be investigated in future studies.

We note that in the hippocampus, dopamine-DARPP32 feedback in the cAMP signalling pathway (P = 1.6 × 10–3) was the most significant pathway identified. Neurons in the midbrain release dopamine, which modulates cAMP (cyclic adenosine 3,5-monophosphate) production by activating dopamine receptors [44]. These results suggest that the antipsychotic effects of olanzapine may involve alterations in gene-specific methylation that leads to the dysregulation of genes involved in dopamine-DARPP32 feedback in the cAMP signalling pathway (Figure 3). This includes several differentially methylated genes such as Drd1/5 and Nos1. The dopamine blockade leads to the progressive reduction of psychosis while its disturbance leads to psychosis [45]. All antipsychotics block post-synaptic D2 receptors [11]. A serotonin-dopamine antagonist was formulated following the synthesis of second-generation antipsychotics [11]. However, patients frequently fail to respond to one antipsychotic but respond to a different drug even if both block the D2 receptor [16]. Also, schizophrenia patients may partially respond to an antipsychotic or do not respond at all [17]. This may be due to several factors, and one possibility would be the delay in the onset of therapeutic actions partly or fully influenced by the downstream effects, such as altered transcription [14, 33]. As such, the differentially methylated genes involved in the dopamine-signalling pathway may stop or reduce transcription and gene expression [14, 21, 33]. In fact, decreased expression of DARPP32 in the prefrontal cortex has been reported in schizophrenia patients [43, 46, 47]. Also, DNA methylation differences have been observed in the dopamine D2 receptor gene within and between pairs of monozygotic twins discordant for psychoses [48] and there is an overwhelming evidence for the involvement of dopamine in psychosis including schizophrenia [19].

Figure 3
figure3

Dopamine-DARPP32 feedback in cAMP signalling is a significant canonical pathway in an olanzapine-treated hippocampus. An asterisk indicates a gene previously implicated in schizophrenia (from Ingenuity Pathway Analysis).

Full size image

Further studies on the effects of drugs may help to identify the genes and pathways that underlie psychosis. For example, a decreased expression of CDC42 was reported in the cerebral cortex of schizophrenic patients in post-mortem studies, and this has been implicated in defects in dendritic spines in cortical neurons in the patients [49]. CDC42 can reorganize septin fibre formation, which is thought to stabilize actin filaments needed for a normal spine shape and synaptic plasticity, as reviewed in Ide and Lewis [49]. However, cautious interpretation of the results is important because actual epigenetic changes in schizophrenic patients may represent changes in methylation status [50].

Our results show that olanzapine caused an increase or a decrease in the methylation of genes previously implicated in schizophrenia (Table 3), which may reflect the fact that olanzapine could alleviate psychiatric symptoms via mechanisms involving DNA methylation. Among the genes that decreased in methylation in the hippocampus is Map6, which is implicated in schizophrenia [51] and is involved in molecular transport, nervous system development and function (Additional file 1: Figure S6). This emphasizes that methylation may serve an intermediary role whose actual effect is realized through gene expression.

Table 3 Genes positively associated with schizophrenia where methylation changed in response to olanzapine
Full size table

Among the genes that showed an increase in methylation in hippocampus was Bdnf, which has been previously implicated in schizophrenia [52]. This corroborates previous findings that showed its regulatory role in the expression of the dopamine D3 receptor gene (DRD3) [53]. The relatively lower methylation and higher expression of BDNF was also observed in schizophrenia patients compared to healthy controls [54].

Our results suggested that the efficacy of olanzapine might be achieved by changes in gene-specific methylation of relevant genes that take part in psychosis-related pathways. A list of 123 such genes that increased in methylation in both brain regions is given in Additional file 9: Table S5. That methylation may serve an intermediary role in modulating gene expression is apparent in the cerebellum, which is dominated by a number of signalling pathways including ephrin receptors and synaptic long-term potentiation (Table 2). Ephrin ligands and receptors guide axons during neural development and regulate neuronal plasticity in adults [55, 56]. Specifically, ephrin plays an important role in the regulation of neuronal migration, which is essential for the development of the nervous system and the proper functioning of the brain [57]. Neuronal cells have ahigher variation in DNA methylation than non-neuronal cells, supporting the idea that the epigenetic status of neuronal cells changes in response to the environment in the brain [58]. Interestingly, DNA methylation was found to be highly heritable and significantly correlated with gene expression in the human brain [59].

Furthermore, the synaptic long-term potentiation (LTP) pathway was one of the top canonical pathways in the cerebellum. Synaptic activity can persistently modify the way a neuron reacts to subsequent inputs by affecting either its intrinsic excitability or its synaptic efficacy, which is enhanced during long-term potentiation [60]. Specifically, in a rat cerebellum, synaptic transmission and granule cell intrinsic excitability are enhanced during LTP [47]. LTP is a well-known model for synaptic plasticity and it is typically induced by high-frequency activation of NMDA receptors at glutamatergic synapses [61]. Such results allow us to postulate that the efficacy of olanzapine may be due to its effect on the regulation of dopamine- DARPP32 feedback in the cAMP signalling pathway in the hippocampus, via DNA methylation.

Further, an atypical antipsychotic induced a restrictive chromatin state in the present study and previous reports [62]. On the other hand, clozapine was found to induce MII1, a mediator of open chromatin [63]. A restrictive chromatin state through DNA methylation has been implicated in psychiatric disorders [64]. Also, olanzapine, unlike clozapine and sulpiride, did not activate brain DNA demethylation in mice [65]. Moreover, atypical antipsychotics might regulate the transcription and function of genes that are related to histone post-transcriptional modifications [62]. Therefore, the mechanisms of actions of olanzapine on the chromatin structure and on any epigenetic machinery need to be studied further.

Adverse consequences of olanzapine treatment may include aberrant methylation

Although the focus of this study was to assess the effect of olanzapine treatment in the hippocampus and cerebellum, we also used the liver as a non-brain tissue sample. In addition to the pathways implicated in schizophrenia and psychosis, olanzapine treatment also affected pathways for lipid metabolism, metabolic diseases and cell deathin the hippocampus; cardiovascular disease and cell signallingin the cerebellum and cardiovascular system function, cell death and survivalin the liver. These effects may reflect specific adverse consequences such as weight gain [29]. Interestingly, olanzapine has been shown to alterlipid metabolism [65]. Further, IL6, which is differentially methylated and is involved in regulating the lipid metabolismpathway, was reported to destabilize atherosclerotic plaques in mice [66]. The Jak/Stat signalling pathway, which is affected by olanzapine treatment, is also known to regulate how muscle mass is lost or gained, which is an essential factor in defining obesity [67]. Furthermore, the pathway for the cellular effects of sildenafil (Viagra), which was identified for the hippocampus, may result in sexual dysfunction, which is a common side effect of many antipsychotics [68]. This reflects the possibility that the epigenomic status of the genes involved in this pathway causes this dysfunction.

It is important to point out that the selected pathways discussed are the most significant pathways identified in this study. The study provides a novel insight into the potential mechanisms in the olanzapine-induced reduction of psychosis and the associated adverse effects. Antipsychotics were previously shown to have modulated promoter methylation and thereby gene expression [33, 34]. We show for the first time that the pathways affected are for the known effects of olanzapine. Also, the effect of this drug on such pathways may involve alterations in gene-specific methylation. Further, the pathways affected are diverse and tissue specific. Thus, the findings in this report support the involvement of epigenetic changes that are known to be reversible and influenced by genetic as well as environmental factors, in neural function. They offer an original insight into any future epigenetic research in psychiatric disorders and potential avenues for personalized medicine.

Conclusions

The known functions of affected genes suggest that the observed epigenetic changes may underlie the amelioration of symptoms and account for certain adverse effects including the metabolic syndrome. The results give a novel insight into the mechanism of action of olanzapine, therapeutic effects and the side effects of antipsychotics.

Methods

Animals

Rats

Adult male Sprague–Dawley rats,12 weeks old weighing 250 to 300 g, were purchased from Charles River, QC, Canada. Upon arrival, the rats were separated into individual cages and housed under controlled humidity and temperature on a 12-hour light/dark cycle (the lights were switched on at 7.00 am). They were fed standard rat chow (LabDiet) and tap water ad libitum. The Institutional Animal Care Committee of the University of Western Ontario approved all animal-related procedures used in this study following the Canadian and National Institute of Health Guides on animal experimentation.

Olanzapine treatment

Before the commencement of olanzapine treatment, the animals were weighed and divided into two treatment groups with comparable mean weights. They were habituated individually for 30 min to an automated open-field activity-monitoring chamber (San Diego Instruments, San Diego, CA, USA), and then subjected to 5 min of tail pinch stress. Their stress-induced locomotor activity was recorded for the next 30 min using the open-field activity chamber. Starting from the following day, the rats received injections of olanzapine (Zyprexa, Lilly, IN, USA; 2.5 mg/kg, intramuscular; n = 8) or vehicle (PBS; n = 8) between 1.30 pm and 3.00 pm daily for 19 days. Although antipsychotic drugs have been administered orally and intramuscularly in rodent studies, we chose the intramuscular route to ensure the rats consistently received the intended dose for the entire duration of therapy.

Phenotypic analysis and tissue collection

Then 48 hours after the last dose of olanzapine or the vehicle (to avoid olanzapine-induced sedation interfering with stress perception and activity), the rats were habituated to the same open-field activity monitoring chamber for 30 min, subjected to a similar 5-min tail pinch as done at the baseline, and monitored to verify whether stress-induced locomotor activity was reduced in olanzapine-treated rats compared to the vehicle-treated group. This paradigm has been widely used to study the therapeutic efficacy of antipsychotic drugs [27, 69]. Subsequently, 24 hours after completion of the stress-induced behavioural assessment (to minimize the potential effect of stress on expected molecular changes), each rat was sacrificed. The rats were decapitated without anaesthesia, brain tissues were micro-dissected promptly in ice-cold PBS and three random biopsy punches through the hippocampus, cerebellum and liver were obtained. These three biopsy punches were considered to be a single sample. The sample from each rat was kept separately and flash-frozen in liquid nitrogen. Genomic DNA was isolated from olanzapine-treated and vehicle-treated samples to analyse the genome-wide methylation using rat methylation arrays.

Assessing genome-wide DNA methylation by MeDIP-chip analysis

Array hybridization

Genomic DNA was isolated from each of the three tissues from two random control samples and two random olanzapine-treated samples. All methylated DNA immunoprecipitation (MeDIP), sample labelling, hybridization, and processing were performed by Arraystar Inc (Rockville, MD, USA). Briefly, isolated genomic DNA was sonicated to generate random fragments of 200 to 1,000 bp. For DNA labelling, the NimbleGen Dual-Color DNA Labeling Kit was used according to the manufacturer’s guideline detailed in the NimbleGen MeDIP-chip protocol (NimbleGen Systems, Inc, Madison, WI, USA). Microarrays were hybridized at 42°C for 16 to 20 h with Cy3/5 labelled DNA in NimbleGen hybridization buffer/hybridization component A in a hybridization chamber (Hybridization System, NimbleGen Systems, Inc, Madison, WI, USA). The methylated DNA was immunoprecipitated using Biomag™ magnetic beads coupled with mouse monoclonal antibodies against 5-methylcytidine. The total input and matched immunoprecipitated DNA were labelled with Cy3- and Cy5-labelled random 9-mers, respectively, and hybridized to NimbleGen RN34 Meth 3×720K CpG plus Promoter arrays. Scanning was performed with the Axon GenePix 4000B microarray scanner.

Data normalization and analysis

Raw data was extracted as pair files using the NimbleScan software (Roche NimbleGen Inc). Median-centring quantile normalization and linear smoothing was performed using the Bioconductor packages Ringo, limma, and MEDME. From the normalized log2 ratio data, a sliding-window peak-finding algorithm provided by NimbleScan v2.5 (Roche NimbleGenInc) was applied to find the enriched peaks with specified parameters (sliding window width: 750 bp; mini probes per peak: 2; P- value minimum cut-off: 2; maximum spacing between nearby probes within peak: 500 bp).

To compare differentially enriched regions between drug-exposed (E) and matched control (C) rats, the log2 ratios were averaged and then used to calculate M' for each probe:

M ' = Average log 2 MeDIPE / InputE - Average log 2 MeDIPC / InputC

The NimbleScan sliding-window peak-finding algorithm was run on this data to find the differential enrichment peaks (DEPs). The differential enrichment peaks, identified by the NimbleScan algorithm, were filtered according to the following criteria: (i) at least one of the two groups had the median value of log2 MeDIP/Input ≥ 0.3 and a median value of M' > 0 within the peak; (ii) at least half of the probes in a peak had amedian value of the coefficient of variability (CV) ≤ 0.8 for both groups.

Using an R script program, a hierarchical clustering analysis was completed. The probe data matrix was obtained using PeakScores from differentially methylated regions selected by DEP analysis. This analysis used PeakScore ≥ 2 to define the DEPs, which is equivalent to the average P ≤ 0.01, for all probes within the peak.

Pathway and bioinformatic analysis of array results

A venn diagram of the genes was used to assess the distribution of genes affected across tissue types [70]. The identified gene promoters with significant alterations to DNA methylation were then subjected to Ingenuity Pathway Analysis (Ingenuity System Inc, CA, USA) [71].

Abbreviations

Bp:

Base pair

CV:

Coefficient of variability

DEP:

Differential enrichment peak

HCP:

High CpG contents

ICP:

Intermediate CpG contents

LCP:

Low CpG contents

LTP:

Long-term potentiation

MeDIP:

Methylated DNA immunoprecipitation

PBS:

Phosphate-buffered saline.

References

  1. 1.

    van Os J, Kapur S: Schizophrenia. Lancet. 2009, 374: 635-644.

  2. 2.

    Fitzgerald PB, Benitez J, Daskalakis JZ, De Castella A, Kulkarni J: The treatment of recurring auditory hallucinations in schizophrenia with rTMS. World J Biol Psychiatry. 2006, 7: 119-122. 10.1080/15622970500474705.

  3. 3.

    Carlsson A, Lindqvist M: Effect of chlorpromazine or haloperidol on formation of 3methoxytyramine and normetanephrine in mouse brain. Acta Pharmacol Toxicol (Copenh). 1963, 20: 140-144.

  4. 4.

    Kane JM: The current status of neuroleptic therapy. J Clin Psychiatry. 1989, 50: 322-328.

  5. 5.

    Lieberman JA, Stroup TS, McEvoy JP, Swartz MS, Rosenheck RA, Perkins DO, Keefe RS, Davis SM, Davis CE, Lebowitz BD, Severe J, Hsiao JK: Effectiveness of antipsychotic drugs in patients with chronic schizophrenia. N Engl J Med. 2005, 353: 1209-1223. 10.1056/NEJMoa051688.

  6. 6.

    Bervoets C, Morrens M, Vansteelandt K, Kok F, de Patoul A, Halkin V, Pitsi D, Constant E, Peuskens J, Sabbe B: Effect of aripiprazole on verbal memory and fluency in schizophrenic patients: results from the ESCAPE study. NS Drugs. 2012, 26: 975-982.

  7. 7.

    Lett TA, Wallace TJ, Chowdhury NI, Tiwari AK, Kennedy JL, Muller DJ: Pharmacogenetics of antipsychotic-induced weight gain: review and clinical implications. Mol Psychiatry. 2012, 17: 242-266. 10.1038/mp.2011.109.

  8. 8.

    Muller DJ, Muglia P, Fortune T, Kennedy JL: Pharmacogeneticsof antipsychotic-induced weight gain. Pharmacol Res. 2004, 49: 309-312. 10.1016/j.phrs.2003.05.001.

  9. 9.

    Muller DJ, Kennedy JL: Genetics of antipsychotic treatment emergent weight gain in schizophrenia. Pharmacogenomics. 2006, 7: 863-887. 10.2217/14622416.7.6.863.

  10. 10.

    Lencz T, Malhotra AK: Pharmacogenetics of antipsychotic-induced side effects. Dialogues Clin Neurosci. 2009, 11: 405-411.

  11. 11.

    Jasovic-Gasic M, Vukovic O, Pantovic M, Cvetic T, Maric-Bojovic N: Antipsychotics – history of development and field of indication, new wine – old glasses. Psychiatr Danub. 2012, 24 (Suppl 3): S342-S344.

  12. 12.

    van Rossum JM: The significance of dopamine-receptor blockade for the mechanism of action of neuroleptic drugs. Arch Int Pharmacodyn Ther. 1966, 160: 492-494.

  13. 13.

    Manschreck TC, Boshes RA: The CATIE schizophrenia trial: results, impact, controversy. Harv Rev Psychiatry. 2007, 15: 245-258. 10.1080/10673220701679838.

  14. 14.

    Razin A, Kantor B: DNA methylation in epigenetic control of gene expression. Prog Mol Subcell Biol. 2005, 38: 151-167. 10.1007/3-540-27310-7_6.

  15. 15.

    Grace AA: Dopamine system dysregulation by the hippocampus: implications for the pathophysiology and treatment of schizophrenia. Neuropharmacology. 2012, 62: 1342-1348. 10.1016/j.neuropharm.2011.05.011.

  16. 16.

    Clark SL, Adkins DE, van den Oord EJ: Analysis of efficacy and side effects in CATIE demonstrates drug response subgroups and potential for personalized medicine. Schizophr Res. 2011, 132: 114-120. 10.1016/j.schres.2011.07.031.

  17. 17.

    Dratcu L, Grandison A, McKay G, Bamidele A, Vasudevan V: Clozapine-resistant psychosis, smoking, and caffeine: managing the neglected effects of substances that our patients consume every day. Am J Ther. 2007, 14: 314-318. 10.1097/01.pap.0000249958.96498.ce.

  18. 18.

    Abi-Dargham A, Martinez D, Mawlawi O, Simpson N, Hwang DR, Slifstein M, Anjilvel S, Pidcock J, Guo NN, Lombardo I, Mann JJ, Van Heertum R, Foged C, Halldin C, Laruelle M: Measurement of striatal and extrastriatal dopamine D1 receptor binding potential with [11C]NNC 112 in humans: Validation and reproducibility. J Cereb Blood Flow Metab. 2000, 20: 225-243.

  19. 19.

    Howes OD, Kapur S: The dopamine hypothesis of schizophrenia: version III – the final common pathway. Schizophr Bull. 2009, 35: 549-562. 10.1093/schbul/sbp006.

  20. 20.

    Kelkar A, Deobagkar D: A novel method to assess the full genome methylation profile using monoclonal antibody combined with the high throughput based microarray approach. Epigenetics. 2009, 4: 415-420. 10.4161/epi.4.6.9768.

  21. 21.

    Razin A, Cedar H: DNA methylation and gene expression. Microbiol Rev. 1991, 55: 451-822.

  22. 22.

    Bredberg A, Bodmer W: Cytostatic drug treatment causes seeding of gene promoter methylation. Eur J Cancer. 2007, 43: 947-954. 10.1016/j.ejca.2006.12.003.

  23. 23.

    Szyf M: Epigenetic therapeutics in autoimmune disease. Clin Rev Allergy Immunol. 2010, 39: 62-77. 10.1007/s12016-009-8172-8.

  24. 24.

    Melas PA, Rogdaki M, Osby U, Schalling M, Lavebratt C, Ekstrom TJ: Epigenetic aberrations in leukocytes of patients with schizophrenia: association of global DNA methylation with antipsychotic drug treatment and disease onset. FASEB J. 2012, 26: 2712-2718. 10.1096/fj.11-202069.

  25. 25.

    Kunii Y, Hyde TM, Ye T, Li C, Kolachana B, Dickinson D, Weinberger DR, Kleinman JE, Lipska BK: Revisiting DARPP-32 in postmortem human brain: changes in schizophrenia and bipolar disorder and genetic associations with t-DARPP-32 expression. Mol Psychiatry. 2013, doi:10.1038/mp.2012.174

  26. 26.

    Bromberg A, Levine J, Nemetz B, Belmaker RH, Agam G: No association between global leukocyte DNA methylation and homocysteine levels in schizophrenia patients. Schizophr Res. 2008, 101: 50-57. 10.1016/j.schres.2008.01.009.

  27. 27.

    Bardgett ME, Humphrey WM, Csernansky JG: The effects of excitotoxic hippocampal lesions in rats on risperidone- and olanzapine-induced locomotor suppression. Neuro Psychopharmacol. 2002, 27: 930-938.

  28. 28.

    Geyer MA, Krebs-Thomson K, Braf DL, Swerdlow NR: Pharmacological studies of prepulse inhibition models of sensorimotor gating deficits in schizophrenia: a decade in review. Psychopharmacology (Berl). 2001, 156: 117-154. 10.1007/s002130100811.

  29. 29.

    Deng C, Lian J, Pai N, Huang XF: Reducing olanzapine-induced weight gain side effect by using betahistine: a study in the rat model. J Psychopharmacol. 2012, 26: 1271-1279. 10.1177/0269881112449396.

  30. 30.

    Schreiner A, Niehaus D, Shuriquie NA, Aadamsoo K, Korcsog P, Salinas R, Theodoropoulou P, Fernández LG, Uçok A, Tessier C, Bergmans P, Hoeben D: Metabolic effects of paliperidone extended release versus oral olanzapine in patients with schizophrenia: a prospective, randomized, controlled trial. J Clin Psychopharmacol. 2012, 32: 449-457. 10.1097/JCP.0b013e31825cccad.

  31. 31.

    Shrivastava A, Johnston M, Terpstra K, Stitt L, Shah N: Atypical antipsychotics usage in long-term follow-up of first episode schizophrenia. Indian J Psychiatry. 2012, 54: 248-252. 10.4103/0019-5545.102425.

  32. 32.

    Gray R: Olanzapine: efficacy in treating the positive and negative symptoms of schizophrenia. Ment Health Care. 1998, 1: 193-194.

  33. 33.

    Feng W, Dong Z, He B, Wang K: Analysis method of epigenetic DNA methylation to dynamically investigate the functional activity of transcription factors in gene expression. BMC Genomics. 2012, 13: 532. 10.1186/1471-2164-13-532.

  34. 34.

    Dong E, Grayson DR, Guidotti A, Costa E: Antipsychotic subtypes can be characterized by differences in their ability to modify GABAergic promoter methylation. Epigenomics. 2009, 1: 201-211. 10.2217/epi.09.2.

  35. 35.

    Ereshefsky L, Riesenman C, Lam YW: Serotonin selective reuptake inhibitor drug interactions and the cytochrome P450 system. J Clin Psychiatry. 1996, 57 (Suppl 8): 17-24.

  36. 36.

    Gray SG, Eriksson T, Ekstrom TJ: Methylation, gene expression and the chromatin connection in cancer (review). Int J Mol Med. 1999, 4: 333-350.

  37. 37.

    Numata S, Ye T, Hyde TM, Guitart-Navarro X, Tao R, Wininger M, Colantuoni C, Weinberger DR, Kleinman JE, Lipska BK: DNA methylation signatures in development and aging of the human prefrontal cortex. Am J Hum Genet. 2012, 90: 260-272. 10.1016/j.ajhg.2011.12.020.

  38. 38.

    Jun H, HussainiS MQ, Rigby MJ, Jang MH: Functional role of adult hippocampal neurogenesis as a therapeutic strategy for mental disorders. Neural Plast. 2012, 2012: 854285-

  39. 39.

    Eisenberg DP, Ianni AM, Wei SM, Kohn PD, Kolachana B, Apud J, Weinberger DR, Berman KF: Brain-derived neurotrophic factor (BDNF) val(66)met polymorphism differentially predicts hippocampal function in medication-free patients with schizophrenia. J Neuropsych Clin Neurosci. 2009, 21: 30-37. 10.1176/appi.neuropsych.21.1.30.

  40. 40.

    Lorincz MC, Dickerson DR, Schmitt M, Groudine M: Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells. Nat Struct Mol Biol. 2004, 11: 1068-1075. 10.1038/nsmb840.

  41. 41.

    Oh YS, Gao P, Lee KW, Ceglia I, Seo JS, Zhang X, Ahn JH, Chait BT, Patel DJ, Kim Y, Greengard P: SMARCA3, a chromatin-remodeling factor, is required for p11-dependent antidepressant action. Cell. 2013, 152: 831-843. 10.1016/j.cell.2013.01.014.

  42. 42.

    Santini MA, Ratner C, Aznar S, Klein AB, Knudsen GM, Mikkelsen JD: Enhanced prefrontal serotonin 2A receptor signaling in the subchronic phencyclidine mouse model of schizophrenia. J Neurosci Res. 2013, 91 (5): 634-641. 10.1002/jnr.23198.

  43. 43.

    Olszewski M, Piasecka J, Goda SA, Kasicki S, Hunt MJ: Antipsychotic compounds differentially modulate high-frequency oscillations in the rat nucleus accumbens: a comparison of first- and second-generation drugs. Int J Neuropsychopharmacol. 2012, 16: 1009-20.

  44. 44.

    Svenningson P, Nishi A, Fisone G, Girault J, Nairn AC, Greengard P: DARPP-32: an integrator of neurotransmission. Annu Rev Pharmacol Toxicol. 2004, 44: 269-296. 10.1146/annurev.pharmtox.44.101802.121415.

  45. 45.

    Ginovart N, Kapur S: Role of dopamine D(2) receptors for antipsychotic activity. Handb Exp Pharmacol. 2012, 212: 27-52. 10.1007/978-3-642-25761-2_2.

  46. 46.

    Albert KA, Hemmings HC, Adamo AI, Potkin SG, Akbarian S, Sandman CA, Cotman CW, Bunney WE, Greengard P: Evidence for decreased DARPP-32 in the prefrontal cortex of patients with schizophrenia. Arch Gen Psychiatry. 2002, 59: 705-712. 10.1001/archpsyc.59.8.705.

  47. 47.

    Li C, Zheng Y, Qin W, Tao R, Pan Y, Xu Y, Li X, Gu N, Feng G, He L: A family-based association study of kinesin heavy chain member 2 gene (KIF2) and schizophrenia. NeurosciLett. 2006, 407: 151-155.

  48. 48.

    Petronis A, Gottesman II, Kan P, Kennedy JL, Basile VS, Paterson AD, Popendikyte V: Monozygotic twins exhibit numerous epigenetic differences: clues to twin discordance?. Schizophr Bull. 2003, 29: 169-178. 10.1093/oxfordjournals.schbul.a006988.

  49. 49.

    Ide M, Lewis DA: Altered cortical CDC42 signaling pathways in schizophrenia: implications for dendritic spine deficits. Biol Psychiatry. 2010, 68: 25-32. 10.1016/j.biopsych.2010.02.016.

  50. 50.

    Tsankova NM, Berton O, Renthal W, Kumar A, Neve RL, Nestler EJ: Sustained hippocampal chromatin regulation in a mouse model of depression and antidepressant action. Nat Neurosci. 2006, 9: 519-525. 10.1038/nn1659.

  51. 51.

    Shimizu H, Iwayama Y, Yamada K, Toyota T, Minabe Y, Nakamura K, Nakajima M, Hattori E, Mori N, Osumi N, Yoshikawa T: Genetic and expression analyses of the STOP (MAP6) gene in schizophrenia. Schizophr Res. 2006, 84: 244-252. 10.1016/j.schres.2006.03.017.

  52. 52.

    Ho BC, Milev P, O’Leary DS, Librant A, Andreasen NC, Wassink TH: Cognitive and magnetic resonance imaging brain morphometric correlates of brain-derived neurotrophic factor Val66Met gene polymorphism in patients with schizophrenia and healthy volunteers. Arch Gen Psychiatry. 2006, 63: 731-740. 10.1001/archpsyc.63.7.731.

  53. 53.

    Maziade M, Martinez M, Rodrigue C, Gauthier B, Tremblay G, Fournier C, Bissonnette L, Simard C, Roy MA, Rouillard E, Mérette C: Childhood/early adolescence-onset and adult-onset schizophrenia. Heterogeneity at the dopamine D3 receptor gene. Br J Psychiatry. 1997, 170: 27-30. 10.1192/bjp.170.1.27.

  54. 54.

    Kordi-Tamandani DM, Sahranavard R, Torkamanzehi A: DNA methylation and expression profiles of the brain-derived neurotrophic factor (BDNF) and dopamine transporter (DAT1) genes in patients with schizophrenia. Mol Biol Rep. 2012, 39: 10889-10893. 10.1007/s11033-012-1986-0.

  55. 55.

    Klein R: Bidirectional modulation of synaptic functions by eph/ephrin signaling. Nat Neurosci. 2009, 12: 15-20. 10.1038/nn.2231.

  56. 56.

    Senturk A, Pfennig S, Weiss A, Burk K, Acker-Palmer A: EphrinBsare essential components of the Reelin pathway to regulate neuronal migration. Nature. 2011, 472: 356-360. 10.1038/nature09874.

  57. 57.

    Iwamoto S, Suzuki T, Sutani A, Kuraki T, Isobe T: A case of atypical drug-induced hypersensitivity syndrome caused by isoniazid. Kekkaku. 2012, 87: 777-782.

  58. 58.

    Zhang D, Cheng L, Badner JA, Chen C, Chen Q, Luo W, Craig DW, Redman M, Gershon ES, Liu C: Genetic control of individual differences in gene-specific methylation in human brain. Am J Hum Genet. 2010, 86: 411-419. 10.1016/j.ajhg.2010.02.005.

  59. 59.

    Armano S, Rossi P, Taglietti V, D’Angelo E: Long-term potentiation of intrinsic excitability at the mossy fiber-granule cell synapse of rat cerebellum. J Neurosci. 2000, 20: 5208-5216.

  60. 60.

    Aizenman CD, Linden DJ: Rapid, synaptically driven increases in the intrinsic excitability of cerebellar deep nuclear neurons. Nat Neurosci. 2000, 3: 109-111. 10.1038/72049.

  61. 61.

    Kurita M, Holloway T, García-Bea A, Kozlenkov A, Friedman AK, Moreno JL, Heshmati M, Golden SA, Kennedy PJ, Takahashi N, Dietz DM, Mocci G, Gabilondo AM, Hanks J, Umali A, Callado LF, Gallitano AL, Neve RL, Shen L, Buxbaum JD, Han MH, Nestler EJ, Meana JJ, Russo SJ, González-Maeso J: HDAC2 regulates atypical antipsychotic responses through the modulation of mGlu2 promoter activity. Nat Neurosci. 2012, 15: 1245-1254. 10.1038/nn.3181.

  62. 62.

    Huang Y, Kang BN, Tian J, Liu Y, Luo HR, Hester L, Snyder SH: The cationic amino acid transporters CAT1 and CAT3 mediate NMDA receptor activation-dependent changes in elaboration of neuronal processes via the mammalian target of rapamycinm TOR pathway. J Neurosci. 2007, 27: 449-458. 10.1523/JNEUROSCI.4489-06.2007.

  63. 63.

    Abdolmaleky HM, Abdolmaleky HM, Smith CL, Faraone SV, Shafa R, Stone W, Glatt SJ, Tsuang MT: Methylomics in psychiatry: modulation of gene–environment interactions may be through DNA methylation. Am J Med Genet B Neuropsychiatr Genet. 2004, 127B: 51-59. 10.1002/ajmg.b.20142.

  64. 64.

    Dong E, Nelson M, Grayson DR, Costa E, Guidotti A: Clozapine and sulpiride but not haloperidol or olanzapine activate brain DNA demethylation. Proc Natl Acad Sci USA. 2008, 105: 13614-13619. 10.1073/pnas.0805493105.

  65. 65.

    Minet-Ringuet J, Even PC, Valet P, Carpéné C, Visentin V, Prévot D, Daviaud D, Quignard-Boulange A, Tomé D, de Beaurepaire R: Alterations of lipid metabolism and gene expression in rat adipocytes during chronic olanzapine treatment. Mol Psychiatry. 2007, 12: 562-571. 10.1038/sj.mp.4001948.

  66. 66.

    Zhang K, Huang XZ, Li XN, Feng M, Li L, Cai XJ: Interleukin 6 destabilizes atherosclerotic plaques by downregulating prolyl-4-hydroxylase alpha1 via a mitogen-activated protein kinase and c-jun pathway. Arch Biochem Biophys. 2012, 528: 127-133. 10.1016/j.abb.2012.09.007.

  67. 67.

    Pijet M, Pijet B, Litwiniuk A, Pajak B, Gajkowska B, Orzechowski A: Leptin impairs myogenesis in C2C12 cells through JAK/STAT and MEK signaling pathways. Cytokine. 2013, 61: 445-54. 10.1016/j.cyto.2012.11.002.

  68. 68.

    Mir A, Shivakumar K, Williamson RJ, McAllister V, O’Keane V, Aitchison KJ: Change in sexual dysfunction with aripiprazole: a switching or add-on study. J Psychopharmacol. 2008, 22: 244-253. 10.1177/0269881107082901.

  69. 69.

    Scornaiencki R, Cantrup R, Rushlow WJ, Rajakumar N: Prefrontal cortical D1 dopamine receptors modulate subcortical D2 dopamine receptor-mediated stress responsiveness. Int J Neuropsychopharmacol. 2009, 12: 1195-2208. 10.1017/S1461145709000121.

  70. 70.

    Pirooznia M, Nagarajan V, Deng Y: GeneVenn - a web application for comparing gene lists using Venn diagrams. Bioinformation. 2007, 1: 420-422. 10.6026/97320630001420.

  71. 71.

    Ingenuity pathway analysis.http://www.ingenuity.com/,

Download references

Acknowledgements

This research was supported by grants from the Canadian Institute of Health Research (SMS) and the Ontario Mental Health Foundation (NR). SMS held a Senior Research Fellowship from the Ontario Mental Health Foundation during the course of this study.

Author information

Affiliations

  1. Molecular Genetics Unit, Department of Biology, The University of Western Ontario, London, ON, N6A 5B7, Canada

    • Melkaye G Melka
    • , Benjamin I Laufer
    • , Patrick McDonald
    • , Christina A Castellani
    •  & Shiva M Singh

  2. Department of Psychiatry, The University of Western Ontario, London, ON, N6A 5B7, Canada

    • Nagalingam Rajakumar
    •  & Richard O’Reilly

  3. Department of Neuroscience, The University of Western Ontario, London, ON, N6A 5B7, Canada

    • Shiva M Singh

  4. Children's Health Research Institute, 800 Commissioners Road East, London, ON, N6C 2V5, Canada

    • Shiva M Singh

Authors

  1. Search for Melkaye G Melka in:

  2. Search for Benjamin I Laufer in:

  3. Search for Patrick McDonald in:

  4. Search for Christina A Castellani in:

  5. Search for Nagalingam Rajakumar in:

  6. Search for Richard O’Reilly in:

  7. Search for Shiva M Singh in:

Corresponding author

Correspondence to Shiva M Singh.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MGM analyze the data and wrote the first draft of the manuscript. BIL helped analyze the data and reviewed the manuscript. PM worked on genomic DNA extraction from the rat brain regions and liver. CAC critically reviewed the manuscript for important intellectual content. NR prepared the animals and provided animal tissues. NR and RO participated in the design of the study and also critically reviewed the manuscript for important intellectual content. SMS conceived of the study and critically reviewed and approved the manuscript. All authors have approved the final version of the manuscript.

Electronic supplementary material

Additional file 1:Supplementary figures.(DOC 2 MB)

Additional file 2: Table S1A: Genes with increased methylation in the hippocampus as a result of olanzapine treatment. (XLS 178 KB)

Additional file 3: Table S1B: Genes with decreased methylation in the hippocampus as a result of olanzapine treatment. (XLS 110 KB)

Additional file 4: Table S2A: Genes with increased methylation in the cerebellum as a result of olanzapine treatment. (XLS 199 KB)

Additional file 5: Table S2B: Genes with decreased methylation in thecerebellum as a result of olanzapine treatment. (XLS 100 KB)

Additional file 6: Table S3A: Genes with increased methylation in the liver as a result of olanzapine treatment. (XLS 202 KB)

Additional file 7: Table S3B: Genes with decreased methylation in the liver as a result of olanzapine treatment. (XLS 96 KB)

Additional file 8: Table S4: Most significant networks identified by pathway analysis for genes that had (a) an increase or (b) a decrease in methylation, in the liver following olanzapine treatment. (DOC 46 KB)

Additional file 9: Table S5: Genes with increased methylation in both the hippocampus and cerebellum following olanzapine treatment. (XLS 24 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Keywords

  • Epigenetic
  • Hippocampus
  • Cerebellum
  • DNA methylation
  • Psychosis
  • Rat

Clinical Epigenetics

ISSN: 1868-7083

Contact us