Temporal trajectories of normal myelination and axonal development assessed by quantitative macromolecular and diffusion MRI: Ultrastructural and immunochemical validation in a rabbit model

INTRODUCTION

Quantitative and non-invasive measures of brain myelination and maturation during development are of great importance to both clinical and translational research communities. While the metrics derived from diffusion tensor imaging, are sensitive to developmental changes and some pathologies, they remain difficult to relate to the actual microstructure of the brain tissue. The advent of advanced model-based microstructural metrics requires histological validation. The purpose of the study was to validate novel, model-based MRI techniques, such as macromolecular proton fraction mapping (MPF) and neurite orientation and dispersion indexing (NODDI), against histologically derived indexes of myelination and microstructural maturation at various stages of development.

METHODS

New Zealand White rabbit kits underwent serial in-vivo MRI examination at postnatal days 1, 5, 11, 18, and 25, and as adults. Multi-shell, diffusion-weighted experiments were processed to fit NODDI model to obtain estimates, intracellular volume fraction (ICVF) and orientation dispersion index (ODI). Macromolecular proton fraction (MPF) maps were obtained from three source (MT-, PD-, and T1-weighted) images. After MRI sessions, a subset of animals was euthanized and regional samples of gray and white matter were taken for western blot analysis, to determine myelin basic protein (MBP), and electron microscopy, to estimate axonal, myelin fractions and g-ratio.

RESULTS

MPF of white matter regions showed a period of fast growth between P5 and P11 in the internal capsule, with a later onset in the corpus callosum. This MPF trajectory was in agreement with levels of myelination in the corresponding brain region, as assessed by western blot and electron microscopy. In the cortex, the greatest increase of MPF occurred between P18 and P26. In contrast, myelin, according to MBP western blot, saw the largest hike between P5 and P11 in the sensorimotor cortex and between P11 and P18 in the frontal cortex, which then seemingly plateaued after P11 and P18 respectively. G-ratio by MRI markers decreased with age in the white matter. However, electron microscopy suggest a relatively stable g-ratio throughout development.

CONCLUSION

Developmental trajectories of MPF accurately reflected regional differences of myelination rate in different cortical regions and white matter tracts. MRI-derived estimation of g-ratio was inaccurate during early development, likely due to the overestimation of axonal volume fraction by NODDI due to the presence of a large proportion of unmyelinated axons.

Longitudinal microstructural MRI markers of demyelination and neurodegeneration in early relapsing-remitting multiple sclerosis: Magnetisation transfer, water diffusion and g-ratio

INTRODUCTION

Quantitative microstructural MRI, such as myelin-sensitive magnetisation transfer ratio (MTR) or saturation (MTsat), axon-sensitive water diffusion Neurite Orientation Dispersion and Density Imaging (NODDI), and the aggregate g-ratio, may provide more specific markers of white matter integrity than conventional MRI for early patient stratification in relapsing-remitting multiple sclerosis (RRMS). The aim of this study was to determine the sensitivity of such markers to longitudinal pathological change within cerebral white matter lesions (WML) and normal-appearing white matter (NAWM) in recently diagnosed RRMS.

METHODS

Seventy-nine people with recently diagnosed RRMS, from the FutureMS longitudinal cohort, were recruited to an extended MRI protocol at baseline and one year later. Twelve healthy volunteers received the same MRI protocol, repeated within two weeks. Ethics approval and written informed consent were obtained. 3T MRI included magnetisation transfer, and multi-shell diffusion-weighted imaging. NAWM and whole brain were segmented from 3D T1-weighted MPRAGE, and WML from T2-weighted FLAIR. MTR, MTsat, NODDI isotropic (ISOVF) and intracellular (ICVF) volume fractions, and g-ratio (calculated from MTsat and NODDI data) were measured within WML and NAWM. Brain parenchymal fraction (BPF) was also calculated. Longitudinal change in BPF and microstructural metrics was assessed with paired t-tests (α = 0.05) and linear mixed models, adjusted for confounding factors with False Discovery Rate (FDR) correction for multiple comparisons. Longitudinal changes were compared with test-retest Bland-Altman limits of agreement from healthy control white matter. The influence of longitudinal change on g-ratio was explored through post-hoc analysis in silico by computing g-ratio with realistic simulated MTsat and NODDI values.

RESULTS

In NAWM, g-ratio and ICVF increased, and MTsat decreased over one year (adjusted mean difference = 0.007, 0.005, and -0.057 respectively, all FDR-corrected p < 0.05). There was no significant change in MTR, ISOVF, or BPF. In WML, MTsat, NODDI ICVF and ISOVF increased over time (adjusted mean difference = 0.083, 0.024 and 0.016, respectively, all FDR-corrected p < 0.05). Group-level longitudinal changes exceeded test-retest limits of agreement for NODDI ISOVF and ICVF in WML only. In silico analysis showed g-ratio may increase due to a decrease in MTsat or ISOVF, or an increase in ICVF.

DISCUSSION

G-ratio and MTsat changes in NAWM over one year may indicate subtle myelin loss in early RRMS, which were not apparent with BPF or NAWM MTR. Increases in NAWM and WML NODDI ICVF were not anticipated, and raise the possibility of axonal swelling or morphological change. Increases in WML MTsat may reflect myelin repair. Changes in NODDI ISOVF are more likely to reflect alterations in water content. Competing MTsat and ICVF changes may account for the absence of g-ratio change in WML. Longitudinal changes in microstructural measures are significant at a group level, however detection in individual patients in early RRMS is limited by technique reproducibility.

CONCLUSION

MTsat and g-ratio are more sensitive than MTR to early pathological changes in RRMS, but complex dependence of g-ratio on NODDI parameters limit the interpretation of aggregate measures in isolation. Improvements in technique reproducibility and validation of MRI biophysical models across a range of pathological tissue states are needed.

Cerebral aggregate g-ratio mapping using magnetic resonance relaxometry and diffusion tensor imaging to investigate sex and age-related differences in white matter microstructure

Axonal demyelination is a cardinal feature of aging and age-related diseases. The g-ratio, mathematically defined as the inner-to-outer diameter of a myelinated axon, is used as a structural index of optimal axonal myelination and has been shown to represent a sensitive imaging biomarker of microstructural integrity. Several magnetic resonance imaging (MRI) methods for whole-brain mapping of aggregate g-ratio have been introduced. Computation of the aggerate g-ratio requires estimates of the myelin volume fraction (MVF) and the axonal volume fraction (AVF). While accurate determinations of MVF and AVF can be obtained through multicomponent relaxometry or diffusion analyses, respectively, these methods require lengthy acquisition times making their implementation challenging in a clinical context. Therefore, any attempt to overcome this drawback is needed. Expanding on our previous work, we introduced a new MRI method for whole-brain mapping of aggregate g-ratio. This new approach is based on the use of a single-shell diffusion for AVF determination, reducing the acquisition time by approximately ~10 min from our recently introduced approach, while offering the possibility to investigate g-ratio differences in previous studies with existing data for MVF mapping and single-shell diffusion data for AVF mapping. Our comparison analysis indicates that our newly derived aggregate g-ratio values were similar to those derived from our previous method, which requires a longer acquisition time. Further, in agreement with our previous observations, we found quadratic U-shaped relationships between aggregate g-ratio and age in this much larger study cohort. However, our results show that sexual dimorphism in g-ratio was not significant in any brain region investigated.

Ultrastructural and morphometric alterations to the peripheral nerve following the administration of immunosuppressive agent tacrolimus (FK506)

Tacrolimus, a widely used immunosuppressive drug for preventing graft rejection following organ transplantation, was reported to develop neurotoxic side effects ranging from mild to severe symptoms in the literature. Rats were randomly divided into three groups as control and 2-week and 3-week treatment groups and received a 2 mg/kg/day tacrolimus by oral gavage. Animals were sacrificed and sciatic nerves obtained from all groups were fixed and processed for light and electron microscopic investigations. The myelinated fiber diameter, axon diameter, G-ratio (axon diameter/myelinated fiber diameter), and myelin thickness were also determined. The data obtained in the control and tacrolimus-treated groups were compared.The control group sciatic nerve fascicles showed normal morphology with myelinated and unmyelinated fibers. Experimental groups exhibited axonal dilatation, irregularly thickened and vacuolated myelin sheaths with separation of myelin layers. The morphometric analysis showed that the myelinated fibers of the 2-week tacrolimus-treated group displayed a moderate increase in the myelin thickness and axon and fiber diameter in comparison with the control and 3-week tacrolimus-treated groups. The G-ratio was found to be in normal range in all groups and there were no statistically significant difference.The present study indicates that the treatment with tacrolimus may produce a mild degenerative change but prolonged drug administration for 3 weeks led to improvement in morphometric and morphologic data and the normal G-ratio values, suggesting that the regeneration capacity of the myelinated fibers maintains their normal function to transmit nerve impulses.

Dendritic Degeneration of Human Auditory Nerve Fibers and Its Impact on the Spiking Pattern Under Regular Conditions and During Cochlear Implant Stimulation

Due to limitations of human in vivo studies, detailed computational models enable understanding the neural signaling in the degenerated auditory system and cochlear implants (CIs). Four human cochleae were used to quantify hearing levels depending on dendritic changes in diameter and myelination thickness from type I of the auditory nerve fibers (ANFs). Type I neurons transmit the auditory information as spiking pattern from the inner hair cells (IHCs) to the cochlear nucleus. The impact of dendrite diameter and degree of myelination on neural signal transmission was simulated for (1) synaptic excitation via IHCs and (2) stimulation from CI electrodes. An accurate three-dimensional human cochlear geometry, along with 30 auditory pathways, mimicked the CI environment. The excitation properties of electrical potential distribution induced by two CI were analyzed. Main findings: (1) The unimodal distribution of control dendrite diameters becomes multimodal for hearing loss cases; a group of thin dendrites with diameters between 0.3 and 1 μm with a peak at 0.5 μm appeared. (2) Postsynaptic currents from IHCs excite such thin dendrites easier and earlier than under control conditions. However, this advantage is lost as their conduction velocity decreases proportionally with the diameter and causes increased spike latency and jitter in soma and axon. Firing probability reduces through the soma passage due to the low intracellular current flow in thin dendrites during spiking. (3) Compared with dendrite diameter, variations in myelin thickness have a small impact on spiking performance. (4) Contrary to synaptic excitation, CIs cause several spike initiation sites in dendrite, soma region, and axon; moreover, fiber excitability reduces with fiber diameter. In a few cases, where weak stimuli elicit spikes of a target neuron (TN) in the axon, dendrite diameter reduction has no effect. However, in many cases, a spike in a TN is first initiated in the dendrite, and consequently, dendrite degeneration demands an increase in threshold currents. (5) Threshold currents of a TN and co-stimulation of degenerated ANFs in other frequency regions depend on the electrode position, including its distance to the outer wall, the cochlear turn, and the three-dimensional pathway of the TN.

Neuroprotective potential of Spirulina platensis on lesioned spinal cord corticospinal tract under experimental conditions in rat models

Spinal cord injury (SCI) results from penetrating or compressive traumatic injury to the spine in humans or by the surgical compression of the spinal cord in experimental animals. In this study, the neuroprotective potential of Spirulina platensis was investigated on ultrastructural and functional recovery of the spinal cord following surgical-induced injury. Twenty-four Sprague-Dawley rats were divided into three groups; sham group, control (trauma) group, and experimental (S. platensis) group (180 mg/kg) of eight rats each. For each group, the rats were then subdivided into two groups to allow measurement at two different timepoints (day 14 and 28) for the microscopic analysis. Rats in the control and experimental S. platensis groups were subjected to partial crush injury at the level of T12 with Inox number 2 modified forceps by compressing on the spinal cord for 30 s. Pairwise comparisons of ultrastructural grading mean scores difference between the control and experimental S. platensis groups reveals that there were significant differences on the axonal ultrastructure, myelin sheath and BBB Score on Day 28; these correlate with the functional locomotor recovery at this timepoint. The results suggest that supplementation with S. platensis induces functional recovery and effective preservation of the spinal cord ultrastructure after SCI. These findings will open new potential avenue for further research into the mechanism of S. platensis-mediated spinal cord repair.

Estimating axon conduction velocity in vivo from microstructural MRI

The conduction velocity (CV) of action potentials along axons is a key neurophysiological property central to neural communication. The ability to estimate CV in humans in vivo from non-invasive MRI methods would therefore represent a significant advance in neuroscience. However, there are two major challenges that this paper aims to address: (1) Much of the complexity of the neurophysiology of action potentials cannot be captured with currently available MRI techniques. Therefore, we seek to establish the variability in CV that can be captured when predicting CV purely from parameters that have been reported to be estimatable from MRI: inner axon diameter (AD) and g-ratio. (2) errors inherent in existing MRI-based biophysical models of tissue will propagate through to estimates of CV, the extent to which is currently unknown. Issue (1) is investigated by performing a sensitivity analysis on a comprehensive model of axon electrophysiology and determining the relative sensitivity to various morphological and electrical parameters. The investigations suggest that 85% of the variance in CV is accounted for by variation in AD and g-ratio. The observed dependency of CV on AD and g-ratio is well characterised by the previously reported model by Rushton. Issue (2) is investigated through simulation of diffusion and relaxometry MRI data for a range of axon morphologies, applying models of restricted diffusion and relaxation processes to derive estimates of axon volume fraction (AVF), AD and g-ratio and estimating CV from the derived parameters. The results show that errors in the AVF have the biggest detrimental impact on estimates of CV, particularly for sparse fibre populations (AVF<0.3). For our equipment set-up and acquisition protocol, CV estimates are most accurate (below 5% error) where AVF is above 0.3, g-ratio is between 0.6 and 0.85 and AD is high (above 4μm). CV estimates are robust to errors in g-ratio estimation but are highly sensitive to errors in AD estimation, particularly where ADs are small. We additionally show CV estimates in human corpus callosum in a small number of subjects. In conclusion, we demonstrate accurate CV estimates are possible in regions of the brain where AD is sufficiently large. Problems with estimating ADs for smaller axons presents a problem for estimating CV across the whole CNS and should be the focus of further study.

Factors Affecting Ultrastructural Quality in the Prefrontal Cortex of the Postmortem Human Brain

Electron microscopy (EM) studies of the postmortem human brain provide a level of resolution essential for understanding brain function in both normal and disease states. However, processes associated with death can impair the cellular and organelle ultrastructural preservation required for quantitative EM studies. Although postmortem interval (PMI), the time between death and preservation of tissue, is thought to be the most influential factor of ultrastructural quality, numerous other factors may also influence tissue preservation. The goal of the present study was to assess the effects of pre- and postmortem factors on multiple components of ultrastructure in the postmortem human prefrontal cortex. Tissue samples from 30 subjects were processed using standard EM histochemistry. The primary dependent measure was number of identifiable neuronal profiles, and secondary measures included presence and/or integrity of synapses, mitochondria, and myelinated axonal fibers. Number of identifiable neuronal profiles was most strongly affected by the interaction of PMI and pH, such that short PMIs and neutral pH values predicted the best preservation. Secondary measures were largely unaffected by pre- and postmortem factors. Together, these data indicate that distinct components of the neuropil are differentially affected by PMI and pH in postmortem human brain.

Influence of diet on axonal damage in the EAE mouse model of multiple sclerosis

Numerous studies have shown that diet influences the development of autoimmune diseases. However, the influence of diet on axonal damage occurring in EAE (experimental autoimmune encephalomyelitis) has not been examined. In the current study we compared changes in axonal damage and myelin thickness in spinal cords of sham- and MOG (myelin oligodendrocyte glycoprotein) peptide-immunized mice kept on a standard mouse chow (Teklad 7012) versus AIN-93 M chow which was developed for improved animal health. Despite that the development of clinical signs was similar in the 2 groups, there were significant differences in axonal caliber and myelin thickness. Following induction of EAE, axonal caliber was significantly reduced in mice fed Teklad diet, but not those fed the AIN diet. Concomitantly, myelin thickness was decreased by EAE in mice fed Teklad, but not AIN diet. Analysis of g-ratios showed that the increase in g-ratio with increasing axonal size was reduced in mice fed AIN diet. These findings demonstrate that differences in axonal pathology occur in the absence of observable differences in clinical signs, and that inter-study comparisons may be confounded by differences in dietary care.

Experimental studies of g-ratio MRI in ex vivo mouse brain

This study aimed to experimentally evaluate a previously proposed MRI method for mapping axonal g-ratio (ratio of axon diameters, measured to the inner and outer boundary of myelin). MRI and electron microscopy were used to study excised and fixed brains of control mice and three mouse models of abnormal white matter. The results showed that g-ratio measured with MRI correlated with histological measures of myelinated axon g-ratio, but with a bias that is likely due to the presence of non-myelinated axons. The results also pointed to cases where the MRI g-ratio model simplifies to be primarily a function of total myelin content.

A comparative study of oculomotor, trochlear and abducens nerves in Arabian foals

We investigated the microscopic structure of transverse sections of the oculomotor, trochlear and abducens nerves of Arabian foals using stereological methods. Bilateral nerve pairs from 2-month-old female Arabian foals were analyzed. The tissues were embedded in plastic blocks, then 1 µm thick sections were cut and stained with osmium tetroxide and methylene blue-azure II. Stereology was performed using light microscopy. Morphometry showed that the right and left pairs of nerves were similar. The transverse sectional areas of the oculomotor, trochlear and abducens nerves were 1.93 ± 0.19 mm², 0.32 ± 0.06 mm² and 0.70 ± 0.08 mm², respectively. The oculomotor nerve exhibited a significantly greater number of myelinated axons (16755 ± 1279) and trochlear (2656 ± 494) and the abducens nerves (4468 ± 447). The ratio of the axon diameter to myelinated nerve fiber diameter was 0.58, 0.55 and 0.55 for the oculomotor, trochlear and abducens nerves, respectively. Of the three nerves studied, the abducens nerve exhibited the greatest nerve fiber area, myelin area, nerve and axon diameters, and myelin thickness. The ratio of small myelinated nerve fibers was greatest in the oculomotor nerve.