Quantitation in MRI : application to ageing and epilepsy

Multi-atlas propagation and label fusion techniques have recently been developed for segmenting the human brain into multiple anatomical regions. In this thesis, I investigate possible adaptations of these current state-of-the-art methods. The aim is to study ageing on the one hand, and on the other hand temporal lobe epilepsy as an example for a neurological disease. Overall effects are a confounding factor in such anatomical analyses. Intracranial volume (ICV) is often preferred to normalize for global effects as it allows to normalize for estimated maximum brain size and is hence independent of global brain volume loss, as seen in ageing and disease. I describe systematic differences in ICV measures obtained at 1.5T versus 3T, and present an automated method of measuring intracranial volume, Reverse MNI Brain Masking (RBM), based on tissue probability maps in MNI standard space. I show that this is comparable to manual measurements and robust against field strength differences. Correct and robust segmentation of target brains which show gross abnormalities, such as ventriculomegaly, is important for the study of ageing and disease. We achieved this with incorporating tissue classification information into the image registration process. The best results in elderly subjects, patients with TLE and healthy controls were achieved using a new approach using multi-atlas propagation with enhanced registration (MAPER). I then applied MAPER to the problem of automatically distinguishing patients with TLE with (TLE-HA) and without (TLE-N) hippocampal atrophy on MRI from controls, and determine the side of seizure onset. MAPER-derived structural volumes were used for a classification step consisting of selecting a set of discriminatory structures and applying support vector machine on the structural volumes as well as morphological similarity information such as volume difference obtained with spectral analysis. Acccuracies were 91-100 %, indicating that the method might be clinically useful. Finally, I used the methods developed in the previous chapters to investigate brain regional volume changes across the human lifespan in over 500 healthy subjects between 20 to 90 years of age, using data from three different scanners (2x 1.5T, 1x 3T), using the IXI database. We were able to confirm several known changes, indicating the veracity of the method. In addition, we describe the first multi-region, whole-brain database of normal ageing

Automatic volumetry can reveal visually undetected disease features on brain MR images in temporal lobe epilepsy

Automatic volumetry can reveal visually undetected disease features on brain MR images in temporal lobe epilepsy

Advanced Magnetic Resonance Imaging and Quantitative Analysis Approaches in Patients with Refractory Focal Epilepsy

Background Epilepsy has a high prevalence of 1%, which makes it the most common serious neurological disorder. The most difficult to treat type of epilepsy is temporal lobe epilepsy (TLE) with its most commonly associated lesion being hippocampal sclerosis (HS). About 30-50% of all patients undergoing resective surgery of epileptogenic tissue continue to have seizures postoperatively. Indication for this type of surgery is only given when lesions are clearly visible on magnetic resonance images (MRI). About 30% of all patients with focal epilepsy do not show an underlying structural lesion upon qualitative neuroradiological MRI assessment (MRI-negative). Objectives The work presented in this thesis uses MRI data to quantitatively investigate structural differences between brains of patients with focal epilepsy and healthy controls using automated imaging preprocessing and analysis methods. Methods All patients studied in this thesis had electrophysiological evidence of focal epilepsy, and underwent routine clinical MRI prior to participation in this study. There were two datasets and both included a cohort of age-matched controls: (i) Patients with TLE and associated HS who later underwent selective amygdalahippocampectomy (cohort 1) and (ii) MRI-negative patients with medically refractory focal epilepsy (cohort 2). The participants received high- resolution routine clinical MRI as well as additional sequences for gray and white matter (GM/WM) structural imaging. A neuroradiologist reviewed all images prior to analysis. Hippocampal subfield volume and automated tractography analysis was performed in patients with TLE and HS and related to post-surgical outcomes, while images of MRI- negative patients were analyzed using voxel-based morphometry (VBM) and manual/automated tractography. All studies were designed to detect quantitative differences between patients and controls, except for the hippocampal subfield analysis as control data was not available and comparisons were limited to patients with persistent postoperative seizures and those without. Results 1. Automated hippocampal subfield analysis (cohort 1): The high-resolution hippocampal subfield segmentation technique cannot establish a link between hippocampal subfield volume loss and post-surgical outcome. Ipsilateral and contralateral hippocampal subfield volumes did not correlate with clinical variables such as duration of epilepsy and age of onset of epilepsy. 2. Automated WM diffusivity analysis (cohort 1): Along-the-tract analysis showed that ipsilateral tracts of patients with right/left TLE and HS were more extensively affected than contralateral tracts and the affected regions within tracts could be specified. The extent of hippocampal atrophy (HA) was not related to (i) the diffusion alterations of temporal lobe tracts or (ii) clinical characteristics of patients, whereas diffusion alterations of ipsilateral temporal lobe tracts were significantly related to age at onset of epilepsy, duration of epilepsy and epilepsy burden.Patients without any postoperative seizure symptoms (excellent outcomes) had more ipsilaterally distributed WM tract diffusion alterations than patients with persistent postoperative seizures (poorer outcomes), who were affected bilaterally. 3. Automated epileptogenic lesion detection (cohort 2): Comparison of individual patients against the controls revealed that focal cortical dysplasia (FCD) can be detected automatically using statistical thresholds. All sites of dysplasia reported at the start of the study were detected using this technique. Two additional sites in two different patients, which had previously escaped neuroradiological assessment, could be identified. When taking these statistical results into account during re-assessment of the dedicated epilepsy research MRI, the expert neuroradiologist was able to confirm these as lesions. 4. Manual and automated WM diffusion tensor imaging (DTI) analysis (cohort 2): The analysis of consistency across approaches revealed a moderate to good agreement between extracted tract shape, morphology and space and a strong correlation between diffusion values extracted with both methods. While whole-tract DTI-metrics determined using Automated Fiber Quantification (AFQ) revealed correlations with clinical variables such as age of onset and duration of epilepsy, these correlations were not found using the manual technique. The manual approach revealed more differences than AFQ in group comparisons of whole-tract DTI-metrics. Along-the-tract analysis provided within AFQ gave a more detailed description of localized diffusivity changes along tracts, which correlated with clinical variables such as age of onset and epilepsy duration. Conclusions While hippocampal subfield volume loss in patients with TLE and HS was not related with any clinical variables or to post-surgical outcomes, WM tract diffusion alterations were more bilaterally distributed in patients with persistent postoperative seizures, compared to patients with excellent outcomes. This may indicate that HS as an initial precipitating injury is not affected by clinical features of the disorder and automated hippocampal subfield mapping based on MRI is not sufficient to stratify patients according to outcome. Presence of persisting seizures may depend on other pathological processes such as seizure propagation through WM tracts and WM integrity. Automated and time-efficient three-dimensional voxel-based analysis may complement conventional visual assessments in patients with MRI-negative focal epilepsy and help to identify FCDs escaping routine neuroradiological assessment. Furthermore, automated along-the-tract analysis may identify widespread abnormal diffusivity and correlations between WM integrity loss and clinical variables in patients with MRI-negative epilepsy. However, automated WM tract analysis may differ from results obtained with manual methods and therefore caution should be exercised when using automated techniques

Imaging of epileptic activity using EEG-correlated functional MRI.

This thesis describes the method of EEG-correlated fMRI and its application to patients with epilepsy. First, an introduction on MRI and functional imaging methods in the field of epilepsy is provided. Then, the present and future role of EEG-correlated fMRI in the investigation of the epilepsies is discussed. The fourth chapter reviews the important practicalities of EEG-correlated fMRI that were addressed in this project. These included patient safety, EEG quality and MRI artifacts during EEG-correlated fMRI. Technical solutions to enable safe, good quality EEG recordings inside the MR scanner are presented, including optimisation of the EEG recording techniques and algorithms for the on-line subtraction of pulse and image artifact. In chapter five, a study applying spike-triggered fMRI to patients with focal epilepsy (n = 24) is presented. Using statistical parametric mapping (SPM), cortical Blood Oxygen Level-Dependent (BOLD) activations corresponding to the presumed generators of the interictal epileptiform discharges (IED) were identified in twelve patients. The results were reproducible in repeated experiments in eight patients. In the remaining patients no significant activation (n = 10) was present or the activation did not correspond to the presumed epileptic focus (n = 2). The clinical implications of this finding are discussed. In a second study it was demonstrated that in selected patients, individual (as opposed to averaged) IED could also be associated with hemodynamic changes detectable with fMRI. Chapter six gives examples of combination of EEG-correlated fMRI with other modalities to obtain complementary information on interictal epileptiform activity and epileptic foci. One study compared spike-triggered fMRI activation maps with EEG source analysis based on 64-channel scalp EEG recordings of interictal spikes using co-registration of both modalities. In all but one patient, source analysis solutions were anatomically concordant with the BOLD activation. Further, the combination of spike- triggered fMRI with diffusion tensor and chemical shift imaging is demonstrated in a patient with localisation-related epilepsy. In chapter seven, applications of EEG-correlated fMRI in different areas of neuroscience are discussed. Finally, the initial imaging findings with the novel technique for the simultaneous and continuous acquisition of fMRI and EEG data are presented as an outlook to future applications of EEG-correlated fMRI. In conclusion, the technical problems of both EEG-triggered fMRI and simultaneous EEG-correlated fMRI are now largely solved. The method has proved useful to provide new insights into the generation of epileptiform activity and other pathological and physiological brain activity. Currently, its utility in clinical epileptology remains unknown

Improving the clinico-radiological association in neurological diseases

Despite the key role of magnetic resonance imaging (MRI) in the diagnosis and monitoring of multiple sclerosis (MS) and cerebral small vessel disease (SVD), the association between clinical and radiological disease manifestations is often only moderate, limiting the use of MRI-derived markers in the clinical routine or as endpoints in clinical trials. In the projects conducted as part of this thesis, we addressed this clinico-radiological gap using two different approaches. Lesion-symptom association: In two voxel-based lesion-symptom mapping studies, we aimed at strengthening lesion-symptom associations by identifying strategic lesion locations. Lesion mapping was performed in two large cohorts: a dataset of 2348 relapsing-remitting MS patients, and a population-based cohort of 1017 elderly subjects. T2-weighted lesion masks were anatomically aligned and a voxel-based statistical approach to relate lesion location to different clinical rating scales was implemented. In the MS lesion mapping, significant associations between white matter (WM) lesion location and several clinical scores were found in periventricular areas. Such lesion clusters appear to be associated with impairment of different physical and cognitive abilities, probably because they affect commissural and long projection fibers. In the SVD lesion mapping, the same WM fibers and the caudate nucleus were identified to significantly relate to the subjects’ cerebrovascular risk profiles, while no other locations were found to be associated with cognitive impairment. Atrophy-symptom association: With the construction of an anatomical physical phantom, we aimed at addressing reliability and robustness of atrophy-symptom associations through the provision of a “ground truth” for atrophy quantification. The built phantom prototype is composed of agar gels doped with MRI and computed tomography (CT) contrast agents, which realistically mimic T1 relaxation times of WM and grey matter (GM) and showing distinguishable attenuation coefficients using CT. Moreover, due to the design of anatomically simulated molds, both WM and GM are characterized by shapes comparable to the human counterpart. In a proof-of-principle study, the designed phantom was used to validate automatic brain tissue quantification by two popular software tools, where “ground truth” volumes were derived from high-resolution CT scans. In general, results from the same software yielded reliable and robust results across scans, while results across software were highly variable reaching volume differences of up to 8%

Quantitative MRI correlates of hippocampal and neocortical pathology in intractable temporal lobe epilepsy

Intractable or drug-resistant epilepsy occurs in over 30% of epilepsy patients, with many of these patients undergoing surgical excision of the affected brain region to achieve seizure control. Advances in MRI have the potential to improve surgical treatment of epilepsy through improved identification and delineation of lesions. However, validation is currently needed to investigate histopathological correlates of these new imaging techniques. The purpose of this work is to investigate histopathological correlates of quantitative relaxometry and DTI from hippocampal and neocortical specimens of intractable TLE patients. To achieve this goal I developed and evaluated a pipeline for histology to in-vivo MRI image registration, which finds dense spatial correspondence between both modalities. This protocol was divided in two steps whereby sparsely sectioned histology from temporal lobe specimens was first registered to the intermediate ex-vivo MRI which is then registered to the in-vivo MRI, completing a pipeline for histology to in-vivo MRI registration. When correlating relaxometry and DTI with neuronal density and morphology in the temporal lobe neocortex, I found T1 to be a predictor of neuronal density in the neocortical GM and demonstrated that employing multi-parametric MRI (combining T1 and FA together) provided a significantly better fit than each parameter alone in predicting density of neurons. This work was the first to relate in-vivo T1 and FA values to the proportion of neurons in GM. When investigating these quantitative multimodal parameters with histological features within the hippocampal subfields, I demonstrated that MD correlates with neuronal density and size, and can act as a marker for neuron integrity within the hippocampus. More importantly, this work was the first to highlight the potential of subfield relaxometry and diffusion parameters (mainly T2 and MD) as well as volumetry in predicting the extent of cell loss per subfield pre-operatively, with a precision so far unachievable. These results suggest that high-resolution quantitative MRI sequences could impact clinical practice for pre-operative evaluation and prediction of surgical outcomes of intractable epilepsy

Cortical Morphology and MRI Signal Intensity Analysis in Paediatric Epilepsy

Epilepsy encompasses a great variety of aetiologies, and as such is not a single disease but a group of diseases characterised by unprovoked seizures.The primary aim of the work presented in this thesis was to use multimodal structural imaging to improve understanding of epilepsy related brain pathology, both the epileptogenic lesions themselves and extralesional pathology, in order to improve pre-surgical planning in medicationresistant epilepsy and improve understanding of the underlying pathogenic mechanisms. The work focuses on 2 epilepsy aetiologies: focal cortical dysplasia (FCD) (chapters 2 and 3) and mesial temporal lobe epilepsy (chapters 4 & 5). Chapter 2 of this thesis develops surface-based, structural MRI post-processing techniques that can be applied to clinical T1 and FLAIR images to complement current MRI-based diagnosis of focal cortical dysplasias. Chapter 3 uses the features developed in Chapter 2 within a machine learning framework to automatically detect FCDs, obtaining 73% sensitivity using a neural network. Chapter 4 develops an in vivo method to explore neocortical gliosis in adults with TLE, while Chapter 5 applies this method to a paediatric cohort. Finally, the concluding chapter discusses contributions, main limitations and outlines options for future research