- 1ENyS (Estudios en Neurociencias y Sistemas Complejos), Consejo Nacional de Investigaciones Científicas y Técnicas, Universidad Nacional Arturo Jauretche y Hospital El Cruce, Florencio Varela, Argentina
- 2Hospital de Pediatría J.P Garrahan, Departamento de Neuroimágenes, Buenos Aires, Argentina
- 3Centro Internacional Franco Argentino de Ciencias de la Información y de Sistemas (CIFASIS) - Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Grupo de Procesamiento de Señales Multimedia - División Neuroimágenes, Universidad Nacional de Rosario, Rosario, Argentina
- 4Instituto de Fisiología y Biofísica B. Houssay (IFIBIO), Consejo Nacional de Investigaciones Científicas y Técnicas, Departamento de Fisiología y Biofísica, Universidad de Buenos Aires, Buenos Aires, Argentina
- 5Hospital J.M Ramos Mejía, Centro de Epilepsia, Buenos Aires, Argentina
- 6Instituto de investigación en Electrónica, Control y Procesamiento de Señales (LEICI), Universidad Nacional de La Plata-Consejo Nacional de Investigaciones Científicas y Técnicas, La Plata, Argentina
- 7Instituto de Ingeniería y Agronomía, Universidad Nacional Arturo Jauretche, Florencio Varela, Argentina
Introduction: Several methods offer free volumetry services for MR data that adequately quantify volume differences in the hippocampus and its subregions. These methods are frequently used to assist in clinical diagnosis of suspected hippocampal sclerosis in temporal lobe epilepsy. A strong association between severity of histopathological anomalies and hippocampal volumes was reported using MR volumetry with a higher diagnostic yield than visual examination alone. Interpretation of volumetry results is challenging due to inherent methodological differences and to the reported variability of hippocampal volume. Furthermore, normal morphometric differences are recognized in diverse populations that may need consideration. To address this concern, we highlighted procedural discrepancies including atlas definition and computation of total intracranial volume that may impact volumetry results. We aimed to quantify diagnostic performance and to propose reference values for hippocampal volume from two well-established techniques: FreeSurfer v.06 and volBrain-HIPS.
Methods: Volumetry measures were calculated using clinical T1 MRI from a local population of 61 healthy controls and 57 epilepsy patients with confirmed unilateral hippocampal sclerosis. We further validated the results by a state-of-the-art machine learning classification algorithm (Random Forest) computing accuracy and feature relevance to distinguish between patients and controls. This validation process was performed using the FreeSurfer dataset alone, considering morphometric values not only from the hippocampus but also from additional non-hippocampal brain regions that could be potentially relevant for group classification. Mean reference values and 95% confidence intervals were calculated for left and right hippocampi along with hippocampal asymmetry degree to test diagnostic accuracy.
Results: Both methods showed excellent classification performance (AUC:> 0.914) with noticeable differences in absolute (cm3) and normalized volumes. Hippocampal asymmetry was the most accurate discriminator from all estimates (AUC:1~0.97). Similar results were achieved in the validation test with an automatic classifier (AUC:>0.960), disclosing hippocampal structures as the most relevant features for group differentiation among other brain regions.
Conclusion: We calculated reference volumetry values from two commonly used methods to accurately identify patients with temporal epilepsy and hippocampal sclerosis. Validation with an automatic classifier confirmed the principal role of the hippocampus and its subregions for diagnosis.
Introduction
Quantification of brain anatomical structures from magnetic resonance images (MR) is being increasingly used to recognize pathologic conditions such as temporal lobe epilepsy. Volumetric estimates of hippocampal size are postulated to be more sensitive than visual assessment alone, and also to improve clinical diagnosis in dementia and epilepsy (1–4).
Temporal lobe epilepsy with hippocampal sclerosis (HS) is one of the most frequent focal epilepsies in adults often refractory to pharmacological treatment; surgical resection is an effective therapeutic option for these patients achieving a seizure-free rate close to 80%.
Patients with temporal epilepsy and HS usually share clinical key features associated with the majority of seizure discharges including characteristic aura, arrest, alteration of consciousness (and amnesia), and automatisms. Relatively typical scalp EEG findings can be seen in the interictal state, at the seizure onset, during the course of the seizure, and postictally.
Hippocampal sclerosis is suspected in epilepsy patients when compatible ictal semiology and scalp EEG findings are found, but definitive diagnosis is established based on characteristic brain MR anomalies. Neuroimaging abnormalities are typically recognized in the hippocampus proper, including atrophy, loss of internal structure, and decreased T1- and increased T2-FLAIR signal intensity in clinical practice (5). Inspection of hippocampal coronal sections allows for a side-by-side comparison of asymmetry in volume, shape, and signal important for clinical diagnosis. Atrophy seems to be the most specific and signal changes the most sensitive biomarker in HS (6). Magnets with high field strengths above 3 T are able to depict subtle blurring of the internal architecture of the hippocampus on T2-weighted images (5). Originally, manual segmentation of the hippocampus based on anatomical knowledge and specific MRI landmarks was used to estimate structural volumes. Previous studies using these methods adequately identified lateralization of seizure origin in the temporal lobe of patients with HS. Earlier reports also documented a strong association between severity of histopathological anomalies and hippocampal volumes with an increased diagnostic yield of MR studies (7, 8)
The recent development of automatic volumetry methods such as FreeSurfer (FS) suite (9) and VolBrain (vB) HIPS (10), among others, makes it possible to account for hippocampal volume differences that may escape visual detection. Several studies validated the utility of hippocampal volumetry for HS detection in temporal epilepsy, mostly based on postoperative correlation or using ex vivo neuroimaging analysis (7–10). The potential of volumetry measures for postsurgical outcome prediction is still modest, with some improvement in reports considering subfields patterns of atrophy (11).
Since numerous publications demonstrate considerable differences in normative brain structural volumes across populations with different genetic backgrounds (12–14), volumetric estimates of the hippocampus in different populations are of particular concern. Previous reports consistently show hippocampal volume differences even when using analogous procedures (15–32). In this regard, normal anatomic variations and differences associated with the implemented methodology need to be considered for the interpretations of clinical conditions. An additional concern is that normative structural data from Latin America populations remains underrepresented.
The main objective of this work is to estimate reference values of sensitivity, specificity, and confidence intervals for classification of a local population of epilepsy patients with unilateral hippocampal sclerosis using two different volumetry approaches.
We analyzed T1 brain MRI volumetry of the hippocampus and hippocampal subfields in a cohort of 61 healthy subjects and in 57 epilepsy patients with confirmed unilateral mesial temporal sclerosis. Anatomical volumes were computed using two well-established automatic methods FS and vB. Recorded values for the hippocampus and subregions are expressed as absolute values (in cm3) and further normalized to brain size, quantified as a percent of total intracranial volume (TIV).
Furthermore, we provide hippocampal and subfield volume distribution for a community-based sample of healthy controls (HC) and evaluate subregion asymmetry differences in HC and between patients. We also compared the degree of asymmetry in left and right HS to investigate its relevance for diagnosis and the presence of distinctive patterns of atrophy at the subregion level.
Finally, a validation process was implemented to explore the contribution of non-hippocampal structures for group classification. This was performed using machine learning techniques, considering only FreeSurfer's morphometric information of whole-brain regions, including anatomical volumes and cortical thickness. Specifically, we used a feature selection technique to obtain the optimal number of features to discriminate between patients and HC, and then we performed three binary classifications for each group using a Monte Carlo cross-validation (MCCV) scheme (33) with a random forest classifier.
Materials and Methods
Participants
Patients were retrospectively enrolled based on medical records from the epilepsy unit between 2014 and 2019 at Nestor Kirchner—El Cruce Hospital at Florencio Varela, Buenos Aires, with a final diagnosis of temporal lobe epilepsy associated with unilateral right (n = 22, 15 females) and left (n = 35, 17 females) hippocampal sclerosis. Diagnosis was established using standardized practices as described in Oddo et al. (34) through clinical examination, assessment of disease history, semiology of seizures along with neuropsychological tests including prolonged video EEG, and compatible findings on 3-T MRI as suggested by ILAE (5). Thirty-one patients (54%) underwent surgical treatment with histopathology confirmation of HS after standard amygdalohippocampectomy with partial temporal lobectomy. The remaining patients are not yet operated but scheduled for surgery. Age- and sex-paired HC (n = 61, 44 females) were recruited mostly from local universities including students and academic personnel.
All participants gave written consent to participate and to make use of medical information for this study. The work described in this paper was carried out in accordance with the code of ethics of the world medical association (Declaration of Helsinki). Research ethics approval was obtained from the Hospital Research Ethics Board at El Cruce Hospital.
Imaging Characteristics and Analysis Methods
Only volumetric T1-weighted images were used in this study. These images were obtained as part of the clinical protocol for epilepsy workout in our institution and were acquired using the same MR unit (Philips Achieva 3T, 8-channel head coil), as recommended on recent specialized guidelines (5). Structural images consist of a 3D T1WI (FFE) sequence, with 180 slices of 1-mm isotropic resolution, TE= 3.3 msec, TR= 2300 msec, TI= 900 msec, flip angle= 9°, and field of view (FOV)= 240 × 240 × 180. Images were exported from the scanner and transformed to Nifti format for further analysis. For the statistical analysis, the same T1 volumetric images were processed using two established and freely available methods used to calculate brain region segmentation and quantification, namely, FreeSurfer Suite v6.0 (FS) working in an offline workstation and VolBrain-Hips 2016 (vB) that provides online services running on remote servers through a website interface.
Both methods offer validated hippocampal and hippocampal subfield segmentation through different approaches, distinct reference atlases, dissimilar processing times, and specific subfield region delineations. Output files and results from both methods were independently reviewed by two experienced neuroradiologists (JPP and GDS) looking for labeling inconsistencies and to assure quality control (no manual correction was performed). (See segmentation details for each method in Figure 1). Full documentation is available for processing details on each software platform, but here we describe a resumed version of each method.
Figure 1. (A) Examples of common subfield's atlas definition using vB and FS. Boxplots represent mean volumes as percent of TIV and whiskers the 95% confidence interval for the HC group. (B) Right hippocampal 3D models for the same subject, constructed using all subfields from both methods in 3D Slicer2; boxplots represent mean hippocampal volumes for left and right hippocampi in HC expressed in mm3 and in percent of TIV. Upper models show anterior–superior view, and lower images represent inferior projections for comparison. Shaded gray-wireframe area embodies whole hippocampus representation created from standard FS segmentation; note reduced size of the vB model. Most noticeable subregion differences are related to the definition of the anterior and posterior extent of CA1 and posterior subiculum; more medially and dorsally extended in vB. Coincidentally, the hippocampal tail, pre-subiculum, and para-subiculum regions defined in FS represents at least partially overlapping areas between methods. Other deep internal hippocampal structures such as the molecular layer, GCMLDG, fissure, and fimbria are individually ascribed only in FS (C). CA4-DG and CA2-CA3 subfields are jointly segmented in vB; CA4 and GCMLDF are grouped together in FS for comparison purposes. Volume differences are probably not only related to atlas definition; both approaches also show methodological discrepancies for intracranial volume computation. *Significant after Bonferroni correction. **Significant uncorrected p < 0.05. Paired-sample T-test for inter-hemispheric comparison of volumes as percent of TIV in HC.
FreeSurfer v6.0
All T1 brain volumes were processed to obtain a complete morphometric description. Cortical reconstruction and volumetric segmentation were performed in each participant's native space on FreeSurfer's1 (v 6.0) image analysis suite.
Briefly, image processing included removal of non-brain tissue using a hybrid watershed/surface deformation procedure, an automatic Talairach transformation, segmentation of the subcortical WM and deep GM volumetric structures (including hippocampus, amygdala, caudate, putamen, and ventricles), intensity normalization, tessellation of the GM–WM boundary, an automatic topology correction, and surface deformation following intensity gradients to optimally place the GM/WM and GM/CSF borders at the location where the greatest shift in intensity defines the transition to the other tissue class (9).
Once the cortical models were complete, a number of deformable procedures were performed for further data processing and analysis, including surface inflation and registration to a spherical atlas—based on individual cortical folding patterns to match cortical geometry across subjects, parcellation of the cerebral cortex into units relative to gyral and sulcal structure, and creation of a variety of surface-based data—including maps of curvature and sulcal depth. These methods use both intensity and continuity information of the entire 3D MR volume from segmentation and deformation procedures to produce representations of cortical thickness, which is calculated as the closest distance from the GM/WM boundary to the GM/CSF boundary at each vertex on the tessellated surface (9). The maps were created using spatial intensity gradients across tissue classes; therefore, they were not simply reliant on absolute signal intensity. Since the ensuing maps were not restricted to the voxel resolution of the original data, they can detect submillimeter differences between groups. Procedures for the measurement of cortical thickness have been validated against histological analysis and manual measurements. FreeSurfer morphometric procedures including principal hippocampal subfields have been demonstrated to show good test–retest reliability across scanner manufacturers and across field strengths (35, 36).
The FreeSurfer v6.0 algorithm follows a generative, parametric approach which focuses on modeling the spatial distribution of the hippocampal subregions and surrounding brain structures, which is learned from labeled training data. FreeSurfer v6.0 is built with a novel atlasing algorithm and ex vivo MRI data from autopsy brains. The segmentation provides 15 different subregions (12 used in for this work), based on the histology and morphometry from Rosene and Van Hoesen (37) and partly also on (38–41). See Figure 1 for details on implemented atlas and segmentation.
The ex vivo imaging protocol yields images with high resolution and signal-to-noise ratio. The segmentation algorithm is similar to Van Leemput (42) which is appropriate for analyzing in vivo MRI scans of different manufacturers using different T1 contrasts.
Compared to other new methods available, FreeSurfer involves a prolonged processing time 8~24 h running on standard single-core systems but also yielding extended quantification of additional brain structures including whole-brain regions beyond hippocampal formations.
We transformed the fixed-width-column plain-text files in which were written down the FreeSurfer outputs to comma separate values (csv) plain-text files which are more suitable to be opened as a Pandas' Dataframe (Python package). To ensure that classifiers did not consider features lacking specific regional information, we eliminated general features like cortical volume, mean cortical thickness, brain volume, and ventricle volume. Finally, to avoid potential biases due to differences among the participants' head size (43), volume measures of each area were normalized as a percentage of the estimated total intracranial volume (eTIV), provided also in FreeSurfer's results.
VolBrain—HIPS
VolBrain is a patch-based segmentation method for high-resolution hippocampus subfields. It has been validated and uses two publicly available segmentation protocols different from FreeSurfer on manually ex vivo segmented datasets (44, 45).
Both hippocampal segmentation protocols are available in volBrain-HIPS; Winterburn atlas disclosing 5 subregions was used for this work because it is more similar to the FreeSurfer v6.0 definition than Kulaga-Yoskovitz. VolBrain-HIPS is based on the combination of MOPAL (46), a multi-contrast extension of the OPAL (47) patch-based label fusion segmentation method and a novel neural network-based error corrector. The method uses an adaptation of MOPAL, a patch matching segmentation method to produce fast and accurate T1 brain segmentations. The method also works on standard MRI acquisition with image resolution of clinical practice as well as on single T1w or single T2w images. The VolBrain approach performs well also on mono-contrast T1w and T2w images as well as under standard resolution images that are upsampled using the LASR (48, 49) super-resolution method. The HIPS method also includes an error corrector post-processing step based on the use of a boosted ensemble of a neural network algorithm that is proposed to minimize systematic segmentation errors at post-processing. It works in a fully automated manner providing accurate results outperforming state-of-the-art methods such as MAGeT (50), ASHS (51), and SurfPatch (52) which usually require extended periods of computing time. VolBrain-HIPS takes <20 min and performs fast segmentation as well as subject-specific library registration that only requires estimating one non-linear registration over small regions to translate the whole library to the case to be segmented.
Finally, an online report is generated and results are plotted as absolute or percent values adjusted for intracranial volume against a normal reference standard for each anatomical region. Segmentation images can also be downloaded for evaluation purposes.
The same T1 volumetric images used for FreeSurfer v6.0 were uploaded to VolBrain-HIPS3 for this analysis, using Winterburn atlas definition for controls and patients (45). The produced final report including absolute values (mm3) and normalized to percent of brain volume were recorded for analysis.
Supervised Classification With a Random Forest Algorithm
As suggested by several previous publications (53–59), the quantification of non-hippocampal volumes in HS patients usually shows widespread modifications, involving the thalamus, amygdala, subcortical temporal white matter, temporal pole, and entorhinal cortex among others.
To study structural changes in the brain without any bias, we used FreeSurfer v6.0 metrics, specifically parcels of cortical thickness and volumes of all the cerebral structures in combination with machine learning methods based on Random Forest Classifiers (RFC) (60). This process was based on the implementation of an automatic classification algorithm to evaluate group discrimination performance considering morphometric contribution of whole-brain structures as independent features, without any a priori consideration. The selection of RFC was based by several premises: (i) We were interested in considering linear and more importantly non-linear relationships between all the features. (ii) As the number of samples was relatively low (although it is high for this type of studies), the parameter tuning should be an optional step. (iii) The interpretability of the relevant features in the classification should be clear. Given these conditions and the experience of the research team, we selected RFC as the best suitable algorithm for the analysis (61–63).
Preprocessed features of cortical–subcortical volumes and cortical thickness normalized to estimated total intracranial volume (eTIV) were analyzed via a progressive feature elimination (PFE) procedure (64) with a Monte Carlo cross-validation scheme (33). Briefly, we performed 30 shuffle-splits in which we randomly selected 80% of the samples (with balanced classes) to train the RFC and the remaining 20% for testing to optimize the accuracy of RFC by varying the number of features from all to a single one according to its classificatory relevance. RFC quantifies a feature's importance depending on how much the average Gini impurity index decreases in the forest due to its use as a node in a tree (65). We used this score to progressively eliminate features by removing the feature with the lowest importance at each iteration. Finally, we kept the N first features in the ranking, where N is the optimal number of features such that using more than N features fails to improve the classifier's performance.
The optimal number of features was selected visually by indicating the minimal quantity at which accuracy became constant. We used this fixed number of features to compute the accuracy, the confusion matrix, and the ROC curve, and to obtain each subject's probability of being in each group (HC, left HS, and right HS).
We implemented this processing framework to perform three classifications: (i) a binary classification to discriminate HC and HS; (ii) a binary classification to discriminate left HS and right HS; and (iii) a multiclass classification to discriminate HC, left HS, and right HS. For each classification, we obtained the optimal number of features, the list of defined features, and the classification performance metrics (accuracy, confusion matrix, and ROC curves). Asymmetry metrics were not included in these analyses given the conceptual basis that RFCs consider the relationship between features, and therefore the asymmetry between hemispheres regions was indirectly taken into account.
These analyses were performed with the RFC implemented in the Python's scikit-learn package, with a fixed number of trees (2000) and the recommended number of features (P) in each split, where P is the square root of the full set of features. The maximum depth in each tree was not restricted a priori, i.e., nodes were expanded until all leaves were pure or until all leaves contained less than two samples.
Statistical Analysis
Results were analyzed independently for each method using the Statistical Package for Social Sciences SPSS (Version 23; IBM, Armonk, New York). Volume mean average and 95% confidence intervals (CI) were calculated for each hemisphere. Receiver operating characteristic analyses were used to obtain optimal sensitivity, specificity, and 95% CI computed for left and right hippocampal sclerosis patients. A normalization of absolute values related to the total TIV was implemented and used for group comparison and correlation tests, since it was previously described as the most significant covariate to be considered (25). Normalization was performed by the following expression:
Both methods implement different atlas definitions and strategies to quantify TIV, thus precluding a direct comparison between absolute values.
Asymmetry degree was analyzed as an independent measure representing the difference between right and left regions divided by their mean (in percent) as implemented in vB and used in previous reports (25). Thus, positive values represent greater volumes on the right side.
Nominal variables were compared using the Chi square test. Paired-sample t-test (right vs. left) and ANCOVA (between groups) were used for normally distributed scalar variables adjusted for age and sex. Correlations were tested using the two-tailed Pearson coefficient controlling for age and sex. Significance level was adjusted for the effect of multiple comparisons using Bonferroni correction when appropriate. To test the difference between HS sides in the group of patients, an ANCOVA test was calculated on z-scores computed for each region using the following formula:
Age, sex, and clinical characteristics of epilepsy were included in the analysis as covariates.
Results
After correction for TIV, no significant correlation was found between age and sex with hippocampal or subfield volumes (p > 0.05) in controls or patients.
Controls and patients were paired according to age and sex, with female prevalence (controls 44f/17m, right HS 15f/7m, and left HS 17f/18m) not reaching significant differences (p.062). Groups were not different in relation to participants' age (p. 495), control subjects with a mean of 32 (18–62y), right HS patients group with 33 (21–64y), and left HS with 34 (19–52y).
Clinical characteristics of epilepsy including seizure frequency, age at onset, and epilepsy evolution time were similar (p > 0.05) in both groups of patients. Right HS patients had 7 (1–30) seizures per month, disease onset at 10 (1–40y), with a duration of 23 (6–40y), and left HS epilepsy patients presented 12 (1–90), 11 (1–32y), and 22 (2–49y), respectively.
No correlation was found between clinical features of epilepsy and hippocampus or subregion volumes.
Hippocampal Results
Estimated hippocampal volume and 95% confidence interval (CI) for controls on the right side were 3,454 cm3 (3.355–3.554)/0.2239% (0.2196–0.2283) for FS, and 2,480 cm3 (2.326–2.490)/0.1750% (0.1713–0.1787) for vB. Results for the left hippocampus were 3,398 cm3 (3.300–3.496)/0.2230% (0.2158–0.2248) for FS, and 2,320 cm3 (2.246–2.394)/0.1686% (0.1653–0.1720) for vB. Volume asymmetry was 1.6% (0.5–2.7) for FS and 3.6% (2.2–5) for vB with significant rightward lateralization (p < 0.003).
Mean ipsilateral hippocampal volume and 95% CI for right HS patients were 2,578 cm3 (2.401–2.755)/0.1743% (0.1603–0.1882) for FS and 1,429 cm3 (1.295–1.563)/0.1073% (0.0979–0.1167) for vB. Left hippocampus volume and 95% CI for the left HS patients were 2,560 cm3 (2.425–2.696)/0.1693% (0.1604–0.1783) for FS and 1,437 cm3 (1.324–1.549)/0.1055% (0.0981–0.1129) for vB.
Hippocampal asymmetry in the right HS group was −27.4% (−31.4/−23.5) for FS and −47% (−53.2/−41.6) for vB. Asymmetry in left HS patients was 33% (28.2/37.7) on FS and 53.2% (48–58.3) for vB. Hippocampal volumes ipsilateral to the HS side were significantly reduced compared with controls and also with the non-lesional side of right and left HS groups (p.000). Additionally, the right hippocampus was greater in left HS patients than in HC (FS, p.022) (see details in Tables 1 and 2).
Hippocampal asymmetry was the most reliable indicator for accurate classification between HC and right and left HS with an AUC:1 for Vb (measured in cm3 and in brain percent), an AUC:0.998 (using cm3), and an AUC:0.977 (in brain percent) based on FS. Optimal sensitivity–specificity was also calculated using hippocampal volumes with elevated accuracy (AUC:0.914 ~ 0.993) for patient classification. Detailed results are specified in Figure 2 and Table 3.
Figure 2. ROC curves: sensitivity and specificity for classification between HC; right (A) and left (B) HS patients. Prediction is calculated for hippocampal volumes in percent of TIV; expressed in cm3 and also using asymmetry differences for both methods. Best accuracy was obtained using asymmetry values; hippocampal volume estimates from vB showed slightly better accuracy than FS.
Table 3. ROC curve: reference values for highest sensitivity and specificity detection for right (A) and left (B) HS; using hippocampal volumes and asymmetry values from both methods.
To specifically account for atrophy differences among HS sides, z-score volumes for each hippocampus were compared, and no significant differences were found (p.692, FS and p.768, vB).
Results for Hippocampal Subfields
The mean volume and 95% CI estimates of hippocampal subfields for HC and patients are detailed in Tables 1 and 2. In the HC group, a significant rightward asymmetry of hippocampal subfields was recognized for CA1, CA2–CA3, and CA4-DG (in vB) and for CA1, CA3, CA4, molecular layer, hippocampal fissure, and GC-ML of DG (in FS). Leftward lateralization was recognized for the subiculum (vB) and pre-subiculum (FS) subregions (see details in Tables 1 and 2).
All subregions on the ipsilateral side of HS patients showed significant volume reduction compared with HC using vB, and most subfields were also reduced considering FS except for the right (p.446) and left (p.140) HATA, right and left fissure (p.1), and right ipsilateral fimbria (p.849).
Most hippocampal subfields contralateral to the sclerotic side in left HS patients revealed greater volumes compared with HC, specifically right CA1(p.048 in vB), CA1(p.002 in FS) and CA3(p.016), CA4 (p.025), HATA (p.009), molecular layer, and GC-ML-DC (p.018) on FS. Only the left subiculum (p.035 vB) of right HS patients was reduced compared with HC.
The only subregion with a significant volume difference between sides of the affected hemisphere in patients was CA2–CA3 (p.024) for the group of right HS patients (observed in vB). Accordingly, ipsilateral to the sclerotic side, CA2–CA3 (vB) and CA3 (FS) subfields in left HS patients were less atrophic than any other cornus ammonia division.
The most atrophic subfield ipsilateral to the sclerotic side for FS were CA4, GCMLDG, and molecular layer, and SLSRSM measured in vB in both right and left HS patients. See details in Figure 3.
Figure 3. Differences in hippocampal subregion atrophy; comparison of Z-scores between HS sides. Mean Z-score volume comparison between left and right HS; obtained from vB (A) and from FS (B). *Significant for ANCOVA test between groups; Bonferroni corrected (p < 0.05) adjusted for age, sex, and epilepsy characteristics. Whiskers represent 95% confidence interval.
Validation With the Automatic (Random Forest) Classifier
Our supervised machine learning validation process disclosed anatomical regions that were restricted to hippocampal subregions as the most relevant features to discriminate between patients and HC. In other terms, non-hippocampal regions were not identified as relevant for the classification.
The classifier was able to discriminate between controls and patients with a high accuracy in the three main classifications we performed: the classification between HC and patients (validation set mean accuracy: 0.907, AUC:0.960), between left and right HS patients (validation set mean accuracy: 0.91, AUC: 0.963), and between the three groups (validation set mean accuracy: 0.857, AUC: 0.960). The most important features and their relevance in each of three classifications are listed in Figure 4.
Figure 4. Validation process; results from the RFC algorithm. (A) Random forest performance. The confusion matrices show the percentage of correct (colored) and incorrect predictions for each class. The value was accumulated over the Monte Carlo cross validation (200-folds). (B) Progressive feature elimination. Random forest mean accuracy over the Monte Carlo cross-validation (30-folds) as a function of the number of features used to train the model. The features are sorted from most to the least important. The dashed lines show the optimal number of features for each classifier. (C) Feature importance for the 20 most important FS metrics. The boxplots show the normalized random forest feature importance distribution over the Monte Carlo cross-validation (200-folds). The colored boxes are the features which were selected by the progressive feature elimination procedure. The feature importance value was normalized with respect to the trivial importance level 1/N, where N is the number of features—that means, at the trivial level all the features have the same importance. Whiskers represent 95% confidence interval; small rhombuses indicate outliers.
Discussion
In this work, we define reference volumetric values and confidence intervals for hippocampus and hippocampal subfields using two commonly available approaches in a small community-based sample of healthy adults from Buenos Aires, Argentina. This is a limited sample but an important contribution to the field due to the scarce research literature on brain morphometric variations available in Latin America (66–68).
Since population variability on brain morphometric estimates are being increasingly reported (14, 69, 70, 70–72), it is important to consider the possibility of innate differences for adequate interpretation of MRI volumetry.
Several methods provide quantification of brain structures by using MRI data, including freely available softwares and online processing services that usually report adjusted values considering intracranial total volume, age, and sex as covariates. Unfortunately, wide variability exists related to the employed methodology that impairs appropriate comparisons of results between different techniques. Results are usually matched against a mixture of publicly available database of normal subjects that may not entirely account for variation among populations. Thus, absence of local references for normal and pathologic hippocampus volumes may also be a challenge for non-neuroimaging experts.
In this work, we report volumes of hippocampal structures and subregions that are specific for two different methods, evaluating patients from Latin America. The proposed reference values are intended to clarify the results obtained using two different methodologies, which are based on unequal anatomical definitions, and therefore the resulting scores cannot be directly used for cross-comparisons (see details in Figure 1).
We calculated mean volumes, confidence intervals, and cutoff estimations to recognize a regional sample of patients with confirmed unilateral mesial sclerosis and temporal lobe epilepsy with high sensitivity and specificity. Hippocampal asymmetry degree was the most accurate measure for classification regardless of the volumetry method used, as previously reported by others (3, 23, 31).
Our results are coincident with previous reports supporting rightward asymmetry for whole hippocampal volume not only in HC but also present in other animal species (73).
Interestingly, as recently reported (74), some hippocampal subregion volumes in our study were leftward lateralized in HC including the subiculum and pre-subiculum, the former based on volBrain and the latter on FreeSurfer. This discrepancy probably represents similar findings observed in overlapping areas related to known differences in atlas definitions (75) (see Figure 1).
Contrary to previous findings (31), our results did not show any significant correlation between hippocampus volume and its subfields with clinical features of epilepsy.
Few studies had focused on assessing subregion atrophy differences between HS sides based on imaging data. We found specific volume reduction of CA2–CA3(vB) in right HS patients with partial preservation in left HS patients. Future investigation using adequate methodology and involving a greater number of participants may confirm our findings. A distinctive pattern of modifications can be expected from left and right HS which are not usually considered on histopathology research, probably supporting differences in functional abilities (76–80).
To our knowledge, only one published study directly addressed asymmetry differences between hippocampal subregions among left and right HS patients using FS v6.0 (81). The authors found reduced contralateral volumes to the side of HS for presubiculum, HATA, and TAIL subfields. Unfortunately, information about known constitutional asymmetries present in HC (74) that could influence the results as in our analysis is not usually considered.
Another recent study used an approach similar to ours (but based on manual segmentation) and found greater (rather than reduced) volume of left subiculum (contrary to our findings) in right HS participants (32). Additionally, the authors also showed significant reduction of ipsilateral CA1 subfield compared against any other subregion on the sclerotic side.
An interesting observation from our analysis is a trend to find larger volumes on mesial–temporal structures contralateral to the side of HS in patients compared with HC. Diverse hippocampal subfields and also the hippocampus (FS) in the right (non-lesional) hemisphere of left HS patients support this assumption showing significant greater volumes compared to the same regions in healthy controls (Table 2). We should stress that in clinical practice the interpretation of hippocampal volumetry alone may not adequately identify some confirmed cases (~10%) with compatible clinical and paraclinical findings of HS which may only show subtle signal intensity changes on T2/FLAIR images (5, 82). Furthermore, it is important to note that a small group (~20%) of confirmed temporal lobe epilepsy patients without abnormal MRI finding will be postoperatively classified as “Gliosis only” without hippocampal sclerosis based on histopathology (83), showing no evidence of neuronal loss nor hippocampal volume reduction.
Supplementary functional imaging examinations are useful for diagnosis in temporal lobe epilepsy with HS and unremarkable MRI findings that may preserve normal hippocampal volumes. Interictal FDG-PET (2-[18F]-fluoro-2-D-deoxyglucose positron emission tomography) is a relatively widely available neuroimaging modality with high sensitivity (~80%) to disclose abnormal cortex hypo-metabolism in temporal lobe epilepsy (84, 85). Importantly, about 20% of patients with confirmed hippocampal sclerosis and normal MRI will show temporal cortex anomalies with reduced 18-FDG uptake (86, 87).
Although great progress has been made in recent years for preoperative diagnosis of HS using non-invasive methods, a considerable group of patients (20~40%) will fail to achieve complete seizure free after surgery (88, 89) following appropriate medical practices in experienced epilepsy centers. A recognized limitation of our study is the absence of histopathology information about recent standardized ILAE classification for HS subtype (ILAE HS I-III) (90) that could allow us to correlate volumetry findings with specific subfield anomalies. Nevertheless, some controversies remain concerning the role of histopathologic classification for predicting clinical evolution in HS patients and also regarding the feasibility of MRI-histopathology correlations, limited by the amount of brain sample available for examination. Additional benefits of MRI volumetry include the ability to examine the entire length of the sclerotic hippocampus and its contralateral homologous and also to consider inherited asymmetries for comparison.
Another caveat of this study is its relatively small sample size and also the uncertainty of segmentation accuracy of automated methods, to quantify structures on atrophic hippocampus. Some studies suggest that manual tracing methods may provide more accurate volumetric measurement than automated segmentation, especially in cases of HS (91, 92). However, validation results from FreeSurfer v6.0 developers indicate that subfield volumes still carry useful information, even when T1 images usually display limited contrast on the internal subregion boundaries (75). Equivalent methodology was also successfully implemented in previous studies on cognitive function and epilepsy (3, 4, 22, 23, 91, 93) with satisfactory results.
Contrary to previous observations supporting a fundamental role for cortical mesial–temporal regions, our machine learning-based validation process using an automatic algorithm failed to identify non-hippocampal structures such as the thalamus, temporal pole, fornix, or mammillary bodies as relevant for group classification. It shall be stressed that the abovementioned structures and others known to be involved in HS patients could falsely not been recognized as important due to a superior performance of hippocampal and subregion metrics in a trade-off between accuracy and number of analyzed features. Moreover, non-hippocampal anomalies preferentially involve white matter tracts (94, 95) and are usually related to prolonged epilepsy duration or high seizure frequency not considered in our validation process.
In conclusion, hippocampal anatomical structures are the most relevant features to recognize HS patients as confirmed by an automatic classification based on RFC. The local reference values proposed for hippocampal volumes and subfields may prove a useful guide for diagnosis in adult patients with temporal lobe epilepsy and suspected HS particularly for non-specialized radiologists.
Providing normal hippocampal reference values are a significant contribution to future studies focusing on regional morphometric variations in Latin America.
Finally, our results are also important for the interpretation of studies reporting hippocampal subfield volumes based on different atlas, which may show noticeable differences even when the same anatomical labels are used (96–98).
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Research and Bio-ethics review board at Hospital El Cruce, Carlos N. Kirschner. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
JP: study design, writing, statistical analysis, and quality revision. SK: study design and manuscript edition. PS: patients selection, clinical data, and follow-up. AN: patients selection, clinical data, and demography curation and analysis. SC: statistical analysis, data pre-processing, and quality revision. AD: data processing, analysis, classification algorithm, manuscript writing, and edition. GP: classification algorithm and manuscript writing and edition. JV: data pre-processing, analysis, quality revision, and manuscript edition. MV-A: data pre-processing, analysis, quality revision, and manuscript edition. PD-K: study design, data processing, analysis, classification algorithm, and manuscript writing and edition. All authors contributed to the article and approved the submitted version.
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.
Acknowledgments
We would like to acknowledge the kind contributions of colleagues from the Epilepsy Division, MRI-EEG technicians, and Administrative Personnel and especially from Prof. Paula Gonzalez at ENyS—El Cruce Hospital.
Footnotes
References
1. Jack CR Jr, Sharbrough FW, Cascino GD, Hirschorn KA, O'Brien PC, Marsh WR. Magnetic resonance image-based hippocampal volumetry: correlation with outcome after temporal lobectomy. Ann Neurol. (1992) 31:138–46. doi: 10.1002/ana.410310204
2. Jack CR Jr, Petersen RC, O'Brien PC, Tangalos EG. MR-based hippocampal volumetry in the diagnosis of Alzheimer's disease. Neurology. (1992) 42:183–8. doi: 10.1212/WNL.42.1.183
3. Louis S, Morita-Sherman M, Jones S, Vegh D, Bingaman W, Blumcke I, et al. Hippocampal sclerosis detection with neuroquant compared with neuroradiologists. AJNR Am J Neuroradiol. (2020) 41:591–7. doi: 10.3174/ajnr.A6454
4. Vos S, Gavin P, Olivia G, Hugh G, Frederik B, Ferran P, et al. Hippocampal profiling: localized magnetic resonance imaging volumetry and T2 relaxometry for hippocampal sclerosis. Epilepsia. (2020) 61:297–309. doi: 10.1111/epi.16416
5. Bernasconi A, Cendes F, Theodore WH, Gill R, Koepp M, Hogan R, et al. Recommendations for the use of structural magnetic resonance imaging in the care of patients with epilepsy: a consensus report from the International League Against Epilepsy Neuroimaging Task Force. Epilepsia. (2019) 60:1054–68. doi: 10.1111/epi.15612
6. Lee DH, Gao FQ, Rogers JM, Gulka I, Mackenzie IR, Parrent AG, et al. MR in temporal lobe epilepsy: analysis with pathologic confirmation. AJNR Am J Neuroradiol. (1998) 19:19–27.
7. Cascino GD, Jack CR Jr, Parisi JE, Sharbrough FW, Hirschorn KA, Meyer FB, et al. Magnetic resonance imaging-based volume studies in temporal lobe epilepsy: pathological correlations. Ann Neurol. (1991) 30:31–6. doi: 10.1002/ana.410300107
8. Kim JH, Tien RD, Felsberg GJ, Osumi AK, Lee N, Friedman AH. Fast spin-echo MR in hippocampal sclerosis: correlation with pathology and surgery. AJNR Am J Neuroradiol. (1995) 16:627–36.
10. Romero JE, Coupé P, Manjón JV. HIPS: a new hippocampus subfield segmentation method. Neuroimage. (2017) 163:286–95. doi: 10.1016/j.neuroimage.2017.09.049
11. Bernhardt BC, Hong SJ, Bernasconi A, Bernasconi N. Magnetic resonance imaging pattern learning in temporal lobe epilepsy: classification and prognostics. Ann Neurol. (2015) 77:436–46. doi: 10.1002/ana.24341
12. Liang P, Shi L, Chen N, Luo Y, Wang X, Liu K, et al. Construction of brain atlases based on a multi-center MRI dataset of 2020 Chinese adults. Sci Rep. (2016) 5:18216. doi: 10.1038/srep18216
13. Isamah N, Faison W, Payne ME, MacFall J, Steffens D, Beyer J, et al. Variability in frontotemporal brain structure: the importance of recruitment of African Americans in neuroscience research. PLoS ONE. (2010) 5:e13642. doi: 10.1371/journal.pone.0013642
14. Pai PP, Mandal PK, Punjabi K, Shukla D, Gael A, Joon S et al. BRAHMA: population specific T1, T2, and FLAIR weighted brain templates and their impact in structural and functional imaging studies. Magn Reson Imaging. (2020) 70:5–21. doi: 10.1016/j.mri.2019.12.009
15. Mohandas AN, Bharath RD, Prathyusha PV, Gupta AK. Hippocampal volumetry: normative data in the Indian population. Ann Indian Acad Neurol. (2014) 17:267–71. doi: 10.4103/0972-2327.138482
16. Jalaluddin WM, Mat Jusoh N, Ali Basahai IA, Abdullah MS, Abdul Karim AH, Gazali AK. Normalised MRI volumetry of the hippocampus among normal malay children and adolescents. Malays J Med Sci. (2013) 20:31–8.
17. Embong MF, Yaacob R, Abdullah MS, Abdul Karim AH, Ghazali AK, Jalaluddin WM. MR volumetry of hippocampus in normal adult malay of age 50 years old and above. Malays J Med Sci. (2013) 20:25–31.
18. Gonçalves-Pereira PM, Oliveira E, Insausti R. Análisis volumétrico mesiotemporal: valores normativos del hipocampo, la amígdala y el córtex entorrinal en la población adulta portuguesa [Quantitative volumetric analysis of the hippocampus, amygdala and entorhinal cortex: normative database for the adult Portuguese population]. Rev Neurol. (2006) 42:713–22. doi: 10.33588/rn.4212.2005352
19. Gaillard WD, Bhatia S, Bookheimer SY, Fazilat S, Sato S, Theodore WH. FDG-PET and volumetric MRI in the evaluation of patients with partial epilepsy. Neurology. (1995) 45:123–6. doi: 10.1212/WNL.45.1.123
20. Szabo CA, Xiong J, Lancaster JL, Rainey L, Fox P. Amygdalar and hippocampal volumetry in control participants: differences regarding handedness. AJNR Am J Neuroradiol. (2001) 22:1342–5.
21. Li YJ, Ga SN, Huo Y, Li SY, Gao XG. Characteristics of hippocampal volumes in healthy Chinese from MRI. Neurol Res. (2007) 29:803–6. doi: 10.1179/016164107X223557
22. Sone D, Sato N, Maikusa N, Ota M, Sumida K, Yokoyama K, et al. Automated subfield volumetric analysis of hippocampus in temporal lobe epilepsy using high-resolution T2-weighed MR imaging. Neuroimage Clin. (2016) 12:57–64. doi: 10.1016/j.nicl.2016.06.008
23. Granados Sánchez AM, Orejuela Zapata JF. Diagnosis of mesial temporal sclerosis: sensitivity, specificity and predictive values of the quantitative analysis of magnetic resonance imaging. Neuroradiol J. (2018) 31:50–9. doi: 10.1177/1971400917731301
24. Nobis L, Manohar SG, Smith SM, Alfaro-Almagro F, Jenkinson M, Mackay C, et al. Hippocampal volume across age: nomograms derived from over 19,700 people in UK Biobank. Neuroimage Clin. (2019) 23:101904. doi: 10.1016/j.nicl.2019.101904
25. Santyr BG, Goubran M, Lau JC, Kwan B, Fateme S, Lee D, et al. Investigation of hippocampal substructures in focal temporal lobe epilepsy with and without hippocampal sclerosis at 7T. J Magn Reson Imaging. (2017) 45:1359–70. doi: 10.1002/jmri.25447
26. Pedraza O, Bowers D, Gilmore R. Asymmetry of the hippocampus and amygdala in MRI volumetric measurements of normal adults. J Int Neuropsychol Soc. (2004) 10:664–78. doi: 10.1017/S1355617704105080
27. Özdemir M, Soysal H, Eraslan Ö, Dilli A. Normative hippocampal volumetric measurements using magnetic resonance imaging. Turk J Med Sci. (2019) 49:1464–70. doi: 10.3906/sag-1903-233
28. Laakso MP, Juottonen K, Partanen K, Vainio P, Soininen H. MRI volumetry of the hippocampus: the effect of slice thickness on volume formation. Magn Reson Imaging. (1997) 15:263–5. doi: 10.1016/S0730-725X(96)00390-6
29. Bhatia S, Bookheimer SY, Gaillard WD, Theodore WH. Measurement of whole temporal lobe and hippocampus for MR volumetry: normative data. Neurology. (1993) 43:2006–10. doi: 10.1212/WNL.43.10.2006
30. Farid N, Girard HM, Kemmotsu N, Smith M, Magda S, Lin W, et al. Temporal lobe epilepsy: quantitative MR volumetry in detection of hippocampal atrophy. Radiology. (2012) 264:542–50. doi: 10.1148/radiol.12112638
31. Kim JB, Suh SI, Kim JH. Volumetric and shape analysis of hippocampal subfields in unilateral mesial temporal lobe epilepsy with hippocampal atrophy. Epilepsy Res. (2015) 117:74–81. doi: 10.1016/j.eplepsyres.2015.09.004
32. Ji C, Zhu L, Chen C, Wang S, Zheng L, Li H. Volumetric changes in hippocampal subregions and memory performance in mesial temporal lobe epilepsy with hippocampal sclerosis. Neurosci Bull. (2018) 34:389–96. doi: 10.1007/s12264-017-0186-2
33. Xu QS, Liang YZ. Monte Carlo cross validation. Chemom Intell Lab Syst. (2001) 56:1–11. doi: 10.1016/S0169-7439(00)00122-2
34. Oddo S, Solis P, Consalvo D, Seaone E, Giagante B, Kochen S, et al. Postoperative neuropsychological outcome in patients with mesial temporal lobe epilepsy in Argentina. Epilepsy Res Treat. (2012) 2012:370351. doi: 10.1155/2012/370351
35. Han X, Jovicich J, Salat D, Avd Kouwe, Quinn BT, Czanner S, et al. Reliability of MRI-derived measurements of human cerebral cortical thickness: the effects of field strength, scanner upgrade and manufacturer. Neuroimage. (2006) 32:180–94. doi: 10.1016/j.neuroimage.2006.02.051
36. Brown EM, Pierce ME, Clark DC, Fischl B, Iglesias J, Milberg W, et al. Test-retest reliability of FreeSurfer automated hippocampal subfield segmentation within and across scanners. Neuroimage. (2020) 210:116563. doi: 10.1016/j.neuroimage.2020.116563
37. Rosene D, Van Hoesen, G. The hippocampal formation of the primate brain: a review of some comparative aspects of cytoarchitecture and connections. In: Jones E, Peters A. (Eds.), Cerebral Cortex. New York, NY: Plenum Press (1987). p. 6, 345–455. doi: 10.1007/978-1-4615-6616-8_9
38. Lorente de No, R. Studies on the structure of the cerebral cortex. II. Continuation of the study of the ammonic system. J Psychol Neurol. (1934) 46:113–77.
39. Insausti R, Amaral, D. Hippocampal formation. Mai J, Paxinos G. (Eds.) The Human Nervous System. (2011). p. 896–942. doi: 10.1016/B978-0-12-374236-0.10024-0
40. Green RC, Mesulam, M. Acetylcholinesterase fiber staining in the human hippo-campus and parahippocampal gyrus. J Comp Neurol. (1988) 273:488–99. doi: 10.1002/cne.902730405
41. Duvernoy, H. The Human Hippocampus. An Atlas of Applied Anatomy. J.F. Bergmann Verlag. Munich. (1988)
42. Van Leemput, K. Encoding probabilistic brain atlases using bayesian inference. Med Imaging IEEE Trans. (2009) 28:822–37. doi: 10.1109/TMI.2008.2010434
43. Whitwell JL, Crum WR, Watt HC, Fox NC. Normalization of cerebral volumes by use of intracranial volume: implications for longitudinal quantitative MR imaging. AJNR Am J Neuroradiol. (2001) 22:1483–9.
44. Yushkevich P, Avants B, Pluta J, Das S, Minkoff D, Mechanic-Hamilton D, et al. A high-resolution computational atlas of the human hippocampus from postmortem magnetic resonance imaging at 9.4 T. NeuroImage. (2009) 44:385–98. doi: 10.1016/j.neuroimage.2008.08.042
45. Winterburn JL, Pruessner JC, Chavez S, Schira M, Lobaugh N, Voineskos A, et al. A novel in vivo atlas of human hippocampal subfields using high-resolution 3 T magnetic resonance imaging. Neuroimage. (2013) 74:254–65. doi: 10.1016/j.neuroimage.2013.02.003
46. Romero JE, Coupe P, Manjón, J.V. High Resolution Hippocampus Subfield Segmentation Using Multispectral Multi-atlas Patch-Based Label Fusion. Patch-MI, LNCS 9993. (2016) 117–24. doi: 10.1007/978-3-319-47118-1_15
47. Giraud R, Ta V.T, Papadakis N, Manjón J.V, Collins D.L, Coupé P. Alzheimer's Disease Neuroimaging Initiative. An optimized patchmatch for multi-scale and multi-feature label fusion. NeuroImage. (2016) 124:770–82. doi: 10.1016/j.neuroimage.2015.07.076
48. Manjón JV, Coupé P, Buades A, Fonov V, Collins DL, Robles, et al. Non-local MRI upsampling. Med Image Anal. (2010) 14:784–92.
49. Coupé P, Manjón JV, Chamberland M, Descoteaux M, Hiba B. Collaborative patch- based super-resolution for diffusion-weighted images. Neuroimage. (2013) 83:245–61. doi: 10.1016/j.neuroimage.2013.06.030
50. Pipitone J, Park MTM, Winterburn J, Lett TA, Lerch JP, Pruessner C, et al. The Alzheimer's disease neuroimaging initiative, 2014 multi- atlas segmentation of the whole hippocampus and subfields using multiple automatically generated templates. NeuroImage. (2014) 101:494–512. doi: 10.1016/j.neuroimage.2014.04.054
51. Yushkevich PA, Pluta JB, Wang H, Xie L, Ding SL, Gertje, et al. Automated volumetry and regional thickness analysis of hippocampal subfields and medial temporal cortical structures in mild cognitive impairment. Hum Brain Mapp. (2015) 36:258–87. doi: 10.1002/hbm.22627
52. Caldairou B, Bernhardt BC, Kulaga-Yoskovitz J, Kim H, Bernasconi N, Bernasconi, et al. A Surface Patch-Based Segmentation Method For Hippocampal Subfields. MICCAI, Part II, LNCS 9901. (2016) 379 - 387. doi: 10.1007/978-3-319-46723-8_44
53. Mahmoudi F, Elisevich K, Bagher-Ebadian H, Nazem-Zadeh M-R, Davoodi-Bojd E, Schwalb JM, et al. Data mining MR image features of select structures for lateralization of mesial temporal lobe epilepsy. PLoS ONE. (2018) 13:e0199137. doi: 10.1371/journal.pone.0199137
54. Mo J, Liu Z, Sun K, Ma Y, Hu W, Zhang C, et al. Automated detection of hippocampal sclerosis using clinically empirical and radiomics features. Epilepsia. (2019) 60:2519–29. doi: 10.1111/epi.16392
55. Rudie JD, Colby JB, Salamon N. Machine learning classification of mesial temporal sclerosis in epilepsy patients. Epilepsy Res. (2015) 117:63–9. doi: 10.1016/j.eplepsyres.2015.09.005
56. Cendes F, Andermann F, Gloor P, Evans A, Jones_Gotman M, Watson C, et al. MRI volumetric measurement of amygdala and hippocampus in temporal lobe epilepsy. Neurology. (1993) 43:719–25. doi: 10.1212/WNL.43.4.719
57. Bernhardt BC, Kim H, Bernasconi N. Patterns of subregional mesiotemporal disease progression in temporal lobe epilepsy. Neurology. (2013) 81:1840–7. doi: 10.1212/01.wnl.0000436069.20513.92
58. Garbelli R, Milesi G, Medici V, Villani F, Didoto G, Delea F, et al. Blurring in patients with temporal lobe epilepsy: clinical, high-field imaging and ultrastructural study. Brain. (2012) 135:2337–49. doi: 10.1093/brain/aws149
59. Kuzniecky R, Bilir E, Gilliam F, Faught E, Martin R, Hugg J. Quantitative MRI in temporal lobe epilepsy: evidence for fornix atrophy. Neurology. (1999) 53:496–501. doi: 10.1212/WNL.53.3.496
60. Donnelly-Kehoe PA, Pascariello GO, Gómez JC. Alzheimers Disease Neuroimaging Initiative. Looking for Alzheimer's Disease morphometric signatures using machine learning techniques. J Neurosci Methods. (2018) 302:24–34. doi: 10.1016/j.jneumeth.2017.11.013
61. Donnelly-Kehoe PA, Pascariello GO, García AM, Hodges J, Miller B, Rosen H, et al. Robust automated computational approach for classifying frontotemporal neurodegeneration: multimodal/multicenter neuroimaging. Alzheimers Dement. (2019) 11:588–98. doi: 10.1016/j.dadm.2019.06.002
62. Misra S, Wu Y. Machine learning assisted segmentation of scanning electron microscopy images of organic-rich shales with feature extraction and feature ranking. Machine Learn Subsurface Characterization. (2019) 289–314. doi: 10.1016/B978-0-12-817736-5.00010-7
63. Palczewska A, Palczewski J, Robinson RM, & Neagu D. Interpreting random forest models using a feature contribution method. In IEEE 14th International Conference on Information Reuse & Integration (IRI). San Francisco, CA: IEEE (2013). p. 112–9. doi: 10.1109/IRI.2013.6642461
65. Finkelsztejn A, Fragoso YD, Bastos EA, Duarte J, Santos Varela J, Houbrechts, et al. Intercontinental validation of brain volume measurements using MSmetrix. Neuroradiol J. (2018) 31:147–9. doi: 10.1177/1971400917721521
66. Curiati PK, Tamashiro JH, Squarzoni P, Duran F, Santos L, Wojngarten M, et al. Brain structural variability due to aging and gender in cognitively healthy Elders: results from the Sáo Paulo Ageing and Health study. AJNR Am J Neuroradiol. (2009) 30:1850–6. doi: 10.3174/ajnr.A1727
67. Squarzoni P, Tamashiro-Duran J, Souza Duran FL, Santos L, Vallada H, Menezes P, et al. Relationship between regional brain volumes and cognitive performance in the healthy aging: an MRI study using voxel-based morphometry. J Alzheimers Dis. (2012) 31:45–58. doi: 10.3233/JAD-2012-111124
68. Tang Y, Hojatkashani C, Dinov I, Sun B, Lingshong F, Lin XW, et al. The construction of a Chinese MRI brain atlas: a morphometric comparison study between Chinese and Caucasian cohorts. Neuroimage. (2010) 51:33–41. doi: 10.1016/j.neuroimage.2010.01.111
69. Rao NP, Jeelani H, Achalia R, Achiala G, Arpitha J, Bharat R, et al. Population differences in brain morphology: need for population specific brain template. Psychiatry Res Neuroimaging. (2017) 265:1–8. doi: 10.1016/j.pscychresns.2017.03.018
70. Yang G, Zhou S, Bozek J, Dong HM, Han M, Zuo XN, et al. Sample sizes and population differences in brain template construction. Neuroimage. (2020) 206:116318. doi: 10.1016/j.neuroimage.2019.116318
71. Sivaswamy J, Thottupattu AJ, Mehta R, Sheelakumari R, Kesavadas C. Construction of Indian human brain atlas. Neurol India. (2019) 67:229–34. doi: 10.4103/0028-3886.253639
72. Lee JS, Lee DS, Kim J, Kyeong K, Kong E, Kong H, et al. Development of Korean standard brain templates. J Korean Med Sci. (2005) 20:483–8. doi: 10.3346/jkms.2005.20.3.483
73. Freeman H, Cantalupo C and Hopkins W. Asymmetries in the hippocampus and amygdala of chimpanzees (Pan troglodytes). Behav Neurosci. (2004) 118:1460–5. doi: 10.1037/0735-7044.118.6.1460
74. Sarica A, Vasta R, Novellino F, Vaccaro M, Cerasa A, Quattrone A, et al. MRI asymmetry index of hippocampal subfields increases through the continuum from the mild cognitive impairment to the Alzheimer's disease. Front Neurosci. (2018) 12:576. doi: 10.3389/fnins.2018.00576
75. Iglesias JE, Augustinack JC, Nguyen K, Player CM, Player A, Wright M et al. A computational atlas of the hippocampal formation using ex vivo, ultra-high resolution MRI: application to adaptive segmentation of in vivo MRI. Neuroimage. (2015) 115:117–37. doi: 10.1016/j.neuroimage.2015.04.042
76. Coras R, Pauli E, Li J, Schwarz M, Rossler K, Buchfelder M, et al. Differential influence of hippocampal subfields to memory formation: insights from patients with temporal lobe epilepsy. Brain. (2014) 137:1945–57. doi: 10.1093/brain/awu100
77. Comper SM, Jardim AP, Corso JT, Botelho-Gaca L, Silva-Noffs M, Lancellotti C, et al. Impact of hippocampal subfield histopathology in episodic memory impairment in mesial temporal lobe epilepsy and hippocampal sclerosis. Epilepsy Behav. (2017) 75:183–9. doi: 10.1016/j.yebeh.2017.08.013
78. Daumas S, Halley H, Francés B, Lassalle JM. Encoding, consolidation, and retrieval of contextual memory: differential involvement of dorsal CA3 and CA1 hippocampal subregions. Learn Mem. (2005) 12:375–82. doi: 10.1101/lm.81905
79. Kesner RP, Lee I, Gilbert, P. A behavioral assessment of hippocampal function based on a subregional analysis. Rev Neurosci. (2004) 15:333–51. doi: 10.1515/REVNEURO.2004.15.5.333
80. Long J, Feng Y, Liao H, Zhou Q, Urbin MA. Motor sequence learning is associated with hippocampal subfield volume in humans with medial temporal lobe epilepsy. Front Hum Neurosci. (2018) 12:367. doi: 10.3389/fnhum.2018.00367
81. Kreilkamp BAK, Weber B, Elkommos SB, Richardson MP, Keller SS. Hippocampal subfield segmentation in temporal lobe epilepsy: relation to outcomes. Acta Neurol Scand. (2018) 137:598–608. doi: 10.1111/ane.12926
82. Briellmann RS, Kalnins RM, Berkovic SF, Jackson GD. Hippocampal pathology in refractory temporal lobe epilepsy: T2-weighted signal change reflects dentate gliosis. Neurology. (2002) 58:265–71. doi: 10.1212/WNL.58.2.265
83. Blumcke I, Pauli E, Clusmann H, Schramm J, Becker A, Elger C et al. A new clinico-pathological classification system for mesial temporal sclerosis. Acta Neuropathol. (2007):113:235–44. doi: 10.1007/s00401-006-0187-0
84. Salanova V, Markand O, Worth R, Smith R, Wellman H, Hutchins, et al. FDG-PET and MRI in temporal lobe epilepsy: relationship to febrile seizures, hippocampal sclerosis and outcome. Acta Neurol Scand. (1998) 97:146–53. doi: 10.1111/j.1600-0404.1998.tb00628.x
85. Guedj E, Bonini F, Gavaret M, Trebuchon A, Aubert S, Boucekine M, et al. 18FDG-PET in different subtypes of temporal lobe epilepsy: SEEG validation and predictive value. Epilepsia. (2015) 56:414–21. doi: 10.1111/epi.12917
86. Henry TR, Roman DD. Presurgical epilepsy localization with interictal cerebral dysfunction. Epilepsy Behav. (2011) 20:194–208. doi: 10.1016/j.yebeh.2010.12.008
87. Yang PF, Pei JS, Zhang HJ, Lin Q, Mei Z, Zhong Z-H, et al. Long-term epilepsy surgery out- comes in patients with PET-positive, MRI-negative temporal lobe epilepsy. Epilepsy Behav. (2014) 41:91–7. doi: 10.1016/j.yebeh.2014.09.054
88. LoPinto-Khoury C, Sperling MR, Skidmore C, Nei M, Evans J, Sharan A, et al. Surgical out- come in PET-positive, MRI-negative patients with temporal lobe epilepsy. Epilepsia. (2012) 53:342–48. doi: 10.1111/j.1528-1167.2011.03359.x
89. McIntosh AM, Wilson SJ, Berkovic SF. Seizure outcome after temporal lobectomy: current research practice and findings. Epilepsia. (2001) 42:1288–307. doi: 10.1046/j.1528-1157.2001.02001.x
90. Blümcke I, Thom M, Aronica E, Armstrong D, Bartolomei F, Bernasconi A, et al. International consensus classification of hippocampal sclerosis in temporal lobe epilepsy: a task force report from the ILAE Commission on Diagnostic Methods. Epilepsia. (2013) 54:1315–29. doi: 10.1111/epi.12220
91. Steve TA, Gargula J, Misaghi E, Nowacki T, Schmitt L, Wheatley B, et al. Hippocampal subfield measurement and ILAE hippocampal sclerosis subtype classification with in vivo 4.7 tesla MRI. Epilepsy Res. (2020) 161:106279. doi: 10.1016/j.eplepsyres.2020.106279
92. Pardoe HR, Pell GS, Abbott DF, Jackson, G.D. Hippocampal volume assessment in temporal lobe epilepsy: how good is automated segmentation? Epilepsia. (2009) 50:2586–92. doi: 10.1111/j.1528-1167.2009.02243.x
93. Winston GP, Cardoso MJ, Williams EJ, Brudett J, Bartlett P, Espak M, et al. Automated hippocampal segmentation in patients with epilepsy: available free online. Epilepsia. (2013) 54:2166–73. doi: 10.1111/epi.12408
94. Besson P, Dinkelacker V, Valabregue R, Thivard L, Leclerc X, Baulac M, et al. Structural connectivity differences in left and right temporal lobe epilepsy. Neuroimage. (2014) 100:135–44. doi: 10.1016/j.neuroimage.2014.04.071
95. Rodríguez-Cruces R, Concha L. White matter in temporal lobe epilepsy: clinico-pathological correlates of water diffusion abnormalities. Quant Imaging Med Surg. (2015) 5:264–78.
96. Zammit AR, Ezzati A, Zimmerman ME, Lipton RB, Lipton ML, Katz MJ. Roles of hippocampal subfields in verbal and visual episodic memory. Behav Brain Res. (2017) 317:157–62. doi: 10.1016/j.bbr.2016.09.038
97. Ono SE, de Carvalho Neto A, Joaquim MJM, Dos Santos GR, de Paola L, Silvado CES. Mesial temporal lobe epilepsy: revisiting the relation of hippocampal volumetry with memory deficits. Epilepsy Behav. (2019) 100:106516. doi: 10.1016/j.yebeh.2019.106516
98. Peixoto-Santos JE, Carvalho LED, Kandratavicius L, Diniz PRB, Scandiuzzi RC, Coras R et al. Manual hippocampal subfield segmentation using high-field MRI: impact of different subfields in hippocampal volume loss of temporal lobe epilepsy patients. Front. Neurol. (2018) 9:927. doi: 10.3389/fneur.2018.00927
Keywords: epilepsy, volumetry, hippocampal sclerosis, random forest classifier, MRI
Citation: Princich JP, Donnelly-Kehoe PA, Deleglise A, Vallejo-Azar MN, Pascariello GO, Seoane P, Veron Do Santos JG, Collavini S, Nasimbera AH and Kochen S (2021) Diagnostic Performance of MRI Volumetry in Epilepsy Patients With Hippocampal Sclerosis Supported Through a Random Forest Automatic Classification Algorithm. Front. Neurol. 12:613967. doi: 10.3389/fneur.2021.613967
Received: 05 October 2020; Accepted: 18 January 2021;
Published: 22 February 2021.
Edited by:
Hans-Peter Müller, University of Ulm, GermanyReviewed by:
Jordi A. Matias-Guiu, Hospital Clínico San Carlos, SpainAlessia Sarica, University of Magna Graecia, Italy
Copyright © 2021 Princich, Donnelly-Kehoe, Deleglise, Vallejo-Azar, Pascariello, Seoane, Veron Do Santos, Collavini, Nasimbera and Kochen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Juan Pablo Princich, jpprincich@gmail.com