Review Article - Imaging in Medicine (2013) Volume 5, Issue 3

Quantitative analysis of structural neuroimaging of mesial temporal lobe epilepsy

Negar Memarian1, Paul M Thompson1,2, Jerome Engel Jr1 and Richard J Staba*1

1Department of Neurology, Reed Neurological Research Center, Suite 2155, University of California, 710 Westwood Plaza, Los Angeles, CA 90095, USA

2Laboratory of Neuroimaging, University of California, Los Angeles, CA, USA

*Corresponding Author:
Richard J Staba
Department of Neurology, Reed Neurological Research Center
Suite 2155, University of California, 710 Westwood Plaza
Los Angeles, CA 90095, USA
Tel:
+1 310 825 8479
Fax:
+1 310 206 8461
E-mail:
rstaba@mednet.ucla.edu

Abstract

Keywords

drug-resistant epilepsy; DTI; mesial temporal lobe epilepsy; morphometry; MRI; quantitative neuroimaging

Epilepsy is a chronic brain disorder characterized by an enduring predisposition to generate spontaneous epileptic seizures. Epilepsy affects nearly 3 million Americans, making it the third most common neurological disorder in the USA. Worldwide, an estimated 50 million individuals are affected by epilepsy, which accounts for 1% of the global burden of disease [1,2]. Between 60 and 80% of people with epilepsy will achieve seizure control with antiseizure drugs [36]. However, these statistics imply that as many as 40% of patients with epilepsy have seizures that are not adequately controlled by antiseizure drugs. The International League Against Epilepsy has proposed that drug-resistant epilepsy is a failure of two tolerated, appropriately chosen antiseizure drug trials to achieve sustained seizure freedom [7]. Research indicates that only a small percentage (<10%) of individuals with epilepsy benefit from subsequent drug trials after failing the first two [8,9]. It is not known why some seizures are or become resistant to medication, but several features are frequently associated with pharmacoresistant seizures [1013]. Among them, the most common pathology and one commonly associated with pharmacoresistant limbic seizures is mesial temporal or hippocampal sclerosis (HS) [14,15]. The classic pattern of HS described by Bratz is associated with significant neuron loss and gliosis in subfield CA1 and prosubiculum [16], as well as the area between blades of dentate gyrus or ‘end folium’ [17]. There is less damage to dentate gyrus granule cells, CA3 and particularly CA2 pyramidal cells, and relative preservation of cells in subicular and parahippocampal gyrus. Patients with unilateral HS, which can be detected using MRI epilepsy protocols, have seizures that arise from or involve the affected mesial temporal lobe (MTL) structures that correspond to characteristic clinical signs and symptoms [18,19]. HS is often associated with widespread bilateral limbic and neocortical disturbances [20]. It may also be present with other lesions including, but not limited to, heterotopia of the temporal lobe, cortical dysplasia, cavernous angioma, tumor, contusion and cerebral infarctions [21,22]. The presence of these features, family history of epilepsy (genetic) and prolonged febrile seizures in infancy are consistent with MTL epilepsy (MTLE) with HS.

Accurate diagnosis is critical as drug-resistant epilepsies, such as MTLE with HS, and those with well-circumscribed epileptogenic lesions, known pathophysiology and predictable natural history can be treated successfully with surgery [8]. Successful surgical outcome – that is, seizure freedom or significant reduction of disabling seizures – depends on accurately delineating the epileptogenic zone that theoretically represents the brain areas necessary and sufficient for generating spontaneous seizures. The epileptogenic zone cannot be measured directly, but is inferred from presurgical diagnostic tests, such as video-EEG monitoring using scalp electrodes or, in some cases, intracranial grid or depth electrodes and neuroimaging. MRI is the preferred imaging modality to identify structural abnormalities responsible for the generation of spontaneous seizures confirmed by electrophysiological studies – that is, epileptogenic lesions. Conventional MRI does not reliably capture epileptogenic abnormalities such as HS or many types of malformations of cortical development; however, detection improves dramatically with an epilepsy MRI protocol and neuroradiologists who are knowledgeable about structural lesions that cause epilepsy [23]. In addition, quantitative analysis of MRI can aid in the detection of structural lesions and, since MRI plays such an important role in the diagnosis and management of epilepsy, some have proposed that MRI should be included in the classification of epilepsy etiologies [24].

Technical advances in neuroimaging have spurred the development of sophisticated analysis techniques. These techniques have been used chiefly in the basic research of epilepsy and have provided new information on structural abnormalities associated with drug-resistant epilepsy and MTLE with HS in particular. Evidence from structural MRI, along with electrophysiological and histological data and postsurgical seizure freedom, indicate that there is not one type of HS, and the ‘hippocampal-centric’ view of MTLE is being revised with greater emphasis on networks that include mesial temporal and extratemporal limbic structures. The purpose of this article is to review several MRI analysis techniques that have been used to quantify the spatial distribution and extent of structural brain abnormalities. The authors discuss the basis of voxel-based morphometry (VBM), surface-based morphometry (SBM) and pattern- based morphometry as well as cortical pattern matching (CPM), and the advantages and disadvantages of each. The focus is on the application of these techniques and detection of MRI structural abnormalities in patient studies of MTLE. Finally, the authors comment on the future of MRI morphometry and implications for basic research and clinical studies of drugresistant epilepsy.

Review of the literature

Quantitative structural imaging of epilepsy

The detection of pathological structural substrates that cause epilepsy has increased significantly with the advent of noninvasive, high-resolution MRI [22]. According to the International League Against Epilepsy Neuroimaging Commission, the essential indications for MRI include partial or secondarily generalized seizures and generalized seizures that do not remit with antiseizure drug treatment, and development of progressive neurological or neuropsychological deficits [25]. Furthermore, the Neuroimaging Commission recommends that imaging of epilepsy should include T1- and T2-weighted sequences to cover the whole brain in a minimum of two orthogonal planes, with a slice thickness of 1.5 mm or less to allow reformatting in any orientation and 3D reconstruction [25]. In patients with HS as the suspected epileptogenic lesion, images are acquired in an oblique coronal plane, perpendicular to the anterior–posterior axis of the hippocampus. The presence of HS on T1-weighted imaging generally appears as shrunken hippocampal gray matter (GM) and often loss of internal architecture ipsilateral to the site of seizure onset. In some cases, GM loss occurs bilaterally, but more extensive alterations are often found in the ipsilateral rather than the contralateral hippocampus. Assisting in the detection of HS are other MRI scans, such as T2-weighted and fluid-attenuated inversion recovery (FLAIR) sequences, which voids the water signal from ventricular and subarachnoid cerebral spinal fluid (CSF) and makes it easier to see signal changes associated with pathology in periventricular GM, such as in the hippocampus. Studies have found patients with histologically verified HS or neocortical lesions identified on conventional MRI correspond to increased signal intensity on T2-weighted and FLAIR sequences [26,27]. These studies and others that correlated MRI measures (e.g., T2 relaxation times) with neuronal and glial cell density suggest that increased signal intensity reflects increased interstitial fluid due to neuron loss and possibly gliosis [2830]. Recent work, however, has not found a correlation between T2 relaxation time or normalized FLAIR intensity and glial immunostaining [30,31].

Advanced MRI analysis techniques that are described in the sections that follow can quantify the extent of GM loss, which is necessary in cases when HS asymmetry is subtle, and can also display the surface contours and spatial patterns of GM loss in hippocampus and other brain structures. Based on imaging studies of MTLE, we now know that GM loss extends beyond the affected hippocampus and includes cortical GM and white matter (WM) alterations in ipsilateral and contralateral temporal lobe and often extratemporal lobe structures.

Voxel-based morphometry techniques

In VBM, whole-brain digital structural images are compared on the basis of properties corresponding with groups of voxels. Voxel-based methods have been implemented in many different ways (e.g., [3235]), but the general premise is to establish voxel-for-voxel correspondence across subjects through nonlinear registration of multiple subjects’ brain images to a standard anatomical template. Once the images are aligned into the same coordinate space, one can perform voxel-wise statistical comparisons of GM and WM volume or concentration across different subject groups [36]. GM volume can be derived from the number of GM voxels in the brain and GM concentration is the amount of GM per unit of intracranial volume. The concentration at a voxel is the volume of the GM in a neighborhood around that voxel divided by the overall volume of the structures in the same region around that voxel [37]. The standard VBM process typically involves four steps: spatial normalization, tissue segmentation, spatial smoothing and statistical analysis [38]. Spatial normalization involves applying a nonlinear registration to each subject’s T1-weighted MR image, so that local areas (voxels corresponding to anatomical brain areas) stretch or compress with respect to each other to match a group template. The deformed image is then segmented into tissue classes (i.e., GM, WM and CSF) based on the intensity in the image as well as a predetermined probability that a particular type of tissue will be found at a given location [37]. The segmented image is then spatially smoothed. The intensity in each voxel of the smoothed image is a local weighted average of GM, WM or CSF from adjacent voxels, which is generally expressed as a GM, WM or CSF concentration [38]. The final step of VBM involves voxel-wise statistical analysis that commonly includes group comparisons or correlations with covariates of interest [33,38].

The initial implementation of VBM was modified to avoid problems that arose when data from multiple subjects were not accurately aligned into the same coordinate space. In some cases, group differences in tissue volumes were inferred when images were not fully aligned across subjects and groups [33]. This occurs because the parameters of the normalization only encode smoothed, low-frequency deformations that may not fully align GM and non-GM structures (e.g., ventricles). Optimized VBM methods often involve an additional step, which multiplies the spatially normalized GM concentration (or other tissue class) by its relative volume before and after spatial normalization. This is referred to as modulation and enables comparisons of voxel-wise GM volume differences between groups [33].

Voxel-based methods are computationally efficient as they work directly on the voxel grid. Furthermore, the smoothing step in VBM methods suppresses confounding variations due to the high variability of gyral anatomy between individuals [38]. Hence, VBM techniques are useful in studies that quantitatively evaluate the spatial distribution and extent of statistically significant GM loss (or gain) in one subject group with respect to another, for example, patients compared with age- and sex-matched controls [36,39]. VBM has been successfully used to characterize structural abnormalities in patients with Parkinson’s disease [40], Alzheimer’s disease [39], schizophrenia [41] and, of course, epilepsy [38].

Studies of MTLE using automated VBM analysis consistently find GM loss in the hippocampus ipsilateral to – and, to a lesser extent, contralateral to – the site of seizure onset in patients with MRI evidence of unilateral HS [4245]. VBM-based GM loss has also been detected in the amygdala, entorhinal cortex and parahippocampal gyri [46,47]. A review of the studies cited above and others indicates that GM reductions extend beyond the sclerotic MTL ipsilateral to seizure onset and often include other limbic structures, such as the thalamus, cingulate gyrus and orbitofrontal cortex, as well as areas of temporal and parietal neocortex, striatum and cerebellum [38]. Comparable results have been found using standard or optimized VBM methods, and optimized VBM may better detect subtle structural abnormalities [46]. Some studies observed that GM loss was greater in patients with MRI evidence of HS than those without [48,49], and VBM was effective in detecting GM alterations associated with cortical dysplasia that can accompany HS [50]. In other studies, GM reductions were more obvious in patients with left compared with right MTLE [5153], although one study detected more prominent GM loss in patients with right-sided MTLE [50]. In addition to GM loss, several studies found WM loss adjacent to MTL structures ipsilateral to seizure onset and in some cases in the contralateral temporal lobe, as well as bilateral loss within the frontal lobe [43,45,5456]. An analysis of GM volume and postsurgical seizure freedom in a cohort of patients with left MTLE and histological evidence for HS found patients who were not seizure free had reduced GM volumes in ipsilateral posterior MTL and the contralateral hippocampus compared with seizure-free patients [57].

Surface-based morphometry techniques

In SBM, morphometric measures are derived from geometric models of the cortical surface [37] or surface-based models of subcortical structures, such as the hippocampus. SBM techniques are used to reconstruct and analyze surfaces that represent structural boundaries within the brain. The major steps involved in SBM of the brain are preprocessing (image registration, intensity normalization, smoothing, brain extraction and brain segmentation), surface extraction, surface inflation or parameterization and surface mapping [58]. In order to align individual brain images and allow for comparisons across subjects, in cortical surface analyses, the brains are mapped to a unit sphere on which their original properties (e.g., cortical thickness) can be compared with each other, and results are mapped back to a reference brain surface (Figure 1). Several implementations of SBM have been made (e.g., [5862]) with the FreeSurfer software from Harvard Martinos Center for Biomedical Imaging (MA, USA) being the most widely used (available through [201]). FreeSurfer is a set of automated tools for the reconstruction of the brain’s cortical surface from structural MRI data. It also allows the overlay of functional MRI data onto the reconstructed surface.

imaging-in-medicine-intensity-normalized

Figure 1. Surface-based morphometry processing steps. (A) High-resolution T1-weighted MRI scan. (B) Computerized reconstruction of the gray/white matter boundary. A smoothed and expanded view of the white matter surface is shown. The image has been intensity normalized, skull-stripped and the cerebellum has been removed. (C) Inflation of the cortical surface to map gyral and sulcal anatomy. (D) Coregistration of the subject’s cortical surface to a common spherical template. This step allows the assessment of cortical tissue properties with respect to a normative database using a common coordinate system. Reproduced with permission from [147]. For color images please see online www.futuremedicine.com/doi/pdf/10.2217/ iim.13.28.

Studies of MTLE using custom [63,64] and commercial SBM-based algorithms, such as FreeSurfer, have found significant hippocampal GM volume loss ipsilateral and, in some cases, contralateral to seizure onset [65,66]. However, in other studies, automated SBM analysis underestimated the extent of hippocampal GM loss compared with manual segmentation methods [67,68]. Patient studies of MTLE with HS found reduced ipsilateral hippocampal GM volumes correlated with cortical GM thinning in the temporal lobe [66,69]. Cortical GM was on average 5–15% thinner in bilateral frontal (prominently within pre- and para-central gyri operculum) and temporal lobes (Heschl’s gyrus) of patients with MTLE compared with controls [70], with somewhat greater thinning in patients with left- versus right-sided MTLE [71]. Furthermore, reduced thalamic GM volume extending from anterior to posterior primarily along the medial surface correlated with reduced hippocampal volume (subfield CA1) and MTL GM thinning [72], as well as bilateral frontal–central and lateral temporal lobe areas [73]. One study evaluated the presurgical presence of MRI hippocampal atrophy and postsurgical seizure freedom and found that, in cases with hippocampal atrophy, patients with postsurgical seizures had ipsilateral temporopolar and bilateral insular atrophy compared with patients who were seizure free [74]. Interestingly, in cases without hippocampal atrophy, postsurgical seizures were associated with ipsilateral posterior lateral temporal and contralateral parietal–occipital lobe atrophy.

Other cortical morphometry techniques

CPM is related to SBM, but differs in that after creating explicit geometric models of the cortex using parametric surfaces, deformation maps are built on the geometric models to spatially align cortical regions across subjects [75,76]. This is done by applying mathematical transformations that align key anatomical landmarks from one dataset to another. CPM techniques usually involve three main steps of cortical parameterization, matching cortical features across individuals and statistical comparisons with map group differences in features such as cortical thickness or GM concentrations. The first step involves creating geometrical models of the cortical surface and includes several steps, such as image segmentation, registration and 3D reconstruction. Matching cortical features across individuals is performed through neuroanatomical labeling (e.g., marking sulcal lines) and warping one brain surface onto another. Finally, statistical testing is performed to investigate the effect of various factors (e.g., disease, aging, treatment and genetics, among others) on various cortical attributes, such as GM thickness or gyrification [77]. These steps are shown in Figure 2 from a study of GM thickness differences between MTLE patients and healthy controls.

imaging-in-medicine-epilepsy-patients

Figure 2. Cortical pattern matching steps (facing page). The images are from the authors’ study of gray matter thickness differences between mesial temporal lobe epilepsy patients and healthy controls currently in progress, which follows the cortical pattern-matching protocol suggested by Thompson et al. [75]. (A) Segmentation of the region of interest from T1-weighted MRI, (B) 3D reconstruction and neuroanatomical labeling (the following sulci are traced on the lateral surface: sylvian fissure, central, postcentral, precentral, superior temporal main body, superior temporal ascending branch, superior temporal posterior branch, intraparietal, primary intermediate, secondary intermediate, transverse, inferior temporal, inferior frontal and superior frontal; on the bottom surface: olfactory, occipital–temporal and collateral; on the medial surface: callosal, callosal inferior, superior rostral, inferior rostral, paracentral, cingulate anterior, cingulate posterior, anterior cingulate outer, parieto–occipital, calcarine anterior, calcarine posterior and subparietal), (C) gray matter thickness measurement, (D) generation of cortical thickness map and (E) generation of statistical maps depicting significant differences in gray matter thickness between patients and controls. For color images please see online www.futuremedicine.com/doi/pdf/10.2217/iim.13.28.

Similar to the studies using VBM and SBM methods described in the preceding sections, CPM techniques have been used to evaluate neocortical GM thickness in patients with MTLE and HS. In one study, GM thickness was reduced by up to 30% and reductions in cortical complexity were found in patients compared with controls [20]. In this latter study, significant GM loss was observed bilaterally in frontal areas (pole, operculum and orbitofrontal areas), lateral temporal and occipital regions independent of the hemisphere of seizure onset. In patients with HS, analyses of hippocampal surface structure found significant alterations on the lateral hippocampal surface corresponding to Sommer’s sector, while patients without HS had surface changes on the superior surface of the hippocampal body, in right-side MTLE, and on the superior surface of the medial hippocampal head region, near the uncinate gyrus, in left-sided MTLE [78]. Another study evaluated the spatial distribution of hippocampal atrophy in relation to different depth electroderecorded ictal EEG onset patterns [79]. In this study, based on EEG morphology, propensity for seizure spread and ipsilateral histological features, it was hypothesized that seizures beginning with hypersynchronous and low-voltage fast EEG patterns reflect different neuronal and anatomical mechanisms of generation. Results of CPM analysis found that patients with hypersynchronous ictal onsets had significant hippocampal atrophy ipsilateral to seizure onset that resembled a classic pattern of HS (sparing of subfield CA2; Figure 3), while patients with lowvoltage fast ictal onset had diffuse ipsilateral hippocampal atrophy and significant atrophy in contralateral hippocampus. Interestingly, in a separate study, surgical patients with MTLE who were not seizure free had atrophy that was more diffuse in the ipsilateral hippocampus and significant atrophy in anterior and lateral areas of the contralateral hippocampus compared with patients who were seizure free [80].

imaging-in-medicine-contour-maps

Figure 3.Structural differences in ipsilateral (epileptogenic) hippocampus of patients with a hypersynchronous and low-voltage fast seizure onset pattern. (A) Superior aspect of ipsilateral (epileptogenic) hippocampus of patients with HYP and LVF seizure onset pattern. (B) Inferior aspect of ipsilateral (epileptogenic) hippocampus of patients with HYP and LVF seizure onset pattern. Surface contour maps depicting areas of local atrophy (P maps) show regions of significant differences in HYP (n = 8) and LVF (n = 9) patients relative to controls. Areas with p < 0.05 indicate significant atrophy. HYP: Hypersynchronous; LVF: Low-voltage fast. Reproduced with permission from [79]. For color images please see online www. futuremedicine.com/doi/pdf/10.2217/iim.13.28.

Due to the complex folding of the human MTL, it is difficult to visualize the hippocampal formation (dentate gyrus, hippocampus subfields CA1-3 and subicular cortex) together with adjacent structures (entorhinal cortex, parahippocampal gyrus and fusiform gyrus). Advances have been made with a novel cortical unfolding technique using high-resolution T2-weighted MRI that increases signal-to-noise ratio and resolution of MTL structures [81,82]. Figure 4 illustrates the main steps in the unfolding procedure that includes the standard manual segmentation of GM, WM and CSF, and an interpolation step to generate a 3D GM ribbon corresponding to MTL structures. The 3D GM ribbon is then computationally unfolded and flattened into 2D and structural boundaries projected onto the flattened image to produce an anatomical map of the adjoining MTL subregions. Similar to other morphometry techniques, GM thickness can be computed for each of the MTL subregions and individual maps of the MTL warped into a common plane with maps from other subjects to compute group averaged GM thickness maps that retain much of the anatomical variability within the subject group. In addition, individual maps provide the capability to localize and characterize functional MRI signal changes simultaneously across MTL subregions [81,83]. Furthermore, MTL maps generated from MRI of presurgical patients can be used to locate the position of multiple microelectrodes extending beyond the distal tip of clinical depth electrodes (Figure 4F; red circles depict location of microelectrodes) [82]. One study applied these techniques to evaluate interictal MTL single neuron firing properties in relation to subregion GM loss computed in patients with respect to controls [84]. Results from this study found significant bilateral MTL GM loss, although, compared with MTL contralateral to seizure onset, there was a stronger correlation between GM loss and increased neuronal firing in the ipsilateral MTL. Adjusting for levels of GM loss revealed significantly higher neuronal firing and bursting in ipsilateral than contralateral MTL, suggesting that synaptic reorganization following cell loss is associated with varying degrees of epileptogenicity.

imaging-in-medicine-temporal-lobe

Figure 4.MRI-based mesial temporal lobe gray matter thickness maps. The unfolding technique for left mesial temporal lobe (MTL) viewed in the coronal plane is demonstrated, but methods were carried out bilaterally in each patient’s MRI. (A) 3-Tesla T2-weighted MRI was acquired perpendicular to the anterior–posterior axis of the hippocampus. (B & C) MTL gray matter (shaded in (C)) was segmented manually by excluding cerebrospinal fluid and white matter (both shaded). (D) MTL gray matter strip was mathematically up-sampled, expanded and smoothed to generate a 3D gray matter ribbon. (E) A rule-based protocol and atlases of histologic and structural hippocampal anatomy were used to draw MTL subregion margins. The blocks beginning from the top right and moving down and across denote margins between CA23DG and CA1, CA1 and Sub, Sub and EC, EC and PRC, and PRC and fusiform gyrus. (F) Computationally unfolded 2D map that contains margins corresponding to MTL subregions. CA23DG: Cornu ammonis subfields 2 and 3/dentate gyrus; CADG: Cornu ammonis/dentate gyrus; EC: Entorhinal cortex; PHG: Parahippocampal gyrus; PRC: Perirhinal cortex; Sub: Subiculum. Reproduced with permission from [84]. For color images please see online www.futuremedicine.com/doi/pdf/10.2217/iim.13.28.

Pattern-based techniques

A more recent approach to study brain structure involves data-driven techniques that use machine learning to extract specific patterns [85], features [86] or objects [87] in MRI. Unlike most morphometry methods that assume one-to-one correspondence of anatomy across subjects, pattern-based methods identify distinctive anatomical patterns (e.g., appearance and geometry) that could be due to disease and, thus, are only present in subsets of subjects. These features can be automatically extracted from a set of subject images and their relationship to diagnostic categories can be modeled using supervised learning. Pattern-based techniques seek to identify instances of the same image feature in different subjects and discover features that tend to co-occur with specific subject groups [86]. Pattern-based techniques have been clinically validated for identification of structural differences between normal and Alzheimer’s disease patients [85,86]. They have also been used for gyrus-based parcellation of the cortical surface [87,88]. However, the pattern-based approaches have not yet been used in epilepsy studies.

Alternative MRI modalities

Lesions may not be obvious on conventional T1-weighted MRI in as many as 20% of patients with pharmacoresistant seizures [89], although alternative techniques, such as T2-weighted MRI, FLAIR and other MRI sequences, including double inversion recovery [90,91] and magnetization transfer ratio [90], can detect structural abnormalities in approximately half of these patients [92]. In addition, a special type of MRI is sensitive to the diffusion of water using pulsed field gradients with a method called diffusionweighted imaging. The extent of water movement freedom is quantified in a parameter called apparent diffusion coefficient (ADC). Several retrospective diffusion-weighted imaging studies found increased ADC values in the temporal lobe including the hippocampus and temporal pole [9396] and posterior aspects of corpus callosum [97]. In many studies, ADC values were appropriately lateralized to the site of seizure onset in a large percentage of patients with MTLE and HS, although measures of ADC alone and some of the MRI methods described above had lower specificity for the seizure onset zone in patients without HS [98100].

A special type of diffusion-weighted imaging is diffusion tensor imaging (DTI), which provides several useful measures of the restricted movement of water in tissue. Parameters derived from DTI include mean diffusivity, which reflects the overall rate of water diffusion in a voxel, and fractional anisotropy (FA), which measures how directionally constrained the diffusion process is, and is an indirect measure of WM (fiber) density or myelination levels. By following the principal directions of water diffusion, tract tracing methods can be used to follow neural pathways throughout the brain, and infer axonal connections. DTI studies have mapped the patterns of connectivity between the parahippocampal gyrus and orbitofrontal areas, as well as direct connectivity between the parahippocampal gyrus and the hippocampus itself (Figure 5A) [101]. The connections between neocortical areas and the hippocampus via the parahippocampal gyrus may provide a structural basis for theoretical models of seizure propagation. In patients with MTLE, voxel-based analysis of DTI has found increased mean diffusivity in the ipsilateral hippocampus, parahippocampal gyrus and thalamus [102104]. Reduced FA was observed in ipsilateral temporal lobe including hippocampus and parahippocampal gyrus, areas of frontal lobe and portions of corpus callosum [105107]. In several studies, reductions in FA were more prominent ipsilateral and, in some cases, contralateral to seizure onset in left compared with right MTLE [105,108,109]. In patients without lesions on MRI, reduced FA was found primarily in ipsilateral temporal lobe, corpus callosum and thalamus [104,110]. DTI studies using tractography found reductions in fornix bilaterally and cingulum [111,112], as well as increased mean diffusivity in ipsilateral WM fibers of the uncinate, arcuate and inferior longitudinal fasciculi, with the latter changes correctly lateralized to the site of seizure onset in 91% of patients with MTLE and hippocampal atrophy [111]. DTI has also been utilized to determine the prospective functional deficit zone in MTLE surgery candidates [113,114].

imaging-in-medicine-extrastriate-occipital

Figure 5.Diffusion tensor imaging for imaging structural connectivity. (A) Averaged connectivity maps from the left parahippocampal gyrus in ten controls, superimposed on the normalized single T1 images from the MRI provided by statistical parametric mapping 99 (sagittal views). These maps are displayed without applying any threshold to the connectivity values and display connectivity between the parahippocampal gyrus and the anterior temporal lobe, orbitofrontal areas, posterior temporal lobe and extrastriate occipital lobe via the lingual and fusiform gyri. (B) Results from probabilistic tractography obtained from different regions of interest averaged over all subjects. Each 3D brain rendering illustrates the voxel-wise probabilistic map of connectivity from one region of interest (top left box: left anterior cingulate; top right box: left inferior orbital region; bottom left box: left hippocampus; and bottom right box: left thalamus). (A) Reproduced with permission from [92] and (B) reproduced with permission from [127]. For color images please see online www.futuremedicine.com/doi/pdf/10.2217/ iim.13.28.

Discussion

None of the structural imaging techniques reviewed here have been used in the clinical presurgical evaluation of epilepsy, although each of the imaging techniques mentioned has strengths and weaknesses that are currently suitable for different research applications. VBM is relatively simple to implement, free and widely used in MRI studies of MTLE. While VBM is generally accurate for measuring GM atrophy in structures below the pial surface (e.g., hippocampus [38,44,115,116], amygdala [117] and thalamus [38,48,118]) and in some studies cingulate gyrus [118,119], it does not consistently detect GM abnormalities in the neocortex. This is due to the wide variations in cortical patterning across subjects, which makes it challenging to automatically align anatomical landmarks and cortical regions across subjects. In addition, it is difficult to detect significant differences in GM associated with neocortical pathology that is subtle or spatially complex [120]. SBM, on the other hand, is well-suited for detecting neocortical GM abnormalities because images are aligned with respect to neocortical folding patterns, and this improves neocortical registration compared with aligning on the basis of voxel intensity [121]. SBM methods also have a somewhat higher reproducibility than VBM and can be more robust in longitudinal studies that measure changes in GM thickness over time [36]. CPM uses high dimensional registration to provide detailed measurements of GM thickness in the hippocampus or neocortex of patients with MTLE. CPM is ideal for capturing small differences in GM thickness in patients without obvious pathology on MRI, although measures of volumetric GM loss can be sensitive to the geometrical calibration of the MRI scanner [121]. Furthermore, CPM can be difficult to implement and CPM algorithms are not as readily available or easy to use as VBM or SBM algorithms. Based on the differences in the design and application between these analysis techniques, it is difficult to specify one as the gold standard and use it in validation studies of other analysis techniques. The most reliable means of validation for these methods are correlation studies that relate measures of GM thickness to histological measures of neuronal and glial cell density.

Results from imaging analysis studies reflect statistically consistent differences across a sample of patients with MTLE typically with respect to control subjects. However, it is not known how accurate group-averaged data correspond to the individual patient, which is the basis for clinical decisions. More studies are required to determine the implication of these findings in the surgical outcome of patients who undergo epilepsy surgery. So far, several studies have shown that VBM in particular may not be a suitable standalone technique for detecting or spatially characterizing focal lesions in individual patients [122124]. Due to the low specificity of the findings, caution is especially needed when evaluating MRI-negative patients for localization of the epileptogenic zone [99]. In a VBM-based analysis of GM in patients with left MTLE and HS and other MRI-negative patients, GM abnormalities were detected in only 20% of patients with HS and 20% of MRInegative patients [124]. More recently, VBM analysis on single case studies that compared single subjects with a control group found a very high false-positive rate [123]. This can be because any significant difference may be driven by individual variability in neuroanatomy rather than the neuropathology of the disease under investigation, or may represent a false positive due to the data being sampled from non-normally distributed populations [123].

Supervised learning algorithms are a potential strategy to incorporate quantitative MRI techniques in surgical decision-making for individual patients. Supervised learning algorithms can be trained by a set of multimodal features (e.g., EEG, MRI and clinical information) for each patient and their corresponding surgery outcome result. These algorithms seek to find patterns between features that are exclusively associated with the given outcome. Usually, such patterns are not obvious from visual inspection of the data. After the algorithm learns from a group of patients, it can then classify a new patient (i.e., predict surgery outcome) based on his/her multimodal features. Examples of this approach in MTLE surgical decision-making can be found in [125] and [126]. One retrospective study that used support vector machine DTI classification achieved high diagnostic accuracy in lateralizing the seizure onset zone in individual subjects [115]. Other studies using novel computer-automated SBM and DTI algorithms have produced results comparable with manual, labor-intensive methods [63,68,108]. Ultimately, these strategies could be used during presurgical evaluation to help plan the placement of intracranial grid or depth electrodes to determine sites of seizure onset or provide information on fiber pathways prior to resective surgery to minimize the risk of cognitive, motor or visual-field impairment [92].

Considering the strengths and weaknesses of the MRI analysis techniques discussed in the preceding paragraphs, it appears that significant structural abnormalities extend beyond the epileptogenic hippocampus in patients with MTLE and HS. In addition, MRI studies have found factors such as hemisphere of seizure onset, duration of epilepsy, age of onset, history of febrile seizures and possibly extent of HS and comorbidity to explain some of the variability in spatial patterns of extra-hippocampal GM loss, as these could be important for identifying brain areas consistently damaged in patients with MTLE [42,59,79,80]. Recent work has used graph theoretical analysis to evaluate anatomical networks and results are consistent with the hypothesis that MTLE is a network disease of reorganized mesial temporal and extratemporal limbic structures (Figure 5B) [127,128]. What is not known is whether all or only some areas of the reorganized network reflect the epileptogenic zone. Remote GM loss found in many of the MRI studies could reflect, for example, deafferentation from the MTL where seizures begin [118], excitotoxic neuron loss from propagated seizure activity, developmental abnormalities that contribute to alterations in extra-hippocampal connectivity and structure [129], GM loss due to an initial injury and/or progression of the damage caused by the initial injury, or chronic exposure to antiseizure drugs [20]. Some of these structural alterations may be epileptogenic, while others may not. The structural analysis techniques mentioned here, such as cortical unfolding or hippocampal CPM, that can localize electrophysiological or functional MRI signals and relate it to measures of GM thickness could help distinguish epileptogenic from nonepileptogenic lesions [81,83,130].

Determining the boundaries of the epileptogenic zone has clear implications in the surgical treatment and prognosis for postsurgical seizure freedom. A review of epilepsy surgery found that the percentage of surgical patients with HS who were seizure free was between 48 and 84% [131]. This latter review and others cite several factors associated with seizure freedom [18], including subtypes of HS. A histological study by Blumcke and colleagues found that patients with less severe HS (e.g., end folium sclerosis) were more likely to have postsurgical seizures compared with patients with the classic pattern of HS who were more likely to be seizure free [132]. A technical explanation for these histological results could be that MRI does not reliably detect the full spatial distribution of damage associated with more subtle forms of HS. Results from a VBM study found residual lesions posterior to ipsilateral hippocampus and/or in the contralateral hippocampus in patients who received anterior temporal resection that were presumed to be the basis for persistent postsurgical seizures [57]. By contrast, a recent randomized trial for early resective surgery found that, in a carefully selected group of patients with unilateral MTLE, neuroimaging indicative of HS, and an absence of MRI and EEG evidence for remote epileptogenic lesions, the percentage of seizure-free patients was 85% at 2 years after surgery [133]. Although histology did not confirm the subtype of HS in this trial, the high percentage of postsurgical seizure freedom is consistent with a homogenous cohort of classical HS.

Conclusion

Basic research studies of drug-resistant epilepsy using semiautomated structural MRI techniques provide noninvasive, efficient and complementary approaches for identifying and quantifying different types of anatomical abnormalities. In cases of MTLE, evidence from structural MRI studies suggest that the brain areas responsible for generating spontaneous seizures include the sclerotic hippocampus, as well as other limbic brain areas. How well spatial patterns of GM loss derived from group-averaged data correspond to an individual subject and whether all GM abnormalities are epileptogenic is not known and complicates the translation of these techniques to clinical epilepsy. We anticipate that MRI and other brain imaging studies, such as DTI and functional MRI using connectivity analysis combined with electrophysiological recordings and clinical history of epilepsy and comorbidity, will provide new information on the spatial distribution of the epileptogenic zone corresponding to different types of MTLE with HS.

Future perspective

Improvements in image acquisition & processing

Improvements in image acquisition technology, primarily improving the spatial resolution and signal-to-noise ratio of images, will be a major area of future work. Enhancing the accuracy and computational speed of quantitative neuroimaging algorithms is also important. With an increasing number of studies using automated SBM, there is a clear need for more accurate segmentation algorithms, particularly for segmenting subcortical regions. Invention of advanced methods for efficient registration of multiple structural imaging modalities, as well as integrating structural and functional imaging data, are crucial for the development of multimodal imaging studies of epilepsy. These improvements will benefit epilepsy research and possibly preoperative and prognosis assessments, as well as image-guided depth electrode implantation and resection surgery.

Multimodal imaging studies of MTLE

Multimodal imaging of MTLE through integration of structural and functional imaging data is another important area of future work. Functional imaging modalities such as PET [134137] and SPECT [138] can be useful in determining the prospective functional deficit zone [92]. Functional MRI [139142] and magnetic source imaging or magnetoencephalography [143,144] can yield more information about the ictal onset zone and the spatial extent of interictal spiking, often termed the ‘irritative zone’ [8]. The combination of functional imaging modalities with structural neuroimaging and techniques that spatially coregister functional and structural signals (e.g., [81,82]) may elucidate the functional role of damaged hippocampal and extra-hippocampal brain areas in the progression of the disease and the surgical outcome of patients.

Constructing a network of the epileptogenic zone

The structural and functional connectivity corresponding to the epileptogenic network is not known. More advanced MRI studies, along with graph theory-based models, are required to shed light on the anatomical and functional connectivity of epileptogenic structures and reveal potential seizure propagation pathways. With studies finding extra-hippocampal neocortical atrophy, further imaging and electrophysiological investigations are crucial to verifying whether those areas reflect damage due to propagated seizure activity, or if they are morphological changes that independently support the generation of spontaneous seizures.

Emerging field of connectomics

The emerging science of connectomics provides an approach to understanding large-scale brain networks in which normal and abnormal brain functions operate [145,146]. Progress in this field might elucidate the yet unknown mechanisms of seizure onset and termination. Furthermore, the combination of connectomics with biologically inspired dynamic models of epilepsy provides a suitable platform to explore the dynamic network nature of epilepsy. This can eventually improve epilepsy surgery by identifying specific connections or network nodes that, if cut or resected, will abolish seizures.

Acknowledgements

The authors would like to thank A Reid for her valuable comments on this manuscript.

image

Financial & competing interests disclosure

This work was supported by the Natural Sciences and Engineering Research Council of Canada (N Memarian), and Grant NS071048 from the NIH (R Staba). The authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript apart from those disclosed.

No writing assistance was utilized in the production of this manuscript.

References

Papers of special note have been highlighted as:
* of interest
* of considerable interest

Website