Tractographie globale sous contraintes anatomiques

This work aims at developing a method inferring white matter fibers reconstructed using a global spin-glass approach constrained by anatomical prior knowledge. Unlike usual methods building fibers independently from one another, our markovian approach reconstructs the whole tractogram in an unique process by minimizing the global energy depending on the spin glass configuration (position, orientation, length and connection(s)) and the match with the local diffusion process in order to increase the robustness and the accuracy of the algorithm and the anatomical reliability of the reconstructed fibers. Thus, the work done during this PhD, along with the development of the global tractography algorithm, consisted in studying the feasibility of the anatomical prior knowledge integration stemming from the T1 weighted MRI and from new diffusion MRI microstructure approaches providing microstructural information of the surrounding tissue. In particular, the algorithm was built to allow a high fiber curvature when getting closer to the cortical ribbon and thus enabling the connection not only at the end of the gyri but also on their sides. The NODDI (Neurite Orientation Dispersion and Density Imaging) model has become more and more popular during the past years thanks to its capability to be used in clinical routine and allows to quantify neurite density and axons angular dispersion. A high dispersion means the existence of different fibers population or a high curvature of a fascicle within a voxel. Thus, the orientation dispersion has been used in our global tractography framework to release the curvature constraint near the cerebral cortex when the angular dispersion is high, allowing fibers to orientate collinear to the local normal to the cortical surface. However, this constraint is removed if the angular dispersion stays low, meaning a low curvature fiber trajectory following the example of the fibers projecting to the end of a gyrus or the U-fibers. The performances of this new tractography approach constrained by anatomical prior knowledge have been evaluated on simulated data, and tested on high resolution post-mortem MRI acquisitions and millimetric resolution in vivo MRI acquisitions. In parallel of this methodological development, a study about local-regional correlations between neurite density and cerebral activation on the cortical surface has been made. This study has been conducted on the healthy volunteers cohort scanned in the frame of the European CONNECT project including anatomical, diffusion and functional data. The anatomical data has been used to extract the pial surface and an individual parcellation on the cortical surface for each volunteer, the diffusion data has been used to evaluate the individual maps of neurite density within the cortical ribbon and the functional data from the BOLD (Blood Oxygen Level Dependent) effect has been used to calculate the individual z-scores of the general linear model for specific contrasts investigating the motor, language and visual networks. A co-localization of neurite density and activation peaks has been observed, which might indicate an increase of the neurite density within functional networks in order to increase its efficiency. This study also corroborates the lateralization of the language functional network and the motor one, in good agreement with the population lateralization, while an increase of the neurite density in the visual cortex has been observed which might be correlated to the results of visuo-spatial attention studies described in the literature on the non-human primate.Ce travail vise au développement d’une méthode d’inférence des fibres de la substance blanche cérébrale fondée sur l’utilisation d’une approche globale de type « verres de spins » sous contraintes anatomiques. Contrairement aux méthodes classiques reconstituant les fibres indépendamment les unes des autres, cette approche markovienne reconstruit l’ensemble des fibres dans un unique processus de minimisation d’une énergie globale dépendant de la configuration des spins (position, orientation, longueur et connexion(s)) et de leur adéquation avec le modèle local du processus de diffusion, afin d'améliorer la robustesse et la réalité anatomique des fibres reconstruites. Le travail mené dans le cadre de cette thèse a donc consisté, en plus du développement de l’algorithme de tractographie, à étudier la possibilité de le contraindre à l’aide d’a priori anatomiques provenant de l’imagerie anatomique pondérée en T1 et des nouvelles approches de microscopie par IRM de diffusion fournissant des informations de nature micro-structurelle sur le tissu. En particulier, l’algorithme a été conçu pour autoriser une forte courbure des fibres à l’approche du ruban cortical et permettre leur connexion au sommet des gyri, mais également sur leurs flancs. Le modèle NODDI (Neurite Orientation Dispersion and Density Imaging) a gagné en popularité au cours des dernières années grâce à sa compatibilité avec une utilisation en routine clinique et permet de quantifier la densité neuritique et la dispersion angulaire des axones. Une forte dispersion traduit l’existence de populations de fibres d’orientations différentes ou une forte courbure d’un même faisceau de fibres au sein d'un voxel. Elle est donc exploitée pour relâcher la contrainte de faible courbure à proximité du cortex cérébral dans notre approche de tractographie globale lorsque cette dispersion angulaire est forte, permettant aux fibres de s'orienter par rapport à la normale locale au cortex. Cette contrainte est en revanche supprimée si la dispersion angulaire reste faible, indiquant une trajectoire à plus faible courbure, à l’instar des fibres se projetant dans le fond du gyrus ou des fibres en U. Les performances de cette nouvelle approche de tractographie sous contraintes anatomiques ont été évaluées à partir de données simulées, et ont été testées sur des données IRM post-mortem de très haute résolution et sur des données IRM in vivo de résolution millimétrique. En parallèle de ce développement méthodologique, une étude des corrélats locaux-régionaux de la densité neuritique et de l’activation cérébrale à la surface du cortex a été réalisée. L'étude a été menée sur la cohorte de sujets sains scannés dans le cadre du projet européen CONNECT dotée de données anatomiques, de diffusion et fonctionnelles reposant sur l’utilisation de paradigmes explorant en particulier les réseaux de la motricité, du langage et de la vision. Les données anatomiques ont permis d’extraire la surface piale et une parcellisation surfacique du cortex de chaque individu, les données de diffusion ont permis l’évaluation des cartographies individuelles de la densité neuritique au sein du ruban cortical et les données fonctionnelles du phénomène BOLD (Blood Oxygen Level Dependent) ont permis le calcul des cartographies individuelles des z-scores du modèle linéaire général pour différents contrastes. Une colocalisation des maxima de la densité neuritique et des pics d'activation a été observée, pouvant être interprétée comme une augmentation de la densité neuritique au sein des réseaux fonctionnels afin d'en améliorer l'efficacité. L’étude a également corroboré la latéralisation du réseau fonctionnel du langage et de la motricité, en accord avec la latéralisation de la population scannée tandis qu'une augmentation de la densité neuritique dans le cortex visuel droit a été observée pouvant être corrélée à des résultats d’étude de l’attention visuo-spatiale reportée dans la littérature chez le primate non-humain