Transgenerational deep sequencing revealed hypermethylation of hippocampal mGluR1 gene with altered mRNA expression of mGluR5 and mGluR3 associated with behavioral changes in Sprague Dawley rats with history of prolonged febrile seizure

Authors: Oluwole Ojo Alese ^aff001; Musa V. Mabandla ^aff001
Authors place of work: Department of Human Physiology, College of Health Sciences, University of Kwazulu-Natal, Durban, South Africa ^aff001
Published in the journal: PLoS ONE 14(11)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0225034

Summary

The impact of febrile seizure has been shown to transcend immediate generation with the alteration of glutamatergic pathway being implicated. However, transgenerational effects of this neurological disorder particularly prolonged febrile seizure (PFS) on neurobehavioral study and methylation profile is unknown. We therefore hypothesized that transgenerational impact of prolonged febrile seizure is dependent on methylation of hippocampal mGluR1 gene. Prolonged febrile seizure was induced on post-natal day (PND) 14, by injecting lipopolysaccharide (LPS; 217μg/kg ip) and kainic acid (KA; 1.83 mg/kg ip). Sucrose preference test (SPT) and Forced swim test (FST) were carried out in the first generation (F₀) of animals at PND37 and PND60. The F₀ rats were decapitated at PND 14, 37 and 60 which corresponded to childhood, adolescent and adulthood respectively and their hippocampal tissue collected. The second generation (F₁) rats were obtained by mating F₀ generation at PND 60 across different groups, F₁ rats were subjected to SPT and FST test on PND 37 only. Decapitation of F₁rats and collection of hippocampal tissues were done on PND 14 and 37. Assessment of mGluR5 and mGluR3 mRNA was done with PCR while mGluR1 methylation profile was assessed with the Quantitative MassARRAY analysis. Results showed that PFS significantly leads to decreased sucrose consumption in the SPT and increased immobility time in the FST in both generations of rats. It also leads to significant decrease in mGluR5 mRNA expression with a resultant increased expression of mGluR3 mRNA expression and hypermethylation of mGluR1 gene across both generations of rats. This study suggested that PFS led to behavioral changes which could be transmitted on to the next generation in rats.

Keywords:

DNA methylation – Epigenetics – Methylation – Rats – Glutamate – Swimming – Hippocampus – Sucrose

Introduction

Febrile seizure is the most prevalent seizure type, affecting 3% to 5% of infants between the ages of 3months and 5years [1,2]. It is caused by infections (including respiratory tract infection, gastroenteritis and otitis media) that trigger the immune system leading to inflammatory responses with a subsequent release of macrophages, neutrophils and pro-inflammatory cytokines [3,4]. Sequel to these, the infections cause systemic fever that subsequently leads to aberrant neuronal excitability in febrile seizure [5,6,7]. However, one important factor that lead to the characterized aberrant neuronal excitability in febrile seizure, is neurotransmitter imbalance between glutamate and γ aminobutyric acid (GABA) [8].

Glutamate is the predominant excitatory neurotransmitter in the brain that has been associated with aberrant signaling and pathology of febrile seizure [7]. Furthermore, glutamate mediates its excitatory effects through two receptors that are ionotropic and metabotropic. Ionotropic glutamate receptors are dependent on ionic gradients while metabotropic glutamate receptors facilitate their action by coupling to secondary G-protein messengers [9]. Among the metabotropic glutamate receptors, the activation of metabotropic glutamate receptor 1 and 5 (mGluR1, mGluR5) enhances presynaptic glutamate release, thus modulating glutamatergic signaling in seizure [7]. Also, the metabotropic glutamate receptor 3 (mGluR3) regulates glutamate release on the presynaptic membrane and glia cell [10]. More importantly, previous studies have revealed up-regulation of the expression of mGluR1 on astrocytes in epilepsy and seizure [7,11,12]. Likewise hypermethylation of mGluR1 gene has been reported in malignant melanoma in mouse [13] and in the hippocampus of prenatal stressed rats offspring with depression [14,15].

Previous studies showed that different genetic alterations may be linked to abnormal behavior and functioning of the developing brain [16,17]. Behavioral and emotional changes have been linked to an early life exposure to unfavorable environmental factors which can be transferred to the offspring [18,19,20]. It has been demonstrated that stress in childhood may result in emotional disturbance in adulthood which could also be transferred across generations [21]. Early childhood PFS exposure in rats has been shown to cause memory deficit in adulthood that can be epigenetically transferred to their offspring [22]. Epigenetic alteration has been attributed to these environmental variations and as such, transferrable to their offspring [23]. Besides, epigenetics is the chromatin network modifications which lead to alteration of gene expression leaving the DNA sequence intact [24]. Epigenetic modifications have been reported to affect regulatory processes involved in learning and memory, synaptic plasticity and central nervous system development [25]. However, such epigenetic modification of mGluR1 in prolonged febrile seizure even across generation of rats is yet to be determined. Hence, we designed the present study to investigate the transgenerational effects of PFS through hippocampal mGluR5and mGluR3 mRNA expression. Also, we examined the transgenerational effects of PFS on neurobehavioral assessments and the methylation profile of mGluR1 gene in Sprague Dawley rats.

Materials and methods

Animals

All experimental procedures and animals used were treated in accordance with the approved guideline by the Animal Ethics Research Committee of the University of KwaZulu-Natal (AREC/042/016D).The adult rats were obtained from the Biomedical Resource Center of the University of KwaZulu-Natal, where they were housed under standard laboratory conditions with a 12 h light/dark cycle (lights on at 06h: 00). The rats had access to food and water ad libitum. A total of 64 pups (56 males and 8 females) obtained from mating female and male adult Sprague- Dawley rats were used in the experiments. The sample size was set according to previous studies where the statistical power was shown [26].

Experimental procedure

On post-natal day (PND) 14, prolonged febrile seizures (PFS) was induced in 32 pups (28 males and 4 females) by injecting lipopolysaccharide (LPS, 217 μg/kg) intraperitoneally (i.p) followed by an injection of kainic acid (KA, 1.83 mg/kg i.p) at 2h 30min later [27]. LPS had been demonstrated to reliably increase body temperature [28]. Temperature was monitored with the use of a digital thermometer before LPS injection and thereafter every 30 min until rats were given KA [29]. The remaining animals were grouped as control and injected with a corresponding volume of normal saline (SAL). Animals were monitored for one hour following the injection of KA and seizure severity was observed as follows: stage 0: no change in behavior; stage 1: chewing; stage 2: gazing and head nodding; stage 3: unilateral forelimb clonus, twitching and scratching; stage 4: rearing with bilateral forelimb clonus; stage 5: widespread muscle spasms, rearing with bilateral forelimb clonus and falling back [30]. All rats induced had at least stage 4 convulsions for a minimum of 30min duration, indicating that the convulsion were visibly obvious to the observer as rearing with bilateral forelimb clonus which is consistent with earlier observation [30]. In the first generation (F₀), 8 males (4 PFS and 4 SAL) rats were decapitated at PND 14 1h after the induction of PFS. The remaining pups were then returned to their mothers and weaned on PND 21.Two sets of 20 male rats from the F₀ generation (10 PFS and 10 SAL in each set) were decapitated after the neurobehavioral test at PND 37 and PND 60 respectively. This timelines of PND 14, 37 and 60 corresponds to the life span of the Sprague-Dawley rats from childhood, through adolescent to adulthood respectively. On PND 60, the second generation (F₁) rats were obtained by mating the remaining 16 rats (8males and 8 females) according to Table 1 below. A total of 16 rats comprising of 4offsprings from each mating group were sacrificed on PND 14 across the various groups while 28 rats were subjected to neurobehavioral test and decapitated on PND 37 across the various groups in F₁.

Table showing the first generation (F<sub>0</sub>) animal grouping for mating. — **Tab. 1. Table showing the first generation (F₀) animal grouping for mating.**

Sucrose preference test (SPT)

The SPT was conducted over a 24 h period to check for anhedonia in the animals [31]. A day prior to the test, (PND 34 and 57)in F₀ rats and PND 34 in F₁ generation, the rats were weighed and placed in separate cages for a 24 h training period. The training phase consists of placing two bottles of known volume of tap water on opposite sides of the cage [32]. After 24 h (PND 35 and 58), one bottle of water was replaced with another containing 5% sucrose solution. For the test phase, the two bottles (one containing tap water and the other 5% sucrose solution) were placed on the left and right sides respectively. After 24 h, the bottles were weighed to determine consumption in grams (converted to ml). The sucrose preference was calculated as the volume of sucrose drank over the total volume of fluid consumed (sucrose + water) x 100. A decreased amount of sucrose solution consumed is suggestive of anhedonia hence depressive-like behavior [31].

Forced swim test (FST)

The forced swim test assesses the despair behaviour in animals [33]. On PND 36 and 59 in F₀ rats and PND 36 in F₁ rats, each rat was placed in a plexiglass cylinder (50 H X30 D in cm) filled with water (24 ± 0.5 ^oC) to a level of 30 cm for 15 min as a pre-test [33]. This pre-test is required to enable quick adoption of an immobile posture by rats on the test day hence an obvious despair behaviour which is expected in depressed rats. The test session (5 min) was conducted on PND 37 and 60in F₀ rats and PND 37 in F₁ rats. The immobile time that the rats spent by floating on water with the head slightly immersed and/or minor movements to maintain the head above water after the 1^st min was video-recorded and documented for each animal [34].

Decapitation

All rats were decapitated using a sharp guillotine. Hippocampal tissues were immediately dissected, weighed and placed in Eppendorf tubes before snap freezing in liquid nitrogen. The tissues were stored in a -80°C bio-freezer until the day of analysis.

Real-time PCR

Twenty (20) micrograms of hippocampal tissue were weighed, homogenized and suspended in 300 μl of RNA lysis buffer (Zymo Research, USA). Thereafter, RNA isolation was carried out as instructed in manufacturer’s protocol (ZR RNA MiniPrep^TM, USA). Purification of RNA isolates was assessed using a NanoDrop (Thermo Scientific USA). Purity of 1.8–2.01 was recommended for use in the construction of cDNA. The cDNA synthesis was carried out using the iScript^TMcDNA Synthesis Kit (BioRad, South Africa) and run through the Thermocycler according to instructions stipulated in the protocol. The Fast start SYBR green kit (Roche Diagnostics, USA) was used in accordance with the manufacturer’s guidelines. The mGluR5 (Grm5) and mGluR3 (Grm3) primers used to assay genes were designed by Inqaba Biotech (Pretoria, South Africa) as shown in Table 2. Primer sequences were reconstituted in RNA nuclease free water according to the manufacturer’s protocol and were added to a master mix comprising of SYBR green dye, nuclease free H₂O and MgCl₂. Thereafter, cDNA was added to the strips and run in the Light cycler 480 (Roche USA) at optimized conditions as instructed in the manufacturer’s protocol. Experiments were performed in triplicate, GAPDH gene was used as the house keeping control and Data were analyzed by comparing C(q) values of the PFS samples to SAL and both normalized to GAPDH using the 2-ΔΔCq comparative method [35]

Quantitative MassARRAY analysis of gene methylation status

DNA extraction

Genomic DNA was extracted from the hippocampal tissue using a DNeasy® Blood and Tissue Kit according to the manufacturer’s instructions (Qiagen, Hilden, Germany). Twenty-five mg of hippocampal tissue of each animal were weighed and placed in a microcentrifuge tube. Then, 180μl of buffer ATL and 20μl of proteinase K was added into the tubes, incubated in a water bath at 56 ^oC and vortexed intermittently until the tissue were lysed completely. Subsequently, 200 μl of AL buffer was added and incubated at 56 ^oC for 10min. Thereafter, 200 μl of 100% ethanol was added and mixed thoroughly. The mixture was then pipetted into the DNeasy Mini spin column placed in 2ml collection tube and centrifuged at 8000 rpm for 1 min (HERMLE Labortechnik GmbH, Germany), the flow-through and the collection tubes were discarded. The Mini spin column was again placed into a new 2 ml collection tube with 500 μl of buffer AW1 added via the column and centrifuge at 8000 rpm for 1min, the flow-through and the collection tubes were discarded. Subsequently, spin column was placed in a new collection tube,500 μl of buffer AW2was added and the mixture was spun for 3min at 14000 rpm with the flow-through and the collection tubes discarded. To elute the DNA, the spin column was placed in 2 ml collection tube and 200 μl of buffer AE was added, incubated at room temperature for 1min and centrifuged for 1min at 8000rpm.This elution step was repeated to enhance the yield. The concentration and purity of the DNA were determined by absorbances at 260–280 nm with a ratio 1.7–2.0 by NanoDropTM 1000 spectrophotometer (Thermo Scientific, Wilmington, USA).

DNA sodium bisulfite conversion

Sodium bisulfate modification of the extracted DNA was performed using an EZ DNA Methylation Kit^TM according to the manufacturer’s instructions (Zymo Research, Orange, CA, USA). Five microliter of M-Dilution buffer was added to each DNA sample of 500ng, adjusted to a total volume to 50μl with distilled water and incubated at 37°C for 15min. Then, 100μl of CT conversion reagent was added to each sample and mixed properly. The samples were incubated in the thermocycler for 20 cycles at 95°C for 30s and 50°C for 15 min, followed by a final holding step at 4°C.A volume of 400 μl of M-Binding buffer was added to a Zymo-Spin^™ IC Column placed in a collection tube, this was then followed by loading the sample from the previous preparation above and mixing thoroughly by inverting the columns severally. The preparation was centrifuged at 8000rpm for 30s and the flow-through discarded, this was followed by adding 200μl of M-Wash buffer to the column and then centrifuged at full speed for 30s. Subsequently, 200μl of M-Desulphonation buffer was added to the column and incubated at room temperature (20–30°C) for 20 min. It was there after centrifuged at full speed for 30s;200μl of M-Wash buffer was added to the column and centrifuged at 14000rpm for 30s. The above step was repeated and the columns were placed into a 1.5ml microcentrifuge tube. Furthermore, 20 μl of pre-warm M-Elution buffer was added directly to the column matrix, incubated at room temperature for 1 min and then centrifuged for 30s at full speed to elute the DNA.

Quantitative mass ARRAY analysis of gene methylation status

The Sequenom MassARRAY platform was used for the quantitative methylation analysis of the upstream sequence of mGluR1 gene promoter. The region analysed and the CpG sites of mGluR1 promoter are shown in Fig 1. The primers used in this study were designed by Inqaba Biotech (Pretoria, South Africa) (forward primer AGGAAGAGAGTTGTTAGGTATTTTGGGTAAAATGG and reverse primer CAGTAATACGACTCACTATAGGGAGAAGGCTAAACCCAAAAATTTAAATACAATTCC). The PCR mixture (10 ng bisulfite-treated DNA, 25 mMdNTP, 5U/μL of PCR enzyme and a 1μM mixture of forward and reverse primers) was pre-heated for 4 min at 94°C and then incubated for 45 cycles of 94°C for 20s, 56°C for 30 s, and 72°C for 60s, followed by 72°C for 3 min. Two microliters of SAP mix containing 1.7μlH₂O and 0.3μl (1.7 U) of shrimp alkaline phosphatase (Sequenom) was added to digest redundant dNTPs with the following program: 37°C for 20 min, 85°C for 5 min, then maintained at 4°C. Five microliters of T Cleavage Transcription/RNase Cocktail, including 0.89 μl 5x T7 polymerase buffer, 0.22 μl T cleavage mix, 3.14 m Mdithiothreitol (DDT) 0.22μl, 22 U of T7 RNA and DNA Polymerase 0.40μ, 0.09 mg/ml RNase A, and 2 μl of product of the PCR/SAP reactions were mixed and incubated under the following conditions: 37°C for 3 h in vitro transcription and RNase A digestion. Then the mixture was further diluted with H₂O to 27 μl, purified with CLEAN resin (Sequenom) and robotically dispensed onto silicon chips preloaded with matrix (SpectroCHIP; Sequenom). The spectra and the methylation values of matrix-associated laser desorption/ ionization time-of-flight mass spectrometries (Sequenom) were collected and analyzed using Epityper software (version 1.0; Sequenom). Inapplicable readings and their corresponding sites were eliminated from analysis. The methylation level was expressed as the percentage of methylated cytosines over the total number of methylated and unmethylated cytosines.

**Fig. 1. Upstream sequence of mGluR1 core promoter region.**
(**Fig 1A**) The chart is the analysis result of Methyl Primer Express; red pillars stand for CpG sites (**Fig 1B**) Another analysis in the same region; blue dots represent valid CpG sites while red dots which are invalid cannot be detected in MassARRAY, valid dots were numbered.

Statistical analysis

The data were analyzed using the Graph Pad Prism (version 5). Results were presented as mean ± SEM. The Shapiro-Wilk test was used to test the distribution of the data and determine whether the data were parametric and non-parametric. Student t-tests were used to compare the genes and neurobehavioral study at different timelines among the first generation (F₀). One-way analysis of variance (ANOVA) with Newman-Keuls post hoc test was used to compare the genes and neurobehavioral characteristics across different groups in the second generation of animals (F₁). A p< 0.05 was considered statistically significant in all analysis.

Results

Sucrose preference test

Sucrose consumption was measured among groups of rats in both F₀ and F₁ generation. In the F₀ generation, rats with PFS showed significantly low sucrose preference when compared to SAL rats on PND 37 (p = 0.0004, t = 4.371) and PND 60 (p = 0.0006, t = 4.371) respectively (Fig 2A and 2B). In Fig 2C, there was a significant difference in PFS effect on the F₁ generation as offspring from both male and female PFS parents (PFSf + PFSm) on PND 37 showed significant low sucrose preference [F (3, 24) = 8.926, p = 0.0004] when compared to offspring from other groups of parents.

Immobility during the forced swim test

The immobile time spent during the forced swim test was assessed in SAL and PFS rats. Fig 3A showed that a prolonged febrile seizure effect was evident at PND 37 in the F₀ generation during the forced swim test as PFS animals spent significant higher immobile time (p<0.0001, t = 5.230) when compared to SAL animals. In Fig 3B,on PND 60, PFS rats in the F₀ generation spent significantly more immobile time when compared to SAL rats(p<0.0001, t = 5.886). There was no significant PFS effect across groups in the F₁ generation of rats [F (3, 24) = 0.2206, p = 0.8811] (Fig 3C).

Metabotropic glutamate receptor 5 (mGluR5) expression in hippocampal tissue

There was a PFS effect in the F₀ and F₁generations of rats as shown in Fig 4. On PND 14, PFS animals showed a significantly decreased hippocampal mGluR5 mRNA expression when compared to the saline groups (p = 0.0272, t = 2.905, Fig 4A). Also, on PND 37, rats with PFS showed significantly decreased expression of mGluR5 mRNA when compared to SAL (p = 0.0126, t = 3.515,Fig 4B). Similarly, this was noticed on PND 60 in the F₀ animals, PFS rats showed decreased expression of mGluR5 mRNA in the hippocampus when compared to SAL animals (p = 0.0186, t = 3.199,Fig 4C). As shown in Fig 4D, there was a significant PFS effect on hippocampal mGluR5 mRNA transferred from the parents to their offspring [F (3, 12) = 5.904, p = 0.0103]. Offspring from adult males with PFS and saline females (PFSm + SALf) showed significantly decreased mGluR5 mRNA expression when compared to offspring from adult male and female saline rats (SALm+SALf) (p<0.05) (Fig 4D). Offspring from adult male with PFS and saline female (PFSm+SALf) showed significantly low mGluR5 expression when compared to offspring from adult male saline and female PFS animals (PFSf + SALm) (p<0.05) (Fig 4D). Offspring from adult male and female rats with PFS (PFSm + PFSf) showed significantly low mGluR5 expression when compared to offspring from adult male and female saline rats (SALm+SALf) (p<0.05) (Fig 4D). Similar observations were made among offspring from adult male and female with PFS as they showed significantly low mGluR5 expression when compared to offspring from adult male saline and female PFS rats (PFSf+SALm) (p<0.05) (Fig 4D).

**Fig. 4. Effects of PFS on the expression of hippocampal mGluR5 mRNA.**
**Fig 4A** represent F₀ generation on PND 14*(SAL vs. PFS, p = 0.0272), **Fig 4B**: F₀ rats at PND 37 *(SAL vs. PFS, p = 0.0126), **Fig 4C**: F₀ at PND 60 *(SAL vs. PFS, p = 0.0188). **Fig 4D** showed the effect of PFS on mGluR5 mRNA in the F₁ generation on PND 37* (PFSm+SALfvs. SALm+SALf, PFSf+SALm p = 0.0103), #(PFSm+PFSf vs. SALm+SALf, PFSf+SALm, p = 0.0103). Results expressed as mean ± SEM (n = 4/group).

Metabotropic glutamate receptor 3 (mGluR3) expression in hippocampal tissue

There was a PFS effect in the F₀ generation of rats as shown in Fig 5. In Fig 5A among the F₀ generation on PND 14, PFS animals had significantly high hippocampal mGluR3 mRNA expression when compared to the SAL groups (p = 0.0343, t = 2.728, Fig 5A). Also, on PND 37 rats with PFS showed significant high expression of mGluR3 mRNA when compared to SAL (p = 0.0008, t = 6.203, Fig 5B). Similar effects were seen on PND 60 in the F₀generation PFS rats as there was increased expression of mGluR3 mRNA in the hippocampus when compared to SAL animals (p = 0.0207, t = 3.114,Fig 5C). As shown in Fig 5D, there was a PFS effect on hippocampal mGluR3 gene transferred from the parents to their offspring [F (3, 12) = 11, 24, p = 0.0008]. Offspring from adult male with PFS and saline female rats (PFSm+SALf) showed significantly higher mGluR3 mRNA expression when compared to offspring from adult male and female saline rats (SALm+SALf) (p<0.05) (Fig 5D). Offspring from adult male and female saline (SALm+SALf) showed significantly low mGluR3 mRNA expression when compared to offspring form adult male saline and female PFS (PFSf+SALm) (p<0.05) (Fig 5D). Offspring from adult male and female with PFS (PFSm+PFSf) showed significantly high mGluR3 mRNA expression when compared to offspring form adult male and female saline rats (SALm+SALf) (p<0.05) (Fig 5D).

**Fig. 5. Effects of PFS on the expression of hippocampal mGluR3mRNA.**
**Fig 5A** represents F₀ generation at PND 14*(SAL vs. PFS, p = 0.0343), **Fig 5B**: F₀ rats at PND 37 ***(SAL vs. PFS, p = 0.0008), **Fig 5C** F₀ at PND 60 *(SAL vs. PFS, p = 0.0207). **Fig 5D** showed the effect of PFS on mGluR3 gene in the F₁ generation at PND 37 * (PFSm+SALf vs. SALm+SALf), ** (PFSf+SALmvs. SALm+SALf), #(PFSm+PFSf vs. SALm+SALf). Results expressed as mean ± SEM (n = 4/group).

Metabotropic glutamate receptor 1 (mGluR1) methylation level in hippocampal tissue

On PND 14, as shown in Fig 6A, PFS led to a significant decreased level of methylation among the F₀ generation when compared to the saline groups at CpG_4 and CpG_5 (p<0.05).It was also observed that among them, there were significant higher levels of methylation at CpG-11 and CpG_15 when compared to the saline group (p< 0.05).No significant difference was found at other CpG units. The methylation level of the CpG regions was also evaluated on PND 37(Fig 6B). Fig 6B shows that in the F₀ generation there were significant increase in the methylation level of animals with prolonged febrile seizure when compared to saline at CpG_2, CpG _9, CpG_12, CpG_13 and CpG_14 (p< 0.05). There are no significant differences seen at other CpG units. Fig 6C shows that there were significant increases in the methylation level of animals in the F₀ generation on PND 60 with prolonged febrile seizure when compared to saline at CpG_9, with p< 0.05 on PND 60. There was no significant difference at other CpG units. Fig 6D shows the comparison between the different groups of offspring in the F₁ generation at the various CpG levels on PND 14. The results however showed that the methylation levels were significantly higher [F (3, 12) = 3,284, p = 0.0584]in offspring of PFSm+PFSf at CpG_8 when compared to offspring from both saline parents(SALm+SALf), but no significant difference in methylation level between the other CpG units in the F₁ generation of rats.

**Fig. 6.**
**Fig 6A:** Effects of PFS on the F₀ generation hippocampal methylation level of mGlR1 gene on PND 14 rats. Data are expressed as mean ±SEM (n = 4/group).*** (CpG _ 4 p = 0.0003), ** (CpG_6 p = 0.004), *(GpG_11 p = 0.02) and *(CpG_15 p = 0.04). **Fig 6B:** Effects of PFS on the hippocampal methylation level of mGlR1 gene in rats on PND 37in the F₀ generation. Data are expressed as mean±SEM (n = 4/group).*** (CpG _ 9 p = 0.0005), ** (CpG_14 p = 0.004), *(GpG_2 p = 0.04) and *(CpG_12; 13 p = 0.03). **Fig 6C:** Effects of PFS on the F₀ generation hippocampal methylation level of mGlR1 gene in rats on PND 60. Data are expressed as mean ±SEM (n = 4/group).* (CpG _ 9 p<0.0005). **Fig 6D:** Effects of PFS on the hippocampal methylation level of mGluR1 gene on F₁ generation rats. *(PFSm+PFSf vs. SALm+SALf, CpG _8 p = 0.0584). Data are expressed as mean ±SEM (n = 4/group).

Discussion

In this study, we examined the effect of neonatal PFS on SPT and FST in rats from adolescence to adulthood in the F₀ generation and at the adolescent age in the F₁ generation. We also looked at the effects of PFS on hippocampal mGluR1gene, mGlur5 mRNA and mGluR3mRNA at different timelines across both F₀ and F₁ generations of rats. We observed that animals with PFS developed anhendonia by consuming less sucrose solution in the SPT and showed despair by exhibiting less mobility time in the FST in the F₀ generation. This may be due to the early life exposure to brain insult caused by the prolonged febrile seizure [36]. This finding is in agreement with Crespo et al. [37], who observed that in their model of PFS, rats exhibited depressive-like behavior at adulthood. Furthermore, we also noticed that offspring from both parents with history of PFS consumed less sucrose solution; this may probably be due to genetic transfer and the level of care received from the dams. This is in agreement with previous study which reported that postpartum depressed mothers had offspring with depressive-like behavior [21]. Interestingly, we observed that there was no change in the offspring immobility time spent within the FST, this may be because loss of hope which is revealed by the FST is one of the late signs of depressive-like behavior [38]. This finding is in accordance with that of Wu et al. [21] who demonstrated that the combined effects of early life maternal separation and maternal stress were only able to increase the immobility time of the offspring in the FST. These depressive-like behaviors has been linked to altered expression of metabotropic receptors (mGluRs) and as such the mGluRs has been suggested as potential targets for drug treatment [39–42]

In this study we observed that mGluR5 mRNA hippocampal expression was down regulated in the PFS group within the F₀ generation across all time lines from childhood to adulthood (PND14-60) when compared to the SAL group. This may be due to the active receptor protein bound during excitation subsequently leading to low expression of the unbounded available receptor site in the hippocampus. Our finding is in agreement with the previous report of Iyo et al. [43], who observed low expression of mGluR5 in the hippocampus of their corticosteroid rat model of depression. Among the F₁ generation, it was observed in our study that offspring born to both parents with history of PFS and those born of PFSm+SALf parents showed decreased expression of hippocampal mGluR5 mRNA with associated decrease sucrose consumption. This may be due to the gender specific ways by which stress affect gene expression in the hippocampus [44]. Our findings ins in agreement with Wang et al, [40] who in their animal model of depression observed that hippocampal mGluR5 mRNA expression was low in the male offspring when compared to their female counterparts.

Metabotropic glutamate receptor 3(mGluR3) are autoreceptors that decrease excitatory glutamatergic transmission [10]. In this study, we observed that expression of hippocampal mGluR3 was upregulated among the PFS rats when compared to SAL animals at PND14, 37 and 60 in the F₀ generation. This is expected since mGluR3 reduces the excitatory effects of glutamate; hence the observed upregulation is due to the active excitatory mGluR5 receptor effect. This finding was corroborated by Stam et al. [45], who reported that in their stressed animals, mGluR3 expression was elevated. Among the F₁ generation, it was observed that offspring born of any parent with history of PFS had increased expression of hippocampal mGluR3 mRNA with associated decrease in sucrose consumption when compared to rats born from male and female saline parents. Similar findings was also noticed by Wang et al. [40] who in their prenatal stressed model observed hippocampus of offspring from stressed adult rats and found elevated mGluR3 expression.

Previous studies have demonstrated the involvement of DNA methylation, as a critical epigenetic change in the control of DNA transcription level [46,47]. While there are diverse consequences of DNA methylation on transcription regulation [48, 49], the alteration in DNA methylation level in a distinct region plays a crucial role in the final reaction of DNA methylation consequent to transcription regulation [50]. In this study, we observed the epigenetic effect of prolonged febrile seizure history at different time lines in rats since factors like age, gender, and the environment have been reported to influence DNA methylation [51].

From our methylation study it was observed that methylation of the mGluR1 gene was not of a particular pattern especially at PND14 this was consistent with findings by Bagot et al. [52], where they observed that Grm1 gene methylation difference which is not specific to a single CpG site. However, it was noticed that PFS resulted more in the hypermethylation of the CpG unit on postnatal day 14, 37 and 60 after the induction of seizure in F₀ rats. This pattern was also observed Bagot et al. [52], who observed hypermethylation of Grm1 gene in offspring of mothers that showed low care (low-LG) when compared to offspring of high-LG mothers. In the F₁ generation, animals born to both parents with history of PFS (PFSm+PFSf),showed increased methylation of mGluR1 when compared to their offspring from both saline parents (SALm+SALf). This may be due to the maternal care and gender specific genetic transfer from parents to offspring. This is in tandem with the findings of Zhang et al. [53], who discovered significant increased methylation level in GAD1 gene promoter in rat hippocampus in response to maternal care. Furthermore, Wu et al. [21], discovered that environmental factors led to epigenetic changes which bring about susceptibility to disease later on in life. Therefore, we observed that PFS history caused decreased sucrose consumption and immobility in rats in both the sucrose preference test and forced swim test respectively, a concomitant increase in methylation of mGluR1 and mRNA expression of mGluR3 with decreased expression of mGluR5 mRNA in the hippocampus. Our findings were consistent with Wang et al, 2015 [40] who in their chronic stress model observed that offspring of mother exposed to chronic mild stress exhibited depressive behavior, decreased mGluR5 hippocampal expression and increased mGluR2/3 expression. Our present observation was also in agreement with Bagot et al, 2012 [14] who showed that poor maternal care which could be as a result of early life time event in rats can lead to behavioral changes and hypermethylation of mGluR1 hippocampal expression in their offspring.

We concluded that PFS may profoundly affect mGluR1 methylation profile, mRNA expression of both mGluR5 and mGluR3 genes which may have led to depressive-like behavior in rats in both the F₀ and F₁ generations. Thus, the metabotropic receptor proteins could be a novel target in drug development towards the treatment of major depression.

Supporting information

S1 Table [pdf]
Showing the synthesis reports of mGluR5 (Grm5) and mGluR3 (Grm3) primers used to assay genes in the hippocampal tissue in this stud.

S2 Table [csv]
Showing the raw data generated from the Mass ARRAY analysis of mGluR1 gene methylation status used in this study.

Zdroje

1. Dube C., Chen K., Eghbal‐Ahmadi M., Brunson K., Soltesz I. and Baram T.Z., 2000. Prolonged febrile seizures in the immature rat model enhance hippocampal excitability long term. Annals of Neurology: Official Journal of the American Neurological Association and the Child Neurology Society, 47(3), pp.336–344.

2. Ackermann S. and Van Toorn R., 2011. Managing first-time seizures and epilepsy in children: A first seizure is a relatively common problem in paediatric general practice. Continuing Medical Education, 29(4).

3. Dubé C., Vezzani A., Behrens M., Bartfai T. and Baram T.Z., 2005. Interleukin‐1β contributes to the generation of experimental febrile seizures. Annals of Neurology: Official Journal of the American Neurological Association and the Child Neurology Society, 57(1), pp.152–155.

4. Cassim S., Qulu L. and Mabandla M.V., 2015. Prenatal stress and early life febrile convulsions compromise hippocampal genes MeCP2/REST function in mid-adolescent life of Sprague-Dawley rats. Neurobiology of learning and memory, 125, pp.195–201. doi: 10.1016/j.nlm.2015.09.002 26358644

5. Qulu L., Daniels W.M. and Mabandla M.V., 2012. Exposure to prenatal stress enhances the development of seizures in young rats. Metabolic brain disease, 27(3), pp.399–404. doi: 10.1007/s11011-012-9300-3 22527993

6. Marchi N., Granata T. and Janigro D., 2014.Inflammatory pathways of seizure disorders.Trends in neurosciences, 37(2), pp.55–65. doi: 10.1016/j.tins.2013.11.002 24355813

7. Barker-Haliski M. and White H.S., 2015.Glutamatergic mechanisms associated with seizures and epilepsy. Cold Spring Harbor perspectives in medicine, 5(8), p.a022863. doi: 10.1101/cshperspect.a022863 26101204

8. Lasoń W., Chlebicka M. and Rejdak K., 2013. Research advances in basic mechanisms of seizures and antiepileptic drug action. Pharmacological Reports, 65(4), pp.787–801. 24145073

9. Danbolt N. C. 2001. Glutamate uptake. Progress in Neurobiology, 65:1, 1–105 doi: 10.1016/s0301-0082(00)00067-8 11369436

10. Sanacora G., Zarate C.A., Krystal J.H. and Manji H.K., 2008.Targeting the glutamatergic system to develop novel, improved therapeutics for mood disorders.Nature reviews Drug discovery, 7(5), p.426. doi: 10.1038/nrd2462 18425072

11. Aronica E., Van Vliet E.A., Mayboroda O.A., Troost D., Da Silva F.H.L. and Gorter J.A., 2000. Upregulation of metabotropic glutamate receptor subtype mGluR3 and mGluR5 in reactive astrocytes in a rat model of mesial temporal lobe epilepsy. European Journal of Neuroscience, 12(7), pp.2333–2344. doi: 10.1046/j.1460-9568.2000.00131.x 10947812

12. Notenboom RG, Hampson DR, Jansen GH, van Rijen PC, van Veelen CW, van Nieuwenhuizen O, de Graan PN. 2006. Up-regulation of hippocampal metabotropic glu- tamate receptor 5 in temporal lobe epilepsy patients. Brain 129: 96–107. doi: 10.1093/brain/awh673 16311265

13. Pollock P.M., Cohen-Solal K., Sood R., Namkoong J., Martino J.J., Koganti A., Zhu H., Robbins C., Makalowska I., Shin S.S. and Marin Y., 2003. Melanoma mouse model implicates metabotropic glutamate signaling in melanocytic neoplasia. Nature genetics, 34(1), p.108. doi: 10.1038/ng1148 12704387

14. Bagot R.C., Tse Y.C., Nguyen H.B., Wong A.S., Meaney M.J. and Wong T.P., 2012. Maternal care influences hippocampal N-methyl-D-aspartate receptor function and dynamic regulation by corticosterone in adulthood. Biological psychiatry, 72(6), pp.491–498. doi: 10.1016/j.biopsych.2012.03.016 22521150

15. Lin T., Dang S., Su Q., Zhang H., Zhang J., Zhang L., Zhang X., Lu Y., Li H. and Zhu Z., 2018. The Impact and Mechanism of Methylated Metabotropic Glutamate Receptors 1 and 5 in the Hippocampus on Depression-Like Behavior in Prenatal Stress Offspring Rats. Journal of clinical medicine, 7(6), p.117.

16. Qureshi I.A. and Mehler M.F., 2010. Epigenetic mechanisms underlying human epileptic disorders and the process of epileptogenesis.Neurobiology of disease, 39(1), pp.53–60. doi: 10.1016/j.nbd.2010.02.005 20188170

17. Lubin F.D., 2012. Epileptogenesis: can the science of epigenetics give us answers?.Epilepsy currents, 12(3), pp.105–110. doi: 10.5698/1535-7511-12.3.105 22690136

18. Sterba S.K., Prinstein M.J. and Cox M.J., 2007. Trajectories of internalizing problems across childhood: Heterogeneity, external validity, and gender differences. Development and psychopathology, 19(2), pp.345–366. doi: 10.1017/S0954579407070174 17459174

19. Arai J.A., Li S., Hartley D.M. and Feig L.A., 2009. Transgenerational rescue of a genetic defect in long-term potentiation and memory formation by juvenile enrichment.Journal of Neuroscience, 29(5), pp.1496–1502. doi: 10.1523/JNEUROSCI.5057-08.2009 19193896

20. Kim H.K., Capaldi D.M., Pears K.C., Kerr D.C. and Owen L.D., 2009. Intergenerational transmission of internalising and externalisingbehaviours across three generations: Gender‐specific pathways. Criminal Behaviour and Mental Health, 19(2), pp.125–141. doi: 10.1002/cbm.708 19274624

21. Wu R., Zhang H., Xue W., Zou Z., Lu C., Xia B., Wang W. and Chen G., 2017.Transgenerational impairment of hippocampal Akt-mTOR signaling and behavioral deficits in the offspring of mice that experience postpartum depression-like illness.Progress in Neuro-Psychopharmacology and Biological Psychiatry, 73, pp.11–18. doi: 10.1016/j.pnpbp.2016.09.008 27693392

22. Dai Y.J., Wu D.C., Feng B., Chen B., Tang Y.S., Jin M.M., Zhao H.W., Dai H.B., Wang Y. and Chen Z., 2019. Prolonged febrile seizures induce inheritable memory deficits in rats through DNA methylation. CNS neuroscience & therapeutics.

23. Bohacek J., Gapp K., Saab B.J. and Mansuy I.M., 2013. Transgenerational epigenetic effects on brain functions.Biological psychiatry, 73(4), pp.313–320. doi: 10.1016/j.biopsych.2012.08.019 23062885

24. Fagiolini M., Jensen C.L. and Champagne F.A., 2009. Epigenetic influences on brain development and plasticity. Current opinion in neurobiology, 19(2), pp.207–212. doi: 10.1016/j.conb.2009.05.009 19545993

25. Kobow K. and Blümcke I., 2011. The methylation hypothesis: do epigenetic chromatin modifications play a role in epileptogenesis?.Epilepsia, 52, pp.15–19. doi: 10.1111/j.1528-1167.2011.03145.x 21732935

26. Eng J., 2003. Sample size estimation: how many individuals should be studied?.Radiology, 227(2), pp.309–313. doi: 10.1148/radiol.2272012051 12732691

27. Rakgantsho C. and Mabandla M.V., 2019. Acetylcholine receptor agonist effect on seizure activity and GABAergic mechanisms involved in prolonged febrile seizure development in an animal model. Brain research bulletin.

28. Heida J.G., Moshé S.L. and Pittman Q.J., 2009. The role of interleukin-1β in febrile seizures. Brain and Development, 31(5), pp.388–393. doi: 10.1016/j.braindev.2008.11.013 19217733

29. Thompson C.I., Brannon A.J. and Heck A.L., 2003. Emotional fever after habituation to the temperature-recording procedure.Physiology & behavior, 80(1), pp.103–108.

30. Racine R.J., 1972. Modification of seizure activity by electrical stimulation: II. Motor seizure.Electroencephalography and clinical neurophysiology, 32(3), pp.281–294. doi: 10.1016/0013-4694(72)90177-0 4110397

31. Slattery D.A., Markou A. and Cryan J.F., 2007. Evaluation of reward processes in an animal model of depression. Psychopharmacology, 190(4), pp.555–568. doi: 10.1007/s00213-006-0630-x 17177055

32. Santiago R. M., Tonin F. S., Barbiero J., Zaminelli T., Boschen S. L., Andreatini R., & Vital M. A. B. F. (2015). The nonsteroidal antiinflammatory drug piroxicam reverses the onset of depressive-like behavior in 6-OHDA animal model of Parkinson’s disease. Neuroscience, 300, 246–253. doi: 10.1016/j.neuroscience.2015.05.030 25999296

33. Yankelevitch-Yahav R., Franko M., Huly A. and Doron R., 2015.The forced swim test as a model of depressive-like behavior.JoVE (Journal of Visualized Experiments), (97), p.e52587.

34. Castagné V., Moser P., Roux S. and Porsolt R.D., 2011. Rodent models of depression: forced swim and tail suspension behavioral despair tests in rats and mice. Current Protocols in Neuroscience, 55(1), pp.8–10.

35. Schmittgen T.D. and Livak K.J., 2008. Analyzing real-time PCR data by the comparative C T method.Nature protocols, 3(6), p.1101. doi: 10.1038/nprot.2008.73 18546601

36. Cukor D. and McGinn L.K., 2006. History of child abuse and severity of adult depression: The mediating role of cognitive schema. Journal of Child Sexual Abuse, 15(3), pp.19–34. doi: 10.1300/J070v15n03_02 16893817

37. Crespo M., León-Navarro D.A. and Martín M., 2018. Early-life hyperthermic seizures upregulate adenosine A2A receptors in the cortex and promote depressive-like behavior in adult rats. Epilepsy & Behavior, 86, pp.173–178.

38. Porsolt R.D., Brossard G., Hautbois C. and Roux S., 2001. Rodent models of depression: forced swimming and tail suspension behavioral despair tests in rats and mice. Current protocols in neuroscience, 14(1), pp.8–10.

39. Lang U.E. and Borgwardt S., 2013. Molecular mechanisms of depression: perspectives on new treatment strategies. Cellular Physiology and Biochemistry, 31(6), pp.761–777. doi: 10.1159/000350094 23735822

40. Wang Y., Ma Y., Hu J., Cheng W., Jiang H., Zhang X., Li M., Ren J. and Li X., 2015. Prenatal chronic mild stress induces depression-like behavior and sex-specific changes in regional glutamate receptor expression patterns in adult rats. Neuroscience, 301, pp.363–374. doi: 10.1016/j.neuroscience.2015.06.008 26071959

41. Murrough J.W., Abdallah C.G. and Mathew S.J., 2017. Targeting glutamate signalling in depression: progress and prospects.Nature Reviews Drug Discovery, 16(7), p.472. doi: 10.1038/nrd.2017.16 28303025

42. Gulyaeva N.V., 2017. Interplay between brain BDNF and glutamatergic systems: a brief state of the evidence and association with the pathogenesis of depression. Biochemistry (Moscow), 82(3), pp.301–307.

43. Iyo A.H., Feyissa A.M., Chandran A., Austin M.C., Regunathan S. and Karolewicz B., 2010. Chronic corticosterone administration down-regulates metabotropic glutamate receptor 5 protein expression in the rat hippocampus. Neuroscience, 169(4), pp.1567–1574. doi: 10.1016/j.neuroscience.2010.06.023 20600666

44. Van den Hove D.L.A., Kenis G., Brass A., Opstelten R., Rutten B.P.F., Bruschettini M., Blanco C.E., Lesch K.P., Steinbusch H.W.M. and Prickaerts J., 2013.Vulnerability versus resilience to prenatal stress in male and female rats; implications from gene expression profiles in the hippocampus and frontal cortex.European Neuropsychopharmacology, 23(10), pp.1226–1246. doi: 10.1016/j.euroneuro.2012.09.011 23199416

45. Stam R., De Lange R.P., Graveland H., Verhave P.S. and Wiegant V.M., 2007.Involvement of group II metabotropic glutamate receptors in stress-induced behavioural sensitization.Psychopharmacology, 191(2), pp.365–375. doi: 10.1007/s00213-006-0659-x 17225168

46. Matouk C.C. and Marsden P.A., 2008. Epigenetic regulation of vascular endothelial gene expression.Circulation research, 102(8), pp.873–887. doi: 10.1161/CIRCRESAHA.107.171025 18436802

47. Yan M.S.C., Matouk C.C. and Marsden P.A., 2010. Epigenetics of the vascular endothelium.Journal of Applied Physiology, 109(3), pp.916–926. doi: 10.1152/japplphysiol.00131.2010 20413423

48. Fan A, Ma K, An X, Ding Y, An P, Song G, Tang L, Zhang S, Zhang P, Tan W, Tang B. Effects of Tet1 Knock-Down on Gene Expression and DNA Methylation in Porcine Induced Pluripotent Stem Cells. Reproduction. 2013 Sep 19:REP-13.

49. Frazier-Wood A.C., Aslibekyan S., Absher D.M., Hopkins P.N., Sha J., Tsai M.Y., Tiwari H.K., Waite L.L., Zhi D. and Arnett D.K., 2014. Methylation at CPT1A locus is associated with lipoprotein subfraction profiles. Journal of lipid research, 55(7), pp.1324–1330. doi: 10.1194/jlr.M048504 24711635

50. Hernandez H.G., Tse M.Y., Pang S.C., Arboleda H. and Forero D.A., 2013. Optimizing methodologies for PCR-based DNA methylation analysis.Biotechniques, 55(4), pp.181–197. doi: 10.2144/000114087 24107250

51. Ollikainen M., Smith K.R., Joo E.J.H., Ng H.K., Andronikos R., Novakovic B., Abdul Aziz N.K., Carlin J.B., Morley R., Saffery R. and Craig J.M., 2010. DNA methylation analysis of multiple tissues from newborn twins reveals both genetic and intrauterine components to variation in the human neonatal epigenome. Human molecular genetics, 19(21), pp.4176–4188. doi: 10.1093/hmg/ddq336 20699328

52. Bagot R.C., Zhang T.Y., Wen X., Nguyen T.T.T., Nguyen H.B., Diorio J., Wong T.P. and Meaney M.J., 2012. Variations in postnatal maternal care and the epigenetic regulation of metabotropic glutamate receptor 1 expression and hippocampal function in the rat. Proceedings of the National Academy of Sciences, 109(Supplement 2), pp.17200–17207.

53. Zhang T.Y., Hellstrom I.C., Bagot R.C., Wen X., Diorio J. and Meaney M.J., 2010.Maternal care and DNA methylation of a glutamic acid decarboxylase 1 promoter in rat hippocampus.Journal of Neuroscience, 30(39), pp.13130–13137. doi: 10.1523/JNEUROSCI.1039-10.2010 20881131