Insights and improvements in correspondence between axonal volume fraction measured with diffusion-weighted MRI and electron microscopy.

Biophysical diffusion-weighted imaging (DWI) models are increasingly used in neuroscience to estimate the axonal water fraction ( f AW ), which in turn is key for noninvasive estimation of the axonal volume fraction ( f A ). These models require thorough validation by comparison with a reference method, for example, electron microscopy (EM). While EM studies often neglect the unmyelinated axons and solely report the fraction of myelinated axons, in DWI both myelinated and unmyelinated axons contribute to the DWI signal. However, DWI models often include simplifications, for example, the neglect of differences in the compartmental relaxation times or fixed diffusivities, which in turn might affect the estimation of f AW . We investigate whether linear calibration parameters (scaling and offset) can improve the comparability between EM- and DWI-based metrics of f A . To this end, we (a) used six DWI models based on the so-called standard model of white matter (WM), including two models with fixed compartmental diffusivities (e.g., neurite orientation dispersion and density imaging, NODDI) and four models that fitted the compartmental diffusivities (e.g., white matter tract integrity, WMTI), and (b) used a multimodal data set including ex vivo diffusion DWI and EM data in mice with a broad dynamic range of fibre volume metrics. We demonstrated that the offset is associated with the volume fraction of unmyelinated axons and the scaling factor is associated with different compartmental T 2 and can substantially enhance the comparability between EM- and DWI-based metrics of f A . We found that DWI models that fitted compartmental diffusivities provided the most accurate estimates of the EM-based f A . Finally, we introduced a more efficient hybrid calibration approach, where only the offset is estimated but the scaling is fixed to a theoretically predicted value. Using this approach, a similar one-to-one correspondence to EM was achieved for WMTI. The method presented can pave the way for use of validated DWI-based models in clinical research and neuroscience.

In Vivo Assessment of Bone Quality Without X-rays.

PURPOSE OF REVIEW: This review summarizes recent advances in the assessment of bone quality using non-X-ray techniques. RECENT FINDINGS: Quantitative ultrasound (QUS) provides multiple measurements of bone characteristics based on the propagation of sound through bone, the attenuation of that sound, and different processing techniques. QUS parameters and model predictions based on backscattered signals can discriminate non-fracture from fracture cases with accuracy comparable to standard bone mineral density (BMD). With advances in magnetic resonance imaging (MRI), bound water and pore water, or a porosity index, can be quantified in several long bones in vivo. Since such imaging-derived measurements correlate with the fracture resistance of bone, they potentially provide new BMD-independent predictors of fracture risk. While numerous measurements of mineral, organic matrix, and bound water by Raman spectroscopy correlate with the strength and toughness of cortical bone, the clinical assessment of person's bone quality using spatially offset Raman spectroscopy (SORS) requires advanced spectral processing techniques that minimize contaminating signals from fat, skin, and blood. Limiting exposure of patients to ionizing radiation, QUS, MRI, and SORS has the potential to improve the assessment of fracture risk and track changes of new therapies that target bone matrix and micro-structure.

Fast bound and pore water mapping of cortical bone with arbitrary slice oriented two-dimensional ultra-short echo time.

PURPOSE: Extend fast, two-dimensional (2D) methods of bound and pore water mapping in bone to arbitrary slice orientation. METHODS: To correct for slice profile artifacts caused by gradient errors of half pulse 2D ultra-short echo time (UTE), we developed a library of predistorted gradient waveforms that can be used to interpolate optimized gradient waveforms for 2D UTE slice selection. We also developed a method to estimate and correct for a bulk phase difference between the two half pulse excitations used for 2D UTE signal excitation. Bound water images were acquired in three healthy subjects with adiabatic inversion recovery prepared 2D UTE, while pore water images were acquired after short-T2 signals were suppressed with double adiabatic inversion recovery preparation. The repeatability of bound and pore water imaging with 2D UTE was tested by repeating acquisitions after repositioning. RESULTS: The library-based interpolation of optimized slice select gradient waveforms combined with the method to estimate bulk phase between two excitations provided compact slice profiles for half pulse excited 2D UTE. Quantitative bound and pore water values were highly repeatable-the pooled SD of bound water across all three subjects was 0.38 mol 1 $$ {}^1 $$ H/L, while pooled SD of pore water was 0.30 mol 1 $$ {}^1 $$ H/L. CONCLUSION: Fast, quantitative, 2D UTE-based bound and pore water images can be acquired at arbitrary oblique orientations after correcting for errors in the slice select gradient waveform and bulk phase shift between the two half acquisitions.

T1 relaxation of bound and pore water in cortical bone.

MRI measures of bound and/or pore water concentration in cortical bone offer potential diagnostics of bone fracture risk. The transverse relaxation characteristics of both bound and pore water are relatively well understood and have been used to design clinical MRI pulse sequences to image each water pool quantitatively. However, these methods are also sensitive to longitudinal relaxation characteristics, which have been less well studied. Here, spectroscopic relaxometry measurements of 31 human cortical bone specimens provided a more detailed picture of T 1 of both bound and pore water. The results included mean, standard deviation, and range of T 1 spectra from both bound and pore water, as well as novel presentations of the 2D T 1 - T 2 distribution of pore water. Importantly, for each sample the pore water T 1 spectrum was found to span more than one order of magnitude and varied substantially across the 31 sample studies. Because many existing methods assume pore water T 1 to be mono-exponential and constant across individuals, the results were used to compute the potential effect neglecting this intra- and intersample T 1 variation on accurate MRI measurement of both bound and pore water concentrations. The greatest effect was found for adiabatic inversion recovery (AIR) based measurements of bound water concentration, which showed an average of 8.8% and as much as 37% error when using a common mono-exponential assumption of pore water T 1 . Despite these errors, the simulated AIR measurements were still moderately well correlated with the bound water concentrations derived from the spectroscopic data.

A hybrid numeric-analytic solution for pulsed CEST.

Chemical exchange saturation transfer (CEST) methods measure the effect of magnetization exchange between solutes and water. While CEST methods are often implemented using a train of off-resonant shaped RF pulses, they are typically analyzed as if the irradiation were continuous. This approximation does not account for exchange of rotated magnetization, unique to pulsed irradiation and exploited by chemical exchange rotation transfer methods. In this work, we derive and test an analytic solution for the steady-state water signal under pulsed irradiation by extending a previous work to include the effects of pulse shape. The solution is largely accurate at all offsets, but this accuracy diminishes at higher exchange rates and when applying pulse shapes with large root-mean-squared to mean ratios (such as multi-lobe sinc pulses).

A simple estimate of axon size with diffusion MRI.

Noninvasive estimation of mean axon diameter presents a new opportunity to explore white matter plasticity, development, and pathology. Several diffusion-weighted MRI (DW-MRI) methods have been proposed to measure the average axon diameter in white matter, but they typically require many diffusion encoding measurements and complicated mathematical models to fit the signal to multiple tissue compartments, including intra- and extra-axonal spaces. Here, Monte Carlo simulations uncovered a straightforward DW-MRI metric of axon diameter: the change in radial apparent diffusion coefficient estimated at different effective diffusion times, ΔD⊥. Simulations indicated that this metric increases monotonically within a relevant range of effective mean axon diameter while being insensitive to changes in extra-axonal volume fraction, axon diameter distribution, g-ratio, and influence of myelin water. Also, a monotonic relationship was found to exist for signals coming from both intra- and extra-axonal compartments. The slope in ΔD⊥ with effective axon diameter increased with the difference in diffusion time of both oscillating and pulsed gradient diffusion sequences.

Glans inflation morphology and female cloaca copulatory interactions of the male American alligator phallus†.

The phallic glans of the American alligator (Alligator mississippiensis) is the distal termination of the semen-conducting sulcus spermaticus and during copulation has the closest, most intimate mechanical interactions with the female urodeum, the middle cloacal chamber that contains the opening to the vaginal passages and oviducts. However, the details of this interface leading to insemination and gamete uptake are unclear. Here, we: (1) histologically characterize the underlying tissue types and morphologically quantify the shape changes associated with glans inflation into the copulatory conformation, (2) digitally reconstruct from MRI the 3D shape of functional tissue compartments, and (3) diffusible iodine-based contrast-enhanced computed tomography image the copulatory fit between male phallus and female cloaca. We discuss these results in relation to tissue type material properties, the transfer on intromittent forces, establishing potential copulatory lock, inflated glans volume scaling with body mass/length, the mechanics of semen targeting and insemination, and potential female cryptic choice impacting multiple clutch paternity. In part, this study further clarifies the phallic morphological variation observed among crocodylians and begins to investigate the role(s) these divergent male forms play during copulation interacting with female cloacal forms to increase reproductive success.

Chemical exchange rotation transfer imaging of phosphocreatine in muscle.

In chemical exchange saturation transfer (CEST) imaging, the signal at 2.6 ppm from the water resonance in muscle has been assigned to phosphocreatine (PCr). However, this signal has limited specificity for PCr since the signal is also sensitive to exchange with protein and macromolecular protons when using some conventional quantification methods, and will vary with changes in the water longitudinal relaxation rate. Correcting for these effects while maintaining reasonable acquisition times is challenging. As an alternative approach to overcome these problems, here we evaluate chemical exchange rotation transfer (CERT) imaging of PCr in muscle at 9.4 T. Specifically, the CERT metric, AREXdouble,cpw at 2.6 ppm, was measured in solutions containing the main muscle metabolites, in tissue homogenates with controlled PCr content, and in vivo in rat leg muscles. PCr dominates CERT metrics around 2.6 ppm (although with nontrivial confounding baseline contributions), indicating that CERT is well-suited to PCr specific imaging, and has the added benefit of requiring a relatively small number of acquisitions.

Diffusion time dependency along the human corpus callosum and exploration of age and sex differences as assessed by oscillating gradient spin-echo diffusion tensor imaging.

Conventional diffusion imaging uses pulsed gradient spin echo (PGSE) waveforms with diffusion times of tens of milliseconds (ms) to infer differences of white matter microstructure. The combined use of these long diffusion times with short diffusion times (<10 â€‹ms) enabled by oscillating gradient spin echo (OGSE) waveforms can enable more sensitivity to changes of restrictive boundaries on the scale of white matter microstructure (e.g. membranes reflecting the axon diameters). Here, PGSE and OGSE images were acquired at 4.7 â€‹T from 20 healthy volunteers aged 20-73 years (10 males). Mean, radial, and axial diffusivity, as well as fractional anisotropy were calculated in the genu, body and splenium of the corpus callosum (CC). Monte Carlo simulations were also conducted to examine the relationship of intra- and extra-axonal radial diffusivity with diffusion time over a range of axon diameters and distributions. The results showed elevated diffusivities with OGSE relative to PGSE in the genu and splenium (but not the body) in both males and females, but the OGSE-PGSE difference was greater in the genu for males. Females showed positive correlations of OGSE-PGSE diffusivity difference with age across the CC, whereas there were no such age correlations in males. Simulations of radial diffusion demonstrated that for axon sizes in human brain both OGSE and PGSE diffusivities were dominated by extra-axonal water, but the OGSE-PGSE difference nonetheless increased with area-weighted outer-axon diameter. Therefore, the lack of OGSE-PGSE difference in the body is not entirely consistent with literature that suggests it is composed predominantly of axons with large diameter. The greater OGSE-PGSE difference in the genu of males could reflect larger axon diameters than females. The OGSE-PGSE difference correlation with age in females could reflect loss of smaller axons at older ages. The use of OGSE with short diffusion times to sample the microstructural scale of restriction implies regional differences of axon diameters along the corpus callosum with preliminary results suggesting a dependence on age and sex.

Slice-selective extended phase graphs in gradient-crushed, transient-state free precession sequences: An application to MR fingerprinting.

PURPOSE: Slice-selective, gradient-crushed, transient-state sequences such as those used in MR fingerprinting (MRF) relaxometry are sensitive to slice profile effects. Whereas balanced steady-state free precession MRF profile effects have been studied, less attention has been given to gradient-crushed MRF forms. Extensions of the extended phase graph (EPG) formalism, called slice-selective EPG (ssEPG), are proposed that model slice profile effects. THEORY AND METHODS: The hard-pulse approximation to slice-selective excitation in the spatial domain is reformulated in k-space. Excitation is modeled by standard EPG shift and transition operators. This ssEPG modeling is validated against Bloch simulations and phantom slice profile measurements. ssEPG relaxometry accuracy and variability are compared with other EPG methods in phantoms and human leg in vivo. The role of ∆B0 interactions with slice profile and gradient crushers is investigated. RESULTS: Simulations and slice profile measurements show that ssEPG can model highly dynamic slice profile effects of gradient-crushed sequences. The MRF ssEPG T2 estimates over 0 < T2 < 100 ms improve accuracy by > 10 ms at some values relative to other modeling approaches. Small deviations in B0 can produce substantial bias in T2 estimations from a range of MRF sequence types, and these effects can be modeled and understood by ssEPG. CONCLUSION: Transient-state, gradient-crushed sequences such as those used in MRF are sensitive to slice profile effects, and these effects depend on RF pulse choice, gradient crusher strength, and ∆B0 . It was found ssEPG was the most accurate EPG-based means to model these effects.

Microscopic susceptibility anisotropy imaging.

PURPOSE: The gradient-echo MR signal in brain white matter depends on the orientation of the fibers with respect to the external magnetic field. To map microstructure-specific magnetic susceptibility in orientationally heterogeneous material, it is thus imperative to regress out unwanted orientation effects. METHODS: This work introduces a novel framework, referred to as microscopic susceptibility anisotropy imaging, that disentangles the 2 principal effects conflated in gradient-echo measurements, (a) the susceptibility properties of tissue microenvironments, especially the myelin microstructure, and (b) the axon orientation distribution relative to the magnetic field. Specifically, we utilize information about the orientational tissue structure inferred from diffusion MRI data to factor out the B0 -direction dependence of the frequency difference signal. RESULTS: A human pilot study at 3 T demonstrates proxy maps of microscopic susceptibility anisotropy unconfounded by fiber crossings and orientation dispersion as well as magnetic field direction. The developed technique requires only a dual-echo gradient-echo scan acquired at 1 or 2 head orientations with respect to the magnetic field and a 2-shell diffusion protocol achievable on standard scanners within practical scan times. CONCLUSIONS: The quantitative recovery of microscopic susceptibility features in the presence of orientational heterogeneity potentially improves the assessment of microstructural tissue integrity.

Relayed nuclear Overhauser enhancement sensitivity to membrane Cho phospholipids.

PURPOSE: Phospholipids are key constituents of cell membranes and serve vital functions in the regulation of cellular processes; thus, a method for in vivo detection and characterization could be valuable for detecting changes in cell membranes that are consequences of either normal or pathological processes. Here, we describe a new method to map the distribution of partially restricted phospholipids in tissues. METHODS: The phospholipids were measured by signal changes caused by relayed nuclear Overhauser enhancement-mediated CEST between the phospholipid Cho headgroup methyl protons and water at around -1.6 ppm from the water resonance. The biophysical basis of this effect was examined by controlled manipulation of head group, chain length, temperature, degree of saturation, and presence of cholesterol. Additional experiments were performed on animal tumor models to evaluate potential applications of this novel signal while correcting for confounding contributions. RESULTS: Negative relayed nuclear Overhauser dips in Z-spectra were measured from reconstituted Cho phospholipids with cholesterol but not for other Cho-containing metabolites or proteins. Significant contrast was found between tumor and contralateral normal tissue signals in animals when comparing both the measured saturation transfer signal and a more specific imaging metric. CONCLUSION: We demonstrated specific relayed nuclear Overhauser effects in partially restricted phospholipid phantoms and similar effects in solid brain tumors after correcting for confounding signal contributions, suggesting possible translational applications of this novel molecular imaging method, which we name restricted phospholipid transfer.

Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry.

PURPOSE: Multi-exponential relaxometry is a powerful tool for characterizing tissue, but generally requires high image signal-to-noise ratio (SNR). This work evaluates the use of principal-component-analysis (PCA) denoising to mitigate these SNR demands and improve the precision of relaxometry measures. METHODS: PCA denoising was evaluated using both simulated and experimental MRI data. Bi-exponential transverse relaxation signals were simulated for a wide range of acquisition and sample parameters, and experimental data were acquired from three excised and fixed mouse brains. In both cases, standard relaxometry analysis was performed on both original and denoised image data, and resulting estimated signal parameters were compared. RESULTS: Denoising reduced the root-mean-square-error of parameters estimated from multi-exponential relaxometry by factors of ≈3×, for typical acquisition and sample parameters. Denoised images and subsequent parameter maps showed little or no signs of spatial artifact or loss of resolution. CONCLUSION: Experimental studies and simulations demonstrate that PCA denoising of MRI relaxometry data is an effective method of improving parameter precision without sacrificing image resolution. This simple yet important processing step thus paves the way for broader applicability of multi-exponential MRI relaxometry.

Recommendations towards standards for quantitative MRI (qMRI) and outstanding needs.

LEVEL OF EVIDENCE: 5 Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2019.

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, T1, T2, and T2*, 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.

Evaluating g-ratio weighted changes in the corpus callosum as a function of age and sex.

Recent years have seen a growing interest in relating MRI measurements to the structural-biophysical properties of white matter fibers. The fiber g-ratio, defined as the ratio between the inner and outer radii of the axon myelin sheath, is an important structural property of white matter, affecting signal conduction. Recently proposed modeling methods that use a combination of quantitative-MRI signals, enable a measurement of the fiber g-ratio in vivo. Here we use an MRI-based g-ratio estimation to observe the variance of the g-ratio within the corpus callosum, and evaluate sex and age related differences. To estimate the g-ratio we used a model (Stikov et al., 2011; Duval et al., 2017) based on two different WM microstructure parameters: the relative amounts of myelin (myelin volume fraction, MVF) and fibers (fiber volume fraction, FVF) in a voxel. We derived the FVF from the fractional anisotropy (FA), and estimated the MVF by using the lipid and macromolecular tissue volume (MTV), calculated from the proton density (Mezer et al., 2013). In comparison to other methods of estimating the MVF, MTV represents a stable parameter with a straightforward route of acquisition. To establish our model, we first compared histological MVF measurements (West et al., 2016) with the MRI derived MTV. We then implemented our model on a large database of 92 subjects (44 males), aged 7 to 81, in order to evaluate age and sex related changes within the corpus callosum. Our results show that the MTV provides a good estimation of MVF for calculating g-ratio, and produced values from the corpus callosum that correspond to those found in animals ex vivo and are close to the theoretical optimum, as well as to published in vivo data. Our results demonstrate that the MTV derived g-ratio provides a simple and reliable in vivo g-ratio-weighted (GR*) measurement in humans. In agreement with theoretical predictions, and unlike other tissue parameters measured with MRI, the g-ratio estimations were found to be relatively stable with age, and we found no support for a significant sexual dimorphism with age.

Myelin volume fraction imaging with MRI.

MRI is a valuable tool to assess myelin during development and demyelinating disease processes. While multiexponential T2 and quantitative magnetization transfer measures correlate with myelin content, neither provides the total myelin volume fraction. In many cases correlative measures are adequate; but to assess microstructure of myelin, (e.g. calculate the g-ratio using MRI), an accurate measure of myelin volume fraction is imperative. Using a volumetric model of white matter, we relate MRI measures of myelin to absolute measures of myelin volume fraction and compare them to quantitative histology. We assess our approach in control mice along with two models of hypomyelination and one model of hypermyelination and find strong agreement between MRI and histology amongst models. This work investigates the sensitivities of MRI myelin measures to changes in axon geometry and displays promise for estimating g-ratio from MRI.

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.

Propagation of error from parameter constraints in quantitative MRI: Example application of multiple spin echo T2 mapping.

PURPOSE: Quantitative MRI may require correcting for nuisance parameters which can or must be constrained to independently measured or assumed values. The noise and/or bias in these constraints propagate to fitted parameters. For example, the case of refocusing pulse flip angle constraint in multiple spin echo T2 mapping is explored. METHODS: An analytical expression for the mean-squared error of a parameter of interest was derived as a function of the accuracy and precision of an independent estimate of a nuisance parameter. The expression was validated by simulations and then used to evaluate the effects of flip angle (θ) constraint on the accuracy and precision of T⁁2 for a variety of multi-echo T2 mapping protocols. RESULTS: Constraining θ improved T⁁2 precision when the θ-map signal-to-noise ratio was greater than approximately one-half that of the first spin echo image. For many practical scenarios, constrained fitting was calculated to reduce not just the variance but the full mean-squared error of T⁁2, for bias in θ⁁≲6%. CONCLUSION: The analytical expression derived in this work can be applied to inform experimental design in quantitative MRI. The example application to T2 mapping provided specific cases, depending on θ⁁ accuracy and precision, in which θ⁁ measurement and constraint would be beneficial to T⁁2 variance or mean-squared error. Magn Reson Med 79:673-682, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Chemical exchange rotation transfer (CERT) on human brain at 3 Tesla.

PURPOSE: To test the ability of a novel pulse sequence applied in vivo at 3 Tesla to separate the contributions to the water signal from amide proton transfer (APT) and relayed nuclear Overhauser enhancement (rNOE) from background direct water saturation and semisolid magnetization transfer (MT). The lack of such signal source isolation has confounded conventional chemical exchange saturation transfer (CEST) imaging. METHODS: We quantified APT and rNOE signals using a chemical exchange rotation transfer (CERT) metric, MTRdouble . A range of duty cycles and average irradiation powers were applied, and results were compared with conventional CEST analyses using asymmetry (MTRasym ) and extrapolated magnetization transfer (EMR). RESULTS: Our results indicate that MTRdouble is more specific than MTRasym and, because it requires as few as 3 data points, is more rapid than methods requiring a complete Z-spectrum, such as EMR. In white matter, APT (1.5 ± 0.5%) and rNOE (2.1 ± 0.7%) were quantified by using MTRdouble with a 30% duty cycle and a 0.5-µT average power. In addition, our results suggest that MTRdouble is insensitive to B0 inhomogeneity, further magnifying its speed advantage over CEST metrics that require a separate B0 measurement. However, MTRdouble still has nontrivial sensitivity to B1 inhomogeneities. CONCLUSION: We demonstrated that MTRdouble is an alternative metric to evaluate APT and rNOE, which is fast, robust to B0 inhomogeneity, and easy to process.