Mapping axon conduction delays in vivo from microstructural MRI

The conduction velocity (CV) of action potentials along axons is a key neurophysiological property central to neural communication. The ability to estimate CV in humans in vivo from non-invasive MRI methods would therefore represent a significant advance in neuroscience. However, there are 2 major challenges that this paper aims to address: (1) much of the complexity of the neurophysiology of action potentials cannot be captured with currently available MRI techniques. Therefore, we seek to establish the variability in CV that can be captured when predicting CV purely from parameters that can be estimated from MRI (axon diameter and g-ratio); and (2) errors inherent in existing MRI-based biophysical models of tissue will propagate through to estimates of CV, the extent to which is currently unknown. Issue (1) is investigated by performing a sensitivity analysis on a comprehensive model of axon electrophysiology and determining the relative sensitivity to various morphological and electrical parameters. The investigations suggest that 89.2 % of the variance in CV is accounted for by variation in AD and g-ratio. The observed dependency of CV on AD and g-ratio is well characterised by a previously reported model by Rushton. Issue (2) is investigated through simulation of diffusion and relaxometry MRI data for a range of axon morphologies, applying models of restricted diffusion and relaxation processes to derive estimates of axon volume fraction (AVF), AD and g-ratio and estimating CV from the derived parameters. The results show that errors in the AVF have the biggest detrimental impact on estimates of CV, particularly for sparse fibre populations (AVF< 0.3). CV estimates are most accurate (below 5% error) where AVF is above 0.3, g-ratio is between 0.6 and 0.85 and AD is below 10 µm. Fortunately, these parameter bounds are typically satisfied by most myelinated axons. In conclusion, we demonstrate that accurate CV estimates can be inferred in axon populations across a range of configurations, except in some exceptional cases or where axonal density is low. As a proof of concept, for the first time, we generated an in vivo map of conduction velocity in the human corpus callosum with estimates consistent with values previously reported from invasive electrophysiology in primates.

[1]  R. S. Roche Is Polyoxymethylene a Cosmopolymer? , 1972, Nature.

[2]  Joseph A. Helpern,et al.  Characterizing microstructure of living tissues with time-dependent diffusion , 2012, 1210.3014.

[3]  D. Attwell,et al.  Node of Ranvier length as a potential regulator of myelinated axon conduction speed , 2017, eLife.

[4]  Julien Cohen-Adad,et al.  In vivo histology of the myelin g-ratio with magnetic resonance imaging , 2015, NeuroImage.

[5]  Zhengyi Yang,et al.  Towards higher sensitivity and stability of axon diameter estimation with diffusion‐weighted MRI , 2016, NMR in biomedicine.

[6]  Andrew L. Alexander,et al.  Mapping an index of the myelin g-ratio in infants using magnetic resonance imaging , 2016, NeuroImage.

[7]  Jelle Veraart,et al.  On the scaling behavior of water diffusion in human brain white matter , 2019, NeuroImage.

[8]  J. Young,et al.  The Rates Of Conduction Of Nerve Fibres Of Various Diameters In Cephalopods , 1938 .

[9]  Stephen M Smith,et al.  Fast robust automated brain extraction , 2002, Human brain mapping.

[10]  Carl-Fredrik Westin,et al.  Resolution limit of cylinder diameter estimation by diffusion MRI: The impact of gradient waveform and orientation dispersion , 2017, NMR in biomedicine.

[11]  Tobias C. Wood QUIT: QUantitative Imaging Tools , 2018, J. Open Source Softw..

[12]  S. Shinomoto,et al.  Can distributed delays perfectly stabilize dynamical networks? , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  S. Waxman Determinants of conduction velocity in myelinated nerve fibers , 1980, Muscle & nerve.

[14]  Kathryn L. West,et al.  Myelin volume fraction imaging with MRI , 2016, NeuroImage.

[15]  A. G. Richardson,et al.  Modelling the effects of electric fields on nerve fibres: Influence of the myelin sheath , 2000, Medical and Biological Engineering and Computing.

[16]  Tim B. Dyrby,et al.  Orientationally invariant indices of axon diameter and density from diffusion MRI , 2010, NeuroImage.

[17]  P. Sterling,et al.  Why Do Axons Differ in Caliber? , 2012, The Journal of Neuroscience.

[18]  V. Kiselev,et al.  Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation , 2016, NMR in biomedicine.

[19]  Carmen Birchmeier,et al.  Neuregulin-1/ErbB Signaling Serves Distinct Functions in Myelination of the Peripheral and Central Nervous System , 2008, Neuron.

[20]  W. Rushton A theory of the effects of fibre size in medullated nerve , 1951, The Journal of physiology.

[21]  H. S. Gasser,et al.  AXON DIAMETERS IN RELATION TO THE SPIKE DIMENSIONS AND THE CONDUCTION VELOCITY IN MAMMALIAN A FIBERS , 1939 .

[22]  S. Deoni,et al.  High‐resolution T1 mapping of the brain at 3T with driven equilibrium single pulse observation of T1 with high‐speed incorporation of RF field inhomogeneities (DESPOT1‐HIFI) , 2007, Journal of magnetic resonance imaging : JMRI.

[23]  Marco De Spirito,et al.  Modifications in solvent clusters embedded along the fibers of a cellulose polymer network cause paper degradation. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  James Kennedy,et al.  Particle swarm optimization , 2002, Proceedings of ICNN'95 - International Conference on Neural Networks.

[25]  J. Rosenbluth,et al.  Permeability of the Paranodal Junction of Myelinated Nerve Fibers , 2010, The Journal of Neuroscience.

[26]  Dmitry S. Novikov,et al.  What dominates the time dependence of diffusion transverse to axons: Intra- or extra-axonal water? , 2017, NeuroImage.

[27]  A. Huxley,et al.  Evidence for saltatory conduction in peripheral myelinated nerve fibres , 1949, The Journal of physiology.

[28]  Dmitry S. Novikov,et al.  Mesoscopic structure of neuronal tracts from time-dependent diffusion , 2015, NeuroImage.

[29]  G. Innocenti Network causality, axonal computations, and Poffenberger , 2017, Experimental Brain Research.

[30]  James L. Sinclair,et al.  Tuning of Ranvier node and internode properties in myelinated axons to adjust action potential timing , 2015, Nature Communications.

[31]  A. Schüz,et al.  Basic Connectivity of the Cerebral Cortex and some Considerations on the Corpus Callosum , 1996, Neuroscience & Biobehavioral Reviews.

[32]  David S. Lorberbaum,et al.  Genetic evidence that Nkx2.2 acts primarily downstream of Neurog3 in pancreatic endocrine lineage development , 2017, eLife.

[33]  J. Hursh CONDUCTION VELOCITY AND DIAMETER OF NERVE FIBERS , 1939 .

[34]  R. Llinás,et al.  1516 On uniform olivocerebellar conduction time that underlies purkinje cell complex spike synchronicity in the rat cerebellum , 1993 .

[35]  Daniel C. Alexander,et al.  NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain , 2012, NeuroImage.

[36]  Daniel R. Berger,et al.  Distinct Profiles of Myelin Distribution Along Single Axons of Pyramidal Neurons in the Neocortex , 2014, Science.

[37]  Jelle Veraart,et al.  In vivo observation and biophysical interpretation of time-dependent diffusion in human white matter , 2016, NeuroImage.

[38]  H. Swadlow,et al.  Characteristics of interhemispheric impulse conduction between prelunate gyri of the rhesus monkey , 1978, Experimental Brain Research.

[39]  S. Deoni,et al.  Transverse relaxation time (T2) mapping in the brain with off‐resonance correction using phase‐cycled steady‐state free precession imaging , 2009, Journal of magnetic resonance imaging : JMRI.

[40]  Deepika Agrawal,et al.  Internodal myelination during development quantitated using X-ray diffraction. , 2009, Journal of structural biology.

[41]  P. Basser,et al.  Axcaliber: A method for measuring axon diameter distribution from diffusion MRI , 2008, Magnetic resonance in medicine.

[42]  R. Caminiti,et al.  The diameter of cortical axons depends both on the area of origin and target. , 2014, Cerebral cortex.

[43]  Derek K. Jones,et al.  Gleaning multicomponent T1 and T2 information from steady‐state imaging data , 2008, Magnetic resonance in medicine.

[44]  M. Bennett,et al.  Relative conduction velocities of small myelinated and non-myelinated fibres in the central nervous system. , 1972, Nature: New biology.

[45]  M. Konishi,et al.  Axonal delay lines for time measurement in the owl's brainstem. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[46]  Péter Buzás,et al.  Neocortical Axon Arbors Trade-off Material and Conduction Delay Conservation , 2010, PLoS Comput. Biol..

[47]  R. Llinás,et al.  Uniform olivocerebellar conduction time underlies Purkinje cell complex spike synchronicity in the rat cerebellum. , 1993, The Journal of physiology.

[48]  Zhengyi Yang,et al.  Parametric Probability Distribution Functions for Axon Diameters of Corpus Callosum , 2016, Front. Neuroanat..

[49]  J W Moore,et al.  Simulations of conduction in uniform myelinated fibers. Relative sensitivity to changes in nodal and internodal parameters. , 1978, Biophysical journal.

[50]  Derek K. Jones,et al.  Precision and Accuracy in Diffusion Tensor Magnetic Resonance Imaging , 2010, Topics in magnetic resonance imaging : TMRI.

[51]  G. Yovel,et al.  In vivo correlation between axon diameter and conduction velocity in the human brain , 2014, Brain Structure and Function.

[52]  Hui Zhang,et al.  PGSE, OGSE, and sensitivity to axon diameter in diffusion MRI: Insight from a simulation study , 2015, Magnetic resonance in medicine.

[53]  R. Caminiti,et al.  Diameter, Length, Speed, and Conduction Delay of Callosal Axons in Macaque Monkeys and Humans: Comparing Data from Histology and Magnetic Resonance Imaging Diffusion Tractography , 2013, The Journal of Neuroscience.

[54]  E. Formisano,et al.  Reproducibility and Reliability of Quantitative and Weighted T1 and T2∗ Mapping for Myelin-Based Cortical Parcellation at 7 Tesla , 2016, Front. Neuroanat..

[55]  Stamatios N. Sotiropoulos,et al.  An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging , 2016, NeuroImage.

[56]  R. S. Smith,et al.  Myelinated nerve fibers: computed effect of myelin thickness on conduction velocity. , 1970, The American journal of physiology.

[57]  Yaniv Assaf,et al.  Composite hindered and restricted model of diffusion (CHARMED) MR imaging of the human brain , 2005, NeuroImage.