Investigating the effect of flow compensation and quantitative susceptibility mapping method on the accuracy of venous susceptibility measurement

Comparing repeatability metrics for quantitative susceptibility mapping in the head and neck

OBJECTIVE

Quantitative susceptibility mapping (QSM) is a technique that has been demonstrated to be highly repeatable in the brain. As QSM is applied to other parts of the body, it is necessary to investigate metrics for quantifying repeatability, to enable optimization of repeatable QSM reconstruction pipelines beyond the brain.

MATERIALS AND METHODS

MRI data were acquired in the head and neck (HN) region in ten healthy volunteers, who underwent six acquisitions across two sessions. QSMs were reconstructed using six representative state-of-the-art techniques. Repeatability of the susceptibility values was compared using voxel-wise metrics (normalized root mean squared error and XSIM) and ROI-based metrics (within-subject and between-subject standard deviation, coefficient of variation (CV), intraclass correlation coefficient (ICC)).

RESULTS

Both within-subject and between-subject variations were smaller than the variation between QSM dipole inversion methods, in most ROIs. autoNDI produced the most repeatable susceptibility values, with ICC > 0.75 in three of six HN ROIs with an average ICC of 0.66 across all ROIs. Joint consideration of standard deviation and ICC offered the best metric of repeatability for comparisons between QSM methods, given typical distributions of positive and negative QSM values.

DISCUSSION

Repeatability of QSM in the HN region is highly dependent on the dipole inversion method chosen, but the most repeatable methods (autoNDI, QSMnet, TFI) are only moderately repeatable in most HN ROIs.

Imaging arterial and venous vessels using Iron Dextran enhanced multi-echo 3D gradient echo MRI at 7T

Iron Dextran is a widely used iron oxide compound to treat iron-deficiency anemia patients in the clinic. Similar to other iron oxide compounds such as Ferumoxytol, it can also be used off-label as an intravascular magnetic resonance imaging (MRI) contrast agent due to its strong iron-induced T2 and T2* shortening effects. In this study, we seek to evaluate the feasibility of using Iron Dextran enhanced multi-echo susceptibility weighted imaging (SWI) MRI at 7T to image arterial and venous blood vessels in the human brain. Phantom experiments were performed to measure the r2* relaxivity for Iron Dextran in blood, based on which the SWI sequence was optimized. Pre- and post-infusion MR images were acquired in human subjects from which maps of arteries and veins were extracted. The post-contrast SWI images showed enhanced susceptibility difference between blood and the surrounding tissue in both arteries and veins. Our results showed that the proposed Iron Dextran enhanced multi-echo SWI approach allowed the visualization of blood vessels with diameters down to ~100 μm, including small blood vessels supplying and draining small brain structures such as the hippocampus. We conclude that Iron Dextran can be an alternative iron-based MRI contrast agent for blood vessel imaging in the human brain.

In vivo mapping of hippocampal venous vasculature and oxygenation using susceptibility imaging at 7T

Mapping the small venous vasculature of the hippocampus in vivo is crucial for understanding how functional changes of hippocampus evolve with age. Oxygen utilization in the hippocampus could serve as a sensitive biomarker for early degenerative changes, surpassing hippocampal tissue atrophy as the main source of information regarding tissue degeneration. Using an ultrahigh field (7T) susceptibility-weighted imaging (SWI) sequence, it is possible to capture oxygen-level dependent contrast of submillimeter-sized vessels. Moreover, the quantitative susceptibility mapping (QSM) results derived from SWI data allow for the simultaneous estimation of venous oxygenation levels, thereby enhancing the understanding of hippocampal function. In this study, we proposed two potential imaging markers in a cohort of 19 healthy volunteers aged between 20 and 74 years. These markers were: 1) hippocampal venous density on SWI images and 2) venous susceptibility (Δχvein) in the hippocampus-associated draining veins (the inferior ventricular veins (IVV) and the basal veins of Rosenthal (BVR) using QSM images). They were chosen specifically to help characterize the oxygen utilization of the human hippocampus and medial temporal lobe (MTL). As part of the analysis, we demonstrated the feasibility of measuring hippocampal venous density and Δχvein in the IVV and BVR at 7T with high spatial resolution (0.25 × 0.25 × 1 mm3). Our results demonstrated the in vivo reconstruction of the hippocampal venous system, providing initial evidence regarding the presence of the venous arch structure within the hippocampus. Furthermore, we evaluated the age effect of the two quantitative estimates and observed a significant increase in Δχvein for the IVV with age (p=0.006, r2 = 0.369). This may suggest the potential application of Δχvein in IVV as a marker for assessing changes in atrophy-related hippocampal oxygen utilization in normal aging and neurodegenerative diseases such as AD and dementia.

Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group

This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.

Modeling venous bias in resting state functional MRI metrics

Resting-state (rs) functional magnetic resonance imaging (fMRI) is used to detect low-frequency fluctuations in the blood oxygen-level dependent (BOLD) signal across brain regions. Correlations between temporal BOLD signal fluctuations are commonly used to infer functional connectivity. However, because BOLD is based on the dilution of deoxyhemoglobin, it is sensitive to veins of all sizes, and its amplitude is biased by draining veins. These biases affect local BOLD signal location and amplitude, and may also influence BOLD-derived connectivity measures, but the magnitude of this venous bias and its relation to vein size and proximity is unknown. Here, veins were identified using high-resolution quantitative susceptibility maps and utilized in a biophysical model to investigate systematic venous biases on common local rsfMRI-derived measures. Specifically, we studied the impact of vein diameter and distance to veins on the amplitude of low-frequency fluctuations (ALFF), fractional ALFF (fALFF), Hurst exponent (HE), regional homogeneity (ReHo), and eigenvector centrality values in the grey matter. Values were higher across all distances in smaller veins, and decreased with increasing vein diameter. Additionally, rsfMRI values associated with larger veins decrease with increasing distance from the veins. ALFF and ReHo were the most biased by veins, while HE and fALFF exhibited the smallest bias. Across all metrics, the amplitude of the bias was limited in voxel-wise data, confirming that venous structure is not the dominant source of contrast in these rsfMRI metrics. Finally, the models presented can be used to correct this venous bias in rsfMRI metrics.

Cerebral oxygen metabolism from MRI susceptibility

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

Multimodal comparisons of QSM and PET in neurodegeneration and aging

Quantitative susceptibility mapping (QSM) has been used to study susceptibility changes that may occur based on tissue composition and mineral deposition. Iron is a primary contributor to changes in magnetic susceptibility and of particular interest in applications of QSM to neurodegeneration and aging. Iron can contribute to neurodegeneration through inflammatory processes and via interaction with aggregation of disease-related proteins. To better understand the local susceptibility changes observed on QSM, its signal has been studied in association with other imaging metrics such as positron emission tomography (PET). The associations of QSM and PET may provide insight into the pathophysiology of disease processes, such as the role of iron in aging and neurodegeneration, and help to determine the diagnostic utility of QSM as an indirect indicator of disease processes typically evaluated with PET. In this review we discuss the proposed mechanisms and summarize prior studies of the associations of QSM and amyloid PET, tau PET, TSPO PET, FDG-PET, 15O-PET, and F-DOPA PET in evaluation of neurologic diseases with a focus on aging and neurodegeneration.

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 fieldB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ with respect to the acquired image volume axes. The direction ofB ^ 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 withB ^ 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 withB ^ 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 withB ^ 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.

Sub-acute Changes on MRI Measures of Cerebral Blood Flow and Venous Oxygen Saturation in Concussed Australian Rules Footballers

Background

Sports-related concussion (SRC) is common in collision sport athletes. There is growing evidence that repetitive SRC can have serious neurological consequences, particularly when the repetitive injuries occur when the brain has yet to fully recover from the initial injury. Hence, there is a need to identify biomarkers that are capable of determining SRC recovery so that they can guide clinical decisions pertaining to return-to-play. Cerebral venous oxygen saturation (SvO₂) and cerebral blood flow (CBF) can be measured using magnetic resonance imaging (MRI) and may provide insights into changing energy demands and recovery following SRC.

Results

In this study we therefore investigated SvO₂ and CBF in a cohort of concussed amateur Australian Football athletes (i.e., Australia’s most participated collision sport). Male and female Australian footballers (n = 13) underwent MRI after being cleared to return to play following a mandatory 13-day recovery period and were compared to a group of control Australian footballers (n = 16) with no recent history of SRC (i.e., > 3 months since last SRC). Despite the concussed Australian footballers being cleared to return to play at the time of MRI, we found evidence of significantly increased susceptibility in the global white matter (p = 0.020) and a trend (F_5,21 = 2.404, p = 0.071) for reduced relative CBF (relCBF) compared to the control group. Further, there was evidence of an interaction between sex and injury in straight sinus susceptibility values (F_1,25 = 3.858, p = 0.061) which were decreased in female SRC athletes (p = 0.053). Of note, there were significant negative correlations between straight sinus susceptibility and relCBF suggesting impaired metabolic function after SRC.

Conclusions

These findings support the use of quantitative susceptibility mapping (QSM) and relCBF as sensitive indicators of SRC, and raise further concerns related to SRC guidelines that allow for return-to-play in less than two weeks.

Multisite reproducibility of quantitative susceptibility mapping and effective transverse relaxation rate in deep gray matter at 3 T using locally optimized sequences in 24 traveling heads

Iron concentration in the human brain plays a crucial role in several neurodegenerative diseases and can be monitored noninvasively using quantitative susceptibility mapping (QSM) and effective transverse relaxation rate (R₂ *) mapping from multiecho T₂ *-weighted images. Large population studies enable better understanding of pathologies and can benefit from pooling multisite data. However, reproducibility may be compromised between sites and studies using different hardware and sequence protocols. This work investigates QSM and R₂ * reproducibility at 3 T using locally optimized sequences from three centers and two vendors, and investigates possible reduction of cross-site variability through postprocessing approaches. Twenty-four healthy subjects traveled between three sites and were scanned twice at each site. Scan-rescan measurements from seven deep gray matter regions were used for assessing within-site and cross-site reproducibility using intraclass correlation coefficient (ICC) and within-subject standard deviation (SDw) measures. In addition, multiple QSM and R₂ * postprocessing options were investigated with the aim to minimize cross-site sequence-related variations, including: mask generation approach, echo-timing selection, harmonizing spatial resolution, field map estimation, susceptibility inversion method, and linear field correction for magnitude images. The same-subject cross-site region of interest measurements for QSM and R₂ * were highly correlated (R²  ≥ 0.94) and reproducible (mean ICC of 0.89 and 0.82 for QSM and R₂ *, respectively). The mean cross-site SDw was 4.16 parts per billion (ppb) for QSM and 1.27 s^-1 for R₂ *. For within-site measurements of QSM and R₂ *, the mean ICC was 0.97 and 0.87 and mean SDw was 2.36 ppb and 0.97 s^-1 , respectively. The precision level is regionally dependent and is reduced in the frontal lobe, near brain edges, and in white matter regions. Cross-site QSM variability (mean SDw) was reduced up to 46% through postprocessing approaches, such as masking out less reliable regions, matching available echo timings and spatial resolution, avoiding the use of the nonconsistent magnitude contrast between scans in field estimation, and minimizing streaking artifacts.

Reduced white matter venous density on MRI is associated with neurodegeneration and cognitive impairment in the elderly

High-resolution susceptibility weighted imaging (SWI) provides unique contrast to small venous vasculature. The conspicuity of these mesoscopic veins, such as deep medullary veins in white matter, is subject to change from SWI venography when venous oxygenation in these veins is altered due to oxygenated blood susceptibility changes. The changes of visualization in small veins shows potential to depict regional changes of oxygen utilization and/or vascular density changes in the aging brain. The goal of this study was to use WM venous density to quantify small vein visibility in WM and investigate its relationship with neurodegenerative features, white matter hyperintensities (WMHs), and cognitive/functional status in elderly subjects (N = 137). WM venous density was significantly associated with neurodegeneration characterized by brain atrophy (β = 0.046± 0.01, p < 0.001), but no significant association was found between WM venous density and WMHs lesion load (p = 0.3963). Further analysis of clinical features revealed a negative trend of WM venous density with the sum-of-boxes of Clinical Dementia Rating and a significant association with category fluency (1-min animal naming). These results suggest that WM venous density on SWI can be used as a sensitive marker to characterize cerebral oxygen metabolism and different stages of cognitive and functional status in neurodegenerative diseases.

A Comparison of MRI Quantitative Susceptibility Mapping and TRUST-Based Measures of Brain Venous Oxygen Saturation in Sickle Cell Anaemia

In recent years, interest has grown in the potential for magnetic resonance imaging (MRI) measures of venous oxygen saturation (Yv) to improve neurological risk prediction. T2-relaxation-under-spin-tagging (TRUST) is an MRI technique which has revealed changes in Yv in patients with sickle cell anemia (SCA). However, prior studies comparing Yv in patients with SCA relative to healthy controls have reported opposing results depending on whether the calibration model, developed to convert blood T2 to Yv, is based on healthy human hemoglobin (HbA), bovine hemoglobin (HbBV) or sickle hemoglobin (HbS). MRI Quantitative Susceptibility Mapping (QSM) is an alternative technique that may hold promise for estimating Yv in SCA as blood magnetic susceptibility is linearly dependent upon Yv, and no significant difference has been found between the magnetic susceptibility of HbA and HbS. Therefore, the aim of this study was to compare estimates of Yv using QSM and TRUST with five published calibration models in healthy controls and patients with SCA. 17 patients with SCA and 13 healthy controls underwent MRI. Susceptibility maps were calculated from a multi-parametric mapping acquisition and Yv was calculated from the mean susceptibility in a region of interest in the superior sagittal sinus. TRUST estimates of T2, within a similar but much smaller region, were converted to Yv using five different calibration models. Correlation and Bland-Altman analyses were performed to compare estimates of Yv between TRUST and QSM methods. For each method, t-tests were also used to explore group-wise differences between patients with SCA and healthy controls. In healthy controls, significant correlations were observed between QSM and TRUST measures of Yv, while in SCA, there were no such correlations. The magnitude and direction of group-wise differences in Yv varied with method. The TRUST-HbBV and QSM methods suggested decreased Yv in SCA relative to healthy controls, while the TRUST-HbS (p < 0.01) and TRUST-HbA models suggested increased Yv in SCA as in previous studies. Further validation of all MRI measures of Yv, relative to ground truth measures such as O15 PET and jugular vein catheterization, is required in SCA before QSM or TRUST methods can be considered for neurological risk prediction.

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 χ 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.

Click here for author‐reader discussions