A tutorial and tool for exploring feature similarity gradients with MRI data
Claude J. Bajada,
Lucas Q. Costa Campos,
Svenja Caspers,
Richard Muscat,
Geoff J.M. Parker,
Matthew A. Lambon Ralph,
Lauren L. Cloutman,
Nelson J. Trujillo-Barreto
Affiliations
Claude J. Bajada
Department of Physiology and Biochemistry, Faculty of Medicine and Surgery, The University of Malta, Malta; Division of Neuroscience & Experimental Psychology, School of Biological Sciences, The University of Manchester, UK; Institute of Neuroscience and Medicine (INM-1), Research Centre Jülich, 52425, Jülich, Germany; Corresponding author. Department of Physiology and Biochemistry, Faculty of Medicine and Surgery, The University of Malta, Msida, Malta
Lucas Q. Costa Campos
Institute of Neuroscience and Medicine (INM-1), Research Centre Jülich, 52425, Jülich, Germany; Institute of Complex Systems and Institute for Advanced Simulation (ICS-2/IAS-2), Research Centre Jülich, 52425, Jülich, Germany
Svenja Caspers
Institute of Neuroscience and Medicine (INM-1), Research Centre Jülich, 52425, Jülich, Germany; Institute for Anatomy I, Medical Faculty, Heinrich-Heine-University Duesseldorf, 40221, Duesseldorf, Germany; JARA-BRAIN, Jülich-Aachen Research Alliance, 52425, Jülich, Germany
Richard Muscat
Department of Physiology and Biochemistry, Faculty of Medicine and Surgery, The University of Malta, Malta
Geoff J.M. Parker
Centre for Medical Image Computing, Department of Computer Science, and Queen Square MS Centre, Department of Neuroinflammation, UCL Institute of Neurology, University College London, UK; Bioxydyn Limited, Manchester, UK
Matthew A. Lambon Ralph
MRC Cognition and Brain Sciences Unit, University of Cambridge, Cambridge, UK
Lauren L. Cloutman
Division of Neuroscience & Experimental Psychology, School of Biological Sciences, The University of Manchester, UK
Nelson J. Trujillo-Barreto
Division of Neuroscience & Experimental Psychology, School of Biological Sciences, The University of Manchester, UK; Corresponding author. Division of Neuroscience and Experimental Psychology, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Academic Health Science Centre, Brunswick Street, Manchester, M13 9PL, UK
There has been an increasing interest in examining organisational principles of the cerebral cortex (and subcortical regions) using different MRI features such as structural or functional connectivity. Despite the widespread interest, introductory tutorials on the underlying technique targeted for the novice neuroimager are sparse in the literature.Articles that investigate various “neural gradients” (for example based on region studied “cortical gradients,” “cerebellar gradients,” “hippocampal gradients” etc … or feature of interest “functional gradients,” “cytoarchitectural gradients,” “myeloarchitectural gradients” etc …) have increased in popularity. Thus, we believe that it is opportune to discuss what is generally meant by “gradient analysis”. We introduce basics concepts in graph theory, such as graphs themselves, the degree matrix, and the adjacency matrix. We discuss how one can think about gradients of feature similarity (the similarity between timeseries in fMRI, or streamline in tractography) using graph theory and we extend this to explore such gradients across the whole MRI scale; from the voxel level to the whole brain level. We proceed to introduce a measure for quantifying the level of similarity in regions of interest. We propose the term “the Vogt-Bailey index” for such quantification to pay homage to our history as a brain mapping community.We run through the techniques on sample datasets including a brain MRI as an example of the application of the techniques on real data and we provide several appendices that expand upon details. To maximise intuition, the appendices contain a didactic example describing how one could use these techniques to solve a particularly pernicious problem that one may encounter at a wedding. Accompanying the article is a tool, available in both MATLAB and Python, that enables readers to perform the analysis described in this article on their own data.We refer readers to the graphical abstract as an overview of the analysis pipeline presented in this work.