Hypothesis and Theory ARTICLE
Proposed evolutionary changes in the role of myelin
- 1The MARCS Institute, Sydney, University of Western Sydney, NSW, Australia
- 2The Blue Brain Project, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
- 3Department of Neurobiology, Hebrew University of Jerusalem, Edmond J. Safra Campus, Jerusalem, Israel
- 4NeuroLinx Research Institute, La Jolla, CA, USA
Myelin is the multi-layered lipid sheet periodically wrapped around neuronal axons. It is most frequently found in vertebrates. Myelin allows for saltatory action potential (AP) conduction along axons. During this form of conduction, the AP travels passively along the myelin-covered part of the axon, and is recharged at the intermittent nodes of Ranvier. Thus, myelin can reduce the energy load needed and/or increase the speed of AP conduction. Myelin first evolved during the Ordovician period. We hypothesize that myelin's first role was mainly energy conservation. During the later “Mesozoic marine revolution,” marine ecosystems changed toward an increase in marine predation pressure. We hypothesize that the main purpose of myelin changed from energy conservation to conduction speed increase during this Mesozoic marine revolution. To test this hypothesis, we optimized models of myelinated axons for a combination of AP conduction velocity and energy efficiency. We demonstrate that there is a trade-off between these objectives. We then compared the simulation results to empirical data and conclude that while the data are consistent with the theory, additional measurements are necessary for a complete evaluation of the proposed hypothesis.
Myelin is a neural cellular specialization unique to vertebrates (with a few interesting exceptions, Davis et al., 1999; Hartline and Colman, 2007). It is derived from glial cells, and is composed of lipid sheaths of cellular membranes wrapped around the axon of a neuron. Myelin improves the conduction of action potentials (APs), large (~100 mV), brief (<2 ms) depolarizations of the neural membrane, which are used to convey information to subsequent downstream neurons. AP conduction is improved by increasing the insulation of the axonal membrane, thereby elevating membrane resistance and decreasing capacitance. This allows the AP to travel along the axon with less loss of current through the membrane. APs propagate passively through myelinated portions of the axons, and are regenerated by active, voltage dependent Na+ and K+ ion-specific conductances at the intermittent, non-myelinated, nodes of Ranvier. This process is called saltatory conduction (Johnston and Wu, 1994; Figure 1).
Figure 1. Myelin sheets engulf axons and can increase conduction speed and energy efficiency of action potential (AP) conduction. (A) An electro-micrograph of a myelin sheath. Image by the Electron Microscopy Facility at Trinity College via Wikipedia. (B) Schematic representation of salutatory AP conduction. Different parameters are optimal for conduction speed (bottom) and energy efficiency (top) of action potential conduction. Our hypothesis is that the Mesozoic marine revolution led to a transition between these different parameter sets.
Improved AP propagation by myelin has two important consequences for the animal: First, AP propagation becomes energetically cheaper, since APs only have to be regenerated at the nodes of Ranvier, not along the whole length of the axon. The cost of AP conduction in unmyelinated axons has been shown to be significant (Crotty et al., 2006). Second, myelin improves the speed of AP conduction (Rushton, 1951; Waxman, 1997). The trade-off between optimizing myelin for either of these advantages is at the center of the hypothesis presented here.
Myelin first evolved in the Ordovician, an epoch starting 488 ma ago (Bullock et al., 1984). While such minute soft tissues like axons never fossilize, we know this from comparing the nervous systems of extant animals. Jawless vertebrates (lampreys and hagfish), which split from the rest of the vertebrate lineages before the Ordovician, lack myelin (Bullock et al., 1984). As an evolutionary innovation, myelination co-occurred with jaws and the main sections of the vertebrate brain (forebrain, midbrain, hindbrain).
In the millions of years following the first appearance of myelin, vertebrates diversified enormously. However, this diversification was curtailed at the end-Permian mass extinction. This was the most severe mass-extinction so far recorded in the history of life on earth. An estimated >95% of marine species, and probably a comparable number of terrestrial species became extinct. The causes of this massive extinction are not completely clear with, ex., vulcanism, oxidation of carbon deposits, global warming, and ocean anoxia among the non-exclusive candidate mechanisms (Hallam and Wignall, 1997). After recovery from this mass-extinction event, the composition of the marine fauna significantly changed. Animals became more mobile, and the predation pressure among the fauna increased. This is inferred from the fossil record, which shows a steep decrease in the number of sessile species and an increase of mobile and actively predating species. In addition, multiple vertebrate and invertebrate lineages increased their predatory adaptions (crushing jaws) and anti-predatory defenses (armament, spines). Bore-holes in invertebrate shells became more common among other indications. This transition is referred to as the “Mesozoic Marine Revolution” (Vermeij, 1977). Long-term, large-scale escalatory trends in evolution have been suggested previously. A change in the role of myelin toward performance-enhancement (increased conduction speed) is in accord with such trends (Vermeij, 1987). The mesozoic marine revolution was probably not an abrupt event (Walker and Brett, 2002), but a qualitative change in marine ecosystem composition undoubtedly occurred. The more mobile, actively moving animals will have had different requirements for their nervous systems. The hypothesis we present here links the observed faunal changes with a proposed repurposing of the primary role of myelin.
Evolution is an optimization process constrained by the properties of biological tissues. A number of properties of the nervous system have also been explained by reference to optimization, such as the wiring of the cortex (Chklovskii and Koulakov, 2004), the scaling of cell numbers in the primate visual system (Stevens, 2001) and the shapes of dendritic trees in the insect visual system (Torben-Nielsen and Stiefel, 2010). These studies first determined where the theoretical optimum for a certain assumed goal lay. Then, they showed that the parameters in the real biological systems agree with the theoretical predictions, and hence these systems are likely to be optimized for the initially assumed goal.
While the optimization process in biological evolution is enormously complex and simulating it is out of question, we can numerically optimize models of myelinated axons. The trajectory toward the optima will be different from biological evolution, but the parameters of optimized myelinated axons will be quite conceivably very similar to those of real axons evolved toward the same set-point.
Our hypothesis is that during the Mesozoic marine revolution, the role of myelin changed from mainly improving the energy efficiency of AP conduction to improving conduction speed. The more actively moving animals, which came to dominate ecosystems, used myelin for increasing the propagation speed of neural signals (Figure 2). This should equally true for predators as well as for the prey in need of improved escape behavior.
Figure 2. The Mesozoic marine revolution and the proposed change in the role of myelin. The cladogram shows the divergence of the vertebrate lineages over geological time (left). During the Mesozoic, starting in the triassic, the predation pressure in marine ecosystems increased (center). (*Predation pressure is estimated from durophagous (shell-crushing) predators, based on Figure 1 of Walker and Brett, 2002). Abundance of a vertebrate or invertebrate durophagous predatory lineage in a geological period was scored with 3 points, presence with 2 points and possible presence with 1 point. According to this estimate, predation pressure increases from a score of 9 in the Permian to a maximum of 20 in the Jurassic. We hypothesize that this upward shift in predatory pressure lead to a change in the structure myelin from being optimal for energy efficiency to being optimal for conduction speed (right).
This hypothesis involves evolutionary changes of non-fossilizing tissues which happened 250 ma ago. Nevertheless, it has testable consequences. Specifically, it predicts that the axonal parameters will be located at different optima for different vertebrate lineages which were already separated before the Mesozoic marine revolution. In contrast, if the Mesozoic marine revolution and the subsequent increase in marine predation pressure did not cause a change in the role of myelin, then the parameters of myelinated axons should be relatively uniform in all vertebrate lineages.
Importantly, evolutionary optimization often does not follow a single objective, but a compromise. In the case investigated here, a trade-off between the usefulness of AP transmission speed and energy efficiency has to be found. To reformulate our hypothesis in this framework, we believe that the set-point of the evolutionary optimization process shifted on the continuum toward transmission speed during the Mesozoic marine revolution, and it shifted to different new set-points in different vertebrate lineages.
In the remainder of the paper, we first describe numerical optimizations of models of myelinated axons, optimized for a range of compromises between transmission speed and energy efficiency. We then look at the empirical literature for support of the hypothesis that the myelinated axons of different vertebrate lineages have parameters corresponding to different optima. Finally, we discuss future experimental work to further test our hypothesis, as well as the implications of our hypothesis for medical research.
We first computed the optimal parameters of myelinated axons for low energy consumption and high conduction speed. Then, we compared the optimization results to the known parameters of vertebrate axons, and conclude with a discussion of myelin evolution.
We have recently implemented the optimization approach outlined above based on applying genetic algorithms to models of neurons (Stiefel and Sejnowski, 2007; Torben-Nielsen and Stiefel, 2009). The details of the implementation of the optimization algorithm can be found in (Torben-Nielsen and Stiefel, 2009).
We performed optimization on the axonal parameters related to myelin to achieve the highest conduction velocity and the most energy-efficient conduction. By optimizing for both objectives simultaneously, we could analyze the trade-off between the axonal parameters and their influence on conductance velocity and energy efficiency.
More specifically, we constructed a simulation of an axon model in which the axonal parameters (Table 1) can be specified. We then ran a simulation in which we generate an AP in the attached soma and measure the conduction velocity. The efficiency is computed from the total amount of conductance, which is a good proxy for the total current flow, and hence, a good indication of energy expenditure.
Table 1. Parameters of the axon model used in the multi-objective optimization for conduction speed and energy usage.
Optimization was done with a multi-objective genetic algorithm (Deb et al., 2002), an optimization technique loosely inspired by survival-of-the-fittest but in no way a copy of biological evolution. For such an optimization algorithm, an initial set of random axonal parameters was generated. Then, models of myelinated axons were generated from the parameter sets, and its performance was assessed according to the optimization objectives (velocity and energy-efficiency). We computed a fitness value jointly based on both objectives. Parameter sets and the models they represent which are better performing, with higher fitness values, were preferentially selected to remain in the pool of candidate models. These parameter sets were altered by mutation (random change of initial parameters) and recombination (reshuffling of parameters).
By iterating this process of assessing, selecting sets which are better than average, and mutating/recombining, the algorithm quickly yielded solutions which were performing the desired optimization criteria. The multi-objective variant of a genetic algorithm used here took two objectives (velocity and energy-efficiency) into account simultaneously. The best performing models displayed one of many optimal combinations of high conduction velocity and high energy-efficiency. The frontier of such models with optimal performance combinations is the Pareto-front.
Energetic cost was defined as the total charge and thus related to the total metabolic cost in one node: more charge offsets the ionic gradient more and thus the pumps needs to work more to achieve balance again (Sengupta et al., 2010). Velocity was computed from the time it takes for an AP to travel from one node to the next divided by the distance between the two nodes. We compute the AP velocity near the end of the axon where the AP shape has stabilized.
We compared theoretically determined optima to empirical data, and a variety of measurements from a wealth of axons of vertebrates. (Additional data points will later be included as it becomes available in order to continue testing the theory).
Types of Data
The data acquired for each species should ideally contain information at least about conduction speed, axonal thickness, myelin thickness, inter-node distance, node length and sodium and potassium concentration at the node. The data should also be available for comparable axonal tracts, across species, from the central and PNSs. Measurements should include classic electron microscopy (for measuring lengths and diameters) and immuno-gold electron microscopy for measuring ion channel densities (Poliak and Peles, 2003). Electrophysiological recordings should provide information about AP conduction speed (Gillespie and Stein, 1983).
Ideally the compared species should be from different marine branches of the vertebrate cladogram (Figure 2). The main groups are the Chondrichthyes (cartilagenous fishes), the Actinopterygii (bony fishes), the Sarcopterygii (coelochants), Dipnoi (lungfishes) and Tetrapoda (land vertebrates). The Mesozoic change of predation pressure on land is less clear than in the ocean, but land vertebrates can at least serve as a control group. As stated above, we hypothesize that these groups will have different axonal/myelin parameter sets. Groups which are descending from a group of common ancestors alive during the Mesozoic share a set of axonal/myelin parameters. Thus, the cladogram constructed from the axon/myelin parameters should allow us to estimate the time of divergence of axonal functioning in vertebrates. Our prediction is that this divergence happened during the early Mesozoic.
The optimization algorithm we used yielded axonal models with parameter sets (within the bounds of their limits, Table 1) having the fastest APs with the highest efficiency possible. The performance of the resulting models in relation to these two objectives is shown in Figure 3. A clear trade-off can be observed between velocity and energy efficiency: some parameter sets give rise to very fast AP propagation (130 m/s) at a relative high cost (2.6, nearly 3x the cost of efficient ones), while energetically efficient parameter sets (0.8) give rise to slower AP propagation (55 m/s). The asymptotic convergence to a minimal energy is defined by the lower limit on the amount of g_na (0.12 mS/cm2) and g_k (0.036 mS/cm2).
Figure 3. Optimization for both high axonal conduction velocity and low energy consumption. Speed and energetic cost for optimized models from four optimization runs (colors). Insets: Na+ and K+ currents, and the excess of Na+ current of the AP in a node at approximately three quarters the length of the axon. The 3 simulated models are optimized for different speed/efficiency trade-offs (taken from the points indicated by arrows).
The AP shape changed slightly and the amplitude of the energetically cheaper ones is lower (almost 20 mV). The axons optimized for fast AP propagation had the highest number of nodes, while the slower had fewer nodes but the exact number varies (10–40). A positive correlation exists between the number of nodes and the AP propagation velocity (Figure 4).
Figure 4. Trends in axonal parameters related to conduction velocity. Left conduction velocity vs. number of nodes. Right: conduction velocity vs. inter node diameter. Higher conduction velocities are achieved by more nodes of Ranvier and by larger axonal diameters.
Both the Chondrichthyes and the Actinopterygii had already diverged into a multitude of lineages before the Mesozoic marine revolution. We hence don't know how many independent lineages separately transitioned the Mesozoic marine revolution. Thus, lineages which have only arisen during the Mesozoic will be of particular interest. Potentially, a single ancestral species later giving rise to these lineages transitioned the Mesozoic marine revolution. An example would be the rays (batoidea) among the Chondrichthyes and the modern bony fishes (teleostei), among the Actinopterygii, two groups which emerged during the Mesozoic (Long, 2010).
Due to several factors involving the difficulty of obtaining specimen and of modifying existing methods for animals of a lesser-used species, modern neurobiology is mostly concentrating on relatively few species (“model organisms”). Thus, the desired data are available from fewer species than is desirable, with a glaring lack of data from aquatic vertebrates. There is some published information along these lines for common model organisms such as rats and cats. We compiled a table (Table 2) with all the empirical data on the aforementioned axonal parameters found in the literature.
Comparison between Numerical Optimizations and Empirical Data
We performed a literature study that is summarized in Table 2 Despite the spareness of the table, one can observe that both a higher number of nodes of Ranvier and larger inter-nodal diameter seem to be associated with higher conduction velocity in the axon. A similar connection can be found in our optimization results: higher propagation velocities are all obtained in models with the maximally allowed inter-node diameter (25 microns, Figure 4, left), and, with a higher number of nodes of Ranvier (Figure 4, right). Even if that the range of allowed inter-nodal diameters is larger (2–25 micron) in our optimization than reasonably found in nature, there is, nevertheless, a clear trend.
The available data is too sparse at this point to convincingly test the hypothesis that the role of myelin changed from energy conservation to propagation speed maximization. The empirical data lacks both complete information (all columns filled in Table 2) from one axonal tract, as well as phylogenetic diversity.
We hypothesize that during the Mesozoic marine revolution, the role of myelin changed from mainly improving the energy efficiency of AP conduction to improving conduction speed. This hypothesis has testable consequences, namely it predicts a phylogenetically biased distribution of axonal/myelin parameters. We conducted numerical optimizations of models of myelinated axons, and attempted to compare them to published information on myelinated axons. Unfortunately, the literature is too sparse, especially when it comes to non-model-organisms, and we cannot reach a definite conclusion about the state of our hypothesis. As of now, it is at least consistent, but not yet strongly supported by the empirical literature on myelinated axons.
Potentially, the peripheral nervous system (PNS) and the central nervous system (CNS) could operate along different speed/efficiency set-points. A number of time-dependent processes, such as the distinct types of oscillations seen in the CNS, could exert unique evolutionary pressures on conduction speed. On the other hand, the longer distances transversed by some PNS axons could exert pressure toward evolving higher conduction speeds as well. Better knowledge of the ontological changes in myelin and axonal conductance changes among the species included in the study would also help refine our understanding. In most cases, this kind of data is non-existent and the preponderance of the data included herein is from adult tissues.
An additional, extremely useful, piece of empirical information would be knowledge of the molecular biology of the ontogeny of myelin across multiple species. A comparison of the proteins involved in ensheathing axons with myelin, such as the oligodendrocyte myelin glycoprotein (Vourc'h and Andres, 2004), across lineages could reveal a time when these systems separated in vertebrates. Again, our hypothesis predicts a separation during the Mesozoic.
Furthermore, an evolutionary divergence of myelin during the Mesozoic marine revolution would have significant consequences for the use of different animal species as models for demyelination diseases. If the model species has myelin located at an optimum different from the human myelin-parameter optimum, its use as model species would be limited.
We hope that this hypothesis-paper serves as a motivation for our fellow neuroscientists to conduct measurements in less frequently investigated vertebrate species. The outcome of the testing of this hypothesis will have significant consequences for our understanding of vertebrate brain evolution. It will tell us if the Mesozoic marine revolution had macro-evolutionary consequences for vertebrate brains.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Boyle, M. E., Berglund, E. O., Murai, K. K., Weber, L., Peles, E., and Ranscht, B. (2001). Contactin orchestrates assembly of the septate-like junctions at the paranode in myelinated peripheral nerve. Neuron 30, 385–397. doi: 10.1016/S0896-6273(01)00296-3
Cullheim, S. (1978). Relations between cell body size, axon diameter and axon conduction velocity of cat sciatic a-motoneurons stained with horseradish peroxidase. Neurosci. Lett. 8, 17–20. doi: 10.1016/0304-3940(78)90090-3
Czarkowska, J., Jankowska, E., and Sybirska, E. (1976). Diameter and internodal length of axons of spinal interneurones excitedby group I afferents in the cat. Brain Res. 118, 119–122. doi: 10.1016/0006-8993(76)90845-3
Deschenes, M., and Landry, P. (1980). Axonal branch diameter and spacing of nodes in the terminal arborization of identified thalamic and cortical neurons. Brain Res. 191, 538–544. doi: 10.1016/0006-8993(80)91302-5
Fried, K., Hildebrand, C., and Erdily, G. (1982). Myelin sheath thickness and nerve fibres in the developing feline inferior alveolar nerve. J. Neurol. Sci. 54, 47–57. doi: 10.1016/0022-510X(82)90217-9
Friede, R. I., Benda, M., Fewitz, A., and Stoll, P. (1984). Is maximum conduction velocity the factor controlling the evolution of nerve structure. J. Neurol. Sci. 63, 369–380. doi: 10.1016/0022-510X(84)90160-6
Gillespie, M. J., and Stein, R. B. (1983). The relationship between axon diameter, myelin thickness and conduction velocity during atrophy of mammalian peripheral nerves. Brain Res. 259, 41–56. doi: 10.1016/0006-8993(83)91065-X
Ibrahim, M., Butt, A. M., and Berry, M. (1995). Relationship between myelin sheath diameter and internodal length in axons of the anterior medullary velum of the adult rat. J. Neurol. Sci. 331, 119–127. doi: 10.1016/0022-510X(95)00174-Z
Mazzatenta, A., Caleo, M., Baldaccini, N., and Maffei, L. (2001). A comparative morphometric analysis of the optic nerve in two cetacean species, the striped dolphin (Stenella coeruleoalba) and fin whale (Balaenoptera physalus). Vis. Neurosci. 18, 319–325. doi: 10.1017/S0952523801182155
Ritchie, J. M., and Rogart, R. B. (1977). Density of sodium channels in mammalian myelinated nerve fibers and nature of the axonal membrane under the myelin sheath. Proc. Natl. Acad. Sci. U.S.A. 74, 211–215. doi: 10.1073/pnas.74.1.211
Sengupta, B., Stemmler, M., Laughlin, S. B., and Niven, J. E. (2010). Action potential energy efficiency varies among neuron types in vertebrates and invertebrates. PLoS Comput. Biol. 6:e1000840. doi: 10.1371/journal.pcbi.1000840
Ulzheimer, J. C., Peles, E., Levinson, S. R., and Martini, R. (2004). Altered expression of ion channel isoforms at the node of Ranvier in P0-deficient myelin mutants. Mol. Cell. Neurosci. 25, 83–94. doi: 10.1016/j.mcn.2003.09.015
Keywords: myelin, evolution, Mesozoic marine revolution, brain energy consumption, axonal conduction speed
Citation: Stiefel KM, Torben-Nielsen B and Coggan JS (2013) Proposed evolutionary changes in the role of myelin. Front. Neurosci. 7:202. doi: 10.3389/fnins.2013.00202
Received: 05 April 2013; Paper pending published: 03 September 2013;
Accepted: 13 October 2013; Published online: 08 November 2013.
Edited by:Axel Cleeremans, Université Libre de Bruxelles, Belgium
Reviewed by:Karol Osipowicz, Jefferson Neuroscience Hospital, USA
Sven G. Meuth, Universitätsklinikum Münster, Germany
Copyright © 2013 Stiefel, Torben-Nielsen and Coggan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Klaus M. Stiefel, The MARCS Institute, University of Western Sydney, Building XB, UWS Penrith Campus, Sydney, NSW, 2751, Australia e-mail: firstname.lastname@example.org