Susceptibility mapping as a means to visualize veins and quantify oxygen saturation

Vessel-specific quantification of cerebral venous oxygenation with velocity-encoding preparation and rapid acquisition

PURPOSE

Non-invasive measurement of cerebral venous oxygenation (Yv) is of critical importance in brain diseases. The present work proposed a fast method to quantify regional Yv map for both large and small veins.

METHODS

A new sequence was developed, referred to as TRU-VERA (T2 relaxation under velocity encoding and rapid acquisition, which isolates blood spins from static tissue with velocity-encoding preparation, modulates the T2 weighting of venous signal with T2-preparation and utilizes a bSSFP readout to achieve fast acquisition with high resolution. The sequence was first optimized to achieve best sensitivity for both large and small veins, and then validated with TRUST (T2 relaxation under spin tagging), TRUPC (T2 relaxation under phase contrast), and accelerated TRUPC MRI. Regional difference of Yv was evaluated, and test-retest reproducibility was examined.

RESULTS

Optimal Venc was determined to be 3 cm/s, while recovery time and balanced SSFP flip angle within reasonable range had minimal effect on SNR efficiency. Venous T2 measured with TRU-VERA was highly correlated with T2 from TRUST (R2 = 0.90), and a conversion equation was established for further calibration to Yv. TRU-VERA sequences showed consistent Yv estimation with TRUPC (R2 = 0.64) and accelerated TRUPC (R2 = 0.79). Coefficient of variation was 0.84% for large veins and 2.49% for small veins, suggesting an excellent test-retest reproducibility.

CONCLUSION

The proposed TRU-VERA sequence is a promising method for vessel-specific oxygenation assessment.

Unveiling metabolic patterns in dementia: Insights from high-resolution quantitative blood-oxygenation-level-dependent MRI

BACKGROUND

Oxygen extraction fraction (OEF) and deoxyhemoglobin (DoHb) levels reflect variations in cerebral oxygen metabolism in demented patients.

PURPOSE

Delineating the metabolic profiles evident throughout different phases of dementia necessitates an integrated analysis of OEF and DoHb levels. This is enabled by leveraging high-resolution quantitative blood oxygenation level dependent (qBOLD) analysis of magnitude images obtained from a multi-echo gradient-echo MRI (mGRE) scan performed on a 3.0 Tesla scanner.

METHODS

Achieving superior spatial resolution in qBOLD necessitates the utilization of an mGRE scan with only four echoes, which in turn limits the number of measurements compared to the parameters within the qBOLD model. Consequently, it becomes imperative to discard non-essential parameters to facilitate further analysis. This process entails transforming the qBOLD model into a format suitable for fitting the log-magnitude difference (L-MDif) profiles of the four echo magnitudes present in each brain voxel. In order to bolster spatial specificity, the log-difference qBOLD model undergoes refinement into a representative form, termed as r-qBOLD, particularly when applied to class-averaged L-MDif signals derived through k-means clustering of L-MDif signals from all brain voxels into a predetermined number of clusters. The agreement between parameters estimated using r-qBOLD for different cluster sizes is validated using Bland-Altman analysis, and the model's goodness-of-fit is evaluated using aχ 2 ${\chi ^2}$ -test. Retrospective MRI data of Alzheimer's disease (AD), mild cognitive impairment (MCI), and non-demented patients without neuropathological disorders, pacemakers, other implants, or psychiatric disorders, who completed a minimum of three visits prior to MRI enrolment, are utilized for the study.

RESULTS

Utilizing a cohort comprising 30 demented patients aged 65-83 years in stages 4-6 representing mild, moderate, and severe stages according to the clinical dementia rating (CDR), matched with an age-matched non-demented control group of 18 individuals, we conducted joint observations of OEF and DoHb levels estimated using r-qBOLD. The observations elucidate metabolic signatures in dementia based on OEF and DoHb levels in each voxel. Our principal findings highlight the significance of spatial patterns of metabolic profiles (metabolic patterns) within two distinct regimes: OEF levels exceeding the normal range (S1-regime), and OEF levels below the normal range (S2-regime). The S1-regime, accompanied by low DoHb levels, predominantly manifests in fronto-parietal and perivascular regions with increase in dementia severity. Conversely, the S2-regime, accompanied by low DoHb levels, is observed in medial temporal (MTL) regions. Other regions with abnormal metabolic patterns included the orbitofrontal cortex (OFC), medial-orbital prefrontal cortex (MOPFC), hypothalamus, ventro-medial prefrontal cortex (VMPFC), and retrosplenial cortex (RSP). Dysfunction in the OFC and MOPFC indicated cognitive and emotional impairment, while hypothalamic involvement potentially indicated preclinical dementia. Reduced metabolic activity in the RSP suggested early-stage AD related functional abnormalities.

CONCLUSIONS

Integrated analysis of OEF and DoHb levels using r-qBOLD reveals distinct metabolic signatures across dementia phases, highlighting regions susceptible to neuronal loss, vascular involvement, and preclinical indicators.

Reduced oxygen extraction fraction in deep cerebral veins associated with cognitive impairment in multiple sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89% vs 72.4 ± 2.23%) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.508, p = 0.031) and the SvO2 (r = -0.498, p = 0.036) exhibited a moderate correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS. .

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.

Strategically Acquired Gradient Echo (STAGE) Imaging, part IV: Constrained Reconstruction of White Noise (CROWN) Processing as a Means to Improve Signal-to-Noise in STAGE Imaging at 3 Tesla

Increasing the signal-to-noise ratio (SNR) has always been of critical importance for magnetic resonance imaging. Although increasing field strength provides a linear increase in SNR, it is more and more costly as field strength increases. Therefore, there is a major effort today to use signal processing methods to improve SNR since it is more efficient and economical. There are a variety of methods to improve SNR such as averaging the data at the expense of imaging time, or collecting the data with a lower resolution, all of these methods, including imaging processing methods, usually come at the expense of loss of image detail or image blurring. Therefore, we developed a new mathematical approach called CROWN (Constrained Reconstruction of White Noise) to enhance SNR without loss of structural detail and without affecting scanning time. In this study, we introduced and tested the concept behind CROWN specifically for STAGE (strategically acquired gradient echo) imaging. The concept itself is presented first, followed by simulations to demonstrate its theoretical effectiveness. Then the SNR improvement on proton spin density (PSD) and R2⁎ maps was investigated using brain STAGE data acquired from 10 healthy controls (HCs) and 10 patients with Parkinson's disease (PD). For the PSD and R2* maps, the SNR and CNR between white matter and gray matter were improved by a factor of 1.87 ± 0.50 and 1.72 ± 0.88, respectively. The white matter hyperintensity lesions in PD patients were more clearly defined after CROWN processing. Using these improved maps, simulated images for any repeat time, echo time or flip angle can be created with improved SNR. The potential applications of this technology are to trade off the increased SNR for higher resolution images and/or faster imaging.

Enhancing Nigrosome-1 Sign Identification via Interpretable AI using True Susceptibility Weighted Imaging

BACKGROUND

Nigrosome 1 (N1), the largest nigrosome region in the ventrolateral area of the substantia nigra pars compacta, is identifiable by the "N1 sign" in long echo time gradient echo MRI. The N1 sign's absence is a vital Parkinson's disease (PD) diagnostic marker. However, it is challenging to visualize and assess the N1 sign in clinical practice.

PURPOSE

To automatically detect the presence or absence of the N1 sign from true susceptibility weighted imaging by using deep-learning method.

STUDY TYPE

Prospective.

POPULATION/SUBJECTS

453 subjects, including 225 PD patients, 120 healthy controls (HCs), and 108 patients with other movement disorders, were prospectively recruited including 227 males and 226 females. They were divided into training, validation, and test cohorts of 289, 73, and 91 cases, respectively.

FIELD STRENGTH/SEQUENCE

3D gradient echo SWI sequence at 3T; 3D multiecho strategically acquired gradient echo imaging at 3T; NM-sensitive 3D gradient echo sequence with MTC pulse at 3T.

ASSESSMENT

A neuroradiologist with 5 years of experience manually delineated substantia nigra regions. Two raters with 2 and 36 years of experience assessed the N1 sign on true susceptibility weighted imaging (tSWI), QSM with high-pass filter, and magnitude data combined with MTC data. We proposed NINet, a neural model, for automatic N1 sign identification in tSWI images.

STATISTICAL TESTS

We compared the performance of NINet to the subjective reference standard using Receiver Operating Characteristic analyses, and a decision curve analysis assessed identification accuracy.

RESULTS

NINet achieved an area under the curve (AUC) of 0.87 (CI: 0.76-0.89) in N1 sign identification, surpassing other models and neuroradiologists. NINet localized the putative N1 sign within tSWI images with 67.3% accuracy.

DATA CONCLUSION

Our proposed NINet model's capability to determine the presence or absence of the N1 sign, along with its localization, holds promise for enhancing diagnostic accuracy when evaluating PD using MR images.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 1.

Reduced Oxygen Extraction Fraction in Deep Cerebral Veins Associated with Cognitive Impairment in Multiple Sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89 % vs 72.4 ± 2.23 %) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.646, p = 0.004) and the SvO2 (r = -0.603, p = 0.008) exhibited a strong correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS.

Simultaneous Arterial and Venous Imaging Using 3D Quantitative Parameter Mapping

PURPOSE

To increase the number of images that can be acquired in MR examinations using quantitative parameters, we developed a method for obtaining arterial and venous images with mapping of proton density (PD), RF inhomogeneity (B1), longitudinal relaxation time (T1), apparent transverse relaxation time (T2*), and magnetic susceptibility through calculation, all with the same spatial resolution.

METHODS

The proposed method uses partially RF-spoiled gradient echo sequences to obtain 3D images of a subject with multiple scan parameters. The PD, B1, T1, T2*, and magnetic susceptibility maps are estimated using the quantification method we previously developed. Arterial images are obtained by adding images using optimized weights to emphasize the arteries. A morphology filter is used to obtain venous images from the magnetic susceptibility maps. For evaluation, images obtained from four out of five healthy volunteers were used to optimize the weights used in the arterial-image calculation, and the optimized weights were applied to the images from the fifth volunteer to obtain an arterial image. Arterial images of the five volunteers were calculated using the leave-one-out method, and the contrast between the arterial and background regions defined using the reference time-of-flight (TOF) method was evaluated using the area under the receiver operation characteristic curve (AUC). The contrast between venous and background regions defined by a reference quantitative susceptibility mapping (QSM) method was also evaluated for the venous image.

RESULTS

The AUC to discriminate blood vessels and background using the proposed method was 0.905 for the arterial image and 0.920 for the venous image.

CONCLUSION

The results indicate that the arterial images and venous images have high signal intensity at the same region as determined from the reference TOF and QSM methods, demonstrating the possibility of acquiring vasculature images with quantitative parameter mapping through calculation in an integrated manner.

Toward vendor-independent measurement of cerebral venous oxygenation: Comparison of TRUST MRI across three major MRI manufacturers and association with end-tidal CO2

Cerebral venous oxygenation (Yv ) is a valuable biomarker for a variety of brain diseases. T2 relaxation under spin tagging (TRUST) MRI is a widely used method for Yv quantification. In this work, there were two main objectives. The first was to evaluate the reproducibility of TRUST Yv measurements across MRI scanners from different vendors. The second was to examine the correlation between Yv and end-tidal CO2 (EtCO2 ) in a multisite, multivendor setting and determine the usefulness of this correlation to account for variations in Yv caused by normal variations and physiological fluctuations. Standardized TRUST pulse sequences were implemented on three scanners from major MRI vendors (GE, Siemens, Philips). These scanners were located at two research institutions. Ten healthy subjects were scanned. On each scanner, the subject underwent two scan sessions, each of which included three TRUST scans, to evaluate the intrasession and intersession reproducibility of Yv . Each scanner was also equipped with a capnograph device to record the EtCO2 of the subject during the MRI scan. We found no significant bias in Yv measurements across the three scanners (P = 0.18). The measured Yv values on the three scanners were also strongly correlated with each other (intraclass correlation coefficients > 0.85, P < 0.001). The intrasession and intersession coefficients of variation of Yv were less than 4% and showed no significant difference among the scanners. In addition, our results revealed that (1) within the same subject, Yv increased with EtCO2 at a rate of 1.24 ± 0.17%/mmHg (P < 0.0001), and (2) across different subjects, individuals with a higher EtCO2 had a higher Yv , at a rate of 0.94 ± 0.36%/mmHg (P = 0.01). These results suggest that (1) the standardized TRUST sequences had similar accuracies and reproducibilities for the quantification of Yv across the scanners, and (2) recording of EtCO2 may be a useful complement to Yv measurement to account for CO2 -related physiological fluctuations in Yv in multisite, multivendor studies.

Brain iron deposition is positively correlated with cognitive impairment in patients with chronic cerebral hypoperfusion: a MRI susceptibility mapping study

AIM

To investigate the relationship of brain iron deposition with cognitive impairment in patients with chronic cerebral hypoperfusion (CHP).

MATERIALS AND METHODS

Brain iron deposition was detected using quantitative susceptibility mapping (QSM), and cognitive function by neuropsychological tests including the Mini Mental State Examination (MMSE), Montreal Cognitive Assessment (MoCA), Activities of Daily Living (ADLs), and verbal fluency tests in a total of 40 participants, 23 with CHP and 17 age- and sex-matched healthy participants without CHP (controls).

RESULTS

The neuropsychological tests revealed that cognitive impairment (p<0.05) and susceptibility values (p<0.05) of the bilateral hippocampus (HP) and substantia nigra (SN) in CHP patients were significantly higher than those of the controls. The susceptibility values of bilateral HP and left putamen correlated closely with the scores of neuropsychological tests in the CHP patients (p<0.05, r2>0.1). The susceptibility values in the left putamen and bilateral HP were significantly higher in CHP patients with mild cognitive impairment (MCI; n=8) than those of CHP patients without MCI (n=15; p<0.05).

CONCLUSIONS

The present findings indicated that brain iron deposition in specific areas may be responsible for the cognitive impairment in CHP patients, and that QSM is a useful tool to determine brain iron, predicting cognitive impairment in CHP patients.

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.

An automatic interpretable deep learning pipeline for accurate Parkinson's disease diagnosis using quantitative susceptibility mapping and T1-weighted images

Parkinson's disease (PD) diagnosis based on magnetic resonance imaging (MRI) is still challenging clinically. Quantitative susceptibility maps (QSM) can potentially provide underlying pathophysiological information by detecting the iron distribution in deep gray matter (DGM) nuclei. We hypothesized that deep learning (DL) could be used to automatically segment all DGM nuclei and use relevant features for a better differentiation between PD and healthy controls (HC). In this study, we proposed a DL-based pipeline for automatic PD diagnosis based on QSM and T1-weighted (T1W) images. This consists of (1) a convolutional neural network model integrated with multiple attention mechanisms which simultaneously segments caudate nucleus, globus pallidus, putamen, red nucleus, and substantia nigra from QSM and T1W images, and (2) an SE-ResNeXt50 model with an anatomical attention mechanism, which uses QSM data and the segmented nuclei to distinguish PD from HC. The mean dice values for segmentation of the five DGM nuclei are all >0.83 in the internal testing cohort, suggesting that the model could segment brain nuclei accurately. The proposed PD diagnosis model achieved area under the the receiver operating characteristic curve (AUCs) of 0.901 and 0.845 on independent internal and external testing cohorts, respectively. Gradient-weighted class activation mapping (Grad-CAM) heatmaps were used to identify contributing nuclei for PD diagnosis on patient level. In conclusion, the proposed approach can potentially be used as an automatic, explainable pipeline for PD diagnosis in a clinical setting.

Toward a realistic in silico abdominal phantom for QSM

PURPOSE

QSM outside the brain has recently gained interest, particularly in the abdominal region. However, the absence of reliable ground truths makes difficult to assess reconstruction algorithms, whose quality is already compromised by additional signal contributions from fat, gases, and different kinds of motion. This work presents a realistic in silico phantom for the development, evaluation and comparison of abdominal QSM reconstruction algorithms.

METHODS

Synthetic susceptibility andR 2 * $$ {R}_2^{\ast } $$ maps were generated by segmenting and postprocessing the abdominal 3T MRI data from a healthy volunteer. Susceptibility andR 2 * $$ {R}_2^{\ast } $$ values in different tissues/organs were assigned according to literature and experimental values and were also provided with realistic textures. The signal was simulated using as input the synthetic QSM andR 2 * $$ {R}_2^{\ast } $$ maps and fat contributions. Three susceptibility scenarios and two acquisition protocols were simulated to compare different reconstruction algorithms.

RESULTS

QSM reconstructions show that the phantom allows to identify the main strengths and limitations of the acquisition approaches and reconstruction algorithms, such as in-phase acquisitions, water-fat separation methods, and QSM dipole inversion algorithms.

CONCLUSION

The phantom showed its potential as a ground truth to evaluate and compare reconstruction pipelines and algorithms. The publicly available source code, designed in a modular framework, allows users to easily modify the susceptibility,R 2 * $$ {R}_2^{\ast } $$ and TEs, and thus creates different abdominal scenarios.

Diagnosing Parkinson's disease by combining neuromelanin and iron imaging features using an automated midbrain template approach

BACKGROUND AND PURPOSE

Early diagnosis of Parkinson's disease (PD) is still a clinical challenge. Most previous studies using manual or semi-automated methods for segmenting the substantia nigra (SN) are time-consuming and, despite raters being well-trained, individual variation can be significant. In this study, we used a template-based, automatic, SN subregion segmentation pipeline to detect the neuromelanin (NM) and iron features in the SN and SN pars compacta (SNpc) derived from a single 3D magnetization transfer contrast (MTC) gradient echo (GRE) sequence in an attempt to develop a comprehensive imaging biomarker that could be used to diagnose PD.

MATERIALS AND METHODS

A total of 100 PD patients and 100 age- and sex-matched healthy controls (HCs) were imaged on a 3T scanner. NM-based SN (SN_NM) boundaries and iron-based SN (SN_QSM) boundaries and their overlap region (representing the SNpc) were delineated automatically using a template-based SN subregion segmentation approach based on quantitative susceptibility mapping (QSM) and NM images derived from the same MTC-GRE sequence. All PD and HC subjects were evaluated for the nigrosome-1 (N1) sign by two raters independently. Receiver Operating Characteristic (ROC) analyses were performed to evaluate the utility of SN_NM volume, SN_QSM volume, SNpc volume and iron content with a variety of thresholds as well as the N1 sign in diagnosing PD. Correlation analyses were performed to study the relationship between these imaging measures and the clinical scales in PD.

RESULTS

In this study, we verified the value of the fully automatic template based midbrain deep gray matter mapping approach in differentiating PD patients from HCs. The automatic segmentation of the SN in PD patients led to satisfactory DICE similarity coefficients and volume ratio (VR) values of 0.81 and 1.17 for the SN_NM, and 0.87 and 1.05 for the SN_QSM, respectively. For the HC group, the average DICE similarity coefficients and VR values were 0.85 and 0.94 for the SN_NM, and 0.87 and 0.96 for the SN_QSM, respectively. The SN_QSM volume tended to decrease with age for both the PD and HC groups but was more severe for the PD group. For diagnosing PD, the N1 sign performed reasonably well by itself (Area Under the Curve (AUC) = 0.783). However, combining the N1 sign with the other quantitative measures (SN_NM volume, SN_QSM volume, SNpc volume and iron content) resulted in an improved diagnosis of PD with an AUC as high as 0.947 (using an SN threshold of 50ppb and an NM threshold of 0.15). Finally, the SN_QSM volume showed a negative correlation with the MDS-UPDRS III (R² = 0.1, p = 0.036) and the Hoehn and Yahr scale (R² = 0.04, p = 0.013) in PD patients.

CONCLUSION

In summary, this fully automatic template based deep gray matter mapping approach performs well in the segmentation of the SN and its subregions for not only HCs but also PD patients with SN degeneration. The combination of the N1 sign with other quantitative measures (SN_NM volume, SN_QSM volume, SNpc volume and iron content) resulted in an AUC of 0.947 and provided a comprehensive set of imaging biomarkers that, potentially, could be used to diagnose PD clinically.

MRI assessment of cerebral oxygen extraction fraction in the medial temporal lobe

The medial temporal lobe (MTL) is a key area implicated in many brain diseases, such as Alzheimer's disease. As a functional biomarker, the oxygen extraction fraction (OEF) of MTL may be more sensitive than structural atrophy of MTL, especially at the early stages of diseases. However, there is a lack of non-invasive techniques to measure MTL-OEF in humans. The goal of this work is to develop an MRI technique to assess MTL-OEF in a clinically practical time without using contrast agents. The proposed method measures venous oxygenation (Yv) in the basal veins of Rosenthal (BVs), which are the major draining veins of the MTL. MTL-OEF can then be estimated as the arterio-venous difference in oxygenation. We developed an MRI sequence, dubbed arterial-suppressed accelerated T2-relaxation-under-phase-contrast (AS-aTRUPC), to quantify the blood T2 of the BVs, which was then converted to Yv through a well-established calibration model. MTL-OEF was calculated as (Ya-Yv)/Ya × 100%, where Ya was the arterial oxygenation. The feasibility of AS-aTRUPC to quantify MTL-OEF was evaluated in 16 healthy adults. The sensitivity of AS-aTRUPC in detecting OEF changes was assessed by a caffeine ingestion (200 mg) challenge. For comparison, T2-relaxation-under-spin-tagging (TRUST) MRI, which is a widely used global OEF technique, was also acquired. The dependence of MTL-OEF on age was examined by including another seven healthy elderly subjects. The results showed that in healthy adults, MTL-OEF of the left and right hemispheres were correlated (P=0.005). MTL-OEF was measured to be 23.9±3.6% (mean±standard deviation) and was significantly lower (P<0.0001) than the OEF of 33.3±2.9% measured in superior sagittal sinus (SSS). After caffeine ingestion, there was an absolute percentage increase of 9.1±4.0% in MTL-OEF. Additionally, OEF in SSS measured with AS-aTRUPC showed a strong correlation with TRUST OEF (intra-class correlation coefficient=0.94 with 95% confidence interval [0.91, 0.96]), with no significant bias (P=0.12). MTL-OEF was found to increase with age (MTL-OEF=20.997+0.100 × age; P=0.02). In conclusion, AS-aTRUPC MRI provides non-invasive assessments of MTL-OEF and may facilitate future clinical applications of MTL-OEF as a disease biomarker.

Quantifying calcium changes in the fetal spine using quantitative susceptibility mapping as extracted from STAGE imaging

OBJECTIVES

To evaluate calcium deposition in the fetal spine in vivo during the second and third trimesters using quantitative susceptibility mapping (QSM).

METHODS

Fifty-four pregnant women in their second and third trimesters underwent a 2D multi-echo STrategically Acquired Gradient Echo (STAGE) MR imaging protocol at 3T covering the fetal spine. The first echo data was used for QSM processing. A linear regression model was used to assess the correlation between magnetic susceptibility and gestational age (GA). A paired sample t-test was used to compare the consistency of QSM measurements from each sequence.

RESULTS

The magnetic susceptibility of the fetal spine decreased linearly with advancing GA, with a slope of -52.3 parts per billion (ppb)/week and a Pearson correlation coefficient (r) of 0.83 (p < 0.001). In 37 subjects for whom the STAGE local QSM data were available from both flip angles, the average magnetic susceptibility values were -1111 ± 278 ppb and -1081 ± 262 ppb for FA = 8° and FA = 40°, respectively. These means were not statistically different according to a paired sample t-test (p = 0.156).

CONCLUSIONS

QSM is a reliable technique for evaluating calcium deposition and bone mineral density of fetal vertebrae. Our results demonstrate an increase in fetal calcium levels as a function of GA. These measures might be able to provide reference values for calcium content in the fetal spine during the second and third trimesters.

KEY POINTS

• Calcium deposition and mineralization in the fetal spine, evaluated by vertebral magnetic susceptibility, increased with advancing gestational age. • Our results provide reference values for calcium content in the fetal spine during the second and third trimesters.

Test-retest analysis of cerebral oxygen extraction estimates in healthy volunteers: comparison of methods based on quantitative susceptibility mapping and dynamic susceptibility contrast magnetic resonance imaging

Background

Estimation of the oxygen extraction fraction (OEF) by quantitative susceptibility mapping (QSM) magnetic resonance imaging (MRI) is promising but requires systematic evaluation. Extraction of OEF-related information from the tissue residue function in dynamic susceptibility contrast MRI (DSC-MRI) has also been proposed. In this study, whole-brain OEF repeatability was investigated, as well as the relationships between QSM-based OEF and DSC-MRI-based parameters, i.e., mean transit time (MTT) and an oxygen extraction index, referred to as apparent OEF (AOEF).

Method

Test-retest data were obtained from 20 healthy volunteers at 3 T. QSM maps were reconstructed from 3D gradient-echo MRI phase data, using morphology-enabled dipole inversion. DSC-MRI was accomplished using gradient-echo MRI at a temporal resolution of 1.24 s.

Results

The whole-brain QSM-based OEF was (40.4±4.8) % and, in combination with a previously published cerebral blood flow (CBF) estimate, this corresponds to a cerebral metabolic rate of oxygen level of CMRO2 = 3.36 ml O2/min/100 g. The intra-class correlation coefficient [ICC(2,1)] for OEF test-retest data was 0.73. The MTT-versus-OEF and AOEF-versus-OEF relationships showed correlation coefficients of 0.61 (p = 0.004) and 0.52 (p = 0.019), respectively.

Discussion

QSM-based OEF showed a convincing absolute level and good test-retest results in terms of the ICC. Moderate to good correlations between QSM-based OEF and DSC-MRI-based parameters were observed. The present results constitute an indicator of the level of robustness that can be achieved without applying extraordinary resources in terms of MRI equipment, imaging protocol, QSM reconstruction, and OEF analysis.

Susceptibility-weighted imaging to evaluate normal and abnormal vertebrae in fetuses:a preliminary study

OBJECTIVE

To evaluate the performance of Susceptibility-weighted imaging (SWI) in visualizing normal and abnormal fetal vertebrae in vivo and in utero.

METHODS

Ninety-seven women with normal fetal vertebrae and 127 women suspected fetal vertebral anomalies on ultrasound were included in our study. SWI, TrueFISP and HASTE of the fetal spine were performed on 1.5-T MRI. The image quality and diagnostic performance between HASTE/TrueFISP and SWI were compared. Pearson correlations to correlate the L1 centrum ossification center (COC) measurements with gestational age (GA) were performed.

RESULTS

The visibility of the fetal vertebral structures on the SWI images (3.58 ± 0.69) was significantly greater than those on the HASTE (1.98 ± 0.51, P < 0.001) and TrueFISP (2.63 ± 0.52, P < 0.001). The diagnostic accuracy of SWI (89.0%) was superior to HASTE/TrueFISP (48.0%) (P < 0.001) and the area under the curve (AUC) for SWI was 0.909 (P < 0.001). The height, transverse, sagittal diameter and area of L1 COC were linearly correlated with GA (all P < 0.001).

CONCLUSION

SWI proved to be a reliable method for depicting fetal vertebral structure and growth, which can significantly improve the diagnostic performance of vertebral anomalies in fetuses. This article is protected by copyright. All rights reserved.

Ex vivo susceptibility-weighted imaging anatomy of canine brain–comparison of imaging and histological sections

Now that access of large domestic mammals to high-field MRI becomes more common, techniques initially implemented for human patients can be used for the structural and functional study of the brain of these animals. Among them, susceptibility-weighted imaging (SWI) is a recent technique obtained from gradient echo (GE) imaging that allow for an excellent anatomical tissue contrast and a non-invasive assessment of brain iron content. The goal of this study was to design an optimal GE SWI imaging protocol to be used in dogs undergoing an MRI examination of the brain in a 3-Tesla scanner. This imaging protocol was applied to ex vivo brains from four dogs. The imaging protocol was validated by visual inspection of the SWI images that provided a high anatomical detail, as demonstrated by their comparison with corresponding microscopic sections. As resolvable brain structures were labeled, this study is the first to provide an anatomic description of SWI images of the canine brain. Once validated in living animals, this GE SWI imaging protocol could be easily included in routine neuroimaging protocols to improve the diagnosis of various intracranial diseases of dogs, or be used in future comparative studies aiming at evaluating brain iron content in animals.

Quantitative susceptibility atlas construction in Montreal Neurological Institute space: towards histological-consistent iron-rich deep brain nucleus subregion identification

Iron-rich deep brain nuclei (DBN) of the human brain are involved in various motoric, emotional and cognitive brain functions. The abnormal iron alterations in the DBN are closely associated with multiple neurological and psychiatric diseases. Quantitative susceptibility mapping (QSM) provides the spatial distribution of the magnetic susceptibility of human brain tissues. Compared to traditional structural imaging, QSM provides superiority for imaging the iron-rich DBN owing to the susceptibility difference existing between brain tissues. In this study, we constructed a Montreal Neurological Institute (MNI) space unbiased QSM human brain atlas via group-wise registration from 100 healthy subjects aged 19-29 years. The atlas construction process was guided by hybrid images that were fused from multi-modal magnetic resonance images (MRI). We named it as Multi-modal-fused magnetic Susceptibility (MuSus-100) atlas. The high-quality susceptibility atlas provides extraordinary image contrast between iron-rich DBN with their surroundings. Parcellation maps of DBN and their subregions that are highly related to neurological and psychiatric pathology were then manually labeled based on the atlas set with the assistance of an image border-enhancement process. Especially, the bilateral thalamus was delineated into 64 detailed subregions referring to the Schaltenbrand-Wahren stereotactic atlas. To our best knowledge, the histological-consistent thalamic nucleus parcellation map is well defined for the first time in the MNI space. Compared with existing atlases that emphasizing DBN parcellation, the newly proposed atlas outperforms on the task of atlas-guided individual brain image DBN segmentation both in accuracy and robustness. Moreover, we applied the proposed DBN parcellation map to conduct detailed identification of the pathology-related iron content alterations in subcortical nuclei for Parkinson's Disease (PD) patients. We envision that the MuSus-100 atlas can play a crucial role in improving the accuracy of DBN segmentation for the research of neurological and psychiatric disease progress and also be helpful for target planning in deep brain stimulation surgery.

Cerebral Venous Oxygen Saturation in Hypoperfusion Regions May Become a New Imaging Indicator to Predict the Clinical Outcome of Stroke

To automatically and quantitatively evaluate the venous oxygen saturation (SvO2) in cerebral ischemic tissues and explore its value in predicting prognosis. A retrospective study was conducted on 48 AIS patients hospitalized in our hospital from 2015−2018. Based on quantitative susceptibility mapping and perfusion-weighted imaging, this paper measured the cerebral SvO2 in hypoperfusion tissues and its change after intraarterial rt-PA treatment. The cerebral SvO2 in different hypoperfusion regions between the favorable and unfavorable clinical outcome groups was analyzed using an independent t-test. Relationships between cerebral SvO2 and clinical scores were determined using the Pearson correlation coefficient. The receiver operating characteristic process was conducted to evaluate the accuracy of cerebral SvO2 in predicting unfavorable clinical outcomes. Cerebral SvO2 in hypoperfusion (Tmax > 4 and 6 s) was significantly different between the two groups at follow-up (p < 0.05). Cerebral SvO2 and its changes before and after treatment were negatively correlated with clinical scores. The positive predictive value, negative predictive value, accuracy, and area under the curve of the cerebral SvO2 were higher than those predicted by the ischemic core. Therefore, the cerebral SvO2 of hypoperfusion regions was a stronger imaging predictor of unfavorable clinical outcomes after stroke.

Quantitative Distribution of Cerebral Venous Oxygen Saturation and Its Prognostic Value in Patients with Acute Ischemic Stroke

This study investigated the quantitative distribution of cerebral venous oxygen saturation (SvO2) based on quantitative sensitivity mapping (QSM) and determined its prognostic value in patients with acute ischemic stroke (AIS). A retrospective study was conducted on 39 hospitalized patients. Reconstructed QSM was used to calculate the cerebral SvO2 of each region of interest (ROI) in the ischemic hemisphere. The intraclass correlation coefficient (ICC) and Bland−Altman analysis were conducted to define the best resolution of the distribution map. The correlation between the cerebral SvO2 in hypoxic regions (SvO2ROI < 0.7) and clinical scores was obtained by Spearman and power analysis. The associations between cerebral SvO2 and unfavorable prognosis were analyzed using multivariate logistic regression. Excellent agreement was found between the cerebral SvO2 in hypoxic regions with a resolution of 7.18 × 7.18 × 1.6 mm3 and asymmetrically prominent cortical veins regions (ICC: 0.879 (admission), ICC: 0.906 (discharge)). The cerebral SvO2 was significantly negative with clinical scores (all |r| > 0.3). The cerebral SvO2 and its changes at discharge were significantly associated with an unfavorable prognosis (OR: 0.812 and 0.866). Therefore, the cerebral SvO2 in hypoxic regions measured by the quantitative distribution map can be used as an indicator for evaluating the early prognosis of AIS.

Quantitative Susceptibility Mapping of Brain Iron Deposition in Patients With Recurrent Depression

OBJECTIVE

Recurrence is the most significant feature of depression and the relationship between iron and recurrent depression is still lack of direct evidence in vivo.

METHODS

Twenty-one patients with depression and twenty control subjects were included. Gradient-recalled echo, T1 and T2 images were acquired using a 3.0T MRI system. After quantitative susceptibility mapping were reconstructed and standardized, a whole-brain and the regions of interest were respectively analyzed.

RESULTS

Significant increases in susceptibility were found in multiple recurrent depression patients, which involved several brain regions (frontal lobes, temporal lobe structures, occipital lobes hippocampal regions, putamen, thalamus, cingulum, and cerebellum). Interestingly, no susceptibility changes after treatment compared to pre-treatment (all p>0.05) and no significant correlation between susceptibility and Hamilton Depression Rating Scale were found. Besides, it was close to significance that those with a higher relapse frequency or a longer mean duration of single episode had a higher susceptibility in the putamen, thalamus, and hippocampus. Further studies showed susceptibility across the putamen (ρ2=0.27, p<0.001), thalamus (ρ2=0.21, p<0.001), and hippocampus (ρ2=0.19, p<0.001) were strongly correlated with total course of disease onset.

CONCLUSION

Brain iron deposition is related to the total course of disease onset, but not the severity of depression, which suggest that brain iron deposition may be a sign of brain damage in multiple recurrent depression.

CAU-Net: A Deep Learning Method for Deep Gray Matter Nuclei Segmentation

The abnormal iron deposition of the deep gray matter nuclei is related to many neurological diseases. With the quantitative susceptibility mapping (QSM) technique, it is possible to quantitatively measure the brain iron content in vivo. To assess the magnetic susceptibility of the deep gray matter nuclei in the QSM, it is mandatory to segment the nuclei of interest first, and many automatic methods have been proposed in the literature. This study proposed a contrast attention U-Net for nuclei segmentation and evaluated its performance on two datasets acquired using different sequences with different parameters from different MRI devices. Experimental results revealed that our proposed method was superior on both datasets over other commonly adopted network structures. The impacts of training and inference strategies were also discussed, which showed that adopting test time augmentation during the inference stage can impose an obvious improvement. At the training stage, our results indicated that sufficient data augmentation, deep supervision, and nonuniform patch sampling contributed significantly to improving the segmentation accuracy, which indicated that appropriate choices of training and inference strategies were at least as important as designing more advanced network structures.

Cerebral oxygen extraction fraction MRI: Techniques and applications

The human brain constitutes 2% of the body's total mass but uses 20% of the oxygen. The rate of the brain's oxygen utilization can be derived from a knowledge of cerebral blood flow and the oxygen extraction fraction (OEF). Therefore, OEF is a key physiological parameter of the brain's function and metabolism. OEF has been suggested to be a useful biomarker in a number of brain diseases. With recent advances in MRI techniques, several MRI-based methods have been developed to measure OEF in the human brain. These MRI OEF techniques are based on the T₂ of blood, the blood signal phase, the magnetic susceptibility of blood-containing voxels, the effect of deoxyhemoglobin on signal behavior in extravascular tissue, and the calibration of the BOLD signal using gas inhalation. Compared to ¹⁵ O PET, which is considered the "gold standard" for OEF measurement, MRI-based techniques are non-invasive, radiation-free, and are more widely available. This article provides a review of these emerging MRI-based OEF techniques. We first briefly introduce the role of OEF in brain oxygen homeostasis. We then review the methodological aspects of different categories of MRI OEF techniques, including their signal mechanisms, acquisition methods, and data analyses. The strengths and limitations of the techniques are discussed. Finally, we review key applications of these techniques in physiological and pathological conditions.

Quantitative susceptibility mapping as a measure of cerebral oxygenation in neonatal piglets

Prominence of cerebral veins using susceptibility weighted magnetic resonance imaging (SWI) has been used as a qualitative indicator of cerebral venous oxygenation (CvO₂). Quantitative susceptibility mapping (QSM) adds more precision to the assessment of CvO₂, but has not been applied to neonatal hypoxic ischemic injury (HII). We proposed to study QSM measures of venous susceptibility and their correlation with direct measures of brain oxygenation and cerebral blood flow (CBF) in the neonatal piglet. The association of QSM intravascular cerebral venous susceptibility, with brain tissue O₂ tension, CBF, cortical tissue oxyhemoglobin saturation, and the partial pressure of oxygen in arterial blood measurement during various oxygenation states was determined by linear regression. Compared to normoxia, venous susceptibility in the straight sinus increased 56.8 ± 25.4% during hypoxia, while decreasing during hyperoxia (23.5 ± 32.9%) and hypercapnia (23.3 ± 73.1%), which was highly correlated to all other measures of oxygenation (p < 0.0001) but did not correlate to CBF (p = 0.82). These findings demonstrate a strong relationship between venous susceptibility and brain tissue O₂ tension. Our results suggest that QSM-derived venous susceptibility is sensitive to cerebral oxygenation status across various oxygenation states.

Evaluation of highly accelerated wave controlled aliasing in parallel imaging (Wave-CAIPI) susceptibility-weighted imaging in the non-sedated pediatric setting: a pilot study

BACKGROUND

Susceptibility-weighted imaging (SWI) is highly sensitive for intracranial hemorrhagic and mineralized lesions but is associated with long scan times. Wave controlled aliasing in parallel imaging (Wave-CAIPI) enables greater acceleration factors and might facilitate broader application of SWI, especially in motion-prone populations.

OBJECTIVE

To compare highly accelerated Wave-CAIPI SWI to standard SWI in the non-sedated pediatric outpatient setting, with respect to the following variables: estimated scan time, image noise, artifacts, visualization of normal anatomy and visualization of pathology.

MATERIALS AND METHODS

Twenty-eight children (11 girls, 17 boys; mean age ± standard deviation [SD] = 128.3±62 months) underwent 3-tesla (T) brain MRI, including standard three-dimensional (3-D) SWI sequence followed by a highly accelerated Wave-CAIPI SWI sequence for each subject. We rated all studies using a predefined 5-point scale and used the Wilcoxon signed rank test to assess the difference for each variable between sequences.

RESULTS

Wave-CAIPI SWI provided a 78% and 67% reduction in estimated scan time using the 32- and 20-channel coils, respectively, corresponding to estimated scan time reductions of 3.5 min and 3 min, respectively. All 28 children were imaged without anesthesia. Inter-reader agreement ranged from fair to substantial (k=0.67 for evaluation of pathology, 0.55 for anatomical contrast, 0.3 for central noise, and 0.71 for artifacts). Image noise was rated higher in the central brain with wave SWI (P<0.01), but not in the peripheral brain. There was no significant difference in the visualization of normal anatomical structures and visualization of pathology between the standard and wave SWI sequences (P=0.77 and P=0.79, respectively).

CONCLUSION

Highly accelerated Wave-CAIPI SWI of the brain can provide similar image quality to standard SWI, with estimated scan time reduction of 3-3.5 min depending on the radiofrequency coil used, with fewer motion artifacts, at a cost of mild but perceptibly increased noise in the central brain.

Motor Band Sign in Motor Neuron Disease: A Marker for Upper Motor Neuron Involvement

BACKGROUND AND OBJECTIVE

The prevalence and role of the motor band sign (MBS) remain unclear in motor neuron disease. We report the frequency of MBS in amyotrophic lateral sclerosis (ALS) and primary lateral sclerosis (PLS), its correlation with clinical upper motor neuron (UMN) signs, and prognostic value in ALS.

METHODS

We conducted a retrospective study of ALS, PLS, and controls with retrievable MRI between 2010 and 2018. We compared the frequencies of MBS across the three groups, and studied correlation between susceptibility-weighted MRI measurements in primary motor cortices and contralateral UMN features. Clinical outcomes were compared between ALS with and without MBS.

RESULTS

Thirteen ALS, 5 PLS, and 10 controls were included (median age 60 years, IQR 54-66 years; 14/28 males). MBS was present in 9/13 (69.2%, 95% CI 38.9-89.6%) and 4/5 (80.0%, 95% CI 29.9-99.0%) of ALS and PLS, respectively, and none in controls. 2/13 (15.4%, 95% CI 2.7-46.3%) ALS and 3/5 (60.0%, 95% CI 17.0-92.7%) PLS had MBS in the absence of corticospinal T2/FLAIR hyperintensity sign. Susceptibility measurements in left motor cortices had a significantly positive correlation with contralateral UMN signs in ALS (τ_b = 0.628, p = 0.03). Similar but nonsignificant trends was observed for right motor cortices in ALS (τ_b = 0.516, p = 0.07). There were no significant differences in mRS at last follow-up, mortality, or time from symptom onset to last follow-up between ALS patients with and without MBS.

CONCLUSIONS

We provide limited evidence that MBS and susceptibility quantification measurements in motor cortices may serve as surrogate markers of UMN involvement in motor neuron disease.

Hybrid data fidelity term approach for quantitative susceptibility mapping

Purpose

Methods

Results

Conclusions

Does the Brush-Sign Reflect Collateral Status and DWI-ASPECTS in Large Vessel Occlusion?

Introduction

The relevance of the brush-sign remained poorly documented in large vessel occlusion (LVO). We aimed to assess the relationship between the brush-sign and collateral status and its potential impact on baseline diffusion-weighted imaging–Alberta Stroke Program Early Computed Tomography Score (DWI-ASPECTS) in acute ischemic stroke (AIS) patients eligible to mechanical thrombectomy (MT).

Methods

Consecutive patients admitted in the Lyon Stroke Center with anterior circulation AIS due to intracranial internal carotid artery (ICA) and/or M1 or M2 segment of the middle cerebral artery (MCA) occlusion eligible for MT were included. The brush-sign was assessed on T2-gradient-echo MRI. Collateral status was assessed on digital subtraction angiography according to the American Society of Interventional and Therapeutic Neuroradiology/Society of Interventional Radiology (ASITN/SIR) score.

Results

In this study, 504 patients were included, among which 171 (33.9%) patients had a brush-sign. Patients with a brush-sign more frequently had a poor collateral status [72 (42.1%) vs. 103 (30.9%); p = 0.017]. In univariable analysis, a DWI-ASPECTS < 7 was associated with a brush sign. Following multivariable analysis, the brush-sign no longer affected DWI-ASPECTS < 7 while the latter remained associated with younger age [odds ratio (OR) 0.97, 95% CI.96–0.99], male sex (OR 1.79, 95% CI 1.08–2.99), a higher National Institutes of Health Stroke Scale (NIHSS) score (OR 1.16, 95% CI 1.1–1.21), a poor collateral status (OR 9.35, 95% CI 5.59-16.02), MCA segment (OR 2.54, 95% CI 1.25–5.38), and intracranial ICA (OR 3.01, 95% CI 1.16–8) occlusion.

Conclusions and Relevance

The brush-sign may be a marker of poor collateral status but did not independently predict a lower DWI-ASPECTS.

Clinical Trial Registration

ClinicalTrials.gov, identifier: NCT04620642.

Three-dimensional simultaneous brain mapping of T1, T2, T2∗ and magnetic susceptibility with MR Multitasking

PURPOSE

To develop a new technique that enables simultaneous quantification of whole-brain T1 , T2 , T 2 ∗ , as well as susceptibility and synthesis of six contrast-weighted images in a single 9.1-minute scan.

METHODS

The technique uses hybrid T2 -prepared inversion-recovery pulse modules and multi-echo gradient-echo readouts to collect k-space data with various T1, T2, and T 2 ∗ weightings. The underlying image is represented as a six-dimensional low-rank tensor consisting of three spatial dimensions and three temporal dimensions corresponding to T1 recovery, T2 decay, and multi-echo behaviors, respectively. Multiparametric maps were fitted from reconstructed image series. The proposed method was validated on phantoms and healthy volunteers, by comparing quantitative measurements against corresponding reference methods. The feasibility of generating six contrast-weighted images was also examined.

RESULTS

High quality, co-registered T1 , T2 , and T 2 ∗ susceptibility maps were generated that closely resembled the reference maps. Phantom measurements showed substantial consistency (R2 > 0.98) with the reference measurements. Despite the significant differences of T1 (p < .001), T2 (p = .002), and T 2 ∗ (p = 0.008) between our method and the references for in vivo studies, excellent agreement was achieved with all intraclass correlation coefficients greater than 0.75. No significant difference was found for susceptibility (p = .900). The framework is also capable of synthesizing six contrast-weighted images.

CONCLUSION

The MR Multitasking-based 3D brain mapping of T1 , T2 , T 2 ∗ , and susceptibility agrees well with the reference and is a promising technique for multicontrast and quantitative imaging.

Penumbra Detection With Oxygen Extraction Fraction Using Magnetic Susceptibility in Patients With Acute Ischemic Stroke

Background

The oxygen extraction fraction (OEF) has been applied to identify ischemic penumbral tissue, but is difficult to use in an urgent care setting. This study aimed to investigate whether an OEF map generated via magnetic resonance quantitative susceptibility mapping (QSM) could help identify the ischemic penumbra in patients with acute ischemic stroke.

Materials and Methods

This prospective imaging study included 21 patients with large anterior circulation vessel occlusion who were admitted <24 h after stroke onset and 21 age-matched healthy controls. We identified the ischemic penumbra as the region with a Tmax of >6 s during dynamic susceptibility contrast-magnetic resonance imaging (DSC-MRI) and calculated the perfusion-core mismatch ratio between the ischemic penumbra and infarct core volumes. The OEF values were measured based on magnetic susceptibility differences between the venous structures and brain tissues using rapid QSM acquisition. Volumes with increased OEF values were compared to the ischemic penumbra volumes using an anatomical template.

Results

Eleven patients had a perfusion-core mismatch ratio of ≥1.8, and reperfusion therapy was recommended. In these patients, the volumes with increased OEF values of >51.5%, which was defined using the anterior circulation territory OEF values from the 21 healthy controls, were positively correlated with the ischemic penumbra volumes (r = 0.636, 95% CI: 0.059 to 0.895, P = 0.035) and inversely correlated with the 30-day change in the National Institutes of Health Stroke Scale scores (r = −0.624, 95% CI: −0.891 to −0.039, P = 0.041).

Conclusion

Tissue volumes with increased OEF values could predict ischemic penumbra volumes based on DSC-MRI, highlighting the potential of the QSM-derived OEF map as a penumbra biomarker to guide treatment selection in patients with acute ischemic stroke.

Comparison of MRI methods for measuring whole‐brain oxygen extraction fraction under different geometric conditions at 7T

Background and Purpose

Methods

Results

Conclusion

Vascular Mapping of the Human Hippocampus Using Ferumoxytol-Enhanced MRI

The hippocampus is a small but complex grey matter structure that plays an important role in spatial and episodic memory and can be affected by a wide range of pathologies including vascular abnormalities. In this work, we introduce the use of Ferumoxytol, an ultra-small superparamagnetic iron oxide (USPIO) agent, to induce susceptibility in the arteries (as well as increase the susceptibility in the veins) to map the hippocampal micro-vasculature and to evaluate the quantitative change in tissue fractional vascular density (FVD), in each of its subfields. A total of 39 healthy subjects (aged 35.4 ± 14.2 years, from 18 to 81 years old) were scanned with a high-resolution (0.22×0.44×1 mm³) dual-echo SWI sequence acquired at four time points during a gradual increase in Ferumoxytol dose (final dose = 4 mg/kg). The volumes of each subfield were obtained automatically from the pre-contrast T₁ -weighted data. The dynamically acquired SWI data were co-registered and adaptively combined to reduce the blooming artifacts from large vessels, preserving the contrast from smaller vessels. The resultant SWI data were used to segment the hippocampal vasculature and to measure the FVD ((volume occupied by vessels)/(total volume)) for each subfield. The hippocampal fissure, along with the fimbria, granular cell layer of the dentate gyrus and cornu ammonis layers (except for CA1), showed higher micro-vascular FVD than the other parts of hippocampus. The CA1 region exhibited a significant correlation with age (R = −0.37, p < 0.05). demonstrating an overall loss of hippocampal vascularity in the normal aging process. Moreover, the vascular density reduction was more prominent than the age correlation with the volume reduction (R = −0.1, p > 0.05) of the CA1 subfield, which would suggest that vascular degeneration may precede tissue atrophy.

Automatic detection of neuromelanin and iron in the midbrain nuclei using a magnetic resonance imaging ‐based brain template

Parkinson disease (PD) is a chronic progressive neurodegenerative disorder characterized pathologically by early loss of neuromelanin (NM) in the substantia nigra pars compacta (SNpc) and increased iron deposition in the substantia nigra (SN). Degeneration of the SN presents as a 50 to 70% loss of pigmented neurons in the ventral lateral tier of the SNpc at the onset of symptoms. Also, using magnetic resonance imaging (MRI), iron deposition and volume changes of the red nucleus (RN), and subthalamic nucleus (STN) have been reported to be associated with disease status and rate of progression. Further, the STN serves as an important target for deep brain stimulation treatment in advanced PD patients. Therefore, an accurate in‐vivo delineation of the SN, its subregions and other midbrain structures such as the RN and STN could be useful to better study iron and NM changes in PD. Our goal was to use an MRI template to create an automatic midbrain deep gray matter nuclei segmentation approach based on iron and NM contrast derived from a single, multiecho magnetization transfer contrast gradient echo (MTC‐GRE) imaging sequence. The short echo TE = 7.5 ms data from a 3D MTC‐GRE sequence was used to find the NM‐rich region, while the second echo TE = 15 ms was used to calculate the quantitative susceptibility map for 87 healthy subjects (mean age ± SD: 63.4 ± 6.2 years old, range: 45–81 years). From these data, we created both NM and iron templates and calculated the boundaries of each midbrain nucleus in template space, mapped these boundaries back to the original space and then fine‐tuned the boundaries in the original space using a dynamic programming algorithm to match the details of each individual's NM and iron features. A dual mapping approach was used to improve the performance of the morphological mapping of the midbrain of any given individual to the template space. A threshold approach was used in the NM‐rich region and susceptibility maps to optimize the DICE similarity coefficients and the volume ratios. The results for the NM of the SN as well as the iron containing SN, STN, and RN all indicate a strong agreement with manually drawn structures. The DICE similarity coefficients and volume ratios for these structures were 0.85, 0.87, 0.75, and 0.92 and 0.93, 0.95, 0.89, 1.05, respectively, before applying any threshold on the data. Using this fully automatic template‐based deep gray matter mapping approach, it is possible to accurately measure the tissue properties such as volumes, iron content, and NM content of the midbrain nuclei.

A novel gradient echo data based vein segmentation algorithm and its application for the detection of regional cerebral differences in venous susceptibility

Accurate segmentation of cerebral venous vasculature from gradient echo data is of central importance in several areas of neuroimaging such as for the susceptibility-based assessment of brain oxygenation or planning of electrode placement in deep brain stimulation. In this study, a vein segmentation algorithm for single- and multi-echo gradient echo data is proposed. First, susceptibility maps, true susceptibility-weighted images, and, in the multi-echo case, R₂^* maps were generated from the gradient echo data. These maps were filtered with an inverted Hamming filter to suppress background contrast as well as artifacts from field inhomogeneities at the brain boundaries. A shearlet-based scale-wise representation was generated to calculate a vesselness function and to generate segmentations based on local thresholding. The accuracy of the proposed algorithm was evaluated for different echo times and image resolutions using a manually generated reference segmentation and two vein segmentation algorithms (Frangi vesselness-based, recursive vesselness filter) as a reference with the Dice and Cohen's coefficients as well as the modified Hausdorff distance. The Frangi-based and recursive vesselness filter methods were significantly outperformed with regard to all error metrics. Applying the algorithm, susceptibility differences likely related to differences in blood oxygenation between superficial and deep venous territories could be demonstrated.

Edge prior guided dictionary learning for quantitative susceptibility mapping reconstruction

BACKGROUND

Compared with conventional magnetic resonance imaging methods, the quantitative magnetic susceptibility mapping (QSM) technique can quantitatively measure the magnetic susceptibility distribution of tissues, which has an important clinical application value in the investigations of brain micro-bleeds, Parkinson's, and liver iron deposition, etc. However, the quantitative susceptibility mapping algorithm is an ill-posed inverse problem due to the near-zero value in the dipole kernel, and high-quality QSM reconstruction with effective streaking artifact suppression remains a challenge. In recent years, the performance of sparse representation has been well validated in improving magnetic resonance image (MRI) reconstruction.

METHODS

In this study, by incorporating feature learning into sparse representation, we propose an edge prior guided dictionary learning-based reconstruction method for the dipole inversion in quantitative susceptibility mapping reconstruction. The structure feature dictionary relies on magnitude images for susceptibility maps have similar structures with magnitude images, and this structure feature dictionary and edge prior information are used in the dipole inversion step.

RESULTS

The performance of the proposed algorithm is assessed through in vivo human brain clinical data, leading to high-quality susceptibility maps with improved streaking artifact suppression, structural recovery, and quantitative metrics.

CONCLUSIONS

The proposed edge prior guided dictionary learning method for dipole inversion in QSM achieves improved performance in streaking artifacts suppression, structural recovery and deep gray matter reconstruction.

Inverse Fourier transformation of combined first order derivative and intensity-curvature functional of magnetic resonance angiography of the human brain

BACKGROUND AND OBJECTIVE

This paper reports a novel image processing technique based on inverse Fourier transformation and its validation procedure.

METHODS

Magnetic Resonance Angiography (MRA) data of the human brain is fitted on a pixel-by-pixel basis with bivariate linear model polynomial function. Polynomial fitting allows the formulation of two measures: the first order derivative (FOD), which is an edge finder, and the intensity-curvature functional (ICF), which is a high pass filter. The calculation of FOD and ICF uses knowledge provided by existing research and is performed through resampling. ICF and FOD are direct Fourier transformed, and their k-space is combined through a nonlinear convolution of terms. The resulting k-space is inverse Fourier transformed so to obtain a novel image called Fourier Convolution Image (FCI).

RESULTS

FCI possesses the characteristics of an edge finder (FOD) and a high pass filter (ICF).

CONCLUSIONS

FC images yield the following properties versus MRA: 1. Change of the contrast; 2. Increased sharpness in the proximity of human brain vessels; 3. Increased visualization of vessel connectivity. The implication of this study is to provide FCI as another viable option for MRA evaluation.

Comparison of ultrafast wave-controlled aliasing in parallel imaging (CAIPI) magnetization-prepared rapid acquisition gradient echo (MP-RAGE) and standard MP-RAGE in non-sedated children: initial clinical experience

BACKGROUND

Fast magnetic resonance imaging (MRI) sequences are advantageous in pediatric imaging as they can lessen child discomfort, decrease motion artifact and improve scanner availability.

OBJECTIVE

To evaluate the feasibility of an ultrafast wave-CAIPI (controlled aliasing in parallel imaging) MP-RAGE (magnetization-prepared rapid gradient echo) sequence for brain imaging of awake pediatric patients.

MATERIALS AND METHODS

Each MRI included a standard MP-RAGE sequence and an ultrafast wave-MP-RAGE sequence. Two neuroradiologists evaluated both sequences in terms of artifacts, noise, anatomical contrast and pathological contrast. A predefined 5-point scale was used by two independent pediatric neuroradiologists. A Wilcoxon signed-rank test was used to evaluate the difference between sequences for each variable.

RESULTS

Twenty-four patients (14 males; mean age: 11.5±4.5 years, range: 1 month to 17.8 years) were included. Wave-CAIPI MP-RAGE provided a 77% reduction in scan time using a 32-channel coil and a 70% reduction using a 20-channel coil. Visualization of the pathology, artifacts and pathological enhancement (including parenchymal, leptomeningeal and dural enhancement) was not significantly different between standard MP-RAGE and wave-CAIPI MP-RAGE (all P>0.05). For central (P<0.001) and peripheral (P<0.001) noise, and the evaluation of the anatomical structures (P<0.001), the observers favored standard MP-RAGE over wave-CAIPI MP-RAGE.

CONCLUSION

Ultrafast brain imaging with wave-CAIPI MP-RAGE is feasible in awake pediatric patients, providing a substantial reduction in scan time at a cost of subjectively increased image noise.

The Frequency and Associated Factors of Asymmetrical Prominent Veins: A Predictor of Unfavorable Outcomes in Patients with Acute Ischemic Stroke

Objectives

The present study is aimed at investigating the frequency and associated factors of asymmetrical prominent veins (APV) in patients with acute ischemic stroke (AIS).

Methods

Consecutive patients with AIS admitted to the Comprehensive Stroke Center of Shanghai Fourth People's Hospital between January 2013 and December 2017 were enrolled. MRI including diffusion-weighted imaging (DWI), perfusion-weighted imaging (PWI), and susceptibility-weighted imaging (SWI) was performed within 12 hours of symptom onset. The volume of asymmetrical prominent veins (APV) was evaluated using the Signal Processing In nuclear magnetic resonance software (SPIN, Detroit, Michigan, USA). Multivariate analysis was used to assess relationships between APV findings and medical history, clinical variables as well as cardio-metabolic indices.

Results

Seventy-six patients met the inclusion criteria. The frequency of APV ≥ 10 mL was 46.05% (35/76). Multivariate analyses showed that proximal artery stenosis or occlusion (≥50%) (P < 0.001, adjusted odds ratio (OR) = 660.0, 95%CI = 57.28-7604.88) and history of atrial fibrillation (P < 0.001, adjusted OR = 10.48, 95%CI = 1.78-61.68) were independent factors associated with high APV (≥10 mL).

Conclusion

Our findings suggest that the frequency of APV ≥ 10 mL is high in patients with AIS within 12 hours of symptom onset. History of atrial fibrillation and severe proximal artery stenosis or occlusion are strong predictors of high APV as calculated by SPIN on the SWI map.

Utility of quantitative susceptibility mapping and diffusion kurtosis imaging in the diagnosis of early Parkinson's disease

•

Putamen susceptibility value was higher in PD than controls one year into diagnosis.

•

Putamen susceptibility value was associated with clinical motor scores in early PD.

•

Mean diffusivity revealed greater cellular loss in the lateral substantial nigra.

•

Putamen and caudate microstructural degradation were driven by radial diffusivity.

•

A composite putamen-caudate DKI-QSM marker classified early PD from controls.

Objective

To investigate the utility of quantitative susceptibility mapping (QSM) and diffusion kurtosis imaging (DKI) as complementary tools in characterizing pathological changes in the deep grey nuclei in early Parkinson’s disease (PD) and their clinical correlates to aid in diagnosis of PD.

Method

Patients with a diagnosis of PD made within a year and age-matched healthy controls were recruited. All participants underwent clinical evaluation using the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS III) and Hoehn & Yahr stage (H&Y), and brain 3 T MRI including QSM and DKI. Regions-of-interest (ROIs) in the caudate nucleus, putamen, globus pallidus, and medial and lateral substantia nigra (SN) were manually drawn to compare the mean susceptibility (representing iron deposition) and DKI indices (representing restricted water diffusion) between PD patients and healthy controls and in correlation with MDS-UPDRS III and H&Y, focusing on susceptibility value, mean diffusivity (MD) and mean kurtosis (MK).

Results

There were forty-seven PD patients (aged 68.7 years, 51% male, disease duration 0.78 years) and 16 healthy controls (aged 67.4 years, 63% male). Susceptibility value was increased in PD in all ROIs except the caudate, and was significantly different after multiple comparison correction in the putamen (PD: 64.75 ppb, HC: 44.61 ppb, p = 0.004). MD was significantly higher in PD in the lateral SN, putamen and caudate, the regions with the lowest susceptibility value. In PD patients, we found significant association between the MDS-UPDRS III score and susceptibility value in the putamen after correcting for age and sex (β = 0.21, p = 0.003). A composite DKI-QSM diagnostic marker based on these findings successfully differentiated the groups (p < 0.0001) and had “good” classification performance (AUC = 0.88).

Conclusions

QSM and DKI are complementary tools allowing a better understanding of the complex contribution of iron deposition and microstructural changes in the pathophysiology of PD.

Quantitative susceptibility mapping in ischemic stroke patients after successful recanalization

Quantitative susceptibility mapping (QSM) is a novel processing method for gradient-echo magnetic resonance imaging (MRI). Higher magnetic susceptibility in cortical veins have been observed on susceptibility maps in the ischemic hemisphere of stroke patients, indicating an increased oxygen extraction fraction (OEF). Our goal was to investigate susceptibility in veins of stroke patients after successful recanalization in order to analyze the value of QSM in predicting tissue prognosis and clinical outcome. We analyzed MR images of 23 patients with stroke due to unilateral middle cerebral artery (MCA)-M1/M2 occlusion acquired 24–72 h after successful thrombectomy. The susceptibilities of veins were obtained from QSM and compared between the stroke territory, the ipsilateral non-ischemic MCA territory and the contralateral MCA territory. As outcome variables, early infarct size and functional disability (modified Rankin Scale, mRS) after 3–5 months was used. The median susceptibility value of cortical veins in the ischemic core was 41% lower compared to the ipsilateral non-ischemic MCA territory and 38% lower than on the contralateral MCA territory. Strikingly, in none of the patients prominent vessels with high susceptibility signal were found after recanalization. Venous susceptibility values within the infarct did not correlate with infarct volume or functional disability after 3–5 months. Low venous susceptibility within the infarct core after successful recanalization of the occluded vessel likely indicates poor oxygen extraction arising from tissue damage. We did not identify peri-infarct tissue with increased susceptibility values as potential surrogate of former penumbral areas. We found no correlation of QSM parameters with infarct size or outcome.

Investigation of the magnetic susceptibility properties of fresh and fixed mouse heart, liver, skeletal muscle and brain tissue

PURPOSE

Several magnetic resonance imaging (MRI) techniques exploit the difference in magnetic susceptibilities between tissues, but systematic measurements of tissue susceptibility are lacking. Furthermore, there is the question as to whether chemical fixation that is used for ex vivo MRI studies, affects the magnetic properties of the tissue. Here, we determined the magnetic susceptibility and water content of fresh and chemically fixed mouse tissue.

METHODS

Mass susceptibility of brain, heart, liver and skeletal muscle samples were determined on a vibrating sample magnetometer at room temperature. Measurements at 50, 125, 200 and 295 K were performed to assess the temperature dependence of susceptibility. Moreover, we measured water content of fresh and fixed samples.

RESULTS

All samples show mass susceptibilities between -0.068 and -1.929 × 10^-8 m³/kg, compared to -9.338 × 10^-9 m³/kg of double distilled water. Heart tissue has a more diamagnetic susceptibility than the other tissues. Compared to fresh tissue, fixed tissue has a less diamagnetic susceptibility. Fixed tissue was not different in water content to fresh tissue and showed no consistent dependence of susceptibility with temperature, whereas fresh tissue shows a decrease to at least 125 K, indicative of a paramagnetic component.

CONCLUSIONS

Biological tissues are diamagnetic in comparison to water, where the heart is more diamagnetic than the other tissues, with paramagnetic contributions. Fixation rendered tissue less diamagnetic compared to fresh tissue. Our measurements revealed differences in tissue susceptibility between VSM and QSM, inviting more research to compare susceptibility-based MRI methods with physical measurements of tissue susceptibility.

Assessing brain iron and volume of subcortical nuclei in idiopathic rapid eye movement sleep behavior disorder

STUDY OBJECTIVES

The relationship of iron with cognitive and motor impairment in idiopathic rapid eye movement sleep behavior disorder (iRBD) remains unknown.

METHODS

Twenty-nine (29) patients and 28 healthy controls (HCs) underwent susceptibility weighted imaging and susceptibility mapping. These images were used to evaluate the nigrosome-1 (N1) sign in the substantia nigra (SN), global and regional high-iron (RII) content and volume of subcortical nuclei.

RESULTS

The number of iRBD patients with N1 loss (12) was significantly higher than HCs (2) (P=0.005). Compared with HCs, the iRBD patients had reduced volume of the right caudate nucleus (RCN) (P<0.05, FDR correction) but no significant changes in global and RII iron of the subcortical nuclei (all P>0.05, FDR correction). Multiple regression analysis revealed that: for cognitive function, the RII iron of the RCN was significantly correlated with visuospatial function and the global iron of the right dentate nucleus (RDN) was correlated with memory function; for motor function, the RII iron of the left DN (LDN) and global iron of the left CN correlated with the Alternate-Tap test (left, average), the global iron of the LDN correlated with the Alternate-Tap test (right), and the global iron of the left GP correlated with the 3-meter Timed Up and Go test (all P<0.05, FDR correction).

CONCLUSIONS

Our exploratory analysis found that iRBD patients had a higher incidence of N1 loss and reduced RCN volume after FDR correction. Cognitive and motor impairment were associated with iron deposition in several cerebral nuclei after FDR correction.

An interictal measurement of cerebral oxygen extraction fraction in MRI-negative refractory epilepsy using quantitative susceptibility mapping

PURPOSE

Oxygen extraction fraction (OEF) can be a factor to identify brain tissue's disability in epileptic patients. This study aimed to assess the OEF's level measurement in refractory epileptic patients (REPs) using a quantitative susceptibility mapping (QSM) method and to determine whether the OEF parameters change.

METHODS

QSM-OEF maps of 26 REPs and 16 healthy subjects were acquired using 3T MRI with a 64-channel coil. Eighteen regions-of-interest (ROIs) were chosen around the cortex in one appropriate slice of the brain and the mean QSM-OEF for each ROI was obtained. The correlations of QSM-OEF among different clinical characteristics of the disease, as well as between the patients and normal subjects, were also investigated.

RESULTS

QSM-OEF was shown to be significantly higher in REPs (44.9 ± 5.8) than that in HS (41.9 ± 6.2) (p < 0.05). Mean QSM-OEF was statistically lower in the ipsilateral side (44.5 ± 6.6) compared to the contralateral side (46.4 ± 6.8) (P < 0.01). QSM-OEF was illustrated to have a strong positive correlation with the attack duration (r = 0.6), and a moderate negative correlation with the attack frequency (r = -0.3). Using an optimized support vector machine algorithm, we could predict the disease in subjects having abnormal OEF values in the brain-selected-ROIs with sensitivity, specificity, AUC, and the precision of 0.96, 1, 0.98, and 1, respectively.

CONCLUSIONS

The results of this study revealed that QSM-OEF of the REPs' brain is higher than that of HS, which indicates that QSM-OEF is associated with disease activity.

Reduced regional cerebral venous oxygen saturation is a risk factor for the cognitive impairment in hemodialysis patients: a quantitative susceptibility mapping study

The purpose of this study was to noninvasively evaluate the changes of regional cerebral venous oxygen saturation (rSvO₂) in hemodialysis patients using quantitative susceptibility mapping (QSM) and investigate the relationship with clinical risk factors and neuropsychological testing. Fifty four (54) hemodialysis patients and 54 age, gender and education matched healthy controls (HCs) were recruited in this prospective study. QSM data were reconstructed from the original phase data of susceptibility weighted imaging to measure the susceptibility of cerebral regional major veins in all subjects and calculate their rSvO₂. The differences in rSvO₂ between hemodialysis patients and HCs were investigated using analysis of covariance adjusting for age and gender as covariates. Stepwise multiple regression and correlation analysis were performed between the cerebral rSvO₂ and clinical factors including neuropsychological testing. The SvO₂ of the bilateral cortical, thalamostriate, septal, cerebral internal and basal veins in hemodialysis patients was significantly lower than that in HCs (p < 0.001, Bonferroni corrected). The cerebral rSvO₂ in all these veins was reduced by 1.67% to 2.30%. The hematocrit, iron, glucose, pre-and post-dialysis diastolic blood pressure (DBP) were independent predictive factors for the cerebral rSvO₂ (all P < 0.05). The Mini-Mental State Examination and Montreal Cognitive Assessment (MoCA) scores were both lower in patients than those in HCs (both P < 0.05). The SvO₂ of the left cerebral internal vein correlated with MoCA scores (r = 0.492; P = 0.02, FDR corrected). In conclusion, our study indicated that the cerebral rSvO₂ was reduced in hemodialysis patients, which was the risk factor for neurocognitive impairment. The hematocrit, iron, glucose, pre-and post-dialysis DBP were independent risk factors for the cerebral rSvO₂.