Abstract
Objectives:
A subset of children with Febrile Status Epilepticus (FSE) are at risk for development of temporal lobe epilepsy later in life. We sought a non-invasive predictive marker of those at risk, that can be identified soon after FSE, within a clinically realistic timeframe.
Methods:
Longitudinal T2-weighted magnetic resonance imaging (T2WI MRI) of rat pups at several time points after experimental (e)FSE was performed on a high-field scanner followed by long-term continuous electroencephalography. In parallel, T2WI MRI scans were performed on a 3.0T clinical scanner. Finally, chronic T2WI MRI signal changes were examined in rats that experienced eFSE and imaged month’s later in adulthood.
Results:
Epilepsy-predicting T2 changes, previously observed at two hours after eFSE, persisted for at least six hours, enabling translation to the clinic. Repeated scans, creating MRI trajectories of T2 relaxation times following eFSE provided improved prediction of epileptogenesis compared with a single MRI scan. Predictive signal changes centered on limbic structures, such as the basolateral and medial amygdala. T2WI MRI changes, originally described on high-field scanners can be also measured on clinical MRI scanners. Chronically elevated T2 relaxation times in hippocampus were observed months after eFSE in rats, as noted for post-FSE changes in children.
Significance:
Early T2WI MRI changes after eFSE provide a strong predictive measure of epileptogenesis following eFSE, both on high-field and clinical MRI scanners. Importantly, the extension of the acute signal changes to at least six hours after the FSE enables its inclusion in clinical studies. Chronic elevations of T2 relaxation times within the hippocampal formation and related structures are common to human and rodent FSE, suggesting that similar processes are involved across species.
Keywords: febrile status epilepticus, febrile seizures, magnetic resonance imaging, epileptogenesis, temporal lobe epilepsy
Introduction
Fever-associated seizures are common, occurring in 2–5% of children.1 Normally they are short and without long-term consequences, but seizures lasting >30 minutes are categorized as febrile status epilepticus (FSE) and are an important risk factor for developing temporal lobe epilepsy (TLE).2–4 Between the FSE and the first spontaneous epileptic seizure, there are years of epileptogenesis called the latent period.4–7 The latent period can last a decade or more, but currently clinicians are not able to predict which children will develop epilepsy following eFSE and which will remain healthy.8,9 A non-invasive technique to predict epileptogenesis would allow clinicians to appropriately counsel and monitor patients and ideally eventually provide a preventative intervention to those at risk before the first spontaneous seizure occurs.
To reach this goal, we previously reported that early magnetic resonance imaging (MRI) changes can predict epileptogenesis in an immature rat model of experimental FSE (eFSE).8 When rat pups underwent high-field 11.7T MRI scans 2h after eFSE, a decrease in T2 relaxation time within the basolateral amygdala (BLA) was predicative of which rats would develop epilepsy. This built on a number of studies in rodents and humans that have found MRI changes in the days, months and years after eFSE, reflecting long-term changes in the brain following a single episode of FSE.10–18 In contrast to the very early decrease (2h post-eFSE) in T2 signal reported in Choy et al., (2014), studies at later time points and on lower magnetic field scanners reported increased T2 relaxation times. For example, the FEBSTAT study, a prospective study of FSE in childhood found that both in an initial MRI scan (within a week of FSE) and follow up scans in the months and years following revealed hyperintensity (increased T2 signal) and volumetric changes in the hippocampus.10,11 The increased T2 relaxation time is consistent with our MRI finding in rats at a month following eFSE.12
This work was undertaken to advance the translatability of the early, predictive MRI signal to the clinic. Specifically, our goals were to: 1) increase the sensitivity and specificity of the prediction of epileptogenesis; 2) increase the time window for imaging evaluation, to allow clinically-relevant interval between the FSE and imaging, and 3) demonstrate that the changes that were observed at 11.7T high-field magnetic field scanner can be translated to the clinically available low-field scanners (i.e. 3.0T).
Methods:
All experimental procedures were approved by University of California, Irvine and Loma Linda University Institutional Animal Care and Use Committees and conformed to National Institutes of Health guidelines. Sprague-Dawley rats were maintained in quiet facilities under controlled temperatures and light-dark cycles. Cages were monitored every 12h and the date of birth was considered postnatal day 0 (P0). On P2, litters were culled to 10–11 pups, to encourage consistent development.
Induction of eFSE:
Experimental FSE was induced as previously described.8,19–21 Briefly, on P10–11, pups were placed in pairs in a 3.0L glass container lined with absorbent paper. Pups were subjected to a continuous stream of warm air until behavioral seizures began (3–5 minutes), characterized by sudden loss of motion (freezing), oral automatisms, and forelimb clonus. Hyperthermia (39.0–41.5°C) was maintained for 40min (Cohort 1) or 60min (other cohorts), an increased duration aiming to generate seizures lasting longer than 50 min, as suggested by human studies10,22. Core temperatures were measured at baseline, seizure onset, and every 2min during hyperthermia. Following eFSE, rats were cooled by placing on a cool metal surface or under running water and allowed to recover for 15 minutes before returning to their home cage.
The cohorts of rats studied were:
Cohort 1: Male only control (n=14) and eFSE rats (n=19) were serially scanned at 2, 18, and 48h at 11.7T, and underwent continuous EEG studies for up to 10 months. Whereas some data from this cohort were reported8, the analyses presented are previously unreported (Figure 1).
Cohort 2: Imaged in vivo at 11.7T 2 and 6h after eFSE (8 controls, 10 eFSE). Both males and females were used and randomly assigned to both groups (Figure 2).
Cohort 3: Imaged in vivo at 3.0T 4 and 48h after eFSE (7 controls, 9 eFSE). Both males and females were used and randomly assigned to both groups (Figure 3).
Cohort 4: Imaged in vivo at 11.7T 2, 48, and 96h (10 controls, 11 eFSE). Both males and females were used and randomly assigned to both groups (Figure 4)
Cohort 5: Imaged ex vivo at 9.4T during early adulthood (3.5m±0.8m; 9 controls, 4 eFSE). Both males and females were used and randomly assigned to both groups (Figures 5, 6)
Figure 1: The trajectory in T2 differences (48h – 2h) following eFSE is a better predictor of epileptogenesis than single early time point alone.
A) Representative pseudocolored 11.7 T2 maps of a control rat, an eFSE rat that did not develop epilepsy, and an eFSE rat that went on to develop epilepsy. B) Whole brain T2 values in individual animals of the three groups decreased significantly in control and eFSE-NoEpi, but not in eFSE-Epi animals, Control, eFSE-NoEpi, and eFSE-Epi (t=7.94, df=12, p<0.001; eFSE-NoEpi t=5.54, df=11, p<0.001; eFSE-Epi t=2.022, df=5, p=0.10) C) ROC curve of the predictive value of the delta T2 between 2 and 48h of the BLA, MeA, and DMThal (BLA: AUC 0.99 ± 0.020, p=0.001; MEA: AUC 1.00 ± 0, p=0.001; DMThal: AUC 0.83 ± 0.098, p<0.05) (inset: original ROC of 2h time point alone in inset, as published in Choy, et al., [2014]; BLA: AUC 0.91 ± 0.08, p= 0.005; MEA: AUC 0.82 ± 0.10, p<0.05; DMThal: AUC = 0.87 ±0.092; p=0.011). The BLA, MEA, and DMThal are able to differentiate between the eFSE-NoEpi and eFSE-Epi groups (D, E, F), while the entorhinal cortex does not (G). One-way ANOVA; Bonferroni Multiple Comparison Test (BLA Ctrl vs eFSE-NoEpi p=0.89; Ctrl vs. eFSE-NoEpi p<0.001; eFSE-NoEpi vs. eFSE-Epi p<0.001) (MeA: Ctrl vs eFSE-NoEpi p>0.99; Ctrl vs. eFSE-NoEpi p<0.001; eFSE-NoEpi vs. eFSE-Epi p=0.001) (DMThal: Ctrl vs eFSE-NoEpi p<0.001; Ctrl vs. eFSE-NoEpi p<0.001; eFSE-NoEpi vs. eFSE-Epi p<0.05) (Entorhinal Cortex: Ctrl vs eFSE-NoEpi p>0.99; Ctrl vs. eFSE-NoEpi p=0.07; eFSE-NoEpi vs. eFSE-Epi p=0.11) *p<0.05, **p<0.01, ***p<0.001
Figure 2: MRI T2 values do not change significantly between 2 and 6 hours.
A) Representative T2 maps of a rat at 2h and 6h after eFSE. The majority of rats remain in the same predictive group compared to controls at 2 and 6h after eFSE (C-E; Paired T-Test: BLA; p=0.45, t=0.78 df=17, correlation r=0.72, p<0.001; MeA: p=0.11, t=1.691, df=19, correlation r=0.81, p<0.001; DMThal: p=0.19, t=1.36, df=15, correlation r=0.89, p<0.001). F, G, H reveal the strong correlations between 2 and 6h time points for the BLA, MeA and DMThal
Figure 3: Immature rat MRI 4 and 48h after eFSE in a human 3.0T scanner is able to differentiate groups of eFSE rats at a whole brain level but not in individual regions.
A) Representative 3.0T T2 maps at 4 and 48h after eFSE. B) Whole brain 3.0T T2 values revealed consistent reductions in all control rats between 4 and 48h, but eFSE rats exhibit increased variability in their trajectories. (Paired T-test: Ctrl: p<0.001, t=17.23, df=5; eFSE p=0.003, t=6.177, df=8) C) Delta T2 (48h-4h) did not differentiate between the eFSE and control groups, but two clear groups within the eFSE emerge (T-Test, Ctrl vs. eFSE: p=0.12). By separating the eFSE animals into those with a similar trajectory as controls, and those with a smaller T2 decrease, the two groups were statistically different from each other and eFSE-Resilient was different from controls (dividing line at control + 2SD; 2206≥19.5) (one-way ANOVA, eFSE-Vulnerable vs. eFSE-Resilient: p<0.01; eFSE-Vulnerable vs. Control p<0.001; eFSE-Resilient vs. Control p>0.99) D-G: In vivo imaging of immature rats in a human 3T scanner did not reveal regional differences between groups (two-way ANOVA, no significant interaction between treatment group and time, significant effect of time. BLA: F(1,14) = 12.62, p<0.01; MEA: F(1,14) = 39.54, p<0.0001). **p<0.01, ***p<0.001
Figure 4: The 48h 11.7T MRI trajectories persist through 96h.
A) Representative 11.7T T2 maps for controls and 48h and 96h after eFSE. B, C) Individual 11.7T whole brain T2 values of eFSE and control rats revealed a decreased rate of change at 96h following eFSE. D, E, F) There was a significant difference in the trajectories from baseline for the basolateral amygdala (BLA) and medial amygdala (MeA) between control rats and rats that underwent eFSE, but no differences in the dorsal medial thalamus (DMThal). D’, E’) For the BLA and MeA, separating the eFSE rats into two groups those that followed a similar trajectory as control rats from 4–48h (eFSE-Resilient; n=6), and those with a small or no decrease (eFSE-Vulnerable; Δ48h-4h > control + 2SD; n=4), revealed two distinct trajectories (BLA: RM two-way ANOVA, Šidák Multiple Comparison test: significant interaction, F(4,34) = 5.46, p<0.01; 48h Ctrl vs. eFSE-Vulnerable p<0.01, eFSE-Vulnerable vs. eFSE-Resilient p <0.001; 96h Ctrl vs. eFSE-Vulnerable p<0.01, Ctrl vs. eFSE-Resilient p<0.05) (MeA: RM two-way ANOVA, Šidák Multiple Comparison test: significant interaction, F(4,34) = 6.60, p<0.001; 48h Ctrl vs. eFSE-Vulnerable p<0.01, eFSE-Vulnerable vs. eFSE-Resilient p <0.01; 96h Ctrl vs. eFSE-Vulnerable p<0.05, Ctrl vs. eFSE-Resilient p<0.001). *p<0.05, **p<0.01, ***p<0.001
Figure 5: MRI of adult rats (ex vivo) that underwent eFSE revealed increased T2 throughout the whole brain.
A) Group averaged images of control and eFSE rats demonstrate the globally elevated T2 values in the eFSE group B) Heat map of brain regions compartmentalized by region type (white matter, cortex, limbic associated regions, other), revealing increased in T2 relaxation times in FSE animals (n=4) relative to controls (n=9), particularly in limbic and cortical regions. A significant difference between the eFSE and control groups was observed when comparing all brain regions (two-way ANOVA, effect of eFSE, F(1,11) = 7.17, p<0.05). ***p<0.001
Figure 6: Limbic regions revealed increased T2 in adult animals that underwent eFSE in early life.
There is a significant T2 increase in the BLA, CA1, and CA3, with an increase trend across all limbic regions as shown in the heatmap (F). (Individual T-Test; BLA: p<0.05, t ratio = 2.63, df = 11; MeA: p=0.13, t ratio = 1.64, df = 11; CA1: p<0.05, t ratio = 2.80, df = 11, CA3; p<0.05, t ratio = 2.25, df = 11; Whole Brain: p = 0.057, t ratio =2.12, df = 11) (Group comparison: two-way ANOVA, effect of eFSE, F(1,11) = 5.87, p<0.05). BLA - Basolateral amygdala, MeA - Medial Amygdala, CA1/3 - cornu ammonis of the hippocampus 1/3; MeDLTh – Medial Dorsal Thalamic Nuclei. *p <0.05. ***p<0.001
In vivo MRI Procedure:
For all MRI studies, rats were lightly anesthetized using 1.5% isoflurane in 100% O2 to minimize motion, and body temperature was maintained at ~37°C with a heated water cushion.
11.7T in vivo MRI:
A single Bruker Avance 11.7T (Bruker Biospin, Billerica, MA) MR scanner (Research Imaging Center, LLU) was used for all 11.7T T2 studies for cohorts 1, 2, and 4. T2-weighted images were acquired using a 2D multi-echo-spin-echo sequence with a Bruker Biospin circular RF coil and the following parameters: Field of View (FOV): 2.3×2.3 cm, Slice Thickness: 0.75 mm, TR: 4697ms; TE: 10.21–100.1ms; inter-TE: 10.21ms, matrix size: 192×192, and number of averages (NA): 2.
3.0T in vivo MRI:
Cohort 3 rats underwent 3.0T T2 imaging on a single Phillips Achieva 3.0T MR scanner (Neuroscience Imaging Center, UCI). T2-weighted images were acquired using a clinical wrist coil a 2D multi-echo-spin-echo sequence with the following parameters: FOV: 2.3×2.3cm; matrix size: 152×153; slice thickness: 1.0 m; slice interval: 0.1mm; TR: 2000 ms; TE 17.20 – 51.6 ms; inter-TE = 17.20 ms; NA: 2.
In vivo MRI Analysis:
Absolute T2 relaxation times (ms) were calculated by log transform followed by linear least-squares fit on a pixel-by-pixel basis and T2 maps were generated using in-house software (MATLAB, Mathworks).
For all in vivo studies, regions were delineated manually without knowledge of treatment group using ImageJ software (versions 1.25–2.0.0). The regions of interest (ROIs) were drawn with a pseudo colored look up table (16 colors) to highlight borders between adjoining regions. All ROIs (shown in Supplemental Figure S1) were delineated by the same investigator (MMC), with a high level of intra-rater reliability (Intraclass Correlation Coefficient; BLA=0.872; MeA=0.832; DMThal=0.978, Whole Brain=0.992). The whole brain ROI was the entire brain on two consecutive slices, with the anterior slice aligning with the anterior BLA. MRI signal changes are often unilateral and always asymmetrical in children after FSE11,23, which is consistent with our previous findings8,14. Thus, we performed separate measurements and analysis of each side for all bilateral structures. Based on the results of the Choy et al, (2014) studies, we a priori chose to analyze data only the side with the lower T2. However, to ensure that there are no differences across time points we present bilateral data in Figure 2. For Cohort 1, the following ROIs were manually drawn: dorsal hippocampus, ventral hippocampus, entorhinal cortex, piriform cortex, cerebellum, medial amygdala, basolateral amygdala, medial thalamus, and corpus callosum. For the remaining in vivo cohorts, only BLA, medial amygdala, and dorsal medial thalamus were delineated.
Ex Vivo 9.4T MRI Acquisition and Analysis:
Cohort 5 rats (~3 months old) were anesthetized with a lethal dose of pentobarbital, perfused with 4% PFA, and brains were removed and post-fixed for 24h in PFA. Post-fixed brains submerged in Fluorinert FC-770 (Sigma-Aldrich) and T2-weighted images (T2WI) were acquired on a 9.4T Bruker Biospin MRI scanner (Billerica, MA, Paravision 5.1) (Experimental Imaging Center, Univ. of Calgary). The 10-echo T2WI images had the following parameters: 50 0.5mm slices, 1.92cm2 FOV, 256×256 matrix, TR: 6500ms, TE: 10ms, NA: 4 resulting in a total scan time of 15 min. Quantitative T2 maps were processed on JIM software (Xinapse Systems Ltd; West Bergholt, Essex; United Kingdom). The Waxholm MRI atlas24 was registered to each individual’s structural image and the Waxholm label atlas (separated by hemisphere in JIM) was transformed to this resulting image using Advanced Normalization Tools (ANTs).25 Native-space T2 values were extracted using the transformed Waxholm labels, and, consistent with acute measures, only the unilateral lower T2 values are presented. To create mean images (Figure 5A), the T2 maps were registered to a representative control rat’s space using FMRIB Software Library’s (FSL) linear registration tool, FLIRT, and averaged by group.
EEG electrode implantation & long-term video-EEG recordings and analysis:
At ~P40, cohort 1 rats had bipolar electrodes (PlasticsOne) implanted bilaterally in the hippocampus (AP:3.3; L:2.3; V: −2.8 mm to bregma), a cortical electrode was placed over the parietal cortex (AP:2; L: −2 mm), and a ground electrode over the cerebellum. EEG recordings were synchronized to video and conducted in freely moving rats beginning 5d after electrode implantation for up to 10 months, progressively increasing monitoring time to optimize seizure detection. Recordings increased from 112h in the first month (15.6% of the time) to 206h in the final month (3–5d segments; 28.6% of the month). Over 37,000h of video EEG were acquired including 595±57h per control rat and 1319±38h per eFSE rat.
EEGs scanned visually for seizures by two experienced investigators who were blinded to group identity and then reanalyzed using a seizure-detection software (LabChart version 7.3; ADInstruments) and concurrent video recordings were analyzed for behavioral manifestations of seizures. These included sudden cessation of activity, facial automatisms, head bobbing, prolonged immobility with staring, alternating or bilateral clonus, rearing and falling.26 Only events with both EEG and behavioral changes that lasted >20s were classified as seizures, and rats with at least one seizure were categorized as epileptic. The analysis and results defining the groups are presented in detail in Choy, et al. (2014).
Statistics:
Statistics and graphs were completed on GraphPad Prism (Version 7.0), except the Intraclass Correlation Coefficients, which were calculated using the VassarStat online calculator.28 Data are presented as box and whisker plots, with bars representing the minimum and maximum, with significance set at p<0.05. Outliers were excluded prior to analysis using the ROUT Test (Q = 1%).27 A list of outliers removed is presented in Supplementary Table S1. To determine whether MRI performed better than chance at predicting epilepsy after eFSE, regional MRI data from eFSE rats underwent receiver operating characteristic (ROC) analysis. Paired T-Test was used to determine significant change in the whole brain T2 relaxation time over 48h, as well as between 2–6 hours for all regions of interest, as well as measure group effects of the adult T2 images. MRI T2 relaxation differences in individual regions were compared using One-way repeated measures ANOVA with Bonferroni correction. To compare the group effects across time, Repeated Measure (RM-)Two-way ANOVA was used with Šidák Multiple Comparison test.
Results:
Dynamic changes over 48 hours in 11.7T T2 relaxation times are a robust predictor of epileptogenesis
In Cohort 1, 6/19 rats developed spontaneous seizures months after eFSE (31.6%)8, outlined in Table 1. As reported previously, T2 signal 2h after eFSE in the BLA and other limbic regions were good predictors of subsequent epileptogenesis. However, the prediction was incomplete, with the BLA predicting epilepsy at 83.3% sensitivity and 76.9% specificity (Figure 1C, inset).
Table 1:
Description of Seizures in Cohort 1
Rat ID | # of Seizures | Mean Seizure Duration (s) | Mean Racine Scale | Seizures/24h Recording |
---|---|---|---|---|
6 | 4 | 64.50 | 2.00 | 0.11 |
7 | 3 | 153.30 | 1.66 | 0.04 |
8 | 8 | 99.20 | 2.13 | 0.16 |
9 | 1 | 75.00 | 1.00 | 0.02 |
17 | 4 | 87.20 | 2.00 | 0.05 |
26 | 3 | 84.00 | 2.33 | 0.06 |
We extended these early studies by analyzing the longitudinal change in MRI signals 2 and 48h after eFSE by computing the difference in T2 relaxation times for each region of interest (Figure 1A). Consistent with normal myelination and maturation,29,30 we found a reduction in T2 relaxation time in the whole brain of both control rats as well as in the eFSE rats that did not develop epilepsy (eFSE-NoEpi) (Paired T-Test: controls; t=7.94, df=12, p<0.001; eFSE-NoEpi: t=5.54, df=11, p<0.001). Surprisingly, the T2 relaxation decrease was blunted across the whole brain of the eFSE rats that developed epilepsy later in life (eFSE-Epi) (Paired T-Test, t=2.02, df=5, p=0.10) (Figure 1B). These findings suggest that a disruption in the developmental reduction of T2 relaxation times throughout the brain may be predictive of epileptogenesis.
In our prior work, the BLA was the strongest predictor of epileptogenesis8. Here, we found that the longitudinal difference in T2 relaxation times from 2 to 48h following eFSE was a better predictor of epilepsy than the BLA at 2h alone. The trajectory of T2 values in the BLA of eFSE-Epi rats increased over the 48h following eFSE, whereas T2 relaxation times for both the controls and eFSE-NoEpi rats decreased (One-way ANOVA, eFSE-Epi vs. controls p<0.001; eFSE-Epi vs. eFSE-NoEpi, p<0.001) (Figure 1D).
Examining the potential value of other limbic structures, we discovered that the T2 difference in medial amygdala (MeA) was an even a better predictor of epilepsy than BLA (Figure 1E). In MeA, like the BLA, the difference in T2 relaxation times (48h - 2h) increased over time in eFSE-Epi compared to either controls or eFSE-NoEpi rats (One-way AVOVA, eFSE-Epi vs. controls, p<0.001; eFSE-Epi vs. eFSE-NoEpi, p=0.001). The dynamic changes of the dorsal medial thalamus (DMThal) also identifies differences following eFSE (Figure 1F, Ctrl vs eFSE-NoEpi p<0.001; Ctrl vs. eFSE-NoEpi p<0.001; eFSE-NoEpi vs. eFSE-Epi p=0.02).
The findings were selective, as these predictive changes were not observed in other limbic regions such as entorhinal cortex (Figure 1G; Ctrl vs eFSE-NoEpi p>0.99; Ctrl vs. eFSE-NoEpi p=0.07; eFSE-NoEpi vs. eFSE-Epi p=0.11). The same pattern of results was found when comparing the T2 values using a two-way ANOVA: the BLA, MeA and DMThal had a significant interaction between time and outcome. Notably, no interaction was found in the entorhinal cortex (Supplemental Figure S2). The predictive power of the T2 trajectory persisted when the individual T2 values were normalized to the whole brain T2 values, establishing that the results were not solely due to whole brain reduction of T2 values. No significant correlation was found when comparing the T2 change within the whole brain, BLA, MEA or MThal with the seizure burden of the individual animals.
The predictive efficacy of the T2 difference in the BLA, MeA and DMThal as markers of epileptogenesis was confirmed by the use of an independent unbiased measure, the receiver operating characteristic (ROC) curve analysis. ROC curve analysis is an unbiased measure of how a successful a test is at predicting the outcomes of a group of subjects. We used it to compare the predictions of epileptogenesis vs. the actual outcomes of the eFSE-NoEpi group and eFSE-Epi group. The curve for a test that has 100% sensitivity and 100% specificity will follow the top left corner and will have an area under the curve (AUC) of 1.00, whereas a test that performs as chance will be at 45° across the graph with an AUC of 0.50 (represented by the grey line). Using the ROC, the predictive value of the T2 differences was very robust (Figure 1C; AUC: BLA 0.99±0.02, p<0.001; MeA 1.00±0.00, p<0.001, DMThal 0.87±0.09, p=0.01).
Longevity of MRI changes of T2 relaxation times after experimental FSE
A challenge to the clinical translation of our original findings in Choy, et al., (2014) was the potential need to image children sustaining FSE within 2–4 hours after the initial insult. Therefore, we imaged a cohort of rats at both 2 and 6h after eFSE (Figure 2A, B). There was a strong correlation within each rat and no effect of time between 2 and 6h at a group level for T2 values in the BLA, MeA, DMThal, and whole brain (Figure 2C-H; Paired T-Test: BLA; p=0.45, t=0.78 df=17, correlation r=0.72, p<0.001; MeA: p=0.11, t=1.691, df=19, correlation r=0.81, p<0.001; DMThal: T-Test p=0.99, t=1.36, df=15, correlation r=0.89, p<0.001; Whole Brain [Supplemental Figure S3], T-Test: p=0.61, t=0.53, df=9, correlation r=0.75, p=0.01). These data demonstrate the stability of the epilepsy-predicting signal between 2h and 6h and indicate that imaging children at 6h after FSE could enable prediction of epilepsy.
Detection of T2 changes after eFSE is feasible in a clinically relevant low-field scanner
We imaged rat pups at 4 and 48h after eFSE on a 3.0T human scanner, the current standard for clinical MRIs (Figure 3A). The whole brain T2 trajectories recapitulated those observed from a high-field scanner (11.7T): controls had a significant and consistent reduction across the two days (Figure 3B; Paired T-Test: Ctrl: p<0.001, t=17.23, df=5). The eFSE group also had a significant decrease, but with increased variance (p<0.001, t=6.177, df=8). When analyzing whole brain T2 changes (48h - 4h), two distinct eFSE groups became apparent, which we termed eFSE-vulnerable and eFSE-resilient (Figure 3B). The eFSE-resilient rats fell within 2SD of the controls, demonstrating the expected reduction of T2 relaxation time similar to the controls (One-way ANOVA: p>0.99). The eFSE-vulnerable rats were characterized by a minimal reduction of T2 and was significantly different from both controls and eFSE-resilient groups (One-way ANOVA, eFSE-Vulnerable vs. eFSE-Resilient: p<0.01; eFSE-Vulnerable vs. Control p<0.001). There was no significant difference between either the T2 relaxation times (Figure 3D, F) or the change over time (Figure 3E, G, Supplemental Figure S4) in the BLA or MeA, likely a result of their small volume and the lower resolution inherent in imaging at 3T (two-way ANOVA, no significant interaction between treatment group and time; significant effect of time for both BLA and MeA. BLA: F(1,14) = 12.62, p<0.01; MeA: F(1,14) = 39.54, p<0.0001).
MRI T2 trajectory from days to months after eFSE
To extend our understanding of the acute and chronic effects of eFSE on T2 relaxation time trajectories, we imaged eFSE rats at additional acute time points and in adulthood. Delineating acute trajectories, (4, 48, and 96h after eFSE, Figure 4A) we found differing developmental patterns for control and eFSE rats across the whole brain and within specific limbic regions. Control rats had consistent reductions in individual T2 whole brain trajectories across 96h, whereas the eFSE group were more variable (Figure 4B). At a group level, there was a significant interaction for the whole brain between the effect of eFSE and imaging time-point (RM-Two-way ANOVA, Significant interaction, F(2,36) = 3.84, p<0.05), specifically at 96h (Šidák Multiple Comparison, p<0.01). Developmental patterns of T2 relaxation times within the BLA and MeA resembled those of whole brain. Notably, patterns within DMThal differed from those of amygdala nuclei: the age-dependent reduction of T2 in controls appeared delayed, commencing mainly after 48h, and the eFSE-induced inflection of the developmental trajectory, found in BLA, MeA and whole brain, was not observed. (Figures 4D-F, Supplemental Figure S5A-D), BLA: Significant interaction F(2,36) = 5.05, p<0.05, 96h Šidák p<0.001; MeA: Significant interaction F(2,36) = 5.40, p<0.01, 96h Šidák p<0.001; DMThal, no interaction, significant effect of time F(2,36) = 22.5, p<0.001).
The large variance within the eFSE group, especially in T2 for BLA and MeA at 48h following eFSE, prompted us to examine if the eFSE group was comprised of eFSE-resilient and vulnerable subsets. We separated the rats that were within 2SD of controls at 48h following eFSE (eFSE-resilient; BLA/MeA, n=6) from those outside of that range (eFSE-vulnerable, BLA/MeA, n=4). For both BLA and MeA, this stratification resulted two distinctive trajectories over the 48h between the eFSE-vulnerable and eFSE-resilient groups which converged at 96h (Figure 4D’, 4E’; Supplemental Figure S5B’-C’). For both the BLA and MeA, interaction was significant between time and treatment (RM-Two-way ANOVA, Šidák significant interaction, BLA: F(4,34) = 5.46, p<0.01; 48h Ctrl vs. eFSE-vulnerable p<0.01, eFSE-vulnerable vs. eFSE-resilient p <0.001; 96h Ctrl vs. eFSE-vulnerable p<0.01, Ctrl vs. eFSE-resilient p<0.05; MeA: F(4,34) = 6.60, p<0.001; 48h Ctrl vs. eFSE-vulnerable p<0.01, eFSE-vulnerable vs. eFSE-resilient p <0.01; 96h Ctrl vs. eFSE-vulnerable p<0.05, Ctrl vs. eFSE-resilient p<0.001). Whereas the basis for these apparent subgroups is unclear, it may suggest that the eFSE-vulnerable group demonstrated an immediate (0–48h) delay in the developmental trajectory of T2 relaxation, whereas a disruption of this T2 developmental trajectory emerged later in the “resilient” group.
The long-term temporal profile of eFSE-induced changes in T2 was examined in adults on a 9.4T scanner (age = 3.5m±0.8m). We found a brain-wide T2 increase in adult eFSE rats compared to controls, as apparent from the group averaged T2 maps (Figure 5A). When whole brain regions were parsed, there was a significant effect of eFSE (Figure 5B, two-way ANOVA, effect of eFSE, F(1,11) = 7.17, p<0.05). T2 increases were observed in eFSE rats in the limbic regions that predicted epileptogenesis (Figure 6A-E), particularly the BLA (T-Test, p<0.05, t=2.63, df=11). Interestingly, similar findings were noted for the CA1 and CA3 of the hippocampus, regions which had not shown acute changes (CA1: p<0.05, t=2.80, df=11, CA3; p<0.05, t=2.25, df=11). In contrast, T2 relaxation times in the MeA, medial dorsal thalamic nuclei, and dentate gyrus were not significantly different between eFSE and control groups. A group comparison of the aggregate BLA, MeA, hippocampal regions, and MeDLThal between eFSE and control rats did provide a robust measure of the consequences of eFSE, revealing a significant effect of eFSE on T2 values (Figure 6F, two-way ANOVA, effect of eFSE, F(1,11) = 5.87, p<0.05).
Discussion:
Here we provide novel information regarding the trajectories of T2 imaging-related brain changes that occur both short- and long-term following eFSE. The principal discoveries are: A) acute longitudinal MR imaging trajectories following eFSE enhance prediction of epileptogenesis over a single MRI scan; B) MRI changes following eFSE persist for at least 6h; C) T2 changes following eFSE can also be measured on clinical MRI scanners; and D) Increased T2, emblematic of post-FSE changes in children, can be observed in eFSE-experiencing rats months after the insult.
Progressive changes in T2 following eFSE offer improved prediction of epileptogenesis
Our work previously discovered high-field T2 signal changes as predictive marker of epileptogenesis immediately (2h) following a single episode of eFSE.8 This was a critical first step in the ability to predict epilepsy before the first spontaneous seizure occurred. While the prediction was robust, the positive predictive value (PPV) was imperfect: one-third of rats that were predicted to become epileptic did not have a spontaneous seizure. By analyzing the change in T2 rather than at a single time point, our new data demonstrates that the repeated longitudinal T2 changes (48h – 2h) increased prediction of epileptogenesis, increasing the PPV to 100% and 85.7% in the MeA and BLA, respectively. Additionally, this work revealed dynamic MRI differences at multiple time points (48 and 96h), reflecting the fluid process of epileptogenesis. We believe that the improved predictive value is based on the ability of serial MR imaging to differentiate how the temporal profile of T2 relaxation values of each individual rat either reflects that normal developmental decrease in T2 relaxation or differs from it, and that it is a lack of developmental T2 reduction that is the marker of early processes of epileptogenesis. Further studies are required to uncover the biological mechanisms that underlie the relationship between “normal” T2 decreases and the flattened trajectory of rats in the early stages of epileptogenesis.
Reconciliation of early T2 reduction and chronic T2 increases following eFSE
There has been a significant body of work that has found T2 increases in the brain (particularly hippocampus) following FSE in both humans and rodents.10–16 Because of this, the results of Choy, et al., (2014) which revealed a decrease in T2 after FSE 2h after eFSE and no changes at 48h were unexpected, though there were two important differences between it and previous work. First, the imaging was completed much earlier after eFSE. Secondly, the rats were imaged on an 11.7T scanner, which has a higher magnetic field and is affected by paramagnetic changes (T2*) which can be measured as decreased T2 values.8
By analyzing the longitudinal change of T2 values between 2 to 96h after eFSE, it became clear that a T2 increase in the days after eFSE was potentially masked by typical development. This developmental reduction in T2 relaxation times was measured across the entire brain and occurs as a result of normal myelination and maturation.29,30 The trajectory is more rapid in infant rats than it is in humans, thus the developmental reduction in infants would not play as important a role clinically. Whereas the underlying cellular mechanisms that lead to FSE-induced increased T2 are not fully understood, they have been postulated to involved edema, based on volumetric human studies,10,15 or gliosis, as our previous work has reported an increased number and activation of astrocytes following eFSE.12,20
Importantly, the trajectory between 4, 48 and 96h after eFSE reconciles the early and late T2 changes. The current findings demonstrate divergent developmental paths of eFSE rats compared with their control littermates. The same regions that have an early T2 decrease predictive of epileptogenesis demonstrate a later T2 increase in adulthood. Interestingly, early imaging of the hippocampus was not predictive of epileptogenesis8, but late increases in T2 are robust in both human and rodent images.11–13,17,18
Clinical translation of predictive T2 signal changes
The present study addresses two main hurdles for translating the predictive MRI signal to clinic: timeline and scanner strength. The original time-point, 2h after FSE, would be difficult to achieve in an emergency department setting, due to access to MRI and time required for consent. Our new findings that the early T2 signal decrease remains stable for up to 6h makes future clinical imaging of children feasible.
Additionally, our initial study was performed on research scanners (4.7T-11.7T), which have higher magnetic fields and smaller coils than the human scanners (1.5T-3.0T) that are typically available. To demonstrate that similar results would be likely found in children, we imaged eFSE rat pups on a 3.0T human scanner. Using a clinical human wrist coil, we observed measurable differences in the brains of rat pups following eFSE, similar to the findings in the high-field scanner. These differences were measurable at the whole brain level, not in specific limbic regions, but this is not too surprising due to resolution differences: the voxel volume at 11.7T was 0.011 mm3, but 0.023 mm3 in the 3.0T. This inherent resolution difference should not play as critical a role when imaging children, as the whole brain of P12 rat (which had measurable differences at 3.0T) is ~1 cm3 in volume, roughly the same volume of an amygdala of a one-year-old infant.31,32 Thus, it is likely the whole brain differences in the rat pup will translate to even stronger measurements when analyzing individual brain regions in children, mirroring the results found in rats on the animal scanners.
In addition to the noted resolution differences, it is important to remember that T2 relaxation properties are modulated by field strengths, age at imaging and acquisition parameters. Human and rodent T2 values are known to decrease with increasing field strength,33,34 (see Choy et al., [2014] for discussion). In the rodent brain, it has been suggested that signal to noise ratio’s may actually decrease as field strength increases,34 thus imaging at the clinically relevant 3T as we did in this study may in fact enhance our findings at this field strength. In addition, the influence of T2 values by field strength is similar in white and gray matter even if acquisition parameters differ.34 Of course, the use of high signal to noise RF coils and purpose-built coils would enhance the acquisition of T2 signals in neonates and adults.
Overall, this study clearly demonstrates the value of quantitative T2 MR imaging as a marker for the alterations in the brain following eFSE which occur at the onset of epileptogenesis. Whereas long term EEG recording was only completed in one group of rats following eFSE, the strength of the longitudinal MRI results in the EEG cohort allows us to predict the epileptogenic group utilizing short term MRI changes in the other cohorts. Moving forward, it will be important to both confirm these findings with a long-term EEG cohort and study the mechanisms underlying the MRI signal changes. This continued research will inform us how epilepsy develops, help us understand how we can prevent it, and move us towards translating these findings to the clinic.
Supplementary Material
Representative regions of interest for the basolateral amygdala (BLA), medial amygdala (MeA), and Dorsal Medial region of the Thalamus (DM Thalamus) shown on a T2 relaxation map.
There is a significant interaction between the effect of time and outcome on T2 relaxation time in the 48h following eFSE A) Basolateral Amygdala: RM two-way ANOVA, Šidák Multiple Comparison test: significant interaction, F(2,29) = 13.92, p<0.001; 2h: Ctrl vs. eFSE-NoEpi p=0.001, Ctrl vs. eFSE-Epi p<0.001, eFSE-NoEpi vs. eFSE-Epi: p=0.04; 48h: Ctrl vs. eFSE-NoEpi p=0.04. B) Medial Amygdala: significant interaction, F(2,29) = 10.36, p<0.001; 2h: Ctrl vs. eFSE-NoEpi p<0.001, Ctrl vs. eFSE-Epi p<0.001; 48h: Ctrl vs. eFSE-NoEpi p<0.01, eFSE-NoEpi vs. eFSE-Epi: p<0.05. C) Dorsal Medial Thalamus: significant interaction, F(2,29) = 14.72, p<0.001; 2h: Ctrl vs. eFSE-NoEpi p<0.001, Ctrl vs. eFSE-Epi p<0.001. D) Entorhinal Cortex: no significant interaction, F(2,29) = 3.13, p>0.05; effect of outcome, p>0.05; effect of time, F(2,29) = 45.81, p<0.001
MRI T2 value across the whole brain does not change between 2h and 6h after eFSE. A) Paired T-Test reveals no difference between groups 2h and 6h after eFSE (p=0.61, t=0.53, df=9). B) There is a strong correlation between the 2h and 6h whole brain data (r=0.75, p=0.01).
There is a significant effect of time but no effect of treatment group on the T2 relaxation time in the BLA (A) or the MeA (B) when measured with 3T MRI (two-way ANOVA, no significant interaction between treatment group and time, significant effect of time for both BLA and MeA. BLA: F(1,14) = 12.62, p<0.01; MeA: F(1,14) = 39.54, p<0.0001).
There is a significant interaction between time and treatment on the T2 relaxation times at 2, 48, and 96h following eFSE. This is apparent both comparing the eFSE and control groups (A-D; BLA: Significant interaction F(2,36) = 5.05, p<0.05, 96h Šidák p<0.001; MeA: Significant interaction F(2,36) = 5.40, p<0.01, 96h Šidák p<0.001; DMThal, no interaction, significant effect of time F(2,36) = 22.5, p<0.001), and when the eFSE rats are broken up into Resilient and Vulnerable groups based upon T2 changes in the first 48h (B’-C’; RM-Two-way ANOVA, Šidák significant interaction, BLA: F(4,34) = 5.46, p<0.01; 48h Ctrl vs. eFSE-vulnerable p<0.01, eFSE-vulnerable vs. eFSE-resilient p <0.001; 96h Ctrl vs. eFSE-vulnerable p<0.01, Ctrl vs. eFSE-resilient p<0.05; MeA: F(4,34) = 6.60, p<0.001; 48h Ctrl vs. eFSE-vulnerable p<0.01, eFSE-vulnerable vs. eFSE-resilient p <0.01; 96h Ctrl vs. eFSE-vulnerable p<0.05, Ctrl vs. eFSE-resilient p<0.001)
Key Points:
Longitudinal magnetic resonance imaging (MRI) trajectories post-febrile status epilepticus (FSE) enhance prediction of epileptogenesis
Predictive MRI changes following FSE persist for at least 6h and can also be measured on clinical scanners
Chronically increased T2, emblematic of post-FSE changes in children, can be observed in FSE-experiencing rats
Acknowledgements:
This work was supported by NIH grants R37/RO1 NS35439; T32 NS45540; and T32 GM008620.
Footnotes
Statement of Ethics: We confirm that we have read the Journal’s position on issues involved in ethical publication and affirm that this report is consistent with those guidelines. The ARRIVE guidelines and Basel Declarations were adhered to, and efforts were taken in experimental designs to allow for replacement, reduction and refinement of the animals that were used and minimize the pain and suffering of those that were included.
Conflict of Interest: None of the authors has any conflict of interest to disclose.
References
- 1.Baram TZ, Shinnar S. Febrile Seizures. San Diego, CA: Academic press; 2002. [Google Scholar]
- 2.Annegers JF, Hauser WA, Shirts SB, et al. Factors prognostic of unprovoked seizures after febrile convulsions. N Engl J Med. 1987;316:493–8. [DOI] [PubMed] [Google Scholar]
- 3.Cendes F, Andermann F, Dubeau F, et al. Early childhood prolonged febrile convulsions, atrophy and sclerosis of mesial structures, and temporal lobe epilepsy: an MRI volumetric study. Neurology. 1993;43:1083–7. [DOI] [PubMed] [Google Scholar]
- 4.Dubé CM, Brewster AL, Richichi C, et al. Fever, febrile seizures and epilepsy. Trends Neurosci. 2007;30:490–6. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 5.Berg AT, Shinnar S. Complex febrile seizures. Epilepsia. 1996; 37:126–33. [DOI] [PubMed] [Google Scholar]
- 6.French JA, Williamson PD, Thadani VM, et al. Characteristics of medial temporal lobe epilepsy: I. Results of history and physical examination. Ann Neurol. 1993;34:774–80. [DOI] [PubMed] [Google Scholar]
- 7.Mathern GW, Babb TL, Vickrey BG, et al. The clinical-pathogenic mechanisms of hippocampal neuron loss and surgical outcomes in temporal lobe epilepsy. Brain. 1995;118:105–18. [DOI] [PubMed] [Google Scholar]
- 8.Choy M, Dubé CM, Patterson K, et al. A novel, noninvasive, predictive epilepsy biomarker with clinical potential. J Neurosci. 2014;34:8672–84. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 9.Patterson KP, Baram TZ, Shinnar S. Origins of Temporal Lobe Epilepsy : Febrile Seizures and Febrile Status Epilepticus. Neurotherapeutics. 2014; 11:242–50. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 10.Lewis DV, Shinnar S, Hesdorffer DC, et al. Hippocampal sclerosis after febrile status epilepticus: the FEBSTAT study. Ann Neurol. 2014;75:178–85. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 11.Shinnar S, Bello JA, Chan S, et al. MRI abnormalities following febrile status epilepticus in children: the FEBSTAT study. Neurology. 2012;79:871–7. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 12.Dubé CM, Ravizza T, Hamamura M, et al. Epileptogenesis provoked by prolonged experimental febrile seizures: mechanisms and biomarkers. J Neurosci. 2010;30:7484–94. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 13.Jansen JFA, Lemmens EMP, Strijkers GJ, et al. Short- and long-term limbic abnormalities after experimental febrile seizures. Neurobiol Dis. 2008;32:293–301. [DOI] [PubMed] [Google Scholar]
- 14.Dubé C, Yu H, Nalcioglu O, et al. Serial MRI after experimental febrile seizures: altered T2 signal without neuronal death. Ann Neurol. 2004;56:709–14. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 15.Scott RC, Gadian DG, King MD, et al. Magnetic resonance imaging findings within 5 days of status epilepticus in childhood. Brain. 2002;125:1951–9. [DOI] [PubMed] [Google Scholar]
- 16.Dubé CM, Zhou J, Hamamura M, et al. Cognitive dysfunction after experimental febrile seizures. Exp Neurol. 2009;215:167–77. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 17.Provenzale JM, Barboriak DP, VanLandingham K, et al. Hippocampal MRI Signal Hyperintensity After Febrile Status Epilepticus Is Predictive of Subsequent Mesial Temporal Sclerosis. Am J Roentgenol. 2008;190:976–83. [DOI] [PubMed] [Google Scholar]
- 18.VanLandingham KE, Heinz ER, Cavazos JE, et al. Magnetic resonance imaging evidence of hippocampal injury after prolonged focal febrile convulsions. Ann Neurol. 1998;43:413–26. [DOI] [PubMed] [Google Scholar]
- 19.Chen K, Baram TZ, Soltesz I. Febrile seizures in the developing brain result in persistent modification of neuronal excitability in limbic circuits. Nat Med. 1999;5:888–94. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 20.Patterson KP, Brennan GP, Curran M, et al. Rapid, Coordinate Inflammatory Responses after Experimental Febrile Status Epilepticus: Implications for Epileptogenesis. eNeuro. 2015;2:1–12. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 21.Koyama R, Tao K, Sasaki T, et al. GABAergic excitation after febrile seizures induces ectopic granule cells and adult epilepsy. Nat Med. 2012;18:1271–8. [DOI] [PubMed] [Google Scholar]
- 22.Hesdorffer DC, Benn EKT, Bagiella E, et al. Distribution of febrile seizure duration and associations with development. Ann Neurol. 2011;70:93–100. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 23.Scott RC, King MD, Gadian DG, et al. Hippocampal abnormalities after prolonged febrile convulsion: A longitudinal MRI study. Brain. 2003;126:2551–7. [DOI] [PubMed] [Google Scholar]
- 24.Papp EA, Leergaard TB, Calabrese E, et al. Waxholm Space atlas of the Sprague Dawley rat brain. Neuroimage. 2014;97:374–86. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 25.Avants BB, Tustison NJ, Song G, et al. A reproducible evaluation of ANTs similarity metric performance in brain image registration. Neuroimage. 2011;54:2033–44. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Racine RJ. Modification of seizure activity by electrical stimulation. II. Motor seizure. Electroencephalogr Clin Neurophysiol.1972; 32:281–94. [DOI] [PubMed] [Google Scholar]
- 27.Motulsky HJ, Brown RE, Barnett V, et al. Detecting outliers when fitting data with nonlinear regression – a new method based on robust nonlinear regression and the false discovery rate. BMC Bioinformatics.2006; 7:123. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28.Correlation for Unordered Pairs: Intraclass Correlation [Internet]. Vassar Stats. 2018. [cited 2018]. Available from: http://vassarstats.net/icc.html
- 29.Matsumae M, Kurita D, Atsumi H, et al. Sequential changes in MR water proton relaxation time detect the process of rat brain myelination during maturation. Mech Ageing Dev. 2001;122:1281–91. [DOI] [PubMed] [Google Scholar]
- 30.Ferrie JC, Barantin L, Saliba E, et al. MR Assessment of the brain maturation during the perinatal period: Quantitative T2 MR study in premature newborns. Magn Reson Imaging. 1999;17:1275–88. [DOI] [PubMed] [Google Scholar]
- 31.Uematsu A, Matsui M, Tanaka C, et al. Developmental Trajectories of Amygdala and Hippocampus from Infancy to Early Adulthood in Healthy Individuals. PLoS One. 2012;7:e46970. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 32.Calabrese E, Badea A, Watson C, et al. A quantitative magnetic resonance histology atlas of postnatal rat brain development with regional estimates of growth and variability. Neuroimage. 2013;71:196–206. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 33.Schenck JF. Imaging of brain iron by magnetic resonance: T2 relaxation at different field strengths. J Neurol Sci. 1995;134 Suppl:10–8. [DOI] [PubMed] [Google Scholar]
- 34.de Graaf RA, Brown PB, McIntyre S, et al. High magnetic field water and metabolite proton T1 and T2 relaxation in rat brain in vivo. Magn Reson Med. 2006;56:386–94. [DOI] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
Representative regions of interest for the basolateral amygdala (BLA), medial amygdala (MeA), and Dorsal Medial region of the Thalamus (DM Thalamus) shown on a T2 relaxation map.
There is a significant interaction between the effect of time and outcome on T2 relaxation time in the 48h following eFSE A) Basolateral Amygdala: RM two-way ANOVA, Šidák Multiple Comparison test: significant interaction, F(2,29) = 13.92, p<0.001; 2h: Ctrl vs. eFSE-NoEpi p=0.001, Ctrl vs. eFSE-Epi p<0.001, eFSE-NoEpi vs. eFSE-Epi: p=0.04; 48h: Ctrl vs. eFSE-NoEpi p=0.04. B) Medial Amygdala: significant interaction, F(2,29) = 10.36, p<0.001; 2h: Ctrl vs. eFSE-NoEpi p<0.001, Ctrl vs. eFSE-Epi p<0.001; 48h: Ctrl vs. eFSE-NoEpi p<0.01, eFSE-NoEpi vs. eFSE-Epi: p<0.05. C) Dorsal Medial Thalamus: significant interaction, F(2,29) = 14.72, p<0.001; 2h: Ctrl vs. eFSE-NoEpi p<0.001, Ctrl vs. eFSE-Epi p<0.001. D) Entorhinal Cortex: no significant interaction, F(2,29) = 3.13, p>0.05; effect of outcome, p>0.05; effect of time, F(2,29) = 45.81, p<0.001
MRI T2 value across the whole brain does not change between 2h and 6h after eFSE. A) Paired T-Test reveals no difference between groups 2h and 6h after eFSE (p=0.61, t=0.53, df=9). B) There is a strong correlation between the 2h and 6h whole brain data (r=0.75, p=0.01).
There is a significant effect of time but no effect of treatment group on the T2 relaxation time in the BLA (A) or the MeA (B) when measured with 3T MRI (two-way ANOVA, no significant interaction between treatment group and time, significant effect of time for both BLA and MeA. BLA: F(1,14) = 12.62, p<0.01; MeA: F(1,14) = 39.54, p<0.0001).
There is a significant interaction between time and treatment on the T2 relaxation times at 2, 48, and 96h following eFSE. This is apparent both comparing the eFSE and control groups (A-D; BLA: Significant interaction F(2,36) = 5.05, p<0.05, 96h Šidák p<0.001; MeA: Significant interaction F(2,36) = 5.40, p<0.01, 96h Šidák p<0.001; DMThal, no interaction, significant effect of time F(2,36) = 22.5, p<0.001), and when the eFSE rats are broken up into Resilient and Vulnerable groups based upon T2 changes in the first 48h (B’-C’; RM-Two-way ANOVA, Šidák significant interaction, BLA: F(4,34) = 5.46, p<0.01; 48h Ctrl vs. eFSE-vulnerable p<0.01, eFSE-vulnerable vs. eFSE-resilient p <0.001; 96h Ctrl vs. eFSE-vulnerable p<0.01, Ctrl vs. eFSE-resilient p<0.05; MeA: F(4,34) = 6.60, p<0.001; 48h Ctrl vs. eFSE-vulnerable p<0.01, eFSE-vulnerable vs. eFSE-resilient p <0.01; 96h Ctrl vs. eFSE-vulnerable p<0.05, Ctrl vs. eFSE-resilient p<0.001)