Elsevier

Neuroscience Letters

Volume 658, 29 September 2017, Pages 97-101
Neuroscience Letters

Research article
What is the optimal distribution of myelin along a single axon?

https://doi.org/10.1016/j.neulet.2017.08.037Get rights and content

Highlights

  • Ion channel distribution determines the optimal pattern of partial myelination.

  • Short node-internode pairs yield fastest conduction if channel number is conserved.

  • Conduction velocity is insensitive to node-internode pair length if channel density is constant.

Abstract

The myelin sheath that insulates some axons in the central nervous system allows for faster signal conduction. Previously, axons were thought to be either unmyelinated or fully myelinated. Recent experimental work has discovered a new pattern of myelination (intermittent myelination) along axons in the mouse brain, in which long unmyelinated axon segments are followed by myelinated segments of comparable length. We use a computational model to explore how myelin distribution (in particular intermittent myelination) affects conduction velocity. We find that although fully myelinated axons minimize conduction velocity, varying the spatial distribution of a fixed amount of myelin along a partially myelinated axon leads to considerable variation in the conduction velocity for action potentials. Whether sodium ion channel number or sodium ion channel density is held constant as the area of the unmyelinated segments increases has a strong influence on the optimal pattern of myelin and the conduction velocity.

Introduction

Some axons in the vertebrate nervous system are wrapped with layers of myelin, which insulate these axons allowing for faster conduction of action potentials. The myelin sheath is produced by Schwann cells in the peripheral nervous system and oligodendrocytes in the central nervous system (CNS) [18], [19]. Unlike Schwann cells, which act on single peripheral axons, oligodendrocytes in the CNS ensheath up to 50 axons [7], allowing them to exert influence on neural processing on a larger scale. An individual oligodendrocyte, or a cluster of neighbouring oligodendrocytes, can have a large number of nearby axons available to myelinate, but experimental data [7] supports the proposition that the choice of which axons to myelinate is not made at random [27]. Moreover, myelination has been shown to be a dynamic process that responds to environmental cues [6], [14].

In a recent study, Tomassy et al. [25] analysed high-resolution maps of myelination by tracing high-throughput electron microscopy reconstructions of single axons of pyramidal neurons in the mouse brain. Analysing neurons in layers II/III from a publicly available dataset of a region of the mouse visual cortex [2], [3], they observed a new pattern of myelination. Historically, axons were thought to be either fully myelinated or unmyelinated. However, when tracing neurons in layers II/III, Tomassy et al. found that 17 out of 22 neurons displayed a pattern of myelination in which myelinated axon segments are interspersed with long unmyelinated segments, and they called this newly identified pattern “intermittent myelination” (IM). In this myelination pattern, unmyelinated sections of these axons were observed to be up to 55 μm long, much longer than typical nodes of Ranvier (approximately 1 μm long).

Since the discovery of intermittent myelination of (excitatory) pyramidal neurons by Tomassy et al. [25], Micheva et al. [15] also observed that the distribution of myelin in layer II/III inhibitory neurons was “patchy” with myelinated segments preferentially located near the cell body. Tomassy et al. [25] noted that neurons in layer II/III are involved in more complex cortical functioning than those found in layers V and VI, where IM was not observed. This raises the possibility that IM may be facilitating more complex neural functioning.

Whilst the evidence that signalling between electrically active axons and oligodendrocyte progenitor cells provides an important cue for inducing OPCs to differentiate into myelinating oligodendrocytes is firmly established, see [16], and the evidence that myelin provides more than just insulation to maximise conduction velocity (CV) continues to mount [9], the mechanisms controlling active myelination remain largely unknown (reviewed in Snaidero and Simons [23]).

The discovery of IM by Tomassy et al. raises several questions. What is the purpose of IM? Does an IM distribution provide any advantages to signal propagation over a fully myelinated pattern? Obviously, the shortest conduction time will be for a fully myelinated axon. However, it is not obvious which distribution of myelin will maximise CV when only a fixed fraction of the length of an axon is to be myelinated. In this instance, the distribution of ion channels also becomes important.

We examine the consequences for CV of the partial myelination of axons. The geometry of our model depends crucially on one parameter, L, which is the length of both the myelinated segments and the unmyelinated segments. This pair is repeated periodically over a fixed length of axon. We only study the effects of myelin distribution on signal CV but we note that there are several other facets of neural information processing that depend on myelin distribution, such as the reliability of action potential (AP) propagation [8] and the energy consumption [11], [20], [21].

Section snippets

Materials and methods

Our modelling was based on the NEURON [4], [5] implementation of a highly influential model of spike initiation in a myelinated mammalian axon by Mainen et al. [12]. We simplify this model by removing the dendrites but keep the ion channel kinetics characteristics, based on recordings of neocortical pyramidal neurons from the rat brain. Mainen et al. [12] used this kinetic model to simulate experimental measurements of rat pyramidal neurons [24].

Our simplification of this model consisted of a

Results

As noted earlier, little is known about the ion channel densities in the extended unmyelinated segments of IM patterned axons. Whilst the axon initial segment and hillock are known to have a high density of channels, the unmyelinated segments discovered by Tomassy et al. [25] extend much further from the soma, highlighting the need for new measurements before a definitive understanding of any biological advantage conferred by IM can be attained. However, our computational modeling enables the

Discussion

Our study reveals considerable differences in CV, with consequent implications on optimal allocation of a fixed total amount of myelin per neuron to maximise CV, depending on the density of sodium channels in the unmyelinated segments. If we assume that the number of sodium channels in the unmyelinated regions remains constant as the common length L of the myelinated and unmyelinated increases, we find that lower values of L lead to faster conduction. If instead, we assume that the density of

Acknowledgements

This work was supported by the Australian Research Council (DP140100339). DMW wishes to thank Dr. James Osborne for programming advice, Dr. Tania Kameneva for suggesting additional references and especially Dr. Tobias Merson for helpful comments.

References (27)

  • A.A. Faisal et al.

    Stochastic simulations on the reliability of action potential propagation in thin axons

    PLoS Comput. Biol.

    (2007)
  • R.D. Fields

    Myelin – more than insulation

    Science

    (2014)
  • S.A. Freeman et al.

    Mechanisms of sodium channel clustering and its influence on axonal impulse conduction

    Cell Mol. Life Sci.

    (2016)
  • Cited by (6)

    • A paradigm shift: Bioengineering meets mechanobiology towards overcoming remyelination failure

      2022, Biomaterials
      Citation Excerpt :

      More complex models are being rapidly developed to include full dendritic arbors with multiple cell types in a 3D environment [194]. Regarding myelination, in silico models were used to explore how OPC proliferation, migration differentiation and death affects the formation of organized arrays in the corpus callosum observed in vivo in early postnatal development [195], to determine how myelin distribution in axons affects nerve impulse conduction velocity [196] and even to assess how recruitment of immune cells to the CNS by chemotaxis induces OL death and eventual demyelination [197]. As computation power grows so will the complexity of the models used, being able to integrate exponentially more information and variables, working towards the ultimate goal of accurately simulating a living organism, from the cell to the whole organism.

    • Neural correlates of symptom severity in obsessive-compulsive disorder using magnetization transfer and diffusion tensor imaging

      2020, Psychiatry Research - Neuroimaging
      Citation Excerpt :

      The regions identified linked to symptom severity accord fairly well with regions previously implicated more generally in the pathophysiology of the disorder (Burguière et al., 2015; Chamberlain et al., 2005; Eng et al., 2015; Frydman et al., 2016; Gan et al., 2017; Graybiel and Rauch, 2000; Menzies et al., 2008a; Gonçalves et al., 2017a,b; Piras et al., 2015; van Velzen et al., 2014; Whiteside et al., 2004). The increased myelination related to OCD severity herein may relate to a decrease conduction velocity along WM axons, according to an “Inverted-U” function (Chomiak and Hu, 2009; Walsh et al., 2017; Wu et al., 2012). Functional MR studies have shown decreased functional connectivity in OCD cohort among numerous brain regions including OFC (Anticevic et al., 2014), cerebellum and occipital regions (Hou et al., 2014), posterior temporal areas (Hou et al., 2014; Zhang et al., 2011), hypoconnectivity between midbrain and ventral striatum as well as between lateral PFC and dorsal striatum (Harrison et al., 2009).

    View full text