Normative volumes and relaxation times at 3T during brain development

While research has unveiled and quantified brain markers of abnormal neurodevelopment, clinicians still work with qualitative metrics for MRI brain investigation. The purpose of the current article is to bridge the knowledge gap between case-control cohort studies and individual patient care. Here, we provide a unique dataset of seventy-three 3-to-17 years-old healthy subjects acquired with a 6-minute MRI protocol encompassing T1 and T2 relaxation quantitative sequence that can be readily implemented in the clinical setting; MP2RAGE for T1 mapping and the prototype sequence GRAPPATINI for T2 mapping. White matter and grey matter volumes were automatically quantified. We further provide normative developmental curves based on these two imaging sequences; T1, T2 and volume normative ranges with respect to age were computed, for each ROI of a pediatric brain atlas. This open-source dataset provides normative values allowing to position individual patients acquired with the same protocol on the brain maturation curve and as such provides potentially useful quantitative biomarkers facilitating precise and personalized care.

Tract-wise microstructural analysis informs on current and future disability in early multiple sclerosis

OBJECTIVES

Microstructural characterization of patients with multiple sclerosis (MS) has been shown to correlate better with disability compared to conventional radiological biomarkers. Quantitative MRI provides effective means to characterize microstructural brain tissue changes both in lesions and normal-appearing brain tissue. However, the impact of the location of microstructural alterations in terms of neuronal pathways has not been thoroughly explored so far. Here, we study the extent and the location of tissue changes probed using quantitative MRI along white matter (WM) tracts extracted from a connectivity atlas.

METHODS

We quantified voxel-wise T1 tissue alterations compared to normative values in a cohort of 99 MS patients. For each WM tract, we extracted metrics reflecting tissue alterations both in lesions and normal-appearing WM and correlated these with cross-sectional disability and disability evolution after 2 years.

RESULTS

In early MS patients, T1 alterations in normal-appearing WM correlated better with disability evolution compared to cross-sectional disability. Further, the presence of lesions in supratentorial tracts was more strongly associated with cross-sectional disability, while microstructural alterations in infratentorial pathways yielded higher correlations with disability evolution. In progressive patients, all major WM pathways contributed similarly to explaining disability, and correlations with disability evolution were generally poor.

CONCLUSIONS

We showed that microstructural changes evaluated in specific WM pathways contribute to explaining future disability in early MS, hence highlighting the potential of tract-wise analyses in monitoring disease progression. Further, the proposed technique allows to estimate WM tract-specific microstructural characteristics in clinically compatible acquisition times, without the need for advanced diffusion imaging.

Submillimeter T1 atlas for subject-specific abnormality detection at 7T

PURPOSE

Studies at 3T have shown that T₁ relaxometry enables characterization of brain tissues at the single-subject level by comparing individual physical properties to a normative atlas. In this work, an atlas of normative T₁ values at 7T is introduced with 0.6 mm isotropic resolution and its clinical potential is explored in comparison to 3T.

METHODS

T₁ maps were acquired in two separate healthy cohorts scanned at 3T and 7T. Using transfer learning, a template-based brain segmentation algorithm was adapted to ultra-high field imaging data. After segmenting brain tissues, volumes were normalized into a common space, and an atlas of normative T₁ values was established by modeling the T₁ inter-subject variability. A method for single-subject comparisons restricted to white matter and subcortical structures was developed by computing Z-scores. The comparison was applied to eight patients scanned at both field strengths for proof of concept.

RESULTS

The proposed method for morphometry delivered segmentation masks without statistically significant differences from those derived with the original pipeline at 3T and achieved accurate segmentation at 7T. The established normative atlas allowed characterizing tissue alterations in single-subject comparisons at 7T, and showed greater anatomical details compared with 3T results.

CONCLUSION

A high-resolution quantitative atlas with an adapted pipeline was introduced and validated. Several case studies on different clinical conditions showed the feasibility, potential and limitations of high-resolution single-subject comparisons based on quantitative MRI atlases. This method in conjunction with 7T higher resolution broadens the range of potential applications of quantitative MRI in clinical practice.

Variability and reproducibility of multi-echo T2 relaxometry: Insights from multi-site, multi-session and multi-subject MRI acquisitions

Quantitative magnetic resonance imaging (qMRI) can increase the specificity and sensitivity of conventional weighted MRI to underlying pathology by comparing meaningful physical or chemical parameters, measured in physical units, with normative values acquired in a healthy population. This study focuses on multi-echo T₂ relaxometry, a qMRI technique that probes the complex tissue microstructure by differentiating compartment-specific T₂ relaxation times. However, estimation methods are still limited by their sensitivity to the underlying noise. Moreover, estimating the model's parameters is challenging because the resulting inverse problem is ill-posed, requiring advanced numerical regularization techniques. As a result, the estimates from distinct regularization strategies are different. In this work, we aimed to investigate the variability and reproducibility of different techniques for estimating the transverse relaxation time of the intra- and extra-cellular space (T2IE) in gray (GM) and white matter (WM) tissue in a clinical setting, using a multi-site, multi-session, and multi-run T₂ relaxometry dataset. To this end, we evaluated three different techniques for estimating the T₂ spectra (two regularized non-negative least squares methods and a machine learning approach). Two independent analyses were performed to study the effect of using raw and denoised data. For both the GM and WM regions, and the raw and denoised data, our results suggest that the principal source of variance is the inter-subject variability, showing a higher coefficient of variation (CoV) than those estimated for the inter-site, inter-session, and inter-run, respectively. For all reconstruction methods studied, the CoV ranged between 0.32 and 1.64%. Interestingly, the inter-session variability was close to the inter-scanner variability with no statistical differences, suggesting that T2IE is a robust parameter that could be employed in multi-site neuroimaging studies. Furthermore, the three tested methods showed consistent results and similar intra-class correlation (ICC), with values superior to 0.7 for most regions. Results from raw data were slightly more reproducible than those from denoised data. The regularized non-negative least squares method based on the L-curve technique produced the best results, with ICC values ranging from 0.72 to 0.92.

Data-driven myelin water imaging based on T1 and T2 relaxometry

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T₁ and T₂ relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T₁ /T₂ data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions.

Evaluating reproducibility and subject-specificity of microstructure-informed connectivity

Tractography enables identifying and evaluating the healthy and diseased brain's white matter pathways from diffusion-weighted magnetic resonance imaging data. As previous evaluation studies have reported significant false-positive estimation biases, recent microstructure-informed tractography algorithms have been introduced to improve the trade-off between specificity and sensitivity. However, a major limitation for characterizing the performance of these techniques is the lack of ground truth brain data. In this study, we compared the performance of two relevant microstructure-informed tractography methods, SIFT2 and COMMIT, by assessing the subject specificity and reproducibility of their derived white matter pathways. Specifically, twenty healthy young subjects were scanned at eight different time points at two different sites. Subject specificity and reproducibility were evaluated using the whole-brain connectomes and a subset of 29 white matter bundles. Our results indicate that although the raw tractograms are more vulnerable to the presence of false-positive connections, they are highly reproducible, suggesting that the estimation bias is subject-specific. This high reproducibility was preserved when microstructure-informed tractography algorithms were used to filter the raw tractograms. Moreover, the resulting track-density images depicted a more uniform coverage of streamlines throughout the white matter, suggesting that these techniques could increase the biological meaning of the estimated fascicles. Notably, we observed an increased subject specificity by employing connectivity pre-processing techniques to reduce the underlaying noise and the data dimensionality (using principal component analysis), highlighting the importance of these tools for future studies. Finally, no strong bias from the scanner site or time between measurements was found. The largest intraindividual variance originated from the sole repetition of data measurements (inter-run).

Periventricular gradient of T1 tissue alterations in multiple sclerosis

•

T₁ relaxation times alterations were assessed based on a study-specific atlas.

•

T₁ alterations depend on distance from lateral ventricles (“periventricular gradient”).

•

Gradient parameters correlate better with disability compared to conventional MRI.

Objective

Pathology in multiple sclerosis is not homogenously distributed. Recently, it has been shown that structures adjacent to CSF are more severely affected. A gradient of brain tissue involvement was shown with more severe pathology in periventricular areas and in proximity to brain surfaces such as the subarachnoid spaces and ependyma, and hence termed the “surface–in” gradient. Here, we study whether (i) the surface-in gradient of periventricular tissue alteration measured by T₁ relaxometry is already present in early multiple sclerosis patients, (ii) how it differs between early and progressive multiple sclerosis patients, and (iii) whether the gradient-derived metrics in normal-appearing white matter and lesions correlate better with physical disability than conventional MRI-based metrics.

Methods

Forty-seven patients with early multiple sclerosis, 52 with progressive multiple sclerosis, and 92 healthy controls were included in the study. Isotropic 3D T₁ relaxometry maps were obtained using the Magnetization-Prepared 2 Rapid Acquisition Gradient Echoes sequence at 3 T. After spatially normalizing the T₁ maps into a study-specific common space, T₁ inter-subject variability within the healthy cohort was modelled voxel-wise, yielding a normative T₁ atlas. Individual comparisons of each multiple sclerosis patient against the atlas were performed by computing z-scores. Equidistant bands of voxels were defined around the ventricles in the supratentorial white matter; the z-scores in these bands were analysed and compared between the early and progressive multiple sclerosis cohorts. Correlations between both conventional and z-score-gradient-derived MRI metrics and the Expanded Disability Status Scale were assessed.

Results

Patients with early and progressive multiple sclerosis demonstrated a periventricular gradient of T₁ relaxation time z-scores. In progressive multiple sclerosis, z-score-derived metrics reflecting the gradient of tissue abnormality in normal-appearing white matter were more strongly correlated with disability (maximal rho = 0.374) than the conventional lesion volume and count (maximal rho = 0.189 and 0.21 respectively). In early multiple sclerosis, the gradient of normal-appearing white matter volume with z-scores > 2 at baseline correlated with clinical disability assessed at two years follow-up.

Conclusion

Our results suggest that the surface-in white matter gradient of tissue alteration is detectable with T1 relaxometry and is already present at clinical disease onset. The periventricular gradients correlate with clinical disability. The periventricular gradient in normal-appearing white matter may thus qualify as a promising biomarker for monitoring of disease activity from an early stage in all phenotypes of multiple sclerosis.

Revisiting the T2 spectrum imaging inverse problem: Bayesian regularized non-negative least squares

Multi-echo T₂ magnetic resonance images contain information about the distribution of T₂ relaxation times of compartmentalized water, from which we can estimate relevant brain tissue properties such as the myelin water fraction (MWF). Regularized non-negative least squares (NNLS) is the tool of choice for estimating non-parametric T₂ spectra. However, the estimation is ill-conditioned, sensitive to noise, and highly affected by the employed regularization weight. The purpose of this study is threefold: first, we want to underline that the apparently innocuous use of two alternative parameterizations for solving the inverse problem, which we called the standard and alternative regularization forms, leads to different solutions; second, to assess the performance of both parameterizations; and third, to propose a new Bayesian regularized NNLS method (BayesReg). The performance of BayesReg was compared with that of two conventional approaches (L-curve and Chi-square (X²) fitting) using both regularization forms. We generated a large dataset of synthetic data, acquired in vivo human brain data in healthy participants for conducting a scan-rescan analysis, and correlated the myelin content derived from histology with the MWF estimated from ex vivo data. Results from synthetic data indicate that BayesReg provides accurate MWF estimates, comparable to those from L-curve and X², and with better overall stability across a wider signal-to-noise range. Notably, we obtained superior results by using the alternative regularization form. The correlations reported in this study are higher than those reported in previous studies employing the same ex vivo and histological data. In human brain data, the estimated maps from L-curve and BayesReg were more reproducible. However, the T₂ spectra produced by BayesReg were less affected by over-smoothing than those from L-curve. These findings suggest that BayesReg is a good alternative for estimating T₂ distributions and MWF maps.

Comparison of non-parametric T2 relaxometry methods for myelin water quantification

Multi-component T₂ relaxometry allows probing tissue microstructure by assessing compartment-specific T₂ relaxation times and water fractions, including the myelin water fraction. Non-negative least squares (NNLS) with zero-order Tikhonov regularization is the conventional method for estimating smooth T₂ distributions. Despite the improved estimation provided by this method compared to non-regularized NNLS, the solution is still sensitive to the underlying noise and the regularization weight. This is especially relevant for clinically achievable signal-to-noise ratios. In the literature of inverse problems, various well-established approaches to promote smooth solutions, including first-order and second-order Tikhonov regularization, and different criteria for estimating the regularization weight have been proposed, such as L-curve, Generalized Cross-Validation, and Chi-square residual fitting. However, quantitative comparisons between the available reconstruction methods for computing the T₂ distribution, and between different approaches for selecting the optimal regularization weight, are lacking. In this study, we implemented and evaluated ten reconstruction algorithms, resulting from the individual combinations of three penalty terms with three criteria to estimate the regularization weight, plus non-regularized NNLS. Their performance was evaluated both in simulated data and real brain MRI data acquired from healthy volunteers through a scan-rescan repeatability analysis. Our findings demonstrate the need for regularization. As a result of this work, we provide a list of recommendations for selecting the optimal reconstruction algorithms based on the acquired data. Moreover, the implemented methods were packaged in a freely distributed toolbox to promote reproducible research, and to facilitate further research and the use of this promising quantitative technique in clinical practice.

Model-informed machine learning for multi-component T2 relaxometry

Recovering the T₂ distribution from multi-echo T₂ magnetic resonance (MR) signals is challenging but has high potential as it provides biomarkers characterizing the tissue micro-structure, such as the myelin water fraction (MWF). In this work, we propose to combine machine learning and aspects of parametric (fitting from the MRI signal using biophysical models) and non-parametric (model-free fitting of the T₂ distribution from the signal) approaches to T₂ relaxometry in brain tissue by using a multi-layer perceptron (MLP) for the distribution reconstruction. For training our network, we construct an extensive synthetic dataset derived from biophysical models in order to constrain the outputs with a priori knowledge of in vivo distributions. The proposed approach, called Model-Informed Machine Learning (MIML), takes as input the MR signal and directly outputs the associated T₂ distribution. We evaluate MIML in comparison to a Gaussian Mixture Fitting (parametric) and Regularized Non-Negative Least Squares algorithms (non-parametric) on synthetic data, an ex vivo scan, and high-resolution scans of healthy subjects and a subject with Multiple Sclerosis. In synthetic data, MIML provides more accurate and noise-robust distributions. In real data, MWF maps derived from MIML exhibit the greatest conformity to anatomical scans, have the highest correlation to a histological map of myelin volume, and the best unambiguous lesion visualization and localization, with superior contrast between lesions and normal appearing tissue. In whole-brain analysis, MIML is 22 to 4980 times faster than the non-parametric and parametric methods, respectively.

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Normal volumetric and T1 relaxation time values at 1.5 T in segmented pediatric brain MRI using a MP2RAGE acquisition

OBJECTIVES

This study introduced a tailored MP2RAGE-based brain acquisition for a comprehensive assessment of the normal maturing brain.

METHODS

Seventy normal patients (35 girls and 35 boys) from 1 to 16 years of age were recruited within a prospective monocentric study conducted from a single University Hospital. Brain MRI examinations were performed at 1.5 T using a 20-channel head coil and an optimized 3D MP2RAGE sequence with a total acquisition time of 6:36 min. Automated 38 region segmentation was performed using the MorphoBox (template registration, bias field correction, brain extraction, and tissue classification) which underwent a major adaptation of three age-group T1-weighted templates. Volumetry and T1 relaxometry reference ranges were established using a logarithmic model and a modified Gompertz growth respectively.

RESULTS

Detailed automated brain segmentation and T1 mapping were successful in all patients. Using these data, an age-dependent model of normal brain maturation with respect to changes in volume and T1 relaxometry was established. After an initial rapid increase until 24 months of life, the total intracranial volume was found to converge towards 1400 mL during adolescence. The expected volumes of white matter (WM) and cortical gray matter (GM) showed a similar trend with age. After an initial major decrease, T1 relaxation times were observed to decrease progressively in all brain structures. The T1 drop in the first year of life was more pronounced in WM (from 1000-1100 to 650-700 ms) than in GM structures.

CONCLUSION

The 3D MP2RAGE sequence allowed to establish brain volume and T1 relaxation time normative ranges in pediatrics.

KEY POINTS

• The 3D MP2RAGE sequence provided a reliable quantitative assessment of brain volumes and T1 relaxation times during childhood. • An age-dependent model of normal brain maturation was established. • The normative ranges enable an objective comparison to a normal cohort, which can be useful to further understand, describe, and identify neurodevelopmental disorders in children.

Fast and high-resolution myelin water imaging: Accelerating multi-echo GRASE with CAIPIRINHA

PURPOSE

Although several MRI methods have been explored to achieve in vivo myelin quantification, imaging the whole brain in clinically acceptable times and sufficiently high resolution remains challenging. To address this problem, this work investigates the acceleration of multi-echo T₂ acquisitions based on the multi-echo gradient and spin echo (GRASE) sequence using CAIPIRINHA undersampling and adapted k-space reordering patterns.

METHODS

A prototype multi-echo GRASE sequence supporting CAIPIRINHA parallel imaging was implemented. Multi-echo T₂ data were acquired from 12 volunteers using the implemented sequence (1.6 × 1.6 × 1.6 mm³ , 84 slices, acquisition time [TA] = 10:30 min) and a multi-echo spin echo (MESE) sequence as reference (1.6 × 1.6 × 3.2 mm³ , single-slice, TA = 5:41 min). Myelin water fraction (MWF) maps derived from both acquisitions were compared via correlation and Bland-Altman analyses. In addition, scan-rescan datasets were acquired to evaluate the repeatability of the derived maps.

RESULTS

Resulting maps from the MESE and multi-echo GRASE sequences were found to be correlated (r = 0.83). The Bland-Altman analysis revealed a mean bias of -0.2% (P = .24) with the limits of agreement ranging from -3.7% to 3.3%. The Pearson's correlation coefficient among MWF values obtained from the scan-rescan datasets was found to be 0.95 and the mean bias equal to 0.11% (P = .32), indicating good repeatability of the retrieved maps.

CONCLUSION

By combining a 3D multi-echo GRASE sequence with CAIPIRINHA sampling, whole-brain MWF maps were obtained in 10:30 min with 1.6 mm isotropic resolution. The good correlation with conventional MESE-based maps demonstrates that the implemented sequence may be a promising alternative to time-consuming MESE acquisitions.

Model-based super-resolution reconstruction of T2 maps

PURPOSE

High-resolution isotropic T₂ mapping of the human brain with multi-echo spin-echo (MESE) acquisitions is challenging. When using a 2D sequence, the resolution is limited by the slice thickness. If used as a 3D acquisition, specific absorption rate limits are easily exceeded due to the high power deposition of nonselective refocusing pulses. A method to reconstruct 1-mm³ isotropic T₂ maps is proposed based on multiple 2D MESE acquisitions. Data were undersampled (10-fold) to compensate for the prolonged scan time stemming from the super-resolution acquisition.

THEORY AND METHODS

The proposed method integrates a classical super-resolution with an iterative model-based approach to reconstruct quantitative maps from a set of undersampled low-resolution data. The method was tested on numerical and multipurpose phantoms, and in vivo data. T₂ values were assessed with a region-of-interest analysis using a single-slice spin-echo and a fully sampled MESE acquisition in a phantom, and a MESE acquisition in healthy volunteers.

RESULTS

Numerical simulations showed that the best trade-off between acceleration and number of low-resolution datasets is 10-fold acceleration with 4 acquisitions (acquisition time = 18 min). The proposed approach showed improved resolution over low-resolution images for both phantom and brain. Region-of-interest analysis of the phantom compartments revealed that at shorter T₂ , the proposed method was comparable with the fully sampled MESE. For the volunteer data, the T₂ values found in the brain structures were consistent across subjects (8.5-13.1 ms standard deviation).

CONCLUSION

The proposed method addresses the inherent limitations associated with high-resolution T₂ mapping and enables the reconstruction of 1 mm³ isotropic relaxation maps with a 10 times faster acquisition.

Quantitative brain relaxation atlases for personalized detection and characterization of brain pathology

PURPOSE

To exploit the improved comparability and hardware independency of quantitative MRI, databases of MR physical parameters in healthy tissue are required, to which tissue properties of patients can be compared. In this work, normative values for longitudinal and transverse relaxation times in the brain were established and tested in single-subject comparisons for detection of abnormal relaxation times.

METHODS

Relaxometry maps of the brain were acquired from 52 healthy volunteers. After spatially normalizing the volumes into a common space, T₁ and T₂ inter-subject variability within the healthy cohort was modeled voxel-wise. A method for a single-subject comparison against the atlases was developed by computing z-scores with respect to the established healthy norms. The comparison was applied to two multiple sclerosis and one clinically isolated syndrome cases for a proof of concept.

RESULTS

The established atlases exhibit a low variation in white matter structures (median RMSE of models equal to 32 ms for T₁ and 4 ms for T₂ ), indicating that relaxation times are in a narrow range for normal tissues. The proposed method for single-subject comparison detected relaxation time deviations from healthy norms in the example patient data sets. Relaxation times were found to be increased in brain lesions (mean z-scores >5). Moreover, subtle and confluent differences (z-scores ~2-4) were observed in clinically plausible regions (between lesions, corpus callosum).

CONCLUSIONS

Brain T₁ and T₂ quantitative norms were derived voxel-wise with low variability in healthy tissue. Example patient deviation maps demonstrated good sensitivity of the atlases for detecting relaxation time alterations.