Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 20 February 2023
Sec. Vascular Physiology
This article is part of the Research Topic Measuring Resting Cerebral Perfusion using Magnetic Resonance Imaging (MRI) View all 14 articles

Imputation models and error analysis for phase contrast MR cerebral blood flow measurements

  • 1Division of Cardiology, Department of Pediatrics, Children’s Hospital Los Angeles, Los Angeles, CA, United States
  • 2Department of Radiology, Children’s Hospital Los Angeles, Los Angeles, CA, United States
  • 3Division of Neonatology, Department of Pediatrics, Children’s Hospital Los Angeles, Keck School of Medicine, Fetal and Neonatal Institute, University of Southern California, Los Angeles, CA, United States

Cerebral blood flow (CBF) supports brain metabolism. Diseases impair CBF, and pharmacological agents modulate CBF. Many techniques measure CBF, but phase contrast (PC) MR imaging through the four arteries supplying the brain is rapid and robust. However, technician error, patient motion, or tortuous vessels degrade quality of the measurements of the internal carotid (ICA) or vertebral (VA) arteries. We hypothesized that total CBF could be imputed from measurements in subsets of these 4 feeding vessels without excessive penalties in accuracy. We analyzed PC MR imaging from 129 patients, artificially excluded 1 or more vessels to simulate degraded imaging quality, and developed models of imputation for the missing data. Our models performed well when at least one ICA was measured, and resulted in R2 values of 0.998–0.990, normalized root mean squared error values of 0.044–0.105, and intra-class correlation coefficient of 0.982–0.935. Thus, these models were comparable or superior to the test-retest variability in CBF measured by PC MR imaging. Our imputation models allow retrospective correction for corrupted blood vessel measurements when measuring CBF and guide prospective CBF acquisitions.

Introduction

The brain requires a significant fraction of the metabolic support of the body, and it lacks the capacity to buffer disruptions in blood supply. Thus, adequate cerebral blood flow (CBF) is vital to brain health (Lassen, 1959) and CBF measurements span many areas of brain research. CBF varies markedly with brain development throughout infancy, childhood, and senescence. (Oshima et al., 2002; Liu et al., 2019; Paniukov et al., 2020). Altered CBF is associated with cardiovascular risk, (Jennings et al., 2013; King et al., 2018), small vessel disease, (Yu et al., 2020), Alzheimer’s dementia, (Leijenaar et al., 2017), white matter lesions, (Hanaoka et al., 2016), anemia, (Borzage et al., 2016), stroke risk in sickle cell disease, (Vernooij et al., 2008; Prohovnik et al., 2009), and higher risk of non-cardiovascular mortality in the elderly. (Sabayan et al., 2013). CBF is a predictive biomarker for patients with severe depression, (Leaver et al., 2019), and CBF is modulated by pharmacological agents including substances of abuse, (Chen et al., 2016), alcohol, (Christie et al., 2008), and anesthetics. (Oshima et al., 2002).

The gold standard for CBF measurements are Kety-Schmidt tracer methods for whole brain measurements and Positron Emitting Tomography (PET) for spatially resolved measurements. However, both are invasive, making them poor candidates for neurovascular screening or research applications. A non-invasive alternative CBF measurement is arterial spin labeling (ASL) magnetic resonance imaging (MRI), which uses magnetically labeled water protons in the blood as endogenous tracers to estimate brain perfusion. However, quantitation of ASL is limited by the measurement or assumption of multiple parameters, including blood-brain partition coefficient, T1 of blood, T1 of brain tissue, labeling efficiency, and arterial transit time. (Alsop et al., 2010). Estimates of each of these key parameters are provided in the literature, however experimental methods demonstrate variability for each: the blood-brain partition coefficient varies with age, underlying pathology and brain region; (Thalman et al., 2019) T1 of blood varies with red cell characteristics and hematocrit; (Lu et al., 2004) labeling efficiency is a confounding variable that depends on the blood velocity and labelling schema; (Robertson et al., 2017) and arterial transit time varies with the health, age (MacIntosh et al., 2015; Dai et al., 2017) and brain region. (Thomas et al., 2006). The reproducibility of ASL is moderate (ICC 0.74–0.78); (Yang et al., 2019) and ASL acquisitions are lengthy, with typical scan times of 5–7 min. As a result, ASL is excellent at determining relative brain perfusion in specific regions, but less suitable for measuring absolute blood flow.

Alternatively, phase contrast (PC) MRI of the internal carotid (ICA) or vertebral (VA) arteries provides rapid and robust quantification of global CBF (ml/min). The measurement yielded by PC MRI can be normalized to brain volume from brain mass estimates derived from anatomic imaging to estimate global brain perfusion in ml/100 g/min. PC MRI is reproduceable between users (ICC 0.97–0.99) and has a low coefficient of variation (4%–9%) for serial measurements. (Koerte et al., 2013; Liu et al., 2014; Sakhare et al., 2019). If phasic variation is unimportant, then accurate CBF can be collected without cardiac gating thereby greatly accelerating the acquisition. Most importantly, PC-MRI does not require modeling or parameter estimation, making it robust across pathological states.

One barrier to deploying PC MRI in clinical practice is the need to minimize contributions from partial-voluming or off-axis flow. The MR operator prevents these issues by carefully localizing the imaging plane orthonormal to the vessel being measured. The four head vessels are not parallel, thus optimal measurement of one of vessel compromises measurements from the others. To overcome this problem, the MR operator can use four PC MR images at the cost of a four-fold increase in time. However, manual optimization of each image remains dependent on operator skill. Automated methods for analyzing the anatomy of patients and planning ideal locations of four PC MR imaging slices have been proposed, but not implemented by MR scanner vendors. (Liu et al., 2014). Alternatively, placing a “best-guess” single PC MR imaging slice to measure CBF is methodologically simple and fast, but is inherently suboptimal and occasionally results in one or more vessels being too oblique for acceptable quantitation.

We postulated that flows in each vessel were sufficiently correlated with one another so that imputation could compensate for measurements corrupted by motion or obliquity, provided that one or more vessels remain measurable. To test this hypothesis, we calculated the internal correlations among the head vessels, developed a table of mathematical models to impute the flow from corrupted arterial flow measurements, and assessed the resulting error introduced by our imputations.

Materials and methods

Patient demographics

This study is a secondary analysis of existing data, which was originally approved by the Children’s Hospital Los Angeles Committee on Clinical Investigations (CCI 11–00083). Informed consent was obtained from N = 129 patients recruited between 2012 and 2017. This cohort included patients with sickle cell disease (N = 55), healthy control patients (N = 42), and patients suffering from various hemoglobinopathies (N = 32). These patients were 23.5 ± 9.7 years (range 9–61 years) old (59M, 70F). Their hematocrits were 32.9% ± 7.2% and 33% (N = 42) were on chronic transfusion. (Borzage et al., 2016).

Image acquisition

The imaging methods are reported elsewhere in detail and summarized here. (Borzage et al., 2016). We obtained all images with a 3T Philips Achieva and eight-element head coil. A magnetic resonance angiogram localized the vessels in the neck, and a PC MR imaging plane was placed approximately 1 cm above the carotid bifurcation. The angiogram was collected in the axial plane with inline reformatting into sagittal and coronal planes to facilitate orthogonal placement of the PC imaging plane. Image parameters for the PC MR examination were as follows: repetition time, 12.3 ms; echo time, 7.5 ms; field of view, 260 mm; thickness, 5 mm; signal averages, 10; acquisition matrix, 204 × 201; reconstruction matrix, 448 × 448; bandwidth, 244 Hz/pixel; and velocity encoding gradient, 200 cm/s. For this retrospective data analysis, we retained only one PC MR image per patient. Figure 1 demonstrates the coronal angiogram showing the carotid and vertebral arteries together with the magnitude and phase images.

FIGURE 1
www.frontiersin.org

FIGURE 1. Left panel: Coronal Angiogram showing the carotid and vertebral arteries. The solid-colored portion of the vessels are above the bifurcation of the common carotid arteries, and below the level of the basilar artery. The green horizontal plane shows the location of the phase contrast slice. Center panel: a magnitude image; and right panel: a phase image, both acquired at the level indicated in the left panel.

Image processing

We performed all phase-contrast image analysis using MATLAB (The MathWorks, Natick, MA). We thresholded the complex difference image to identify moving voxels (defined as greater than the mean plus two standard deviations of stationary voxels sampled from a non-vascular region). We fit the phase differences of stationary voxels using a two-dimensional second-order polynomial to remove the background phase. We identified vessel boundaries using a Canny edge-detector of the complex difference image, dilating the edge by a single voxel, and excluding any stationary voxels. We calculated the blood flow in each artery by summing the blood velocities (cm/s) within the vessel multiplied by the voxel area (cm2). When the automatic edge detection failed (<5% of the time), the carotid or vertebral artery boundaries were identified manually by an MR researcher (JCW) with 22 years of experience analyzing PC MR images.

Modeling cerebral blood flow

We identified all 16 possible scenarios wherein combinations of two, one, or zero ICAs or VAs might be able to be analyzed from an image (Table 1). We synthesized each of these 16 scenarios from each PC MR image to simulate the effects of a sub-optimal image which failed to assess the four vessels to be analyzed. We applied a standard least-squares model to model the total cerebral blood flow (CBF) as a function of the vessels able to be analyzed. In scenarios where two ICAs or two VAs were analyzed, we simplified and reduced the degrees of freedom in the models by calculating the total anterior (sum of ICAs) or posterior (sum of VAs) flow. We evaluated the quality of each model using root mean squared error, intra-class correlation coefficients (ICC), R2 statistic. The ICC (Lassen, 1959; Liu et al., 2019) is a two-way random, single measures absolute agreement between model 0 (gold standard) and models 1-8, calculated using MATLAB. We also performed Bland-Altman analyses of all imputed models compared to total CBF. We calculated the biases as the mean difference between model 0 and each other model, and the 95% limits of agreement as twice the standard deviation of the differences of individual measurements between model 0 and each other model.

TABLE 1
www.frontiersin.org

TABLE 1. Models for computing cerebral blood flow.

Results

We report the total cerebral blood flow (933.7 ± 297.9), anterior circulation (sum of ICAs, 652.6 ± 209.2), and posterior circulation (sum of VAs, 281.0 ± 106.2). We also report the flow in the individual arteries: left ICA (328.1 ± 110.5), right ICA (324.9 ± 105.7), left VA (145.9 ± 76.8), and right VA (134.9 ± 53.1), all values are reported in units of ml/min as mean ± standard deviation. (Figure 1, Left). We also report the ratio of flow in individual arteries versus total cerebral blood flow, which are: left ICA (0.352 ± 0.039), right ICA (0.348 ± 0.039), left VA (0.154 ± 0.049), and right VA (0.147 ± 0.044), all values are unitless, and reported as mean ± standard deviation. (Figure 1, Right).

We evaluated all four ratios (left ICA-CBF, right ICA-CBF, left VA-CBF, right VA-CBF) with an ANOVA, and tested all six pairwise comparisons with Tukey-Kramer HSD. Only two pairs of ratios were not different: the left ICA-CBF versus right ICA-CBF ratios (p = 0.549), and the left VA-CBF versus right VA-CBF ratios (p = 0.549). Thus, we simplified from 16 models to eight using symmetry of the left and right CBF ratios. All other pairs of ratios were different (p < 0.0001).

We created eight statistical models of CBF using standard least squares simple linear regression. We computed the performance statistics for the models, including R2value, root mean square error (RMSE), normalized RMSE and ICC to evaluate the performance and the reliability of the models (Table 1). These models are numbered from the case with all arteries present (model 0), in ascending order of increasing expected error to the model with the most error because there are no arteries present (model 8). As anticipated from model 1 to model 8, the error in imputation (RMSE) increased and the R2 decreased. All these models were statistically significant (p < 0.0001), except when there were no usable blood vessels (model 8).

We performed Bland-Altman analyses of models 1-8 compared to total CBF (model 0), demonstrated in Figure 2. Models 1-3 and 5 are all statistically unbiased (p > 0.05), models 4 and 6 had small biases (1.71% and 3.33%, respectively). Models 1-4 had narrow limits of agreement, with standard deviations of 5.43%–8.78%. Models 5-6 and had wider limits of agreement, with standard deviations of 12.31%–15.12%. Model 7, which was derived from a measurement in one VA, and it performed exceptionally poorly (Bias −9.57%, standard deviation 28.53%). Thus, model 7 had similar performance to model 8 (Bias −9.57% standard deviation 33.69%) wherein flow is assumed equal to the population mean.

FIGURE 2
www.frontiersin.org

FIGURE 2. Bland-Altman analysis of total CBF quantified by PC in four vessels (model 0) versus various models of one or more missing vessels. Blue line demonstrates bias, dashed blue line demonstrates insignificant (p > 0.05) bias, and red dash lines demonstrate 95% confidence intervals. These results suggest that it is important to successfully capture at least one ICA or both VAs to impute CBF measurements.

Discussion

Our models were suitable to impute CBF for all scenarios wherein at least one ICA was usable (Model 1–5). We do not recommend attempting to impute CBF without any ICA measurements (models 6–8) because of their poor performance. We included models 6 and 7 because they were valid correlates of total cerebral blood flow, and model 8 (mean CBF) for completeness. We recommend the reader interpret our model performances in context of the accuracy and precision needed for each specific use case for measuring CBF. In any situation in which these models are used they add error versus the ideal situation wherein all four vessels are perfectly measured in one PC MR image. Thus, if any of the models are used, the error they introduce should be considered and the modeled results should be held with appropriate uncertainty.

The best comparisons of our model performance would be versus the performance of other CBF imputation models that use PC MR data. However, because our motivation for these models was the lack of any extant model we cannot ascertain how our models perform versus other approaches. As an imperfect alternative, we provide the reader with examples of alternative studies of CBF correlations using different measurement and validation approaches as context for our model performance. The correlation of arterial spin labeling versus 15O positron emission tomography (R2 0.47 total CBF and R2 0.48 voxelwise blood flow); (Heijtel et al., 2014) and the intraclass correlation of test-retest measurements with PC MR with no repositioning of slice or patient (ICC 0.79). (Sakhare et al., 2019). provide context to judge our imputation performance versus other methods and errors commonly encountered by those who measure CBF.

Imputation error increased more when an ICA was corrupted than a VA. This was anticipated because the ICA makes a greater contribution to CBF than a VA. Thus, a single PC MR slice placement that optimizes the measurement of the two ICAs may be preferable to a single PC MR slice placement that attempts to optimize both the ICAs and VAs. The ability to accurately estimate total CBF based only on carotid measurements (one or both) is perhaps not surprising based on prior work using transcranial Doppler. While losing one carotid compromised CBF accuracy, knowledge of the posterior flow adequately compensated for this loss, reflecting the strong conservation of flow balance between the anterior and posterior circulations.

We suggest that in time-limited scanning scenarios, that the following PC MR localization be used. (A1) Acquire a low resolution, highly accelerated MR angiography scout image with sagittal and coronal maximum intensity projection (MIP) images; this can typically be performed in approximately 30 s (A2) Use the sagittal MRA image to optimize the anterior-posterior pitch of the PC MR plane for the ICAs. (A3) Use the coronal image to optimize the left-right roll of the PC MR plane for the ICAs, (A4) ensure the PC MR slice is above the external carotid artery bifurcation, and (A5) below the junction of the VAs into the basilar artery. Alternatively, if there is not adequate time to acquire the MR angiography image, it is also quite possible to use the ubiquitous initial anatomical survey image; (B1) place the PC MR imaging plane at the level of the C2 cervical vertebra, and (B2) angle it to be perpendicular to the spinal cord. Either approach will intersect all the great vessels of the neck and therefore create the images suitable for our imputation models.

FIGURE 3
www.frontiersin.org

FIGURE 3. Left panel: measurements of flow in units of ml per minute. From left to right the figures shows total, anterior circulation, left internal carotid, right internal carotid, posterior circulation, left vertebral, right vertebra. Right panel: ratios of flow in individual arteries versus total cerebral blood flow. From left to right the figure shows left internal carotid, right internal carotid, left vertebral, right vertebra. The figures demonstrate visually (1) the range of total cerebral blood flow, (2) the contribution from individual arteries is predominantly from the anterior circulation, that flow through (3) carotid and (4) vertebral arteries are symmetric on the left and right sides.

Limitations

Our models are only exploiting image-based information, however we know that CBF is associated with hematocrit, sex and age. (Borzage et al., 2016; Bush et al., 2016). Thus, excluding this data limits our modeling but it also allows our approach to be suitable for scenarios in which this information is either never collected, or removed to anonymize the datasets. We assumed fixed ratios between vessel flow and total CBF. This assumption is implicit on the understanding that the flow would be proportional to the volume of tissue perfused by these vessels, and that those tissues would be present in a fixed volumetric ratio. However, the ratio of these structures may change in development or senescence which limits our model application to the developing brain and provides an opportunity for further development. (Bethlehem et al., 2021). Our model assumed the PC MR slice contains the ICAs and VAs, however with an exceptionally poor localization the image might include the common carotid arteries or basilar artery. If these are incorrectly identified as ICA or VA arteries, then it would cause overestimation of the total CBF. One opportunity for future research is developing models to impute cerebral blood flow from images of common carotid or basilar arteries. Our results demonstrate flow in the two ICAs and two VAs are statistically equivalent. However, we did not include patients with systemic or cerebrovascular disease, nor patients with profound variance in the anatomy of their arteries (e.g. patients who actually lack an ICA or VA), thus our data demonstrating lateral symmetry might not generalize to populations with vascular disease. However, patients with known vascular disease or profoundly abnormal vessels may benefit from the more time-insensitive option of ASL MR. In contrast the PC MR approach is suitable for large studies and population-based screening. Moreover, our observed inter-vessel relationships could potentially be useful in recognizing deranged blood flow distribution in conditions such as steno-occlusive disease.

We did not explore other approaches to acquisition or analysis of the PC MR image data. Our secondary analysis was unable to change the prior approach to acquisition and exploring image processing approaches is beyond the scope of this project. Moreover, our imputation results are based on vascular physiology, not imaging technology. Thus, our results remain valid even with improved PC MR data acquisitions or image analyses. Our dimensionless approach means that even if our imaging methods are biased compared to approaches taken by others, they can use our imputation models. If our image acquisition methods have higher variance compared to approaches by others, our modeling error will be conservative and overestimate the variance compared to those their improved methods.

Conclusion

Phase contrast MR is an efficient and effective way to assess total cerebral blood flow. Our results indicate that using our imputation models based on at least one ICA flow measurement provides lower variance in results than ASL, and higher intraclass correlation than test-retest PC MR. Therefore, our methods are an important set of equations describing vascular physiology. Our equations enable a new approach for dealing with real-world data and make it easier to use PC MR to obtain measurements of CBF in large numbers of patients or volunteers.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Childrenʼs Hospital Los Angeles Committee on Clinical Investigations (CCI 11-00083). Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

Research reported in this publication was supported by The National Heart, Lung, and Blood Institute of the National Institutes of Health under awards 5K25HL153954, 1U01-HL-117718-01 and Minority Supplement to Grant 1U01-HL-117718-01.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Alsop, D. C., Dai, W., Grossman, M., and Detre, J. A. (2010). Arterial spin labeling blood flow MRI: Its role in the early characterization of alzheimer's disease. J. Alzheimers Dis. 20 (3), 871–880. doi:10.3233/JAD-2010-091699

PubMed Abstract | CrossRef Full Text | Google Scholar

Bethlehem, R. A. I., Seidlitz, J., White, S. R., Vogel, J. W., Anderson, K. M., Adamson, C., et al. 2021.

Borzage, M. T., Bush, A. M., Choi, S., Nederveen, A. J., Václavů, L., Coates, T. D., et al. (2016). Predictors of cerebral blood flow in patients with and without anemia. J. Appl. Physiol. (1985) 120 (8), 976–981. doi:10.1152/japplphysiol.00994.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Bush, A. M., Borzage, M. T., Choi, S., Václavů, L., Tamrazi, B., Nederveen, A. J., et al. (2016). Determinants of resting cerebral blood flow in sickle cell disease. Am. J. Hematol. 91 (9), 912–917. doi:10.1002/ajh.24441

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Liu, P., Volkow, N. D., Pan, Y., and Du, C. (2016). Cocaine attenuates blood flow but not neuronal responses to stimulation while preserving neurovascular coupling for resting brain activity. Mol. Psychiatry 21 (10), 1408–1416. doi:10.1038/mp.2015.185

PubMed Abstract | CrossRef Full Text | Google Scholar

Christie, I. C., Price, J., Edwards, L., Muldoon, M., Meltzer, C. C., and Jennings, J. R. (2008). Alcohol consumption and cerebral blood flow among older adults. Alcohol 42 (4), 269–275. doi:10.1016/j.alcohol.2008.03.132

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, W., Fong, T., Jones, R. N., Marcantonio, E., Schmitt, E., Inouye, S. K., et al. (2017). Effects of arterial transit delay on cerebral blood flow quantification using arterial spin labeling in an elderly cohort. J. Magn. Reson Imaging 45 (2), 472–481. doi:10.1002/jmri.25367

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanaoka, T., Kimura, N., Aso, Y., Takemaru, M., Kimura, Y., Ishibashi, M., et al. (2016). Relationship between white matter lesions and regional cerebral blood flow changes during longitudinal follow up in Alzheimer's disease. Geriatr. Gerontol. Int. 16 (7), 836–842. doi:10.1111/ggi.12563

PubMed Abstract | CrossRef Full Text | Google Scholar

Heijtel, D. F., Mutsaerts, H. J., Bakker, E., Schober, P., Stevens, M. F., Petersen, E. T., et al. (2014). Accuracy and precision of pseudo-continuous arterial spin labeling perfusion during baseline and hypercapnia: A head-to-head comparison with 1⁵O H₂O positron emission tomography. Neuroimage 92, 182–192. doi:10.1016/j.neuroimage.2014.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Jennings, J. R., Heim, A. F., Kuan, D. C., Gianaros, P. J., Muldoon, M. F., and Manuck, S. B. (2013). Use of total cerebral blood flow as an imaging biomarker of known cardiovascular risks. Stroke 44 (9), 2480–2485. doi:10.1161/STROKEAHA.113.001716

PubMed Abstract | CrossRef Full Text | Google Scholar

King, K. S., Sheng, M., Liu, P., Maroules, C. D., Rubin, C. D., Peshock, R. M., et al. (2018). Detrimental effect of systemic vascular risk factors on brain hemodynamic function assessed with MRI. Neuroradiol. J. 31 (3), 253–261. doi:10.1177/1971400917750375

PubMed Abstract | CrossRef Full Text | Google Scholar

Koerte, I., Haberl, C., Schmidt, M., Pomschar, A., Lee, S., Rapp, P., et al. (2013). Inter- and intra-rater reliability of blood and cerebrospinal fluid flow quantification by phase-contrast MRI. J. Magn. Reson Imaging 38 (3), 655–662. doi:10.1002/jmri.24013

PubMed Abstract | CrossRef Full Text | Google Scholar

Lassen, N. A. (1959). Cerebral blood flow and oxygen consumption in man. Physiol. Rev. 39 (2), 183–238. doi:10.1152/physrev.1959.39.2.183

PubMed Abstract | CrossRef Full Text | Google Scholar

Leaver, A. M., Vasavada, M., Joshi, S. H., Wade, B., Woods, R. P., Espinoza, R., et al. (2019). Mechanisms of antidepressant response to electroconvulsive therapy studied with perfusion magnetic resonance imaging. Biol. Psychiatry 85 (6), 466–476. doi:10.1016/j.biopsych.2018.09.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Leijenaar, J. F., van Maurik, I. S., Kuijer, J. P. A., van der Flier, W. M., Scheltens, P., Barkhof, F., et al. (2017). Lower cerebral blood flow in subjects with Alzheimer's dementia, mild cognitive impairment, and subjective cognitive decline using two-dimensional phase-contrast magnetic resonance imaging. Alzheimers Dement. (Amst). 9, 76–83. doi:10.1016/j.dadm.2017.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Lu, H., Filbey, F. M., Pinkham, A. E., McAdams, C. J., Adinoff, B., et al. (2014). Automatic and reproducible positioning of phase-contrast MRI for the quantification of global cerebral blood flow. PLoS One 9 (5), e95721. doi:10.1371/journal.pone.0095721

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Qi, Y., Lin, Z., Guo, Q., Wang, X., and Lu, H. (2019). Assessment of cerebral blood flow in neonates and infants: A phase-contrast MRI study. Neuroimage 185, 926–933. doi:10.1016/j.neuroimage.2018.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, H., Clingman, C., Golay, X., and van Zijl, P. C. (2004). Determining the longitudinal relaxation time (T1) of blood at 3.0 Tesla. Magn. Reson Med. 52 (3), 679–682. doi:10.1002/mrm.20178

PubMed Abstract | CrossRef Full Text | Google Scholar

MacIntosh, B. J., Swardfager, W., Robertson, A. D., Tchistiakova, E., Saleem, M., Oh, P. I., et al. (2015). Regional cerebral arterial transit time hemodynamics correlate with vascular risk factors and cognitive function in men with coronary artery disease. AJNR Am. J. Neuroradiol. 36 (2), 295–301. doi:10.3174/ajnr.A4094

PubMed Abstract | CrossRef Full Text | Google Scholar

Oshima, T., Karasawa, F., and Satoh, T. (2002). Effects of propofol on cerebral blood flow and the metabolic rate of oxygen in humans. Acta Anaesthesiol. Scand. 46 (7), 831–835. doi:10.1034/j.1399-6576.2002.460713.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Paniukov, D., Lebel, R. M., Giesbrecht, G., and Lebel, C. (2020). Cerebral blood flow increases across early childhood. Neuroimage 204, 116224. doi:10.1016/j.neuroimage.2019.116224

PubMed Abstract | CrossRef Full Text | Google Scholar

Prohovnik, I., Hurlet-Jensen, A., Adams, R., De Vivo, D., and Pavlakis, S. G. (2009). Hemodynamic etiology of elevated flow velocity and stroke in sickle-cell disease. J. Cereb. Blood Flow. Metab. 29 (4), 803–810. doi:10.1038/jcbfm.2009.6

PubMed Abstract | CrossRef Full Text | Google Scholar

Robertson, A. D., Matta, G., Basile, V. S., Black, S. E., Macgowan, C. K., Detre, J. A., et al. (2017). Temporal and spatial variances in arterial spin-labeling are inversely related to large-artery blood velocity. AJNR Am. J. Neuroradiol. 38 (8), 1555–1561. doi:10.3174/ajnr.A5257

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabayan, B., van der Grond, J., Westendorp, R. G., Jukema, J. W., Ford, I., Buckley, B. M., et al. (2013). Total cerebral blood flow and mortality in old age: A 12-year follow-up study. Neurology 81 (22), 1922–1929. doi:10.1212/01.wnl.0000436618.48402.da

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakhare, A. R., Barisano, G., and Pa, J. (2019). Assessing test-retest reliability of phase contrast MRI for measuring cerebrospinal fluid and cerebral blood flow dynamics. Magn. Reson Med. 82 (2), 658–670. doi:10.1002/mrm.27752

PubMed Abstract | CrossRef Full Text | Google Scholar

Thalman, S. W., Powell, D. K., Ubele, M., Norris, C. M., Head, E., and Lin, A. L. (2019). Brain-blood partition coefficient and cerebral blood flow in canines using calibrated short TR recovery (CaSTRR) correction method. Front. Neurosci. 13, 1189. doi:10.3389/fnins.2019.01189

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, D. L., Lythgoe, M. F., van der Weerd, L., Ordidge, R. J., and Gadian, D. G. (2006). Regional variation of cerebral blood flow and arterial transit time in the normal and hypoperfused rat brain measured using continuous arterial spin labeling MRI. J. Cereb. Blood Flow. Metab. 26 (2), 274–282. doi:10.1038/sj.jcbfm.9600185

PubMed Abstract | CrossRef Full Text | Google Scholar

Vernooij, M. W., van der Lugt, A., Ikram, M. A., Wielopolski, P. A., Vrooman, H. A., Hofman, A., et al. (2008). Total cerebral blood flow and total brain perfusion in the general population: The rotterdam scan study. J. Cereb. Blood Flow. Metab. 28 (2), 412–419. doi:10.1038/sj.jcbfm.9600526

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, F. N., Xu, S., Spaeth, A., Galli, O., Zhao, K., Fang, Z., et al. (2019). Test-retest reliability of cerebral blood flow for assessing brain function at rest and during a vigilance task. Neuroimage 193, 157–166. doi:10.1016/j.neuroimage.2019.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, C., Lu, W., Qiu, J., Wang, F., Li, J., and Wang, L. (2020). Alterations of the whole cerebral blood flow in patients with different total cerebral small vessel disease burden. Front. Aging Neurosci. 12, 175. doi:10.3389/fnagi.2020.00175

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: magnetic resonance imaging (MRI), cerebral blood flow (CBF), phase contrast (PC), internal carotid artery (ICA), vertebral artery

Citation: Shah P, Doyle E, Wood JC and Borzage MT (2023) Imputation models and error analysis for phase contrast MR cerebral blood flow measurements. Front. Physiol. 14:1096297. doi: 10.3389/fphys.2023.1096297

Received: 11 November 2022; Accepted: 24 January 2023;
Published: 20 February 2023.

Edited by:

Alex Bhogal, Utrecht University, Netherlands

Reviewed by:

Mikio Hiura, Aomori University, Japan
Zhensen Chen, Fudan University, China
Xiang Kong, Nanjing General Hospital of Nanjing Military Command, China

Copyright © 2023 Shah, Doyle, Wood and Borzage. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Matthew T. Borzage, Borzage@usc.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.