Abstract
Silent cerebral infarcts (SCI) are common in patients with sickle cell disease (SCD) and are thought to be caused by a mismatch between oxygen delivery and consumption. Functional cerebrovascular shunting is defined as reduced oxygen offloading due to the rapid transit of blood through the capillaries caused by increased flow and has been suggested as a potential mechanism underlying reduced oxygenation and SCI. We investigated the venous arterial spin labeling signal (VS) in the sagittal sinus as a proxy biomarker of cerebral functional shunting, and its association with hemodynamic imaging and hematological laboratory parameters. We included 28 children and 38 adults with SCD, and ten healthy racematched adult controls. VS, cerebral blood flow (CBF), velocity in the brain feeding arteries, oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2) were measured before and after acetazolamide administration. VS was higher in patients with SCD compared to controls (P<0.01) and was increased after acetazolamide administration in all groups (P<0.01). VS was primarily predicted by CBF (P<0.01), but CBF-corrected VS was also associated with decreased CMRO2 (P<0.01). Additionally, higher disease severity defined by low hemoglobin and increased hemolysis was associated with higher CBF-corrected VS. Finally, CMRO2 was negatively correlated with fetal hemoglobin, and positively correlated with lactate dehydrogenase, which could be explained by changes in oxygen affinity. These findings provide evidence for cerebral functional shunting and encourage future studies investigating the potential link to aberrant capillary exchange in SCD.
Introduction
Sickle cell disease (SCD) is characterized by chronic hemolytic anemia resulting in organ damage including silent cerebral infarcts (SCI).1,2 SCI are associated with cognitive impairment at an early stage in life, which can result in under- and unemployment as well as lower quality of life.2,3 SCI are primarily localized in deep white matter and are hypothesized to be a result of impaired oxygen delivery.4–6 In order to maintain brain oxygenation, cerebral blood flow (CBF) is increased in patients with SCD.7–9 However, despite the preservation of normal global oxygen delivery at rest,10 the incidence of SCI continues to rise throughout the lifetime of these patients and is prevalent in more than 50% of adult SCD patients by the age of 32 years.11 One potential pathophysiological mechanism that has been proposed suggests that SCI result from localized impairments in oxygen delivery caused by cerebral functional shunting.10,12,13 Cerebral functional shunting is defined as the rapid transit of blood through the brain capillaries as a result of increased flow, limiting offloading of oxygen to the tissue.13 However, functional shunting is difficult to assess and more insight into this process is required to understand the underlying cerebral hemodynamics in SCD.
Arterial spin labeling (ASL) magnetic resonance imaging (MRI) is widely used in research settings as a non-invasive technique for cerebral perfusion measurements in patients with SCD.9,14,15 ASL applies radiofrequency pulses in the brain feeding arteries that change the magnetization of blood, also known as labeling. Following labeling, the blood travels to the brain and a labeled image is acquired. The same procedure is repeated without labeling to acquire a control image. By subtracting the labeled image from the control image, a quantitative perfusion map is obtained. A requirement of ASL is that images are acquired after a predefined delay (post-label delay), allowing the labeled blood to travel to the capillary bed, and for the labeled spins to enter the tissue. The difference between labeled spins detected in arterial and venous circulations represents spins that have exchanged with the brain parenchyma.16 According to the Renkin-Crone model, which relates blood flow to extraction fraction, the amount of unexchanged labeled spins in the venous outflow, i.e., the venous signal (VS), will increase approximately proportionally with CBF.17,18 However, in the presence of cerebral functional shunting, a larger volume of the labeled blood will pass unexchanged through the capillaries and may arrive on the venous side at the time of imaging. Indeed, a previous study demonstrated that, with longer post-label delays, the VS in the sagittal sinus can be used to inform on the exchange of spins at the capillary level.16 As such, VS may provide additional hemodynamic information about oxygen utilization at the capillary level, which may be affected in patients with SCD.
Support for this hypothesis comes from prior work in SCD patients.13,19,20 Juttukonda et al.13,20 observed the presence of ASL signal in the venous sinuses at a regular post-label delay (1,900 ms), which was associated with increased CBF and elevated velocities in the carotid arteries. When administered to adult patients with SCD, acetazolamide (ACZ; a compound that induces vasodilation and increased CBF) led to observed reductions in oxygen extraction fraction (OEF) as well as the cerebral metabolic rate of oxygen (CMRO2),10 indicating that oxygen consumption can actually deteriorate despite increasing oxygen delivery. Moreover, in SCD patients with SCI, an inverse relation between VS scores and OEF was found, providing support for a mechanism involving functional shunting.20 However, whether VS relates to cerebral functional shunting and could therefore be a marker of capillary oxygen exchange efficiency in SCD patients, remains to be investigated.
Therefore, in the present study, we investigated cerebral functional shunting in patients with SCD, by examining the relationship between VS and cerebrovascular imaging parameters of perfusion, cerebral oxygen extraction and metabolism as well as laboratory measures. In order to directly probe the VS-CBF relationship, we studied the imaging parameters both before and after ACZ administration in adult participants.
Methods
The data were obtained from two studies6,21 (one adult and one pediatric study) that were approved by the medical ethics committee of Academic Medical Center in Amsterdam and performed in accordance with the Declaration of Helsinki. Written informed consent was obtained from all adults and from parents/legal guardians of participants in the pediatric study cohort. Sixty-six patients with SCD (59 HbSS and 7 HbSβ0thalassemia) were recruited from hematology outpatient clinics.6,21 Additionally, ten healthy adult race-matched controls (8 HbAA and 2 HbAS) were recruited. Exclusion criteria were MRI contraindications, history of cerebral pathology, sickle cell crisis at inclusion, hospitalization 1 month prior to the study day, pregnancy, and ACZ contraindications for adult participants.
Data acquisition
Images were acquired on two 3.0 T MRI systems (Philips Intera and Ingenia, Philips Healthcare, Best, the Netherlands) with an 8-channel head coil for pediatric patients and a 32-channel head coil for adults. Prior to MRI, blood was drawn to quantify hematological laboratory parameters (for details, see the Online Supplementary Appendix). Anatomical sequences included a 2D T2-weighted and a 2D fluid-attenuated inversion recovery (FLAIR) scan for children and 3D FLAIR for adult participants were used for lesion and gray matter segmentation, and registration purposes.
For VS and CBF measurements, pseudo-continuous arterial spin labeling (pCASL) was used (Table 1). Adult participants received intravenous ACZ (16 mg/kg with a maximum of 1,400 mg) 10 minutes prior to the second ASL scan.
In order to obtain venous saturation for OEF measurements in adult participants’ pre- and post-ACZ, T2-prepared tissue Relaxation with an Inversion Recovery MRI (T2-TRIR) sequence22 was used. Finally, to obtain velocity measurements in the brain feeding arteries, a 2D phase-contrast single-shot gradient-echo T1 FFE sequence was acquired both pre- and post-ACZ (Table 1) (for more details on the acquisition parameters, see the Online Supplementary Appendix).
Data analysis
VS was assessed both qualitatively and quantitatively. Qualitatively, ASL images were visually inspected by an observer (LA) for the presence or absence of the venous ASL signal. Quantitatively, we scaled the perfusion-weighted signal (control-label) from the ASL images, by the group-based T1 of blood (patients = 1,818 ms23 and controls = 1,650 ms24), subject-specific estimates of labeling efficiency based on velocity,25 and M0 in ExploreASL.26 The superior sagittal sinus and straight sinus were manually segmented in three group-average (children and adults with SCD and healthy control group) images of all participants in standard MNI space27 (Figure 1). The segmented average image of each group was resampled to the native ASL space of each individual and used as an ASL mask from which the average VS is measured. In case the segmented venous sinuses were not correctly aligned after resampling, manual correction was applied to only include the VS. In addition, subject-specific gray matter masks were used in the native space to measure the mean signal in the GM. For the comparison between children and adults, the proportion of signal in the venous regions of interest (ROI) relative to GM was calculated in the VS images (VS/GM). This ratio is referred to as VGR henceforth. Taking the ratio rather than VS cancels out global physiological perfusion confounders and imaging acquisition differences and allows a comparison between the pediatric and adult group. However, the VS was used to study the associations between VS and other hemodynamic parameters within each group.
Table 1.Scan parameters of functional scans in children and adult subjects.
Figure 1.An example of the semi-automatic segmentation method of the sagital and straight sinuses. (A) Sagittal group-average image in standard space. (B) Manual segmentation of the venous signal (VS) in image (A). (C) Image in native arterial spin labeling (ASL) space of a representative patient. (D) Red overlay of the segmented VS on the ASL image.
CBF, OEF and CMRO2 were calculated as described in the previous studies10,21 and are illustrated in Figure 2 and described in the Online Supplementary Appendix, together with the analysis of the velocity in the brain feeding arteries. White matter lesions were delineated on FLAIR images using manual segmentation and quantified as previously described.6,21
Figure 2.Visual representation of the image analysis steps. Left: the arterial spin labeling (ASL) signal was obtained by subtracting the control images from the labeled images. Cerebral blood flow (CBF) was then quantified using the dual compartment model and the obtained CBF image was multiplied by the gray matter (GM) mask (created from anatomical images not shown here) to obtain the GM CBF image. CBF images were also used for venous signal assessment. Right: T2-prepared tissue relaxation with an inversion recovery magnetic resonance imaging (T2-TRIR) images were used to obtain the T1 and T2 values of venous blood. An automatic localizer tool was used to detect the region of interest in the posterior part of the brain (in red). This tool searched for high-intensity signals in the last 7 phases of the magnitude reconstructed T2-TRIR data to detect the sagittal sinus. Subsequently, a fit was initiated in the sagittal sinus voxels (shown in green on the mask image). The fitted T1 values were used for CBF quantification in adults and the fitted T2 values were converted to venous oxygen saturation using the hemoglobin A (HbA) and HbS calibration models. Subsequently, the venous oxygen saturation and arterial oxygen saturation (assumed to be 0.98) were used to calculate the oxygen extraction fraction (OEF). CMRO2 was obtained by multiplying the OEF by the CBF and oxygen-carrying capacity (CaO2). For details, see the Online Supplementary Appendix. ATT: arterial transit time.
Statistical analyses
Statistical analyses were performed in SPSS v26 (IBM, NY, USA). P<0.05 was considered statistically significant. Independent t-tests and paired-sample t-tests (or non-parametric alternatives in case of non-normality) were performed to test differences between the groups and to test the statistical differences of VGR before and after ACZ administration. We performed the correlation analysis to assess the relationship between VS and CBF. In order to test if VS contained additional hemodynamic information independently from CBF, the residuals of the regression between VS and CBF (VSCBF) were used as the dependent variable in stepwise multiple linear regression analyses against age, sex, participant groups, hydroxyl-urea, hemoglobin (Hb; as a measure of anemia) and lactate dehydrogenase (LDH) (as a measure of hemolysis) for the baseline data. In order to test the role of oxygen metabolism and ACZ on VS in adult participants, linear mixed-effects modeling was performed, in which OEF, CMRO2 and ACZ condition were added as additional predictors.
In order to further explore the relationship between parameters of oxygen metabolism and hemodynamic and laboratory parameters in adult SCD patients, we performed linear mixed-effects modeling with CMRO2, corrected for VSCBF as the dependent variable. The predictors tested in the model were CBF, VS, Hb, LDH, HbF, and ACZ condition.
Table 2.Demographic, clinical and imaging summary of participants.
Results
Demographic, clinical and imaging parameters of 28 pediatric patients with SCD, 38 adult patients and ten healthy controls are presented in Table 2. The Online Supplementary Table S1 shows the available sample size for each included parameter. Of this sample, nine (32%) pediatric and 15 (39%) adult patients were using hydroxyurea. In the adult SCD group, three (8%) patients received regular blood exchange transfusions and were studied 3-28 days since their last transfusion.
Venous arterial spin labeling signal
Venous ASL signal was observed in 27 (96%) pediatric patients, 36 (95%) adult patients and one (10%) healthy control. After ACZ administration, VS was observed in all adult participants (Figure 3A). VGR was significantly higher in patients compared to controls (Z=-4.24, P<0.01), but not between pediatric and adult SCD patients (t=-1.78, P=0.08) (Figure 3B). After ACZ administration, VGR increased in adult patients (t=-6.10, P<0.01) and controls (Z=2.80, P<0.01) compared to baseline VGR (Figure 3B). No significant VGR differences were found between patients receiving hydroxyurea and those not receiving this intervention.
Associations of venous signal with other parameters
We observed a positive relationship between baseline CBF and baseline VS (R2=0.59; F(1,73)=104.4, P<0.001; Figure 4A). After ACZ administration, these associations remained significant (R2=0.57; F(1,45)=59.4, P<0.001; Figure 4B). No significant association was found between VS and lesion volume in patients with SCD (β=-0.002 P=0.89). Subsequently, we tested associations with baseline VSCBF in all participants. Stepwise multiple regression analysis demonstrated significant associations with Hb (β=-0.26, P<0.001) and participant group (β=0.55, P=0.004; total model R2=0.23). Subsequent analyses demonstrated a negative correlation between Hb and VSCBF in the combined patient groups (β=-0.45, P<0.001), but not in the control group (β=0.02, P=0.85) (Figure 5A). After splitting the patient groups, the association between Hb and VSCBF remained significant in both adults (β=-0.43, P<0.001) and children (β=-0.45, P=0.02) with SCD.
For adult participants, we added OEF, CMRO2 and ACZ condition as additional parameters and used a linear mixed model to accommodate the repeated measures dependencies. The strongest predictor of VSCBF was CMRO2 (β=-0.79, F(1,81)=-24.5, P<0.001), demonstrating that subjects with lower CMRO2 showed higher VS, independent of CBF (Figure 5B). Subsequently, participant group (adult SCD vs. controls) in combination with either Hb (model 1) or LDH (model 2) were significant additional predictors (model 1: CMRO2 β=-0.75, F(1,82)=19.8, P<0.001 / group β=-69.6, F(1,72)=11.5, P=0.001 / Hb β=-11.6, F(1,78)=11.7, P=0.001; model 2: CMRO2 β=-0.92, F(1,78)=27.8, P<0.001 / group β=-40.1, F(1,78)=5.4, P=0.023 / LDH β=0.11, F(1,79)=10.0, P=0.002). When splitting the group into SCD adults and controls, we observed that in addition to CMRO2, Hb and LDH were significant independent predictors of VSCBF in patients (but not in controls) with Hb being a negative predictor (β=-14.2, P<0.001) and LDH being a positive predictor (β=0.13, P=0.002). We repeated our statistical analyses excluding the three patients receiving exchange transfusion and found comparable results.
Figure 3.Venous signal across groups. (A) Group average axial and sagittal arterial spin labeling (ASL) images in Montreal Neurological Institute (MNI) space before and after acetazolamide (ACZ) administration. In healthy controls, venous signal (VS) is minimal at baseline (pre-ACZ) but appears after ACZ administration. (B) Dot plots and the mean with the standard deviation of VS to gray matter signal ratio (VGR) between the groups.
In order to assess further contributing factors to CMRO2 in SCD patients, we explored the relationships between CMRO2 and the various hemodynamic and hematological markers. Linear mixed modeling demonstrated that HbF and LDH were significant predictors of CMRO2 when corrected for VSCBF. Higher HbF was associated with lower CMRO2 (β=-1.4, F(1,56)=23.4, P<0.001) and higher LDH was associated with higher CMRO2 (β=0.05, F(1,61)=8.5, P=0.005). As HbF is increased by hydroxyurea, we further split the data into hydroxyurea and no hydroxyurea groups to investigate its effects. The associations between HbF and CMRO2, and LDH and CMRO2 were similar across both groups.
Discussion
The purpose of this study was to investigate cerebral functional shunting in patients with SCD, by means of exploring the relationship between the VS intensity in ASL images and cerebral circulation and oxygen metabolism. Our results confirmed prior observations of higher VS in patients with SCD compared to healthy controls, although no differences between pediatric and adult patients were observed. Our VS data showed a strong association with CBF, confirming that venous outflow is indeed strongly dependent on the inflow to the brain. A key finding of this study was that the CBF-corrected VS was negatively associated with CMRO2, demonstrating that VS contains important information about microvascular oxygen offloading, providing evidence for functional shunting. Additionally, higher disease severity, as reflected by more severe anemia and hemolysis was, independently of CBF, associated with higher VS. Finally, in addition to potential cerebral shunting mechanisms in SCD patients, we found that CMRO2 was further negatively correlated with HbF and positively correlated with LDH.
The oxygen-carrying capacity of blood in patients with SCD is reduced and as a result, compensatory increases of CBF are observed.7– 9 Previous studies using both continuous and categorical estimates of VS13,28,29 observed that higher perfusion in SCD patients results in the presence of venous ASL signal at standard post-label delays, and the results of the current study corroborate these findings. In addition to CBF, microvascular permeability affects the fraction of labels arriving in the venous sinuses. In case of endothelial dysfunction, the permeability-surface area product (PS), defined as the flow of molecules through the capillary membranes in a certain volume of tissue, may be altered. At physiologically plausible CBF levels and in subjects with intact microvasculature, labeled spins enter the tissue and barely any VS is observed (Figure 6A). In case of increased CBF (e.g., after ACZ administration) and intact PS, VS can be observed in the brain (Figure 6B). However, at elevated CBF with low PS, we expect more labeled spins to pass unexchanged into the cerebral venous circulation (Figure 6C), a process we refer to as microvascular shunting. Importantly, our current findings provide empirical evidence for shunting, by showing that higher VS (independently of CBF) is associated with lower CMRO2. Our results are concordant with decreased OEF values observed in SCD patients with categorically increased venous signal scores.20 Taken together, these data suggest that the venous signal, beyond that expected for a given CBF, is a biomarker of cerebral functional shunting, which might be a result of PS reduction.
The inverse relation between CBF and the parameters OEF and CMRO2 in chronic anemia might seem paradoxical, but can be explained by physiological mechanisms at the capillary level. Low Hb levels trigger reciprocal increases in CBF to preserve resting oxygen delivery. However, elevated CBF results in decreased capillary transit times, potentially promoting trans-capillary pressure, and increasing heterogeneity in transit times across capillaries in the vascular bed.30,31 As oxygen extraction is dependent on the time the blood resides in the capillary network, elevated CBF can result in lower OEF and CMRO2. However, it is important to consider our findings with respect to the existing body of literature. Studies using T2-relaxation-under-spin-tagging (TRUST) have shown increased, unaffected or decreased OEF values in patients with SCD compared to controls,13,19,20,32,33 depending on the calibration model used for the SCD patients. Inconsistent OEF values were also found using other techniques, with a previous PET study showing no differences in OEF and CMRO2 between patients and controls,7 whereas studies applying susceptibility MRI techniques reported both increased or decreased OEF values in SCD patients.34–36 Our findings cannot be fully attributed to microvascular shunting per se, as we did not assess shunting locally at the capillary level, but rather used VS as a proxy. Nevertheless, although anatomical shunts in this population have been reported,37,38 microvascular shunting seems the most likely explanation, given the relation between VS and CMRO2 in the current study, and a previous report showing a reduction of VS after blood transfusion.39 Interestingly, a very recent study using multi-TI ASL showed that despite lower bolus arrival times in the tissue, much longer arrival times were found in the sagittal sinus of patients with SCD compared to controls, suggesting that hyper-perfusion in the arterial tree may be accompanied by altered capillary and/or venular micro-circulatory flow patterns, potentially due to microvascular resistance.29 Future studies are needed to further elucidate these processes in SCD at the capillary level.
Figure 4.Association between venous signal and cerebral blood flow. Scatterplots of massociations between venous signal (VS) and gray matter cerebral blood flow (CBF) before (left) and after (right) acetazolamide (ACZ). SCD: sickle cell disease administration.
Increased VS beyond that expected for a given CBF (VSCBF) might represent a reduction in PS that may be pathological. We demonstrated an inverse relationship of VSCBF with Hb in both pediatric and adult patients. In adults, VSCBF was also highly correlated with LDH. This may indicate that SCD patients with higher disease severity, i.e., more anemia and hemolysis, have lower functional cerebral microvascular surface area. Further studies are needed to confirm this hypothesis and to investigate if such changes represent permanent damage or a reversible physiologic phenotype.
We demonstrated that besides VS, CMRO2 was also influenced by HbF and LDH. Whereas HbF reduces HbS polymerization, which diminishes the complications of SCD, we found a negative association of HbF with CMRO2. This negative association is likely a result of the higher oxygen-binding affinity of HbF compared to HbA,40 potentially resulting in lower oxygen extraction and CMRO2 for the same CBF. Notably, the HbF effect was not driven specifically by patients on hydroxyurea. The positive association between LDH and CMRO2 most likely reflects the impact of dense red blood cells (DRBC) on the Hb-dissociation curve. DRBC are an important biomarker of SCD severity because of their predisposition for polymerization and hemolysis. Dense cells have markedly low oxygen affinity41 which would facilitate oxygen extraction. Indeed, another study demonstrated that, on multivariate analysis, LDH was the strongest predictor of DRBC, with biliru-bin and HbF levels also retained in the model.42
Figure 5.Predictors of VSCBF. Left: scatterplot of hemoglobin (Hb) and VSCBF (the residual of the regression between venous signal [VS] and cerebral blood flow [CBF]) for all pre-acetazolamide (pre-ACZ) data of all subject groups. The association between Hb and VSCBF was only significant in the patients, but not in controls. Right: scatterplot of cerebral metabolic rate of oxygen (CMRO2) and VSCBF for the pre-ACZ and post-ACZ data in all adult subjects, showing a significant correlation across all subjects. Please note that the values for VSCBF differ for both regressions as the pre-ACZ data were analyzed using multiple linear regression, whereas the pre-ACZ vs. post-ACZ data were analyzed using linear mixed models to accommodate the repeated measures variation.
VSCBF was comparable between pediatric and adult SCD patients, despite differences in absolute CBF. On the one hand, one might postulate higher VSCBF in children because of their increased cerebral metabolic rates and greater risk for ischemic injury.43 Alternatively, one might expect progressive microvascular damage and loss of cer-ebrovascular cross-sectional area with increasing age, similar to microvascular disease in other organs in SCD patients. A preliminary report assessing blood-brain barrier permeability demonstrated that children with SCD have higher PS than controls, despite lower extraction, suggesting increased functional capillary exchange area.44 However, this remains to be confirmed in future studies that are encouraged to incorporate OEF and CMRO2 measurements from pediatric patients as well as age-matched healthy controls. This is important because SCI are also present in pediatric patients with SCD and it is, therefore, important to understand potential age-dependent underlying mechanisms.4,45
We were unable to demonstrate any link between VSCBF and white matter damage, however, the lesion load was relatively low, and no patients with overt stroke were included. Previous studies incorporating these samples also showed no association between lesion volume and hemodynamic parameters.10,15 Although evidence from neurodegenerative disorders and vascular disease suggest that there is a link between hypoxia and white matter lesions,46,47 it remains to be determined how lower OEF and CMRO2, and potentially functional shunting mechanisms are related to such lesions in patients with chronic anemia. A previous study20 demonstrated a relationship between VS and OEF only in patients with a significant lesion load, but not in patients with overt stroke. However, they also did not correlate VS to lesion volumes.
In this study, T2-TRIR instead of TRUST was used to obtain venous oxygenation estimates for OEF calculation. T2-TRIR enables T1 measurements in addition to T2 measurements. Given this advantage and a previous study showing OEF values in the same range as the values obtained from TRUST,22 we opted for T2-TRIR here. However, a recent study from our group48 demonstrated that the OEF values obtained from T2-TRIR were significantly lower compared to OEF values obtained from the TRUST sequence (despite comparable reproducibility), which may explain our lower OEF and CMRO2 values compared to previous studies using the TRUST sequence.20,33 Additionally, we here used the HbS calibration model introduced by Bush et al.19 to obtain the venous oxygenation in patients with SCD (although see the Online Supplementary Appendix, showing comparable findings using the HbA model for completeness). Recently, a novel sickle cell specific calibration model, the Bush-Li model,49 was introduced, based on a larger range of hematocrit. However, they did not show statistical superiority compared to the currently used Bush model, and therefore, we do not expect that our conclusions would be different using the Bush-Li model. An important strength of this study is the fact that we used a continuous estimate of VS, which allowed us to conduct regression analyses to investigate the relation between VS, hemodynamic measures and hematological parameters. However, our sample size is relatively small for multiple regression analyses, which is why we have carefully chosen our included predictors. Future studies should confirm our regression analyses using larger sample sizes, to further explore the presence of functional shunting and the associations between shunting and the parameters of oxygen metabolism. Furthermore, our pediatric sample did not include age-matched healthy controls. However, previous studies have already shown that CBF is higher in children with SCD compared to healthy children.9,15 In addition, Wu et al.28 showed that the signal in the sagittal sinus was lower in healthy children compared to children with SCD. Furthermore, the MRI scanner, the head coil and multiple scan parameters differed between the pediatric and adult patients, which precluded direct comparison of CBF and lesion volume estimates. Individual M0 scans were not acquired in our pediatric cohort, because the scanning was done prior to the M0 consensus recommendations in 2014.24 Instead, we used a single fixed M0 value to quantify CBF. In a post hoc analysis in adults, we found comparable coefficients of variation for GM CBF obtained from a group average M0 and from a subject-specific M0. Moreover, using a group average M0 did not change the associations between VS and the parameter GM CBF and the markers Hb and LDH. Nevertheless, using a single M0 value does not account for slight spatial differences in magnetization, which is why a subject-specific M0 was used in our adult participants and why this is also recommended for future studies. Additionally, the T2-TRIR scan was not available for the pediatric study, and therefore OEF and CMRO2 estimates could not be obtained in pediatric patients. In the adults, only a small proportion of patients were on chronic exchange transfusion, and therefore its effect could not be estimated in this study. Although hydroxyurea did not appear to affect the relationships found in this study, future studies in larger samples could shed light on the effects of hydroxyurea treatment on functional shunting. Our control group included two sickle cell trait participants whose hematological parameters and T2 values were in the same range as those of the HbAA controls, in line with previous studies.19,50 This suggests that they were appropriately included as control subjects.
Figure 6.Schematic hypothetical outline of venous signal in the absence and presence of shunting. Arterially labeled spins enter brain parenchyma at a rate determined by the product of vascular surface area and intrinsic water permeability (permeability-surface area product). The re-entry of the label was previously shown to be negligible (Lin et al. 2018) and is therefore not displayed. (A) In the brain of a healthy control with normal cerebral blood flow (CBF) levels, (almost) no venous signal (VS) is observed. (B) In the brain of a healthy control after acetazolamide administration with increased CBF and intact permeability-surface area product (PS), VS signal can be observed as flow predominates the microvascular uptake. (C) In a patient with sickle cell disease (SCD) with the presence of shunting, blood preferentially passes through short, low resistance capillary pathways leading to functional loss of capillary surface area, further increasing VS and impairing oxygen unloading. HC: healthy control, ACZ: acetazolamide.
In summary, our findings suggest that higher CBF-corrected venous signal may reflect the loss of capillary exchange area as shown by its relationship to CMRO2, providing evidence for cerebral functional shunting. Moreover, we found that higher disease severity is related to higher VS and that CMRO2 is additionally influenced by HbF and LDH. These findings indicate that the venous signal on ASL images can be considered as a complementary biomarker of cerebral perfusion and oxygen metabolism.
Footnotes
- Received October 19, 2021
- Accepted May 5, 2022
Correspondence
Disclosures
BB has financial relationship with GBT, Novartis, Sanquin, Novo nordisk, Celgene, CSL Behring, CHIESI and Bluebird Bio. LV has a financial relationship with Philips Healthcare
Contributions
LAH was involved in the study conception and design, analysis and interpretation of the data and manuscript drafting; LV and BJB were involved in data acquisition, analysis and interpretation of the data and manuscript revision; JCW was involved in analysis and interpretation of the data and manuscript revision; AJN and HJMMM were involved in study conception and design, analysis and interpretation of the data and manuscript revision; AS was involved in analysis and interpretation of the data and manuscript drafting..
Data-sharing statement
Data available on request from the authors.
Funding
The authors disclose receipt of the following financial support for the research, authorship, and/or publication of this article: National Heart Lung and Blood Institute (1RO1HL136484-A1).
Acknowledgements
The authors would like to thank all the participants for their time and efort to take part in this study and they are grateful to Dr. Jan Petr for his advice regarding ExploreASL.
References
- Rees DC, Williams TN, Gladwin MT. Sickle-cell disease. Lancet. 2010; 376(9757):2018-2031. https://doi.org/10.1016/S0140-6736(10)61029-X|PubMed|Google Scholar
- Debaun MR, Kirkham FJ. Central nervous system complications and management in sickle cell disease. Blood. 2016; 127(7):829-838. https://doi.org/10.1182/blood-2015-09-618579|PubMed|Google Scholar
- Strouse JJ, Jordan LC, Lanzkron S. The excess burden of stroke in hospitalized adults with sickle cell disease. Am J Hematol. 2009; 84(9):548-552. https://doi.org/10.1002/ajh.21476|PubMed|PubMed Central|Google Scholar
- Ford AL, Ragan DK, Fellah S. Silent infarcts in sickle cell disease occur in the border zone region and are associated with low cerebral blood flow. Blood. 2018; 132(16):1714-1723. https://doi.org/10.1182/blood-2018-04-841247|PubMed|PubMed Central|Google Scholar
- Bernaudin F, Verlhac S, Arnaud C. Chronic acute anemia and extracranial internal carotid stenosis are risk factors for silent cerebral infarcts in sickle cell anemia. Blood. 2015; 125(10):1653-1661. https://doi.org/10.1182/blood-2014-09-599852|PubMed|Google Scholar
- van der Land V, Mutsaerts HJMM, Engelen M. Risk factor analysis of cerebral white matter hyperintensities in children with sickle cell disease. Br J Haematol. 2016; 172(2):274-284. https://doi.org/10.1111/bjh.13819|PubMed|Google Scholar
- Herold S, Brozovic M, Gibbs J. Measurement of regional cerebral blood flow, blood volume and oxygen metabolism in patients with sickle cell disease using positron emission tomography. Stroke. 1986; 17(4):692-698. https://doi.org/10.1161/01.STR.17.4.692|PubMed|Google Scholar
- Numaguchi Y, Haller JS, Humbert JR. Cerebral blood flow mapping using stable Xenon-enhanced CT in sickle cell cerebrovascular disease. Neuroradiology. 1990; 32(4):289-295. https://doi.org/10.1007/BF00593047|PubMed|Google Scholar
- Oguz KK, Golay X, Pizzini FB. Sickle cell disease: continuous arterial spin-labeling perfusion MR imaging in children. Radiology. 2003; 227(2):567-574. https://doi.org/10.1148/radiol.2272020903|PubMed|Google Scholar
- Václavů L, Petr J, Petersen ET. Cerebral oxygen metabolism in adults with sickle cell disease. Am J Hematol. 2020; 95(4):401-412. https://doi.org/10.1002/ajh.25727|PubMed|PubMed Central|Google Scholar
- Kassim AA, Pruthi S, Day M. Silent cerebral infarcts and cerebral aneurysms are prevalent in adults with sickle cell anemia. Blood. 2016; 127(16):2038-2040. https://doi.org/10.1182/blood-2016-01-694562|PubMed|Google Scholar
- Bush A, Chai Y, Choi SY. Pseudo continuous arterial spin labeling quantification in anemic subjects with hyperemic cerebral blood flow. Magn Reson Imaging. 2018; 47:137-146. https://doi.org/10.1016/j.mri.2017.12.011|PubMed|PubMed Central|Google Scholar
- Juttukonda MR, Donahue MJ, Davis LT. Preliminary evidence for cerebral capillary shunting in adults with sickle cell anemia. J Cereb Blood Flow Metab. 2019; 39(6):1099-1110. https://doi.org/10.1177/0271678X17746808|PubMed|PubMed Central|Google Scholar
- Van Den Tweel XW, Nederveen AJ, Majoie CBLM. Cerebral blood flow measurement in children with sickle cell disease using continuous arterial spin labeling at 3.0-tesla MRI. Stroke. 2009; 40(3):795-800. https://doi.org/10.1161/STROKEAHA.108.523308|PubMed|Google Scholar
- Gevers S, Nederveen AJ, Fijnvandraat K. Arterial spin labeling measurement of cerebral perfusion in children with sickle cell disease. J Magn Reson Imaging. 2012; 35(4):779-787. https://doi.org/10.1002/jmri.23505|PubMed|Google Scholar
- Lin Z, Li Y, Su P. Non-contrast MR imaging of blood-brain barrier permeability to water. Magn Reson Med. 2018; 80(4):1507-1520. https://doi.org/10.1002/mrm.27141|PubMed|PubMed Central|Google Scholar
- Renkin EM. Transport of potassium-42 from blood to tissue in isolated mammalian skeletal muscles. Am J Physiol. 1959; 197:1205-1210. https://doi.org/10.1152/ajplegacy.1959.197.6.1205|PubMed|Google Scholar
- Crone C. The permeability of capillaries in various organs as determined by use of the ‘indicator diffusion’ method. Acta Physiol Scand. 1963; 58:292-305. https://doi.org/10.1111/j.1748-1716.1963.tb02652.x|PubMed|Google Scholar
- Bush AM, Coates TD, Wood JC. Diminished cerebral oxygen extraction and metabolic rate in sickle cell disease using T2 relaxation under spin tagging MRI. Magn Reson Med. 2018; 80(1):294-303. https://doi.org/10.1002/mrm.27015|PubMed|PubMed Central|Google Scholar
- 20.Juttukonda MR, Donahue MJ, Waddle SL. Reduced oxygen extraction efficiency in sickle cell anemia patients with evidence of cerebral capillary shunting. J Cereb Blood Flow Metab. 2020; 41(3):546-560. https://doi.org/10.1177/0271678X20913123|PubMed|PubMed Central|Google Scholar
- Václavů L, Meynart BN, Mutsaerts HJMM. Hemodynamic provocation with acetazolamide shows impaired cerebrovascular reserve in adults with sickle cell disease. Haematologica. 2019; 104(4):690-699. https://doi.org/10.3324/haematol.2018.206094|PubMed|PubMed Central|Google Scholar
- De Vis JB, Petersen ET, Alderliesten T. Non-invasive MRI measurements of venous oxygenation, oxygen extraction fraction and oxygen consumption in neonates. Neuroimage. 2014; 95:185-192. https://doi.org/10.1016/j.neuroimage.2014.03.060|PubMed|Google Scholar
- Vaclavu L, Van Der Land V, Heijtel DFR. In vivo T1 of blood measurements in children with sickle cell disease improve cerebral blood flow quantification from arterial spin-labeling MRI. Am J Neuroradiol. 2016; 37(9):1727-1732. https://doi.org/10.3174/ajnr.A4793|PubMed|PubMed Central|Google Scholar
- Alsop DC, Detre JA, Golay X. 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. Magn Reson Med. 2015; 73(1):102-116. https://doi.org/10.1002/mrm.25197|PubMed|PubMed Central|Google Scholar
- Zhao MY, Václavů L, Petersen ET. Quantification of cerebral perfusion and cerebrovascular reserve using Turbo-UASAR arterial spin labeling MRI. Magn Reson Med. 2020; 83(2):731-748. https://doi.org/10.1002/mrm.27956|PubMed|PubMed Central|Google Scholar
- Mutsaerts H, Petr J, Groot P. ExploreASL: an image processing pipeline for multi-center ASL perfusion MRI studies. Neuroimage. 2020; 219:117031. https://doi.org/10.1016/j.neuroimage.2020.117031|PubMed|Google Scholar
- Evans AC, Janke AL, Collins DL. Brain templates and atlases. Neuroimage. 2012; 62(2):911-922. https://doi.org/10.1016/j.neuroimage.2012.01.024|PubMed|Google Scholar
- Wu WC, St Lawrence KS, Licht DJ. Quantification issues in arterial spin labeling perfusion magnetic resonance imaging. Top Magn Reson Imaging. 2010; 21(2):65-73. https://doi.org/10.1097/RMR.0b013e31821e570a|PubMed|Google Scholar
- Stotesbury H, Hales PW, Koelbel M. Venous cerebral blood flow quantification and cognition in patients with sickle cell anemia. J Cereb Blood Flow Metab. 2022; 42(6):1061-1077. https://doi.org/10.1177/0271678X211072391|PubMed|PubMed Central|Google Scholar
- Jespersen SN, Østergaard L. The roles of cerebral blood flow, capillary transit time heterogeneity, and oxygen tension in brain oxygenation and metabolism. J Cereb Blood Flow Metab. 2012; 32(2):264-277. https://doi.org/10.1038/jcbfm.2011.153|PubMed|PubMed Central|Google Scholar
- Sakadžić S, Mandeville ET, Gagnon L. Large arteriolar component of oxygen delivery implies a safe margin of oxygen supply to cerebral tissue. Nat Commun. 2014; 5:5734. https://doi.org/10.1038/ncomms6734|PubMed|PubMed Central|Google Scholar
- Jordan LC, Gindville MC, Scott AO. Non-invasive imaging of oxygen extraction fraction in adults with sickle cell anaemia. Brain. 2016; 139(Pt 3):738-750. https://doi.org/10.1093/brain/awv397|PubMed|PubMed Central|Google Scholar
- Vu C, Bush A, Choi S. Reduced global cerebral oxygen metabolic rate in sickle cell disease and chronic anemias. Am J Hematol. 2021; 96(8):901-913. https://doi.org/10.1002/ajh.26203|PubMed|PubMed Central|Google Scholar
- Guilliams KP, Fields ME, Ragan DK. Red cell exchange transfusions lower cerebral blood flow and oxygen extraction fraction in pediatric sickle cell anemia. Blood. 2018; 131(9):1012-1021. https://doi.org/10.1182/blood-2017-06-789842|PubMed|PubMed Central|Google Scholar
- Croal PL, Leung J, Phillips CL. Quantification of pathophysiological alterations in venous oxygen saturation: A comparison of global MR susceptometry techniques. Magn Reson Imaging. 2019; 58:18-23. https://doi.org/10.1016/j.mri.2019.01.008|PubMed|Google Scholar
- Wang Y, Fellah S, Fields ME. Cerebral oxygen metabolic stress, microstructural injury, and infarction in adults with sickle cell disease. Neurology. 2021; 97(9):e902-e912. https://doi.org/10.1212/WNL.0000000000012404|PubMed|PubMed Central|Google Scholar
- Hambley BC, Rahman RA, Reback M. Intracardiac or intrapulmonary shunts were present in at least 35% of adults with homozygous sickle cell disease followed in an outpatient clinic. Haematologica. 2019; 104(1):e1-e3. https://doi.org/10.3324/haematol.2018.196212|PubMed|PubMed Central|Google Scholar
- Dowling MM, Quinn CT, Ramaciotti C. Increased prevalence of potential right-to-left shunting in children with sickle cell anaemia and stroke. Br J Haematol. 2017; 176(2):300-308. https://doi.org/10.1111/bjh.14391|PubMed|PubMed Central|Google Scholar
- DeBeer T, Jordan LC, Lee CA. Evidence of transfusion-induced reductions in cerebral capillary shunting in sickle cell disease. Am J Hematol. 2020; 95(9):E228-E230. https://doi.org/10.1002/ajh.25863|PubMed|PubMed Central|Google Scholar
- Young RC, Rachal RE, Del Pilar Aguinaga M. Automated oxyhemoglobin dissociation curve construction to assess sickle cell anemia therapy. J Natl Med Assoc. 2000; 92(9):430-435. Google Scholar
- Di Liberto G, Kiger L, Marden MC. Dense red blood cell and oxygen desaturation in sickle-cell disease. Am J Hematol. 2016; 91(10):1008-1013. https://doi.org/10.1002/ajh.24467|PubMed|Google Scholar
- Bartolucci P, Brugnara C, Teixeira-Pinto A. Erythrocyte density in sickle cell syndromes is associated with specific clinical manifestations and hemolysis. Blood. 2012; 120(15):3136-3141. https://doi.org/10.1182/blood-2012-04-424184|PubMed|Google Scholar
- DeBaun MR, Sarnaik SA, Rodeghier MJ. Associated risk factors for silent cerebral infarcts in sickle cell anemia: Low baseline hemoglobin, sex, and relative high systolic blood pressure. Blood. 2012; 119(16):3684-3690. https://doi.org/10.1182/blood-2011-05-349621|PubMed|PubMed Central|Google Scholar
- Lin Z, Lance E, Li Y. Impaired blood-brain barrier function in pediatric sickle cell disease.Google Scholar
- DeBaun MR, Armstrong FD, McKinstry RC. Silent cerebral infarcts: a review on a prevalent and progressive cause of neurologic injury in sickle cell anemia. Blood. 2012; 119(20):4587-4596. https://doi.org/10.1182/blood-2011-02-272682|PubMed|PubMed Central|Google Scholar
- Yatawara C, Lee D, Ng KP. Mechanisms linking white matter lesions, tract integrity, and depression in alzheimer disease. Am J Geriatr Psychiatry. 2019; 27(9):948-959. https://doi.org/10.1016/j.jagp.2019.04.004|PubMed|Google Scholar
- Van Dijk EJ, Prins ND, Vrooman HA. Progression of cerebral small vessel disease in relation to risk factors and cognitive consequences: Rotterdam scan study. Stroke. 2008; 39(10):2712-2719. https://doi.org/10.1161/STROKEAHA.107.513176|PubMed|Google Scholar
- Baas K, Coolen B, Petersen E. Comparative analysis of blood T2 values measured by T2-RIR and TRUST. J Magn Reson Imaging. 2022; 56(2):516-526. https://doi.org/10.1002/jmri.28066|PubMed|Google Scholar
- Bush A, Vu C, Choi S. Calibration of T2 oximetry MRI for subjects with sickle cell disease. Magn Reson Med. 2021; 86(2):1019-1028. https://doi.org/10.1002/mrm.28757|PubMed|PubMed Central|Google Scholar
- Eaton W, Hofrichter J. Hemoglobin S gelation and sickle cell disease. Blood. 1987; 70(5):1245-1266. https://doi.org/10.1182/blood.V70.5.1245.bloodjournal7051245|Google Scholar
Data Supplements
Figures & Tables
Table 1.Scan parameters of functional scans in children and adult subjects.
Figure 1.An example of the semi-automatic segmentation method of the sagital and straight sinuses. (A) Sagittal group-average image in standard space. (B) Manual segmentation of the venous signal (VS) in image (A). (C) Image in native arterial spin labeling (ASL) space of a representative patient. (D) Red overlay of the segmented VS on the ASL image.
Figure 2.Visual representation of the image analysis steps. Left: the arterial spin labeling (ASL) signal was obtained by subtracting the control images from the labeled images. Cerebral blood flow (CBF) was then quantified using the dual compartment model and the obtained CBF image was multiplied by the gray matter (GM) mask (created from anatomical images not shown here) to obtain the GM CBF image. CBF images were also used for venous signal assessment. Right: T2-prepared tissue relaxation with an inversion recovery magnetic resonance imaging (T2-TRIR) images were used to obtain the T1 and T2 values of venous blood. An automatic localizer tool was used to detect the region of interest in the posterior part of the brain (in red). This tool searched for high-intensity signals in the last 7 phases of the magnitude reconstructed T2-TRIR data to detect the sagittal sinus. Subsequently, a fit was initiated in the sagittal sinus voxels (shown in green on the mask image). The fitted T1 values were used for CBF quantification in adults and the fitted T2 values were converted to venous oxygen saturation using the hemoglobin A (HbA) and HbS calibration models. Subsequently, the venous oxygen saturation and arterial oxygen saturation (assumed to be 0.98) were used to calculate the oxygen extraction fraction (OEF). CMRO2 was obtained by multiplying the OEF by the CBF and oxygen-carrying capacity (CaO2). For details, see the Online Supplementary Appendix. ATT: arterial transit time.
Table 2.Demographic, clinical and imaging summary of participants.
Figure 3.Venous signal across groups. (A) Group average axial and sagittal arterial spin labeling (ASL) images in Montreal Neurological Institute (MNI) space before and after acetazolamide (ACZ) administration. In healthy controls, venous signal (VS) is minimal at baseline (pre-ACZ) but appears after ACZ administration. (B) Dot plots and the mean with the standard deviation of VS to gray matter signal ratio (VGR) between the groups.
Figure 4.Association between venous signal and cerebral blood flow. Scatterplots of massociations between venous signal (VS) and gray matter cerebral blood flow (CBF) before (left) and after (right) acetazolamide (ACZ). SCD: sickle cell disease administration.
Figure 5.Predictors of VSCBF. Left: scatterplot of hemoglobin (Hb) and VSCBF (the residual of the regression between venous signal [VS] and cerebral blood flow [CBF]) for all pre-acetazolamide (pre-ACZ) data of all subject groups. The association between Hb and VSCBF was only significant in the patients, but not in controls. Right: scatterplot of cerebral metabolic rate of oxygen (CMRO2) and VSCBF for the pre-ACZ and post-ACZ data in all adult subjects, showing a significant correlation across all subjects. Please note that the values for VSCBF differ for both regressions as the pre-ACZ data were analyzed using multiple linear regression, whereas the pre-ACZ vs. post-ACZ data were analyzed using linear mixed models to accommodate the repeated measures variation.
Figure 6.Schematic hypothetical outline of venous signal in the absence and presence of shunting. Arterially labeled spins enter brain parenchyma at a rate determined by the product of vascular surface area and intrinsic water permeability (permeability-surface area product). The re-entry of the label was previously shown to be negligible (Lin et al. 2018) and is therefore not displayed. (A) In the brain of a healthy control with normal cerebral blood flow (CBF) levels, (almost) no venous signal (VS) is observed. (B) In the brain of a healthy control after acetazolamide administration with increased CBF and intact permeability-surface area product (PS), VS signal can be observed as flow predominates the microvascular uptake. (C) In a patient with sickle cell disease (SCD) with the presence of shunting, blood preferentially passes through short, low resistance capillary pathways leading to functional loss of capillary surface area, further increasing VS and impairing oxygen unloading. HC: healthy control, ACZ: acetazolamide.
Article Information

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.