Skip to main content

ORIGINAL RESEARCH article

Front. Neurosci., 05 November 2019
Sec. Brain Imaging Methods
This article is part of the Research Topic Metabolic and Vascular Imaging Biomarkers for Brain Aging and Alzheimer’s Disease View all 20 articles

Brain–Blood Partition Coefficient and Cerebral Blood Flow in Canines Using Calibrated Short TR Recovery (CaSTRR) Correction Method

\r\nScott W. Thalman,Scott W. Thalman1,2David K. Powell,David K. Powell1,3Margo UbeleMargo Ubele2Christopher M. Norris,Christopher M. Norris2,4Elizabeth Head,Elizabeth Head5,6Ai-Ling Lin,,,*Ai-Ling Lin1,2,4,7*
  • 1F. Joseph Halcomb III, Department of Biomedical Engineering, University of Kentucky, Lexington, KY, United States
  • 2Sanders–Brown Center on Aging, University of Kentucky, Lexington, KY, United States
  • 3Magnetic Resonance Imaging and Spectroscopy Center, University of Kentucky, Lexington, KY, United States
  • 4Department of Pharmacology and Nutritional Sciences, University of Kentucky, Lexington, KY, United States
  • 5Department of Pathology and Laboratory Medicine, University of California, Irvine, Irvine, CA, United States
  • 6University of California Irvine Institute for Memory Impairments and Neurological Disorders (UCI MIND), University of California, Irvine, Irvine, CA, United States
  • 7Department of Neuroscience, University of Kentucky, Lexington, KY, United States

The brain–blood partition coefficient (BBPC) is necessary for quantifying cerebral blood flow (CBF) when using tracer based techniques like arterial spin labeling (ASL). A recent improvement to traditional MRI measurements of BBPC, called calibrated short TR recovery (CaSTRR), has demonstrated a significant reduction in acquisition time for BBPC maps in mice. In this study CaSTRR is applied to a cohort of healthy canines (n = 17, age = 5.0 – 8.0 years) using a protocol suited for application in humans at 3T. The imaging protocol included CaSTRR for BBPC maps, pseudo-continuous ASL for CBF maps, and high resolution anatomical images. The standard CaSTRR method of normalizing BBPC to gadolinium-doped deuterium oxide phantoms was also compared to normalization using hematocrit (Hct) as a proxy value for blood water content. The results show that CaSTRR is able to produce high quality BBPC maps with a 4 min acquisition time. The BBPC maps demonstrate significantly higher BBPC in gray matter (0.83 ± 0.05 mL/g) than in white matter (0.78 ± 0.04 mL/g, p = 0.006). Maps of CBF acquired with pCASL demonstrate a negative correlation between gray matter perfusion and age (p = 0.003). Voxel-wise correction for BBPC is also shown to improve contrast to noise ratio between gray and white matter in CBF maps. A novel aspect of the study was to show that that BBPC measurements can be calculated based on the known Hct of the blood sample placed in scanner. We found a strong correlation (R2 = 0.81 in gray matter, R2 = 0.59 in white matter) established between BBPC maps normalized to the doped phantoms and BBPC maps normalized using Hct. This obviates the need for doped water phantoms which simplifies both the acquisition protocol and the post-processing methods. Together this suggests that CaSTRR represents a feasible, rapid method to account for BBPC variability when quantifying CBF. As canines have been used widely for aging and Alzheimer’s disease studies, the CaSTRR method established in the animals may further improve CBF measurements and advance our understanding of cerebrovascular changes in aging and neurodegeneration.

Introduction

When using tracer-based techniques like arterial spin labeling (ASL) to quantify cerebral blood flow (CBF), it is necessary to determine the partition coefficient of the tracer between the perfused tissue and the arterial blood. ASL is a non-invasive, quantitative magnetic resonance imaging (MRI) technique that uses magnetically labeled protons in the water molecules of the blood as the tracer (Williams et al., 1992; Petcharunpaisan et al., 2010; Alsop et al., 2015). So in the case of ASL, the relevant partition coefficient is the brain–blood partition coefficient of water (BBPC) which is the ratio of the solubility of water in brain tissue to the solubility of water in the blood. The BBPC is tissue specific and varies with age, species, pathology, and brain region (Herscovitch and Raichle, 1985; Kudomi et al., 2005; Leithner et al., 2010; Hirata et al., 2011). This means that BBPC should be determined experimentally for each subject.

However, the standard practice in ASL studies is to assume a constant BBPC value of 0.9 mL/g for all regions of the brain regardless of the known variability of this parameter (Herscovitch and Raichle, 1985; Alsop et al., 2015). This assumption is made because the previously published MRI methods to experimentally determine BBPC required prohibitively long acquisition times and ASL studies were generally focused on gray matter perfusion where BBPC variability was assumed to be small (Roberts et al., 1996; Leithner et al., 2010). A recent study in mice at 7T reported an 87% reduction in the acquisition time for BBPC maps using an MRI technique called calibrated short TR recovery (CaSTRR) (Thalman et al., 2019). Like previous methods, CaSTRR determines relative proton density using a series of gradient echo acquisitions with varying repetition times (TR) and then calibrates the proton density map using a set of deuterium doped phantoms which provide an absolute scale of water content. The method is accelerated in CaSTRR by using shorter TR values and using gadolinium doped water phantoms to acquire similar quality BBPC maps in a fraction of the time.

The goal of this study is to apply the CaSTRR technique to a cohort of healthy canines using a protocol suited for application in humans at 3T. To do so we acquired BBPC images using a CaSTRR protocol adapted for use on a 3T human scanner. We then acquired CBF maps using pseudo-continuous ASL (pCASL) to assess the effect of BBPC correction on CBF maps, and high resolution anatomical images using magnetization prepared rapid acquisition gradient echo (MPRAGE) to facilitate segmentation and coregistration. Finally, we compare two methods of normalizing the proton density maps using the doped water phantoms and using blood water content derived from Hct values.

Materials and Methods

Animals

All animal experiments were performed in accordance with NIH guidelines and approved by the University of Kentucky Institutional Animal Care and Use Committee (approval number #2017–2680). Middle aged beagles (n = 17, age = 5.0–8.0 years, male = 24%) were acquired as part of a longitudinal study on aging and Alzheimer’s disease. The scans in this report represent pretreatment observations, and all animals were healthy at the time of their scans. The animals were anesthetized during the MRI procedure using 3–4 mg/Kg propofol for induction and 1–4% isoflurane mixed with air for maintenance. Respiratory rate, heart rate, body temperature, and blood pressure were monitored and maintained throughout the procedure. Two 5 mL vials of blood were drawn from the jugular vein using ethylenediaminetetraacetate (EDTA) treated vials. One of these was placed in the scanner with the animal according to the CaSTRR protocol, and the other was sent for laboratory analysis including Hct (ANTECH Diagnostics, Louisville, KY, United States).

Scanning Procedure

Magnetic resonance imaging experiments were performed using a 3T Siemens Prisma scanner (Siemens, Erlangen, Germany) at the MRI and Spectroscopy Center at the University of Kentucky. The animal was placed prone with their head resting in a 155 mm diameter, 15 channel transmit/receive birdcage coil commonly used for scanning human knees. The doped water phantoms along with blood sample were centered on the top of the head. The CaSTRR, pCASL, and MPRAGE acquisitions were all acquired in a single scanning session.

Calibrated Short TR Recovery Imaging

For the CaSTRR proton density measurements a series of 2-D image stacks were acquired using a phase-spoiled, fast low-angle shot gradient echo (FLASH-GRE) sequence with varying repetition times (TR = 125, 250, 500, 1000, and 2000 ms) (Thalman et al., 2019). The shortest possible echo time (TE = 1.9 ms) was used to minimize T2 decay. Image matrix parameters were: field of view = 135 × 124 mm, matrix = 96 × 88, in-plane resolution = 1.4 × 1.4 mm, slice thickness = 3 mm, number of slices = 30, flip angle = 90°, acquisition time = 4 min, labeling offset = 12 mm (see Figure 1A). A B1 mapping was done to confirm accuracy and homogeneity of the B1 field. A simulation with a saline phantom was used to simulate the gradient echo signal over a range of flip angles (see Supplementary Figure S1).

FIGURE 1
www.frontiersin.org

Figure 1. An explanation of CaSTRR and pCASL methods used in this study. CaSTRR utilizes a series of FLASH-GRE acquisitions with varying TR which include a blood sample and gadolinium-doped deuterium samples placed on the head (A). An exponential regression is fit to the signal recovery curve for each voxel (B) yielding a map of relative proton density values (C). The relative proton density values are calibrated using either water content estimated from hematocrit (Hct) (D) or the scale of water content present in the phantoms (E). Uncorrected CBF maps are derived from pCASL acquisitions (F) and are corrected on a voxel-wise basis using the BBPC map normalized to the phantoms to create a corrected map of CBF (G).

A Qualitative proton density map was calculated for each subject in a voxel-wise manner by fitting the signal recovery curve to the mono-exponential equation S = M0[1-e^(-TR/T1)] to yield a map of M0 in arbitrary units. Next a Bayesian bias field correction (Figure 2) was applied to the M0 maps to account for inhomogeneity in the receiver coil profile (Iglesiast et al., 2016). The low spatial frequency bias field was calculated using 4th order polynomials and six Gaussian components. To avoid artificially attenuating the higher signal water phantoms, the blood and water phantoms were excluded when calculating the smooth bias field, and the correction was then applied to the entire volume (see Figures 1B,C).

FIGURE 2
www.frontiersin.org

Figure 2. A representative images of M0 as calculated by exponential regression (A), signal inhomogeneity determined by Bayesian bias field correction (B), and corrected M0 (C).

The calculation of M0 values by voxel-wise exponential regression on the signal recovery curves resulted in a range of M0 values on an arbitrary scale. Due to the qualitative nature of MRI signal measurement, the average value of M0 in brain tissue varied from subject to subject on this arbitrary scale. Therefore, the M0 maps were calibrated to absolute water content in two ways. The first follows the previously published CaSTRR technique which uses a series of five phantoms containing mixtures of deuterium oxide and distilled water at 40, 30, 20, 10, and 0% (Thalman et al., 2019). These phantoms were also doped with 0.18 mM gadobutrol (Gadavist, Bayer Healthcare Pharmaceuticals, Whippany, NJ, United States) to reduce the longitudinal relaxation rate (T1) to be similar to the T1 of brain tissue (≈ 1.2 s). Because the deuterium oxide does not produce signal in MRI this creates a scale of known water concentration from 60–100% water. Each voxel in the image can then be normalized to this internal scale to yield a measure of absolute water content. It was noted, however, that in subjects where the arbitrary M0 values of brain tissue were high, the M0 values of the doped water phantoms tended to be higher than expected. This meant that while the M0 value of tissues increased linearly with the overall image intensity, the M0 values of the phantoms exhibited a quadratic increase over the range of arbitrary M0 values measured. The result was that the ratio of average M0,tissue to average M0,phantoms used in image normalization was not constant from subject to subject but had a strong negative linear relationship to the phantom signal intensity. To correct for the relationship, linear regression of this ratio against the average M0 value of the phantoms was used to quantify the quadratic error term for each subject. This negative error term was then subtracted from the blood and tissue M0 values prior to normalization. Finally, a linear regression was calculated based on the average M0 value in each phantom and its known water content, and every voxel in the image was normalized to the resultant equation (see Figure 1E).

The second method of calibrating M0 maps utilized the arterial Hct value to determine the absolute water content of the blood sample. Water content was determined according to the equation WCblood = −0.271Hct + 0.912 (Lijnema et al., 1993). The M0 image was then normalized such that the average water content in the blood sample matched the average value calculated according to hematocrit.

Brain–blood partition coefficient maps were generated for both normalization methods by comparing the measured water content of each voxel in the brain to the average water content of the blood using the equation BBPC = WCbrain/(WCblood1.04 g/mL) (see Figure 1D).

Cerebral Blood Flow and Anatomical Imaging

The CBF maps were acquired using a pCASL sequence with a three-dimensional gradient and spin echo (GRASE) readout (Alsop et al., 2015). The acquisition parameters were as follows: TR/TE/TI = 3200/16/1400 ms, field of view = 270 × 250 × 90 mm, matrix = 96 × 88, resolution = 3 × 3 × 3 mm, acquisition time = 6:15 min (see Figure 1F).

Because the original CBF maps were created with an assumed BBPC value of 0.9 mL/g, corrected maps were generated by dividing the entire map by this value and then multiplying by the BBPC map in a voxel-wise manner. BBPC correction was performed using the BBPC maps generated using the gadolinium doped water phantoms (see Figure 1G).

Contrast to noise ratio (CNR) was calculated according to the equation CNR = (Meangray – Meanwhite)/Pooled Standard Deviation (Cohen, 1988).

Anatomical images were acquired using a high resolution T1 weighted magnetization prepared rapid acquisition gradient echo (MPRAGE) sequence as recommended for optimal use with FreeSurfer (MGH, 2009). Scan parameters were: TR/TE = 1690/2.56 ms, flip angle = 12°, field of view = 320 × 320 × 160 mm, matrix = 256 × 256 × 128, resolution = 0.8 × 0.8 × 0.8 mm, acquisition time = 10:49 min.

Image Analysis

All images were coregistered by first resampling the anatomical volumes to match the slice thickness of the CaSTRR and pCASL acquisitions. The CaSTRR and pCASL volumes were then registered to the anatomical using an intensity-based registration algorithm in Matlab (Mathworks, Natick, MA, United States) (Styner et al., 2000). The brain region of interest was extracted manually, and then segmented into gray and white regions of interest using an expectation maximization algorithm with classes for gray matter, white matter, and cerebrospinal fluid (Wells et al., 1996). To avoid partial volume effects, the gray and white matter regions of interest in each slice were eroded by two pixels. Due to image inhomogeneities in the MPRAGE acquisitions in some animals, the segmentation algorithm often failed to adequately segment gray and white matter regions in the most rostral and caudal sections of the brain. While the BBPC maps did not display these inhomogeneities, the regions of interest were not reliable for analysis. For this reason, we chose to include only the centermost 10 slices in our analysis. These regions of interest were then applied to the BBPC and CBF maps. All image analysis was performed in Matlab (Mathworks, Natick, MA, United States).

Statistical Analysis

All statistical analysis was performed in Matlab. Data is expressed as mean ± standard deviation. Gray and white matter comparisons were assessed using three-way analysis of covariance with age, sex, and tissue type as independent variables. Linear regressions against age were also controlled for sex. Values of p < 0.05 were considered statistically significant.

Results

BBPC Is Higher in Gray Matter Than in White Matter

The first comparison was drawn on BBPC maps generated using the previously published CaSTRR method of normalizing to the gadolinium doped water phantoms (see Figure 1E). The average BBPC in gray matter was 0.83 ± 0.05 mL/g which is 5.6% higher than the BBPC in white matter (0.78 ± 0.04 mL/g, p = 0.007) (see Figure 3B). When plotted against age, neither the gray nor the white matter regions of interest demonstrated a significant correlation over the age range studied (gray matter: p = 0.645, white matter: p = 0.483 (see Figure 3A).

FIGURE 3
www.frontiersin.org

Figure 3. Average BBPC in gray matter (blue) and white matter (red) regions of interest plotted against age (A) and as group averages (B). Regions corresponding to the same animal are connected with a vertical dotted line. No significant linear correlation to age was found in either region. Average BBPC in gray matter was which is 5.6% higher than in white matter (BBPCgray = 0.83 ± 0.05 mL/g, BBPCwhite 0.78 ± 0.04 mL/g, ∗∗ indicates p < 0.01, error bars represent 1 standard deviation).

Gray Matter CBF Is Negatively Correlated With Age

Next, gray and white matter perfusion were compared in uncorrected CBF maps (see Figure 1F) and again in corrected CBF maps (see Figure 1G) which used the standard CaSTRR derived BBPC maps for correction (see Figure 1E). Gray matter has 45% higher CBF than white matter in the uncorrected CBF maps (CBFgray = 81 ± 12 mL/100g/min, CBFwhite = 56 ± 12 mL/100 g/min, p < 0.001) and 53% higher CBF in the maps corrected for BBPC (cCBFgray = 73 ± 13 mL/100 g/min, cCBFwhite = 49 ± 11 mL/100 g/min, p < 0.01) (see Figure 4B). Gray matter demonstrated a significant negative correlation with age with a reduction of 7.5 ± 2.1 mL/100 g/min each year or 9% of the average perfusion (CBFgray = 128 – 7.5 Age mL/100 g/min, p = 0.003). The corrected CBF maps also revealed a reduction of 6.6 ± 2.6 mL/100 g/min/year (cCBFgray = 117 – 6.6 Age mL/100 g/min, p = 0.02), but this relationship was not significantly different in the corrected maps than the uncorrected (p = 0.81). While there appears to be a downward trend in white matter perfusion with age, this correlation was not statistically significant in the uncorrected CBF maps (p = 0.20) or the corrected maps (p = 0.33) (see Figure 4A).

FIGURE 4
www.frontiersin.org

Figure 4. Gray and white matter perfusion plotted against age (A). Plotted points represent uncorrected average CBF for each animal. Gray matter CBF demonstrates a negative linear correlation with age in both uncorrected CBF maps (CBFgray = 128 – 7.5 Age mL/100 g/min) and maps corrected using measured BBPC (cCBFgray = 117 – 6.6 Age mL/100 g/min). Linear regression was not significant in the white matter region in either case. When considered as group averages (B), the gray matter has significantly higher CBF in both uncorrected (CBFgray = 81 ± 12 mL/100 g/min, CBFwhite = 56 ± 12 mL/100 g/min), and corrected maps (cCBFgray = 73 ± 13 mL/100 g/min, cCBFwhite = 49 ± 11 mL/100 g/min). indicates p < 0.05, ∗∗ indicates p < 0.01, ∗∗∗ indicates p < 0.001, error bars represent 1 standard deviation.

BBPC Correction Improved Contrast to Noise Ratio in CBF Maps

Next, the CNR of corrected CBF maps (see Figure 1G) was compared to uncorrected CBF maps (see Figure 1F). On average BBPC correction improved CNR between gray and white matter regions of the CBF maps by 3.6% (95% confidence interval = 0.6 – 6.5%). The average uncorrected CNR was 0.81 and the average corrected CNR was 0.84.

BBPC Values Generated Using Hematocrit to Estimate Water Content Agree With Those Generated Using Water Phantoms

The final comparison was between maps normalized using the doped water phantoms (see Figure 1E) to ones normalized using Hct derived blood water content (see Figure 1D), we observed positive correlations between the BBPC values generated using these two methods. The Pearson correlation was R2 = 0.81 for gray matter indicating strong correlation between these measures in this region (see Figure 5A). Due to higher variability in the white matter regions the correlation was moderate in white matter (R2 = 0.59) (see Figure 5C). The measured BBPC values were slightly lower in maps normalized to hematocrit, though not statistically different, and Bland-Altman analysis demonstrates no significant bias in either region of interest (see Figures 5B,D). Again the BBPC in gray matter was 5.9% higher than in white matter when using Hct to normalize (BBPCgray = 0.81 ± 0.06 mL/g, BBPCwhite = 0.77 ± 0.05 mL/g, p = 0.02).

FIGURE 5
www.frontiersin.org

Figure 5. Correlation of BBPC values in maps normalized according to hematocrit (Hct) with values in maps normalized to the gadolinium-doped phantoms. Correlation is very strong for the gray matter BBPC values (R2 = 0.81) (A) and moderately strong for white matter (R2 = 0.59) (C). Bland-Altman analysis demonstrates no significant bias in either gray (B) or white matter values (D).

Discussion

The CaSTRR technique represents a significant improvement in the acquisition speed of BBPC maps. A previously published report measuring BBPC using a 1.5T human scanner acquired a single slice of the BBPC map in approximately 30 min (Roberts et al., 1996). In this study we were able to produce BBPC maps of sufficient quality to perform voxel-wise correction of CBF maps with coverage of the entire brain using a CaSTRR acquisition protocol that only required 4 min of scan time. This is an improvement over the reported CaSTRR technique in mice which required 17 min of scan time due to the much higher resolution requirement of scanning small animals at 7T (Thalman et al., 2019). A 4 min scan time is comparable to the acquisition time of the pCASL technique for CBF. The experiment was also done with commercially available equipment and pulse sequences that are directly applicable for use with a human subject. Furthermore, the CaSTRR scans were performed with the same birdcage receive coil that was used for the pCASL acquisitions. This was not possible when scanning mice at 7T and represents a significant advantage of scanning at lower fields. The greatly reduced scan time and ready availability of equipment and pulse sequences demonstrate that CaSTRR is a feasible approach to correct CBF maps using empirically measured BBPC values instead of assuming a constant value for all tissue types, pathologies, ages, and species. This technique has the potential for rapid translation to use in human studies.

While BBPC has not been previously reported in canines, our reported values of 0.83 ± 0.05 mL/g in gray matter and 0.78 ± 0.04 mL/g in white matter are lower than published reports in humans, non-human primates, and mice (Herscovitch and Raichle, 1985; Kudomi et al., 2005; Leithner et al., 2010). One possible reason for this is the temperature discrepancy between the blood sample and the brain tissue. When measured at room temperature instead of physiologic temperature, the proton density of the blood could be overestimated by as much as 5% causing BBPC to be underestimated by the same amount (Tofts, 2003). The amount of inter-species variability in BBPC values is further evidence for the importance of empirical BBPC correction when quantifying CBF.

The correlation between the two methods of normalization represents a distinct advantage of this study. Our results suggest that the CaSTRR technique can be further simplified by omitting the gadolinium doped water phantoms. One of the difficulties of this study was the non-linear signal increase observed in the water phantoms. This is possibly due to the pre-scan normalization algorithm of the scanner used for this study and/or reduced T2 decay in the water phantoms. These effects would likely be specific to a given scanner and would need to be determined empirically. However, the correlation between BBPC values derived using the water phantoms with those derived using Hct indicate that future studies using CaSTRR could rely solely on the water content determined by the hematocrit. This would also obviate the need for correcting the non-linear signal increases observed in the water phantoms as done in this study.

Another significant advantage of this study is the use of Bayesian bias field correction to account for inhomogeneities in the receiver coil sensitivity profile. The CaSTRR method described in mice assumed a sufficiently homogenous profile in the birdcage receive coil, but observed significant signal loss near the large ear canals of the mouse (Thalman et al., 2019). Other BBPC studies attempted to correct for bias field using a separate measurement on a uniform phantom (Roberts et al., 1996; Leithner et al., 2010), but it is unlikely that the receiver profile would be the same when measuring the non-uniform tissue of a live subject.

Arterial spin labeling has an inherently low signal to noise ratio because it is a subtractive technique. Including an empirical measurement of BBPC to the quantification of CBF will increase noise, as we observed in the greater variance of CBF values in the corrected maps. However, there was an improvement in contrast to noise between areas with significantly different BBPC values despite this addition of noise. This could become important when studying perfusion in models of pathology where the tissue composition is likely to change.

Canines have been widely used for aging and neurodegeneration studies (Martin et al., 2011; Vite and Head, 2014; Triani et al., 2018). There are many examples of pathologies that could affect water balance in the brain. Brain edema caused by ischemia (Rosenberg, 1999), infection (Niemoller and Tauber, 1989), or trauma (Winkler et al., 2016) can cause localized increases in free water and potentially affect the BBPC. Another important field of interest where ASL is commonly used is the study of Alzheimer’s disease (AD). The deposition of hydrophobic amyloid-β protein, which is a hallmark of AD pathology, may reduce the BBPC in regions of plaque development (Aleksis et al., 2017). AD occurs in the context of aging and typically causes pronounced brain atrophy (Bobinski et al., 1999). Both of these could result in reduced brain water content and therefore reduced BBPC. Both brain volume (Duning et al., 2005) and Hct (Jimenez et al., 1999) can also change significantly with a subject’s level of hydration. So while our result showed that BBPC correction did not affect the observed relationship between CBF and age between the ages of 5–8 year, it is possible that BBPC correction could improve sensitivity in studies where BBPC is expected to change significantly. In future studies, CaSTRR imaging could be used to study how BBPC changes in canines with pathology and could also be used to account for water balance effects when measuring perfusion.

We acknowledge that in addition to CaSTRR, other efforts have also been made to improve the BBPC measurement, CBF quantification, and gray-white CBF contrast. It has been reported that a uniform, brain-tissue-type-dependent magnetization image could be generated using a sensitivity calibration (Dai et al., 2011). Another study showed that BBPC can be improved by exploiting the partial-volume data to adjust the ratio between BBPC and the proton density-weighted image (Ahlgren et al., 2018). Furthermore, gray-white matter CBF can be enhanced with background suppression methods (van Osch et al., 2009). Here we provide another approach that can quantify BBPC rapidly, and improve CBF quantification and gray-white CBF contrast using Hct calibration.

Conclusion

In conclusion, this study demonstrates the feasibility of CaSTRR as a method to correct CBF measurements for regional and inter-subject variability in BBPC. Further, we demonstrated that the correction can be achieved using Hct calibration. The developed CaSTRR method has potential contribution for future translational studies in aging and neurodegeneration.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Ethics Statement

The animal study was reviewed and approved by the Institutional Animal Care and Use Committee of the University of Kentucky.

Author Contributions

ST was responsible for the experimental and scanning protocol design, analysis software development, image acquisition, data and statistical analyses, and manuscript preparations. DP contributed to the sequence development, scanning protocol design, technical support, and manuscript editing. MU was responsible for the animal handling, image acquisition, laboratory samples, and manuscript editing. CN and EH contributed to the project design, animal acquisition, data interpretation, and manuscript editing. A-LL was the primary investigator and contributed to the project design, interpretation of the results, and manuscript preparation.

Funding

This research was supported by the National Institute of Health (NIH) (Grant Numbers R01AG054459, RF1AG062480, R01AG056998, and T32AG057461). The 7T ClinScan small animal MRI scanner of the University of Kentucky was funded by the S10 NIH Shared Instrumentation Program Grant (Grant Number 1S10RR029541-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.

Acknowledgments

The authors would like to thank Kathy Boaz, Stephanie Krumholz, and Beverly Meacham for their assistance in animal handling and MRI experiments.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2019.01189/full#supplementary-material

FIGURE S1 | B1 mapping. The B1 mapping was aquired using a Tfl_b1map sequence with a saline phantom (salinity is the same as mammalian extracellular space). Right image proportional to flip angle, showing maximum flip angle variation of 4% over a 10 cm phantom. The FOV is comparible to beagle’s brain, which is approximately 4 cm × 5 cm × 7 cm. (Left) The magnitude image. (Right) A flip angle map where the intensity is linearly proportional to the flip angle. The examples are shown in two ROIs: 83.8 and 89.5 degrees, respectively.

References

Ahlgren, A., Wirestam, R., Knutsson, L., and Petersen, E. T. (2018). Improved calculation of the equilibrium magnetization of arterial blood in arterial spin labeling. Magn. Reson. Med. 80, 2223–2231. doi: 10.1002/mrm.27193

PubMed Abstract | CrossRef Full Text | Google Scholar

Aleksis, R., Oleskovs, F., Jaudzems, K., Pahnke, J., and Biverstal, H. (2017). Structural studies of amyloid-beta peptides: unlocking the mechanism of aggregation and the associated toxicity. Biochimie 140, 176–192. doi: 10.1016/j.biochi.2017.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Alsop, D. C., Detre, J. A., Golay, X., Gunther, M., Hendrikse, J., Hernandez-Garcia, L., et al. (2015). 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. 73, 102–116. doi: 10.1002/mrm.25197

PubMed Abstract | CrossRef Full Text | Google Scholar

Bobinski, M., De Leon, M. J., Wegiel, J., Desanti, S., Convit, A., Saint Louis, L. A., et al. (1999). The histological validation of post mortem magnetic resonance imaging-determined hippocampal volume in Alzheimer’s disease. Neuroscience 95, 721–725. doi: 10.1016/s0306-4522(99)00476-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, J. (1988). Statistical Power Analysis for the Behavioral Sciences. Hillsdale, NJ: Lawrence Erlbaum Associates.

Google Scholar

Dai, W., Robson, P. M., Shankaranarayanan, A., and Alsop, D. C. (2011). Sensitivity calibration with a uniform magnetization image to improve arterial spin labeling perfusion quantification. Magn. Reson. Med. 66, 1590–1600. doi: 10.1002/mrm.22954

PubMed Abstract | CrossRef Full Text | Google Scholar

Duning, T., Kloska, S., Steinstrater, O., Kugel, H., Heindel, W., and Knecht, S. (2005). Dehydration confounds the assessment of brain atrophy. Neurology 64, 548–550. doi: 10.1212/01.wnl.0000150542.16969.cc

PubMed Abstract | CrossRef Full Text | Google Scholar

Herscovitch, P., and Raichle, M. E. (1985). What is the correct value for the brain–blood partition coefficient for water? J. Cereb. Blood Flow Metab. 5, 65–69. doi: 10.1038/jcbfm.1985.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirata, K., Hattori, N., Katoh, C., Shiga, T., Kuroda, S., Kubo, N., et al. (2011). Regional partition coefficient of water in patients with cerebrovascular disease and its effect on rCBF assessment. Nucl. Med. Commun. 32, 63–70. doi: 10.1097/MNM.0b013e3283412106

PubMed Abstract | CrossRef Full Text | Google Scholar

Iglesiast, J. E., Paz-Alonso, P., Lerma-Usabiaga, G., Insausti, R., Miller, K., and Caballero-Gaudes, C. (2016). “Simultaneous Bayesian correction of slab boundary artifacts and bias field for high resolution ex vivo MRI,” in Proceedings of the 13th International Symposium on Biomedical Imaging (ISBI), (Piscataway, NJ: IEEE), 1384–1387.

Google Scholar

Jimenez, C., Melin, B., Koulmann, N., Allevard, A. M., Launay, J. C., and Savourey, G. (1999). Plasma volume changes during and after acute variations of body hydration level in humans. Eur. J. Appl. Physiol. Occup. Physiol. 80, 1–8. doi: 10.1007/s004210050550

PubMed Abstract | CrossRef Full Text | Google Scholar

Kudomi, N., Hayashi, T., Teramoto, N., Watabe, H., Kawachi, N., Ohta, Y., et al. (2005). Rapid quantitative measurement of CMRO(2) and CBF by dual administration of (15)O-labeled oxygen and water during a single PET scan-a validation study and error analysis in anesthetized monkeys. J. Cereb. Blood Flow Metab. 25, 1209–1224. doi: 10.1038/sj.jcbfm.9600118

PubMed Abstract | CrossRef Full Text | Google Scholar

Leithner, C., Muller, S., Fuchtemeier, M., Lindauer, U., Dirnagl, U., and Royl, G. (2010). Determination of the brain-blood partition coefficient for water in mice using MRI. J. Cereb. Blood Flow Metab. 30, 1821–1824. doi: 10.1038/jcbfm.2010.160

PubMed Abstract | CrossRef Full Text | Google Scholar

Lijnema, T. H., Huizenga, J. R., Jager, J., Mackor, A. J., and Gips, C. H. (1993). Gravimetric determination of the water concentration in whole blood, plasma and erythrocytes and correlations with hematological and clinicochemical parameters. Clin. Chim. Acta 214, 129–138. doi: 10.1016/0009-8981(93)90105-d

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, S. B., Dowling, A. L., and Head, E. (2011). Therapeutic interventions targeting beta amyloid pathogenesis in an aging dog model. Curr. Neuropharmacol. 9, 651–661. doi: 10.2174/157015911798376217

PubMed Abstract | CrossRef Full Text | Google Scholar

MGH, (2009). FreeSurfer Suggested Morphometry Protocols. Available at: http://surfer.nmr.mgh.harvard.edu/fswiki/FreeSurferWiki?action=AttachFile&do=get&target=FreeSurfer_Suggested_Morphometry_Protocols.pdf (accessed April 29, 2019).

Google Scholar

Niemoller, U. M., and Tauber, M. G. (1989). Brain edema and increased intracranial pressure in the pathophysiology of bacterial meningitis. Eur. J. Clin. Microbiol. Infect. Dis. 8, 109–117. doi: 10.1007/bf01963892

PubMed Abstract | CrossRef Full Text | Google Scholar

Petcharunpaisan, S., Ramalho, J., and Castillo, M. (2010). Arterial spin labeling in neuroimaging. World J. Radiol. 2, 384–398. doi: 10.4329/wjr.v2.i10.384

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberts, D. A., Rizi, R., Lenkinski, R. E., and Leigh, J. S. Jr. (1996). Magnetic resonance imaging of the brain: blood partition coefficient for water: application to spin-tagging measurement of perfusion. J. Magn. Reson. Imaging 6, 363–366. doi: 10.1002/jmri.1880060217

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosenberg, G. A. (1999). Ischemic brain edema. Prog. Cardiovasc. Dis. 42, 209–216.

PubMed Abstract | Google Scholar

Styner, M., Brechbuhler, C., Szckely, G., and Gerig, G. (2000). Parametric estimate of intensity inhomogeneities applied to MRI. IEEE Trans. Med. Imaging 19, 153–165. doi: 10.1109/42.845174

PubMed Abstract | CrossRef Full Text | Google Scholar

Thalman, S. W., Powell, D. K., and Lin, A. L. (2019). Novel calibrated short TR recovery (CaSTRR) method for brain-blood partition coefficient correction enhances gray-white matter contrast in blood flow measurements in mice. Front. Neurosci. 13:308. doi: 10.3389/fnins.2019.00308

PubMed Abstract | CrossRef Full Text | Google Scholar

Tofts, P. S. (ed.) (2003). “PD: Proton Density of Tissue Water,” in Quantitative MRI of the Brain, (Chichester: John Wiley).

Google Scholar

Triani, F., Tramutola, A., Di Domenico, F., Sharma, N., Butterfield, D. A., Head, E., et al. (2018). Biliverdin reductase-A impairment links brain insulin resistance with increased Abeta production in an animal model of aging: implications for Alzheimer disease. Biochim. Biophys. Acta Mol. Basis Dis. 1864, 3181–3194. doi: 10.1016/j.bbadis.2018.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

van Osch, M. J., Teeuwisse, W. M., Van Walderveen, M. A., Hendrikse, J., Kies, D. A., and Van Buchem, M. A. (2009). Can arterial spin labeling detect white matter perfusion signal? Magn. Reson. Med. 62, 165–173. doi: 10.1002/mrm.22002

PubMed Abstract | CrossRef Full Text | Google Scholar

Vite, C. H., and Head, E. (2014). Aging in the canine and feline brain. Vet. Clin. North Am. Small Anim. Pract. 44, 1113–1129. doi: 10.1016/j.cvsm.2014.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Wells, W. M., Grimson, W. E. L., Kikinis, R., and Jolesz, F. A. (1996). Adaptive segmentation of MRI data. IEEE Trans. Med. Imaging 15, 429–442. doi: 10.1109/42.511747

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams, D. S., Detre, J. A., Leigh, J. S., and Koretsky, A. P. (1992). Magnetic resonance imaging of perfusion using spin inversion of arterial water. Proc. Natl Acad. Sci. U.S.A. 89, 212–216.

PubMed Abstract | Google Scholar

Winkler, E. A., Minter, D., Yue, J. K., and Manley, G. T. (2016). Cerebral edema in traumatic brain injury: pathophysiology and prospective therapeutic targets. Neurosurg. Clin. N. Am. 27, 473–488. doi: 10.1016/j.nec.2016.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: cerebral blood flow, brain–blood partition coefficient, calibrated short TR recovery, arterial spin labeling, perfusion, magnetic resonance imaging, canines

Citation: Thalman SW, Powell DK, Ubele M, Norris CM, 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

Received: 18 June 2019; Accepted: 21 October 2019;
Published: 05 November 2019.

Edited by:

Federico Giove, Centro Fermi – Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi, Italy

Reviewed by:

Eric R. Muir, Stony Brook University, United States
Danny J. J. Wang, University of Southern California, United States

Copyright © 2019 Thalman, Powell, Ubele, Norris, Head and Lin. 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: Ai-Ling Lin, ailing.lin@uky.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.