Multi-component T1 relaxation and magnetisation transfer in peripheral nerve

Magnetization transfer ratio of the sciatic nerve differs between patients in type 1 and type 2 diabetes

BACKGROUND

Previous studies on magnetic resonance neurography (MRN) found different patterns of structural nerve damage in type 1 diabetes (T1D) and type 2 diabetes (T2D). Magnetization transfer ratio (MTR) is a quantitative technique to analyze the macromolecular tissue composition. We compared MTR values of the sciatic nerve in patients with T1D, T2D, and healthy controls (HC).

METHODS

3-T MRN of the right sciatic nerve at thigh level was performed in 14 HC, 10 patients with T1D (3 with diabetic neuropathy), and 28 patients with T2D (10 with diabetic neuropathy). Results were subsequently correlated with clinical and electrophysiological data.

RESULTS

The sciatic nerve's MTR was lower in patients with T2D (0.211 ± 0.07, mean ± standard deviation) compared to patients with T1D (T1D 0.285 ± 0.03; p = 0.015) and HC (0.269 ± 0.05; p = 0.039). In patients with T1D, sciatic MTR correlated positively with tibial nerve conduction velocity (NCV; r = 0.71; p = 0.021) and negatively with hemoglobin A1c (r =  - 0.63; p < 0.050). In patients with T2D, we found negative correlations of sciatic nerve's MTR peroneal NCV (r =  - 0.44; p = 0.031) which remained significant after partial correlation analysis controlled for age and body mass index (r = 0.51; p = 0.016).

CONCLUSIONS

Lower MTR values of the sciatic nerve in T2D compared to T1D and HC and diametrical correlations of MTR values with NCV in T1D and T2D indicate that there are different macromolecular changes and pathophysiological pathways underlying the development of neuropathic nerve damage in T1D and T2D.

TRIAL REGISTRATION

https://classic.

CLINICALTRIALS

gov/ct2/show/NCT03022721 . 16 January 2017.

RELEVANCE STATEMENT

Magnetization transfer ratio imaging may serve as a non-invasive imaging method to monitor the diseases progress and to encode the pathophysiology of nerve damage in patients with type 1 and type 2 diabetes.

KEY POINTS

• Magnetization transfer imaging detects distinct macromolecular nerve lesion patterns in diabetes patients. • Magnetization transfer ratio was lower in type 2 diabetes compared to type 1 diabetes. • Different pathophysiological mechanisms drive nerve damage in type 1 and 2 diabetes.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Infantile status epilepticus disrupts myelin development

Temporal lobe epilepsy (TLE) is the most prevalent type of epilepsy in adults; it often starts in infancy or early childhood. Although TLE is primarily considered to be a grey matter pathology, a growing body of evidence links this disease with white matter abnormalities. In this study, we explore the impact of TLE onset and progression in the immature brain on white matter integrity and development utilising the rat model of Li-pilocarpine-induced TLE at the 12th postnatal day (P). Diffusion tensor imaging (DTI) and Black-Gold II histology uncovered disruptions in major white matter tracks (corpus callosum, internal and external capsules, and deep cerebral white matter) spreading through the whole brain at P28. These abnormalities were mostly not present any longer at three months after TLE induction, with only limited abnormalities detectable in the external capsule and deep cerebral white matter. Relaxation Along a Fictitious Field in the rotating frame of rank 4 indicated that white matter changes observed at both timepoints, P28 and P72, are consistent with decreased myelin content. The animals affected by TLE-induced white matter abnormalities exhibited increased functional connectivity between the thalamus and medial prefrontal and somatosensory cortex in adulthood. Furthermore, histological analyses of additional animal groups at P15 and P18 showed only mild changes in white matter integrity, suggesting a gradual age-dependent impact of TLE progression. Taken together, TLE progression in the immature brain distorts white matter development with a peak around postnatal day 28, followed by substantial recovery in adulthood. This developmental delay might give rise to cognitive and behavioural comorbidities typical for early-onset TLE.

Magnetic Resonance Neurography: Improved Diagnosis of Peripheral Neuropathies

Peripheral neuropathies account for the most frequent disorders seen by neurologists, and causes are manifold. The traditional diagnostic gold-standard consists of clinical neurologic examinations supplemented by nerve conduction studies. Due to well-known limitations of standard diagnostics and atypical clinical presentations, establishing the correct diagnosis can be challenging but is critical for appropriate therapies. Magnetic resonance neurography (MRN) is a relatively novel technique that was developed for the high-resolution imaging of the peripheral nervous system. In focal neuropathies, whether traumatic or due to nerve entrapment, MRN has improved the diagnostic accuracy by directly visualizing underlying nerve lesions and providing information on the exact lesion localization, extension, and spatial distribution, thereby assisting surgical planning. Notably, the differentiation between distally located, complete cross-sectional nerve lesions, and more proximally located lesions involving only certain fascicles within a nerve can hold difficulties that MRN can overcome, when basic technical requirements to achieve sufficient spatial resolution are implemented. Typical MRN-specific pitfalls are essential to understand in order to prevent overdiagnosing neuropathies. Heavily T2-weighted sequences with fat saturation are the most established sequences for MRN. Newer techniques, such as T2-relaxometry, magnetization transfer contrast imaging, and diffusion tensor imaging, allow the quantification of nerve lesions and have become increasingly important, especially when evaluating diffuse, non-focal neuropathies. Innovative studies in hereditary, metabolic or inflammatory polyneuropathies, and motor neuron diseases have contributed to a better understanding of the underlying pathomechanism. New imaging biomarkers might be used for an earlier diagnosis and monitoring of structural nerve injury under causative treatments in the future.

Whole-Brain Imaging of Subvoxel T1-Diffusion Correlation Spectra in Human Subjects

T1 relaxation and water mobility generate eloquent MRI tissue contrasts with great diagnostic value in many neuroradiological applications. However, conventional methods do not adequately quantify the microscopic heterogeneity of these important biophysical properties within a voxel, and therefore have limited biological specificity. We describe a new correlation spectroscopic (CS) MRI method for measuring how T1 and mean diffusivity (MD) co-vary in microscopic tissue environments. We develop a clinical pulse sequence that combines inversion recovery (IR) with single-shot isotropic diffusion encoding (IDE) to efficiently acquire whole-brain MRIs with a wide range of joint T1-MD weightings. Unlike conventional diffusion encoding, the IDE preparation ensures that all subvoxel water pools are weighted by their MDs regardless of the sizes, shapes, and orientations of their corresponding microscopic diffusion tensors. Accordingly, IR-IDE measurements are well-suited for model-free, quantitative spectroscopic analysis of microscopic water pools. Using numerical simulations, phantom experiments, and data from healthy volunteers we demonstrate how IR-IDE MRIs can be processed to reconstruct maps of two-dimensional joint probability density functions, i.e., correlation spectra, of subvoxel T1-MD values. In vivo T1-MD spectra show distinct cerebrospinal fluid and parenchymal tissue components specific to white matter, cortical gray matter, basal ganglia, and myelinated fiber pathways, suggesting the potential for improved biological specificity. The one-dimensional marginal distributions derived from the T1-MD correlation spectra agree well with results from other relaxation spectroscopic and quantitative MRI studies, validating the T1-MD contrast encoding and the spectral reconstruction. Mapping subvoxel T1-diffusion correlations in patient populations may provide a more nuanced, comprehensive, sensitive, and specific neuroradiological assessment of the non-specific changes seen on fluid-attenuated inversion recovery (FLAIR) and diffusion-weighted MRIs (DWIs) in cancer, ischemic stroke, or brain injury.

Reliability and reproducibility of sciatic nerve magnetization transfer imaging and T2 relaxometry

Objectives

To assess the interreader and test-retest reliability of magnetization transfer imaging (MTI) and T2 relaxometry in sciatic nerve MR neurography (MRN).

Materials and methods

In this prospective study, 21 healthy volunteers were examined three times on separate days by a standardized MRN protocol at 3 Tesla, consisting of an MTI sequence, a multi-echo T2 relaxometry sequence, and a high-resolution T2-weighted sequence. Magnetization transfer ratio (MTR), T2 relaxation time, and proton spin density (PSD) of the sciatic nerve were assessed by two independent observers, and both interreader and test-retest reliability for all readout parameters were reported by intraclass correlation coefficients (ICCs) and standard error of measurement (SEM).

Results

For the sciatic nerve, overall mean ± standard deviation MTR was 26.75 ± 3.5%, T2 was 64.54 ± 8.2 ms, and PSD was 340.93 ± 78.8. ICCs ranged between 0.81 (MTR) and 0.94 (PSD) for interreader reliability and between 0.75 (MTR) and 0.94 (PSD) for test-retest reliability. SEM for interreader reliability was 1.7% for MTR, 2.67 ms for T2, and 21.3 for PSD. SEM for test-retest reliability was 1.7% for MTR, 2.66 ms for T2, and 20.1 for PSD.

Conclusions

MTI and T2 relaxometry of the sciatic nerve are reliable and reproducible. The values of measurement imprecision reported here may serve as a guide for correct interpretation of quantitative MRN biomarkers in future studies.

Key Points

• Magnetization transfer imaging (MTI) and T2 relaxometry of the sciatic nerve are reliable and reproducible.

• The imprecision that is unavoidably associated with different scans or different readers can be estimated by the here presented SEM values for the biomarkers T2, PSD, and MTR.

• These values may serve as a guide for correct interpretation of quantitative MRN biomarkers in future studies and possible clinical applications.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-021-08072-9.

Multidimensional correlation MRI

Multidimensional correlation spectroscopy is emerging as a novel MRI modality that is well suited for microstructure and microdynamic imaging studies, especially of biological specimens. Conventional MRI methods only provide voxel-averaged and mostly macroscopically averaged information; these methods cannot disentangle intra-voxel heterogeneity on the basis of both water mobility and local chemical interactions. By correlating multiple MR contrast mechanisms and processing the data in an integrated manner, correlation spectroscopy is able to resolve the distribution of water populations according to their chemical and physical interactions with the environment. The use of a non-parametric, phenomenological representation of the multidimensional MR signal makes no assumptions about tissue structure, thereby allowing the study of microscopic structure and composition of complex heterogeneous biological systems. However, until recently, vast data requirements have confined these types of measurement to non-localized NMR applications and prevented them from being widely and successfully used in conjunction with imaging. Recent groundbreaking advancements have allowed this powerful NMR methodology to be migrated to MRI, initiating its emergence as a promising imaging approach. This review is not intended to cover the entire field of multidimensional MR; instead, it focuses on pioneering imaging applications and the challenges involved. In addition, the background and motivation that have led to multidimensional correlation MR development are discussed, along with the basic underlying mathematical concepts. The goal of the present work is to provide the reader with a fundamental understanding of the techniques developed and their potential benefits, and to provide guidance to help refine future applications of this technology.

Multidimensional correlation spectroscopic imaging of exponential decays: From theoretical principles to in vivo human applications

Multiexponential modeling of relaxation or diffusion MR signal decays is a popular approach for estimating and spatially mapping different microstructural tissue compartments. While this approach can be quite powerful, it is also limited by the fact that one-dimensional multiexponential modeling is an ill-posed inverse problem with substantial ambiguities. In this article, we present an overview of a recent multidimensional correlation spectroscopic imaging approach to this problem. This approach helps to alleviate ill-posedness by making advantageous use of multidimensional contrast encoding (e.g., 2D diffusion-relaxation encoding or 2D relaxation-relaxation encoding) combined with a regularized spatial-spectral estimation procedure. Theoretical calculations, simulations, and experimental results are used to illustrate the benefits of this approach relative to classical methods. In addition, we demonstrate an initial proof-of-principle application of this kind of approach to in vivo human MRI experiments.

Non-negative least squares computation for in vivo myelin mapping using simulated multi-echo spin-echo T2 decay data

Multi-compartment T₂ mapping has gained particular relevance for the study of myelin water in the brain. As a facilitator of rapid saltatory axonal signal transmission, myelin is a cornerstone indicator of white matter development and function. Regularized non-negative least squares fitting of multi-echo T₂ data has been widely employed for the computation of the myelin water fraction (MWF), and the obtained MWF maps have been histopathologically validated. MWF measurements depend upon the quality of the data acquisition, B₁⁺ homogeneity and a range of fitting parameters. In this special issue article, we discuss the relevance of these factors for the accurate computation of multi-compartment T₂ and MWF maps. We generated multi-echo spin-echo T₂ decay curves following the Carr-Purcell-Meiboom-Gill approach for various myelin concentrations and myelin T₂ scenarios by simulating the evolution of the magnetization vector between echoes based on the Bloch equations. We demonstrated that noise and imperfect refocusing flip angles yield systematic underestimations in MWF and intra-/extracellular water geometric mean T₂ (gmT₂ ). MWF estimates were more stable than myelin water gmT₂ time across different settings of the T₂ analysis. We observed that the lower limit of the T₂ distribution grid should be slightly shorter than TE₁ . Both TE₁ and the acquisition echo spacing also have to be sufficiently short to capture the rapidly decaying myelin water T₂ signal. Among all parameters of interest, the estimated MWF and intra-/extracellular water gmT₂ differed by approximately 0.13-4 percentage points and 3-4 ms, respectively, from the true values, with larger deviations observed in the presence of greater B₁⁺ inhomogeneities and at lower signal-to-noise ratio. Tailoring acquisition strategies may allow us to better characterize the T₂ distribution, including the myelin water, in vivo.

Magnetization transfer ratio: a quantitative imaging biomarker for 5q spinal muscular atrophy

BACKGROUND AND PURPOSE

We quantified peripheral nerve lesions in adults with 5q-linked spinal muscular atrophy (SMA) type 3 by analysing the magnetization transfer ratio (MTR) of the sciatic nerve, and tested its potential as a novel biomarker for macromolecular changes.

METHODS

Eighteen adults with SMA 3 (50% SMA 3a, 50% SMA 3b) and 18 age-/sex-matched healthy controls prospectively underwent magnetization transfer contrast imaging in a 3-Tesla magnetic resonance scanner. Two axial three-dimensional gradient echo sequences, with and without an off-resonance saturation rapid frequency pulse, were performed at the right distal thigh. Sciatic nerve regions of interest were manually traced on 10 consecutive axial slices in the images generated without off-resonance saturation, and then transferred to corresponding slices generated by the sequence with the off-resonance saturation pulse. Subsequently, MTR and cross-sectional areas (CSAs) of the sciatic nerve were analysed. In addition, detailed neurologic, physiotherapeutic and electrophysiologic examinations were conducted in all patients.

RESULTS

Sciatic nerve MTR and CSA reliably differentiated between healthy controls and SMA 3, 3a or 3b. MTR was lower in the SMA 3 (P < 0.0001), SMA 3a (P < 0.0001) and SMA 3b groups (P = 0.0020) than in respective controls. In patients with SMA 3, MTR correlated with all clinical scores, and arm nerve compound motor action potentials (CMAPs). CSA was lower in the SMA 3 (P < 0.0001), SMA 3a (P < 0.0001) and SMA 3b groups (P = 0.0006) than in controls, but did not correlate with clinical scores or electrophysiologic results.

CONCLUSIONS

Magnetization transfer ratio is a novel imaging marker that quantifies macromolecular nerve changes in SMA 3, and positively correlates with clinical scores and CMAPs.

Maturation of the Human Cerebral Cortex During Adolescence: Myelin or Dendritic Arbor?

Previous in vivo studies revealed robust age-related variations in structural properties of the human cerebral cortex during adolescence. Neurobiology underlying these maturational phenomena is largely unknown. Here we employ a virtual-histology approach to gain insights into processes associated with inter-regional variations in cortical microstructure and its maturation, as indexed by magnetization transfer ratio (MTR). Inter-regional variations in MTR correlate with inter-regional variations in expression of genes specific to pyramidal cells (CA1) and ependymal cells; enrichment analyses indicate involvement of these genes in dendritic growth. On the other hand, inter-regional variations in the change of MTR during adolescence correlate with inter-regional profiles of oligodendrocyte-specific gene expression. Complemented by a quantitative hypothetical model of the contribution of surfaces associated with dendritic arbor (1631 m2) and myelin (48 m2), these findings suggest that MTR signals are driven mainly by macromolecules associated with dendritic arbor while maturational changes in the MTR signal are associated with myelination.

Magnetization transfer ratio quantifies polyneuropathy in hereditary transthyretin amyloidosis

Objective

To quantify peripheral nerve lesions in symptomatic and asymptomatic hereditary transthyretin amyloidosis with polyneuropathy (ATTRv‐PNP) by analyzing the magnetization transfer ratio (MTR) of the sciatic nerve, and to test its potential as a novel biomarker for macromolecular changes.

Methods

Twenty‐five patients with symptomatic ATTRv‐PNP, 30 asymptomatic carriers of the mutant transthyretin gene (mutTTR), and 20 age‐/sex‐matched healthy controls prospectively underwent magnetization transfer contrast imaging at 3 Tesla. Two axial three‐dimensional gradient echo sequences with and without an off‐resonance saturation rapid frequency pulse were conducted at the right distal thigh. Sciatic nerve regions of interest were manually drawn on 10 consecutive axial slices in the images without off‐resonance saturation, and then transferred to the corresponding slices that were generated by the sequence with the off‐resonance saturation pulse. Subsequently, the MTR and cross‐sectional area (CSA) of the sciatic nerve were evaluated. Detailed neurologic and electrophysiologic examinations were conducted in all ATTRv‐PNP patients and mutTTR‐carriers.

Results

Sciatic nerve MTR and CSA reliably differentiated between ATTRv‐PNP, mutTTR‐carriers, and controls. MTR was lower in ATTRv‐PNP (26.4 ± 0.7; P < 0.0001) and in mutTTR‐carriers (32.6 ± 0.8; P = 0.0005) versus controls (39.4 ± 2.1), and was also lower in ATTRv‐PNP versus mutTTR‐carriers (P = 0.0009). MTR correlated negatively with the NIS‐LL and positively with CMAPs and SNAPs. CSA was higher in ATTRv‐PNP (34.3 ± 1.7 mm³) versus mutTTR‐carriers (26.0 ± 1.1 mm³; P = 0.0005) and versus controls (20.4 ± 1.2 mm³; P < 0.0001). CSA was also higher in mutTTR‐carriers versus controls.

Interpretation

MTR is a novel imaging marker that can quantify macromolecular changes in ATTRv‐PNP and differentiate between symptomatic ATTRv‐PNP and asymptomatic mutTTR‐carriers and correlates with electrophysiology.

Dipolar induced spin-lattice relaxation in the myelin sheath: A molecular dynamics study

Interactions between hydrogen protons of water molecules and macromolecules within the myelin sheath surrounding the axons are a major factor influencing the magnetic resonance (MR) contrast in white matter (WM) regions. In past decades, several studies have investigated the underlying effects and reported a wide range of R1 rates for the myelin associated compartments at different field strengths. However, it was also shown that the experimental quantification of the compartment-specific R1 rates is associated with large uncertainties. The current study therefore investigates the longitudinal relaxation rates within the myelin sheath using a molecular dynamic (MD) simulation. For this purpose, a realistic molecular model of the myelin sheath was employed to determine the dipole-dipole induced R1 relaxation rate of the hydrogen protons at clinically relevant field strengths. The results obtained clearly reflect the spatial heterogeneity of R1 with a increased relaxivity of myelin water due to a reduced molecular mobility near the membrane surface. Moreover, the calculated R1 rates for both myelin water and macromolecules are in excellent agreement with experimental findings from the literature at different field strengths.

Magnetic Resonance Imaging as a Biomarker in Rodent Peripheral Nerve Injury Models Reveals an Age-Related Impairment of Nerve Regeneration

Assessment of myelin integrity in peripheral nerve injuries and pathologies has largely been limited to post-mortem analysis owing to the difficulty in obtaining biopsies without affecting nerve function. This is further encumbered  by the small size of the tissue and its location. Therefore, the development of robust, non-invasive methods is highly attractive. In this study, we used magnetic resonance imaging (MRI) techniques, including magnetization transfer ratio (MTR), to longitudinally and non-invasively characterize both the sciatic nerve crush and lysolecithin (LCP) demyelination models of peripheral nerve injury in rodents. Electrophysiological, gene expression and histological assessments complemented the extensive MRI analyses in young and aged animals. In the nerve crush model, MTR analysis indicated a slower recovery in regions distal to the site of injury in aged animals, as well as incomplete recovery at six weeks post-crush when analyzing across the entire nerve surface. Similar regional impairments were also found in the LCP demyelination model. This research underlines the power of MTR for the study of peripheral nerve injury in small tissues such as the sciatic nerve of rodents and contributes new knowledge to the effect of aging on recovery after injury. A particular advantage of the approach is the translational potential to human neuropathies.

Two-dimensional magnetization-transfer - CPMG MRI reveals tract-specific signatures in fixed rat spinal cord

Multiexponential T₂ (MET₂) Relaxometry and Magnetization Transfer (MT) are among the most promising MRI-derived techniques for white matter (WM) characterization. Both techniques are shown to have histologically correlated sensitivity to myelin, but these correlations are not fully understood. Furthermore, MET₂ and MT report on different WM features, thus they can be considered specific to different (patho)physiological states. Two-dimensional studies potentially resolving interactions, such as those commonly used in NMR, have been rarely performed in this context. Here, we investigated how off-resonance irradiation affects different MET₂ components in fixed rat spinal cord white matter at 16.4 T. These 2D MT-MET₂ experiments reveal that MT affects both short and long T₂ components in a tract-specific fashion. The spatially distinct signal modulations enhanced contrast between microstructurally-distinct spinal cord tracts. Two hypotheses to explain these findings were proposed: either selective elimination of a short T₂ component through pre-saturation combines with intercompartmental water exchange effects occurring on the irradiation timescale; or, other macromolecular species that exist within the tissue - other than myelin - such as neurofilaments, may be involved in the apparent microstructural segregation of the spinal cord (SC) from MET₂. Though further investigation is required to elucidate the underlying mechanism, this phenomenon adds a new dimension for WM characterization.

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

On the accuracy and precision of PLANET for multiparametric MRI using phase-cycled bSSFP imaging

Purpose

In this work we demonstrate how sequence parameter settings influence the accuracy and precision in T₁, T₂, and off‐resonance maps obtained with the PLANET method for a single‐component signal model. In addition, the performance of the method for the particular case of a two‐component relaxation model for white matter tissue was assessed.

Methods

Numerical simulations were performed to investigate the influence of sequence parameter settings on the accuracy and precision in the estimated parameters for a single‐component model, as well as for a two‐component white matter model. Phantom and in vivo experiments were performed for validation. In addition, the effects of Gibbs ringing were investigated.

Results

By making a proper choice for sequence parameter settings, accurate and precise parameter estimation can be achieved for a single‐component signal model over a wide range of relaxation times at realistic SNR levels. Due to the presence of a second myelin‐related signal component in white matter, an underestimation of approximately 30% in T₁ and T₂ was observed, predicted by simulations and confirmed by measurements. Gibbs ringing artifacts correction improved the precision and accuracy of the parameter estimates.

Conclusion

For a single‐component signal model there is a broad “sweet spot” of sequence parameter combinations for which a high accuracy and precision in the parameter estimates is achieved over a wide range of relaxation times. For a multicomponent signal model, the single‐component PLANET reconstruction results in systematic errors in the parameter estimates as expected.

Conventional and quantitative MRI in a novel feline model of demyelination and endogenous remyelination

Background

The feeding of irradiated food to healthy adult cats results in widespread, noninflammatory demyelination of the central nervous system (CNS); a return to a normal diet results in endogenous remyelination with functional recovery. This recently discovered, reversible disease might provide a compelling clinical neuroimaging model system for the development and testing of myelin‐directed MRI methods as well as future remyelination therapies.

Purpose

Identify the noninvasive imaging characteristics of this new disease model and determine whether it features measurable changes on conventional and quantitative MRI.

Study Type

Pilot study.

Animal Model

Ten adult cats at various stages of demyelinating disease induced by an irradiated diet (35–55 kGy), and during recovery following a return to a normal diet.

Field Strength/Sequence

Conventional (T₂‐weighted) and quantitative (diffusion tensor, magnetization transfer) at 3T.

Assessment

MRI of the brain, optic nerves, and cervical spinal cord; a subset of diseased cats was euthanized for comparative histopathology.

Statistical Tests

Descriptive statistics.

Results

Disease produced T₂ prolongation, progressing from patchy to diffuse throughout most of the cerebral white matter (eventually involving U‐fibers) and spinal cord (primarily dorsal columns, reminiscent of subacute combined degeneration but without evidence of B12 deficiency). Magnetization transfer parameters decreased by 50–53% in cerebral white matter and by 25–30% in optic nerves and spinal cord dorsal columns. Fractional diffusion anisotropy decreased by up to 20% in pyramidal tracts, primarily driven by increased radial diffusivity consistent with axon preservation. Histopathology showed scattered myelin vacuolation of major white matter tracts as well as many thin myelin sheaths consistent with remyelination in the recovery phase, which was detectable on magnetization transfer imaging.

Data Conclusion

Feline irradiated diet‐induced demyelination features noninvasively imageable and quantifiable demyelination and remyelination of the CNS. It is therefore a compelling clinical neuroimaging model system.

Level of Evidence: 4

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2019;49:1304–1311.

Microstructural imaging in the spinal cord and validation strategies

In vivo histology using magnetic resonance imaging (MRI) is a newly emerging research field that aims to non-invasively characterize tissue microstructure. The implications of in vivo histology are many, from discovering novel biomarkers to studying human development, to providing tools for disease diagnosis and monitoring the effects of novel treatments on tissue. This review focuses on quantitative MRI (qMRI) techniques that are used to map spinal cord microstructure. Opening with a rationale for non-invasive imaging of the spinal cord, this article continues with a brief overview of the existing MRI techniques for axon and myelin imaging, followed by the specific challenges and potential solutions for acquiring and processing such data. The final part of this review focuses on histological validation, with suggested tissue preparation, acquisition and processing protocols for large-scale microscopy.

MRI gradient-echo phase contrast of the brain at ultra-short TE with off-resonance saturation

Larmor-frequency shift or image phase measured by gradient-echo sequences has provided a new source of MRI contrast. This contrast is being used to study both the structure and function of the brain. So far, phase images of the brain have been largely obtained at long echo times as maximum phase signal-to-noise ratio (SNR) is achieved at TE = T2* (∼40 ms at 3T). The structures of the brain, however, are compartmentalized and complex with a wide range of signal relaxation times. At such long TE, the short-T2 components are largely attenuated and contribute minimally to phase contrast. The purpose of this study was to determine whether proton gradient-echo images of the brain exhibit phase contrast at ultra-short TE (UTE). Our data showed that UTE images acquired at 7 T without off-resonance saturation do not contain significant phase contrast between gray and white matter. However, UTE images of the brain can attain strong phase contrast even at a nominal TE of 106 μs by using off-resonance RF saturation pulses, which provide direct saturation of ultra-short-T2 components and indirect saturation of longer-T2 components via magnetization transfer. In addition, phase contrast between gray and white matter acquired at UTE with off-resonance saturation is reversed compared to that of the long-T2 signals acquired at long TEs. This finding opens up a potential new way to manipulate image phase contrast of the brain. By accessing short and ultra-short-T2 species, MRI phase images may further improve the characterization of tissue microstructure in the brain.

Inferring brain tissue composition and microstructure via MR relaxometry

MRI relaxometry is sensitive to a variety of tissue characteristics in a complex manner, which makes it both attractive and challenging for characterizing tissue. This article reviews the most common water proton relaxometry measures, T₁, T₂, and T₂^*, and reports on their development and current potential to probe the composition and microstructure of brain tissue. The development of these relaxometry measures is challenged by the need for suitably accurate tissue models, as well as robust acquisition and analysis methodologies. MRI relaxometry has been established as a tool for characterizing neural tissue, particular with respect to myelination, and the potential for further development exists.

Diffusion-relaxation correlation spectroscopic imaging: A multidimensional approach for probing microstructure

PURPOSE

To propose and evaluate a novel multidimensional approach for imaging subvoxel tissue compartments called Diffusion-Relaxation Correlation Spectroscopic Imaging.

THEORY AND METHODS

Multiexponential modeling of MR diffusion or relaxation data is commonly used to infer the many different microscopic tissue compartments that contribute signal to macroscopic MR imaging voxels. However, multiexponential estimation is known to be difficult and ill-posed. Observing that this ill-posedness is theoretically reduced in higher dimensions, diffusion-relaxation correlation spectroscopic imaging uses a novel multidimensional imaging experiment that jointly encodes diffusion and relaxation information, and then uses a novel constrained reconstruction technique to generate a multidimensional diffusion-relaxation correlation spectrum for every voxel. The peaks of the multidimensional spectrum are expected to correspond to the distinct tissue microenvironments that are present within each macroscopic imaging voxel.

RESULTS

Using numerical simulations, experiment data from a custom-built phantom, and experiment data from a mouse model of traumatic spinal cord injury, diffusion-relaxation correlation spectroscopic imaging is demonstrated to provide substantially better multicompartment resolving power compared to conventional diffusion- and relaxation-based methods.

CONCLUSION

The diffusion-relaxation correlation spectroscopic imaging approach provides powerful new capabilities for resolving the different components of multicompartment tissue models, and can be leveraged to significantly expand the insights provided by MRI in studies of tissue microstructure. Magn Reson Med 78:2236-2249, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Lysophosphatidyl Choline Induced Demyelination in Rat Probed by Relaxation along a Fictitious Field in High Rank Rotating Frame

In this work a new MRI modality entitled Relaxation Along a Fictitious Field in the rotating frame of rank 4 (RAFF4) was evaluated in its ability to detect lower myelin content in lysophosphatidyl choline (LPC)-induced demyelinating lesions. The lesions were induced in two areas of the rat brain with either uniform or complex fiber orientations, i.e., in the corpus callosum (cc) and dorsal tegmental tract (dtg), respectively. RAFF4 showed excellent ability to detect demyelinated lesions and good correlation with myelin content in both brain areas. In comparison, diffusion tensor imaging metrices, fractional anisotropy, mean diffusivity and axonal and radial diffusivity, and magnetization transfer (MT) metrices, longitudinal relaxation during off-resonance irradiation and MT ratio, either failed to detect demyelination in dtg or showed lower correlation with myelin density quantified from gold chloride stained histological sections. Good specifity of RAFF4 to myelin was confirmed by its low correlation with cell density assesed from Nissl stained sections as well as its lack of sensitivity to pH changes in the physiological range as tested in heat denaturated bovine serum albumin phantoms. The excellent ability of RAFF4 to detect myelin content and its insensitivity to fiber orientation distribution, gliosis and pH, together with low specific absorption rate, demonstrates the promise of rotating frame of rank n (RAFFn) as a valuable MRI technique for non-invasive imaging of demyelinating lesions.

Fast, accurate 2D-MR relaxation exchange spectroscopy (REXSY): Beyond compressed sensing

Previously, we showed that compressive or compressed sensing (CS) can be used to reduce significantly the data required to obtain 2D-NMR relaxation and diffusion spectra when they are sparse or well localized. In some cases, an order of magnitude fewer uniformly sampled data were required to reconstruct 2D-MR spectra of comparable quality. Nonetheless, this acceleration may still not be sufficient to make 2D-MR spectroscopy practicable for many important applications, such as studying time-varying exchange processes in swelling gels or drying paints, in living tissue in response to various biological or biochemical challenges, and particularly for in vivo MRI applications. A recently introduced framework, marginal distributions constrained optimization (MADCO), tremendously accelerates such 2D acquisitions by using a priori obtained 1D marginal distribution as powerful constraints when 2D spectra are reconstructed. Here we exploit one important intrinsic property of the 2D-MR relaxation exchange spectra: the fact that the 1D marginal distributions of each 2D-MR relaxation exchange spectrum in both dimensions are equal and can be rapidly estimated from a single Carr-Purcell-Meiboom-Gill (CPMG) or inversion recovery prepared CPMG measurement. We extend the MADCO framework by further proposing to use the 1D marginal distributions to inform the subsequent 2D data-sampling scheme, concentrating measurements where spectral peaks are present and reducing them where they are not. In this way we achieve compression or acceleration that is an order of magnitude greater than that in our previous CS method while providing data in reconstructed 2D-MR spectral maps of comparable quality, demonstrated using several simulated and real 2D T2 - T2 experimental data. This method, which can be called "informed compressed sensing," is extendable to other 2D- and even ND-MR exchange spectroscopy.

Use of marginal distributions constrained optimization (MADCO) for accelerated 2D MRI relaxometry and diffusometry

Measuring multidimensional (e.g., 2D) relaxation spectra in NMR and MRI clinical applications is a holy grail of the porous media and biomedical MR communities. The main bottleneck is the inversion of Fredholm integrals of the first kind, an ill-conditioned problem requiring large amounts of data to stabilize a solution. We suggest a novel experimental design and processing framework to accelerate and improve the reconstruction of such 2D spectra that uses a priori information from the 1D projections of spectra, or marginal distributions. These 1D marginal distributions provide powerful constraints when 2D spectra are reconstructed, and their estimation requires an order of magnitude less data than a conventional 2D approach. This marginal distributions constrained optimization (MADCO) methodology is demonstrated here with a polyvinylpyrrolidone-water phantom that has 3 distinct peaks in the 2D D-T1 space. The stability, sensitivity to experimental parameters, and accuracy of this new approach are compared with conventional methods by serially subsampling the full data set. While the conventional, unconstrained approach performed poorly, the new method had proven to be highly accurate and robust, only requiring a fraction of the data. Additionally, synthetic T1-T2 data are presented to explore the effects of noise on the estimations, and the performance of the proposed method with a smooth and realistic 2D spectrum. The proposed framework is quite general and can also be used with a variety of 2D MRI experiments (D-T2,T1-T2,D-D, etc.), making these potentially feasible for preclinical and even clinical applications for the first time.

Differences in iron and manganese concentration may confound the measurement of myelin from R1 and R2 relaxation rates in studies of dysmyelination

A model of dysmyelination, the Long Evans Shaker (les) rat, was used to study the contribution of myelin to MR tissue properties in white matter. A large region of white matter was identified in the deep cerebellum and was used for measurements of the MR relaxation rate constants, R1  = 1/T1 and R2  = 1/T2 , at 7 T. In this study, R1 of the les deep cerebellar white matter was found to be 0.55 ± 0.08 s (-1) and R2 was found to be 15 ± 1 s(-1) , revealing significantly lower R1 and R2 in les white matter relative to wild-type (wt: R1  = 0.69 ± 0.05 s(-1) and R2  = 18 ± 1 s(-1) ). These deviated from the expected ΔR1 and ΔR2 values, given a complete lack of myelin in the les white matter, derived from the literature using values of myelin relaxivity, and we suspect that metals could play a significant role. The absolute concentrations of the paramagnetic transition metals iron (Fe) and manganese (Mn) were measured by a micro-synchrotron radiation X-ray fluorescence (μSRXRF) technique, with significantly greater Fe and Mn in les white matter than in wt (in units of μg [metal]/g [wet weight tissue]: les: Fe concentration,19 ± 1; Mn concentration, 0.71 ± 0.04; wt: Fe concentration,10 ± 1; Mn concentration, 0.47 ± 0.04). These changes in Fe and Mn could explain the deviations in R1 and R2 from the expected values in white matter. Although it was found that the influence of myelin still dominates R1 and R2 in wt rats, there were non-negligible changes in the contribution of the metals to relaxation. Although there are already problems with the estimation of myelin from R1 and R2 changes in disease models with pathology that also affects the relaxation rate constants, this study points to a specific pitfall in the estimation of changes in myelin in diseases or models with disrupted concentrations of paramagnetic transition metals. Copyright © 2016 John Wiley & Sons, Ltd.

Modeling T(1) and T(2) relaxation in bovine white matter

The fundamental basis of T1 and T2 contrast in brain MRI is not well understood; recent literature contains conflicting views on the nature of relaxation in white matter (WM). We investigated the effects of inversion pulse bandwidth on measurements of T1 and T2 in WM. Hybrid inversion-recovery/Carr-Purcell-Meiboom-Gill experiments with broad or narrow bandwidth inversion pulses were applied to bovine WM in vitro. Data were analysed with the commonly used 1D-non-negative least squares (NNLS) algorithm, a 2D-NNLS algorithm, and a four-pool model which was based upon microscopically distinguishable WM compartments (myelin non-aqueous protons, myelin water, non-myelin non-aqueous protons and intra/extracellular water) and incorporated magnetization exchange between adjacent compartments. 1D-NNLS showed that different T2 components had different T1 behaviours and yielded dissimilar results for the two inversion conditions. 2D-NNLS revealed significantly more complicated T1/T2 distributions for narrow bandwidth than for broad bandwidth inversion pulses. The four-pool model fits allow physical interpretation of the parameters, fit better than the NNLS techniques, and fits results from both inversion conditions using the same parameters. The results demonstrate that exchange cannot be neglected when analysing experimental inversion recovery data from WM, in part because it can introduce exponential components having negative amplitude coefficients that cannot be correctly modeled with nonnegative fitting techniques. While assignment of an individual T1 to one particular pool is not possible, the results suggest that under carefully controlled experimental conditions the amplitude of an apparent short T1 component might be used to quantify myelin water.

In vivo mapping of human spinal cord microstructure at 300mT/m

The ability to characterize white matter microstructure non-invasively has important applications for the diagnosis and follow-up of several neurological diseases. There exists a family of diffusion MRI techniques, such as AxCaliber, that provide indices of axon microstructure, such as axon diameter and density. However, to obtain accurate measurements of axons with small diameters (<5μm), these techniques require strong gradients, i.e. an order of magnitude higher than the 40-80mT/m currently available in clinical systems. In this study we acquired AxCaliber diffusion data at a variety of different q-values and diffusion times in the spinal cord of five healthy subjects using a 300mT/m whole body gradient system. Acquisition and processing were optimized using state-of-the-art methods (e.g., 64-channel coil, template-based analysis). Results consistently show an average axon diameter of 4.5+/-1.1μm in the spinal cord white matter. Diameters ranged from 3.0μm (gracilis) to 5.9μm (spinocerebellar tracts). Values were similar across laterality (left-right), but statistically different across spinal cord pathways (p<10(-5)). The observed trends are similar to those observed in animal histology. This study shows, for the first time, in vivo mapping of axon diameter in the spinal cord at 300mT/m, thus creating opportunities for applications in spinal cord diseases.

MRI relaxation in the presence of fictitious fields correlates with myelin content in normal rat brain

PURPOSE

Brain myelin plays an important role in normal brain function. Demyelination is involved in many degenerative brain diseases, thus quantitative imaging of myelin has been under active investigation. In previous work, we demonstrated the capability of the method known as Relaxation Along a Fictitious Field (RAFF) in the rotating frame of rank n (RAFFn) to provide image contrast between white and gray matter in human and rat brains. Here, we provide evidence pointing to myelin being the major source of this contrast.

METHODS

RAFFn relaxation time constant (TRAFFn) was mapped in rat brain ex vivo. TRAFFn was quantified in 12 different brain areas. TRAFFn values were compared with multiple other MRI metrics (T1, T2 , continuous wave T1ρ, adiabatic T1ρ and T2ρ, magnetization transfer ratio), and with histologic measurements of cell density, myelin and iron content.

RESULTS

Highest contrast between white and grey matter was obtained with TRAFFn in the rotating frames of ranks n = 4 and 5. TRAFFn values correlated strongly with myelin content, whereas no associations between TRAFFn and iron content or cell density were found.

CONCLUSION

TRAFFn with n = 4 or 5 provides a high sensitivity for selective myelin mapping in the rat brain.

MRI-based myelin water imaging: A technical review

Multiexponential T2 relaxation time measurement in the central nervous system shows a component that originates from water trapped between the lipid bilayers of myelin. This myelin water component is of significant interest as it provides a myelin-specific MRI signal of value in assessing myelin changes in cerebral white matter in vivo. In this article, the various acquisition and analysis strategies proposed to date for myelin water imaging are reviewed and research conducted into their validity and clinical applicability is presented. Comparisons between the imaging methods are made with a discussion regarding potential difficulties and model limitations.

Investigation of outer cortical magnetisation transfer ratio abnormalities in multiple sclerosis clinical subgroups

BACKGROUND

Pathological abnormalities including demyelination and neuronal loss are reported in the outer cortex in multiple sclerosis (MS).

OBJECTIVE

We investigated for in vivo evidence of outer cortical abnormalities by measuring the magnetisation transfer ratio (MTR) in MS patients of different subgroups.

METHODS

Forty-four relapsing-remitting (RR) (mean age 41.9 years, median Expanded Disability Status Scale (EDSS) 2.0), 25 secondary progressive (SP) (54.1 years, EDSS 6.5) and 19 primary progressive (PP) (53.1 years, EDSS 6.0) MS patients and 35 healthy control subjects (mean age 39.2 years) were studied. Three-dimensional (3D) 1×1×1mm(3) T1-weighted images and MTR data were acquired. The cortex was segmented, then subdivided into outer and inner bands, and MTR values were calculated for each band.

RESULTS

In a pairwise analysis, mean outer cortical MTR was lower than mean inner cortical MTR in all MS groups and controls (p<0.001). Compared with controls, outer cortical MTR was decreased in SPMS (p<0.001) and RRMS (p<0.01), but not PPMS. Outer cortical MTR was lower in SPMS than PPMS (p<0.01) and RRMS (p<0.01).

CONCLUSIONS

Lower outer than inner cortical MTR in healthy controls may reflect differences in myelin content. The lowest outer cortical MTR was seen in SPMS and is consistent with more extensive outer cortical (including subpial) pathology, such as demyelination and neuronal loss, as observed in post-mortem studies of SPMS patients.

Improved intravoxel incoherent motion analysis of diffusion weighted imaging by data driven Bayesian modeling

In addition to the diffusion coefficient, fitting the intravoxel incoherent motion model to multiple b-value diffusion-weighted MR data gives pseudo-diffusion measures associated with rapid signal attenuation at low b-values that are of use in the assessment of a number of pathologies. When summary measures are required, such as the average parameter for a region of interest, least-squares based methods give adequate estimation accuracy. However, using least-squares methods for pixel-wise fitting typically gives noisy estimates, especially for the pseudo-diffusion parameters, which limits the applicability of the approach for assessing spatial features and heterogeneity. In this article, a Bayesian approach using a shrinkage prior model is proposed and is shown to substantially reduce estimation uncertainty so that spatial features in the parameters maps are more clearly apparent. The Bayesian approach has no user-defined parameters, so measures of parameter variation (heterogeneity) over regions of interest are determined by the data alone, whereas it is shown that for the least-squares estimates, measures of variation are essentially determined by user-defined constraints on the parameters. Use of a Bayesian shrinkage prior approach is, therefore, recommended for intravoxel incoherent motion modeling.

Solving 2D Fredholm Integral from Incomplete Measurements Using Compressive Sensing

We present an algorithm to solve the two-dimensional Fredholm integral of the first kind with tensor product structure from a limited number of measurements, with the goal of using this method to speed up nuclear magnetic resonance spectroscopy. This is done by incorporating compressive sensing-type arguments to fill in missing measurements, using a priori knowledge of the structure of the data. In the first step we recover a compressed data matrix from measurements that form a tight frame, and establish that these measurements satisfy the restricted isometry property. Recovery can be done from as few as 10% of the total measurements. In the second and third steps, we solve the zeroth-order regularization minimization problem using the Venkataramanan-Song-Hürlimann algorithm. We demonstrate the performance of this algorithm on simulated data and show that our approach is a realistic approach to speeding up the data acquisition.

The role of tissue microstructure and water exchange in biophysical modelling of diffusion in white matter

Biophysical models that describe the outcome of white matter diffusion MRI experiments have various degrees of complexity. While the simplest models assume equal-sized and parallel axons, more elaborate ones may include distributions of axon diameters and axonal orientation dispersions. These microstructural features can be inferred from diffusion-weighted signal attenuation curves by solving an inverse problem, validated in several Monte Carlo simulation studies. Model development has been paralleled by microscopy studies of the microstructure of excised and fixed nerves, confirming that axon diameter estimates from diffusion measurements agree with those from microscopy. However, results obtained in vivo are less conclusive. For example, the amount of slowly diffusing water is lower than expected, and the diffusion-encoded signal is apparently insensitive to diffusion time variations, contrary to what may be expected. Recent understandings of the resolution limit in diffusion MRI, the rate of water exchange, and the presence of microscopic axonal undulation and axonal orientation dispersions may, however, explain such apparent contradictions. Knowledge of the effects of biophysical mechanisms on water diffusion in tissue can be used to predict the outcome of diffusion tensor imaging (DTI) and of diffusion kurtosis imaging (DKI) studies. Alterations of DTI or DKI parameters found in studies of pathologies such as ischemic stroke can thus be compared with those predicted by modelling. Observations in agreement with the predictions strengthen the credibility of biophysical models; those in disagreement could provide clues of how to improve them. DKI is particularly suited for this purpose; it is performed using higher b-values than DTI, and thus carries more information about the tissue microstructure. The purpose of this review is to provide an update on the current understanding of how various properties of the tissue microstructure and the rate of water exchange between microenvironments are reflected in diffusion MRI measurements. We focus on the use of biophysical models for extracting tissue-specific parameters from data obtained with single PGSE sequences on clinical MRI scanners, but results obtained with animal MRI scanners are also considered. While modelling of white matter is the central theme, experiments on model systems that highlight important aspects of the biophysical models are also reviewed.

Characterizing inter-compartmental water exchange in myelinated tissue using relaxation exchange spectroscopy

PURPOSE

To investigate inter-compartmental water exchange in two model myelinated tissues ex vivo using relaxation exchange spectroscopy.

METHODS

Building upon a previously developed theoretical framework, a three-compartment (myelin, intra-axonal, and extra-axonal water) model of the inversion-recovery prepared relaxation exchange spectroscopy signal was applied in excised rat optic nerve and frog sciatic nerve samples to estimate the water residence time constants in myelin (τmyelin ).

RESULTS

In the rat optic nerve samples, τmyelin = 138 ± 15 ms (mean ± standard deviation) was estimated. In sciatic nerve, which possesses thicker myelin sheaths than optic nerve, a much longer τmyelin = 2046 ± 140 ms was observed.

CONCLUSION

Consistent with previous studies in rat spinal cord, the extrapolation of exchange rates in optic nerve to in vivo conditions indicates that τmyelin < 100 ms. This suggests that there is a significant effect of inter-compartmental water exchange on the transverse relaxation of water protons in white matter. The much longer τmyelin values in sciatic nerve supports the postulate that the inter-compartmental water exchange rate is mediated by myelin thickness. Together, these findings point to the potential for MRI methods to probe variations in myelin thickness in white matter.

Quantitative MR imaging of two-pool magnetization transfer model parameters in myelin mutant shaking pup

Magnetization transfer (MT) imaging quantitatively assesses cerebral white matter disease through its sensitivity to macromolecule-bound protons including those associated with myelin proteins and lipid bilayers. However, traditional MT contrast measured by the MT ratio (MTR) lacks pathologic specificity as demyelination, axon loss, inflammation and edema all impact MTR, directly and/or indirectly through multiple covariances among imaging parameters (particularly MTR with T(1)) and tissue features (e.g. axon loss with demyelination). In this study, more complex modeling of MT phenomena ("quantitative" MT or qMT) was applied to a less complex disease model (the myelin mutant shaking [sh] pup, featuring hypomyelination but neither inflammation nor axon loss) in order to eliminate the covariances on both sides of the MR-pathology "equation" and characterize these important relationships free from the usual confounds. qMT measurements were acquired longitudinally in 6 sh pups and 4 age-matched controls ranging from 3 to 21 months of age and compared with histology. The qMT parameter, bound pool fraction (f), was the most distinctive between diseased and control animals; both f and longitudinal relaxation rate R(1) tracked myelination with normal aging, whereas MTR did not--presumably owing to counterbalancing MT and R(1) effects. qMT imaging provides a more accurate and potentially more specific non-invasive tissue characterization.

Origins of the ultrashort-T2 1H NMR signals in myelinated nerve: a direct measure of myelin content?

Recently developed MRI techniques have enabled clinical imaging of short-lived (1)H NMR signals with T(2) < 1 ms. Using these techniques, novel signal enhancement has been observed in myelinated tissues, although the source of this enhancement has not been identified. Herein, we report studies of the nature and origins of ultrashort T(2) (uT(2)) signals (50 μs < T(2) < 1 ms) from amphibian and mammalian myelinated nerves. NMR measurements and comparisons with myelin phantoms and expected myelin components indicate that these uT(2) signals arise predominantly from methylene (1)H on/in the myelin membranes, which suggests that direct measurement of uT(2) signals can be used as a new means for quantitative myelin mapping.

Diffusion effects on longitudinal relaxation in poorly mixed compartments

Diffusion of spins between physical or virtual, communicating compartments having different states of longitudinal magnetization leads to diffusion-driven longitudinal relaxation. Herein, in two model systems, the effects of diffusion-driven longitudinal relaxation are explored experimentally and analyzed quantitatively. In the first case, longitudinal relaxation in a single slice of a water phantom is monitored spectroscopically as a function of slice thickness. In the second case, mimicking vascular flow/diffusion effects, longitudinal relaxation is monitored in a two-compartment, semi-permeable fiber phantom. In both cases, apparent longitudinal relaxation, though clearly multi-exponential, is well-modeled as bi-exponential.

Compartment-specific enhancement of white matter and nerve ex vivo using chromium

Chromium--Cr(VI) in the form of potassium dichromate--has been shown to specifically enhance white matter signal. The proposed mechanism for this enhancement is reduction of diamagnetic Cr(VI) to paramagnetic chromium species by oxidizable myelin lipids. The purpose of the study herein was to better understand the microanatomical basis of this enhancement (i.e., the relative enhancement of myelin, intra-axonal, and extra-axonal water). Toward this end, integrated T(1)-T(2) measurements were performed in potassium dichromate loaded (hereafter referred to as chromated) rat brains, rat optic nerve samples, and frog sciatic nerve samples ex vivo. In control optic nerve and white matter, two T(1)-T(2) components were resolved, representing myelin and nonmyelin water (intra- and extra-axonal water). Following chromation, three T(1)-T(2) components were resolved in these same tissues. Results from similar measurements in sciatic nerve-all three components are resolvable in control and chromated samples-and quantitative histologic analysis suggest that this additional T(1)-T(2) component is due to a splitting of the nonmyelin water component into intra- and extra-axonal water components. This compartment-specific enhancement may provide unique contrast for MR histology, as well as allow one to probe the compartmental basis of various contrast mechanisms in neural tissue.

AnalyzeNNLS: magnetic resonance multiexponential decay image analysis

Exponential decays are fundamental to magnetic resonance imaging, yet adequately sampling and analyzing multiexponential decays is rarely attempted. The advantage of multiexponential analysis is the quantification of sub-voxel structure caused by water compartmentalization, with application as a non-invasive imaging biomarker for myelin. We have developed AnalyzeNNLS, software designed specifically for multiexponential decay image analysis that has a user-friendly graphical user interface and can analyze data from many MR manufacturers. AnalyzeNNLS is a simple, platform independent analysis tool that was created using the extensive mathematical and visualization libraries in Matlab, and released as open source code allowing scientists to evaluate, scrutinize, improve, and expand.

Multiexponential T2, magnetization transfer, and quantitative histology in white matter tracts of rat spinal cord

Quantitative MRI measures of multiexponential T(2) relaxation and magnetization transfer were acquired from six samples of excised and fixed rat spinal cord and compared with quantitative histology. MRI and histology data were analyzed from six white matter tracts, each of which possessed unique microanatomic characteristics (axon diameter and myelin thickness, in particular) but a relatively constant volume fraction of myelin. The results indicated that multiexponential T(2) relaxation characteristics varied substantially with variation of microanatomy, while the magnetization transfer characteristics remained close to constant. The most-often-cited multiexponential T(2) relaxation metric, myelin water fraction, varied by almost a factor of 2 between two regions with myelin volume fractions that differed by only approximately 12%. Based on the quantitative histology, the proposed explanation for this variation was intercompartmental water exchange, which caused the underestimation of myelin water fraction and T(2) values and is, presumably, a greater factor in white matter regions where axons are small and myelin is thin. In contrast to the multiexponential T(2) relaxation observations, magnetization transfer metrics were relatively constant across white matter tracts and concluded to be relatively insensitive to intercompartmental water exchange.

Sciatic nerve injury model in the axolotl: functional, electrophysiological, and radiographic outcomes

OBJECT

The 2 aims of this study were as follows: 1) to establish outcome measures of nerve regeneration in an axolotl model of peripheral nerve injury; and 2) to define the timing and completeness of reinnervation in the axolotl following different types of sciatic nerve injury.

METHODS

The sciatic nerves in 36 axolotls were exposed bilaterally in 3 groups containing 12 animals each: Group 1, left side sham, right side crush; Group 2, left side sham, right side nerve resected and proximal stump buried; and Group 3 left side cut and sutured, right side cut and sutured with tibial and peroneal divisions reversed. Outcome measures included the following: 1) an axolotl sciatic functional index (ASFI) derived from video swim analysis; 2) motor latencies; and 3) MR imaging evaluation of nerve and muscle edema.

RESULTS

For crush injuries, the ASFI returned to baseline by 2 weeks, as did MR imaging parameters and motor latencies. For buried nerves, the ASFI returned to 20% below baseline by 8 weeks, with motor evoked potentials present. On MR imaging, nerve edema peaked at 3 days postintervention and gradually normalized over 12 weeks, whereas muscle denervation was present until a gradual decrease was seen between 4 and 12 weeks. For cut nerves, the ASFI returned to 20% below baseline by Week 4, where it plateaued. Motor evoked potentials were observed at 2-4 weeks, but with an increased latency until Week 6, and MR imaging analysis revealed muscle denervation for 4 weeks.

CONCLUSIONS

Multiple outcome measures in which an axolotl model of peripheral nerve injury is used have been established. Based on historical controls, recovery after nerve injury appears to occur earlier and is more complete than in rodents. Further investigation using this model as a successful "blueprint" for nerve regeneration in humans is warranted.

The MT pool size ratio and the DTI radial diffusivity may reflect the myelination in shiverer and control mice

A quantitative magnetization transfer (qMT) technique was employed to quantify the ratio of the sizes of the bound and free water proton pools in ex vivo mouse brains. The goal was to determine the pool size ratio sensitivity to myelin. Fixed brains from both shiverer mice and control littermates were imaged. The pool size ratio in the corpus callosum of shiverer mice was substantially lower than that in the control mice, while there was no distinguishable difference in the pool size ratio in the gray matter. These results correlate with diffusion tensor imaging (DTI) derived radial diffusivity which previously was shown to reflect myelin integrity in this animal model. Histological study reveals the presence of myelin in control mice white matter and the absence of myelin in shiverer mice white matter, supporting the qMT and DTI results. Our findings support the view that qMT may be used for estimating myelin integrity.

Gleaning multicomponent T1 and T2 information from steady-state imaging data

The driven-equilibrium single-pulse observation of T(1) (DESPOT1) and T(2) (DESPOT2) are rapid, accurate, and precise methods for voxelwise determination of the longitudinal and transverse relaxation times. A limitation of the methods, however, is the inherent assumption of single-component relaxation. In a variety of biological tissues, in particular human white matter (WM) and gray matter (GM), the relaxation has been shown to be more completely characterized by a summation of two or more relaxation components, or species, each believed to be associated with unique microanatomical domains or water pools. Unfortunately, characterization of these components on a voxelwise, whole-brain basis has traditionally been hindered by impractical acquisition times. In this work we extend the conventional DESPOT1 and DESPOT2 approaches to include multicomponent relaxation analysis. Following numerical analysis of the new technique, renamed multicomponent driven equilibrium single pulse observation of T(1)/T(2) (mcDESPOT), whole-brain multicomponent T(1) and T(2) quantification is demonstrated in vivo with clinically realistic times of between 16 and 30 min. Results obtained from four healthy individuals and two primary progressive multiple sclerosis (MS) patients demonstrate the future potential of the approach for identifying and assessing tissue changes associated with several neurodegenerative conditions, in particular those associated with WM.

Magnetization transfer induced biexponential longitudinal relaxation

Longitudinal relaxation of brain water (1)H magnetization in mammalian brain in vivo is typically analyzed on a per-voxel basis using a monoexponential model, thereby assigning a single relaxation time constant to all (1)H magnetization within a given voxel. This approach was tested by obtaining inversion recovery (IR) data from gray matter of rats at 64 exponentially spaced recovery times. Using Bayesian probability for model selection, brain water data were best represented by a biexponential function characterized by fast and slow relaxation components. At 4.7T, the amplitude fraction of the rapidly relaxing component is 3.4% +/- 0.7% with a rate constant of 44 +/- 12 s(-1) (mean +/- SD; 174 voxels from four rats). The rate constant of the slow relaxing component is 0.66 +/- 0.04 s(-1). At 11.7T, the corresponding values are 6.9% +/- 0.9%, 19 +/- 5 s(-1), and 0.48 +/- 0.02 s(-1) (151 voxels from four rats). Several putative mechanisms for biexponential relaxation behavior were evaluated, and magnetization transfer (MT) between bulk water protons and nonaqueous protons was determined to be the source of biexponential longitudinal relaxation. MR methods requiring accurate quantification of longitudinal relaxation may need to take this effect explicitly into account.