Signatures of microstructure in gradient-echo and spin-echo signals.

PURPOSE: To determine whether the spatial scale and magnetic susceptibility of microstructure can be evaluated robustly from the decay of gradient-echo and spin-echo signals. THEORY AND METHODS: Gradient-echo and spin-echo images were acquired from suspensions of spherical polystyrene microbeads of 10, 20, and 40 µm nominal diameter. The sizes of the beads and their magnetic susceptibility relative to the medium were estimated from the signal decay curves, using a lookup table generated from Monte Carlo simulations and an analytic model based on the Gaussian phase approximation. RESULTS: Fitting Monte Carlo predictions to spin-echo data yielded acceptable estimates of microstructural parameters for the 20 and 40 µm microbeads. Using gradient-echo data, the Monte Carlo lookup table provided satisfactory parameter estimates for the 20 µm beads but unstable results for the diameter of the largest beads. Neither spin-echo nor gradient-echo data allowed accurate parameter estimation for the smallest beads. The analytic model performed poorly over all bead sizes. CONCLUSIONS: Microstructural sources of magnetic susceptibility produce distinctive non-exponential signatures in the decay of gradient-echo and spin-echo signals. However, inverting the problem to extract microstructural parameters from the signals is nontrivial and, in certain regimes, ill-conditioned. For microstructure with small characteristic length scales, parameter estimation is hampered by the difficulty of acquiring accurate data at very short echo times. For microstructure with large characteristic lengths, the gradient-echo signal approaches the static-dephasing regime, where it becomes insensitive to size. Applicability of the analytic model was further limited by failure of the Gaussian phase approximation for all but the smallest beads.

The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy.

The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.

Measuring water exchange on a preclinical MRI system using filter exchange and diffusion time dependent kurtosis imaging.

PURPOSE: Filter exchange imaging (FEXI) and diffusion time (t)-dependent diffusion kurtosis imaging (DKI(t)) are both sensitive to water exchange between tissue compartments. The restrictive effects of tissue microstructure, however, introduce bias to the exchange rate obtained by these two methods, as their interpretation conventionally rely on the Kärger model of barrier limited exchange between Gaussian compartments. Here, we investigated whether FEXI and DKI(t) can provide comparable exchange rates in ex vivo mouse brains. THEORY AND METHODS: FEXI and DKI(t) data were acquired from ex vivo mouse brains on a preclinical MRI system. Phase cycling and negative slice prewinder gradients were used to minimize the interferences from imaging gradients. RESULTS: In the corpus callosum, apparent exchange rate (AXR) from FEXI correlated with the exchange rate (the inverse of exchange time, 1/τex ) from DKI(t) along the radial direction. In comparison, discrepancies between FEXI and DKI(t) were found in the cortex due to low filter efficiency and confounding effects from tissue microstructure. CONCLUSION: The results suggest that FEXI and DKI(t) are sensitive to the same exchange processes in white matter when separated from restrictive effects of microstructure. The complex microstructure in gray matter, with potential exchange among multiple compartments and confounding effects of microstructure, still pose a challenge for FEXI and DKI(t).

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange.

Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging anisotropic compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along neuronal and glial processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Kärger model of two exchanging compartments. Using datasets acquired at multiple diffusion weightings (b) and diffusion times (t) in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the gray matter, compared to the other two possible contributions. Our results for the diffusion time window 20-45 ms show minimal diffusivity time-dependence and more pronounced kurtosis decay with time, which is well fit by the exchange model. Moreover, we observe lower signal for longer diffusion times at high b. In light of these observations, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for gray matter microstructure mapping. We finally highlight multi-b multi-t acquisition protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 15 - 60 ms in the rat cortex and hippocampus in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions. This suggests water exchange as an essential component for interpreting diffusion MRI measurements in gray matter.

Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems.

Estimating intra- and extra-axonal microstructure parameters, such as volume fractions and diffusivities, has been one of the major efforts in brain microstructure imaging with MRI. The Standard Model (SM) of diffusion in white matter has unified various modeling approaches based on impermeable narrow cylinders embedded in locally anisotropic extra-axonal space. However, estimating the SM parameters from a set of conventional diffusion MRI (dMRI) measurements is ill-conditioned. Multidimensional dMRI helps resolve the estimation degeneracies, but there remains a need for clinically feasible acquisitions that yield robust parameter maps. Here we find optimal multidimensional protocols by minimizing the mean-squared error of machine learning-based SM parameter estimates for two 3T scanners with corresponding gradient strengths of 40and80mT/m. We assess intra-scanner and inter-scanner repeatability for 15-minute optimal protocols by scanning 20 healthy volunteers twice on both scanners. The coefficients of variation all SM parameters except free water fraction are ≲10% voxelwise and 1-4% for their region-averaged values. As the achieved SM reproducibility outcomes are similar to those of conventional diffusion tensor imaging, our results enable robust in vivo mapping of white matter microstructure in neuroscience research and in the clinic.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome.

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Improved Task-based Functional MRI Language Mapping in Patients with Brain Tumors through Marchenko-Pastur Principal Component Analysis Denoising.

Background Functional MRI improves preoperative planning in patients with brain tumors, but task-correlated signal intensity changes are only 2%-3% above baseline. This makes accurate functional mapping challenging. Marchenko-Pastur principal component analysis (MP-PCA) provides a novel strategy to separate functional MRI signal from noise without requiring user input or prior data representation. Purpose To determine whether MP-PCA denoising improves activation magnitude for task-based functional MRI language mapping in patients with brain tumors. Materials and Methods In this Health Insurance Portability and Accountability Act-compliant study, MP-PCA performance was first evaluated by using simulated functional MRI data with a known ground truth. Right-handed, left-language-dominant patients with brain tumors who successfully performed verb generation, sentence completion, and finger tapping functional MRI tasks were retrospectively identified between January 2017 and August 2018. On the group level, for each task, histograms of z scores for original and MP-PCA denoised data were extracted from relevant regions and contralateral homologs were seeded by a neuroradiologist blinded to functional MRI findings. Z scores were compared with paired two-sided t tests, and distributions were compared with effect size measurements and the Kolmogorov-Smirnov test. The number of voxels with a z score greater than 3 was used to measure task sensitivity relative to task duration. Results Twenty-three patients (mean age ± standard deviation, 43 years ± 18; 13 women) were evaluated. MP-PCA denoising led to a higher median z score of task-based functional MRI voxel activation in left hemisphere cortical regions for verb generation (from 3.8 ± 1.0 to 4.5 ± 1.4; P < .001), sentence completion (from 3.7 ± 1.0 to 4.3 ± 1.4; P < .001), and finger tapping (from 6.9 ± 2.4 to 7.9 ± 2.9; P < .001). Median z scores did not improve in contralateral homolog regions for verb generation (from -2.7 ± 0.54 to -2.5 ± 0.40; P = .90), sentence completion (from -2.3 ± 0.21 to -2.4 ± 0.37; P = .39), or finger tapping (from -2.3 ± 1.20 to -2.7 ± 1.40; P = .07). Individual functional MRI task durations could be truncated by at least 40% after MP-PCA without degradation of clinically relevant correlations between functional cortex and functional MRI tasks. Conclusion Denoising with Marchenko-Pastur principal component analysis led to higher task correlations in relevant cortical regions during functional MRI language mapping in patients with brain tumors. © RSNA, 2020 Online supplemental material is available for this article.

Removal of partial Fourier-induced Gibbs (RPG) ringing artifacts in MRI.

PURPOSE: To investigate and remove Gibbs-ringing artifacts caused by partial Fourier (PF) acquisition and zero filling interpolation in MRI data. THEORY AND METHODS: Gibbs ringing of fully sampled data, leading to oscillations around tissue boundaries, is caused by the symmetric truncation of k-space. Such ringing can be removed by conventional methods, with the local subvoxel shifts method being the state-of-the-art. However, the asymmetric truncation of k-space in routinely used PF acquisitions leads to additional ringings of wider intervals in the PF sampling dimension that cannot be corrected solely based on magnitude images reconstructed via zero filling. Here, we develop a pipeline for the Removal of PF-induced Gibbs ringing (RPG) to remove ringing patterns of different periods by applying the conventional method twice. The proposed pipeline is validated on numerical phantoms, demonstrated on in vivo diffusion MRI measurements, and compared with the conventional method and neural network-based approach. RESULTS: For PF = 7/8 and 6/8, Gibbs-ringings and subsequent bias in diffusion metrics induced by PF acquisition and zero filling are robustly removed by using the proposed RPG pipeline. For PF = 5/8, however, ringing removal via RPG leads to excessive image blurring due to the interplay of image phase and convolution kernel. CONCLUSIONS: RPG corrects Gibbs-ringing artifacts in magnitude images of PF acquired data and reduces the bias in quantitative MR metrics. Considering the benefit of PF acquisition and the feasibility of ringing removal, we suggest applying PF = 6/8 when PF acquisition is necessary.

Training a neural network for Gibbs and noise removal in diffusion MRI.

PURPOSE: To develop and evaluate a neural network-based method for Gibbs artifact and noise removal. METHODS: A convolutional neural network (CNN) was designed for artifact removal in diffusion-weighted imaging data. Two implementations were considered: one for magnitude images and one for complex images. Both models were based on the same encoder-decoder structure and were trained by simulating MRI acquisitions on synthetic non-MRI images. RESULTS: Both machine learning methods were able to mitigate artifacts in diffusion-weighted images and diffusion parameter maps. The CNN for complex images was also able to reduce artifacts in partial Fourier acquisitions. CONCLUSIONS: The proposed CNNs extend the ability of artifact correction in diffusion MRI. The machine learning method described here can be applied on each imaging slice independently, allowing it to be used flexibly in clinical applications.

Assessment of myofiber microstructure changes due to atrophy and recovery with time-dependent diffusion MRI.

Current clinical MRI evaluation of musculature largely focuses on nonquantitative assessments (including T1-, T2- and PD-weighted images), which may vary greatly between imaging systems and readers. This work aims to determine the efficacy of a quantitative approach to study the microstructure of muscles at the cellular level with the random permeable barrier model (RPBM) applied to time-dependent diffusion tensor imaging (DTI) for varying diffusion time. Patients (N = 15, eight males and seven females) with atrophied calf muscles due to immobilization of one leg in a nonweight-bearing cast, were enrolled after providing informed consent. Their calf muscles were imaged with stimulated echo diffusion for DTI, T1-mapping and RPBM modeling. Specifically, After cast removal, both calf muscles (atrophied and contralateral control leg) were imaged with MRI for all patients, with follow-up scans to monitor recovery of the atrophied leg for six patients after 4 and 8 weeks. We compare RPBM-derived microstructural metrics: myofiber diameter, a, and sarcolemma permeability, κ, along with macroscopic anatomical parameters (muscle cross-sectional area, fiber orientation, <θ>, and T1 relaxation). ROC analysis was used to compare parameters between control and atrophied muscle, while the Friedman test was used to evaluate the atrophied muscle longitudinally. We found that the RPBM framework enables measurement of microstructural parameters from diffusion time-dependent DTI, of which the myofiber diameter is a stronger predictor of intramuscular morphological changes than either macroscopic (anatomical) measurements or empirical diffusion parameters. This work demonstrates the potential of RPBM to assess pathological changes in musculature that seem undetectable with standard diffusion and anatomical MRI.

Measurement of cellular-interstitial water exchange time in tumors based on diffusion-time-dependent diffusional kurtosis imaging.

PURPOSE: To assess the feasibility of using diffusion-time-dependent diffusional kurtosis imaging (tDKI) to measure cellular-interstitial water exchange time (τex ) in tumors, both in animals and in humans. METHODS: Preclinical tDKI studies at 7 T were performed with the GL261 glioma model and the 4T1 mammary tumor model injected into the mouse brain. Clinical studies were performed at 3 T with women who had biopsy-proven invasive ductal carcinoma. tDKI measurement was conducted using a diffusion-weighted STEAM pulse sequence with multiple diffusion times (20-800 ms) at a fixed echo time, while keeping the b-values the same (0-3000 s/mm2 ) by adjusting the diffusion gradient strength. The tDKI data at each diffusion time t were used for a weighted linear least-squares fit method to estimate the diffusion-time-dependent diffusivity, D(t), and diffusional kurtosis, K(t). RESULTS: Both preclinical and clinical studies showed that, when diffusion time t ≥ 200 ms, D(t) did not have a noticeable change while K(t) decreased monotonically with increasing diffusion time in tumors and t ≥ 100 ms for the cortical ribbon of the mouse brain. The estimated τex averaged median and interquartile range (IQR) of GL261 and 4T1 tumors were 93 (IQR = 89) ms and 68 (78) ms, respectively. For the cortical ribbon, the estimated τex averaged median and IQR were 41 (34) ms for C57BL/6 and 30 (17) ms for BALB/c. For invasive ductal carcinoma, the estimated τex median and IQR of the two breast cancers were 70 (94) and 106 (92) ms. CONCLUSION: The results of this proof-of-concept study substantiate the feasibility of using tDKI to measure cellular-interstitial water exchange time without using an exogenous contrast agent.

In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter.

The dependence of the diffusion MRI signal on the diffusion time t is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity D(t) and mean kurtosis K(t) in cortical gray matter and in 25 â€‹gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for t=21.2-100 â€‹ms, and is characterized by a power-law tail ∼t-ϑ with dynamical exponent ϑ. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using the effective medium theory formalism, we derive an exact relation between the power-law tails in D(t) and K(t). The estimated dynamical exponent ϑ≃1/2 in both D(t) and K(t) is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the ϑ=1/2 exponent. In contrast, the Kärger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from K(t) falls below our measured t-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the t-1/2 tail in D(t) and K(t) altogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.

The impact of realistic axonal shape on axon diameter estimation using diffusion MRI.

To study axonal microstructure with diffusion MRI, axons are typically modeled as straight impermeable cylinders, whereby the transverse diffusion MRI signal can be made sensitive to the cylinder's inner diameter. However, the shape of a real axon varies along the axon direction, which couples the longitudinal and transverse diffusion of the overall axon direction. Here we develop a theory of the intra-axonal diffusion MRI signal based on coarse-graining of the axonal shape by 3-dimensional diffusion. We demonstrate how the estimate of the inner diameter is confounded by the diameter variations (beading), and by the local variations in direction (undulations) along the axon. We analytically relate diffusion MRI metrics, such as time-dependent radial diffusivity D⊥(t)and kurtosis K⊥(t),to the axonal shape, and validate our theory using Monte Carlo simulations in synthetic undulating axons with randomly positioned beads, and in realistic axons reconstructed from electron microscopy images of mouse brain white matter. We show that (i) In the narrow pulse limit, the inner diameter from D⊥(t)is overestimated by about twofold due to a combination of axon caliber variations and undulations (each contributing a comparable effect size); (ii) The narrow-pulse kurtosis K⊥|t→∞deviates from that in an ideal cylinder due to caliber variations; we also numerically calculate the fourth-order cumulant for an ideal cylinder in the wide pulse limit, which is relevant for inner diameter overestimation; (iii) In the wide pulse limit, the axon diameter overestimation is mainly due to undulations at low diffusion weightings b; and (iv) The effect of undulations can be considerably reduced by directional averaging of high-b signals, with the apparent inner diameter given by a combination of the axon caliber (dominated by the thickest axons), caliber variations, and the residual contribution of undulations.

Retrieving neuronal orientations using 3D scanning SAXS and comparison with diffusion MRI.

While diffusion MRI (dMRI) is currently the method of choice to non-invasively probe tissue microstructure and study structural connectivity in the brain, its spatial resolution is limited and its results need structural validation. Current ex vivo methods employed to provide 3D fiber orientations have limitations, including tissue-distorting sample preparation, small field of view or inability to quantify 3D fiber orientation distributions. 3D fiber orientation in tissue sections can be obtained from 3D scanning small-angle X-ray scattering (3D sSAXS) by analyzing the anisotropy of scattering signals. Here we adapt the 3D sSAXS method for use in brain tissue, exploiting the high sensitivity of the SAXS signal to the ordered molecular structure of myelin. We extend the characterization of anisotropy from vectors to tensors, employ the Funk-Radon-Transform for converting scattering information to real space fiber orientations, and demonstrate the feasibility of the method in thin sections of mouse brain with minimal sample preparation. We obtain a second rank tensor representing the fiber orientation distribution function (fODF) for every voxel, thereby generating fODF maps. Finally, we illustrate the potential of 3D sSAXS by comparing the result with diffusion MRI fiber orientations in the same mouse brain. We show a remarkably good correspondence, considering the orthogonality of the two methods, i.e. the different physical processes underlying the two signals. 3D sSAXS can serve as validation method for microstructural MRI, and can provide novel microstructural insights for the nervous system, given the method's orthogonality to dMRI, high sensitivity to myelin sheath's orientation and abundance, and the possibility to extract myelin-specific signal and to perform micrometer-resolution scanning.

Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising.

Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.

Effect of intravoxel incoherent motion on diffusion parameters in normal brain.

At very low diffusion weighting the diffusion MRI signal is affected by intravoxel incoherent motion (IVIM) caused by dephasing of magnetization due to incoherent blood flow in capillaries or other sources of microcirculation. While IVIM measurements at low diffusion weightings have been frequently used to investigate perfusion in the body as well as in malignant tissue, the effect and origin of IVIM in normal brain tissue is not completely established. We investigated the IVIM effect on the brain diffusion MRI signal in a cohort of 137 radiologically-normal patients (62 male; mean age = 50.2 ±â€¯17.8, range = 18 to 94). We compared the diffusion tensor parameters estimated from a mono-exponential fit at b = 0 and 1000 s/mm2 versus at b = 250 and 1000 s/mm2. The asymptotic fitting method allowed for quantitative assessment of the IVIM signal fraction f* in specific brain tissue and regions. Our results show a mean (median) percent difference in the mean diffusivity of about 4.5 (4.9)% in white matter (WM), about 7.8 (8.7)% in cortical gray matter (GM), and 4.3 (4.2)% in thalamus. Corresponding perfusion fraction f* was estimated to be 0.033 (0.032) in WM, 0.066 (0.065) in cortical GM, and 0.033 (0.030) in the thalamus. The effect of f* with respect to age was found to be significant in cortical GM (Pearson correlation ρ â€‹= â€‹0.35, p â€‹= â€‹3*10-5) and the thalamus (Pearson correlation ρ = 0.20, p = 0.022) with an average increase in f* of 5.17*10-4/year and 3.61*10-4/year, respectively. Significant correlations between f* and age were not observed for WM, and corollary analysis revealed no effect of gender on f*. Possible origins of the IVIM effect in normal brain tissue are discussed.

On the scaling behavior of water diffusion in human brain white matter.

Development of therapies for neurological disorders depends on our ability to non-invasively diagnose and monitor the progression of underlying pathologies at the cellular level. Physics and physiology limit the resolution of human MRI to be orders of magnitude coarser than cell dimensions. Here we identify and quantify the MRI signal coming from within micrometer-thin axons in human white matter tracts in vivo, by utilizing the sensitivity of diffusion MRI to Brownian motion of water molecules restricted by cell walls. We study a specific power-law scaling of the diffusion MRI signal with the diffusion weighting, predicted for water confined to narrow axons, and quantify axonal water fraction and orientation dispersion.

Genomic kinship construction to enhance genetic analyses in the human connectome project data.

Imaging genetic analyses quantify genetic control over quantitative measurements of brain structure and function using coefficients of relationship (CR) that code the degree of shared genetics between subjects. CR can be inferred through self-reported relatedness or calculated empirically using genome-wide SNP scans. We hypothesized that empirical CR provides a more accurate assessment of shared genetics than self-reported relatedness. We tested this in 1,046 participants of the Human Connectome Project (HCP) (480 M/566 F) recruited from the Missouri twin registry. We calculated the heritability for 17 quantitative traits drawn from four categories (brain diffusion and structure, cognition, and body physiology) documented by the HCP. We compared the heritability and genetic correlation estimates calculated using self-reported and empirical CR methods Kinship-based INference for GWAS (KING) and weighted allelic correlation (WAC). The polygenetic nature of traits was assessed by calculating the empirical CR from chromosomal SNP sets. The heritability estimates based on whole-genome empirical CR were higher but remained significantly correlated (r ∼0.9) with those obtained using self-reported values. Population stratification in the HCP sample has likely influenced the empirical CR calculations and biased heritability estimates. Heritability values calculated using empirical CR for chromosomal SNP sets were significantly correlated with the chromosomal length (r 0.7) suggesting a polygenic nature for these traits. The chromosomal heritability patterns were correlated among traits from the same knowledge domains; among traits with significant genetic correlations; and among traits sharing biological processes, without being genetically related. The pedigree structures generated in our analyses are available online as a web-based calculator (www.solar-eclipse-genetics.org/HCP).

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation.

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

TE dependent Diffusion Imaging (TEdDI) distinguishes between compartmental T2 relaxation times.

Biophysical modeling of macroscopic diffusion-weighted MRI signal in terms of microscopic cellular parameters holds the promise of quantifying the integrity of white matter. Unfortunately, even fairly simple multi-compartment models of proton diffusion in the white matter do not provide a unique, biophysically plausible solution. Here we report a nontrivial diffusion MRI signal dependence on echo time (TE) in human white matter in vivo. We demonstrate that such TE dependence originates from compartment-specific T2 values and that it is a promising "orthogonal measure" able to break the degeneracy in parameter estimation, and to yield important relaxation metrics robustly. We thereby enable the precise estimation of the intra- and extra-axonal water T2 relaxation times, which is precluded by a limited signal-to-noise ratio when using multi-echo relaxometry alone.