Name: how to measure cortical folding from mr images: a step-by-step tutorial to compute local gyrification index
Uploaded: 2012-01-02T13:00:00
Description: Watch this Scientific Journal Video about How to Measure Cortical Folding from MR Images: a Step-by-Step Tutorial to Compute Local Gyrification Index at JoVE.com

This research was supported by the National Center of Competence in Research (NCCR) "SYNAPSY - The Synaptic Bases of Mental Diseases" financed by the Swiss National Science Foundation (n&deg; 51AU40_125759). Development of the local Gyrification Index was supported by grants from the Swiss National Research Fund to Dr. Marie Schaer (323500-111165) and to Dr. Stephan Eliez (3200-063135.00/1, 3232-063134.00/1, PP0033-102864 and 32473B-121996) and by the Center for Biomedical Imaging (CIBM) of the Geneva-Lausanne Universities and the EPFL, as well as the foundations Leenaards and Louis-Jeantet. Support for the development of FreeSurfer software was provided in part by the National Center for Research Resources (P41-RR14075, and the NCRR BIRN Morphometric Project BIRN002, U24 RR021382), the National Institute for Biomedical Imaging and Bioengineering (R01 EB001550, R01EB006758), the National Institute for Neurological Disorders and Stroke (R01 NS052585-01) as well as the Mental Illness and Neuroscience Discovery (MIND) Institute, and is part of the National Alliance for Medical Image Computing (NAMIC), funded by the National Institutes of Health through the NIH Roadmap for Medical Research, Grant U54 EB005149. Additional support was provided by The Autism &amp; Dyslexia Project funded by the Ellison Medical Foundation.

1. Reconstruct the 3D cortical surfaces
This first part of the protocol uses the standard FreeSurfer pipeline as described in the Wiki (<a href="http://surfer.nmr.mgh.harvard.edu/fswiki">http://surfer.nmr.mgh.harvard.edu/fswiki</a>). Note that the commands detailed here describe one way of achieving the cortical surface reconstructions, but equivalent commands may also be used.
<ol class="lalpha">
 <li>Import the raw MRI DICOM into FreeSurfer and verify the quality of the i...

<ol>
	<li>Schaer, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+surface-based+approach+to+quantify+local+cortical+gyrification.">A surface-based approach to quantify local cortical gyrification.</a> IEEE. Trans. Med. Imaging. 27, 161-170 (2008).</li><li>Mangin, J. F., Jouvent, E., Cachia, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In-vivo+measurement+of+cortical+morphology:+means+and+meanings.">In-vivo measurement of cortical morphology: means and meanings.</a> Curr. Opin. Neurol. 23, 359-367 (2010).</li><li>Mangin, J. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+framework+to+study+the+cortical+folding+patterns.">A framework to study the cortical folding patterns.</a> Neuroimage. 23, S129-S138 (2004).</li><li>Luders, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+curvature-based+approach+to+estimate+local+gyrification+on+the+cortical+surface.">A curvature-based approach to estimate local gyrification on the cortical surface.</a> Neuroimage. 29, 1224-1230 (2006).</li><li>Zilles, K., Armstrong, E., Schleicher, A., Kretschmann, H. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+human+pattern+of+gyrification+in+the+cerebral+cortex.">The human pattern of gyrification in the cerebral cortex.</a> Anat. Embryol. (Berl). 179, 173-179 (1988).</li><li>Rakic, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Specification+of+cerebral+cortical+areas.">Specification of cerebral cortical areas.</a> Science. 241, 170-176 (1988).</li><li>Toro, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Brain+size+and+folding+of+the+human+cerebral+cortex.">Brain size and folding of the human cerebral cortex.</a> Cereb. Cortex. 18, 2352-2357 (2008).</li><li>Fischl, B., Sereno, M. I., Dale, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Cortical+surface-based+analysis.+II:+Inflation,+flattening,+and+a+surface-based+coordinate+system.">Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system.</a> Neuroimage. 9, 195-207 (1999).</li><li>Dale, A. M., Fischl, B., Sereno, M. I. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Cortical+surface-based+analysis.+I.+Segmentation+and+surface+reconstruction.">Cortical surface-based analysis. I. Segmentation and surface reconstruction.</a> Neuroimage. 9, 179-194 (1999).</li><li>Genovese, C. R., Lazar, N. A., Nichols, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Thresholding+of+statistical+maps+in+functional+neuroimaging+using+the+false+discovery+rate.">Thresholding of statistical maps in functional neuroimaging using the false discovery rate.</a> Neuroimage. 15, 870-878 (2002).</li><li>Desikan, R. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=An+automated+labeling+system+for+subdividing+the+human+cerebral+cortex+on+MRI+scans+into+gyral+based+regions+of+interest.">An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest.</a> Neuroimage. 31, 968-980 (2006).</li><li>Schaer, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Congenital+heart+disease+affects+local+gyrification+in+22q11.2+deletion+syndrome.">Congenital heart disease affects local gyrification in 22q11.2 deletion syndrome.</a> Dev. Med. Child. Neurol. 51, 746-753 (2009).</li><li>Palaniyappan, L., Mallikarjun, P., Joseph, V., White, T. P., Liddle, P. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Folding+of+the+Prefrontal+Cortex+in+Schizophrenia:+Regional+Differences+in+Gyrification.">Folding of the Prefrontal Cortex in Schizophrenia: Regional Differences in Gyrification.</a> Biol. Psychiatry. , (2011).</li><li>Zhang, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Decreased+gyrification+in+major+depressive+disorder.">Decreased gyrification in major depressive disorder.</a> Neuroreport. 20, 378-380 (2009).</li><li>Juranek, J., Salman, M. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Anomalous+development+of+brain+structure+and+function+in+spina+bifida+myelomeningocele.">Anomalous development of brain structure and function in spina bifida myelomeningocele.</a> Dev. Disabil. Res. Rev. 16, 23-30 (2010).</li><li>Zhang, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Reduced+cortical+folding+in+mental+retardation.">Reduced cortical folding in mental retardation.</a> AJNR. Am. J. Neuroradiol. 31, 1063-1067 (2010).</li><li>Kuperberg, G. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Regionally+localized+thinning+of+the+cerebral+cortex+in+schizophrenia.">Regionally localized thinning of the cerebral cortex in schizophrenia.</a> Archives of general psychiatry. 60, 878-888 (2003).</li><li>Milad, M. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Thickness+of+ventromedial+prefrontal+cortex+in+humans+is+correlated+with+extinction+memory.">Thickness of ventromedial prefrontal cortex in humans is correlated with extinction memory.</a> Proc. Natl. Acad. Sci. U. S. A. 102, 10706-10711 (2005).</li><li>Rauch, S. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+magnetic+resonance+imaging+study+of+cortical+thickness+in+animal+phobia.">A magnetic resonance imaging study of cortical thickness in animal phobia.</a> Biol. Psychiatry. 55, 946-952 (2004).</li><li>Fjell, A. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Selective+increase+of+cortical+thickness+in+high-performing+elderly--structural+indices+of+optimal+cognitive+aging.">Selective increase of cortical thickness in high-performing elderly--structural indices of optimal cognitive aging.</a> Neuroimage. 29, 984-994 (2006).</li><li>Walhovd, K. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Regional+cortical+thickness+matters+in+recall+after+months+more+than+minutes.">Regional cortical thickness matters in recall after months more than minutes.</a> Neuroimage. 31, 1343-1351 (2006).</li><li>Gold, B. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Differing+neuropsychological+and+neuroanatomical+correlates+of+abnormal+reading+in+early-stage+semantic+dementia+and+dementia+of+the+Alzheimer+type.">Differing neuropsychological and neuroanatomical correlates of abnormal reading in early-stage semantic dementia and dementia of the Alzheimer type.</a> Neuropsychologia. 43, 833-846 (2005).</li><li>Salat, D. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Thinning+of+the+cerebral+cortex+in+aging.">Thinning of the cerebral cortex in aging.</a> Cereb. Cortex. 14, 721-730 (2004).</li><li>Schaer, M., Eliez, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Contribution+of+structural+brain+imaging+to+our+understanding+of+cortical+development+process.">Contribution of structural brain imaging to our understanding of cortical development process.</a> European Psychiatry Reviews. 2, 13-16 (2009).</li><li>Shaw, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Neurodevelopmental+trajectories+of+the+human+cerebral+cortex.">Neurodevelopmental trajectories of the human cerebral cortex.</a> J. Neurosci. 28, 3586-3594 (2008).</li></ol>

The protocol above describes how to measure local Gyrification Index based on cerebral T1-weighted MRI and conduct statistical group comparisons. Our method has been specifically designed to localize early disruption in the cortical expansion process and as such is of particular interest in many neurodevelopmental or psychiatric conditions. Examples of group comparisons in clinical samples can be found in publications by our group1,12 or by others13-16. The process is fully automated and requ...

Material: a Unix or Mac workstation with a processor of 2GHz or faster and a minimum of 4GB of RAM, with FreeSurfer installed (<a href="http://surfer.nmr.mgh.harvard.edu/fswiki">http://surfer.nmr.mgh.harvard.edu/fswiki</a>, preferably the latest version, but no older than version 4.0.3). In order to compute the local Gyrification Index, MATLAB is also required (<a href="http://www.mathworks.ch/">http://www.mathworks.com/</a>) along with the Image Processing Toolbox. 
Data: A sample of good quality (high-resolution, high contrast) cerebral MRI T1-weighted dataset. Your group of subjects must be preferably matched for age and gender. Given the normal inter-individual variability in cerebral morphology, the number of subjects in each group should be sufficient to identify an existing group difference (the more - the better). A reasonable minimum sample size would be around 20 subjects per group (although you can probably go for less if the intensity of changes is large and if your groups are tightly matched for gender and age).
<table border="1">
<tbody>
 <tr>
 <td>Name of the equipment</td>
 <td>Company</td>
 <td>Catalogue number</td>
 <td>Comments</td>
 </tr>
 <tr>
 <td>FreeSurfer </td>
 <td>Martinos Center for Biomedical Imaging, MGH</td>
 <td>&nbsp;</td>
 <td>Version newer than 4.0.3</td>
 </tr>
 <tr>
 <td>Matlab</td>
 <td>Mathworks</td>
 <td>&nbsp;</td>
 <td>Image Processing Toolbox</td>
 </tr>
 </tbody>
</table>

how to measure cortical folding from mr images: a step-by-step tutorial to compute local gyrification index

Cortical folding (gyrification) is determined during the first months of life, so that adverse events occurring during this period leave traces that will be identifiable at any age. As recently reviewed by Mangin and colleagues2, several methods exist to quantify different characteristics of gyrification. For instance, sulcal morphometry can be used to measure shape descriptors such as the depth, length or indices of inter-hemispheric asymmetry3. These geometrical properties have the advantage of being easy to interpret. However, sulcal morphometry tightly relies on the accurate identification of a given set of sulci and hence provides a fragmented description of gyrification. A more fine-grained quantification of gyrification can be achieved with curvature-based measurements, where smoothed absolute mean curvature is typically computed at thousands of points over the cortical surface4. The curvature is however not straightforward to comprehend, as it remains unclear if there is any direct relationship between the curvedness and a biologically meaningful correlate such as cortical volume or surface. To address the diverse issues raised by the measurement of cortical folding, we previously developed an algorithm to quantify local gyrification with an exquisite spatial resolution and of simple interpretation. Our method is inspired of the Gyrification Index5, a method originally used in comparative neuroanatomy to evaluate the cortical folding differences across species. In our implementation, which we name local Gyrification Index (lGI1), we measure the amount of cortex buried within the sulcal folds as compared with the amount of visible cortex in circular regions of interest. Given that the cortex grows primarily through radial expansion6, our method was specifically designed to identify early defects of cortical development.
In this article, we detail the computation of local Gyrification Index, which is now freely distributed as a part of the FreeSurfer Software (<a href="http://surfer.nmr.mgh.harvard.edu/">http://surfer.nmr.mgh.harvard.edu/</a>, Martinos Center for Biomedical Imaging, Massachusetts General Hospital). FreeSurfer provides a set of automated reconstruction tools of the brain's cortical surface from structural MRI data. The cortical surface extracted in the native space of the images with sub-millimeter accuracy is then further used for the creation of an outer surface, which will serve as a basis for the lGI calculation. A circular region of interest is then delineated on the outer surface, and its corresponding region of interest on the cortical surface is identified using a matching algorithm as described in our validation study1. This process is repeatedly iterated with largely overlapping regions of interest, resulting in cortical maps of gyrification for subsequent statistical comparisons (Fig. 1). Of note, another measurement of local gyrification with a similar inspiration was proposed by Toro and colleagues7, where the folding index at each point is computed as the ratio of the cortical area contained in a sphere divided by the area of a disc with the same radius. The two implementations differ in that the one by Toro et al. is based on Euclidian distances and thus considers discontinuous patches of cortical area, whereas ours uses a strict geodesic algorithm and include only the continuous patch of cortical area opening at the brain surface in a circular region of interest.

Watch this Scientific Journal Video about How to Measure Cortical Folding from MR Images: a Step-by-Step Tutorial to Compute Local Gyrification Index at JoVE.com

How to Measure Cortical Folding from MR Images: a Step-by-Step Tutorial to Compute Local Gyrification Index

Measuring gyrification (cortical folding) at any age represents a window into early brain development. Hence, we previously developed an algorithm to measure local gyrification at thousands of points over the hemisphere1. In this paper, we detail the computation of this local gyrification index.

Cortical folding (gyrification) is determined during the first months of life, so that adverse events occurring during this period leave traces that will ...

Research

JoVE Journal

Medicine

Intrauterine Telemetry to Measure Mouse Contractile Pressure In Vivo

Intrauterine Telemetry to Measure Mouse Contractile Pressure In Vivo

A complex integration of molecular and electrical signals is needed to transform a quiescent uterus into a contractile organ at the end of pregnancy. ...

Step By Step: Microsurgical training method combining two nonliving animal models

The learning of microsurgical techniques and the maintenance of microsurgical skills have been traditionally based on the use of living animals, mainly ...

Assessing Cortical Cerebral Microinfarcts on High Resolution MR Images

Cerebral microinfarcts are frequent findings in the post-mortem human brain, and are related to cognitive decline and dementia. Due to their small sizes ...

Heterotopic Renal Autotransplantation in a Porcine Model: A Step-by-Step Protocol

Kidney transplantation is the treatment of choice for patients suffering from end-stage renal disease. It offers better life expectancy and higher quality ...

A Step by Step Protocol for Subretinal Surgery in Rabbits

Age related macular degeneration (AMD), retinitis pigmentosa, and other RPE related diseases are the most common causes for irreversible loss of vision in ...

A Clinical Trial Assessing the Safety, Efficacy, and Delivery of Olive-Oil-Based Three-Chamber Bags for Parenteral Nutrition

Limited evidence exists to precisely estimate efficacy and safety differences between parenteral nutrition (PN) prepared using olive-oil-based ...

Conducting Maximal and Submaximal Endurance Exercise Testing to Measure Physiological and Biological Responses to Acute Exercise in Humans

Regular physical activity has a positive effect on human health, but the mechanisms controlling these effects remain unclear. The physiologic and biologic ...

Obtaining Quality Extended Field-of-View Ultrasound Images of Skeletal Muscle to Measure Muscle Fascicle Length

Muscle fascicle length, which is commonly measured in vivo using traditional ultrasound, is an important parameter defining a muscle&#8217;s force ...

Step-by-Step Stapedotomy through Transcanal Exclusive Endoscopic Approach

In recent years there has been an increasing trend in the use of the endoscope to treat a variety of middle ear pathologies, including otosclerosis. ...

Improved Renal Denervation Mitigated Hypertension Induced by Angiotensin II Infusion

The benefits of renal sympathetic denervation (RDN) on blood pressure have been proved in a large number of clinical trials in recent years. However, the ...