Axonal Velocity Distributions in Neural Field Equations

By modelling the average activity of large neuronal populations, continuum mean field models (MFMs) have become an increasingly important theoretical tool for understanding the emergent activity of cortical tissue. In order to be computationally tractable, long-range propagation of activity in MFMs is often approximated with partial differential equations (PDEs). However, PDE approximations in current use correspond to underlying axonal velocity distributions incompatible with experimental measurements. In order to rectify this deficiency, we here introduce novel propagation PDEs that give rise to smooth unimodal distributions of axonal conduction velocities. We also argue that velocities estimated from fibre diameters in slice and from latency measurements, respectively, relate quite differently to such distributions, a significant point for any phenomenological description. Our PDEs are then successfully fit to fibre diameter data from human corpus callosum and rat subcortical white matter. This allows for the first time to simulate long-range conduction in the mammalian brain with realistic, convenient PDEs. Furthermore, the obtained results suggest that the propagation of activity in rat and human differs significantly beyond mere scaling. The dynamical consequences of our new formulation are investigated in the context of a well known neural field model. On the basis of Turing instability analyses, we conclude that pattern formation is more easily initiated using our more realistic propagator. By increasing characteristic conduction velocities, a smooth transition can occur from self-sustaining bulk oscillations to travelling waves of various wavelengths, which may influence axonal growth during development. Our analytic results are also corroborated numerically using simulations on a large spatial grid. Thus we provide here a comprehensive analysis of empirically constrained activity propagation in the context of MFMs, which will allow more realistic studies of mammalian brain activity in the future.

[1]  Ingo Bojak,et al.  Self-organized 40Hz synchronization in a physiological theory of EEG , 2007, Neurocomputing.

[2]  Ben H. Jansen,et al.  Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns , 1995, Biological Cybernetics.

[3]  W. Freeman Mass action in the nervous system : examination of the neurophysiological basis of adaptive behavior through the EEG , 1975 .

[4]  P. Robinson,et al.  Mechanisms of cortical electrical activity and emergence of gamma rhythm. , 2000, Journal of theoretical biology.

[5]  Karl J. Friston,et al.  Dynamic causal modelling , 2003, NeuroImage.

[6]  D. Alistair Steyn-Ross,et al.  Modeling brain activation patterns for the default and cognitive states , 2009, NeuroImage.

[7]  S. E. Ahmed,et al.  Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference , 2008, Technometrics.

[8]  Peter A. Robinson,et al.  Unified neurophysical model of EEG spectra and evoked potentials , 2002, Biological Cybernetics.

[9]  Fabrice Wendling,et al.  Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals , 2000, Biological Cybernetics.

[10]  H. Haken,et al.  Field Theory of Electromagnetic Brain Activity. , 1996, Physical review letters.

[11]  J A Kelso,et al.  Spatiotemporal pattern formation in neural systems with heterogeneous connection topologies. , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[12]  Nicholas V Swindale Neural Synchrony, Axonal Path Lengths, and General Anesthesia: A Hypothesis , 2003, The Neuroscientist : a review journal bringing neurobiology, neurology and psychiatry.

[13]  P. Bressloff,et al.  Neural pattern formation in networks with dendritic structure , 1998 .

[14]  Stiliyan Kalitzin,et al.  Dynamical diseases of brain systems: different routes to epileptic seizures , 2003, IEEE Transactions on Biomedical Engineering.

[15]  J. J. Wright,et al.  EEG simulation: variation of spectral envelope, pulse synchrony and ≈40 Hz oscillation , 1997, Biological Cybernetics.

[16]  F. H. Lopes da Silva,et al.  Models of neuronal populations: the basic mechanisms of rhythmicity. , 1976, Progress in brain research.

[17]  D. Liley,et al.  Theoretical electroencephalogram stationary spectrum for a white-noise-driven cortex: evidence for a general anesthetic-induced phase transition. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[18]  O. Sporns,et al.  Mapping the Structural Core of Human Cerebral Cortex , 2008, PLoS biology.

[19]  Stephen G. Waxman,et al.  Ultrastructure of visual callosal axons in the rabbit , 1976, Experimental Neurology.

[20]  D. Liley,et al.  Drug-induced modification of the system properties associated with spontaneous human electroencephalographic activity. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  M. Young,et al.  Advanced database methodology for the Collation of Connectivity data on the Macaque brain (CoCoMac). , 2001, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[22]  Karl J. Friston,et al.  The Dynamic Brain: From Spiking Neurons to Neural Masses and Cortical Fields , 2008, PLoS Comput. Biol..

[23]  S. Waxman The Axon : structure, function, and pathophysiology , 1995 .

[24]  D. Liley,et al.  Understanding the Transition to Seizure by Modeling the Epileptiform Activity of General Anesthetic Agents , 2005, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[25]  Jean-Christophe Houzel,et al.  Visual inter-hemispheric processing: Constraints and potentialities set by axonal morphology , 1999, Journal of Physiology - Paris.

[26]  I. Tasaki,et al.  A study of conduction velocity in nonmyelinated nerve fibers. , 1977, Biophysical journal.

[27]  Fumitaka Kimura,et al.  Myelination and Isochronicity in Neural Networks , 2009, Front. Neuroanat..

[28]  P. Robinson,et al.  Dynamics of large-scale brain activity in normal arousal states and epileptic seizures. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[29]  J. Cowan,et al.  A mathematical theory of the functional dynamics of cortical and thalamic nervous tissue , 1973, Kybernetik.

[30]  B. Rosner,et al.  Inhalational anesthetics and conduction velocity of human peripheral nerve. , 1971, Electroencephalography and clinical neurophysiology.

[31]  Axel Hutt,et al.  Spontaneous and evoked activity in extended neural populations with gamma-distributed spatial interactions and transmission delay , 2007 .

[32]  J. Sleigh,et al.  The Sleep Cycle Modelled as a Cortical Phase Transition , 2005, Journal of biological physics.

[33]  M. Kramer,et al.  Pathological pattern formation and cortical propagation of epileptic seizures , 2005, Journal of The Royal Society Interface.

[34]  Karl J. Friston,et al.  Dynamic causal modeling of evoked responses in EEG and MEG , 2006, NeuroImage.

[35]  M. Bruce MacIver,et al.  Halothane depresses action potential conduction in hippocampal axons , 1998, Brain Research.

[36]  L. Kristiansson,et al.  Performance of a model for a local neuron population , 1978, Biological Cybernetics.

[37]  P. Rakić,et al.  Cytological and quantitative characteristics of four cerebral commissures in the rhesus monkey , 1990, The Journal of comparative neurology.

[38]  D. Liley,et al.  Modeling the effects of anesthesia on the electroencephalogram. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[39]  Axel Hutt,et al.  Effects of distributed transmission speeds on propagating activity in neural populations. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[40]  Mathew P. Dafilis,et al.  A spatially continuous mean field theory of electrocortical activity , 2002, Network.

[41]  P. Robinson Patchy propagators, brain dynamics, and the generation of spatially structured gamma oscillations. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  Rolf Kötter,et al.  Online retrieval, processing, and visualization of primate connectivity data from the CoCoMac Database , 2007, Neuroinformatics.

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

[44]  P. Szeszko,et al.  MRI atlas of human white matter , 2006 .

[45]  Egon Wanke,et al.  Mapping brains without coordinates , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[46]  H. Swadlow,et al.  Variations in conduction velocity and excitability following single and multiple impulses of visual callosal axons in the rabbit , 1976, Experimental Neurology.

[47]  D. Liley,et al.  Modeling electrocortical activity through improved local approximations of integral neural field equations. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  J. Bullier,et al.  Feedforward and feedback connections between areas V1 and V2 of the monkey have similar rapid conduction velocities. , 2001, Journal of neurophysiology.

[49]  I. A. Boyd,et al.  Scaling factor relating conduction velocity and diameter for myelinated afferent nerve fibres in the cat hind limb. , 1979, The Journal of physiology.

[50]  S. Wakana,et al.  MRI Atlas of Human White Matter , 2005 .

[51]  James J. Wright,et al.  Intracortical connectivity of pyramidal and stellate cells: estimates of synaptic densities and coupling symmetry , 1994 .

[52]  J. Ritchie,et al.  Physiology of axons , 1995 .

[53]  K. H. Lee,et al.  Correlation of cell body size, axon size, and signal conduction velocity for individually labelled dorsal root ganglion cells in the cat , 1986, The Journal of comparative neurology.

[54]  J. Cowan,et al.  Excitatory and inhibitory interactions in localized populations of model neurons. , 1972, Biophysical journal.

[55]  P. Robinson,et al.  A Quantitative Model of Sleep-Wake Dynamics Based on the Physiology of the Brainstem Ascending Arousal System , 2007, Journal of biological rhythms.

[56]  James J. Wright,et al.  Propagation and stability of waves of electrical activity in the cerebral cortex , 1997 .

[57]  Stephen G. Waxman,et al.  Unmyelinated and myelinated axon membrane from rat corpus callosum: differences in macromolecular structure , 1988, Brain Research.

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

[59]  M. Cross,et al.  Pattern formation outside of equilibrium , 1993 .

[60]  A. L. Muler,et al.  [Electrical properties of anisotropic neuromuscular syncytia. II. Distribution of a flat front of excitation]. , 1977, Biofizika.

[61]  A. M. Turing,et al.  The chemical basis of morphogenesis , 1952, Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences.

[62]  Karl J. Friston,et al.  Dynamic causal modelling of distributed electromagnetic responses , 2009, NeuroImage.

[63]  I. Tasaki,et al.  On the cable theory of nerve conduction , 2002, Bulletin of mathematical biology.

[64]  F. H. Lopes da Silva,et al.  Model of brain rhythmic activity , 1974, Kybernetik.

[65]  D. Oorschot,et al.  The number, size, and type of axons in rat subcortical white matter on left and right sides: A stereological, ultrastructural study , 2003, Journal of neurocytology.

[66]  E. Gordon,et al.  Stimulant drug action in attention deficit hyperactivity disorder (ADHD): inference of neurophysiological mechanisms via quantitative modelling , 2005, Clinical Neurophysiology.