K-Surfer: a KNIME extension for the management and analysis of Human brain MRI FreeSurfer/FSL data

Human brain imaging techniques, such as Magnetic Resonance Imaging (MRI) or Diffusion Tensor Imaging (DTI), have been established as scientific and diagnostic tools and their adoption is growing in popularity. Statistical methods, machine learning and data mining algorithms have successfully been adopted to extract predictive and descriptive models from neuroimage data. However, the knowledge discovery process typically requires also the adoption of pre-processing, post-processing and visualisation techniques in complex data workflows. Currently, a main problem for the integrated preprocessing and mining of MRI data is the lack of comprehensive platforms able to avoid the manual invocation of preprocessing and mining tools, that yields to an error-prone and inefficient process. In this work we present K-Surfer, a novel plug-in of the Konstanz Information Miner (KNIME) workbench, that automatizes the preprocessing of brain images and leverages the mining capabilities of KNIME in an integrated way. K-Surfer supports the importing, filtering, merging and pre-processing of neuroimage data from FreeSurfer, a tool for human brain MRI feature extraction and interpretation. K-Surfer automatizes the steps for importing FreeSurfer data, reducing time costs, eliminating human errors and enabling the design of complex analytics workflow for neuroimage data by leveraging the rich functionalities available in the KNIME workbench.

A new combined surface and volume registration

3D registration of brain MRI data is vital for many medical imaging applications. However, purely intensitybased approaches for inter-subject matching of brain structure are generally inaccurate in cortical regions, due to the highly complex network of sulci and gyri, which vary widely across subjects. Here we combine a surfacebased cortical registration with a 3D fluid one for the first time, enabling precise matching of cortical folds, but allowing large deformations in the enclosed brain volume, which guarantee diffeomorphisms. This greatly improves the matching of anatomy in cortical areas. The cortices are segmented and registered with the software Freesurfer. The deformation field is initially extended to the full 3D brain volume using a 3D harmonic mapping that preserves the matching between cortical surfaces. Finally, these deformation fields are used to initialize a 3D Riemannian fluid registration algorithm, that improves the alignment of subcortical brain regions. We validate this method on an MRI dataset from 92 healthy adult twins. Results are compared to those based on volumetric registration without surface constraints; the resulting mean templates resolve consistent anatomical features both subcortically and at the cortex, suggesting that the approach is well-suited for cross-subject integration of functional and anatomic data.

Cortical thickness changes in the non-lesioned hemisphere associated with non-paretic arm immobilization in modified CI therapy

Recent evidence suggests that immobilization of the upper limb for 2–3 weeks induces changes in cortical thickness as well as motor performance. In constraint induced (CI) therapy, one of the most effective interventions for hemiplegia, the non-paretic arm is constrained to enforce the use of the paretic arm in the home setting. With the present study we aimed to explore whether non-paretic arm immobilization in CI therapy induces structural changes in the non-lesioned hemisphere, and how these changes are related to treatment benefit. 31 patients with chronic hemiparesis participated in CI therapy with (N = 14) and without (N = 17) constraint. Motor ability scores were acquired before and after treatment. Diffusion tensor imaging (DTI) data was obtained prior to treatment. Cortical thickness was measured with the Freesurfer software. In both groups cortical thickness in the contralesional primary somatosensory cortex increased and motor function improved with the intervention. However the cortical thickness change was not associated with the magnitude of motor function improvement. Moreover, the treatment effect and the cortical thickness change were not significantly different between the constraint and the non-constraint groups. There was no correlation between fractional anisotropy changes in the non-lesioned hemisphere and treatment outcome. CI therapy induced cortical thickness changes in contralesional sensorimotor regions, but this effect does not appear to be driven by the immobilization of the non-paretic arm, as indicated by the absence of differences between the constraint and the non-constraint groups. Our data does not suggest that the arm immobilization used in CI therapy is associated with noticeable cortical thinning.

Cortical correlates of affective syndrome in dementia due to Alzheimer’s disease

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Investigation of the cingulate cortex in idiopathic generalized epilepsy

Studies using quantitative neuroimaging have shown subtle abnormalities in patients with idiopathic generalized epilepsy (IGE). These findings have several locations, but the midline parasagittal structures are most commonly implicated. The cingulate cortex is related and may be involved. The objective of the current investigation was to perform a comprehensive analysis of the cingulate cortex using multiple quantitative structural neuroimaging techniques. Thirty-two patients (18 women, 30 ± 10 years) and 36 controls (18 women, 32 ± 11 years) were imaged by 3 Tesla magnetic resonance imaging (MRI). A volumetric three-dimensional (3D) sequence was acquired and used for this investigation. Regions-of-interest were selected and voxel-based morphometry (VBM) analyses compared the cingulate cortex of the two groups using Statistical Parametric Mapping (SPM8) and VBM8 software. Cortical analyses of the cingulate gyrus was performed using Freesurfer. Images were submitted to automatic processing using built-in routines and recommendations. Structural parameters were extracted for individual analyses, and comparisons between groups were restricted to the cingulate gyrus. Finally, shape analyses was performed on the anterior rostral, anterior caudal, posterior, and isthmus cingulate using spherical harmonic description (SPHARM). VBM analyses of cingulate gyrus showed areas of gray matter atrophy, mainly in the anterior cingulate gyrus (972 mm(3) ) and the isthmus (168 mm(3) ). Individual analyses of the cingulate cortex were similar between patients with IGE and controls. Surface-based comparisons revealed abnormalities located mainly in the posterior cingulate cortex (718.12 mm(2) ). Shape analyses demonstrated a predominance of anterior and posterior cingulate abnormalities. This study suggests that patients with IGE have structural abnormalities in the cingulate gyrus mainly localized at the anterior and posterior portions. This finding is subtle and variable among patients.

Confronto di dati volumetrici di neuroimaging RM, acquisiti con differenti sonde a radiofrequenza

Il lavoro di questa tesi si propone di esaminare i dati di immagini cerebrali ricevuti da due differentii bobine RF per l'imaging in risonanza magnetica. La strumentazione utilizzata é: un apparecchio di RMN che utilizza uno scanner a 1.5T (GE Medical System signa HDx 15) e due differenti bobine a radiofrequenza: 1) 8-channel brain phased array coil GE (1.5T HD 8 Channel High Resolution Head Array per il GE HDx MR System); 2) GE Quad HEAD Birdcage Coil. I software utilizzati invece sono stati quattro, due per la segmentazione e la parcellizzazione dei dati dalle acquisizioni di MRI (FSL e Freesurfer), e due per l'analisi dei dati (SPSS e Microsoft Excel). I tool utilizzati di FSL sono stati SIENA, per un'indagine sulla variazione percentile della totalitá del volume cerebrale; SIENAX invece permette una segmentazione di volumi di 6 sotto regioni: sostanza grigia (GREY), sostanza bianca (WHITE), totalitá del cervello (BRAIN), liquor cerebrospinale dei ventricoli (vcsf), scaling volumetrico (vscaling), sostanza grigia periferica (pgrey). Freesurfer invece parcellizza la corteccia cerebrale e segmenta le zone della sostanza grigia profonda cerebrale. Lo scopo ultimo era quello di analizzare come la differenza di segnale acquisito dalle due bobine variasse i risultati delle analisi volumetriche delle immagini e vedere se il t-test evidenziasse variazioni sistematiche. Questa analisi aveva come scopo quello di determinare la possibilità di confrontare i soggetti i cui dati sono stati ottenuti con due bobine differenti. I dati analizzati non hanno evidenziato particolari differenze tra le due bobine, se non per i valori del liquor dei ventricoli che sono risultati variare, verosimilmente, per i processi fisiologici di atrofia della sostanza grigia cerebrale che avvengono nel tempo intercorso tra un'acquisizione e l'altra.

Delay of cortical thinning in very preterm born children

BACKGROUND: Cortical gray matter thinning occurs during childhood due to pruning of inefficient synaptic connections and an increase in myelination. Preterms show alterations in brain structure, with prolonged maturation of the frontal lobes, smaller cortical volumes and reduced white matter volume. These findings give rise to the question if there is a differential influence of age on cortical thinning in preterms compared to controls. AIMS: To investigate the relationship between age and cortical thinning in school-aged preterms compared to controls. STUDY DESIGN AND OUTCOME MEASURES: The automated surface reconstruction software FreeSurfer was applied to obtain measurements of cortical thickness based on T1-weighted MRI images. SUBJECTS: Forty-one preterms (<32weeks gestational age and/or <1500g birth weight) and 30 controls were included in the study (7-12years). RESULTS: In preterms, age correlated negatively with cortical thickness in right frontal, parietal and inferior temporal regions. Furthermore, young preterms showed a thicker cortex compared to old preterms in bilateral frontal, parietal and temporal regions. In controls, age was not associated with cortical thickness. CONCLUSION: In preterms, cortical thinning still seems to occur between the age of 7 and 12years, mainly in frontal and parietal areas whereas in controls, a substantial part of cortical thinning appears to be completed before they reach the age of 7years. These data indicate slower cortical thinning in preterms than in controls.

Delay of cortical thinning in very preterm born children

Objective: Cortical gray matter thinning takes place during childhood due to pruning of inefficient synaptic connections and an increase in myelination. Alterations in brain structure occur in very preterm born children with prolonged maturation of the frontal lobes and smaller cortical and white matter volume. These findings give rise to the question if age affects cortical thinning differently in very preterm born children compared to controls. The aim of the present study was to investigate the relationship between age and cortical thickness in very preterm born children when compared to controls. Participants and Methods: Forty-one very preterm born children (<32 weeks gestational age and/or < 1500 gram birth weight) and 30term born controls were included in the study (7-12 years). The automated surface reconstruction software FreeSurfer was applied to obtain measurements of cortical thickness based on T1-weighted MRI images. Results: Cortical thickness was lower in bilateral frontal and left parietal regions and higher in left temporal gyri in very preterm born children compared to controls. However, these differences depended on age. In very preterm born children, age correlated negatively with cortical thickness in right frontal, parietal and inferior temporal regions. Accordingly, cortical thickness was higher in young compared to old very preterm born children in bilateral frontal, parietal and temporal regions. In controls, age was not associated with cortical thickness. Conclusions: In very preterm born children, cortical thinning still occurs between the age of 7 and 12 years, mainly in frontal and parietal areas. In controls, however, a substantial part of cortical thinning appears to be completed in these regions before they reach the age of 7 years. These data indicate a delay in cortical thinning in very preterm born children.

Cortical thickness decreases with age in very preterm born school-children but not in term born controls

Background: Cortical gray matter thinning occurs during childhood due to pruning of inefficient synaptic connections and an increase in myelination. Preterms show alterations in brain structure, with prolonged maturation of the frontal lobes, smaller cortical volumes and reduced white matter volume. These findings give rise to the question if there is a differential influence of age on cortical thinning in preterms compared to controls. Aims: To investigate the relationship between age and cortical thickness in preterms when compared to controls. Study design and outcome measures: The automated surface reconstruction software FreeSurfer was applied to obtain measurements of cortical thickness based on T1-weighted MRI images. Subjects: Forty-one preterms (< 32 weeks gestational age and/or < 1500 gram birth weight) and 30 controls were included in the study (7-12 years). Results: Cortical thickness was lower in bilateral frontal and left parietal regions and higher in left temporal gyri in preterms compared to controls. However, these differences depended on age. In preterms, age correlated negatively with cortical thickness in right frontal, parietal and inferior temporal regions. Accordingly, cortical thickness was higher in young compared to old preterms in bilateral frontal, parietal and temporal regions. In controls, age was not associated with cortical thickness. Conclusion: In preterms, cortical thinning still seems to occur between the age of 7 and 12 years, mainly in frontal and parietal areas whereas in controls, a substantial part of cortical thinning appears to be completed before they reach the age of 7 years. These data indicate slower cortical thinning in preterms than in controls.

Herramienta de modelado disfuncional tridimensional basado en estudios de neuroimagen

El modelado disfuncional basado en estudios de neuroimagen mejora la comprensión de los cambios estructurales provocados ante la presencia de lesiones cerebrales. Actualmente, existen numerosas herramientas para el análisis y procesado de estudios de neuroimagen. Algunas de ellas, como el 3D Slicer, BrainVoyager y el FreeSurfer permiten la creación y navegación sobre modelos tridimensionales cerebrales sin alteraciones estructurales. Sin embargo, no se han detectado herramientas que permitan modelar tridimensionalmente lesiones a partir de estudios de neuroimagen, concretamente de estudios de resonancia magnética. El objetivo de este trabajo es el diseño de una metodología que permite la creación de este tipo de modelos y su visualización y navegación.

Inter-site and inter-scanner diffusion MRI data harmonization.

We propose a novel method to harmonize diffusion MRI data acquired from multiple sites and scanners, which is imperative for joint analysis of the data to significantly increase sample size and statistical power of neuroimaging studies. Our method incorporates the following main novelties: i) we take into account the scanner-dependent spatial variability of the diffusion signal in different parts of the brain; ii) our method is independent of compartmental modeling of diffusion (e.g., tensor, and intra/extra cellular compartments) and the acquired signal itself is corrected for scanner related differences; and iii) inter-subject variability as measured by the coefficient of variation is maintained at each site. We represent the signal in a basis of spherical harmonics and compute several rotation invariant spherical harmonic features to estimate a region and tissue specific linear mapping between the signal from different sites (and scanners). We validate our method on diffusion data acquired from seven different sites (including two GE, three Philips, and two Siemens scanners) on a group of age-matched healthy subjects. Since the extracted rotation invariant spherical harmonic features depend on the accuracy of the brain parcellation provided by Freesurfer, we propose a feature based refinement of the original parcellation such that it better characterizes the anatomy and provides robust linear mappings to harmonize the dMRI data. We demonstrate the efficacy of our method by statistically comparing diffusion measures such as fractional anisotropy, mean diffusivity and generalized fractional anisotropy across multiple sites before and after data harmonization. We also show results using tract-based spatial statistics before and after harmonization for independent validation of the proposed methodology. Our experimental results demonstrate that, for nearly identical acquisition protocol across sites, scanner-specific differences can be accurately removed using the proposed method.

Facteurs de risque de développement de troubles intériorisés : études en Imagerie par Résonance Magnétique structurelle

Plusieurs facteurs de risque de développement de troubles intériorisés, tels que les troubles d’anxiété et de l’humeur, ont été identifiés dans la littérature. Les deux plus importants facteurs de risques regroupent l’adversité vécue durant l’enfance (par exemple la maltraitance) et le risque parental (c’est-à-dire la présence d’un trouble intériorisé chez l’un ou les deux parents). Ces facteurs de risque ont été liés à des changements neuroanatomiques similaires à ceux observés en lien avec les troubles intériorisés. Ainsi, en présence de ces facteurs de risque, des anomalies anatomiques pourraient laisser présager l’apparition prochaine d’une symptomatologie de troubles intériorisés chez des individus encore asymptomatiques. Chez les quelques populations de jeunes investiguées, les participants présentaient des comorbidités et/ou étaient sous médication, ce qui rend difficile l’interprétation des atteintes cérébrales observées. Ce travail de thèse s’est intéressé aux liens entre ces deux facteurs de risque et les substrats neuroanatomiques associés à chacun d’eux, chez des adolescents asymptomatiques et n’étant sous aucune médication. Une première étude a examiné le lien entre le niveau de pratiques parentales coercitives et le niveau de symptômes d’anxiété, mesurés de manière longitudinale depuis la naissance, et les différences neuroanatomiques observées à l’adolescence (voir Chapitre 2). Une deuxième étude a examiné le lien entre le risque parental de développer des troubles d’anxiété et les différences neuroanatomiques observées à l’adolescence (voir Chapitre 3). Une troisième étude s’est intéressée au lien entre le risque parental de développer un trouble de dépression ou un trouble bipolaire et les différences neuroanatomiques observées à l’adolescence (voir Chapitre 4). Les résultats démontrent des différences de volume et/ou d’épaisseur corticale dans plusieurs structures clés impliquées dans le traitement et la régulation des émotions. C’est le cas du cortex préfrontal, de l’amygdale, de l’hippocampe et du striatum. Ces résultats suggèrent que certaines des différences neuroanatomiques observées dans les troubles intériorisés peuvent être présentes avant que le trouble ne se manifeste, et représenter des marqueurs neuronaux du risque de développer le trouble. Les implications théoriques et les limites de ces trois études sont finalement discutées dans le Chapitre 5.

Facteurs de risque de développement de troubles intériorisés : études en Imagerie par Résonance Magnétique structurelle

Plusieurs facteurs de risque de développement de troubles intériorisés, tels que les troubles d’anxiété et de l’humeur, ont été identifiés dans la littérature. Les deux plus importants facteurs de risques regroupent l’adversité vécue durant l’enfance (par exemple la maltraitance) et le risque parental (c’est-à-dire la présence d’un trouble intériorisé chez l’un ou les deux parents). Ces facteurs de risque ont été liés à des changements neuroanatomiques similaires à ceux observés en lien avec les troubles intériorisés. Ainsi, en présence de ces facteurs de risque, des anomalies anatomiques pourraient laisser présager l’apparition prochaine d’une symptomatologie de troubles intériorisés chez des individus encore asymptomatiques. Chez les quelques populations de jeunes investiguées, les participants présentaient des comorbidités et/ou étaient sous médication, ce qui rend difficile l’interprétation des atteintes cérébrales observées. Ce travail de thèse s’est intéressé aux liens entre ces deux facteurs de risque et les substrats neuroanatomiques associés à chacun d’eux, chez des adolescents asymptomatiques et n’étant sous aucune médication. Une première étude a examiné le lien entre le niveau de pratiques parentales coercitives et le niveau de symptômes d’anxiété, mesurés de manière longitudinale depuis la naissance, et les différences neuroanatomiques observées à l’adolescence (voir Chapitre 2). Une deuxième étude a examiné le lien entre le risque parental de développer des troubles d’anxiété et les différences neuroanatomiques observées à l’adolescence (voir Chapitre 3). Une troisième étude s’est intéressée au lien entre le risque parental de développer un trouble de dépression ou un trouble bipolaire et les différences neuroanatomiques observées à l’adolescence (voir Chapitre 4). Les résultats démontrent des différences de volume et/ou d’épaisseur corticale dans plusieurs structures clés impliquées dans le traitement et la régulation des émotions. C’est le cas du cortex préfrontal, de l’amygdale, de l’hippocampe et du striatum. Ces résultats suggèrent que certaines des différences neuroanatomiques observées dans les troubles intériorisés peuvent être présentes avant que le trouble ne se manifeste, et représenter des marqueurs neuronaux du risque de développer le trouble. Les implications théoriques et les limites de ces trois études sont finalement discutées dans le Chapitre 5.

Brain development during adolescence: A mixed-longitudinal investigation of cortical thickness, surface area, and volume.

What we know about cortical development during adolescence largely stems from analyses of cross-sectional or cohort-sequential samples, with few studies investigating brain development using a longitudinal design. Further, cortical volume is a product of two evolutionarily and genetically distinct features of the cortex - thickness and surface area, and few studies have investigated development of these three characteristics within the same sample. The current study examined maturation of cortical thickness, surface area and volume during adolescence, as well as sex differences in development, using a mixed longitudinal design. 192 MRI scans were obtained from 90 healthy (i.e., free from lifetime psychopathology) adolescents (11-20 years) at three time points (with different MRI scanners used at time 1 compared to 2 and 3). Developmental trajectories were estimated using linear mixed models. Non-linear increases were present across most of the cortex for surface area. In comparison, thickness and volume were both characterised by a combination of non-linear decreasing and increasing trajectories. While sex differences in volume and surface area were observed across time, no differences in thickness were identified. Furthermore, few regions exhibited sex differences in the cortical development. Our findings clearly illustrate that volume is a product of surface area and thickness, with each exhibiting differential patterns of development during adolescence, particularly in regions known to contribute to the development of social-cognition and behavioral regulation. These findings suggest that thickness and surface area may be driven by different underlying mechanisms, with each measure potentially providing independent information about brain development. Hum Brain Mapp, 2016. © 2016 Wiley Periodicals, Inc.