SPM to the heart: mapping of 4D continuous velocities for motion abnormality quantification

International audienceThis paper proposes to apply parallel transport and statistical atlas techniques to quantify 4D myocardial motion abnormalities. We take advantage of our previous work on cardiac motion , which provided a continuous spatiotemporal representation of velocities, to interpolate and reorient cardiac motion fields to an unbiased reference space. Abnormal motion is quantified using SPM analysis on the velocity fields, which includes a correction based on random field theory to compensate for the spatial smoothness of the velocity fields. This paper first introduces the imaging pipeline for constructing a continuous 4D velocity atlas. This atlas is then applied to quantify abnormal motion patterns in heart failure patients

Development of an MRI Template and Analysis Pipeline for the Spinal Cord and Application in Patients with Spinal Cord Injury

La moelle épinière est un organe fondamental du corps humain. Étant le lien entre le cerveau et le système nerveux périphérique, endommager la moelle épinière, que ce soit suite à un trauma ou une maladie neurodégénérative, a des conséquences graves sur la qualité de vie des patients. En effet, les maladies et traumatismes touchant la moelle épinière peuvent affecter l’intégrité des neurones et provoquer des troubles neurologiques et/ou des handicaps fonctionnels. Bien que de nombreuses voies thérapeutiques pour traiter les lésions de la moelle épinière existent, la connaissance de l’étendue des dégâts causés par ces lésions est primordiale pour améliorer l’efficacité de leur traitement et les décisions cliniques associées. L’imagerie par résonance magnétique (IRM) a démontré un grand potentiel pour le diagnostic et pronostic des maladies neurodégénératives et traumas de la moelle épinière. Plus particulièrement, l’analyse par template de données IRM du cerveau, couplée à des outils de traitement d’images automatisés, a permis une meilleure compréhension des mécanismes sous-jacents de maladies comme l’Alzheimer et la Sclérose en Plaques. Extraire automatiquement des informations pertinentes d’images IRM au sein de régions spécifiques de la moelle épinière présente toutefois de plus grands défis que dans le cerveau. Il n’existe en effet qu’un nombre limité de template de la moelle épinière dans la littérature, et aucun ne couvre toute la moelle épinière ou n’est lié à un template existant du cerveau. Ce manque de template et d’outils automatisés rend difficile la tenue de larges études d’analyse de la moelle épinière sur des populations variées. L’objectif de ce projet est donc de proposer un nouveau template IRM couvrant toute la moelle épinière, recalé avec un template existant du cerveau, et intégrant des atlas de la structure interne de la moelle épinière (e.g., matière blanche et grise, tracts de la matière blanche). Ce template doit venir avec une série d’outils automatisés permettant l’extraction d’information IRM au sein de régions spécifiques de la moelle épinière. La question générale de recherche de ce projet est donc « Comment créer un template générique de la moelle épinière, qui permettrait l’analyse non biaisée et reproductible de données IRM de la moelle épinière ? ». Plusieurs contributions originales ont été proposées pour répondre à cette question et vont être décrites dans les prochains paragraphes. La première contribution de ce projet est le développement du logiciel Spinal Cord Toolbox (SCT). SCT est un logiciel open-source de traitement d’images IRM multi-parametrique de la moelle épinière (De Leener, Lévy, et al., 2016). Ce logiciel intègre notamment des outils pour la détection et la segmentation automatique de la moelle épinière et de sa structure interne (i.e., matière blanche et matière grise), l’identification et la labellisation des niveaux vertébraux, le recalage d’images IRM multimodales sur un template générique de la moelle épinière (précédemment le template MNI-Poly-AMU, maintenant le template PAM50, proposé içi). En se basant sur un atlas de la moelle, SCT intègre également des outils pour extraire des données IRM de régions spécifiques de la moelle épinière, comme la matière blanche et grise et les tracts de la matière blanche, ainsi que sur des niveaux vertébraux spécifiques. D’autres outils additionnels ont aussi été proposés, comme des outils de correction de mouvement et de traitement basiques d’images appliqués le long de la moelle épinière. Chaque outil intégré à SCT a été validé sur un jeu de données multimodales. La deuxième contribution de ce projet est le développement d’une nouvelle méthode de recalage d’images IRM de la moelle épinière (De Leener, Mangeat, et al., 2017). Cette méthode a été développée pour un usage particulier : le redressement d’images IRM de la moelle épinière, mais peut également être utilisé pour recaler plusieurs images de la moelle épinière entre elles, tout en tenant compte de la distribution vertébrale de chaque sujet. La méthode proposée se base sur une approximation globale de la courbure de la moelle épinière dans l’espace et sur la résolution analytique des champs de déformation entre les deux images. La validation de cette nouvelle méthode a été réalisée sur une population de sujets sains et de patients touchés par une compression de la moelle épinière. La contribution majeure de ce projet est le développement d’un système de création de template IRM de la moelle épinière et la proposition du template PAM50 comme template de référence pour les études d’analyse par template de données IRM de la moelle épinière. Le template PAM50 a été créé à partir d’images IRM tiré de 50 sujets sains, et a été généré en utilisant le redressement d’images présenté ci-dessus et une méthode de recalage d’images itératif non linéaire, après plusieurs étapes de prétraitement d’images. Ces étapes de prétraitement incluent la segmentation automatique de la moelle épinière, l’extraction manuelle du bord antérieur du tronc cérébral, la détection et l’identification des disques intervertébraux, et la normalisation d’intensité le long de la moelle. Suite au prétraitement, la ligne centrale moyenne de la moelle et la distribution vertébrale ont été calculées sur la population entière de sujets et une image initiale de template a été générée. Après avoir recalé toutes les images sur ce template initial, le template PAM50 a été créé en utilisant un processus itératif de recalage d’image, utilisé pour générer des templates de cerveau. Le PAM50 couvre le tronc cérébral et la moelle épinière en entier, est disponible pour les contrastes IRM pondérés en T1, T2 et T2*, et intègre des cartes probabilistes et atlas de la structure interne de la moelle épinière. De plus, le PAM50 a été recalé sur le template ICBM152 du cerveau, permettant ainsi la tenue d’analyse par template simultanément dans le cerveau et dans la moelle épinière. Finalement, plusieurs résultats complémentaires ont été présentés dans cette dissertation. Premièrement, une étude de validation de la répétabilité et reproductibilité de mesures de l’aire de section de la moelle épinière a été menée sur une population de patients touchés par la sclérose en plaques. Les résultats démontrent une haute fiabilité des mesures ainsi que la possibilité de détecter des changements très subtiles de l’aire de section transverse de la moelle, importants pour mesurer l’atrophie de la moelle épinière précoce due à des maladies neurodégénératives comme la sclérose en plaques. Deuxièmement, un nouveau biomarqueur IRM des lésions de la moelle épinière a été proposé, en collaboration avec Allan Martin, de l’Université de Toronto. Ce biomarqueur, calculé à partir du ratio d’intensité entre la matière blanche et grise sur des images IRM pondérées en T2*, utilise directement les développements proposés dans ce projet, notamment en utilisant le recalage du template de la moelle épinière et les méthodes de segmentation de la moelle. La faisabilité d’extraire des mesures de données IRM multiparamétrique dans des régions spécifiques de la moelle épinière a également été démontrée, permettant d’améliorer le diagnostic et pronostic de lésions et compression de la moelle épinière. Finalement, une nouvelle méthode d’extraction de la morphométrie de la moelle épinière a été proposée et utilisée sur une population de patients touchés par une compression asymptomatique de la moelle épinière, démontrant de grandes capacités de diagnostic (> 99%). Le développement du template PAM50 comble le manque de template de la moelle épinière dans la littérature mais présente cependant plusieurs limitations. En effet, le template proposé se base sur une population de 50 sujets sains et jeunes (âge moyen = 27 +- 6.5) et est donc biaisée vers cette population particulière. Adapter les analyses par template pour un autre type de population (âge, race ou maladie différente) peut être réalisé directement sur les méthodes d’analyse mais aussi sur le template en lui-même. Tous le code pour générer le template a en effet été mis en ligne (https://github.com/neuropoly/template) pour permettre à tout groupe de recherche de développer son propre template. Une autre limitation de ce projet est le choix d’un système de coordonnées basé sur la position des vertèbres. En effet, les vertèbres ne représentent pas complètement le caractère fonctionnel de la moelle épinière, à cause de la différence entre les niveaux vertébraux et spinaux. Le développement d’un système de coordonnées spinal, bien que difficile à caractériser dans des images IRM, serait plus approprié pour l’analyse fonctionnelle de la moelle épinière. Finalement, il existe encore de nombreux défis pour automatiser l’ensemble des outils développés dans ce projet et les rendre robuste pour la majorité des contrastes et champs de vue utilisés en IRM conventionnel et clinique. Ce projet a présenté plusieurs développements importants pour l’analyse de données IRM de la moelle épinière. De nombreuses améliorations du travail présenté sont cependant requises pour amener ces outils dans un contexte clinique et pour permettre d’améliorer notre compréhension des maladies affectant la moelle épinière. Les applications cliniques requièrent notamment l’amélioration de la robustesse et de l’automatisation des méthodes d’analyse d’images proposées. La caractérisation de la structure interne de la moelle épinière, incluant la matière blanche et la matière grise, présente en effet de grands défis, compte tenu de la qualité et la résolution des images IRM standard acquises en clinique. Les outils développés et validés au cours de ce projet ont un grand potentiel pour la compréhension et la caractérisation des maladies affectant la moelle épinière et aura un impact significatif sur la communauté de la neuroimagerie.----------ABSTRACT The spinal cord plays a fundamental role in the human body, as part of the central nervous system and being the vector between the brain and the peripheral nervous system. Damaging the spinal cord, through traumatic injuries or neurodegenerative diseases, can significantly affect the quality of life of patients. Indeed, spinal cord injuries and diseases can affect the integrity of neurons, and induce neurological impairments and/or functional disabilities. While various treatment procedures exist, assessing the extent of damages and understanding the underlying mechanisms of diseases would improve treatment efficiency and clinical decisions. Over the last decades, magnetic resonance imaging (MRI) has demonstrated a high potential for the diagnosis and prognosis of spinal cord injury and neurodegenerative diseases. Particularly, template-based analysis of brain MRI data has been very helpful for the understanding of neurological diseases, using automated analysis of large groups of patients. However, extracting MRI information within specific regions of the spinal cord with minimum bias and using automated tools is still a challenge. Indeed, only a limited number of MRI template of the spinal cord exists, and none covers the full spinal cord, thereby preventing large multi-centric template-based analysis of the spinal cord. Moreover, no template integrates both the spinal cord and the brain region, thereby preventing simultaneous cerebrospinal studies. The objective of this project was to propose a new MRI template of the full spinal cord, which allows simultaneous brain and spinal cord studies, that integrates atlases of the spinal cord internal structures (e.g., white and gray matter, white matter pathways) and that comes with tools for extracting information within these subregions. More particularly, the general research question of the project was “How to create generic MRI templates of the spinal cord that would enable unbiased and reproducible template-based analysis of spinal cord MRI data?”. Several original contributions have been made to answer this question and to enable template-based analysis of spinal cord MRI data. The first contribution was the development of the Spinal Cord Toolbox (SCT), a comprehensive and open-source software for processing multi-parametric MRI data of the spinal cord (De Leener, Lévy, et al., 2016). SCT includes tools for the automatic segmentation of the spinal cord and its internal structure (white and gray matter), vertebral labeling, registration of multimodal MRI data (structural and non-structural) on a spinal cord MRI template (initially the MNI-Poly-AMU template, later the PAM50 template), co-registration of spinal cord MRI images, as well as the robust extraction of MRI metric within specific regions of the spinal cord (i.e., white and gray matter, white matter tracts, gray matter subregions) and specific vertebral levels using a spinal cord atlas (Lévy et al., 2015). Additional tools include robust motion correction and image processing along the spinal cord. Each tool included in SCT has been validated on a multimodal dataset. The second contribution of this project was the development of a novel registration method dedicated to spinal cord images, with an interest in the straightening of the spinal cord, while preserving its topology (De Leener, Mangeat et al., 2017). This method is based on the global approximation of the spinal cord and the analytical computation of deformation fields perpendicular to the centerline. Validation included calculation of distance measurements after straightening on a population of healthy subjects and patients with spinal cord compression. The major contribution of this project was the development of a framework for generating MRI template of the spinal cord and the PAM50 template, an unbiased and symmetrical MRI template of the brainstem and full spinal cord. Based on 50 healthy subjects, the PAM50 template was generated using an iterative nonlinear registration process, after applying normalization and straightening of all images. Pre-processing included segmentation of the spinal cord, manual delineation of the brainstem anterior edge, detection and identification of intervertebral disks, and normalization of intensity along the spinal cord. Next, the average centerline and vertebral distribution was computed to create an initial straight template space. Then, all images were registered to the initial template space and an iterative nonlinear registration framework was applied to create the final symmetrical template. The PAM50 covers the brainstem and the full spinal cord, from C1 to L2, is available for T1-, T2- and T2*-weighted contrasts, and includes probabilistic maps of the white and the gray matter and atlases of the white matter pathways and gray matter subregions. Additionally, the PAM50 template has been merged with the ICBM152 brain template, thereby allowing for simultaneous cerebrospinal template-based analysis. Finally, several complementary results, focused on clinical validation and applications, are presented. First, a reproducibility and repeatability study of cross-sectional area measurements using SCT (De Leener, Granberg, Fink, Stikov, & Cohen-Adad, 2017) was performed on a Multiple Sclerosis population (n=9). The results demonstrated the high reproducibility and repeatability of SCT and its ability to detect very subtle atrophy of the spinal cord. Second, a novel biomarker of spinal cord injury has been proposed. Based on the T2*-weighted intensity ratio between the white and the gray matter, this new biomarker is computed by registering MRI images with the PAM50 template and extracting metrics using probabilistic atlases. Additionally, the feasibility of extracting multiparametric MRI metrics from subregions of the spinal cord has been demonstrated and the diagnostic potential of this approach has been assessed on a degenerative cervical myelopathy (DCM) population. Finally, a method for extracting shape morphometrics along the spinal cord has been proposed, including spinal cord flattening, indentation and torsion. These metrics demonstrated high capabilities for the diagnostic of asymptomatic spinal cord compression (AUC=99.8% for flattening, 99.3% for indentation, and 98.4% for torsion). The development of the PAM50 template enables unbiased template-based analysis of the spinal cord. However, the PAM50 template has several limitations. Indeed, the proposed template has been generated with multimodal MRI images from 50 healthy and young individuals (age = 27+/- 6.5 y.o.). Therefore, the template is specific to this particular population and could not be directly usable for age- or disease-specific populations. One solution is to open-source the templategeneration code so that research groups can generate and use their own spinal cord MRI template. The code is available on https://github.com/neuropoly/template. While this project introduced a generic referential coordinate system, based on vertebral levels and the pontomedullary junction as origin, one limitation is the choice of this coordinate system. Another coordinate system, based spinal segments would be more suitable for functional analysis. However, the acquisition of MRI images with high enough resolution to delineate the spinal roots is still challenging. Finally, several challenges in the automation of spinal cord MRI processing remains, including the robust detection and identification of vertebral levels, particularly in case of small fields-of-view. This project introduced key developments for the analysis of spinal cord MRI data. Many more developments are still required to bring them into clinics and to improve our understanding of diseases affecting the spinal cord. Indeed, clinical applications require the improvement of the robustness and the automation of the proposed processing and analysis tools. Particularly, the detection and segmentation of spinal cord structures, including vertebral labeling and white/gray matter segmentation, is still challenging, given the lowest quality and resolution of standard clinical MRI acquisition. The tools developed and validated here have the potential to improve our understanding and the characterization of diseases affecting the spinal cord and will have a significant impact on the neuroimaging community

Recent Advances in Signal Processing

The signal processing task is a very critical issue in the majority of new technological inventions and challenges in a variety of applications in both science and engineering fields. Classical signal processing techniques have largely worked with mathematical models that are linear, local, stationary, and Gaussian. They have always favored closed-form tractability over real-world accuracy. These constraints were imposed by the lack of powerful computing tools. During the last few decades, signal processing theories, developments, and applications have matured rapidly and now include tools from many areas of mathematics, computer science, physics, and engineering. This book is targeted primarily toward both students and researchers who want to be exposed to a wide variety of signal processing techniques and algorithms. It includes 27 chapters that can be categorized into five different areas depending on the application at hand. These five categories are ordered to address image processing, speech processing, communication systems, time-series analysis, and educational packages respectively. The book has the advantage of providing a collection of applications that are completely independent and self-contained; thus, the interested reader can choose any chapter and skip to another without losing continuity

Machine learning approaches to model cardiac shape in large-scale imaging studies

Recent improvements in non-invasive imaging, together with the introduction of fully-automated segmentation algorithms and big data analytics, has paved the way for large-scale population-based imaging studies. These studies promise to increase our understanding of a large number of medical conditions, including cardiovascular diseases. However, analysis of cardiac shape in such studies is often limited to simple morphometric indices, ignoring large part of the information available in medical images. Discovery of new biomarkers by machine learning has recently gained traction, but often lacks interpretability. The research presented in this thesis aimed at developing novel explainable machine learning and computational methods capable of better summarizing shape variability, to better inform association and predictive clinical models in large-scale imaging studies. A powerful and flexible framework to model the relationship between three-dimensional (3D) cardiac atlases, encoding multiple phenotypic traits, and genetic variables is first presented. The proposed approach enables the detection of regional phenotype-genotype associations that would be otherwise neglected by conventional association analysis. Three learning-based systems based on deep generative models are then proposed. In the first model, I propose a classifier of cardiac shapes which exploits task-specific generative shape features, and it is designed to enable the visualisation of the anatomical effect these features encode in 3D, making the classification task transparent. The second approach models a database of anatomical shapes via a hierarchy of conditional latent variables and it is capable of detecting, quantifying and visualising onto a template shape the most discriminative anatomical features that characterize distinct clinical conditions. Finally, a preliminary analysis of a deep learning system capable of reconstructing 3D high-resolution cardiac segmentations from a sparse set of 2D views segmentations is reported. This thesis demonstrates that machine learning approaches can facilitate high-throughput analysis of normal and pathological anatomy and of its determinants without losing clinical interpretability.Open Acces

Learning Biosignals with Deep Learning

The healthcare system, which is ubiquitously recognized as one of the most influential system in society, is facing new challenges since the start of the decade.The myriad of physiological data generated by individuals, namely in the healthcare system, is generating a burden on physicians, losing effectiveness on the collection of patient data. Information systems and, in particular, novel deep learning (DL) algorithms have been prompting a way to take this problem. This thesis has the aim to have an impact in biosignal research and industry by presenting DL solutions that could empower this field. For this purpose an extensive study of how to incorporate and implement Convolutional Neural Networks (CNN), Recursive Neural Networks (RNN) and Fully Connected Networks in biosignal studies is discussed. Different architecture configurations were explored for signal processing and decision making and were implemented in three different scenarios: (1) Biosignal learning and synthesis; (2) Electrocardiogram (ECG) biometric systems, and; (3) Electrocardiogram (ECG) anomaly detection systems. In (1) a RNN-based architecture was able to replicate autonomously three types of biosignals with a high degree of confidence. As for (2) three CNN-based architectures, and a RNN-based architecture (same used in (1)) were used for both biometric identification, reaching values above 90% for electrode-base datasets (Fantasia, ECG-ID and MIT-BIH) and 75% for off-person dataset (CYBHi), and biometric authentication, achieving Equal Error Rates (EER) of near 0% for Fantasia and MIT-BIH and bellow 4% for CYBHi. As for (3) the abstraction of healthy clean the ECG signal and detection of its deviation was made and tested in two different scenarios: presence of noise using autoencoder and fully-connected network (reaching 99% accuracy for binary classification and 71% for multi-class), and; arrhythmia events by including a RNN to the previous architecture (57% accuracy and 61% sensitivity). In sum, these systems are shown to be capable of producing novel results. The incorporation of several AI systems into one could provide to be the next generation of preventive medicine, as the machines have access to different physiological and anatomical states, it could produce more informed solutions for the issues that one may face in the future increasing the performance of autonomous preventing systems that could be used in every-day life in remote places where the access to medicine is limited. These systems will also help the study of the signal behaviour and how they are made in real life context as explainable AI could trigger this perception and link the inner states of a network with the biological traits.O sistema de saúde, que é ubiquamente reconhecido como um dos sistemas mais influentes da sociedade, enfrenta novos desafios desde o ínicio da década. A miríade de dados fisiológicos gerados por indíviduos, nomeadamente no sistema de saúde, está a gerar um fardo para os médicos, perdendo a eficiência no conjunto dos dados do paciente. Os sistemas de informação e, mais espcificamente, da inovação de algoritmos de aprendizagem profunda (DL) têm sido usados na procura de uma solução para este problema. Esta tese tem o objetivo de ter um impacto na pesquisa e na indústria de biosinais, apresentando soluções de DL que poderiam melhorar esta área de investigação. Para esse fim, é discutido um extenso estudo de como incorporar e implementar redes neurais convolucionais (CNN), redes neurais recursivas (RNN) e redes totalmente conectadas para o estudo de biosinais. Diferentes arquiteturas foram exploradas para processamento e tomada de decisão de sinais e foram implementadas em três cenários diferentes: (1) Aprendizagem e síntese de biosinais; (2) sistemas biométricos com o uso de eletrocardiograma (ECG), e; (3) Sistema de detecção de anomalias no ECG. Em (1) uma arquitetura baseada na RNN foi capaz de replicar autonomamente três tipos de sinais biológicos com um alto grau de confiança. Quanto a (2) três arquiteturas baseadas em CNN e uma arquitetura baseada em RNN (a mesma usada em (1)) foram usadas para ambas as identificações, atingindo valores acima de 90 % para conjuntos de dados à base de eletrodos (Fantasia, ECG-ID e MIT -BIH) e 75 % para o conjunto de dados fora da pessoa (CYBHi) e autenticação, atingindo taxas de erro iguais (EER) de quase 0 % para Fantasia e MIT-BIH e abaixo de 4 % para CYBHi. Quanto a (3) a abstração de sinais limpos e assimptomáticos de ECG e a detecção do seu desvio foram feitas e testadas em dois cenários diferentes: na presença de ruído usando um autocodificador e uma rede totalmente conectada (atingindo 99 % de precisão na classificação binária e 71 % na multi-classe), e; eventos de arritmia incluindo um RNN na arquitetura anterior (57 % de precisão e 61 % de sensibilidade). Em suma, esses sistemas são mais uma vez demonstrados como capazes de produzir resultados inovadores. A incorporação de vários sistemas de inteligência artificial em um unico sistema pederá desencadear a próxima geração de medicina preventiva. Os algoritmos ao terem acesso a diferentes estados fisiológicos e anatómicos, podem produzir soluções mais informadas para os problemas que se possam enfrentar no futuro, aumentando o desempenho de sistemas autónomos de prevenção que poderiam ser usados na vida quotidiana, nomeadamente em locais remotos onde o acesso à medicinas é limitado. Estes sistemas também ajudarão o estudo do comportamento do sinal e como eles são feitos no contexto da vida real, pois a IA explicável pode desencadear essa percepção e vincular os estados internos de uma rede às características biológicas

Brain connectivity mapping with diffusion MRI across individuals and species

The human brain is a highly complex organ that integrates functionally specialised subunits. Underpinning this complexity and functional specialisation is a network of structural connections, which may be probed using diffusion tractography, a unique, powerful and non-invasive MRI technique. Estimates of brain connectivity derived through diffusion tractography allow for explorations of how the brain’s functional subunits are inter-linked to subsequently produce experiences and behaviour. This thesis develops new diffusion tractography methodology for mapping brain connectivity, both across individuals and also across species; and explores frameworks for discovering associations of such brain connectivity features with behavioural traits. We build upon the hypothesis that connectional patterns can probe regions of functional equivalence across brains. To test this hypothesis we develop standardised and automated frameworks for mapping these patterns in very diverse brains, such as from human and non-human primates. We develop protocols to extract homologous fibre bundles across two species (human and macaque monkeys). We demonstrate robustness and generalisability of these protocols, but also their ability to capture individual variability. We also present investigations into how structural connectivity profiles may be used to inform us of how functionally-related features can be linked across different brains. Further, we explore how fully data-driven tractography techniques may be utilised for similar purposes, opening the door for future work on data-driven connectivity mapping. Subsequently, we explore how such individual variability in features that probe brain organisation are associated with differences in human behaviour. One approach to performing such explorations is the use of powerful multivariate statisitical techniques, such as canonical correlation analysis (CCA). After identifying issues in out-of-sample replication using multi-modal connectivity information, we perform comprehensive explorations into the robustness of such techniques and devise a generative model for forward predictions, demonstrating significant challlenges and limitations in their current applications. Specifically, we predict that the stability and generalisability of these techniques requires an order of magnitude more subjects than typically used to avoid overfitting and mis-interpretation of results. Using population-level data from the UK Biobank and confirmations from independent imaging modalities from the Human Connectome Project, we validate this prediction and demonstrate the direct link of CCA stability and generalisability with the number of subjects used per considered feature

Brain connectivity mapping with diffusion MRI across individuals and species

The human brain is a highly complex organ that integrates functionally specialised subunits. Underpinning this complexity and functional specialisation is a network of structural connections, which may be probed using diffusion tractography, a unique, powerful and non-invasive MRI technique. Estimates of brain connectivity derived through diffusion tractography allow for explorations of how the brain’s functional subunits are inter-linked to subsequently produce experiences and behaviour. This thesis develops new diffusion tractography methodology for mapping brain connectivity, both across individuals and also across species; and explores frameworks for discovering associations of such brain connectivity features with behavioural traits. We build upon the hypothesis that connectional patterns can probe regions of functional equivalence across brains. To test this hypothesis we develop standardised and automated frameworks for mapping these patterns in very diverse brains, such as from human and non-human primates. We develop protocols to extract homologous fibre bundles across two species (human and macaque monkeys). We demonstrate robustness and generalisability of these protocols, but also their ability to capture individual variability. We also present investigations into how structural connectivity profiles may be used to inform us of how functionally-related features can be linked across different brains. Further, we explore how fully data-driven tractography techniques may be utilised for similar purposes, opening the door for future work on data-driven connectivity mapping. Subsequently, we explore how such individual variability in features that probe brain organisation are associated with differences in human behaviour. One approach to performing such explorations is the use of powerful multivariate statisitical techniques, such as canonical correlation analysis (CCA). After identifying issues in out-of-sample replication using multi-modal connectivity information, we perform comprehensive explorations into the robustness of such techniques and devise a generative model for forward predictions, demonstrating significant challlenges and limitations in their current applications. Specifically, we predict that the stability and generalisability of these techniques requires an order of magnitude more subjects than typically used to avoid overfitting and mis-interpretation of results. Using population-level data from the UK Biobank and confirmations from independent imaging modalities from the Human Connectome Project, we validate this prediction and demonstrate the direct link of CCA stability and generalisability with the number of subjects used per considered feature