Modeling of microstructural kinematics during simple elongation of central nervous system tissue

Damage to axons and glial cells in the central nervous system (CNS) white matter is a nearly universal feature of traumatic brain injury, yet it is not clear how the tissue mechanical deformations are transferred to the cellular components of the CNS. Defining how cellular deformations relate to the applied tissue deformation field can both highlight cellular populations at risk for mechanical injury, and define the fraction of cells in a specific population that will exhibit damage. In this investigation, microstructurally based models of CNS white matter were developed and tested against measured transformations of the CNS tissue microstructure under simple elongation. Results show that axons in the unstretched optic nerves were significantly wavy or undulated, where the measured axonal path length was greater than the end-to-end distance of the axon. The average undulation parameter--defined as the true axonal length divided by the end-to-end length--was 1.13. In stretched nerves, mean axonal undulations decreased with increasing applied stretch ratio (lambda)--the mean undulation values decreased to 1.06 at lambda = 1.06, 1.04 at lambda = 1.12, and 1.02 at lambda = 1.25. A model describing the gradual coupling, or tethering, of the axons to the surrounding glial cells best fit the experimental data. These modeling efforts indicate the fraction of the axonal and glial populations experiencing deformation increases with applied elongation, consistent with the observation that both axonal and glial cell injury increases at higher levels of white matter injury. Ultimately, these results can be used in conjunction with computational simulations of traumatic brain injury to aid in establishing the relative risk of cellular structures in the CNS white matter to mechanical injury.     

Computational modeling of axonal microtubule bundles under tension

Microtubule bundles cross-linked by tau protein serve a variety of neurological functions including maintaining mechanical integrity of the axon, promoting axonal growth, and facilitating cargo transport. It has been observed that axonal damage in traumatic brain injury leads to bundle disorientation, loss of axonal viability, and cognitive impairment. This study investigates the initial mechanical response of axonal microtubule bundles under uniaxial tension using a discrete bead-spring representation. Mechanisms of failure due to traumatic stretch loading and their impact on the mechanical response and stability are also characterized. This study indicates that cross-linked axonal microtubule bundles in tension display stiffening behavior similar to a power-law relationship from nonaffine network deformations. Stretching of cross-links and microtubule bending were the primary deformation modes at low stresses. Microtubule stretch was negligible up to tensile stresses of ～1 MPa. Bundle failure occurred by failure of cross-links leading to pull-out of microtubules and loss of bundle integrity. This may explain the elongation, undulation, and delayed elasticity of axons following traumatic stretch loading. More extensively cross-linked bundles withstood higher tensile stresses before failing. The bundle mechanical behavior uncovered by these computational techniques should guide future experiments on stretch-injured axons.     

Diffusible signals and fasciculated growth in reticulospinal axon pathfinding in the hindbrain

We have addressed the control of longitudinal axon pathfinding in the developing hindbrain, including the caudal projections of reticular and raphe neurons. To test potential sources of guidance signals, we assessed axon outgrowth from embryonic rat hindbrain explants cultured in collagen gels at a distance from explants of midbrain-hindbrain boundary (isthmus), caudal hindbrain, or cervical spinal cord. Our results showed that the isthmus inhibited caudally directed axon outgrowth by 80% relative to controls, whereas rostrally directed axon outgrowth was unaffected. Moreover, caudal hindbrain or cervical spinal cord explants did not inhibit caudal axons. Immunohistochemistry for reticular and raphe neuronal markers indicated that the caudal, but not the rostral projections of these neuronal subpopulations were inhibited by isthmic explants. Companion studies in chick embryos showed that, when the hindbrain was surgically separated from the isthmus, caudal reticulospinal axon projections failed to form and that descending pioneer axons of the medial longitudinal fasciculus (MLF) play an important role in the caudal reticulospinal projection. Taken together, these results suggest that diffusible chemorepellent or nonpermissive signals from the isthmus and substrate-anchored signals on the pioneer MLF axons are involved in the caudal direction of reticulospinal projections and might influence other longitudinal axon projections in the brainstem.     

Changes within maturing neurons limit axonal regeneration in the developing spinal cord

Embryonic birds and mammals display a remarkable ability to regenerate axons after spinal injury, but then lose this ability during a discrete developmental transition. To explain this transition, previous research has emphasized the emergence of myelin and other inhibitory factors in the environment of the spinal cord. However, research in other CNS tracts suggests an important role for neuron-intrinsic limitations to axon regeneration. Here we re-examine this issue quantitatively in the hindbrain-spinal projection of the embryonic chick. Using heterochronic cocultures we show that maturation of the spinal cord environment causes a 55% reduction in axon regeneration, while maturation of hindbrain neurons causes a 90% reduction. We further show that young neurons transplanted in vivo into older spinal cord can regenerate axons into myelinated white matter, while older axons regenerate poorly and have reduced growth cone motility on a variety of growth-permissive ligands in vitro, including laminin, L1, and N-cadherin. Finally, we use video analysis of living growth cones to directly document an age-dependent decline in the motility of brainstem axons. These data show that developmental changes in both the spinal cord environment and in brainstem neurons can reduce regeneration, but that the effect of the environment is only partial, while changes in neurons by themselves cause a nearly complete reduction in regeneration. We conclude that maturational events within neurons are a primary cause for the failure of axon regeneration in the spinal cord.     

Micromechanical analysis of brain’s diffuse axonal injury

Computational models are important tools which help researchers understand traumatic brain injury (TBI). A mechanistic multi-scale numerical approach is introduced to quantify diffuse axonal injury (DAI), the most important mechanism of TBI, induced by a mechanical insult at micro-scale regions of the white matter or voxels where fiber orientations are the same. Using the mechanical properties of a single axon with a viscoelastic constitutive relation and its functional failure in terms of electrophysiological impairment, a numerical 2D micro-level lattice method is implemented to directly analyze the percentage of injured axons in a voxel containing a bundle of axons all with the same orientation under biaxial stretches. Reference micro-injury maps are then developed with the input parameters based on the principal strain or stretch values and their direction with respect to axons, which provide the percentage of injured axons in the voxel of interest as the output. The methodology is independent of any statistical analyses of the accident data and medical reports to derive probabilistic injury risk curves for DAI. Avoiding a structurally detailed full finite element head model, this study proposes a micro-mechanical approach which considers the anatomical structure of neural axons in the white matter together with their mechanical properties using a numerical lattice method to analyze the brain’s diffuse axonal injury. This work has the potential to help develop safer prevention tools and more effective diagnosis methods for DAI.     

Pathological changes in the white matter after spinal contusion injury in the rat

It has been shown previously that after spinal cord injury, the loss of grey matter is relatively faster than loss of white matter suggesting interventions to save white matter tracts offer better therapeutic possibilities. Loss of white matter in and around the injury site is believed to be the main underlying cause for the subsequent loss of neurological functions. In this study we used a series of techniques, including estimations of the number of axons with pathology, immunohistochemistry and mapping of distribution of pathological axons, to better understand the temporal and spatial pathological events in white matter following contusion injury to the rat spinal cord. There was an initial rapid loss of axons with no detectable further loss beyond 1 week after injury. Immunoreactivity for CNPase indicated that changes to oligodendrocytes are rapid, extending to several millimetres away from injury site and preceding much of the axonal loss, giving early prediction of the final volume of white matter that survived. It seems that in juvenile rats the myelination of axons in white matter tracts continues for some time, which has an important bearing on interpretation of our, and previous, studies. The amount of myelin debris and axon pathology progressively decreased with time but could still be observed at 10 weeks after injury, especially at more distant rostral and caudal levels from the injury site. This study provides new methods to assess injuries to spinal cord and indicates that early interventions are needed for the successful sparing of white matter tracts following injury.     

Correlations between tissue-level stresses and strains and cellular damage within the guinea pig spinal cord white matter

Strain magnitude, strain rate, axon location, axon size, and the local tissue stress state have been proposed as the mechanisms governing primary cellular damage within the spinal cord parenchyma during slow compression injury. However, the mechanism of axon injury has yet to be fully elucidated. The objective of this study was to correlate cellular damage within the guinea pig spinal cord white matter, quantified by a horseradish peroxidase (HRP) exclusion test, with tissue-level stresses and strains using a combined experimental and computational approach. Force-deformation curves were acquired by transversely compressing strips of guinea pig spinal cord white matter at a quasi-static rate. Hyperelastic material parameters, derived from a Mooney-Rivlin constitutive law, were varied within a nonlinear, plane strain finite element model of the white matter strips until the computational force-deformation curve converged to the experimental results. In addition, white matter strips were subjected to nominal compression levels of 25%, 50%, 70%, and 90% to assess axonal damage by quantifying HRP uptake. HRP uptake density increased with tissue depth and with increased nominal compression. Using linear and nonlinear regression analyses, the strongest correlations with HRP uptake density were found for groups of tissue-level stresses and groups of log-transformed tissue-level strains.     

Identified motor terminals in Drosophila larvae show distinct differences in morphology and physiology

In Drosophila, the type I motor terminals innervating the larval ventral longitudinal muscle fibers 6 and 7 have been the most popular preparation for combining synaptic studies with genetics. We have further characterized the normal morphological and physiological properties of these motor terminals and the influence of muscle size on terminal morphology. Using dye-injection and physiological techniques, we show that the two axons supplying these terminals have different innervation patterns: axon 1 innervates only muscle fibers 6 and 7, whereas axon 2 innervates all of the ventral longitudinal muscle fibers. This difference in innervation pattern allows the two axons to be reliably identified. The terminals formed by axons 1 and 2 on muscle fibers 6 and 7 have the same number of branches; however, axon 2 terminals are approximately 30% longer than axon 1 terminals, resulting in a corresponding greater number of boutons for axon 2. The axon 1 boutons are approximately 30% wider than the axon 2 boutons. The excitatory postsynaptic potential (EPSP) produced by axon 1 is generally smaller than that produced by axon 2, although the size distributions show considerable overlap. Consistent with vertebrate studies, there is a correlation between muscle fiber size and terminal size. For a single axon, terminal area and length, the number of terminal branches, and the number of boutons are all correlated with muscle fiber size, but bouton size is not. During prolonged repetitive stimulation, axon 2 motor terminals show synaptic depression, whereas axon 1 EPSPs facilitate. The response to repetitive stimulation appears to be similar at all motor terminals of an axon.     

Commissural axon pathfinding on the contralateral side of the floor plate: a role for B-class ephrins in specifying the dorsoventral position of longitudinally projecting commissural axons.

In both invertebrate and lower vertebrate species, decussated commissural axons travel away from the midline and assume positions within distinct longitudinal tracts. We demonstrate that in the developing chick and mouse spinal cord, most dorsally situated commissural neuron populations extend axons across the ventral midline and through the ventral white matter along an arcuate trajectory on the contralateral side of the floor plate. Within the dorsal (chick) and intermediate (mouse) marginal zone, commissural axons turn at a conserved boundary of transmembrane ephrin expression, adjacent to which they form a discrete ascending fiber tract. In vitro perturbation of endogenous EphB-ephrinB interactions results in the failure of commissural axons to turn at the appropriate dorsoventral position on the contralateral side of the spinal cord; consequently, axons inappropriately invade more dorsal regions of B-class ephrin expression in the dorsal spinal cord. Taken together, these observations suggest that B-class ephrins act locally during a late phase of commissural axon pathfinding to specify the dorsoventral position at which decussated commissural axons turn into the longitudinal axis.     

The number, size, and type of axons in rat subcortical white matter on left and right sides: A stereological, ultrastructural study

Abundant evidence indicates important functional differences between the two cerebral hemispheres of humans, although the cellular basis of these differences is unknown. A recent hypothesis proposes that these functional differences depend on differences between sides in the “repertoire” of axonal conduction delays for cortico-cortical axons. In morphological terms this corresponds to differences in caliber, or proportion, of myelinated versus unmyelinated axons. Several behavioural studies have indicated that cerebral asymmetry occurs in rodents, in which rigorous morphological analysis is possible. The hypothesis was therefore tested for the first time in adult male Wistar rats, using transmission electron microscopy and stereological methods. Subcortical white matter was compared between left and right sides in three regions (frontal, parietal, and occipital). The average caliber and numerical density of unmyelinated and myelinated axons was compared between sides and between regions. All data were corrected for shrinkage. No significant differences between sides were found in the average caliber of either type of axon in any region. The numerical density of either type of axon also yielded no significant differences between sides in any region. Significant differences were evident between regions in both caliber and numerical density of the two axonal types, and these quantitative data are reported. The proportion of unmyelinated axons in the lateral white matter was also higher than in previous studies of hemispheric white matter that studied the corpus callosum. The present study provides no evidence supporting the hypothesis that functional hemispheric specialization is due to differences in axonal number, caliber or type.     

Lateral neuron--glia interactions steer the response of axons to the Robo code

Glia are required for axon pathfinding along longitudinal trajectories, but it is unknown how this relates to the molecular paradigm of axon guidance across the midline. Most interneuron axons in bilateral organisms cross the midline only once. Preventing them from recrossing the midline requires the expression of Robo receptors on the axons. These sense the repulsive signal Slit, which is produced by the midline. The lateral positioning of longitudinal axons depends on the response to Slit by the combination of Robo receptors expressed by the axons, on selective fasciculation, and on longitudinal (lateral) glia. Here, we analyse how longitudinal glia influence reading of the 'Robo code' by axons. We show that whereas loss of robo1 alone only affects the most medial axons, loss of both glial cells missing (gcm) and robo1 causes a severe midline collapse of longitudinal axons, similar to that caused by the loss of multiple Robo receptors. Furthermore, whereas ectopic expression of robo2 is sufficient to displace the medial MP2 axons along a more lateral trajectory, this does not occur in gcm-robo1 double-mutant embryos, where axons either do not extend at all or they misroute exiting the CNS. Hence, lateral neuron-glia interactions steer the response of axons to the Robo code.     

Evidence for Cooperative Selection of Axons for Myelination by Adjacent Oligodendrocytes in the Optic Nerve

The cellular mechanisms that regulate the topographic arrangement of myelin internodes along axons remain largely uncharacterized. Recent clonal analysis of oligodendrocyte morphologies in the mouse optic nerve revealed that adjacent oligodendrocytes frequently formed adjacent internodes on one or more axons in common, whereas oligodendrocytes in the optic nerve were never observed to myelinate the same axon more than once. By modelling the process of axonal selection at the single cell level, we demonstrate that internode length and primary process length constrain the capacity of oligodendrocytes to myelinate the same axon more than once. On the other hand, probabilistic analysis reveals that the observed juxtaposition of myelin internodes among common sets of axons by adjacent oligodendrocytes is highly unlikely to occur by chance. Our analysis may reveal a hitherto unknown level of communication between adjacent oligodendrocytes in the selection of axons for myelination. Together, our analyses provide novel insights into the mechanisms that define the spatial organization of myelin internodes within white matter at the single cell level.     

Isolation and characterization of axolemma-enriched fractions from discrete areas of bovine CNS

Myelinated axons were isolated by flotation from bovine pons, middle cerebellar peduncle, cervical spinal cord and three regions of the subcortical white matter. The myelinated axons were osmotically and mechanically shocked, followed by fractionation on a linear 15% sucrose to 45% sucrose density gradient. Axolemma-enriched fractions (AEF) found in the 28% to 32% sucrose region of the gradient from brainstem and cord white matter had high acetylcholinesterase (AChE) while little or nil AChE activity was found in corresponding AEF derived from the subcortical white matter. Morphologically, the subcortical white matter from all regions contained a heterogeneous population of well-myelinated to thinly myelinated axons, while brainstem and cord regions contained a more homogeneous population of well-myelinated axons. Histochemical analysis of AChE localized this enzyme to axonal elements. The AEF derived from any white matter source had similar polypeptide compositions. AEF derived from subcortical white matter contained two-fold more myelin basic protein and a three-fold greater content of 2 3 cyclic nucleotide 3 phosphodiesterase (CNP) compared with AEF derived from well myelinated white matter. We conclude that the purity of the AEF is related to the degree of myelination of the white matter from which the AEF is derived. Homogeneously well myelinated white matter (pons, cerebellar peduncle, cervical spinal cord) yields the highest purity AEF, as judged by the low CNP and myelin basic protein content and highest enrichment in AChE specific activity.     

Mineralized collagen fibril network spatial arrangement influences cortical bone fracture behavior

Bone is a hierarchical material exhibiting different fracture mechanisms at each length scale. At the submicroscale, the bone is composed of mineralized collagen fibrils (MCF). At this scale, the fracture processes in cortical bone have not been extensively studied in the literature. In this study, the influence of MCF size and orientation on the fracture behavior of bone under both transverse and longitudinal loading was investigated using novel 3D models of MCF networks with explicit representation of extra-fibrillar matrix. The simulation results showed that separation between MCFs was the main cause of damage and failure under transverse loading whereas under longitudinal loading, the main damage and failure mechanism was MCF rupture. When the MCF network was loaded in the transverse direction the mechanical properties increased as the orientation of fibrils deviated farther from the main fibril orientation whereas the opposite trend was observed under longitudinal loading. The fracture energy was much larger in longitudinal than transverse loading. MCF diameter variation did not affect the mechanical properties under longitudinal loading but led to higher mechanical properties with increasing MCF diameter under transverse loading. The new modeling framework established in this study generate unique information on the effect of MCF network spatial arrangement on the fracture behavior of bone at the submicroscale which is not currently possible to measure via experiments. This unique information may improve the understanding of how structural alterations at the submicroscale due to disease, age-related changes, and treatments affect the fracture processes at larger length scales.     

Synaptic connections of physiologically identified geniculocortical axons in kitten cortical area 17.

Single geniculocortical axons were recorded in the cortical white matter of kittens and adult cats by using micropipettes filled with horseradish peroxidase (HRP). Of 41 axons recovered in 4-5 week old kittens, three well-filled axons arborized in area 17; the remainder were incomplete or arborized in area 18. One axon had Y-like physiological properties, two were X-like. They were recovered from two 34-day-old kittens. All three axons formed clustered arborizations, mainly in layer 4A. Electron microscopic (EM) analysis of 50 boutons from kitten and 38 boutons from adult controls revealed that the boutons from kitten made synapses more frequently on spines (91% of targets) than did the boutons from the adult (71%). One X-like axon in kitten also had a collateral projection that made synapses in layer 1; this has not been seen in adult cats. In overall extent, the axons from kitten fell within the adult range.     

Synaptic specificity in the first instar cockroach: patterns of monosynaptic input from filiform hair afferents to giant interneurons

Summary Direct evidence for monosynaptic connections between filiform hair sensory axons and giant interneurons (GIs) in the first instar cockroach, Periplaneta americana, was obtained using intracellular recording and HRP injection followed by electron microscopy. GIs 1–6 all receive monosynaptic input from at least one filiform afferent axon. GI1, GI2 and GI5 receive input only from the medial (M) axon, while GI3, GI4 and GI6 receive input from both M and lateral (L) axons. The dendrites of GI3 and GI6 which are contralateral to the cell bodies receive input from both axons whereas the smaller ipsilateral dendritic fields have synapses only from the L axon. GI5 has M axon input only onto its contralateral dendrites. In 50% of preparations GI7 receives weak input from the ipsilateral L axon. There is no obvious relationship between the morphology of the giant interneurons and the pattern of input they receive from the filiform afferents.Abbreviations GI giant interneuron - HRP horseradish peroxidase - L lateral axon - M medial axon     

Alpha 1E subunit of the R-type calcium channel is associated with myelinogenesis

During myelinogenesis, we found an exceedingly strong, transient expression of the α_1E gene for the R-type voltage-gated calcium channel in CNS white matter. This immunoreactivity appeared in glial cells along specific pathways of the brainstem, cerebellum, and telencephalon. The reactivity followed a wave that progressed from the brainstem at P5, to the cerebellar peduncles by P8, the arbor vitae by P14, and the granular layer by P17. The reactivity-peaked about 3–4 days later and decreased gradually to become negligible in all areas before adulthood. Ultrastructural analysis confirmed that α_1E immunoreactivity was located in oligodendroglial somata, their projections, paranodal wraps and loose myelin sheaths. There was a distinct association of the channel protein reactivity on oligodendroglial membranes in contact with the axon. We propose that glial projections, contacting axons, sense axonal firing through small K⁺ currents and open the high voltage R-type calcium channels to signal myelination.     

Genes necessary for directed axonal elongation or fasciculation in C. elegans.

The outgrowth of single axons through different cellular environments requires distinct sets of genes in the nematode C. elegans. Three genes are required for the pioneering circumferential outgrowth of identified motor neuron axons between the lateral hypodermal cell membrane and the basal lamina. Three other genes are required for the longitudinal outgrowth of these axons along preexisting axon bundles as well as for the fasciculation of axons within these neuron bundles. Five additional genes are required for circumferential outgrowth, longitudinal outgrowth, and fasciculation; mutations in three of these genes disrupt axon ultrastructure, suggesting that they function in axon formation rather than in axon guidance.     

A Novel Approach for Studying the Physiology and Pathophysiology of Myelinated and Non-Myelinated Axons in the CNS White Matter

Advances in brain connectomics set the need for detailed knowledge of functional properties of myelinated and non-myelinated (if present) axons in specific white matter pathways. The corpus callosum (CC), a major white matter structure interconnecting brain hemispheres, is extensively used for studying CNS axonal function. Unlike another widely used CNS white matter preparation, the optic nerve where all axons are myelinated, the CC contains also a large population of non-myelinated axons, making it particularly useful for studying both types of axons. Electrophysiological studies of optic nerve use suction electrodes on nerve ends to stimulate and record compound action potentials (CAPs) that adequately represent its axonal population, whereas CC studies use microelectrodes (MEs), recording from a limited area within the CC. Here we introduce a novel robust isolated "whole" CC preparation comparable to optic nerve. Unlike ME recordings where the CC CAP peaks representing myelinated and non-myelinated axons vary broadly in size, "whole" CC CAPs show stable reproducible ratios of these two main peaks, and also reveal a third peak, suggesting a distinct group of smaller caliber non-myelinated axons. We provide detailed characterization of "whole" CC CAPs and conduction velocities of myelinated and non-myelinated axons along the rostro-caudal axis of CC body and show advantages of this preparation for comparing axonal function in wild type and dysmyelinated shiverer mice, studying the effects of temperature dependence, bath-applied drugs and ischemia modeled by oxygen-glucose deprivation. Due to the isolation from gray matter, our approach allows for studying CC axonal function without possible "contamination" by reverberating signals from gray matter. Our analysis of "whole" CC CAPs revealed higher complexity of myelinated and non-myelinated axonal populations, not noticed earlier. This preparation may have a broad range of applications as a robust model for studying myelinated and non-myelinated axons of the CNS in various experimental models.     

Longitudinal axons are guided by Slit/Robo signals from the floor plate

Longitudinal axons grow long distances along precise pathways to connect major CNS regions. However, during embryonic development, it remains largely undefined how the first longitudinal axons choose specific positions and grow along them. Here, we review recent evidence identifying a critical role for Slit/Robo signals to guide pioneer longitudinal axons in the embryonic brain stem. These studies indicate that Slit/Robo signals from the floor plate have dual functions: to repel longitudinal axons away from the ventral midline, and also to maintain straight longitudinal growth. These dual functions likely cooperate with other guidance cues to establish the major longitudinal tracts in the brain.Key words: Slit, Robo, longitudinal axon, hindbrain, axon guidance