The Application of Magnetic Susceptibility Separation for Measuring Cerebral Oxygenation in Preterm Neonates

¹ Integrated Sciences, The University of British Columbia, Vancouver, BC, Canada
² BC Children’s Hospital Research Institute, The University of British Columbia, Vancouver, BC, Canada
³ Pediatrics, The University of British Columbia, Vancouver, BC, Canada
⁴ Physics and Astronomy, The University of British Columbia, Vancouver, BC, Canada

^✉ Correspondence: Alexander Mark Weber <aweber@bcchr.ca>

Impact Statement

• This study evaluated the use of QSM and its paramagnetic components to measure cerebral oxygenation in neonates.
• By comparing susceptibility-derived oxygen saturation (SvO₂) in the superior sagittal sinus (SSS) and central cerebral veins (CCV), it adds to the field of neonatal cerebral oxygenation measurement.
• Decomposing QSM into paramagnetic components shows potential for improving SvO₂ accuracy, particularly in the SSS, though variability remains a challenge.
• The results suggest no significant oxygenation difference between the SSS and CCV, contrasting with previous findings, indicating a need for further research on neonatal venous oxygenation.

Category of Study: basic science

Abstract:

Background: Quantitative susceptibility mapping (QSM), a magnetic resonance imaging (MRI) modality sensitive to deoxyhemoglobin, is a promising method for measuring cerebral oxygenation in human neonates. Paramagnetic sources, like deoxyhemoglobin, however, can be obscured by diamagnetic sources such as water and myelin. This study evaluated whether QSM images, or isolated paramagnetic components, are more accurate for measuring oxygenation of cerebral veins of preterm neonates, and explored oxygenation differences between the major cerebral veins.

Methods: 19 preterm neonates were scanned on at term equivalent age on a 3T MRI using a multi-echo susceptibility-weighted imaging sequence. Susceptibility values were calculated from QSM images to determine oxygen saturation (SvO₂) in the superior sagittal sinus (SSS) and central cerebral veins (CCV). The paramagnetic components of QSM images were isolated, and SvO₂ values were recalculated.

Results: The mean SvO₂ values from QSM were 72.4% (SD, 3.4%) for the SSS and and 68.7% (SD, 3.5%) for the CCV. SvO₂ values for paramagnetic components were 58.1% (SD, 7.3%) for the SSS and 57.7% (SD, 7.0%) for the CCV.

Conclusion: While paramagnetic component decomposition yielded SSS values closer to those found in the literature, it increased variability. No significant oxygenation differences were found between the SSS and CCV, contrasting with prior studies.

Introduction

With advances in neonatal medical care, more infants born preterm are surviving into childhood¹. These children are at high risk of acquiring adverse neurodevelopmental outcomes when compared to their term-born peers². Irregularities in early cerebral oxygen levels have been identified as a potential source of such delays, where too little oxygen provided during NICU care can result in white matter injury, while too much oxygen can result in reduced cortical connectivity³. As such, being able to precisely, accurately, and non-invasively measure cerebral oxygenation is necessary for understanding and improving neurodevelopmental outcomes in preterm neonates.

Unfortunately, there exist many challenges in measuring cerebral oxygen metabolism in neonates. Cerebral metabolic rate of oxygen (CMRO₂) using oxygen-15 positron emission tomography (PET)⁴, has been measured in infants⁵, and is considered the gold standard. However, this method is invasive, requiring ionizing radiation, which limits its suitability for neonates. A less invasive option for evaluating brain hemodynamics is near-infrared spectroscopy (NIRS), which uses the attenuation of near-infrared light (~650–950 nm) as it passes through biological tissue⁶. Deoxygenated and oxygenated hemoglobin absorb this light differently, allowing NIRS to estimate changes in deoxyhemoglobin and oxyhemoglobin⁷ and thus provide an estimate of cerebral venous oxygen saturation (SvO₂). While NIRS offers the advantage of being non-invasive and continuous bedside monitoring, it is limited to regional assessments where the probe is placed and is sensitive only to superficial brain tissue due to the shallow penetration depth of near-infrared light⁸.

For the preceding reasons, non-invasive MRI-based techniques are actively being explored to assess regional and whole-brain blood oxygenation. While MRI-based methods have been developed for adults^9–11, their application in neonates is only beginning to be explored^12–16. This delay is likely due to the unique challenges posed by neonates, including their smaller anatomies, distinct hemodynamic profiles, susceptibility to motion artifacts, and the difficulties associated with recruiting this population for research. These methods have almost all relied on T2 relaxation to estimate CSvO₂^12–14,16 with the exception of¹⁵ which used susceptometry¹⁷. One limitation of these T2 relaxation methods, however, is the fact that SvO₂ is often measured using a single imaging slice, averaging values across several voxels, and only in the superior sagittal sinus (SSS). In the case of¹⁵, they obtained regional and whole-brain data, but with thick slices (5mm), and still only estimated CSvO₂ in the SSS. An alternative MRI method using quantitative susceptibility mapping (QSM) has been proposed, which can measure SvO₂ regionally and across the whole-brain at high resolution (< 1mm³ per voxel)¹⁸. However, this method left room for improvement, as it removed the SSS (averaging CSvO₂ across the internal veins), and required an arbitrary threshold value of 0.15 ppm in order to acquire realistic results¹⁸. Furthermore, QSM tends to underestimates parametric components due to the inclusion of diamagnetic tissue, and vice versa, as the opposing magnetic susceptibilities effectively subtract from one another¹⁹.

In the present study, we set out to determine whether decomposing the QSM image into its paramagnetic and diamagnetic components would allow for a more accurate assessment of SvO₂ in the central cerebral veins (CCV) of a cohort of preterm neonates. We also had a secondary aim of preserving the SSS vessel in our QSM images and using this data to determine whether a difference in oxygenation existed between the SSS and the CCV.

Methods

The study was approved by the Clinical Research Ethics Board at the University of British Columbia and Children’s & Women’s Hospital (H21-00655) and written informed consent was obtained from the parent/guardian for each infant.

Study population

Participant data comes from a previous study²⁰ (under review). Participants consisted of preterm neonates born between 25- and 31-weeks gestational age (GA) who were admitted to the level III NICU at BC Women’s Hospital. Recruitment took place over a span of one year, from February 2021 to January 2022, facilitated by a dedicated research nurse. Parents of eligible infants were approached by the research nurse prior to discharge from the NICU to explain the study objectives and seek their consent for participation. Infants meeting the criteria for inclusion were scanned for the study if they had already been discharged from the NICU, were in stable condition, and had reached a term equivalent age of 37 to 44 weeks GA. However, certain exclusion criteria were applied to ensure the homogeneity and integrity of the study sample: infants were excluded if there was clinical evidence of a congenital malformation or syndrome, a TORCH infection, or ultrasound evidence of large parenchymal hemorrhagic infarction (>2 cm, Grade 4 intraventricular hemorrhage).

Image acquisition

MR imaging was performed on a 3.0 Tesla General Electric Discovery MR750 scanner (scanner software version DV26.0_R03) equipped with a SREE Medical Systems (Cleveland, OH) single-channel neonatal head coil (Table 1). The scans were conducted at the BC Children’s MRI Research Facility. Prior to the scanning procedure, subjects were carefully prepared by a research nurse. Swaddling and feeding were used to ensure the comfort and cooperation of the subjects during the scan. Importantly, no sedatives or invasive markers were utilized throughout the procedure. Subjects were placed within a specially designed SREE Medical Systems MRI compatible incubator, which facilitated both safety and motion minimization. Molded foam was strategically positioned around the head and body within the incubator to further restrict subject movement. To protect against potential hearing damage, ear plugs were employed during the scanning process. Additionally, a pulse oximeter was affixed to the subject’s foot to monitor arterial oxygen saturation and heart rate throughout the scan.

Table 1. Technical parameters for MR imaging pulse sequences

	T1w	T2w	pcASL	SWI
Sequence	3D FSPGR	3D CUBE	Multi-shot 3D fast spin-echo	3D spoiled GRE flow-compensated
Acquisition plane	Coronal	Sagittal	Axial	Axial
Phase-encoding direction	Left-Right	Posterior-Anterior	Posterior-Anterior	Left-Right
TR (ms)	7.74	2,300	4,680	30.9
TE (ms)	2.97	66.29	10.55	5 echoes; first echo: 5; echo spacing: 5.24
Flip angle	12°	90°	111°	20°
FOV (cm)	20	20	24	25
Acquisition matrix	512 x 512	256 x 256	128 x 128	256 x 256
In-plane resolution (mm)	0.39 x 0.39	0.78 x 0.78	1.875 x 1.875	0.977 x 0.977
Slice thickness (mm)	1	1	4	2, reconstructed to 1 with zero filling (ZIP2)
Number of slices	126	106	50	92
Additional parameters	n/a	n/a	1,450 ms label period; 2,025 ms pulse label; 24 control-label pairs	n/a
Scan duration	4 min 39 s	5 min 1 s	5 min 26 s	5 min 29 s
T1w = T1-weighted; T2w = T2-weighted; pcASL = pseudo-continuous arterial spin labelling; SWI = susceptibility weighted imaging; FSPG = fast spoiled gradient echo; CUBE = General Electric name of sequence, not an acronym; GRE = gradient echo; ZIP2 = through-plane zero filling interpolation

Source: Article Notebook

The MRI scan protocol comprised of the following sequences (plane of acquisition in parentheses): a T1-weighted scan (coronal), a T2-weighted scan (sagittal), a pseudo-continuous arterial spin labeling (ASL) scan²¹ (axial), a multi-echo susceptibility-weighted imaging scan²² (axial), and a diffusion-weighted imaging (DWI) spin-echo echo planar imaging (EPI) sequence (axial). The DWI sequence was not used for the present study.

Image analysis

The raw DICOM files acquired from the scanning procedure were converted to NIfTI (Neuroimaging Informatics Technology Initiative) format using Chris Rorden’s dcmniix tool²³. SWI magnitude data files were then used to create subject-specific brain masks that would not erode the SSS during QSM processing, an issue faced by our group in the past¹⁸. A step-by-step summary of the pipeline used is shown in Figure 1.

Figure 1. Pipeline for generating subject-specific brain masks that include the superior sagittal sinus (SSS). Initial steps involved (1) creating brain a mask from the magnitude of the fifth echo of the susceptibility weighted scan. Subsequently, the brain mask is dilated and then (2) utilized in conjunction with a quantitative susceptibility mapping (QSM) script to generate a preliminary QSM image. Further refinement involved (3) segmenting the SSS from the QSM image manually to create a tissue mask of the SSS region. Finally, (4) the vascular mask of the SSS is integrated with the initial brain mask, forming the comprehensive brain mask essential for obtaining susceptibility data that includes the SSS.

Source: Python Figures

First, the fifth echo (TE = 25.96 ms) SWI magnitude file was processed using FSL’s (v. 6.0.7.3)²⁴ fslroi, fslmaths, and bet²⁵ to create a preliminary brain mask, similar to our previous efforts, which does not contain the SSS. The last echo is used to generate the brain mask as it reliably removes artifacts from air-tissue and bone-tissue interfaces, e.g., sinuses, without the need for manual erosion. The reason for this is at longer echo times, tissues with rapid signal decay (such as bone, air, and sinuses) lose their MRI signal due to dephasing caused by greater magnetic field inhomogeneities. Fslroi was used to isolate the fifth echo of the magnitude data, which was then squared using fslmaths and the option -sqr. Squaring the magnitude image was found to dramatically improve subsequent brain extraction. The resulting image was then used to create the preliminary brain mask using bet with the options -m and -R. The former flag generated a binary brain mask, while the latter performed a more robust brain centre estimation. The brain mask was then dilated by 7 voxels using fslmaths and the options -kernel boxv and -dilM in order for the dilated mask to contain the SSS (along with unwanted tissue as well). This mask was then used, along with the phase images, in a MATLAB script for QSM calculation from Christian Kames²⁶ to produce a preliminary QSM image that contained the SSS, albeit with fairly low signal-to-noise ratio and other unwanted tissue. Given the high contrast in voxel intensity between the SSS and surrounding tissue, the select by intensity tool in FSLeyes²⁷ was then used to segment the SSS from the QSM image and create a 3D mask of the selected region. Using fslmaths and the options -add and -bin, the SSS mask was then combined with the original brain mask of the fifth echo. This resulted in a brain mask that contained only brain and SSS signal. Finally, this mask was used in a final QSM post-processing step to create a QSM image that includes the SSS while maintaining a high signal-to-noise ratio, making it suitable to obtain accurate susceptibility values.

STI Suite (v. 3.0)²⁸, was used to process the final QSM images as it produced the images with the least amount of artifacts (based on a visual assessment by the authors) without eroding the SSS. The finalized brain mask along with all five echoes of the magnitude and phase images were used in STI Suite along with the following parameters: 0.9766 x 0.9766 x 1 mm³ voxel size, 5 ms TE1, 5.24 ms \(\Delta\)TE, and 77.4 ms sum TE, B0 strength = 3, and B0 direction = [0, 0, 1]. The 3D GRE data option was selected for the phase processing stage, and STAR-QSM was selected for the QSM stage. STAR-QSM outputs a single QSM map for each echo (i.e. five total). The last three echoes of the QSM maps were then averaged to create the final QSM image using fslmaths, as the accumulation of phase due to susceptibility is small in early echoes and artifacts dominate the phase²⁹. Finally, the `select by intensity’ tool in FSLeyes was then used to semi-automatically make vascular masks of the SSS and CCV from each subject’s QSM image (Figure 2). The vascular masks were used to calculate the mean susceptibility of each subject’s SSS and CCV from their QSM image with fslstats.

Figure 2. Sample venous masks. A sample superior sagittal sinus (red) and central cerebral vein mask (blue) displayed in coronal, sagittal, and axial view. The QSM image is used as the underlay. y, x and z values represent the slice number in each plane (coronal, sagittal, and axial, respectively).

To isolate the paramagnetic component of subjects’ QSM data, the \(\chi\)-separation toolbox³⁰ from the Laboratory for Imaging Science and Technology was used. All five echoes of each subject’s magnitude and phase SWI data were used along with the following parameters: 0.9766 x 0.9766 x 1 mm³ voxel size; TE (ms) = [5, 10.24, 15.48, 20.72, 26.96]; \(\Delta\)TE (ms) = 5.24; B0 strength = 3; B0 direction = [0, 0, 1]. The \(\chi\)-separation toolbox outputs a single negative (diamagnetic), positive (paramagnetic), and total \(\chi\) map (QSM). The mean susceptibility of each subject’s SSS and CCV in their paramagnetic maps was calculated with the same vascular masks used for the QSM images. Sample images showing the magnitude, final QSM, and final paramagnetic component images are shown in Figure 3.

Figure 3. An example of subject imaging data. A sample coronal, sagittal, and axial slice is displayed for each image. (a) The 1st echo of the magnitude susceptibility weighted imaging sequence; (b) the final quantitative susceptibility mapping image; and (c) the paramagnetic component isolated from the quantitative susceptibility map. The color bars in (b) and (c) indicates the range of susceptibility \(\chi\) values in parts per million. y, x and z values represent the slice number in each plane (coronal, sagittal, and axial, respectively).

Once the mean susceptibility values of the SSS and CCV were obtained from the subjects’ QSM images or paramagnetic maps, venous oxygen saturation (SvO₂) was calculated with the following equation³¹:

\[ SvO_{2} = 1 - \frac{\Delta \chi _{blood} - (\Delta \chi _{oxy} * Hct)}{\Delta \chi _{do} * Hct} \tag{1}\]

where \(\Delta \chi _{blood}\) is the vessel’s measured susceptibility, \(\Delta \chi _{oxy}\) is the constant representing the susceptibility changes of oxygenated blood relation to water, \(\Delta \chi _{do}\) is the difference in susceptibility between oxygenated and deoxygenated blood, and Hct is the subject’s hematocrit. \(\Delta \chi _{oxy}\) was -0.21 * 4\(\pi\) ppm as per³² and³³, while \(\Delta \chi _{do}\) was -0.03 * 4\(\pi\) ppm as per³⁴. Subjects’ Hct for the day of the scan was calculated using a four-parameter Weibull function with previously measured values while still in the NICU.

Statistical analysis

Statistical analysis was performed using R and RStudio (v. 2023.09.1 Build 494)^35,36. Mean and standard deviation values are reported for most statistics, unless specified otherwise. A paired Student’s t-test was used to determine statistical significance (p <0.05) between two parameters (e.g. \(\chi\) values between venous structures).

Results

A total sample size of 19 infants were scanned, with a mean (\(\pm\) standard deviation) gestational age of 28.8 \(\pm\) 1.68 weeks and a mean birth weight of 1276.05 \(\pm\) 294.87 grams. A comprehensive summary of neonatal characteristics, including additional demographic and clinical data, is provided in Table 2 for reference.

Table 2. Demographic and clinical characteristic of the study sample.

Variable	Subject data (n = 19)
Gestational age, weeks (mean ± SD)	28.8 ± 1.68
Corrected gestational age on scan day, weeks (mean ± SD)	40.36 ± 1.4
Number of male neonates (%)	10 (52.63)
Birth weight, g (mean ± SD)	1276.05 ± 294.87
Weight on scan day, g (mean ± SD)	3396.58 ± 597.72
Days spent in NICU (median, IQR)	53, 23
Days on ventilation (median, IQR)	31, 28.5
SD = standard deviation; IQR = inter quartile range

The mean SvO₂ values for the SSS and the CCV were found to be 0.72 \(\pm\) 0.03% and 0.69 \(\pm\) 0.03%, respectively, when determined from the QSM data. When determined from the paramagnetic map, the mean SvO₂ values for the SSS and the CCV were found to be 0.58 \(\pm\) 0.07% and 0.58 \(\pm\) 0.07%, respectively. A summary of the measured physiological parameters, including the chi values used to calculate SvO₂, can found in Table 3.

Table 3. Summary of acquired physiological parameters. Mean \(\pm\) SD is shown for chi and SvO₂ values. The P-value and 95% confidence interval (CI) were obtained through the comparison of values between QSM and paramagnetic maps; (n=19).

Region	Measure	QSM	Paramagnetic map	p-value	95% CI
SSS	Chi (ppm)	0.1 ± 0.02	0.21 ± 0.05	2.84e-11	-0.13, -0.09
SSS	SvO₂ (%)	72.46 ± 3.43	58.14 ± 7.3	6.12e-10	0.12, 0.17
CCV	Chi (ppm)	0.13 ± 0.02	0.22 ± 0.05	6.25e-09	-0.1, -0.07
CCV	SvO₂ (%)	68.71 ± 3.46	57.69 ± 6.97	2.16e-09	0.09, 0.13
QSM = quantitative susceptibility mapping; CI = confidence interval; SSS = superior sagitall sinus; CCV = central cerebral vein

Region-specific \(\chi\) and SvO₂ values acquired from QSM were compared to values acquired from paramagnetic maps. In both the SSS and CCV, it was found that a significant difference existed between values acquired (\(\chi\) and SvO₂) from QSM and paramagnetic maps (p < 0.05). Violin plots of the comparisons are shown in Figure 4.

Figure 4. Vein-specific susceptibility and oxygen saturation values by method. (A, B) contains violin plots comparing subject chi (ppm) acquired from the cerebral central veins; (C, D) contains violin plots comparing subject SvO2 (%) acquired from the superior sagittal sinus. Raw data points from paramagnetic maps are shown as filled green circles and raw data points from QSM are shown as filled orange circles. Each line connects the raw data points of a single subject. (***) indicates P<0.05.

Source: gavin_thesis_R

The acquired \(\chi\) and SvO₂ values were additionally compared between veins. In data created from QSM, a significant difference was found between the CCV and SSS in mean \(\chi\) (p < 0.05; 95% CI [0.017, 0.04]) and mean SvO₂ (p < 0.05; 95% CI [-0.052, -0.023]). In data acquired from paramagnetic maps, no significant difference was observed between the CCV and the SSS in either mean \(\chi\) (p = 0.711; 95% CI [-0.02, 0.029]) or mean SvO₂ (p = 0.752; 95% CI [-0.034, 0.029]). A summary of these comparisons is represented in Figure 5.

Figure 5. Inter-venous comparisons of susceptibility and oxygen saturation. Violin plots comparing (A, C) chi (ppm) and (B, D) SvO2 (%) between the CCV and the SSS. Panels (A) and (B) used data acquired from QSM, and its raw data points are shown as filled orange circles. Panels (C) and (D) used data acquired from paramagnetic maps, and its raw data points are shown as filled green circles. Each line connects the raw data points of a single subject. (***) indicates p<0.05; (n.s.) indicates no significant difference.

Discussion

The primary objective of the present study was to assess whether the application of magnetic susceptibility separation to neonatal QSM data could provide more accurate SvO₂ measurements without the need for an arbitrary threshold value. To the best of our knowledge, we are the first to test this in a neonatal cohort, as susceptibility separation has been typically evaluated as a method of imaging myelin and brain iron in adult subjects^30,37. Our results showed that the SvO₂ values of the SSS and CCV obtained from susceptibility separation are significantly lower than the respective SvO₂ values obtained from QSM alone. When our results were compared to the literature (see below), we found that our SSS SvO₂ data from susceptibility separation agreed well with the findings of other studies measuring SvO₂ of the SSS in similar subject populations. Conversely, the paramagnetic CCV SvO₂ data saw less agreement with the existing literature than the corresponding data from QSM. However, there is reason to believe our paramagnetic CCV values may be accurate given their similarity to the paramagnetic SSS values and the limitations of the two studies that observed CCV SvO₂. Additionally, it is important to note that our SvO₂ measurements from susceptibility separation had greater variance than our measurements from QSM, indicating a limitation that should be addressed in future research. Overall, the present work demonstrates the promise of susceptibility separation as an MRI post-processing technique that can measure the oxygenation of the cerebral veins of infant subjects, a useful marker of regional oxygen consumption in the brain.

Comparison with literature values

To evaluate the validity of our results, we compared the mean SvO₂ values we obtained through QSM and susceptibility separation to the mean SvO₂ values found by MRI studies investigating the oxygenation of the SSS or the CCV. Notably, the number of studies that measure the SvO₂ of the CCV, or any of its individual veins, in infants is fairly lower than the number of studies investigating the oxygenation of the SSS. Our comparison is summarized in Table 4.

Table 4. Cerebral oxygenation values of neonates and fetuses in the literature.

Region	Study	Subjects	Method	SvO2 (%)
Whole-brain	Skov et al. (1993)	Preterm neonates (n=9)	NIRS	53.4 ± 15.4
Whole-brain	Skov et al. (1993)	Asphyxiated term neonates (n=10)	NIRS	67.3 ± 9.4
Whole-brain	Altman et al. (1993)	Preterm and term neonates with HIE (n=11)	PET	21.6 ± 21.0
SSS	Gou et al. (2024)	Healthy neonates (n=37)	MRI: TRUST	66.7 ± 4.9
SSS	Jiang et al. (2019)	Healthy neonates (n=4)	MRI: aTRUPC	64.8 ± 2.9
SSS	De Vis et al. (2014)	PT-TEA neonates (n=18)	MRI: T2-TRIR	52.0 ± 12.0
SSS	Yadav et al. (2020)	Late third trimester fetuses (n=33)	MRI: Susceptometry	58.6 ± 4.8
SSS	This study	PT-TEA neonates n=19	MRI: QSM	72.46 ± 3.43
SSS	This study	PT-TEA neonates n=19	MRI: χ-separation	58.14 ± 7.3
CCV	Weber et al. (2021)	Preterm neonates with HIE (n=8)	MRI: QSM	72.2 ± 6.0
CCV	Weber et al. (2021)	Healthy neonates (n=8)	MRI: QSM	73.6 ± 2.8
CCV	Jiang et al. (2019)	Healthy neonates (n=4)	MRI: aTRUPC	70.2 ± 3.3
CCV	This study	PT-TEA neonates n=19	MRI: QSM	68.71 ± 0.03
CCV	This study	PT-TEA neonates n=19	MRI: χ-separation	57.69 ± 6.97
PT-TEA = born preterm and scanned at term-equivalent age; late third trimester = >35 weeks gestational age; HIE = hypoxic-ischemic encephalopathy; TRUST = T2-relaxation-under-spin tagging; aTRUPC = accelerated T2-relaxation-under-phase-contrast; T2-TRIR = T2-tissue-relaxation-inversion-recovery; SWI = susceptibility weighted imaging.

As shown in Table 4, the infants observed in MRI studies investigating cerebral vein oxygenation noticeably differ in clinical status, with three studies involving healthy neonates^16,18,38, three studies (including the present study) involving preterm neonates^12,18, and one study involving late third trimester fetuses³². In the studies involving healthy neonates, the SvO₂ of the SSS fell within the range of 64.8% – 66.6%^16,38, while the SvO₂ of the CCV fell within the range of 70.2% - 73.6%^16,18. Notably, the SvO₂ value of the SSS we obtained from susceptibility separation (58.14%) was closest to values obtained from the studies involving late third trimester fetuses³⁹ or pre-term neonates¹², each finding an SSS SvO₂ value of 58.6% and 52.0%, respectively. It is important to note the difference in MRI modalities used to obtain these values. For their study,³⁹ used MR susceptometry, which involves measuring the difference in phase between the chosen vessel and its background in imaging data from an SWI scanning sequence³⁹. In¹², the authors used T2-TRIR, which allowed them to determine the transverse relaxation rate of blood within the vessel, which can be used alongside hematocrit data to estimate SvO₂. Additionally, the GA of infants scanned in our study ranged between 37 and 44 weeks, while the GA of the fetuses scanned in³⁹ was ≥35 weeks and the GA of infants scanned in¹² ranged between 38 and 40 weeks. As such, our SSS SvO₂ values found through susceptibility separation show promise given their similarity to the SvO₂ values found by³⁹ and¹², two studies that involved comparable subject populations and used considerably different methods.

Conversely, the SvO₂ value of the CCV we obtained through QSM (68.71%) was closest to values from similar studies in the literature. In¹⁸, QSM was used to measure an SvO₂ of 71.5% in preterm neonates with HIE and an SvO₂ of 73.6% in healthy neonates. In their study,¹⁶ also involved healthy neonates and obtained an SvO₂ of 70.2% through an accelerated TRUPC sequence. In contrast, the SvO₂ of the CCV we obtained through susceptibility separation was 57.69%. This disparity from the literature, however, may not undermine the value we obtained, as the study design of¹⁸ and¹⁶ may prevent their values from being representative of the study demographic. In¹⁸, the authors utilized an arbitrary 0.15 ppm threshold and included all \(\chi\) values above 0.15 when measuring the mean \(\chi\) of the CCV, which potentially led to the introduction of \(\chi\) from veins outside the CCV. In¹⁶, the authors acquired their data from 4 subjects, a notably small sample size. Given the limitations of the existing literature and the similarity of the mean paramagnetic CCV SvO₂ value (57.69%) to the mean paramagnetic SSS SvO₂ value (58.14%), it is plausible that susceptibility separation provides more accurate measurements of oxygenation in both cortical and subcortical veins. One reason for this is due to its ability to mitigate partial volume effects, which are likely to contaminate other methods resulting in inaccurate CSvO₂ values³⁰.

Another notable distinction between our findings and those of the existing literature was that we observed no significant oxygenation difference between the SSS and the CCV when \(\chi\) was derived from paramagnetic maps.¹⁶, the only other study that also measured SvO₂ in both the SSS and CCV, observed significantly lower oxygenation in the SSS (64.8%) when compared to the CCV (70.2%). Given the small sample size utilized by¹⁶, it is difficult to ascertain whether this is generalizable to all neonates.

Limitations and future directions

This study has a few limitations that should be considered for future research. Firstly, only 19 infants were recruited for scanning. Given the emotional toll placed on parents when their child is born preterm, it is understandable that they may show reluctance in consenting to further testing that is not medically necessary. Obtaining a larger sample size in future studies, however, may provide greater insight into the efficacy of susceptibility separation. Secondly, this study did not include healthy neonates born at term, resulting in a lack of a control cohort. This is because recruiting healthy controls when there is no contraindication is very difficult. The addition of such a group may provide further validity to any findings and may reveal potential differences in cerebral oxygen consumption between term and preterm neonates. Finally, future studies should consider the use of multi-echo T2 imaging data when performing the decomposition of QSM images. The toolbox applied by this study for QSM decomposition³⁰ utilizes R2 data, which can be obtained from multi-echo T2 imaging. Our study protocol involved the collection of multi-echo SWI imaging data, and as such, we could only use R2^* data to perform the decomposition. Furthermore, this may account for the reduced precision of SvO₂ values obtained through susceptibility separation.

Conclusion

This study aimed to evaluate how the use of susceptibility separation on preterm neonatal QSM images could be used in determining the oxygenation of cerebral venous vessels. We compared venous specific SvO₂ values obtained from QSM images and their respective paramagnetic components to SvO₂ values from neonatal MRI studies. We found that susceptibility separation provided SvO₂ values of the SSS that were comparable to values found in the literature, providing evidence that this processing technique may be a valid tool for measuring regional cerebral oxygen consumption. Additionally, we were able to simultaneously measure SvO₂ in both the SSS and CCV, which permitted us to observe no difference in oxygenation between the two vessels when considering data from isolated paramagnetic components. Ultimately, we hope our work inspires future studies that seek to explore and improve the capabilities of magnetic susceptibility separation, culminating in the development of a tool for clinicians and researchers alike.

Data availability

The manuscript was written in a ‘reproducible manner’. The entire manuscript, including statistics reported, figures, and tables, can be reproduced here: weberlab.github.io/Chisep_CSVO2_Manuscript

Unfortunately, we can not upload our MRI images to an open repository as we did not obtain permission in our consent forms.

References

1.

McKenzie, K., Lynch, E. & Msall, M. E. Scaffolding Parenting and Health Development for Preterm Flourishing Across the Life Course. Pediatrics 149, e2021053509K (2022).

2.

Twilhaar, E. S. et al. Cognitive Outcomes of Children Born Extremely or Very Preterm Since the 1990s and Associated Risk Factors: A Meta-analysis and Meta-regression. JAMA Pediatrics 172, 361 (2018).

3.

Rantakari, K. et al. Early oxygen levels contribute to brain injury in extremely preterm infants. Pediatric Research 90, 131–139 (2021).

4.

Mintun, M. A., Raichie, M. E., Martin, W. R. W. & Herscovitch, P. Brain oxygen utilization measured with O-15 radiotracers and positron emission tomography. Journal of Nuclear Medicine 25, 177–187 (1984).

5.

Altman, D. I. et al. Cerebral blood flow requirement for brain viability in newborn infants is lower than in adults. Annals of Neurology 24, 218–226 (1988).

6.

Skov, L., Pryds, O., Greisen, G. & Lou, H. Estimation of cerebral venous saturation in newborn infants by near infrared spectroscopy. Pediatric Research 33, 52–55 (1993).

7.

Wray, S., Cope, M., Delpy, D. T., Wyatt, J. S. & Reynolds, E. O. R. Characterization of the near infrared absorption spectra of cytochrome Aa3 and haemoglobin for the non-invasive monitoring of cerebral oxygenation. Biochimica et Biophysica Acta (BBA) - Bioenergetics 933, 184–192 (1988).

8.

Boas, D. A., Dale, A. M. & Franceschini, M. A. Diffuse optical imaging of brain activation: Approaches to optimizing image sensitivity, resolution, and accuracy. NeuroImage 23, S275–S288 (2004).

9.

Jain, V. et al. Rapid magnetic resonance measurement of global cerebral metabolic rate of oxygen consumption in humans during rest and hypercapnia. Journal of Cerebral Blood Flow & Metabolism 31, 1504–1512 (2011).

10.

Lu, H. & Ge, Y. Quantitative evaluation of oxygenation in venous vessels using T2-Relaxation-Under-Spin-Tagging MRI. Magnetic Resonance in Medicine 60, 357–363 (2008).

11.

Xu, F., Ge, Y. & Lu, H. Noninvasive quantification of whole-brain cerebral metabolic rate of oxygen (CMRO2) by MRI. Magnetic Resonance in Medicine 62, 141–148 (2009).

12.

De Vis, J. B. et al. Non-invasive MRI measurements of venous oxygenation, oxygen extraction fraction and oxygen consumption in neonates. NeuroImage 95, 185–192 (2014).

13.

Liu, P. et al. Quantitative assessment of global cerebral metabolic rate of oxygen (CMRO2) in neonates using MRI. NMR in biomedicine 27, 332–340 (2014).

14.

Qi, Y., Liu, P., Lin, Z., Lu, H. & Wang, X. Hemodynamic and Metabolic Assessment of Neonates With Punctate White Matter Lesions Using Phase-Contrast MRI and T2-Relaxation-Under-Spin-Tagging (TRUST) MRI. Frontiers in Physiology 9, 233 (2018).

15.

Jain, V. et al. Cerebral oxygen metabolism in neonates with congenital heart disease quantified by MRI and optics. Journal of Cerebral Blood Flow and Metabolism: Official Journal of the International Society of Cerebral Blood Flow and Metabolism 34, 380–388 (2014).

16.

Jiang, D. et al. Vessel-specific Quantification of Neonatal Cerebral Venous Oxygenation. Magnetic resonance in medicine 82, 1129–1139 (2019).

17.

Jain, V., Langham, M. C. & Wehrli, F. W. MRI estimation of global brain oxygen consumption rate. Journal of Cerebral Blood Flow and Metabolism: Official Journal of the International Society of Cerebral Blood Flow and Metabolism 30, 1598–1607 (2010).

18.

Weber, A. M., Zhang, Y., Kames, C. & Rauscher, A. Quantitative Susceptibility Mapping of Venous Vessels in Neonates with Perinatal Asphyxia. American Journal of Neuroradiology 42, 1327–1333 (2021).

19.

Kim, W. et al. \(\chi\)-Separation Imaging for Diagnosis of Multiple Sclerosis versus Neuromyelitis Optica Spectrum Disorder. Radiology 307, e220941 (2023).

20.

Zhu, A. et al. Regional cerebral metabolic rate of oxygen and levels of respiratory support in preterm neonates. AJNR: American Journal of Neuroradiology (2024).

21.

Alsop, D. C. et al. Recommended implementation of arterial spin-labeled perfusion MRI for clinical applications: A consensus of the ISMRM perfusion study group and the European consortium for ASL in dementia: Recommended Implementation of ASL for Clinical Applications. Magnetic Resonance in Medicine 73, 102–116 (2015).

22.

Denk, C. & Rauscher, A. Susceptibility weighted imaging with multiple echoes. Journal of magnetic resonance imaging: JMRI 31, 185–191 (2010).

23.

Li, X., Morgan, P. S., Ashburner, J., Smith, J. & Rorden, C. The first step for neuroimaging data analysis: DICOM to NIfTI conversion. Journal of Neuroscience Methods 264, 47–56 (2016).

24.

Woolrich, M. W. et al. Bayesian analysis of neuroimaging data in FSL. NeuroImage 45, S173–S186 (2009).

25.

Smith, S. M. Fast robust automated brain extraction. Human Brain Mapping 17, 143–155 (2002).

26.

Kames, C., Wiggermann, V. & Rauscher, A. Rapid two-step dipole inversion for susceptibility mapping with sparsity priors. NeuroImage 167, 276–283 (2018).

27.

McCarthy, P. FSLeyes. (2023) doi:10.5281/zenodo.8376979.

28.

Li, W., Avram, A. V., Wu, B., Xiao, X. & Liu, C. Integrated Laplacian-based phase unwrapping and background phase removal for quantitative susceptibility mapping. NMR in Biomedicine 27, 219–227 (2014).

29.

Zhang, Y., Rauscher, A., Kames, C. & Weber, A. M. Quantitative Analysis of Punctate White Matter Lesions in Neonates Using Quantitative Susceptibility Mapping and R2* Relaxation. AJNR. American journal of neuroradiology 40, 1221–1226 (2019).

30.

Shin, H.-G. et al. \(\chi\)-separation: Magnetic susceptibility source separation toward iron and myelin mapping in the brain. NeuroImage 240, 118371 (2021).

31.

Berg, R. C., Preibisch, C., Thomas, D. L., Shmueli, K. & Biondetti, E. Investigating the Effect of Flow Compensation and Quantitative Susceptibility Mapping Method on the Accuracy of Venous Susceptibility Measurement. 2021.04.14.439812 (2021) doi:10.1101/2021.04.14.439812.

32.

Portnoy, S., Milligan, N., Seed, M., Sled, J. G. & Macgowan, C. K. Human umbilical cord blood relaxation times and susceptibility at 3 T: Human Umbilical Cord Blood Relaxation Times and Susceptibility at 3 T. Magnetic Resonance in Medicine 79, 3194–3206 (2018).

33.

Sedlacik, J., Rauscher, A. & Reichenbach, J. R. Obtaining blood oxygenation levels from MR signal behavior in the presence of single venous vessels. Magnetic Resonance in Medicine 58, 1035–1044 (2007).

34.

Weisskoff, R. M. & Kiihne, S. MRI susceptometry: Image-based measurement of absolute susceptibility of MR contrast agents and human blood: COMMUNICATIONS. Magnetic Resonance in Medicine 24, 375–383 (1992).

35.

R Core Team. R: A Language and Environment for Statistical Computing. (2022).

36.

RStudio Team. RStudio: Integrated Development Environment for R.

37.

Ahmed, M. et al. The diamagnetic component map from quantitative susceptibility mapping (QSM) source separation reveals pathological alteration in Alzheimer’s disease-driven neurodegeneration. NeuroImage 280, 120357 (2023).

38.

Gou, Y. et al. Automatic Rejection based on Tissue Signal (ARTS) for motion-corrected quantification of cerebral venous oxygenation in neonates and older adults. Magnetic Resonance Imaging 105, 92–99 (2024).

39.

Yadav, B. K. et al. Imaging putative foetal cerebral blood oxygenation using susceptibility weighted imaging (SWI). European Radiology 28, 1884–1890 (2018).

Acknowledgments

We wish to acknowledge the work of and thank Victoria Tapics (Research Nurse); Vicki Goh (Research Nurse); Chacko Anil (Neonatologist); and Michael A. Sargent (Radiologist).

Author contributions

TGC wrote the original draft, performed the formal analysis, and contributed to methodology, validation, and visualization. AR helped with writing, reviewing & editing. REG helped with writing, reviewing & editing, and with initial funding acquisition. AMW was involved in project administration, supervision, validation, visualization, resources, methodology, formal analysis, funding acquisition, writing, reviewing & editing, conceptualization, data curation, and investigation.

Funding

Authors AMW and REG were co-primary applicants for a BC Children’s Hospital Research Institute - Brain, Behaviour and Development Catalyst Grant ($20,000). AMW was supported by an establishment award from BCCHRI. Scanning was partly funded through a special award to AMW from BCCHRI .

Competing interests

The authors have no competing interests to declare.

Consent statement

The study was approved by the Clinical Research Ethics Board at the University of British Columbia and Children’s & Women’s Hospital (H21-00655) and written informed consent was obtained from the parent/guardian for each infant.