ORIGINAL RESEARCH article

Front. Neurol., 08 May 2019

Sec. Epilepsy

Volume 10 - 2019 | https://doi.org/10.3389/fneur.2019.00488

Relationship Between Seizure Frequency and Functional Abnormalities in Limbic Network of Medial Temporal Lobe Epilepsy

  • 1. Department of Neurology, Mayo Clinic, Rochester, MN, United States

  • 2. Department of Neurologic Surgery, Mayo Clinic, Rochester, MN, United States

  • 3. Department of Radiology, Mayo Clinic, Rochester, MN, United States

  • 4. Department of Psychiatry and Psychology, Mayo Clinic, Rochester, MN, United States

Abstract

Background: We compared resting-state functional connectivity (RSFC) among limbic and temporal lobe regions between patients with medial temporal lobe epilepsy (mTLE) and healthy control subjects to identify imaging evidence of functional networks related to seizure frequency, age of seizure onset, and duration of epilepsy.

Methods: Twelve patients with drug-resistant, unilateral medial temporal lobe epilepsy and 12 healthy control subjects matched for age, sex, and handedness participated in the imaging experiments. We used network-based statistics to compare functional connectivity graphs in patients with mTLE and healthy controls to investigate the relationship between functional connectivity abnormalities and seizure frequency.

Results: Among mTLE patients, we found functional network abnormalities throughout the limbic system, but primarily in the hemisphere ipsilateral to the seizure focus. The RSFCs between ipsilateral hypothalamus and ventral anterior cingulate cortex and between ipsilateral subiculum and contralateral posterior cingulate cortex were highly correlated with seizure frequency.

Discussion: These findings suggest that in mTLE, changes in limbic networks ipsilateral to the epileptic focus are common. The pathological changes in connectivity between cingulate cortex, hypothalamus and subiculum ipsilateral to the seizure focus were correlated with increased seizure frequency.

Introduction

The limbic system of the brain is a complex network of structurally and functionally linked anatomic regions (). Among its components are areas in anteromedial temporal lobe including the amygdala, hippocampal complex, and entorhinal cortex, as well as the cingulate cortex. These areas are interconnected with each other via the thalamus (particularly the anterior thalamic and medial dorsal nuclear groups as well as the midline thalamic nuclear group), and project to the hypothalamus and midbrain, and anteromedial temporal lobe. The limbic system is an important component of the systems subserving emotion, behavior, and memory. In patients with medial temporal lobe epilepsy (mTLE), perturbations of the limbic system often result in in debilitating comorbidities and functional impairments in addition to the direct consequences of seizures.

The role of the limbic system in epilepsy is well-described, with mTLE frequently involving limbic structures and the pathological hallmark of hippocampal sclerosis (). The limbic system's physiologic ability to produce and propagate synchronized activity during normal cognition makes the limbic system an ideal environment for the propagation of pathological synchronization during a seizure (). Given the role of the limbic system in long-term recall and in emotion, it has been postulated that the high prevalence of memory, mood, and affective symptoms among patients with mTLE correlates with a disruption of normal limbic function ().

Networks throughout the brain are altered in epilepsy. Network changes have been associated with cognitive decline, seizure onset zone locations, and surgical outcomes (, ). For example, increased thalamic “hubness” (the importance of a node in a network) prior to anterior temporal lobectomy surgery is associated with increased risk of seizure recurrence after surgery (). Effective connectivity inferred by dynamical causal modeling has suggested that strengthened connectivity between the hippocampus and parahippocampal gyrus is associated with poor seizure control in TLE (). Despite promising associations between limbic circuitry and clinical features, the contribution of limbic connectivity to actual seizure burden remains unclear.

To identify resting-state functional connectivity (RS FC) features associated with seizure frequency, we assessed limbic and temporal RS FC in mTLE patients and healthy controls. We used network-based statistics to identify functional connectivity abnormalities in the limbic system, and to find a relationship between those abnormalities and seizure frequency. By identifying the networks underlying clinical presentations of epilepsy, we may better target neuromodulatory approaches to the circuit components either potentiating or regulating epileptic networks.

Materials and Methods

Participants and Clinical Scores

Consecutive mTLE patients who underwent comprehensive epilepsy evaluations including EEG monitoring at the Mayo Clinic Epilepsy Center were identified from an epilepsy research database. Twelve patients (five females) with unilateral mTLE (Table 1) and 12 healthy control subjects (Table 2), matched to the mTLE subjects by age, gender, and handedness, were studied. The control subjects were free of neurological and psychological disease. All participants provided written, informed consent in accordance with research protocols approved by the institutional review board of Mayo Clinic. We collected the numbers of seizures per month during the 3 months preceding the MRI scan as a quantitative score of disease burden. Structural MRI data were reviewed for all subjects to assess for neuroradiological abnormalities. A board certified psychiatrist (PEC) retrospectively reviewed all clinical records to ascertain if subjects had a co-occurring psychiatric disorder. In mTLE group, there were 8 subjects with neuroradiological abnormalities and 9 subjects with psychiatric disorders.

Table 1

Subject index (N = 12)Video-EEG diagnosis; seizure onset zoneHandednessAge range (year)Epilepsy duration (year)Seizure frequency (number of seizures per months)Neuroradiological abnormality in the structural MRI assessmentPsychiatric symptomPsychotropic medication
01Right mTLERight60–69500.1NoneDepression, AnxietyCitalopram
02Left mTLERight30–39168NoneRemote history of alcohol and cannabis abuse (8–10 years before presentation)None
03Left mTLERight40–49428Stable non-specific foci of increased T2/FLAIR signal in the subcortical left frontal lobe white matterNoneNone
04Right mTLERight20–292830Right mesial temporal sclerosisAnxietyCitalopram, sertraline
05Left mTLERight10–19164NoneInsomniaClonazepam
06Left mTLERight30–393224Left hippocampal atrophyDepression, ADHDMixed amphetamine (Adderall XR)
07Left mTLERight20–29123Mild leukoaraiosisAnxietyNone
08Left mTLELeft20–2972T2 hyperintensities suggestive of migraineDepression, nicotine-use disorderCitalopram
09Left mTLERight20–29120NoneNicotine-use disorder (smokeless tobacco)None
10Right mTLERight50–5923Right temporal encephaloceleNoneNone
11Left mTLELeft20–29221Left hippocampal atrophyNoneNone
12Left mTLERight30–39616Nonspecific T2 hyperintensitiesAnxietyAlprazolam

Demographic information for participants with medial temporal lobe epilepsy (mTLE).

Twelve healthy control subjects were matched for age, sex, and handedness (N = 12, mean age difference = 0.00; p > 0.99). The indirectly identifiable patient data (gender, exact age, exact seizure onset age) were removed according to editorial guidance.

Table 2

Subject index (N = 12)HandednessAge range (year)Neuroradiological abnormality in the structural MRI assessmentPsychiatric symptomPsychotropic medication
01Right50–59NoneNoneNone
02Right30–39NoneNoneNone
03Right40–49NoneNoneNone
04Right20–29NoneNoneNone
05Right10–19Possible pineal cystNoneNone
06Right30–39NoneNoneNone
07Right20–29NoneNoneNone
08Right20–29NoneNoneNone
09Right20–29NoneNoneNone
10Right50–59NoneNoneNone
11Left20–29White matter hyperintensitiesNoneNone
12Left20–29NoneNoneNone

Demographic information for healthy control subjects (five females).

Image Acquisition

Anatomical and functional MRI images were acquired in all subjects on a Siemens 3T Magnetom Skyra system using a 32-channel array head coil and tetrahedron-shaped foam pads to minimize head movement. High-resolution structural whole-brain images were acquired using a T1-weighted sequence with 0.5 × 0.5 × 1.2 mm3 resolution, TR = 2.3 s, TI = 0.9 s, TE = 1.96 ms, and FA = 9°. Subjects were instructed to keep their eyes open during which resting-state (RS) functional magnetic resonance imaging (FMRI) data were acquired by using a gradient echo-planar sequence sensitive to blood oxygenation level-dependent contrast with 3.28 × 3.28 × 3.3 mm3 resolution, 50 slices, TR = 2.9 s, TE = 30 ms, FA = 90°, and total acquisition time of 464 s.

Preprocessing of Imaging Data

Preprocessing of all imaging data was conducted using the Analysis of Functional NeuroImages (AFNI) software package (http://afni.nimh.nih.gov). The RS FMRI data were preprocessed and denoised by the standard protocol of the AFNI package (). Using the robust non-linear warping function of the AFNI package, all T1 anatomy data were registered to the MNI152-T1-2009c atlas of the Montreal Neurological Institute and linearly resampled in the 1 mm isocubic grid space. The EPI data were aligned to the T1 images, and then warped to the template brain space along with their T1 images. The registration results for the image data of all subjects were visually inspected for subcortical and cortical brain structures1. Note that the imaging data of three mTLE patients with seizure foci in the right hemisphere were left-right flipped before the preprocessing, to match the ipsi- and contra-lateral concept of the analysis. The terms “ipsilateral” and “contralateral” hemispheres stand, respectively for the hemisphere of seizure onset and the hemisphere opposite the side of seizure onset.

Functional Connectivity Analysis

For the functional network analysis (, ), masks for regions-of-interest (ROIs) in the limbic system and temporal lobe structures were defined by multiple atlases (), with reference to existing studies of animal seizure models and human patients (, ). The entire list of ROIs and the atlas details are presented in Figure 1A and Table 3. The RS FMRI time series were separately averaged in each ROI mask, and then a Pearson correlation matrix between those was calculated as the network data of each individual subject (by a permutation test for 10,000 iterations with random network extents). A total of 253 RSFCs between ROI pairs were calculated for each subject. The group difference graph between healthy control and mTLE groups were also determined by a two-sample t-test with a threshold level at the family-wise-error-corrected p < 0.01 ().

Figure 1

Table 3

ROI indexRegionBase atlasRemarks
HIPPOCAMPAL FORMATION
1Hippocampus (Hp)FS-HpCA1, CA2, CA3, and CA4 were merged into one mask.
2Subiculum (Sb)FS-HpSubiculum, parasubiculum, and subsubiculum were merged.
3Hippocampal tailFS-Hp
DIENCEPHALON
4Anterior nuclei of thalamusMorel-ThAnterior dorsal, medial, and ventral nuclei were merged.
5Mediodorsal nucleus of thalamus (ThMD)Morel-ThMagnocellular and parvocellular mediodorsal nuclei were merged.
6Hypothalamus (Hth)FS-asegHypothalamus is included in the inferior diencephalon mask.
SUBCORTEX
7AmygdalaFS-aseg
VENTROMEDIAL PREFRONTAL CORTEX (VmPFC)
8Brodmann area 11mMNI-VmPFC
9Brodmann area 14 (BA14)MNI-VmPFCBA14c, BA14m, BA14r, and BA14rr were merged.
CINGULATE CORTEX
10Dorsal anterior (dACC)MNI-VmPFCBrodmann area 24
11Ventral anterior (vACC)MNI-VmPFCBrodmann area 32
12Subgenual anterior (BA25)MNI-VmPFCBrodmann area 25
13Anterior middleFS-a2009s
14Posterior middleFS-a2009s
15Dorsal posterior (dPCC)FS-a2009s
16Ventral posterior (vPCC)FS-a2009s
17Subcallosal gyrusFS-a2009s
OTHER CORTICAL AREAS
18Parahippocampal gyrus (Php)FS-a2009s
19InsulaFS-a2009s
20Temporal poleFS-a2009s
21Superior temporal gyrus/sulcusFS-a2009s
22Middle temporal gyrusFS-a2009s
23Inferior temporal gyrus/sulcusFS-a2009s

Definition of regions-of-interest (ROIs) in limbic structures and seizure-relevant areas for the network-level analysis.

FS-HP, -aseg, -a2009s are FreeSurfer templates for hippocampal subfields, automatic segmentation, and cortical parcellation, respectively. Morel-Th is the reconstructed template image of Morel histological thalamic atlas, and MNI-VmPFC is the atlas for ventromedial prefrontal cortex of the Montreal Neurological Institute.

Correlation Analysis With Clinical Scores

We selected ROI pairs with significant group difference in their RSFC, and then performed a linear regression analysis to find correlations between functional connectivity of those selected ROI pairs and the seizure frequency for mTLE patients. R2 values were obtained by the linear regression function of MATLAB (version 2016a, The MathWorks, Inc., Natick, MA.).

Results

Functional Connectivity Difference Between mTLE and Healthy Control Groups

A total of 253 connections between ROI pairs were considered for group comparisons, and significant connectivity at the threshold level of pFWE−corrected < 0.001 was observed in 17 and 13% of the total ROI pairs in the healthy control and mTLE group, respectively (Figures 1A,B). The mTLE patients showed decreased functional connectivity in the 17 ROI pairs of the limbic system, relative to the healthy control group, by a two-sample t-test with the default threshold level of the NBS tool (pFWE−corrected < 0.01). Specifically, the mTLE groups showed diminished intra-limbic connectivity between ipsilateral hippocampus, subiculum, hypothalamus, mediodorsal thalamus, ventromedial prefrontal cortex (BA14), bilateral parahippocampal gyri, ventral anterior, dorsal and ventral posterior cingulate cortices, and contralateral dorsal anterior cingulate cortex (Figures 1C). The functional connectivity between limbic structures and temporal lobe regions showed less change compared to intra-limbic connectivity changes.

Correlations Between RSFC and Clinical Scores

Two out of 17 limbic ROI pairs, had a significant correlation between their RSFCs and seizure frequency at the threshold level of p < 0.05: ipsilateral ventral anterior cingulate cortex and hypothalamus (R2 = 0.391, p = 0.030) and ipsilateral subiculum and contralateral ventral posterior cingulate cortex (R2 = 0.362, p = 0.039) (Figure 1C). The results of group comparisons and their correlation with seizure frequency are summarized in the Table 4 and Figure 1D.

Table 4

Node pairt-valuep-value (two-tailed)Correlation with seizure frequency
Hippocampus (ipsilateral)- vACC (contralateral)−5.60< 0.0001Not significant
- vPCC (ipsilateral)−5.46< 0.0001Not significant
- vPCC (contralateral)−5.37< 0.0001Not significant
- dPCC (ipsilateral)−6.41< 0.0001Not significant
- dPCC (contralateral)−4.110.0005Not significant
- Ventromedial prefrontal cortex (BA14, ipsilateral)−3.960.0007Not significant
- Parahippocampal gyrus (contralateral)−3.300.0033Not significant
Subiculum (ipsilateral)- vACC (contralateral)−3.640.0014Not significant
- vPCC (ipsilateral)−3.650.0014Not significant
- vPCC (contralateral)−3.770.0011R2 = 0.391, p = 0.030
- dPCC (ipsilateral)−4.310.0003Not significant
Hypothalamus (ipsilateral)- vACC (ipsilateral)−3.720.0012R2 = 0.362, p = 0.039
- dACC (contralateral)−3.670.0013Not significant
Mediodorsal Thalamus (ipsilateral)- dACC (contralateral)−3.940.0007Not significant
Parahippocampal gyrus (ipsilateral)- vACC (contralateral)−3.110.0051Not significant
- dACC (contralateral)−4.82< 0.0001Not significant
Parahippocampal gyrus (contralateral)- vACC (contralateral)−3.360.0028Not significant

The ROI pairs (circles) in limbic system, showing significant difference in functional connectivity between mTLE and HC groups.

Discussion

In this study, we compared RS FC in patients with mTLE and healthy controls. Using network-based statistics, we identified connectivity changes associated with mTLE and with seizure frequency in mTLE. RS FC within the limbic system decreased in mTLE, more so than limbic—temporal connectivity. These network abnormalities primarily affected the hemisphere ipsilateral to the seizure focus. Seizure frequency correlated with RS FC between ipsilateral hypothalamus and ventral anterior cingulate cortex and between ipsilateral subiculum and contralateral posterior cingulate cortex.

Our findings are consistent with the observed decrease in functional connectivity between limbic areas in mTLE patients, which is particularly pronounced in the hemisphere ipsilateral to the seizure focus (, ). It is possible that the decrease in functional connectivity reflects underlying limbic pathology in the mTLE group through hypothalamus-to-mediodorsal-thalamus and hippocampus/subiculum-to-cingulate-cortex pathways (, ). These changes would be consistent with the affective, memory, and cognitive/behavioral symptoms experienced by mTLE patients ().

It is also possible that the decrease in limbic RS FC is compensatory. It may be speculated that decreasing functional connectivity ipsilateral to the seizure focus helps limit the participation of the abnormal temporal lobe in limbic circuitry, thereby preserving some measure of normal limbic function. The findings reported here for FMRI connectivity are consistent with studies using intracranial electrocorticography, which show that brain regions generating seizures are functionally isolated from surrounding brain regions (, ).

It is unclear why connectivity pairs in the cingulate cortex and hypothalamus and in the cingulate cortex and subiculum emerged as associated with seizure frequency. Previous network studies have focused on the clinical outcome of seizure freedom and not seizure frequency. One study of seizure propagation networks showed that seizure freedom associated with decreased functional connectivity in bilateral midline structures including the precuneus and midcingulate (). This finding, together with our observation that select cingulate connectivity features are associated with seizure frequency, may implicate midline connectivity as an overall indicator of limbic network health. Importantly, without larger studies pooling subjects with similar disease courses, this proposal is purely speculative.

With the limbic system acting as a centerpoint for a range of cognitive and seizure-related processes, broader network changes may help explain the role of specific limbic regions in epilepsy (). A pilot study of functional connectivity in deep brain stimulation (DBS) for epilepsy has suggested that patients who respond to DBS at the anterior nucleus of the thalamus have increased thalamic connectivity to the DMN (). Intrinsic DMN connectivity increases prior to inter-ictal epileptiform discharges in TLE before returning to its baseline, decreased connectivity after the discharge (). Further study of how paired connectivity features interact with other resting-state functional networks is warranted.

The primary limitation of this study was the small sample size. We did not conduct structured psychiatric interviews for this study and only retrospectively collected psychiatric histories from the medical record. Future efforts building on the present work will include prospective and comprehensive structured psychiatric assessments. A deeper understanding of the cause and nature of limbic connectivity in patients with mTLE may prove important in treating mTLE and the comorbidities associated with limbic network dysfunction. Also, we flipped the EPI data for left-handed subjects, but there are reports showing that mTLE patients have different behaviors regarding to the anatomical and functional features on the laterality (). To observe more details on the laterality of mTLE connectivity, the further investigation with larger cohorts is required.

In conclusion, we identified functional connectivity changes associated with mTLE and cingulate and limbic network connectivity abnormalities that correlated with increased seizure frequency.

Statements

Ethics statement

All participants provided written, informed consent in accordance with research protocols approved by the institutional review board of Mayo Clinic.

Author contributions

HJ, KW, IB, and GW: design and conceptualization of the study. HJ, DK-J, PC, EB, and GW: analysis and interpretation of the data. HJ, IB, and DK-J: drafting the manuscript for intellectual content. PC, DJ, EB, GW, and IB: revising the manuscript for intellectual content. KW and DJ: major role in the acquisition of data.

Acknowledgments

The Thalamus Atlas is copyrighted by © University of Zurich and ETH Zurich, Axel Krauth, Rémi Blanc, Alejandra Poveda, Daniel Jeanmonod, Anne Morel, Gábor Székely. The authors also acknowledge the Mayo High Performance Research Computing Facility (RCF) at Mayo Clinic in Rochester, MN, for providing computational resources.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Footnotes

1.^Aside from the visual inspection to assess thalamic atrophy, recent studies using quantitative MRI analyses have reported global as well as segmental thalamic atrophy in TLE [see () for review].”

References

Summary

Keywords

medial temporal lobe epilepsy, partial seizure, functional magnetic resonance imaging, network-based statistics, limbic system

Citation

Jo HJ, Kenney-Jung DL, Balzekas I, Welker KM, Jones DT, Croarkin PE, Benarroch EE and Worrell GA (2019) Relationship Between Seizure Frequency and Functional Abnormalities in Limbic Network of Medial Temporal Lobe Epilepsy. Front. Neurol. 10:488. doi: 10.3389/fneur.2019.00488

Received

06 November 2018

Accepted

23 April 2019

Published

08 May 2019

Volume

10 - 2019

Edited by

Andrea Romigi, Mediterranean Neurological Institute (IRCCS), Italy

Reviewed by

Marino M. Bianchin, Federal University of Rio Grande do Sul, Brazil; Luiz Eduardo Betting, São Paulo State University, Brazil

Updates

Copyright

*Correspondence: Hang Joon Jo Gregory A. Worrell

This article was submitted to Epilepsy, a section of the journal Frontiers in Neurology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics