In Vivo Imaging of Inflammation and InfectionView this Special Issue
Research Article | Open Access
Lars Jødal, Svend B. Jensen, Ole L. Nielsen, Pia Afzelius, Per Borghammer, Aage K. O. Alstrup, Søren B. Hansen, "Kinetic Modelling of Infection Tracers [18F]FDG, [68Ga]Ga-Citrate, [11C]Methionine, and [11C]Donepezil in a Porcine Osteomyelitis Model", Contrast Media & Molecular Imaging, vol. 2017, Article ID 9256858, 18 pages, 2017. https://doi.org/10.1155/2017/9256858
Kinetic Modelling of Infection Tracers [18F]FDG, [68Ga]Ga-Citrate, [11C]Methionine, and [11C]Donepezil in a Porcine Osteomyelitis Model
Introduction. Positron emission tomography (PET) is increasingly applied for infection imaging using [18F]FDG as tracer, but uptake is unspecific. The present study compares the kinetics of [18F]FDG and three other PET tracers with relevance for infection imaging. Methods. A juvenile porcine osteomyelitis model was used. Eleven pigs underwent PET/CT with 60-minute dynamic PET imaging of [18F]FDG, [68Ga]Ga-citrate, [11C]methionine, and/or [11C]donepezil, along with blood sampling. For infectious lesions, kinetic modelling with one- and two-tissue-compartment models was conducted for each tracer. Results. Irreversible uptake was found for [18F]FDG and [68Ga]Ga-citrate; reversible uptake was found for [11C]methionine (two-tissue model) and [11C]donepezil (one-tissue model). The uptake rate for [68Ga]Ga-citrate was slow and diffusion-limited. For the other tracers, the uptake rate was primarily determined by perfusion (flow-limited uptake). Net uptake rate for [18F]FDG and distribution volume for [11C]methionine were significantly higher for infectious lesions than for correspondingly noninfected tissue. For [11C]donepezil in pigs, labelled metabolite products appeared to be important for the analysis. Conclusions. The kinetics of the four studied tracers in infection was characterized. For clinical applications, [18F]FDG remains the first-choice PET tracer. [11C]methionine may have a potential for detecting soft tissue infections. [68Ga]Ga-citrate and [11C]donepezil were not found useful for imaging of osteomyelitis.
Positron emission tomography (PET) allows imaging of molecular uptake but is dependent on the availability of tracers with uptake related to the investigated disease. For infection imaging, [18F]FDG is useful , but a drawback is the nonspecific uptake of [18F]FDG in metabolically active tissues (brain, muscles, etc.) .
We previously reported on static PET imaging in a porcine osteomyelitis model with a series of tracers that have been proposed for infection imaging [3–5]. In those reports, we concluded that the performance of [18F]FDG for locating infection was superior to the other tracers studied. However, our static imaging shows only the uptake at a given (typically late) time interval.
In the present paper, we elaborate on those studies [3–5] by including data from dynamic PET imaging of the same animals (based on the same tracer injections), using kinetic analysis to compare [18F]FDG with [68Ga]Ga-citrate, [11C]methionine (L-[11C-methyl]methionine), and [11C]donepezil ([5-11C-methoxy]donepezil) as infection tracers. Dynamic imaging allows for a more detailed study of the uptake and release of tracers, which may contribute to a better understanding of the underlying physiology and also help in determining the optimal time for static imaging. Blood perfusion results from dynamic imaging of [15O]water in the same animals have already been reported .
Gamma-camera imaging (scintigraphy or SPECT) of infection is sometimes performed with [67Ga]Ga-citrate. However, image quality is suboptimal and the radiation dose from 67Ga is high. The positron-emitter 68Ga allows PET imaging, which has a better spatial resolution and a higher detection efficiency than standard gamma-camera imaging. Furthermore, the shorter half-life of 68Ga (67.7 min versus 78 hours for 67Ga) reduces the radiation dose. In the body, the Ga3+ ion acts as an analogue of iron (Fe3+), associating with transferrin; for more details, see . Preliminary studies on [68Ga]Ga-citrate in rats  and patients [9, 10] reported promising results.
Methionine is a naturally occurring essential amino acid. It can be labelled with the positron-emitter 11C to obtain the chemically identical PET tracer [11C]methionine. Methionine (natural or labelled) is transported into cells via the L-type amino acid transporter 1. It is crucial for the formation of proteins and is involved in the synthesis of phospholipids. The accumulation of methionine therefore reflects amino acid transport and protein synthesis. During cell replication, the demand for essential amino acids increases, as do protein and phospholipid synthesis. Although primarily used to detect malignant tumours , methionine is known from brain studies to accumulate in inflammatory tissue [12, 13]. Furthermore, a high uptake of [11C]methionine has been observed after acute myocardial infarction, indicating that this tracer may be a good marker of inflammatory reactions during the phase of tissue repair .
Donepezil is a reversible inhibitor of acetylcholinesterase (AChE), the break-down enzyme of acetylcholine (ACh) . ACh and AChE are primarily known as constituents in cholinergic neuronal signalling pathways, but immune cells have been shown to use ACh as a paracrine signalling molecule , and AChE and cholinergic receptors are upregulated in immune cells when exposed to phytohaemagglutinin . It was recently shown that the AChE ligand [11C]donepezil accumulates at sites of bacterial infection, suggesting that cholinergic PET imaging may have potential in the detection of inflammation and infections .
Therefore, we investigated the uptake processes of these PET tracers in order to determine their clinical potential in infection imaging in relation to both osteomyelitis and soft tissue infection. We hypothesized that categorization of tracers into those having either reversible or irreversible uptake in infectious lesions is possible and perfusion effects on tracer accumulation can be estimated. Finally, we wanted to estimate the optimal time points for static imaging of the tracers.
2. Materials and Methods
2.1. Porcine Osteomyelitis Protocol
The protocol for inducing osteomyelitis in domestic pigs has been described in detail elsewhere [3, 19–21]. Briefly, osteomyelitis was induced in the right hind limb of juvenile Danish Landrace × Yorkshire female pigs by intra-arterial injection of S. aureus (porcine strain S54F9). The injection was performed into the right femoral artery, allowing for selective infection of this hind limb while leaving the contralateral left hind limb as a noninfected control. Osteomyelitis was allowed to develop for one week, after which the pig was scanned and then euthanized. If a pig reached predefined humane endpoints, it was then euthanized (and not scanned). A refinement of the model, combining reduced body weight of the pigs (from 40 kg to 20 kg) with administration of penicillin at the onset of the first clinical signs of disease, proved effective in reducing systemic infection .
We attempted to prolong the infection period from one to two weeks, for the purpose of developing more chronic lesions, and managed to scan one pig two weeks after inoculation (pig number 5 of the 11 scanned pigs described below). However, 4 out of 5 pigs planned for scanning with the prolonged protocol had to be prematurely euthanized due to humane endpoints, so we returned to the subacute one-week pig protocol.
The animal protocol was approved by the Danish Animal Experimental Board, journal number 2012-15-2934-00123, and all procedures followed the European Directive 2010/63/EU on the protection of animals used for scientific purposes.
2.2. Animals and Lesions
Eleven juvenile female domestic pigs were scanned. Pigs number 1–4 had body weights of 39–42 kg (initial model) and pigs number 5–11 had body weights of 19–23 kg (refined model). Blood perfusion in the lesions in these pigs has been described previously , including 17 osteomyelitic (OM) lesions and 8 soft tissue (ST) lesions, with volume of interest (VOI) drawing primarily based on computed tomography (CT) scans. Other OM lesions were also found in pedal bones but were too small for robust volume of interest (VOI) drawing.
In the present paper, the kinetics of [18F]FDG, [68Ga]Ga-citrate, [11C]methionine, and [11C]donepezil were studied and modelled in the same 17 OM lesions and 8 ST lesions, with VOIs redrawn on the CTs from the dynamic PET/CT scans of these tracers. Not all tracers were used in all pigs, and the limitations set by the PET scanner field of view (FOV, further described below) resulted in not all lesions being dynamically scanned in all pigs. An overview is given in Table 1.
| as presented in the 15Owater paper . 15Owater paper presents an exhaustive list of the lesions studied. In pig number 3, inoculation failed to produce infection . soft tissue lesions within the FOV (pig number 4).|
We have previously reported some data on the characterization of the lesions evolving from S. aureus inoculation in these pigs, including the methods used for bacteriological and immunohistochemical (IHC) identification of bacteria [3–5]. Generally, the lesions observed in the pigs were considered to be caused by infection associated with the inoculated S. aureus if the bacterial culture and/or S. aureus specific IHC staining confirmed the presence of S. aureus in one or several of the lesions within each individual pig (results not shown).
2.3. Dynamic PET Scans
PET scans were carried out at the Department of Nuclear Medicine & PET Centre (Aarhus University Hospital), and at Department of Nuclear Medicine (Aalborg University Hospital), with the pigs transported from one hospital to the other. The PET scans at Aarhus included the 11C-labelled tracers in pigs number 1–10 and [18F]FDG in pig number 11. The PET scans at Aalborg included [68Ga]Ga-citrate in pigs number 1–5 and [18F]FDG in pigs number 1–10.
At Aarhus, the PET data were acquired on a Biograph TruePoint 64 PET/CT scanner (Siemens, Erlangen, Germany). The scan field covered 21 cm in the axial direction and was positioned over the pelvic region and the hind limbs. The images were reconstructed with an OSEM algorithm with resolution recovery (TrueX, Siemens). The reconstruction parameters were 6 iterations, 21 subsets, 336 × 336 matrix in 109 slices, voxel size 2 × 2 × 2 mm3, and a 2 mm Gaussian filter. The spatial resolution of the reconstructed images was approximately 4 mm.
At Aalborg, the PET data were acquired on a GE VCT Discovery 64 PET/CT scanner (GE Healthcare, USA). The scan field covered 15 cm in the axial direction and was positioned over the pelvic region and the hind limbs. The images were reconstructed with an OSEM algorithm without resolution recovery (3D Vue Point, GE). The reconstruction parameters were 2 iterations, 28 subsets, 128 × 128 matrix in 47 slices, voxel size 5.5 × 5.5 × 3.3 mm3, and a 6 mm Gaussian filter.
On both scanners, image reconstruction included decay-correction to the start of scanning and attenuation-correction based on CT scanning.
For all tracers, the pigs were dynamically PET scanned for 60 minutes in 23 frames: 8 × 15 s, 4 × 30 s, 2 × 60 s, 2 × 120 s, 4 × 300 s, and 3 × 600 s. For [68Ga]Ga-citrate, the animals were scanned for an additional 6 × 600 s (i.e., 120 minutes total scan time in 29 frames). After each of the dynamic scans, the pigs were statically scanned; the static PET/CT scans have been described in previous papers [3–5] and will not be further discussed here.
The pigs were scanned in dorsal recumbency (supine position). The hind limbs were positioned for the entire hind limbs and pelvis to be within the axial field of view (FOV) of the scanner. This was, however, not always possible; in particular, the 15 cm FOV at the Aalborg scanner (versus 21 cm FOV at the Aarhus scanner) was a limitation. To optimize the fixation position to the scanner FOV, custom-made fixation devices were used for pigs number 6–11 (Figure 1).
2.4. Blood Samples
Blood samples were drawn from the carotid artery at predetermined time points as listed below. The samples were manually drawn, and small variations occurred; the actual time of each sample was recorded. The samples were centrifuged to obtain plasma samples, which were counted in calibrated gamma-counters.
At Aarhus, plasma samples were counted in a Packard Cobra gamma counter. An energy window from 400 to 1400 keV was used. No signs of interference between consecutively administered radionuclides were seen in the plasma curves.
At Aalborg, plasma samples and full-blood samples were counted in a Wizard 2480 gamma counter (PerkinElmer, Turku, Finland). To avoid interference from remnants of 111In from 111In-leukocytes (used in the same animals ), an energy window from 450 to 1200 keV was used .
In the 40 kg pigs (pigs number 1–4), 37 blood samples were drawn per tracer: every 5 seconds for 1 minute (12 samples), at 70, 80, 90, 100, 120, 140, 160, 180, 210, 240, 270, and 300 seconds (12 samples), and at 6, 7, 8, 9, 10, 15, 20, 25, 30, 35, 40, 50, and 60 minutes (13 samples). For [68Ga]Ga-citrate, blood samples were also drawn at 75, 90, 105, and 120 minutes (four samples), for a total of 41 blood samples.
In the 20 kg pigs (pigs number 5–11), the number of blood samples was reduced to 26 samples per tracer: every 5 seconds for 50 seconds (10 samples), at 60, 80, 100, 120, 150, 180, 240, and 300 seconds (8 samples), and at 6, 8, 10, 15, 20, 30, 40, and 55 minutes (8 samples).
For [11C]methionine and [11C]donepezil, additional blood samples were drawn for metabolite analysis: at 2, 5, 10, 15, 25, and 40 minutes in the 40 kg pigs, and at 2, 5, 10, 20, 30, 40, and 55 minutes in the 20 kg pigs.
2.5. Input Function
The blood plasma rather than full blood was considered to be the reference fluid for tracer delivery (corresponding to equilibration between blood cells and plasma being slow compared to single-passage time of blood in the tissue). Accordingly, all input functions were based on plasma samples.
For both [18F]FDG and [68Ga]Ga-citrate, the decay-corrected plasma sample data were used as the input function. No metabolite correction was performed for these two tracers because metabolite products are not expected to be found in the blood. [18F]FDG is phosphorylated within the cells, but the resulting (radioactive) metabolite is trapped within the cells . Regarding [68Ga]Ga-citrate, the Ga-citrate complex quickly dissociates into Ga3+ and citrate3- within the blood, but the gallium ion attaches to transferrin (making [68Ga]Ga-transferrin the actual tracer). Therefore, free gallium is not found in the blood .
For [11C]methionine and [11C]donepezil, metabolite correction was performed. A fractionated HPLC analysis was used to separate the metabolites from the parent tracer, and counting was used to determine the fraction of activity representing the parent tracer. Based on the obtained data points, Hill-type fraction curves were fitted:where is the sampling time (seconds postinjection). The function starts at (thus assuming no metabolism before injection) and has an asymptotic value . The parameters , , and were fitted for both tracers in the individual pigs.
In the following text, the uncorrected input function will denote the activity concentration (decay-corrected Bq/mL) from plasma samples and the metabolite-corrected input function will denote times the uncorrected input function.
2.6. Determination of the Delay-Correction
The measured input function can be biased by delay and dispersion effects due to the differences in distance between the blood sampling site (carotid artery) and the infection sites in the hind limbs. Such a delay can be determined by applying a series of possible delay-correction values, fitting the data with each value, and selecting the delay-correction resulting in the best fit .
We applied this procedure for each tracer in each animal, using input function offsets from −60 s to +60 s in 1-second steps. This was done using the uncorrected plasma data for the input function, and fitting a reversible two-tissue compartment model (rev2TCM in Figure 2; the models are further discussed below) to the first 300 s of the full field-of-view data. Using full FOV data ensures that the statistical noise is low. Restricting this fit to the early data has several advantages. First, this focuses on the part of the study where the input function changes fast and the delay therefore is important. Second, this means that physiological differences in uptake (e.g., bladder versus nonbladder, infection versus noninfection) will only have had little time to manifest; thus the full FOV data will be dominated by the bolus passage in this anatomical part of the animal, rather than by a mix of physiologies. Third, the possible metabolism effects will not yet dominate the input function, for which reason the uncorrected input function can be considered representative, even if the tracer over time is metabolized.
2.7. Kinetic Models
Relatively little literature exists on tracer kinetics in infections for the tracers investigated in this work. Rather than imposing a specific model on these tracers, we examined the applicability of three different models for each tracer (Figure 2).
Physiologically, the 1TCM corresponds to the tracer entering and leaving the tissue with no binding or other specific uptake. The irr2TCM describes uptake and irreversible trapping of the tracer or metabolites (e.g., phosphorylated [18F]FDG). The rev2TCM corresponds to uptake in the tissue followed by reversible binding (or reversible metabolism) of the tracer.
In all of the models, the rate constant describes the first-pass uptake of tracer, equal to the product of perfusion and the first-pass extraction fraction. We report in units of mL/min/100 cm3, that is, mL uptake per minute per 100 cm3 of tissue. The other rate constants (, , and , all with unit min−1) describe how rapid the concentration in a compartment changes due to a given process (excretion, binding, and metabolism in the tissue). The blood fraction is the fraction of measured PET signal that originates from the blood in the vascular bed.
For models with irreversible uptake, the net uptake rate (same unit as ) represents the effective irreversible uptake from the input. Whereas represents the rate of immediate (first-pass) uptake of tracer, can similarly be interpreted as the rate of long-term uptake. For the irr2TCM, the theoretical net uptake rate isThe net uptake rate may also be determined as the slope of a Patlak plot [26, 27].
For models with reversible uptake, the distribution volume (DV, unit mL/cm3) is the ratio of the tissue concentration to the “input” concentration, once a steady-state has been reached, that is, the volume of “input” needed to account for the activity in 1 cm3 of tissue . Although not being an uptake rate, DV may be used as an indicator of the degree of long-term uptake. For the 1TCM and rev2TCM, the theoretical distribution volumes areThe distribution volume of reversible uptake may also be determined as the slope of a Logan plot [29, 30].
2.8. Weighting of the PET Data in Modelling
Theoretically, least-squares fitting is optimal with weights proportional to where is the variance of the noise. With counts (Poisson distributed) during a frame length , the count rate has variance:After decay-correction:where is the decay-correction factor, calculated from the radionuclide half-life and frame time interval. For weighting purposes, it can be well approximated by using the mid-time of the frame:Seemingly, optimal weighting should be for non-decay-corrected data and ) for decay-corrected data.
However, the count rate (or the activity concentration) is known only from a measurement that includes noise. In a simulation study, Thiele and Buchert  found that noise in the weighting factors can severely degrade parameter estimation and should therefore be avoided. Consistent with the weighting that gave the best results in that study, we used the following noise-free weighting factors for the decay-corrected PET data:These weights correctly include the effects of decay and the large differences in frame length (from 15 to 600 seconds). The weights are noise-free, at the cost of ignoring the statistical effects from variation in tracer concentration due to kinetics.
This approximation (weights not fully reflecting 1/σ2) is unlikely to be a problem. A simulation study by Yaqub et al.  found that kinetic modelling of PET data was reasonably robust against some misrepresentation of the variance in (noise-free) weighting, with only severe misrepresentation being a problem.
For each tracer in each model (1TCM, irr2TCM, and rev2TCM), the fitted parameters were determined using least-squares fitting with the weighting described above. Additionally, Patlak plots and Logan plots were computed, based on the data from 10 minutes postinjection (p.i.) and onwards.
For [11C]methionine and [11C]donepezil, this procedure was performed twice: using an uncorrected input function and using a metabolite-corrected input function. For these two tracers, modelling was restricted to data from the first 40 minutes (out of 60 minutes), as metabolite data were in many cases incomplete for later frames.
Modelling was performed using software acquired from the Turku PET Centre website . The parameter was fitted as the ratio / (corresponding to the distribution volume of the first compartment). In the rev2TCM, the parameter was fitted as the ratio / (corresponding to the binding potential BP if the first and second compartments are considered to represent unspecific and specific uptake, resp.).
In addition to visual inspection of the fits, the three models were compared with the corrected Akaike Information Criterion (), which rewards a good fit but punishes the use of a model with many fitting parameters. For a given data set, favours the model resulting in the lowest value [34, 35].
For the determination of uptake as reversible or irreversible, the Patlak plot also was considered. If the uptake is reversible (i.e., not irreversible), the Patlak plot will eventually approach a constant value. Therefore, the linearity of the Patlak plot with a nonzero slope can be used as a test for irreversible uptake.
All these results are based on plasma input functions. To compare with perfusion, previously published results for blood perfusion in the same animals  were transformed into plasma perfusion by the formulawith haematocrit being measured at a time roughly corresponding to the time of the [15O]water PET.
3.1. Metabolite Correction
For both [11C]methionine and [11C]donepezil, the “uncorrected” (i.e., not corrected for metabolism, but corrected for physical decay) plasma curves showed an unexpected tendency of slightly rising values after typically 20 minutes. As an example, see the data for [11C]methionine in Figure 5 (the data for [11C]donepezil were similar). Possible reasons like unresolved background counts, either from other tracers in the multitracer study or from the surroundings, were investigated, but no sign of any such problems was found; for example, blood samples taken before the arrival of the bolus injection were far below the level that would cause this background. Therefore, we conclude that the curves correctly represent the activity concentration in the plasma. After metabolite correction, the curves decreased as expected (Figure 5).
3.2. Delay-Correction of Input Function
Each input function () was individually delay-corrected (but with a common delay for all lesions in the same scan). The mean ± SD of the corrections was seconds. A negative correction corresponds to the tracer arriving earlier to the scanned tissue (PET data) than to the site of blood sampling (plasma input data). All individual corrections were numerically smaller than the initial PET frame length (15 seconds).
The overall results for the different models are summarized in Table 2, with elaborating comments given here.
Despite the lowest values having been found for rev2TCM, the irr2TCM appeared visually to give a reasonable fit for the investigated time range. Example fits are shown in Figure 6. The Patlak plots (not shown) were linear with nonzero slopes, indicating the presence of irreversible uptake. Also, the majority of fitted ratios / were above 2 (median value ~4 for all VOIs, ~5 if restricted to infected side); that is, overall was considerably lower than . For these reasons and because irr2TCM has fewer parameters than rev2TCM, we pragmatically chose to base the further analysis of [18F]FDG uptake on the irr2TCM. The values of were overall very similar for irr2TCM and rev2TCM.
Overall, the values indicated nearly equal quality of fits for irr2TCM and rev2TCM. Example fits are shown in Figure 7. The Patlak plots were linear with nonzero slopes. Further analysis will assume the simplest model: irr2TCM.
For both [18F]FDG and [68Ga]Ga-citrate, good correspondence was observed between the slope of the Patlak plot and calculated from the irr2TCM parameters.
The values sometimes favoured rev2TCM and sometimes favoured irr2TCM, but, from mean the and median values, rev2TCM was favoured. Visually, the irr2TCM fit showed a problematic upward trend in the late part of the fits (Figure 8). Despite the extent of metabolism during the study (Figure 3), only a small difference was observed between using the uncorrected or the metabolite-corrected input function with the rev2TCM: the values were comparable, the fits were visually very similar, and the fitted rate parameters , /, and were very similar. Only the / values differed markedly, being lower for the corrected than the uncorrected input function. However, the ratio is important for the distribution volume (see (4)), and the further analysis of the [11C]methionine data will assume the rev2TCM model with a metabolite-corrected input function.
Both the plots and the values unequivocally favoured fits with a metabolite-corrected input function (Figure 9). Within these, the values variably favoured each of the three models but with no obvious pattern (e.g., not distinguishing infected versus noninfected tissue, bone versus soft tissue, or high versus low ). Visually, however, the 1TCM (with metabolite-corrected input) fits well in all cases. A typical fit for the 1TCM is seen in Figure 9(b). Further analysis of [11C]donepezil will be based on the 1TCM with a metabolite-corrected input function.
3.4. Perfusion and First-Pass Uptake Rate
The parameter represents the product of perfusion and the extraction fraction for the tracer. Thus, plotting as a function of perfusion gives a measure of extraction; if the extraction fraction is independent of perfusion, the plot will show proportionality between and perfusion. The blood perfusion of these lesions in these animals has been described previously, based on [15O]water PET scans . These data were transformed into plasma perfusion according to (9).
Figure 10 shows the values of plotted as a function of plasma perfusion. Notably, [68Ga]Ga-citrate shows only small uptake compared to perfusion, that is, a small extraction fraction.
The paradoxical plot for [11C]donepezil, showing values that are significantly higher than perfusion (corresponding to >100% extraction), is not a result of using the simple 1TCM rather than one of the 2TCM models. Overall, the data from these models (not shown) were very similar, in some cases even higher, resulting in very similar plots (not shown). For further explanations, see Discussion.
For [18F]FDG, [11C]methionine, and [11C]donepezil, a paired t-test showed a higher in the infected (right) side than in the noninfected (left) side (). For [68Ga]Ga-citrate, no significant difference was found ().
3.5. Irreversible Net Uptake Rate () for [18F]FDG and [68Ga]Ga-Citrate
For both [18F]FDG and [68Ga]Ga-citrate, good agreement was observed between calculated from the Patlak plots and from the irr2TCM model (data not shown).
Figure 11 shows as a function of plasma perfusion. Figure 12 compares in the infected versus corresponding noninfected positions. For [18F]FDG, values were significantly higher in the infected lesions than in the corresponding noninfected positions (, paired t-test). For [68Ga]Ga-citrate, the difference was not statistically significant ().
3.6. Distribution Volume (DV) of Reversible Tracers
For [11C]donepezil fitted with the 1TCM, good agreement was observed between DV calculated from the Logan plots and from (3). For [11C]methionine fitted with the rev2TCM, the agreement was not as good between DV from the Logan plots and from (4); however, the lack of agreement was due to the rev2TCM giving unrealistically high values in some cases (e.g., >100 mL/cm3), corresponding to cases with high values of /. Excluding these cases, good agreement was observed in the DV calculations.
Overall, the slope of the Logan plot was used as a robust measure of DV. Figure 13 shows DV as a function of plasma perfusion. Figure 14 compares DV in the infected and corresponding noninfected locations. For [11C]methionine, DV was significantly higher in the infected locations (), while the minor difference seen for [11C]donepezil was not significant ().
Four very different PET tracers were studied for their potential as markers of infection in a porcine model of osteomyelitis (including associated soft tissue lesions). As noted, the results from static imaging have been reported earlier [3–5], but dynamic imaging can provide more information on the uptake processes, which is not available when the tracer concentration is measured at only one time interval.
The modelling of PET data is restricted by the length of the acquisition, the number of data points, and statistical noise. Therefore, some level of pragmatism is needed when setting up or choosing a model. Strictly irreversible uptake (in the sense of the molecules staying in the body for life) is rare, but the efflux level from a compartment may be practically zero relative to the length of the PET acquisition. In the present study, we attempted to find a level that makes the models useful for providing information about the uptake process while accepting a level of pragmatism to distinguish “practically irreversible” from “practically reversible” uptake.
For the kinetic modelling of [18F]FDG, the two classical models are the irr2TCM by Sokoloff et al.  and the rev2TCM by Phelps et al. . In both models, the second tissue compartment represents the metabolite product [18F]FDG-6-phosphate, and in the latter model represents dephosphorylation back to [18F]FDG. Relative to infection, the kinetic modelling of [18F]FDG uptake appears only to have been performed in lung studies (reviewed in ) and in a single study of an acute viral infection . All of these studies focus on irreversible uptake models, generally the irr2TCM, although some lung studies included a separate reversible compartment for uptake in pulmonary oedema [38, 40].
In our porcine osteomyelitis model, the uptake of [18F]FDG was found to be (practically) irreversible, with a reasonable fit by the irr2TCM. Accordingly, the level of uptake was evaluated based on the irreversible net uptake rate, , in almost all cases showing elevated uptake in lesions compared to healthy tissue, for both OM and ST lesions (Figure 12). The correlation with perfusion seen for both first-pass uptake (Figure 10) and net uptake (Figure 11) indicates uptake of [18F]FDG to be more flow-limited than diffusion-limited.
Dynamic PET studies of [68Ga]Ga-citrate and [68Ga]Ga-transferrin have been published before [8, 41, 42], but the present study appears to be the first to include kinetic modelling. As noted by Kumar and Boddeti , no literature is available on the early imaging times of [67Ga]Ga-citrate SPECT. We found the uptake of 68Ga to be well described by the irr2TCM (Figure 7). However, first-pass uptake () of 68Ga was small compared to that of the other tracers, indicating a very small extraction fraction with little dependence on perfusion (Figures 10 and 11), that is, diffusion-limited.
The physiological reason for the slow uptake of 68Ga may be related to the binding of gallium to a large protein (transferrin). Uptake will either require extravasation of the large [68Ga]Ga-transferrin complex or require a two-step process, such as the release of 68Ga from the protein before the uptake or the uptake of [68Ga]Ga-transferrin by leukocytes which then enter tissue . Therefore, even though bacteria may show increased uptake of gallium due to its chemical similarities with iron , the overall uptake mechanism appears to be quite slow.
For static imaging, slow uptake favours late imaging. Physically, the half-life of 68Ga restricts imaging to a few hours after injection. In a study of lung lesions, Vorster et al.  recommended imaging to start no later than 120 minutes p.i.
Compared to other 68Ga infection studies [8–10, 41], the results in the porcine osteomyelitis model (present paper and [3, 4]) appear disappointing. At the basic level, the tracer is a Ga3+ ion (with chemical similarities to the Fe3+ ion), which makes a species difference unlikely. A different reason may be indicated by Figure 12: maybe the tracer is able to differentiate infected soft tissue from healthy tissue (despite the slow uptake rate) but is not suitable for bone infections. Mäkinen et al.  did find uptake in bone lesions in a rat model, but noted as a limitation of the study that their model “perhaps best simulates osteomyelitis arising from grossly contaminated long-bone fractures.” In contrast, the porcine osteomyelitis model represents haematogenous osteomyelitis without bone trauma.
The modelling of [11C]methionine uptake required a rev2TCM. For comparison, Fischman et al.  described the muscle uptake of [11C]methionine with a 2TCM, where the second compartment represented the incorporation of [11C]methionine in tissue proteins, and they assumed that the degradation rate of labelled protein could be ignored (corresponding to = 0 in our notation). That is, they suggested an irr2TCM rather than rev2TCM. Our finding of rev2TCM as preferable thus corresponds to protein degradation being nonnegligible.
Despite the considerable metabolism of [11C]methionine during the acquisition time (Figure 3), practically only the rate constant depended on whether modelling was based on the uncorrected or the metabolite-corrected input function (with higher values in the latter case). It appears that the metabolite products have kinetics quite similar to the original molecule. The correlation with perfusion seen for both first-pass uptake (Figure 10) and distribution volume (Figure 13) indicates uptake of [11C]methionine to be flow-limited.
The unexpected rise in the plasma activity curves after ~20 minutes (Figure 5) might be explained by a heavy uptake by metabolizing organs (the liver), followed by a later release of radioactive metabolite products to the blood pool.
Generally, the distribution volume (DV) for [11C]methionine was higher in the infected tissue than in the noninfected tissue (Figure 14). This difference between infected and noninfected tissue could point to [11C]methionine having a role in infection imaging (regardless of the role of perfusion in causing the distinction). In line with the results from static imaging [3, 5], the distinction between infected and noninfected tissue appeared more clear for soft tissue than for bone (Figure 14). As indicated by the two already mentioned case reports [12, 13], [11C]methionine could be useful for brain infection imaging, where [18F]FDG suffers from the high physiological uptake in healthy brain tissue.
The uptake curves in both infected and noninfected tissues appear quite stable after approximately 10–15 minutes (Figure 8), indicating that static imaging could be performed starting at this time.
Metabolite-corrected [11C]donepezil could be modelled with the 1TCM model. However, the parameter was systematically higher than plasma (and blood) perfusion, paradoxically indicating an extraction fraction above 100% (Figure 10). We consider this to be an indication that at least one of the radioactive metabolite products of [11C]donepezil has marked uptake along with the main tracer. Indeed, Funaki et al.  reported that affinity of the M1 metabolite for AChE is almost as high as the affinity of donepezil for AChE; the M1 metabolite is radioactive when the parent tracer is [5-11C-methoxy]donepezil. In a steady-state study with a regular administration of donepezil, Meier-Davis et al.  found M1 to be relatively more prominent in minipigs than in humans and rats.
Meier-Davis et al.  also found the overall level of metabolites to be higher in the pigs. This species difference may explain why we found relatively fast metabolization of [11C]donepezil (Figure 4), in contrast to the human study by Hiraoka et al.  who saw only minor metabolism of [11C]donepezil (>85% remaining after 30 minutes) and therefore did not need metabolite correction. Also, Hiraoka et al. found the rev2TCM to be unequivocally better than the 1TCM, while our study finds the distinction less clear. As noted, however, their input functions were not corrected for metabolites.
Regarding [11C]donepezil as an infection tracer, our data did show an overall higher first-pass uptake in the infected lesions and a strong correlation with perfusion (Figure 10), while the DV was only slightly higher in the infected sites than in the corresponding noninfected locations (Figure 14), and the difference was not statistically significant.
The pig model was developed as a model for osteomyelitis, for which reason only relatively few soft tissue lesions were available, limiting the scope of the study as a general infection study. For [68Ga]Ga-citrate, only relatively limited data were available.
[18F]FDG was reasonably well described by the irr2TCM (irreversible uptake, three rate constants) for the 60-minute length studied, and, for both bone and soft tissue, [18F]FDG showed increased uptake in infected tissue (Figure 12). The correlation with perfusion indicated that the tracer is mainly flow-limited.
[68Ga]Ga-citrate was also well described by irrTCM but showed very little or very slow uptake, which was a limitation for infection imaging. The difference between infected and noninfected sites appeared to be higher in soft tissue rather than bone lesions, but too little data were available to draw a conclusion (Figure 12). Uptake was slow and diffusion-limited. To allow time for uptake, “late” imaging is preferable, but not so late that the radionuclide has decayed. Imaging at 120 minutes p.i. appears a good compromise.
[11C]methionine needed a rev2TCM (reversible uptake, four rate constants) for modelling. Despite considerable metabolism during the 40 minutes modelled, the , , and rate constants were only slightly affected if an uncorrected input function was used, that is, only was markedly affected. In a majority of cases, the uptake (measured as the distribution volume) was elevated in the infected tissue compared to the noninfected tissue (Figure 14), but the difference was less than for [18F]FDG (measured as the net uptake rate, Figure 12). Uptake appeared flow-limited. Based on the activity curves, imaging at ~15 minutes p.i. appears favourable.
[11C]donepezil could be modelled with a 1TCM (reversible uptake, two rate constants) but required metabolite correction—at least in this juvenile, porcine model. The uptake of labelled metabolite products appeared to be nonnegligible. Based on the present study, the uptake of [11C]donepezil in osteomyelitis seems to depend more on perfusion (flow-limited) than on differences between infected and noninfected tissues. For soft tissue infection, too few data were available to draw a conclusion.
Overall, among the studied PET tracers [18F]FDG showed optimal characteristics for the detection of infectious foci. [68Ga]Ga-citrate and [11C]donepezil were not found to be useful for imaging of osteomyelitis. For soft tissue, [11C]methionine and perhaps [68Ga]Ga-citrate may be applicable to quantify different aspects of inflammatory or infectious processes, while too few soft tissue data on [11C]donepezil were available to draw any conclusions.
The funder had no influence on the design, analysis, interpretation, or the writing of the manuscript.
Conflicts of Interest
The authors declare no conflicts of interest.
The authors are grateful for the support provided by the technical staff at Aarhus University Hospital, Aalborg University Hospital, and Copenhagen University. Vesa Oikonen is acknowledged for stimulating discussions and helping with the modelling software. This work was supported by the Danish Council for Independent Research, Technology and Production Sciences, Grant no. 0602-01911B (11-107077).
- S. Basu, T. Chryssikos, S. Moghadam-Kia, H. Zhuang, D. A. Torigian, and A. Alavi, “Positron emission tomography as a diagnostic tool in infection: present role and future possibilities,” Seminars in Nuclear Medicine, vol. 39, no. 1, pp. 36–51, 2009.
- S. Ahmad Sarji, “Physiological uptake in FDG PET simulating disease,” Biomedical Imaging and Intervention Journal, vol. 2, no. 4, 2006.
- OL. Nielsen, P. Afzelius, D. Bender et al., “Comparison of autologous 111In-leukocytes, 18F-FDG, 11C-methionine, 11C-PK11195 and 68Ga-citrate for diagnostic nuclear imaging in a juvenile porcine haematogenous staphylococcus aureus osteomyelitis model,” American Journal of Nuclear Medicine and Molecular Imaging, vol. 5, no. 2, pp. 169–182, 2015.
- P. Afzelius, O. L. Nielsen, A. K. O. Alstrup et al., “Biodistribution of the radionuclides 18F-FDG, 11C-methionine, 11C-PK11195, and 68Ga-citrate in domestic juvenile female pigs and morphological and molecular imaging of the tracers in hematogenously disseminated Staphylococcus aureus lesions,” American Journal of Nuclear Medicine and Molecular Imaging, vol. 6, no. 1, pp. 42–58, 2016.
- P. Afzelius, A. K. O. Alstrup, H. C. Schønheyder et al., “Utility of 11C-methionine and 11C-donepezil for imaging of Staphylococcus aureus induced osteomyelitis in a juvenile porcine model: comparison to autologous 111In-labelled leukocytes, , and 18F-FDG,” American Journal of Nuclear Medicine and Molecular Imaging, vol. 6, no. 6, pp. 286–300, 2016.
- L. Jødal, O. L. Nielsen, P. Afzelius, A. K. O. Alstrup, and S. B. Hansen, “Blood perfusion in osteomyelitis studied with [15O]water PET in a juvenile porcine model,” EJNMMI Research, vol. 7, no. 1, article no. 4, 2017.
- S. B. Jensen, K. M. Nielsen, D. Mewis, and J. Kaufmann, “Fast and simple one-step preparation of 68Ga citrate for routine clinical PET,” Nuclear Medicine Communications, vol. 34, no. 8, pp. 806–812, 2013.
- T. J. Mäkinen, P. Lankinen, T. Pöyhönen, J. Jalava, H. T. Aro, and A. Roivainen, “Comparison of 18F-FDG and 68Ga PET imaging in the assessment of experimental osteomyelitis due to Staphylococcus aureus,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 32, no. 11, pp. 1259–1268, 2005.
- C. Nanni, C. Errani, and L. Boriani, “68Ga-citrate PET/CT for evaluating patients with infections of the bone: preliminary results,” Journal of Nuclear Medicine, vol. 51, no. 12, pp. 1932–1936, 2010.
- M. Vorster, A. Maes, A. Jacobs et al., “Evaluating the possible role of 68Ga-citrate PET/CT in the characterization of indeterminate lung lesions,” Annals of Nuclear Medicine, vol. 28, no. 6, pp. 523–530, 2014, http://link.springer.com/10.1007/s12149-014-0842-9.
- S. Zhao, Y. Kuge, M. Kohanawa et al., “Usefulness of 11C-methionine for differentiating tumors from granulomas in experimental rat models: a comparison with 18F-FDG and 18F-FLT,” Journal of Nuclear Medicine, vol. 49, no. 1, pp. 135–141, 2008.
- K. Hirata, T. Shiga, N. Fujima et al., “11C-methionine positron emission tomography may monitor the activity of encephalitis,” Acta Radiologica, vol. 53, no. 10, pp. 1155–1157, 2012.
- Y. Maeda, H. Oguni, Y. Saitou et al., “Rasmussen syndrome: multifocal spread of inflammation suggested from MRI and PET findings,” Epilepsia, vol. 44, no. 8, pp. 1118–1121, 2003.
- M. Morooka, K. Kubota, H. Kadowaki et al., “11C-methionine PET of acute myocardial infarction,” Journal of Nuclear Medicine, vol. 50, no. 8, pp. 1283–1287, 2009.
- H. Sugimoto, H. Ogura, Y. Arai, Y. Iimura, and Y. Yamanishi, “Research and development of donepezil hydrochloride, a new type of acetylcholinesterase inhibitor,” Japanese Journal of Pharmacology, vol. 89, no. 1, pp. 7–20, 2002.
- K. Kawashima, T. Fujii, Y. Moriwaki, and H. Misawa, “Critical roles of acetylcholine and the muscarinic and nicotinic acetylcholine receptors in the regulation of immune function,” Life Sciences, vol. 91, no. 21-22, pp. 1027–1032, 2012.
- T. Fujii, Y. Watanabe, K. Fujimoto, and K. Kawashima, “Expression of acetylcholine in lymphocytes and modulation of an independent lymphocytic cholinergic activity by immunological stimulation,” Biogenic Amines, vol. 17, no. 4-6, pp. 373–386, 2003.
- N. P. Jørgensen, A. K. O. Alstrup, F. V. Mortensen et al., “Cholinergic PET imaging in infections and inflammation using 11C-donepezil and 18F-FEOBV,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 44, no. 3, pp. 449–458, 2017.
- L. K. Johansen, J. Koch, D. Frees et al., “Pathology and Biofilm Formation in a Porcine Model of Staphylococcal Osteomyelitis,” Journal of Comparative Pathology, vol. 147, no. 2-3, pp. 343–353, 2012.
- L. K. Johansen, E. L. Svalastoga, D. Frees et al., “A new technique for modeling of hematogenous osteomyelitis in pigs: Inoculation into femoral artery,” Journal of Investigative Surgery, vol. 26, no. 3, pp. 149–153, 2013.
- A. K. O. Alstrup, K. M. Nielsen, H. C. Schønheyder et al., “Refinement of a hematogenous localized osteomyelitis model in pigs,” Scandinavian Journal of Laboratory Animal Science, vol. 42, no. 3, pp. 1–4, 2016.
- L. Jødal, S. B. Hansen, and S. B. Jensen, “Impact of contamination with long-lived radionuclides on PET kinetics modelling in multitracer studies,” Nuclear Medicine Communications, vol. 37, no. 8, pp. 818–824, 2016.
- P. Som, H. L. Atkins, D. Bandoypadhyay et al., “A fluorinated glucose analog, 2-fluoro-2-deoxy-D-glucose (F-18): Nontoxic tracer for rapid tumor detection,” Journal of Nuclear Medicine, vol. 21, no. 7, pp. 670–675, 1980.
- M.-F. Tsan, “Mechanism of gallium-67 accumulation in inflammatory lesions,” Journal of Nuclear Medicine, vol. 26, no. 1, pp. 88–92, 1985.
- J. Van den Hoff, W. Burchert, W. Muller-Schauenburg, G.-J. Meyer, and H. Hundeshagen, “Accurate local blood flow measurements with dynamic PET: Fast determination of input function delay and dispersion by multilinear minimization,” Journal of Nuclear Medicine, vol. 34, no. 10, pp. 1770–1777, 1993.
- C. S. Patlak, R. G. Blasberg, and J. D. Fenstermacher, “Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data,” Journal of Cerebral Blood Flow and Metabolism, vol. 3, no. 1, pp. 1–7, 1983.
- C. S. Patlak and R. G. Blasberg, “Graphical Evaluation of Blood-to-Brain Transfer Constants from Multiple-Time Uptake Data. Generalizations,” Journal of Cerebral Blood Flow & Metabolism, vol. 5, no. 4, pp. 584–590, 2016.
- R. B. Innis, V. J. Cunningham, and J. Delforge, “Consensus nomenclature for in vivo imaging of reversibly binding radioligands,” Journal of Cerebral Blood Flow & Metabolism, vol. 27, no. 9, pp. 1533–1539, 2007.
- J. Logan, J. S. Fowler, N. D. Volkow et al., “Graphical analysis of reversible radioligand binding from time-activity measurements applied to [N-11C-methyl]-(-)-cocaine PET studies in human subjects,” Journal of Cerebral Blood Flow & Metabolism, vol. 10, no. 5, pp. 740–747, 1990.
- J. Logan, “Graphical analysis of PET data applied to reversible and irreversible tracers,” Nuclear Medicine and Biology, vol. 27, no. 7, pp. 661–670, 2000.
- F. Thiele and R. Buchert, “Evaluation of non-uniform weighting in non-linear regression for pharmacokinetic neuroreceptor modelling,” Nuclear Medicine Communications, vol. 29, no. 2, pp. 179–188, 2008.
- M. Yaqub, R. Boellaard, M. A. Kropholler, and A. A. Lammertsma, “Optimization algorithms and weighting factors for analysis of dynamic PET studies,” Physics in Medicine and Biology, vol. 51, no. 17, article no. 007, pp. 4217–4232, 2006.
- TPC. List of applications in TPCCLIB. Turku PET Centre web site. Available from: http://www.turkupetcentre.net/programs/doc/index.html.
- H. Akaike, “A new look at the statistical model identification,” IEEE Transactions on Automatic Control, vol. 19, pp. 716–723, 1974.
- K. P. Burnham and D. R. Anderson, “Multimodel inference: understanding AIC and BIC in model selection,” Sociological Methods and Research, vol. 33, no. 2, pp. 261–304, 2004.
- L. Sokoloff, M. Reivich, and C. Kennedy, “The [14C]deoxyglucose method for the measurement of local cerebral glucose utilization: theory, procedure, and normal values in the conscious and anesthetized albino rat,” Journal of Neurochemistry, vol. 28, no. 5, pp. 897–916, 1977.
- M. E. Phelps, S. C. Huang, E. J. Hoffman, C. Selin, L. Sokoloff, and D. E. Kuhl, “Tomographic measurement of local cerebral glucose metabolic rate in humans with (F-18)2-fluoro-2-deoxy-D-glucose: validation of method,” Annals of Neurology, vol. 6, no. 5, pp. 371–388, 1979.
- T. Schroeder, M. F. Vidal Melo, and J. G. Venegas, “Analysis of 2-[Fluorine-18]-Fluoro-2-deoxy-D-glucose uptake kinetics in PET studies of pulmonary inflammation,” Academic Radiology, vol. 18, no. 4, pp. 418–423, 2011.
- S. Chefer, D. Thomasson, J. Seidel et al., “Modeling [18F]-FDG lymphoid tissue kinetics to characterize nonhuman primate immune response to Middle East respiratory syndrome-coronavirus aerosol challenge,” EJNMMI Research, vol. 5, no. 1, article no. 65, pp. 1–11, 2015, http://www.ejnmmires.com/content/5/1/65.
- T. Schroeder, M. F. Vidal Melo, G. Musch, R. S. Harris, J. G. Venegas, and T. Winkler, “Modeling Pulmonary Kinetics of 2-Deoxy-2-[18F]fluoro-d-glucose During Acute Lung Injury,” Academic Radiology, vol. 15, no. 6, pp. 763–775, 2008.
- V. Kumar, D. K. Boddeti, S. G. Evans, F. Roesch, and R. Howman-Giles, “Potential use of 68Ga-apo-transferrin as a PET imaging agent for detecting Staphylococcus aureus infection,” Nuclear Medicine and Biology, vol. 38, no. 3, pp. 393–398, 2011.
- J. T. Thackeray, J. P. Bankstahl, Y. Wang et al., “Targeting post-infarct inflammation by PET imaging: comparison of 68Ga-citrate and 68Ga-DOTATATE with 18F-FDG in a mouse model,” European Journal of Nuclear Medicine and Molecular Imaging, vol. 42, no. 2, pp. 317–327, 2014.
- V. Kumar and D. K. Boddeti, “68Ga-radiopharmaceuticals for PET imaging of infection and inflammation,” in Theranostics, Gallium-68, and Other Radionuclides, R. P. Baum and F. Rösch, Eds., pp. 189–219, Springer, Berlin, Heidelberg, 2013, http://link.springer.com/10.1007/978-3-642-27994-2_11.
- E. K. J. Pauwels, V. R. McCready, J. H. M. B. Stoot, and D. F. P. Van Deurzen, “The mechanism of accumulation of tumour-localising radiopharmaceuticals,” European Journal of Nuclear Medicine, vol. 25, no. 3, pp. 277–305, 1998.
- A. J. Fischman, Y. M. Yu, E. Livni et al., “Muscle protein synthesis by positron-emission tomography with L-[methyl- 11C]methionine in adult humans,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 22, pp. 12793–12798, 1998.
- Y. Funaki, M. Kato, R. Iwata et al., “Evaluation of the binding characteristics of [5-11C-methoxyl]donepezil in the rat brain for in vivo visualization of acetylcholinesterase,” Journal Pharmacological Sciences, vol. 91, no. 2, pp. 105–112, 2003.
- S. R. Davis and R. Murgasova, “Comparison of metabolism of donepezil in rat, mini-pig and human, following oral and transdermal administration, and in an in vitro model of human epidermis,” Journal of Drug Metabolism & Toxicology, vol. 3, no. 4, 2012.
- K. Hiraoka, N. Okamura, Y. Funaki et al., “Quantitative analysis of donepezil binding to acetylcholinesterase using positron emission tomography and [5-11C-methoxy]donepezil,” NeuroImage, vol. 46, no. 3, pp. 616–623, 2009.
Copyright © 2017 Lars Jødal et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.