Skip to main content

ORIGINAL RESEARCH article

Front. Med., 07 April 2021
Sec. Nuclear Medicine

Model Corrected Blood Input Function to Compute Cerebral FDG Uptake Rates From Dynamic Total-Body PET Images of Rats in vivo

  • 1Department of Radiology and Medical Imaging, University of Virginia, Charlottesville, VA, United States
  • 2Department of Experimental Physiology and Pathophysiology, Medical University of Białystok, Białystok, Poland
  • 3Preclinical Imaging Division, Bruker Biospin, Billerica, MA, United States
  • 4Department of Biomedical Engineering, University of Virginia, Charlottesville, VA, United States
  • 5Cardiovascular Research Center, University of Virginia, Charlottesville, VA, United States

Recently, we developed a three-compartment dual-output model that incorporates spillover (SP) and partial volume (PV) corrections to simultaneously estimate the kinetic parameters and model-corrected blood input function (MCIF) from dynamic 2-[18F] fluoro-2-deoxy-D-glucose positron emission tomography (FDG PET) images of mouse heart in vivo. In this study, we further optimized this model and utilized the estimated MCIF to compute cerebral FDG uptake rates, Ki, from dynamic total-body FDG PET images of control Wistar–Kyoto (WKY) rats and compared to those derived from arterial blood sampling in vivo. Dynamic FDG PET scans of WKY rats (n = 5), fasted for 6 h, were performed using the Albira Si Trimodal PET/SPECT/CT imager for 60 min. Arterial blood samples were collected for the entire imaging duration and then fitted to a seven-parameter function. The 60-min list mode PET data, corrected for attenuation, scatter, randoms, and decay, were reconstructed into 23 time bins. A 15-parameter dual-output model with SP and PV corrections was optimized with two cost functions to compute MCIF. A four-parameter compartment model was then used to compute cerebral Ki. The computed area under the curve (AUC) and Ki were compared to that derived from arterial blood samples. Experimental and computed AUCs were 1,893.53 ± 195.39 kBq min/cc and 1,792.65 ± 155.84 kBq min/cc, respectively (p = 0.76). Bland–Altman analysis of experimental vs. computed Ki for 35 cerebral regions in WKY rats revealed a mean difference of 0.0029 min−1 (~13.5%). Direct (AUC) and indirect (Ki) comparisons of model computations with arterial blood sampling were performed in WKY rats. AUC and the downstream cerebral FDG uptake rates compared well with that obtained using arterial blood samples. Experimental vs. computed cerebral Ki for the four super regions including cerebellum, frontal cortex, hippocampus, and striatum indicated no significant differences.

Introduction

Noninvasive determination of blood input function from dynamic 2-[18F] fluoro-2-deoxy-D-glucose positron emission tomography (FDG PET) (Table 1) images of rodents including mice and rats has been challenging. With the availability of whole-body rodent scanners, obtaining image-derived input function (IDIF) from different sources including the inferior vena cava (IVC), left ventricular blood pool (LVBP), and the carotid arteries is possible. Due to the limited spatial resolution of the rodent scanners and the size of the regions (~1–3 mm), IDIF derived from these regions are susceptible to incomplete radioactive recovery or partial volume (PV) averaging effects and hence spillover (SP) of radioactivity into surrounding regions of interest (ROIs). Recent work by others (1) and from our laboratory has primarily utilized the LVBP (2) and the IVC (3) to derive the IDIF. The IDIF derived from the IVC is underestimated due to PV effects, while the IDIF derived from the LVBP is overestimated due to SP effects. Recent work from our laboratory has optimized a dual-output model for the myocardial tissue and the LVBP in a three-compartment kinetic model to simultaneously estimate model-corrected blood input function (MCIF) and the kinetic model parameters with SP and PV corrections to compute the rate of myocardial FDG uptake, Ki, from dynamic FDG PET images of control BL/6 mouse hearts in vivo (2). Recently, this model was adapted to compute myocardial Ki from total-body FDG PET images of control Wistar–Kyoto (WKY) and experimental spontaneously hypertensive rats (SHR) in vivo (4).

TABLE 1
www.frontiersin.org

Table 1. Abbreviations.

In this study, we further optimized the dual-output model with two cost functions to compute cerebral Ki with MCIF derived from the LV blood pool and compared to that derived from arterial blood sampling in WKY rats in vivo.

Materials and Methods

The image and blood sampling data used to support the findings of this study are available from the corresponding author upon request.

Five (n = 5) male WKY rats, purchased from Charles River (Kingston, NY), were housed under controlled conditions (temperature 21 ± 1°C, humidity 60 ± 10%, 12-h light/12-h dark cycle, and free access to standard rat chow and water). Animal experiments were approved by the Institutional Animal Care and Use Committee of the University of Virginia.

Dynamic FDG PET Imaging

Five male WKY rats, 3–5 months of age, were imaged using a state-of-the-art Albira Si Trimodal PET/SPECT/CT scanner (5). The Albira PET imager is a three-ring scanner with an axial field of view (FOV) of 150 mm and trans-axial FOV of 80 mm, thereby enabling dynamic PET imaging of both the heart and the brain, in the same FOV, of 3–5-month-old rats weighing on an average of 300 g. The volumetric spatial resolution of the PET system ranges from 0.41 to 0.87 mm3 at the axial center with a sensitivity of 11% at the center of the FOV. The noise equivalent count rate (NECR) for the rat phantom was measured to be 240 kcps at 23 MBq (6). Dynamic FDG PET imaging of the rats, fasted for 5–6 h, was performed using a similar protocol as described recently in rats (3, 4). Briefly, a 60-min list mode (LM) acquisition, under anesthesia, was initiated, followed by intravenous injection of ~200–300 μCi of FDG via a tail vein catheter over a period of ~10 s. A three-bed CT scan followed the emission scan for attenuation correction. A small-animal gating and monitoring system (Small Animal Instruments, Inc., model 1025L for PET) was used for continuously monitoring respiration and core body temperature (using a rectal probe). The LM data were then sorted into 23 time bins [frames, time (s): 11,8; 1,12; 2,60; 1,180; 8,400], and the sinograms were reconstructed with attenuation, scatter, randoms, and decay corrections using maximum-likelihood expectation maximization (MLEM) algorithm with six iterations and 0.75-mm isotropic voxel resolution (3).

Arterial Blood Sampling

Arterial blood sampling was performed in WKY rats following a similar protocol developed in mice (2). Twelve arterial blood samples (~0.2 ml) were collected for the entire imaging duration from the same rats more frequently at the early time points (10, 20, 36, and 49 s, 1, 1.2, 2.2, and 3.3 min) and at lower frequency at the later time points (5, 10, 30, and 60 min post FDG administration). Inactin hydrate (thiobutabarbital sodium salt hydrate) was used for our studies, as a longer duration of anesthesia was needed to perform catheterization for arterial blood sampling along with dynamic FDG PET imaging followed by CT. An incision was made under Inactin anesthesia (100 mg/kg body weight, injected intraperitoneally) in the throat area, and the carotid artery was exposed by blunt dissection. The artery was tied off above the insertion site with 4–0 silk ligature and temporarily occluded distal to the insertion site using pressure from an untied ligature. A catheter was formed from PE-50 tubing and placed in the carotid artery. The catheter was ligated with silk in three places and the pressure was removed. Blood draws were followed by flushing the catheter with a small volume of heparinized saline. Each whole-blood sample was weighed, and FDG activity was measured in a Hidex Automatic γ-Counter. The arterial blood sample data were first fitted to a seven-parameter function as described (2) to serve as the ground truth. The seven-parameter FDG blood input model has been validated by others and us for quantifying rates of FDG uptake in rodent hearts and brains. This model accurately captures the wash-in and wash-out of FDG kinetics in blood (7). The peak of this blood model was forced through the maximum tissue activity value obtained from the image-derived myocardium time activity curves (TACs) for the entire dynamic range (0–60 min) to compensate for artificially low early blood sample values as a result of limitations in blood sample timing. This was then utilized to compute experimental cerebral Ki using a four-parameter compartment model as described (3).

Cerebral Ki Modeling

Brain TACs for four super regions (cerebellum, striatum, cortex, and hippocampus) were generated through the application of the W. Schiffer T2 rat brain atlas onto dynamic PET images using PMOD (pixel-wise modeling, pmod.com). Volumes and average TACs were computed for a total of 58 volumes of interest (VOIs) consisting of 26 regions split into two different hemispheres [left (L)/right (R)] and an additional six regions without lateral split. Thirty-five of these VOIs were binned and combined to generate volumes and average TACs for the four super regions. The 58 regions generated directly from the atlas, as well as the four super regions, were analyzed in Matlab (MathWorks Inc., Natick, MA). Kinetic modeling for all regions was completed using a four-parameter compartment model (3).

Modified Objective Function

A 15-parameter dual-output model (2) with SP and PV corrections optimized the following objective functions:

O1(p)=i=1n[(ModelIDIF,i             PETIDIF,i)2+(Modelmyo,iPETmyo,i)2]    (1)
O2(p)=[(ModelPeakIDIFPETPeakIDIF)2             +(ModelPeakmyoPETPeakmyo)2]    (2)
O(p)=O1(p)+  O2(p),    (3)

where ModelIDIF and Modelmyo are model output equations (Eqs. 4 and 5) and PETIDIF and PETmyo are image-derived blood and myocardium TACs, respectively, as described (2).

ModelIDIF,i=tbitei[SmbCT(t)+rbCa(t)]dttei-tbi;    (4)
Modelmyo,i=tbitei[rmCT(t)+SbmCa(t)]dttei-tbi.    (5)

Smb and Sbm are SP contamination factors from the myocardium to the blood and vice versa; rb and rm are recovery coefficients (RCs) for blood pool and myocardium, respectively, and tb and te are the beginning and end of a time frame. CT(t), the model tissue, was obtained by solving FDG transport differential equations from blood to tissue spaces (2). Ca(t) is the seven-parameter model blood for FDG transport as described (2, 7).

The second objective function, O2(p), minimizes the square of the difference between the model and image-derived blood and myocardium peak values. The ModelPeak in Eq. (2) was computed from the model equations for the IDIF (Model_IDIF) and myocardium (Model_myo) (Eqs. 4 and 5, respectively). The PETPeak values were derived from the dynamic PET images for both the blood and the myocardium for each rat. Optimization of the objective function, O(p), using nonlinear regression analysis written in Matlab, resulted in the estimation of MCIF, which was then used to compute cerebral Ki, defined as:

Ki=K1k3(k2+k3),    (6)

using a four-parameter compartment model,

Cm(t)=1(t2t1)t1t2{(1TBV)[K1k3k2+k30TCp(u)du             +K1k2k2+k30TCp(u).e(k2+k3)(Tu)du]+TBV.Cp(T)} dT,    (7)

as described (3).

The computed area under the curve (AUC under MCIF) and Ki were compared to those derived from arterial blood samples. We refer to AUC as direct and downstream Ki as indirect comparisons.

Statistical Analysis

All data were reported as mean ± SE. Bland–Altman analysis was used to assess the agreement between experimental and computed cerebral Ki for 35 VOIs. Computed and measured AUC and downstream Ki of the four super regions were compared using paired Student's t-tests. A p-value of <0.05 was considered statistically significant.

Results

In Figure 1, we show arterial blood sampling data in WKY rats. A sketch of the experimental arterial catheterization setup is shown in Figure 1A. In Figure 1B, example picture of carotid artery catheterization is shown in a WKY rat. Average activity (with standard error) normalized by the injected dose and the weight of the animal for five WKY rats are shown in Figure 1C.

FIGURE 1
www.frontiersin.org

Figure 1. Arterial Blood Sampling in Wistar–Kyoto (WKY) rats. (A) Sketch showing catheterization of the carotid artery. (B) Picture of a WKY rat during the experimental procedure. (C) Arterial Blood Sample measurements averaged over n = 5 WKY rats with standard error as a function of image acquisition time.

Example dynamic FDG PET images co-registered with CT over a period of 60 min are shown in Figures 2A–E. The images show that FDG traverses from the vascular space at the early time points (Figures 2A,B) at ~1 min, to the extravascular spaces (Figure 2C) at ~5 min, and gets trapped as FDG-6-phosphatase in the tissue space including the heart and the brain at the late time points (Figures 2D,E) between 30 and 56 min post FDG administration.

FIGURE 2
www.frontiersin.org

Figure 2. Time-resolved total-body 2-[18F] fluoro-2-deoxy-D-glucose positron emission tomography (FDG PET) images. (A–E) Example PET/CT rendering cut across showing uptake in Wistar–Kyoto (WKY) rat myocardium and brain over a period of 60 min.

Initial runs of the optimization using just the first objective term, O1(p) (Eq. 1), showed fits with relatively poor matching of the peak values at the early time points for both tissue (myocardium) and blood data. We hypothesized that this may have been due to the fact that the points are all weighted equally under O1(p), which meant that the matching of peak values may have been sacrificed by the algorithm in favor of later time point values, which are more numerous. Knowing that first-pass behavior is a crucial component of dynamic tracer analysis, we introduced a second objective function, O2(p) (Eq. 2), to create models that more accurately reflect the behavior of this key system component. By introducing O2(p), we weigh the early time points more heavily than others, which gives more equal distribution of emphasis during optimization. This resulted in more accurate fitting of peak values, without significantly lowering the late time point fits as illustrated in Figure 3A (peak blood fit) and Figure 3B (peak tissue fit). Since the blood peaks at the earlier time points (decays with time), forcing the maximum of the model peak to the maximum of the image-derived peak results in more accurate fits at the same time point. However, since the myocardium TACs derived from the images have a lower peak uptake (SP contamination from the blood to the myocardium) at the early time points and higher uptake 10–15 min and beyond (trapping of FDG into the myocardium), forcing the maximum of the Model_myo to the maximum of the image-derived values within the first 10 min results in a slight mismatch in the peak fits for the myocardium.

FIGURE 3
www.frontiersin.org

Figure 3. Peak fits with objective function O2(p). (A) Peak blood fit with and without the second objective function, O2(p) Equation (2). (B) Peak tissue fit with and without O2(p).

In Figures 4A–C, we show example ROIs and model calculations to compute MCIF. Example ROIs in the myocardium and the LVBP at the last time frame are shown in Figures 4A,B. TACs for the myocardium and LV blood pool obtained from dynamic FDG PET images of the heart and model calculations in a dual-output model to compute MCIF along with arterial blood samples for the entire imaging duration of 60 min are shown in Figure 4C. The AUC (direct) for MCIF averaged over five WKY rats was computed to be 1,792.65 ± 348.48 kBq min/cc and that for the arterial blood samples was measured to be 1,893.53 ± 436.91 kBq min/cc (p = 0.76). Although the blood inputs were not significantly different, point by point, the MCIF was consistently lower than the experimental blood input by 6.18% on average relative to the experimental curve.

FIGURE 4
www.frontiersin.org

Figure 4. Model-corrected blood input function (MCIF). (A) Representative 2-[18F] fluoro-2-deoxy-D-glucose positron emission tomography (FDG PET) images with regions of interest drawn in the myocardium and (B) left ventricular blood in the last time frame. (C) Dual-output model (Model IDIF, Model Myo) fitted to image-derived blood (PET IDIF) and myocardium (PET Myo) time activity curves (TACs) to generate the MCIF. For comparison, blood samples derived from arterial blood sampling are shown.

In Figure 5A, the W. Schiffer T2 rat brain atlas fused with a coronal slice of a late time point dynamic PET image is shown for an example WKY rat. Bland–Altman analysis was used to compare the two methods, including tests for similarity and model bias. Bland–Altman analysis between computed and measured Ki for the 35 VOIs contributing to the four super regions for five WKY rats (175 data points) revealed a mean difference of 0.0029 min−1 (~13.5%). The precision (standard deviation of differences) was 0.0103 min−1. All differences were inside the lower and upper limits of agreement (mean ± 1.96 SD), which were −0.0172 and 0.0230, respectively (Figure 5B). The plot appeared to show a trend of increasing error with increasing magnitude of Ki, which necessitated testing for bias to ensure there were no inherent biases beyond the level of the insignificantly underestimated blood input function. The 95% confidence interval limits for mean and two standard deviation-based agreement limits revealed that the mean difference was likely due to a systematic bias from MCIF underestimating the blood sample inputs, with the line of equality falling outside the confidence interval for the mean percentage difference; however, 98.8% of measurements remained within the limits of agreement despite this bias. The four super regions (formed by combining the above 35 VOIs), namely, cerebellum, striatum, cortex, and hippocampus are outlined for a WKY rat in Figure 5C. Example TACs obtained from the dynamic FDG PET images of the four super regions, model calculations, and MCIF are shown in Figure 5D. The RC in Eq. (5) is the term (1–TBV). In Figure 5E, we plot RC as a function of volume, where we show that the model computations result in radioactivity recovery between 0.97 and 1 for all the brain structures (58 data points). RCs for the four super regions (combining 35 of the 58 data points) were close to 1 (data not shown). Computed Ki and measured Ki for the four super regions are shown in Figure 5F. Although the underestimation of the blood input led to a small systematic bias factor, leading to slight overestimation of Ki by the model, no differences for any region were statistically significant. Experimental and computed Ki for the cerebellum were 0.024 ± 0.001 and 0.027 ± 0.01 min−1, respectively (p = 0.67); for the hippocampus were 0.025 ± 0.001 and 0.029 ± 0.01 min−1, respectively (p = 0.47); for the frontal cortex were 0.019 ± 0.001 and 0.022±0.004 min−1, respectively (p = 0.60); and for the striatum were 0.027 ± 0.001 and 0.032±0.01 min−1, respectively (p = 0.45), thereby indicating no significant differences.

FIGURE 5
www.frontiersin.org

Figure 5. Computed cerebral 2-[18F] fluoro-2-deoxy-D-glucose (FDG) uptake rates in Wistar–Kyoto (WKY) rats. (A) Example brain PET image of WKY rat co-registered onto W. Schiffer T2 rat brain atlas showing 35 volumes of interest. (B) Bland–Altman plot comparing computed model and experimental Ki. (C) Example brain PET image of WKY rat co-registered onto W. Schiffer T2 rat brain atlas combining the 35 volumes of interest (VOIs) into four regions of interest. (D) Four-parameter compartment model fitted to time activity curves (TACs) derived from the dynamic PET data for the regions indicated in panel C. (E) Recovery coefficient (RC) of all the 58 VOIs derived from the atlas for an example WKY rat. (F) Comparison between experimental and computed Ki for the cerebellum, striatum, hippocampus, and frontal cortex.

Discussion

Arterial blood sampling, although invasive in nature, is the gold standard for obtaining the blood input function to compute myocardial or cerebral FDG uptake rates in dynamic FDG PET scans of rodents (8). Several image-derived methods have been developed over the last several years for noninvasive determination of the blood input function. One such method is the hybrid method (9) that relies on image-derived sampling from the LVBP at the early time points due to the rapid change in FDG metabolism and invasive arterial or venous blood sampling at the late time points. Another method (1012) employed factor analysis techniques using principle component analysis (PCA) to segment out LVBP from the myocardium and the surrounding background to compute image-derived blood input function based on time-dependent pixels in dynamic FDG PET scans. However, the former method overestimates and the latter method underestimates the computed myocardial Ki due to significant SP-out of radioactivity from the blood to the tissue at the early time points and SP-in from the tissue to the blood at the late time points, respectively. Recent work by others and from our laboratory utilized the IVC as a possible source of deriving the blood input function for computing cerebral (13) and myocardial (3) Ki. Both these methods require either arterial blood sampling for dispersion correction (13) or structural information for PV corrections (3).

Recently, a 15-parameter dual output model was developed in a three-compartment model that simultaneously estimates the MCIF derived from the LVBP with SP and PV corrections and the kinetic parameters for determination of myocardial and cerebral Ki from dynamic FDG PET images of mouse and rat brain and heart in vivo (1, 2, 4). The optimization method here used “fmincon” in Matlab, which utilizes a deterministic algorithm based on interior-reflective Newton (IRN) algorithm to search for a local minimum. One method (1) utilizes a blood sample in the optimization of their objective function, while our methods optimize the same objective function (Eq. 1) on high-resolution ordered subset expectation maximization (OSEM) dynamic gated images of mouse hearts (2) and high-resolution MLEM dynamic images of rat hearts (3, 4), which possibly obviates the need for any blood samples in the optimization. In this study, we used this and further minimized the differences between the model and image-derived peak values (Eq. 2) to estimate MCIF without the need for any arterial blood samples as optimization parameters. This MCIF was further utilized to compute cerebral FDG Ki and also compared to that derived using arterial blood samples in WKY rats. The multiparameter nature of the model may introduce uncertainties in the computed results using nonlinear regression. The initial guesses and bounds for the PV averaging coefficients (Eqs. 4 and 5) were determined from the phantom experiments (RC plot) performed previously in our laboratory (2, 3). The dimensions of the myocardial wall and the left ventricle were based on our recent MRI measurements in WKY rats (4). We kept the bounds for SP factors (Eqs. 4 and 5) open between 0 and 1, as they are a function of time. As for the model blood input function, one side of the bounds was determined by the distribution of the input function and the other side was left wide open as described (2, 14).

It has been recently shown that for structures ranging between 4 and 5 mm and beyond the RC for six MLEM iterations is between 0.6 and 0.85. The RC plateaus out between 0.8 and 0.95 for similar structures at 10 iterations (15). Additionally, we perform partial volume (1-TBV) and spillover (TBV) corrections during modeling as shown in Eq. (5), resulting in radioactivity recovery between 0.97 and 1 for all the brain structures (Figure 5E). Since the brain structures in the study are >2 times the full width half maximum (FWHM) spatial resolution of the scanner, we expect 100% recovery of the radioactivity. However, due to minimal TBV (SP contamination in the tissue), the radioactivity recoveries (1-TBV) for all the brain structures are close to 100% (97% and greater). This model however assumes that the SP and PV (TBV and 1-TBV) add up to 1. A future computation will consider these as independent parameters. We also find for example WKY rat data that the downstream Ki with SP and PV corrections for an image reconstructed with six iterations vs. 24 iterations are very similar (data not shown), indicating stability of the computed values.

A recent review article (16) has identified WKY and Sprague–Dawley (SD) rats as commonly used strains in both cardiac PET and SPECT imaging studies. Previous studies have computed and measured cerebral Ki in SD rats (1, 11). Our computed and measured cerebral Ki in WKY rats are higher than those from SD rats in one case (1) and lower than those in the other (11), although both these studies performed their arterial blood sampling and imaging under isoflurane anesthesia under fed conditions (unfasted rats). Anesthesia plays an important role in tissue FDG uptake in a rodent PET scan especially in the myocardium (four-fold higher uptake for isoflurane vs. ketamine/xylazine). The effect of anesthesia on rodent cerebral FDG uptake, however, is only ~1.3-fold higher using isoflurane compared to that using ketamine/xylazine (17). Fasting duration also affects tissue uptake in an FDG PET scan (18). Whereas, our studies fasted the WKY rats for ~5–6 h, the other studies in SD rats were in fed state (unfasted). Additionally, one of the studies performed attenuation correction using a Co-57 transmission scan (1) compared to CT-based attenuation correction in our study, whereas the other study did not perform any attenuation correction (11). Thus, the differences in the measured cerebral Ki between the earlier studies and ours may possibly be attributed to a combination of the differences in the mode of attenuation correction (transmission based vs. CT based vs. no attenuation correction), anesthetic used (isoflurane vs. Inactin), fasting state (fed vs. unfed), and strain of the rats (WKY vs. SD). A limitation of this study is the sample size. Since the studies were performed in control WKY rats, we do not expect much variation in the metabolic changes in these rats. However, it should be noted that this study is a validation of the methodology for computing the rate of cerebral FDG uptake using input function from the LVBP with SP and PV corrections for the first time in vivo.

Conclusions

Direct and indirect comparisons of model computations with arterial blood sampling were performed for the first time in WKY rats. MCIF computed from the LVBP compared well with arterial blood samples. Computed cerebellum, frontal cortex, hippocampus, and striatal FDG uptake rates agreed well with those obtained from arterial blood samples.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

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

Author Contributions

JM performed modeling of image data. BK conceptualized the study, performed PET imaging, and modeling of image data. BK is the corresponding author. JL performed PET imaging and arterial blood sampling. CM performed image reconstruction. KM performed PET imaging and arterial blood sampling. SP performed PET imaging. VS performed modeling of image data. All authors contributed to the article and approved the submitted version.

Funding

The research in this study was supported in part by NIH grants R01 HL123627 and Radiology start-up funds at UVA (to BK).

Conflict of Interest

CM was employed by company Bruker Biospin.

The remaining 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

We thank Jeremy Gatesman in the Center for Comparative Medicine at the University of Virginia for tail vein catheterizations during PET imaging studies. Imaging data were acquired through the University of Virginia Molecular Imaging Core Lab, with NIH S10OD021672 funding for the Albira Si trimodal scanner. 18F-FDG for imaging was produced through the University of Virginia Radiochemistry Core Lab, with NIH high end funding for the Siemens Eclipse™ HP 11 MeV Cyclotron. FDG activity in blood samples were measured through the University of Virginia Radiochemistry Core Lab housing the Hidex Automatic γ-Counter.

References

1. Fang YH, Muzic RF Jr. Spillover and partial-volume correction for image-derived input functions for small-animal 18F-FDG PET studies. J Nucl Med. (2008). 49:606–14. doi: 10.2967/jnumed.107.047613

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Zhong M, Kundu BK. Optimization of a model corrected blood input function from dynamic FDG-PET images of small animal heart in vivo. IEEE Trans Nucl Sci. (2013) 60:3417–22. doi: 10.1109/TNS.2013.2269032

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Huang Q, Massey JC, Mińczuk K, Li J, Kundu BK. Non-invasive determination of blood input function to compute rate of myocardial glucose uptake from dynamic FDG PET images of rat heart in vivo: comparative study between the inferior vena cava and the left ventricular blood pool with spill over and partial volume corrections. Phys Med Biol. (2019) 64:165010. doi: 10.1088/1361-6560/ab3238

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Li J, Kemp B, Howell N, Massey J, Mińczuk K, Huang Q, et al. Metabolic changes in spontaneously hypertensive rat hearts precede cardiac dysfunction, and left ventricular hypertrophy. J Am Heart Assoc. (2019) 8:e010926. doi: 10.1161/JAHA.118.010926

PubMed Abstract | CrossRef Full Text | Google Scholar

5. González AJ, Aguilar A, Conde P, Hernández L, Moliner L, Vidal LF, et al. A PET design based on SiPM and monolithic LYSO crystals: performance evaluation. IEEE Trans Nucl Sci. (2016) 63:2471–7. doi: 10.1109/TNS.2016.2522179

CrossRef Full Text | Google Scholar

6. Gsell W, Molinos C, Correcher C, Belderbos S, Wouters J, Junge S, et al. Characterization of a preclinical PET insert in a 7 Tesla MRI scanner: beyond NEMA testing. Phys Med Biol. (2020) 65:245016. doi: 10.1088/1361-6560/aba08c

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Feng D, Huang SC, Wang X. Models for computer simulation studies of input functions for tracer kinetic modeling with positron emission tomography. Int J Biomed Comp. (1993) 32:95–110. doi: 10.1016/0020-7101(93)90049-C

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Laforest R, Sharp TL, Engelbach JA, Fettig NM, Herrero P, Kim J, et al. Measurement of input functions in rodents: challenges and solutions. Nucl Med Biol. (2005) 32:679–85. doi: 10.1016/j.nucmedbio.2005.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Shoghi KI, Welch MJ. Hybrid image and blood sampling input function for quantification of small animal dynamic PET data. Nucl Med Biol. (2007) 34:989–94. doi: 10.1016/j.nucmedbio.2007.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Kim J, Herrero P, Sharp T, Laforest R, Rowland DJ, Tai YC, et al. Minimally invasive method of determining blood input function from PET images in rodents. J Nucl Med. (2006). 47:330–6.

PubMed Abstract | Google Scholar

11. Fakhri GE, Sitek A, Guerin B, Kijewski MF, Di Carli MF, Moore SC. Quantitative dynamic cardiac 82Rb PET using generalized factor and compartment analyses. J Nucl Med. (2005) 46:1264–71.

PubMed Abstract | Google Scholar

12. Huang Q, Li J, Taegtmeyer H, Kundu BK. Blood input function from dynamic FDG PET images of Wistar Kyoto rat heart to compute rate of myocardial glucose utilization in vivo: comparison between region of interest and factor analysis methods. J Nucl Med Abstract. (2018) 59:367.

Google Scholar

13. Lanz B, Poitry-Yamate C, Gruetter R. Image-derived input function from the vena cava for 18F-FDG PET studies in rats and mice. J Nucl Med. (2014) 55:1380–8. doi: 10.2967/jnumed.113.127381

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhong M, Alonso CE, Taegtmeyer H, Kundu BK. Quantitative PET imaging detects early metabolic remodeling in a mouse model of pressure-overload left ventricular hypertrophy in vivo. J Nucl Med. (2013) 54:609–15. doi: 10.2967/jnumed.112.108092

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Attarwala AA, Karanja YW, Hardiansyah D, Romanó C, Roscher M, Wängler B, et al. Investigation of the imaging characteristics of the ALBIRA II small animal PET system for 18F, 68Ga and 64Cu. Z Med Phys. (2017) 27:132–44. doi: 10.1016/j.zemedi.2017.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Cicone F, Viertl D, Quintela Pousa AM, Denoël T, Gnesin S, Scopinaro F, et al. Cardiac radionuclide imaging in rodents: a review of methods, results and factors at play. Front Med. (2017) 4:35. doi: 10.3389/fmed.2017.00035

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Toyama H, Ichise M, Liow J-S, Vines DC, Seneca NM, Modell KJ, et al. Evaluation of anesthesia effects on [18F]FDG uptake in mouse brain and heart using small animal PET. Nucl Med Biol. (2004) 31:251–6. doi: 10.1016/S0969-8051(03)00124-0

CrossRef Full Text | Google Scholar

18. Lee KH, Ko BH, Paik JY, Jung KH, Choe YS, Choi Y, et al. Effects of anesthetic agents and fasting duration on 18F-FDG biodistribution and insulin levels in tumor-bearing mice. J Nucl Med. (2005) 46:1531–6.

PubMed Abstract | Google Scholar

Keywords: arterial blood sampling, dual output model, cerebral fluoro-2-deoxy-D-glucose uptake rate, Wistar–Kyoto rat, dynamic fluoro-2-deoxy-D-glucose positron emission tomography

Citation: Massey JC, Seshadri V, Paul S, Mińczuk K, Molinos C, Li J and Kundu BK (2021) Model Corrected Blood Input Function to Compute Cerebral FDG Uptake Rates From Dynamic Total-Body PET Images of Rats in vivo. Front. Med. 8:618645. doi: 10.3389/fmed.2021.618645

Received: 17 October 2020; Accepted: 15 February 2021;
Published: 07 April 2021.

Edited by:

Francesco Cicone, University of Catanzaro, Italy

Reviewed by:

David John Macfarlane, The University of Queensland, Australia
Jing Wu, Beijing Normal University, China

Copyright © 2021 Massey, Seshadri, Paul, Mińczuk, Molinos, Li and Kundu. 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: Bijoy K. Kundu, bkk5a@virginia.edu

This author has contributed equally to this work

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.