Estimating axon conduction velocity 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 two 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 have been reported to be estimatable from MRI: inner axon diameter (AD) and g-ratio. (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 85% 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 the 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). For our equipment set-up and acquisition protocol, 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 high (above 4μm). CV estimates are robust to errors in g-ratio estimation but are highly sensitive to errors in AD estimation, particularly where ADs are small. We additionally show CV estimates in human corpus callosum in a small number of subjects. In conclusion, we demonstrate accurate CV estimates are possible in regions of the brain where AD is sufficiently large. Problems with estimating ADs for smaller axons presents a problem for estimating CV across the whole CNS and should be the focus of further study.

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

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

[3]  G. Schwarz Estimating the Dimension of a Model , 1978 .

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

[5]  Sebastien Ourselin,et al.  The importance of correcting for signal drift in diffusion MRI , 2017, Magnetic resonance in medicine.

[6]  Siawoosh Mohammadi,et al.  Four in vivo g‐ratio‐weighted imaging methods: Comparability and repeatability at the group level , 2018, Human brain mapping.

[7]  A. Scheibel,et al.  Fiber composition of the human corpus callosum , 1992, Brain Research.

[8]  D J Tolhurst,et al.  Effect of myelination on the conduction velocity of optic nerve fibres , 1992, Ophthalmic & physiological optics : the journal of the British College of Ophthalmic Opticians.

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

[10]  D. Barazany,et al.  AxCaliber 3D , 2010 .

[11]  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.

[12]  R. Friede,et al.  The precise geometry of large internodes , 1980, Journal of the Neurological Sciences.

[13]  Guy Eyal,et al.  Unique membrane properties and enhanced signal processing in human neocortical neurons , 2016, eLife.

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

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

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

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

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

[19]  P. V. van Zijl,et al.  Evaluation of restricted diffusion in cylinders. Phosphocreatine in rabbit leg muscle. , 1994, Journal of magnetic resonance. Series B.

[20]  Jean-Philippe Thiran,et al.  COMMIT: Convex Optimization Modeling for Microstructure Informed Tractography , 2015, IEEE Transactions on Medical Imaging.

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

[22]  J. Young,et al.  Internode length and fibre diameter in developing and regenerating nerves. , 1948, Journal of anatomy.

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

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

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

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

[27]  Rainer Goebel,et al.  Robust and fast nonlinear optimization of diffusion MRI microstructure models , 2017, NeuroImage.

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

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

[30]  Julien Cohen-Adad,et al.  Promise and pitfalls of g-ratio estimation with MRI , 2017, NeuroImage.

[31]  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.

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

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

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

[35]  M. Neal,et al.  Autoradiographic localization of 3 H-GABA in rat retina. , 1972, Nature: New biology.

[36]  Richard Bowtell,et al.  Microstructural imaging of the human brain with a ‘super-scanner’: 10 key advantages of ultra-strong gradients for diffusion MRI , 2018, NeuroImage.

[37]  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.

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

[39]  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.

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

[41]  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.

[42]  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.

[43]  M. Berry,et al.  Relationship between myelin sheath diameter and internodal length in axons of the anterior medullary velum of the adult rat , 1995, Journal of the Neurological Sciences.

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

[45]  R. Yarmukhamedov,et al.  Comments to the paper , 2010 .

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

[47]  M. J. D. Powell,et al.  An efficient method for finding the minimum of a function of several variables without calculating derivatives , 1964, Comput. J..

[48]  K. Smith,et al.  Internodal myelin volume and axon surface area A relationship determining myelin thickness? , 1982, Journal of the Neurological Sciences.

[49]  S. Berman,et al.  Modeling conduction delays in the corpus callosum using MRI-measured g-ratio , 2019, NeuroImage.

[50]  R. Caminiti,et al.  Comments on the paper by Horowitz et al. (2014) , 2014, Brain Structure and Function.

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

[52]  G. B. Pike,et al.  Quantitative imaging of magnetization transfer exchange and relaxation properties in vivo using MRI , 2001, Magnetic resonance in medicine.

[53]  Joseph V. Hajnal,et al.  Inherent and unpredictable bias in multi-component DESPOT myelin water fraction estimation , 2019, NeuroImage.

[54]  W. Blakemore,et al.  The relationship between internodal length and fibre diameter in the spinal cord of the cat , 1980, Journal of the Neurological Sciences.

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

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

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

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

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

[60]  A. Vizoso The relationship between internodal length and growth in human nerves. , 1950, Journal of anatomy.

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

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

[63]  Derek K. Jones,et al.  Including diffusion time dependence in the extra-axonal space improves in vivo estimates of axonal diameter and density in human white matter , 2016, NeuroImage.

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

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

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

[67]  Els Fieremans,et al.  Revealing mesoscopic structural universality with diffusion , 2014, Proceedings of the National Academy of Sciences.

[68]  P. Batchelor,et al.  International Society for Magnetic Resonance in Medicine , 1997 .

[69]  Hui Zhang,et al.  Whole brain g-ratio mapping using myelin water imaging (MWI) and neurite orientation dispersion and density imaging (NODDI) , 2017, NeuroImage.

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

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

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

[73]  Mark Jenkinson,et al.  The minimal preprocessing pipelines for the Human Connectome Project , 2013, NeuroImage.

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

[75]  Christopher L Lankford,et al.  On the inherent precision of mcDESPOT , 2013, Magnetic resonance in medicine.

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

[77]  Bibek Dhital,et al.  Gibbs‐ringing artifact removal based on local subvoxel‐shifts , 2015, Magnetic resonance in medicine.

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

[79]  Nikos K. Logothetis,et al.  Distribution of axon diameters in cortical white matter: an electron-microscopic study on three human brains and a macaque , 2014, Biological Cybernetics.

[80]  S. Nelson,et al.  Distinct Profiles of Myelin Distribution Along Single Axons of Pyramidal Neurons in the Neocortex , 2014 .

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

[82]  Julien Cohen-Adad,et al.  g-Ratio weighted imaging of the human spinal cord in vivo , 2017, NeuroImage.

[83]  R. Caminiti,et al.  Areal differences in diameter and length of corticofugal projections. , 2012, Cerebral Cortex.

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

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