Axon kinematics change during growth and development

The biomechanical implications of neck position in cervical contusion animal models of SCI

Large animal contusion models of spinal cord injury are an essential precursor to developing and evaluating treatment options for human spinal cord injury. Reducing variability in these experiments has been a recent focus as it increases the sensitivity with which treatment effects can be detected while simultaneously decreasing the number of animals required in a study. Here, we conducted a detailed review to explore if head and neck positioning in a cervical contusion model of spinal cord injury could be a factor impacting the biomechanics of a spinal cord injury, and thus, the resulting outcomes. By reviewing existing literature, we found evidence that animal head/neck positioning affects the exposed level of the spinal cord, morphology of the spinal cord, tissue mechanics and as a result the biomechanics of a cervical spinal cord injury. We posited that neck position could be a hidden factor contributing to variability. Our results indicate that neck positioning is an important factor in studying biomechanics, and that reporting these values can improve inter-study consistency and comparability and that further work needs to be done to standardize positioning for cervical spinal cord contusion injury models.

White Matter Tract-Oriented Deformation Is Dependent on Real-Time Axonal Fiber Orientation

Traumatic axonal injury (TAI) is a critical public health issue with its pathogenesis remaining largely elusive. Finite element (FE) head models are promising tools to bridge the gap between mechanical insult, localized brain response, and resultant injury. In particular, the FE-derived deformation along the direction of white matter (WM) tracts (i.e., tract-oriented strain) has been shown to be an appropriate predictor for TAI. The evolution of fiber orientation in time during the impact and its potential influence on the tract-oriented strain remains unknown, however. To address this question, the present study leveraged an embedded element approach to track real-time fiber orientation during impacts. A new scheme to calculate the tract-oriented strain was proposed by projecting the strain tensors from pre-computed simulations along the temporal fiber direction instead of its static counterpart directly obtained from diffuse tensor imaging. The results revealed that incorporating the real-time fiber orientation not only altered the direction but also amplified the magnitude of the tract-oriented strain, resulting in a generally more extended distribution and a larger volume ratio of WM exposed to high deformation along fiber tracts. These effects were exacerbated with the impact severities characterized by the acceleration magnitudes. Results of this study provide insights into how best to incorporate fiber orientation in head injury models and derive the WM tract-oriented deformation from computational simulations, which is important for furthering our understanding of the underlying mechanisms of TAI.

On the Transversely Isotropic, Hyperelastic Response of Central Nervous System White Matter Using a Hybrid Approach

A numerical and experimental hybrid approach is developed to study the constitutive behavior of the central nervous system white matter. A published transversely isotropic hyperelastic strain energy function is reviewed and used to determine stress–strain relationships for three idealized, simple loading scenarios. The proposed constitutive model is simplified to a three-parameter hyperelastic model by assuming the white matter's incompressibility. Due to a lack of experimental data in all three loading scenarios, a finite element model that accounts for microstructural axons and their kinematics is developed to simulate behaviors in simple shear loading scenarios to supplement existing uniaxial tensile test data. The parameters of the transversely isotropic hyperelastic material model are determined regressively using the hybrid data. The results highlight that a hybrid numerical virtual test coupled with experimental data, can determine the transversely isotropic hyperelastic model. It is noted that the model is not limited to small strains and can be applied to large deformations.

Mechanical regulation of oligodendrocyte biology

Oligodendrocytes (OL) are a subset of glial cells in the central nervous system (CNS) comprising the brain and spinal cord. The CNS environment is defined by complex biochemical and biophysical cues during development and response to injury or disease. In the last decade, significant progress has been made in understanding some of the key biophysical factors in the CNS that modulate OL biology, including their key role in myelination of neurons. Taken together, those studies offer translational implications for remyelination therapies, pharmacological research, identification of novel drug targets, and improvements in methods to generate human oligodendrocyte progenitor cells (OPCs) and OLs from donor stem cells in vitro. This review summarizes current knowledge of how various physical and mechanical cues affect OL biology and its implications for disease, therapeutic approaches, and generation of human OPCs and OLs.

Axons Embedded in a Tissue May Withstand Larger Deformations Than Isolated Axons Before Mechanoporation Occurs

Diffuse axonal injury (DAI) is the pathological consequence of traumatic brain injury (TBI) that most of all requires a multiscale approach in order to be, first, understood and then possibly prevented. While in fact the mechanical insult usually happens at the head (or macro) level, the consequences affect structures at the cellular (or microlevel). The quest for axonal injury tolerances has so far been addressed both with experimental and computational approaches. On one hand, the experimental approach presents challenges connected to both temporal and spatial resolution in the identification of a clear axonal injury trigger after the application of a mechanical load. On the other hand, computational approaches usually consider axons as homogeneous entities and therefore are unable to make inferences about their viability, which is thought to depend on subcellular damages. Here, we propose a computational multiscale approach to investigate the onset of axonal injury in two typical experimental scenarios. We simulated single-cell and tissue stretch injury using a composite finite element axonal model in isolation and embedded in a matrix, respectively. Inferences on axonal damage are based on the comparison between axolemma strains and previously established mechanoporation thresholds. Our results show that, axons embedded in a tissue could withstand higher deformations than isolated axons before mechanoporation occurred and this is exacerbated by the increase in strain rate from 1/s to 10/s.

Diffusion imaging of reversible and irreversible microstructural changes within the corticospinal tract in idiopathic normal pressure hydrocephalus

The symptoms of idiopathic normal pressure hydrocephalus (iNPH) can be improved by shunt surgery, but prediction of treatment outcome is not established. We investigated changes of the corticospinal tract (CST) in iNPH before and after shunt surgery by using diffusion microstructural imaging, which infers more specific tissue properties than conventional diffusion tensor imaging. Two biophysical models were used: neurite orientation dispersion and density imaging (NODDI) and white matter tract integrity (WMTI). In both methods, the orientational coherence within the CSTs was higher in patients than in controls, and some normalization occurred after the surgery in patients, indicating axon stretching and recovery. The estimated axon density was lower in patients than in controls but remained unchanged after the surgery, suggesting its potential as a marker for irreversible neuronal damage. In a Monte-Carlo simulation that represented model axons as undulating cylinders, both NODDI and WMTI separated the effects of axon density and undulation. Thus, diffusion MRI may distinguish between reversible and irreversible microstructural changes in iNPH. Our findings constitute a step towards a quantitative image biomarker that reflects pathological process and treatment outcomes of iNPH.

•

NODDI and WMTI provide markers of reversible and irreversible changes in iNPH.

•

Measures of axon orientation indicated recovery from stretching after surgery.

•

Axon density remained low after surgery, suggesting chronic neuronal damage.

•

Axon stretching in simulations differentially affected diffusion metrics.

Low-intensity pulsed ultrasound upregulates pro-myelination indicators of Schwann cells enhanced by co-culture with adipose-derived stem cells

OBJECTIVES

Peripheral nerve injuries are a common occurrence, resulting in considerable patient suffering; it also represents a major economic burden on society. To improve treatment options following peripheral nerve injuries, scientists aim to find a way to promote Schwann cell (SC) myelination to help nerves to carry out their functions effectively. In this study, we investigated myelination ability of SCs, regulated by co-culture with adipose-derived stem cells (ASCs) or low-intensity pulsed ultrasound (LIPUS), and synergistic effects of combined treatments.

MATERIALS AND METHODS

Schwann cells were co-cultured with or without ASCs, and either left untreated or treated with LIPUS for 10 min/d for 1, 4 or 7 days. Effects of LIPUS and ASC co-culture on pro-myelination indicators of SCs were analysed by real-time PCR (RT-PCR), Western blotting and immunofluorescence staining (IF).

RESULTS

Our results indicate that ASC-SC co-culture and LIPUS, together or individually, promoted mRNA levels of epidermal growth factor receptor 3 (EGFR3/ErbB3), neuregulin1 (NRG1), early growth response protein 2 (Egr2/Krox20) and myelin basic protein (MBP), with corresponding increases in protein levels of ErbB3, NRG1 and Krox20. Interestingly, combination of ASC-SC co-culture and LIPUS displayed the most remarkable effects.

CONCLUSION

We demonstrated that ASCs upregulated pro-myelination indicators of SCs by indirect contact (through co-culture) and that effects could be potentiated by LIPUS. We conclude that LIPUS, as a mechanical stress, may have potential in nerve regeneration with potential clinical relevance.

CNS cell distribution and axon orientation determine local spinal cord mechanical properties

Mechanical signaling plays an important role in cell physiology and pathology. Many cell types, including neurons and glial cells, respond to the mechanical properties of their environment. Yet, for spinal cord tissue, data on tissue stiffness are sparse. To investigate the regional and direction-dependent mechanical properties of spinal cord tissue at a spatial resolution relevant to individual cells, we conducted atomic force microscopy (AFM) indentation and tensile measurements on acutely isolated mouse spinal cord tissue sectioned along the three major anatomical planes, and correlated local mechanical properties with the underlying cellular structures. Stiffness maps revealed that gray matter is significantly stiffer than white matter irrespective of directionality (transverse, coronal, and sagittal planes) and force direction (compression or tension) (K_g= ∼130 Pa vs. K_w= ∼70 Pa); both matters stiffened with increasing strain. When all data were pooled for each plane, gray matter behaved like an isotropic material under compression; however, subregions of the gray matter were rather heterogeneous and anisotropic. For example, in sagittal sections the dorsal horn was significantly stiffer than the ventral horn. In contrast, white matter behaved transversely isotropic, with the elastic stiffness along the craniocaudal (i.e., longitudinal) axis being lower than perpendicular to it. The stiffness distributions we found under compression strongly correlated with the orientation of axons, the areas of cell nuclei, and cellular in plane proximity. Based on these morphological parameters, we developed a phenomenological model to estimate local mechanical properties of central nervous system (CNS) tissue. Our study may thus ultimately help predicting local tissue stiffness, and hence cell behavior in response to mechanical signaling under physiological and pathological conditions, purely based on histological data.

Characterization of the three-dimensional kinematic behavior of axons in central nervous system white matter

Traumatic injury to axons in white matter of the brain and spinal cord occurs primarily via tensile stretch. During injury, the stress and strain experienced at the tissue level is transferred to the microscopic axons. How this transfer occurs, and the primary constituents dictating this transfer must be better understood to develop more accurate multi-scale models of injury. Previous studies have characterized axon tortuosity and kinematic behavior in 2-dimensions (2-D), where axons have been modeled to exhibit non-affine (discrete), affine (composite-like), or switching behavior. In this study, we characterize axon tortuosity and model axon kinematic behavior in 3-dimensions (3-D). Embryonic chick spinal cords at different development stages were excised and stretched. Cords were then fixed, transversely sectioned, stained, and imaged. 3-D axon tortuosity was measured from confocal images using a custom-built MATLAB script. 2-D kinematic models previously described in Bain et al. (J Biomech Eng 125(6):798, 2003) were extended, re-derived, and validated for the 3-D case. Results showed that 3-D tortuosity decreased with stretch, exhibiting similar trends with changes in development as observed in the 2-D studies. Kinematic parameters also displayed similar general trends. Axons demonstrated more affine behavior with increasing stretch and development. In comparison with 2-D results, a smaller percentage of the populations of 3-D axons were predicted to follow pure non-affine behavior. The data and kinematic models presented herein can be incorporated into multi-scale CNS injury models, which can advance the accuracy of the models and improve the potential to identify axonal injury thresholds.

Finite Element Modeling of CNS White Matter Kinematics: Use of a 3D RVE to Determine Material Properties

Axonal injury represents a critical target area for the prevention and treatment of traumatic brain and spinal cord injuries. Finite element (FE) models of the head and/or brain are often used to predict brain injury caused by external mechanical loadings, such as explosive waves and direct impact. The accuracy of these numerical models depends on correctly determining the material properties and on the precise depiction of the tissues' microstructure (microscopic level). Moreover, since the axonal microstructure for specific regions of the brain white matter is locally oriented, the stress, and strain fields are highly anisotropic and axon orientation dependent. Additionally, mechanical strain has been identified as the proximal cause of axonal injury, which further demonstrates the importance of this multi-scale relationship. In this study, our previously developed FE and kinematic axonal models are coupled and applied to a pseudo 3-dimensional representative volume element of central nervous system white matter to investigate the multi-scale mechanical behavior. An inverse FE procedure was developed to identify material parameters of spinal cord white matter by combining the results of uniaxial testing with FE modeling. A satisfactory balance between simulation and experiment was achieved via optimization by minimizing the squared error between the simulated and experimental force-stretch curve. The combination of experimental testing and FE analysis provides a useful analysis tool for soft biological tissues in general, and specifically enables evaluations of the axonal response to tissue-level loading and subsequent predictions of axonal damage.

The DIADEM metric: comparing multiple reconstructions of the same neuron

Digital reconstructions of neuronal morphology are used to study neuron function, development, and responses to various conditions. Although many measures exist to analyze differences between neurons, none is particularly suitable to compare the same arborizing structure over time (morphological change) or reconstructed by different people and/or software (morphological error). The metric introduced for the DIADEM (DIgital reconstruction of Axonal and DEndritic Morphology) Challenge quantifies the similarity between two reconstructions of the same neuron by matching the locations of bifurcations and terminations as well as their topology between the two reconstructed arbors. The DIADEM metric was specifically designed to capture the most critical aspects in automating neuronal reconstructions, and can function in feedback loops during algorithm development. During the Challenge, the metric scored the automated reconstructions of best-performing algorithms against manually traced gold standards over a representative data set collection. The metric was compared with direct quality assessments by neuronal reconstruction experts and with clocked human tracing time saved by automation. The results indicate that relevant morphological features were properly quantified in spite of subjectivity in the underlying image data and varying research goals. The DIADEM metric is freely released open source ( http://diademchallenge.org ) as a flexible instrument to measure morphological error or change in high-throughput reconstruction projects.

The importance of axonal undulation in diffusion MR measurements: a Monte Carlo simulation study

Many axons follow wave-like undulating courses. This is a general feature of extracranial nerve segments, but is also found in some intracranial nervous tissue. The importance of axonal undulation has previously been considered, for example, in the context of biomechanics, where it has been shown that posture affects undulation properties. However, the importance of axonal undulation in the context of diffusion MR measurements has not been investigated. Using an analytical model and Monte Carlo simulations of water diffusion, this study compared undulating and straight axons in terms of diffusion propagators, diffusion-weighted signal intensities and parameters derived from diffusion tensor imaging, such as the mean diffusivity (MD), the eigenvalues and the fractional anisotropy (FA). All parameters were strongly affected by the presence of undulation. The diffusivity perpendicular to the undulating axons increased with the undulation amplitude, thus resembling that of straight axons with larger diameters. Consequently, models assuming straight axons for the estimation of the axon diameter from diffusion MR measurements might overestimate the diameter if undulation is present. FA decreased from approximately 0.7 to 0.5 when axonal undulation was introduced into the simulation model structure. Our results indicate that axonal undulation may play a role in diffusion measurements when investigating, for example, the optic and sciatic nerves and the spinal cord. The simulations also demonstrate that the stretching or compression of neuronal tissue comprising undulating axons alters the observed water diffusivity, suggesting that posture may be of importance for the outcome of diffusion MRI measurements.

Probing the influence of myelin and glia on the tensile properties of the spinal cord

Although glia have been historically classified as the structurally supporting cells of the central nervous system, their role in tissue mechanics is still largely unstudied. The influence of myelin and glia on the mechanical properties of spinal cord tissue was examined by testing embryonic day 18 chick embryo spinal cords in uniaxial tension following disruption of the glial matrix using either ethidium bromide (EB) or an antibody against galactocerebroside (alphaGalC) in the presence of complement. Demyelination was confirmed by myelin basic protein immunoreactivity and quantified using osmium tetroxide staining. A substantial loss of astrocytes and oligodendrocytes concurrent with demyelination was observed following EB injection but not alphaGalC injection. No morphological changes were observed following injection of saline or IgG with complement as controls for EB and alphaGalC. Demyelinated spinal cords demonstrated significantly lower stiffness and ultimate tensile stress than myelinated spinal cords. No significant differences were observed in the tensile response between the two demyelinating protocols. The results demonstrate that the glial matrix provides significant mechanical support to the spinal cord, and suggests that myelin and cellular coupling of axons via the glial matrix in large part dictates the tensile response of the tissue.