Investigating the effect of oblique image acquisition on the accuracy of QSM and a robust tilt correction method.

PURPOSE: Quantitative susceptibility mapping (QSM) is used increasingly for clinical research where oblique image acquisition is commonplace, but its effects on QSM accuracy are not well understood. THEORY AND METHODS: The QSM processing pipeline involves defining the unit magnetic dipole kernel, which requires knowledge of the direction of the main magnetic field B ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ with respect to the acquired image volume axes. The direction of B ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ is dependent on the axis and angle of rotation in oblique acquisition. Using both a numerical brain phantom and in vivo acquisitions in 5 healthy volunteers, we analyzed the effects of oblique acquisition on magnetic susceptibility maps. We compared three tilt-correction schemes at each step in the QSM pipeline: phase unwrapping, background field removal and susceptibility calculation, using the RMS error and QSM-tuned structural similarity index. RESULTS: Rotation of wrapped phase images gave severe artifacts. Background field removal with projection onto dipole fields gave the most accurate susceptibilities when the field map was first rotated into alignment with B ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ . Laplacian boundary value and variable-kernel sophisticated harmonic artifact reduction for phase data background field removal methods gave accurate results without tilt correction. For susceptibility calculation, thresholded k-space division, iterative Tikhonov regularization, and weighted linear total variation regularization, all performed most accurately when local field maps were rotated into alignment with B ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ before susceptibility calculation. CONCLUSION: For accurate QSM, oblique acquisition must be taken into account. Rotation of images into alignment with B ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ should be carried out after phase unwrapping and before background-field removal. We provide open-source tilt-correction code to incorporate easily into existing pipelines: https://github.com/o-snow/QSM_TiltCorrection.git.

Multi-echo quantitative susceptibility mapping: how to combine echoes for accuracy and precision at 3 Tesla.

PURPOSE: To compare different multi-echo combination methods for MRI QSM. Given the current lack of consensus, we aimed to elucidate how to optimally combine multi-echo gradient-recalled echo signal phase information, either before or after applying Laplacian-base methods (LBMs) for phase unwrapping or background field removal. METHODS: Multi-echo gradient-recalled echo data were simulated in a numerical head phantom, and multi-echo gradient-recalled echo images were acquired at 3 Tesla in 10 healthy volunteers. To enable image-based estimation of gradient-recalled echo signal noise, 5 volunteers were scanned twice in the same session without repositioning. Five QSM processing pipelines were designed: 1 applied nonlinear phase fitting over TEs before LBMs; 2 applied LBMs to the TE-dependent phase and then combined multiple TEs via either TE-weighted or SNR-weighted averaging; and 2 calculated TE-dependent susceptibility maps via either multi-step or single-step QSM and then combined multiple TEs via magnitude-weighted averaging. Results from different pipelines were compared using visual inspection; summary statistics of susceptibility in deep gray matter, white matter, and venous regions; phase noise maps (error propagation theory); and, in the healthy volunteers, regional fixed bias analysis (Bland-Altman) and regional differences between the means (nonparametric tests). RESULTS: Nonlinearly fitting the multi-echo phase over TEs before applying LBMs provided the highest regional accuracy of χ $$ \chi $$ and the lowest phase noise propagation compared to averaging the LBM-processed TE-dependent phase. This result was especially pertinent in high-susceptibility venous regions. CONCLUSION: For multi-echo QSM, we recommend combining the signal phase by nonlinear fitting before applying LBMs.

An optimized and highly repeatable MRI acquisition and processing pipeline for quantitative susceptibility mapping in the head-and-neck region.

PURPOSE: Quantitative Susceptibility Mapping (QSM) is an emerging technique sensitive to disease-related changes including oxygenation. It is extensively used in brain studies and has increasing clinical applications outside the brain. Here we present the first MRI acquisition protocol and QSM pipeline optimized for the head-and-neck region together with a repeatability analysis performed in healthy volunteers. METHODS: We investigated both the intrasession and the intersession repeatability of the optimized method in 10 subjects. We also implemented two, Tikhonov-regularisation-based susceptibility calculation techniques that were found to have higher contrast-to-noise than existing methods in the head-and-neck region. Repeatability was evaluated by calculating the distributions of susceptibility differences between repeated scans and the corresponding minimum detectable effect sizes (MDEs). RESULTS: Deep brain regions had higher QSM repeatability than neck regions. As expected, intrasession repeatability was generally better than intersession repeatability. Susceptibility maps calculated using projection onto dipole fields for background field removal were more repeatable than using the Laplacian boundary value method in the head-and-neck region. Small (short-axis diameter <5 mm) lymph nodes had the lowest repeatability (MDE = 0.27 ppm) as imperfect segmentation included some of the surrounding paramagnetic fatty fascia, highlighting the importance of accurate region delineation. MDEs in the larger lymph nodes (0.16 ppm), submandibular glands (0.10 ppm), and especially the parotid glands (0.06 ppm) were much lower, comparable to those of the brain regions. CONCLUSIONS: The high repeatability of the acquisition and pipeline optimized for QSM will facilitate clinical studies in the head-and-neck region.

Investigating the accuracy and precision of TE-dependent versus multi-echo QSM using Laplacian-based methods at 3 T.

PURPOSE: Multi-echo gradient-recalled echo acquisitions for QSM enable optimizing the SNR for several tissue types through multi-echo (TE) combination or investigating temporal variations in the susceptibility (potentially reflecting tissue microstructure) by calculating one QSM image at each TE (TE-dependent QSM). In contrast with multi-echo QSM, applying Laplacian-based methods (LBMs) for phase unwrapping and background field removal to single TEs could introduce nonlinear temporal variations (independent of tissue microstructure) into the measured susceptibility. Here, we aimed to compare the effect of LBMs on the QSM susceptibilities in TE-dependent versus multi-echo QSM. METHODS: TE-dependent recalled echo data simulated in a numerical head phantom and gradient-recalled echo images acquired at 3 T in 10 healthy volunteers. Several QSM pipelines were tested, including four distinct LBMs: sophisticated harmonic artifact reduction for phase data (SHARP), variable-radius sophisticated harmonic artifact reduction for phase data (V-SHARP), Laplacian boundary value background field removal (LBV), and one-step total generalized variation (TGV). Results from distinct pipelines were compared using visual inspection, summary statistics of susceptibility in deep gray matter/white matter/venous regions of interest, and, in the healthy volunteers, regional susceptibility bias analysis and nonparametric tests. RESULTS: Multi-echo versus TE-dependent QSM had higher regional accuracy, especially in high-susceptibility regions and at shorter TEs. Everywhere except in the veins, a processing pipeline incorporating TGV provided the most temporally stable TE-dependent QSM results with an accuracy similar to multi-echo QSM. CONCLUSIONS: For TE-dependent QSM, carefully choosing LBMs can minimize the introduction of LBM-related nonlinear temporal susceptibility variations.

The effect of low resolution and coverage on the accuracy of susceptibility mapping.

PURPOSE: Quantitative susceptibility mapping (QSM) has found increasing clinical applications. However, to reduce scan time, clinical acquisitions often use reduced resolution and coverage, particularly in the through-slice dimension. The effect of these factors on QSM has begun to be assessed using only balloon phantoms and downsampled brain images. Here, we investigate the effects (and their sources) of low resolution or coverage on QSM using both simulated and acquired images. METHODS: Brain images were acquired at 1 mm isotropic resolution and full brain coverage, and low resolution (up to 6 mm slice thickness) or coverage (down to 20 mm) in 5 healthy volunteers. Images at reduced resolution or coverage were also simulated in these volunteers and in a new, anthropomorphic, numerical phantom. Mean susceptibilities in 5 brain regions, including white matter, were investigated over varying resolution and coverage. RESULTS: The susceptibility map contrast decreased with increasing slice thickness and spacing, and with decreasing coverage below ~40 mm for 2 different QSM pipelines. Our simulations showed that calculated susceptibility values were erroneous at low resolution or very low coverage, because of insufficient sampling and overattenuation of the susceptibility-induced field perturbations. Susceptibility maps calculated from simulated and acquired images showed similar behavior. CONCLUSIONS: Both low resolution and low coverage lead to loss of contrast and errors in susceptibility maps. The widespread clinical practice of using low resolution and coverage does not provide accurate susceptibility maps. Simulations in images of healthy volunteers and in a new, anthropomorphic numerical phantom were able to accurately model low-resolution and low-coverage acquisitions.

Association of bone mineral density and fat fraction with magnetic susceptibility in inflamed trabecular bone.

PURPOSE: To evaluate the relationship between bone mineral density (BMD) and magnetic susceptibility, and between proton density fat fraction and susceptibility, in inflamed trabecular bone. METHODS: Two different phantoms modeling the fat fraction (FF) and BMD values of healthy bone marrow and disease states were scanned using a multiecho gradient echo acquisition at 3T. After correction for fat-water chemical shift, susceptibility mapping was performed, and susceptibility measurements were compared with BMD and FF values using linear regression. Patients with spondyloarthritis were scanned using the same protocol, and susceptibility values were calculated in areas of inflamed bone (edema) and fat metaplasia, both before and after accounting for the contribution of fat to the total susceptibility. RESULTS: Susceptibility values in the phantoms were accurately described by a 2D linear function, with a negative correlation between BMD and susceptibility and a positive correlation between FF and susceptibility (adjusted R2 = 0.77; P = 3·10-5 ). In patients, significant differences in susceptibility were observed between fat metaplasia and normal marrow, but these differences were eliminated by removing the fat contribution to the total susceptibility. CONCLUSIONS: BMD and proton density fat fraction both influence the total susceptibility of bone marrow and failure to account for the fat contribution could lead to errors in BMD quantification. We propose a method for removing the fat contribution from the total susceptibility, based on the observed linear relationship between susceptibility and FF. In inflamed bone, the overall increase in susceptibility in areas of fat metaplasia is at least partly due to increased fat content.

A comprehensive protocol for quantitative magnetic resonance imaging of the brain at 3 Tesla.

Quantitative MRI (qMRI) has been shown to be clinically useful for numerous applications in the brain and body. The development of rapid, accurate, and reproducible qMRI techniques offers access to new multiparametric data, which can provide a comprehensive view of tissue pathology. This work introduces a multiparametric qMRI protocol along with full postprocessing pipelines, optimized for brain imaging at 3 Tesla and using state-of-the-art qMRI tools. The total scan time is under 50 minutes and includes eight pulse-sequences, which produce range of quantitative maps including T1, T2, and T2* relaxation times, magnetic susceptibility, water and macromolecular tissue fractions, mean diffusivity and fractional anisotropy, magnetization transfer ratio (MTR), and inhomogeneous MTR. Practical tips and limitations of using the protocol are also provided and discussed. Application of the protocol is presented on a cohort of 28 healthy volunteers and 12 brain regions-of-interest (ROIs). Quantitative values agreed with previously reported values. Statistical analysis revealed low variability of qMRI parameters across subjects, which, compared to intra-ROI variability, was x4.1 ± 0.9 times higher on average. Significant and positive linear relationship was found between right and left hemispheres' values for all parameters and ROIs with Pearson correlation coefficients of r>0.89 (P<0.001), and mean slope of 0.95 ± 0.04. Finally, scan-rescan stability demonstrated high reproducibility of the measured parameters across ROIs and volunteers, with close-to-zero mean difference and without correlation between the mean and difference values (across map types, mean P value was 0.48 ± 0.27). The entire quantitative data and postprocessing scripts described in the manuscript are publicly available under dedicated GitHub and Figshare repositories. The quantitative maps produced by the presented protocol can promote longitudinal and multi-center studies, and improve the biological interpretability of qMRI by integrating multiple metrics that can reveal information, which is not apparent when examined using only a single contrast mechanism.

Investigating the relationship between thalamic iron concentration and disease severity in secondary progressive multiple sclerosis using quantitative susceptibility mapping: Cross-sectional analysis from the MS-STAT2 randomised controlled trial.

Background: Deep grey matter pathology is a key driver of disability worsening in people with multiple sclerosis. Quantitative susceptibility mapping (QSM) is an advanced magnetic resonance imaging (MRI) technique which quantifies local magnetic susceptibility from variations in phase produced by changes in the local magnetic field. In the deep grey matter, susceptibility has previously been validated against tissue iron concentration. However, it currently remains unknown whether susceptibility is abnormal in older progressive MS cohorts, and whether it correlates with disability. Objectives: To investigate differences in mean regional susceptibility in deep grey matter between people with secondary progressive multiple sclerosis (SPMS) and healthy controls; to examine in patients the relationships between deep grey matter susceptibility and clinical and imaging measures of disease severity. Methods: Baseline data from a subgroup of the MS-STAT2 trial (simvastatin vs. placebo in SPMS, NCT03387670) were included. The subgroup underwent clinical assessments and an advanced MRI protocol at 3T. A cohort of age-matched healthy controls underwent the same MRI protocol. Susceptibility maps were reconstructed using a robust QSM pipeline from multi-echo 3D gradient-echo sequence. Regions of interest (ROIs) in the thalamus, globus pallidus and putamen were segmented from 3D T1-weighted images, and lesions segmented from 3D fluid-attenuated inversion recovery images. Linear regression was used to compare susceptibility from ROIs between patients and controls, adjusting for age and sex. Where significant differences were found, we further examined the associations between ROI susceptibility and clinical and imaging measures of MS severity. Results: 149 SPMS (77% female; mean age: 53 yrs; median Expanded Disability Status Scale (EDSS): 6.0 [interquartile range 4.5-6.0]) and 33 controls (52% female, mean age: 57) were included.Thalamic susceptibility was significantly lower in SPMS compared to controls: mean (SD) 28.6 (12.8) parts per billion (ppb) in SPMS vs. 39.2 (12.7) ppb in controls; regression coefficient: -12.0 [95% confidence interval: -17.0 to -7.1], p < 0.001. In contrast, globus pallidus and putamen susceptibility were similar between both groups.In SPMS, a 10 ppb lower thalamic susceptibility was associated with a +0.13 [+0.01 to +0.24] point higher EDSS (p < 0.05), a -2.4 [-3.8 to -1.0] point lower symbol digit modality test (SDMT, p = 0.001), and a -2.4 [-3.7 to -1.1] point lower Sloan low contrast acuity, 2.5% (p < 0.01).Lower thalamic susceptibility was also strongly associated with a higher T2 lesion volume (T2LV, p < 0.001) and lower normalised whole brain, deep grey matter and thalamic volumes (all p < 0.001). Conclusions: The reduced thalamic susceptibility found in SPMS compared to controls suggests that thalamic iron concentrations are lower at this advanced stage of the disease. The observed relationships between lower thalamic susceptibility and more severe physical, cognitive and visual disability suggests that reductions in thalamic iron may correlate with important mechanisms of clinical disease progression. Such mechanisms appear to intimately link reductions in thalamic iron with higher T2LV and the development of thalamic atrophy, encouraging further research into QSM-derived thalamic susceptibility as a biomarker of disease severity in SPMS.

Feasibility of smartphone colorimetry of the face as an anaemia screening tool for infants and young children in Ghana.

BACKGROUND: Anaemia affects approximately a quarter of the global population. When anaemia occurs during childhood, it can increase susceptibility to infectious diseases and impair cognitive development. This research uses smartphone-based colorimetry to develop a non-invasive technique for screening for anaemia in a previously understudied population of infants and young children in Ghana. METHODS: We propose a colorimetric algorithm for screening for anaemia which uses a novel combination of three regions of interest: the lower eyelid (palpebral conjunctiva), the sclera, and the mucosal membrane adjacent to the lower lip. These regions are chosen to have minimal skin pigmentation occluding the blood chromaticity. As part of the algorithm development, different methods were compared for (1) accounting for varying ambient lighting, and (2) choosing a chromaticity metric for each region of interest. In comparison to some prior work, no specialist hardware (such as a colour reference card) is required for image acquisition. RESULTS: Sixty-two patients under 4 years of age were recruited as a convenience clinical sample in Korle Bu Teaching Hospital, Ghana. Forty-three of these had quality images for all regions of interest. Using a naïve Bayes classifier, this method was capable of screening for anaemia (<11.0g/dL haemoglobin concentration) vs healthy blood haemoglobin concentration (≥11.0g/dL) with a sensitivity of 92.9% (95% CI 66.1% to 99.8%), a specificity of 89.7% (72.7% to 97.8%) when acting on unseen data, using only an affordable smartphone and no additional hardware. CONCLUSION: These results add to the body of evidence suggesting that smartphone colorimetry is likely to be a useful tool for making anaemia screening more widely available. However, there remains no consensus on the optimal method for image preprocessing or feature extraction, especially across diverse patient populations.

SEGUE: A Speedy rEgion-Growing Algorithm for Unwrapping Estimated Phase.

Recent magnetic resonance imaging (MRI) techniques, such as quantitative magnetic susceptibility mapping, employ the signal phase to reveal disease-related changes in tissue composition, including iron or calcium content. The MRI phase is also routinely used in functional and diffusion MRI for distortion correction. However, phase images are wrapped into a range of 2π radians. Phase region expanding labeller for unwrapping discrete estimates (PRELUDE) is the gold standard method for robust, spatial, 3-D, MRI phase unwrapping. Unfortunately, PRELUDE's computation time can reach 15 min for a severely wrapped brain image and nearly 10 h to unwrap a full head-and-neck image on a standard PC. In this paper, we develop a Speedy rEgion-Growing algorithm for Unwrapping Estimated phase (SEGUE) based on similar principles to PRELUDE, implemented with additional methods for acceleration. We compared PRELUDE and SEGUE in numerical phantoms, and using in vivo images of the brain, head and neck, and pelvis acquired in 4-5 healthy volunteers and at 4-6 echo times. To overcome chemical-shift-induced errors within the head and neck, and pelvic images, we also investigated applying both techniques within fat and water masks separately. SEGUE provided almost identical unwrapped phase maps to the gold standard PRELUDE. SEGUE was (1.5 to 70 times) faster than PRELUDE, especially in severely wrapped images at later echoes and in the head and neck, and pelvic images. Applying these techniques within fat and water masks separately removed chemical-shift-induced errors successfully. SEGUE's MATLAB implementation is available for download. SEGUE is a general unwrapping algorithm not specific to MRI, and therefore could be used in images acquired with other modalities.