- Original research
- Open Access
- Published:

# Reproducible quantification of cardiac sympathetic innervation using graphical modeling of carbon-11-meta-hydroxyephedrine kinetics with dynamic PET-CT imaging

*EJNMMI Research*
**volume 8**, Article number: 63 (2018)

## Abstract

### Background

Graphical methods of radiotracer kinetic modeling in PET are ideal for parametric imaging and data quality assurance but can suffer from noise bias. This study compared the Logan and Multilinear Analysis-1 (MA1) graphical models to the standard one-tissue-compartment (1TC) model, including correction for partial-volume effects, in dynamic PET-CT studies of myocardial sympathetic innervation in the left ventricle (LV) using [^{11}C]HED.

### Methods

Test and retest [^{11}C]HED PET imaging (47 ± 22 days apart) was performed in 18 subjects with heart failure symptoms. Myocardial tissue volume of distribution (V_{T}) was estimated using Logan and MA1 graphical methods and compared to the 1TC standard model values using intraclass correlation (ICC) and Bland-Altman analysis of the non-parametric reproducibility coefficient (NPC).

### Results

A modeling start-time of *t** = 5 min gave the best fit for both Logan and MA1 (R^{2} = 0.95) methods. Logan slightly underestimated V_{T} relative to 1TC (*p* = 0.002), whereas MA1 did not (*p* = 0.96). Both the MA1 and Logan models exhibited good-to-excellent agreement with the 1TC (MA1-1TC ICC = 0.96; Logan-1TC ICC = 0.93) with no significant differences in NPC between the two comparisons (*p* = 0.92). All methods exhibited good-to-excellent test-retest repeatability with no significant differences in NPC (*p* = 0.57).

### Conclusions

Logan and MA1 models exhibited similar agreement and variability compared to the 1TC for modeling of [^{11}C]HED kinetics. Using *t** = 5 min and partial-volume correction produced accurate estimates of V_{T} as an index of myocardial sympathetic innervation.

## Background

Developed as a positron emission tomography (PET) imaging agent to target the cardiac sympathetic nervous system, carbon-11-labeled meta-hydroxyephedrine ([^{11}C]HED) is a norepinephrine analog that is taken up by nerve terminal varicosities in the myocardium, and used to assess sympathetic nerve function [1]. Since its genesis, it has been the cornerstone PET tracer for cardiac sympathetic innervation, employed in determination of neuronal-based defects leading to improved diagnosis and prognosis for pathologies such as heart failure, arrhythmia, and cardiomyopathy, in which cardiac neuronal function is often compromised, leading to decreased catecholamine sensitivity and lowered beta adrenergic receptor density [1]. Using PET [^{11}C]HED imaging of cardiac tissues, the volume of distribution (V_{T}) of the injected radiotracer is an invaluable metric that quantifies the uptake and retention of tracer, providing an index of sympathetic nerve density and reuptake-1 transporter activity. For cardiac PET applications especially, V_{T} and other kinetic modeling parameters measured in the myocardium may be used to aid in the diagnosis of various innervation and perfusion-based pathologies.

In PET imaging studies, V_{T} is defined as the equilibrium ratio of tracer concentration in tissue to that of unmetabolized parent tracer in plasma, but this direct measurement is typically not feasible due to the long time needed to reach equilibrium. Alternatively, kinetic modeling is commonly used to determine V_{T} from a significantly shorter temporal sample following tracer injection [2]. While the physiological kinetics of [^{11}C]HED may be modeled using a two-tissue-compartment model, the one-tissue-compartment (1TC) model has been shown to provide a robust representation with optimal clinical reproducibility in myocardial uptake studies, without sacrificing the accuracy of V_{T} quantification [3]. Two graphical methods reported in the literature for kinetic modeling of reversible-binding tracers are the Logan [4] and Multilinear Analysis-1 (MA1) models [5], which are both computationally simpler than non-graphical (compartmental) methods [6], while being able to provide visual representations of kinetic parameters. The Logan method has been established as the standard graphical model to estimate V_{T} in a wide range of PET applications in the brain and heart, while MA1 was proposed as an alternative numerical formulation to estimate V_{T} with lower noise bias compared to Logan estimates [5]. Although [^{11}C]HED is a widely used tracer, a comprehensive evaluation of the performance of the graphical and non-graphical methods to quantify its kinetics has not been performed. Furthermore, the effects of partial-volume losses on quantification of V_{T} have not been well defined in the context of graphical kinetic modeling in the heart, where the effects of blood-pool spillover and motion are more apparent compared to the brain. The goal of this study was to determine a method of partial-volume correction applicable to graphical kinetic modeling and to compare the Logan and MA1 models to the standard 1TC kinetic model for accurate quantification of myocardial sympathetic innervation using dynamic [^{11}C]HED PET-CT studies.

## Methods

### Patient study design

Twenty-three heart failure patients were recruited as control subjects for a previous study (PET-OSA: NCT00756366) investigating the effects of continuous positive airway pressure (CPAP) on sympathetic nerve function and cardiac energetics in heart failure patients with obstructive sleep apnea (OSA) [7]. These control patients had the same inclusion and exclusion criteria as the PET-OSA study, except they did not have OSA. Patient demographics were collected at baseline and follow-up visits.

Three of the 23 patients were missing baseline or follow-up PET scans and were excluded. Two other patients were also excluded: one with atrial fibrillation at baseline that was treated before the follow-up scan and the other with uncorrectable severe motion artifact, leaving *N* = 18 subjects included in the final analysis. All patients provided written informed consent, according to the research protocol approved by the Human Research Ethics Board at the University of Ottawa Heart Institute.

### [^{11}C]HED PET imaging

[^{11}C]HED was synthesized from [^{11}C]methyl-iodide and metaraminol-free base, with the use of standard methods for high purity and specific activity [8]. Images were obtained at baseline and follow-up (47 ± 22 days apart) on the ECAT-ART PET (Siemens/CTI, Knoxville, TN) or Discovery RX PET-VCT (GE Healthcare, Waukesha, WI) scanner, with ECG, heart rate, and blood pressure monitored at regular intervals. A transmission scan for attenuation correction was performed using Cs-137 isotope or X-ray CT [9], immediately after which 10–15 mCi (370–550 MBq) of [^{11}C]HED was injected over 30 s and a dynamic PET series was acquired over a 40-min period (10 × 10 s; 1 × 60 s; 5 × 100 s; 3 × 180 s; 4 × 300 s) [10]. Image reconstruction was performed using filtered-back-projection with a 12-mm Hann filter and all corrections enabled for quantification of radioactivity concentration [11].

### Tracer kinetic modeling

#### Blood metabolites correction

Quantitative analysis of [^{11}C]HED kinetics requires correction for radiolabeled metabolites that accumulate over time in the bloodstream, which are not present in the myocardium [12]. The arterial whole-blood tracer concentration C_{WB}(*t*) is typically measured using an image-derived region of interest (ROI) placed in the LV cavity and must be differentiated from the unchanged parent tracer concentration in plasma C_{p}(t) as defined using the standardized nomenclature of Innis et al. [13]. The relation between C_{p}(*t*) and C_{WB}(*t*) was characterized as a time-varying function of plasma-to-whole-blood and unchanged parent-to-metabolized radiotracer in the bloodstream and expressed as a combined parent fraction in plasma pfp(t) function (Additional file 1: Figure S1) derived from studies performed previously in humans [3]:

#### Compartment modeling and partial volume correction

The tracer volume of distribution V_{T} in the myocardium is defined as the ratio of the concentration of tracer in tissue divided by the concentration of tracer in arterial plasma, after the system has reached equilibrium (at *t* ≥ T_{E}) [13].

Where C_{T} is the concentration of tracer in myocardial tissue and C_{p} is the concentration of tracer in plasma. For a reversible 1TC kinetic model, the rate-of-change of tracer concentration in myocardial tissue is defined using the rate of influx from arterial plasma-to-tissue (*K*_{1}) and the rate of efflux from the tissue compartment (*k*_{2}) according to Eq. 3.

At equilibrium (*t* ≥ T_{E}), the rate-of-change of tracer concentration in tissue is equal to zero (\( \frac{d{C}_T(t)}{dt}=0\Big) \) [14]. Combining (2) and (3), the volume of distribution may be expressed as:

This simple derivation is applied widely in the analysis of neuro-PET imaging studies. However, in cardiac PET applications, additional image blurring due to cardiac contractile and respiratory motion makes it difficult to deduce the exact boundaries of myocardial tissue based on the measured ROI. There is also a 10–15% fraction of blood volume within normal myocardial tissue that must be considered. These effects may be lessened with modern PET-CT hybrid scanners with improved spatial and contrast resolution, but for cardiac imaging, these effects are still pronounced, necessitating implementation of partial-volume corrections [15, 16]. With partial-volume spillover considered, our model of the imaging process becomes:

where C_{ROI}(*t*) is the measured tracer concentration in the PET myocardial image ROI. F_{WB} is the fraction of whole-blood signal C_{WB}(*t*) contained in the measured ROI curve due to imaging spillover effects and anatomical blood volume in the myocardial tissue [11]. C_{T}(*t*) is the tracer concentration in the myocardial tissue (excluding blood), and RC is the partial-volume recovery coefficient describing the fractional underestimation of C_{T}(*t*) due to limited spatial resolution and myocardial motion blurring. In this study, the value of RC was estimated regionally as 1 − F_{WB,} according to the method of Hutchins et al. [15] used commonly in the compartmental analysis of cardiac PET dynamic imaging studies. Blood spillover from the right ventricle cavity to the interventricular septum was not modeled explicitly. Spillover from the myocardium to blood-pool was not corrected, which might affect C_{WB}(*t*) at later time points. Isolating for C_{T}(*t*) we have:

Substituting C_{T} from (5) into the definition of \( {V}_T=\frac{C_T\left({T}_E\right)}{C_p\left({T}_E\right)} \) in (1) and assuming *t* ≥ *T*_{E}, we obtain:

C_{WB} may be expressed in terms of C_{p} from (1), only considering *t* ≥ T_{E}:

Maintaining the same logic as (2) and distinguishing between V_{T} for the volume of distribution that corresponds to the true myocardial tissue compartment (C_{T}) and V_{ROI} for the volume of distribution that corresponds to the measured region of interest (C_{ROI}), V_{ROI} may be expressed as the ratio: \( \frac{C_{ROI}\left({T}_E\right)}{C_p\left({T}_E\right)}={V}_{ROI} \). This can hence be substituted into (9); then, C_{p}(t) may also be canceled from the second term of (9), yielding:

where pfp(*T*_{E}) represents the equilibrium value of pfp(*t*). Then, (10) can be rearranged to isolate V_{T} as:

From (11), we propose that V_{T} may be estimated from V_{ROI}, with plasma-to-whole-blood and metabolite corrections as well as partial-volume effects considered explicitly.

### Graphical kinetic modeling

The Logan model (12) was derived from the first-order differential equations for general compartment models. Its purpose was to create a graphical method of quantifying V_{T}, while making kinetic modeling more mathematically and computationally simple and robust. The base equation, adapted from the original formulation [5] to fit the reversible 1TC model in the context of cardiac PET, is:

C_{ROI}(*t*) and C_{p}(*t*) time-activity curves are used as measured input data. At a certain time (*t**), the intercept term (Int) will become a constant value [17], at which point the equation becomes a linear system where the slope represents the volume of distribution in the ROI. Since the measured tissue curve C_{ROI}(*t*) is subject to blood spillover and partial-volume losses, only V_{ROI} may be obtained from the graphical model directly. Previous applications of this model have estimated V_{ROI} without explicit correction factors for partial-volume effects, which is required for cardiac applications. In our proposed model, Eq. (11) may be used to determine V_{T} from the slope determined by the Logan model.

The second graphical method investigated is the MA1 model, originally formulated as a more numerically stable alternative to the Logan model [6]:

As with the Logan model, only *T* > *t** are used for MA1 analysis. MA1 is a multilinear equation with two independent variables, and the corresponding Logan slope V_{ROI} is equal to the negative ratio of the two coefficients, such that:

V_{T}may again be determined from the V_{ROI} value estimated using MA1, according to the relation defined in (11).

#### Determination of *t** for graphical models

The estimation start-time (*t**) was varied systematically from 1.5 to 20 min for a subset of five [^{11}C]HED studies to determine the optimal value to be used for the main analysis. Goodness-of-fit was evaluated on the Logan plot as the Pearson correlation (*r*^{2}) of the points from *t** to 40 min, indicating the subset of points best described by a line. Since the *r*^{2} is not effective to assess goodness-of-fit of the near-horizontal fitted plane on the MA1 plots, an alternative metric was computed using the relative standard error of the estimate (rSEE) as 1 − SEE/mean. The optimal *t** was determined by comparing V_{T} values from the graphical methods to the 1TC model standard. Then, all subsequent analysis was performed using the same start-time for both Logan and MA1 models.

#### PET image analysis

The compartmental and graphical analysis models were implemented in the FlowQuant® analysis software (University of Ottawa Heart Institute, ON). The operator reliability of this automated software has been reported previously [18]. Briefly, the left ventricle (LV) myocardium was segmented automatically and partitioned into voxels using a 2D polar-map representation, with each voxel representing a transmural sub-region of the LV myocardial tissue. The arterial whole-blood (WB) ROI was positioned automatically at the center of the left atrioventricular valve plane. Time-activity curves were generated based on measured tracer activity in the LV cavity C_{WB}(*t*) and myocardial tissue C_{T}(t) ROIs, as input to the tracer kinetic models.

In each polar-map voxel, the 1TC model rate constants K_{1} and k_{2}, as well as V_{T} and the blood spillover fraction F_{WB}, were estimated using weighted least-squares regression, according to Eqs. (3), (4), and (5). The Logan and MA1 graphical models in Eqs. (12) and (13) were used to calculate LV polar-maps of V_{ROI}. Scan-specific spillover values were calculated as the polar-map median F_{WB} and the corresponding partial-volume recovery coefficient RC (1 − F_{WB}), which were then used to estimate V_{T} from the graphical model estimates of V_{ROI} according to Eq. (11). Image and data analyses were performed using MATLAB 2013b (The Mathworks, Natick, MA).

### Statistical analysis

LV median V_{T} values obtained from the 1TC, Logan, and MA1 methods were tabulated. Inter-model and test-retest mean effects were evaluated using two-way repeated measures ANOVA. Bland-Altman analyses and Intra-class correlation (ICC) were employed to evaluate the inter-model (MA1 vs 1TC, and Logan vs 1TC) and test-retest (baseline vs follow-up) reliability [19, 20]. Absolute-agreement ICC with two-way mixed effects was used for the inter-model reproducibility and test-retest repeatability [21]. To correct for skew in the V_{T} distributions, V_{T} values were logarithmically transformed before the ANOVA and ICC analyses. ICC values were categorized as: ICC > 0.90 excellent, > 0.75 very good, > 0.40 good, and ≤ 0.40 poor [22]. The limits-of-agreement of repeated measures were estimated using the following: (i) median difference ± non-parametric repeatability coefficient (NPC = 1.45 × IQR) to account for the variable effect of outliers and (ii) mean difference ± coefficient-of-repeatability (CR = 1.96 × SD). √(3/N) × SD_{difference}(*t*_{95%, n − 1}) was used to calculate the 95% confidence intervals on the limits-of-agreement, where *N* is the number of pairs being analyzed [20]. Differences in V_{T} values were divided by the mean V_{T} to account for the increased variability of differences associated with increased mean V_{T}. The NPC was also reported as it is a more robust measure of repeatability [23]. Non-parametric Levene’s test was used to assess the equality of variance between groups. Bias in the Bland-Altman plots was assessed using the one-sample Wilcoxon Signed Ranked test against zero. A 2-tailed *p* value < 0.05 was considered statistically significant for all tests. Statistical analyses were performed using Excel 2016 (Microsoft) and SPSS 20.0 (IBM).

## Results

### Patient demographics

Baseline patient demographics are listed in Table 1. Most patients in this study were male (66.5%) with mean age of 66.5 ± 9.3 years. The majority (83.3%) were classified as having NYHA Class II heart failure and were taking one or more cardiac medications. The demographics were stable at follow-up compared those reported at baseline.

### Adjustment of start-time (*t**) for graphical models

Table 2 shows the goodness-of-fit metrics for the five patients randomly sampled from the entire cohort used for this study. These ranged from 0.80 to 0.99 across all scans and *t** values of 1.5–20 min, as summarized in Fig. 1. Corresponding V_{T} values ranged from 9 to 21 mL/cm^{3} for Logan and 10–28 mL/cm^{3} for MA1. The Logan V_{T} increased systematically up to approximately 5–15 min, interpreted as the start of the steady-state (linear) phase. There was no *t** value with Logan V_{T} estimates equal to the 1TC reference value (20 mL/cm^{3}); therefore, *t** = 5 min was selected as the optimal start-time based on agreement of the MA1 V_{T} values with the reference 1TC model. This value of *t** also demonstrated the highest Logan *r*^{2} value (0.96), suggesting the best fit of a line was obtained for the points starting at 5 min. MA1 plots exhibited a steady increase in goodness-of-fit up to 5 min (0.97) with relatively little improvement at later start times.

### PET image analysis

Figure 2 shows V_{T} polar-maps from a single patient scan using all 3 models, including graphical representations of the Logan and MA1 plots as well as the 1TC modeling results in Fig. 3. V_{T} polar-maps were found to show very similar spatial distributions for all three kinetic models, as expected.

### Comparison of Logan and MA1 versus 1TC

At baseline, V_{T} values were 20 ± 8 mL/cm^{3} for 1TC, 17 ± 8.0 mL/cm^{3} for Logan, and 20 ± 16 mL/cm^{3} for MA1, as shown in Table 3. At follow-up, V_{T} values were 21 ± 11 mL/cm^{3} for 1TC, 19 ± 12 mL/cm^{3} for Logan, and 23 ± 16 mL/cm^{3} for MA1. Intra-model comparison of the V_{T} values at baseline vs. follow-up revealed that there was no significant difference between baseline and follow-up V_{T} values for any of the models (*p* = 0.379). However, on average, V_{T} values generated by Logan were 15% lower than those generated by MA1 (*p* = 0.002) and 12% lower than the 1TC values (*p* = 0.002). V_{T} values generated by MA1 were not significantly different from those of 1TC (*p* = 0.958).

To evaluate reproducibility between the three models, the 1TC model was taken as the reference standard for comparison with the Logan and MA1 models. Of the 36 scans compared in the inter-model analysis (Table 4, Fig. 4), the Logan-vs-1TC and MA1-vs-1TC comparisons exhibited similar reproducibility with NPC ~ 26.5%. However, the V_{T} values generated from the Logan model were systematically lower than those generated from the 1TC model (median bias = − 14.5% and mean bias = − 16.3%, *p* < 0.001), but there was no systematic difference in V_{T} when comparing MA1-vs-1TC models, *p* = 0.2). There was excellent agreement between MA1-vs-1TC values (ICC = 0.955, 95% CI [0.915, 0.977]) and good-to-excellent agreement between Logan-vs-1TC (ICC = 0.928, 95% CI [0.432, 0.978]). There was no difference in reproducibility between the MA1-vs-1TC and Logan-vs-1TC NPC values (nonparametric Levene’s test, *p* = 0.915).

### Test-retest repeatability of kinetic models

All models demonstrated very good repeatability (Table 5) with consistent ICC values = 0.837–0.852. The mean test-retest differences were all < 2% without any systematic bias observed between baseline and follow-up (Fig. 5), but this could be the result of relatively small sample size with fewer points (*N* = 18) compared to the inter-method analysis (*N* = 36). There was no difference in the test-retest reproducibility (NPC) values between the three methods (non-parametric Levene’s test, *p* = 0.57).

## Discussion

In an effort to improve and expand the use of kinetic modeling in cardiac PET studies of sympathetic innervation, we sought to evaluate multiple kinetic models for the analysis of [^{11}C]HED studies. This was achieved by comparing the inter-method differences in V_{T} quantified by the Logan and MA1 graphical models compared to the reference 1TC model in a sample of heart failure patients and assessing the test-retest repeatability between baseline and follow-up scans. HED PET is often used to evaluate therapy or disease progression in heart failure patients; therefore, evaluation of the test-retest repeatability is most relevant in this same population, as opposed to healthy normal subjects who generally have lower sympathetic tone. The patients’ heart failure symptoms and medications were stable over the test-retest interval; therefore, any impact on the repeatability data should be minimal.

The MA1 model exhibited excellent agreement with 1TC, the Logan model exhibited good-to-excellent agreement with 1TC, and all models had good-to-excellent test-retest repeatability. Logan V_{T} values were significantly lower than MA1 and 1TC V_{T} values, while MA1 V_{T} values were not significantly different from those obtained using the 1TC model (Table 3). While 1TC is the reference standard kinetic model in this instance, graphical models such as the Logan and MA1 are computationally simpler alternatives that allow for linearized visualization and analysis of tracer kinetic data. Our findings support the reliable use of both graphical analysis methods in addition to the standard 1TC model for tracer kinetic analysis of V_{T}. These findings agree with previous studies using other PET tracers that compared various graphical models, including the Logan method, finding the results to be in agreement with standard compartment models, but computationally simpler, and potentially more robust [24,25,26,27].

In the present cardiac PET study, partial volume and spillover corrections were critical to implement into the graphical modeling calculations to avoid misinterpretation. The commonly used Logan and MA1 methods (Eqs. 12 and 13) only estimate the volume of distribution in the PET image region (V_{ROI}) as opposed to the myocardial tissue of interest (V_{T}). Compared to PET measurements in other organ systems such as the brain, in cardiac studies, the measured ROI region contains much more spillover of blood signal within and adjacent to the myocardial tissues. Our implementation of a partial-volume correction method based on estimated recovery coefficients and whole-blood spillover fractions allowed accurate measurement of myocardial V_{T} values using Logan and MA1 graphical models on a scan-specific basis. In this validation study, F_{WB} was estimated first using the 1TC with spillover model, and then used to calculate the corresponding RC values for consistent partial-volume and spillover correction of the graphical model V_{ROI} estimates. It is clear that independent estimates of RC and F_{WB} are required to determine V_{T} from V_{ROI} as shown in Eq. 11; therefore, any error in the estimation of these correction factors in practice will be propagated directly into the corresponding values of V_{T}. In the present study, the average F_{WB} value was 0.37 ± 0.07, which could be used to estimate RC and hence V_{T} in similar patient population studies with minimal added variability.

We investigated the effect of varying *t** on the graphical model results (Table 2), which quantified V_{T} using the plotted values at *t* ≥ *t**. It has been reported that *t** may be deduced directly from kinetic modeling data for some tracers [5], but the method we presented used a simpler and systematic approach to determine the *t** which produced the same V_{T} values on average compared to the MA1 plots. This approach is beneficial for tracers for which it is more difficult to estimate *t** directly from the study data, such as those with relatively slower kinetics [28]. It also removes the need to estimate *t** for each individual scan, which may be subject to variable noise effects. We propose *t** = 5 min as an effective start-time for cardiac studies employing [^{11}C]HED as it also gave the highest quality of linear fit (*r*^{2} > 0.95) using the Logan model, in addition to MA1 estimates of V_{T} that were equal to the 1TC reference value on average. This start time was shown with our comparison of the three models to be robust, producing results for V_{T} with excellent goodness-of-fit to the graphical models and inter-method agreement. It is worth noting that a slightly later start time of 10–15 min may have provided Logan V_{T} values that correspond better with 1TC and MA1 (Fig. 1), but at the cost of a lower quality fit of the linear model and wider variability due to fewer fitted points.

Interestingly, the V_{T} values determined by Logan were significantly lower than those determined by both MA1 and 1TC, while V_{T} values determined by MA1 did not show a significant difference to those obtained from 1TC. More precisely, Logan exhibited a greater negative bias where V_{T} was underestimated relative to 1TC, whereas a bias was not present between MA1 and 1TC (Table 4). In a similar kinetic model comparison using [^{18}F]FCWAY and [^{11}C]MDL neurological tracers, Ichise et al. [6] demonstrated that the MA1 model generated higher V_{T} estimates than Logan, and that MA1 exhibited less bias compared to Logan at multiple imaging noise levels. Our results are consistent with these findings, affirming the original report of MA1 as a method to reduce the magnitude of bias induced by noise when using the Logan model [6]. Although Logan seemed to underestimate V_{T} in our study population, it should be realized that the median bias of − 14.5% relative to the 1TC gold standard did not greatly affect the inter-model reproducibility of the models, which exhibited good to excellent agreement despite the bias that was present.

The use of [^{11}C]HED to examine sympathetic function in cardiac PET is becoming increasingly widespread. Recently, it has been shown to be a powerful diagnostic and prognostic tool for patients with heart failure, arrhythmias, flow-innervation mismatches, and microvascular dysfunction in both infarcted and non-infarcted tissues [1, 29,30,31,32,33]. This field continues to be improved and shows promise for a wider variety of applications [34]. As cardiac innervation tracers increase in prevalence, the optimization and validation of kinetic modeling techniques becomes more important; extensions of the current study may be anticipated, such as those investigating the use of a two-tissue-compartment model to quantify cardiac NET re-uptake function more specifically. Moreover, comparisons of multiple kinetic modeling options, in particular those of a graphical nature as presented here, are possible with other cardiac innervation-based tracers such as the [^{18}F]-labeled sympathetic innervation tracers MFBG, MHPG, LMI1195, etc., for more detailed evaluation of their kinetics [35, 36].

A few limitations were present in this study. The current study is a retrospective, single-center study that examined stable heart failure patients only from the PET-OSA trial. The results may be limited by the relatively small sample size (*N* = 18). Larger prospective studies would be beneficial to further validate the performance of the kinetic models as proposed.

## Conclusion

A start time of 5 min was found to provide the best fit for Logan and MA1 models. The MA1-1TC comparison demonstrated excellent agreement while Logan-1TC and test-retest comparisons demonstrated good-to-excellent agreement when quantifying V_{T} with partial volume correction. Although Logan underestimated V_{T} due to the recognized noise bias, Logan and MA1 both exhibited similar test-retest variability, suggesting that they may be used in addition to 1TC in the modeling of [^{11}C]HED kinetics, with benefits of greater computational simplicity and the ability to mathematically visualize kinetic parameters for better quality assurance.

## References

- 1.
Boschi S, Lodi F, Boschi L, Nanni C, Chondrogiannis S, Colletti PM, et al. 11C-Meta-Hydroxyephedrine. Clin Nucl Med. 2015;40(2):e96–103. Available from: https://www.ncbi.nlm.nih.gov/pubmed/24999701.

- 2.
Logan J, Fowler JS, Volkow ND, Wang GJ, Ding YS, Alexoff DL. Distribution volume ratios without blood sampling from graphical analysis of PET data. J Cereb Blood Flow Metab. 1996;16(5):834–40.

- 3.
Harms HJ, Haan SD, Knaapen P, Allart CP, Rijnierse MT, Schuit RC, et al. Quantification of [11C]-meta-hydroxyephedrine uptake in human myocardium. EJNMMI Res. 2014;4:52.

- 4.
Logan J, Fowler JS, Volkow ND, Wolf AP, Dewey SL, Schlyer DJ, et al. Graphical analysis of reversible radioligand binding from time-activity measurements applied to [N-11C-methyl]-(−)-cocaine PET studies in human subjects. J Cereb Blood Flow Metab. 1990;10(5):740–7. Available from: http://dx.doi.org/10.1038/jcbfm.1990.127

- 5.
Ichise M, Toyama H, Innis RB, Carson RE. Strategies to improve neuroreceptor parameter estimation by linear regression analysis. J Cereb Blood Flow Metab. 2002;22(10):1271–81. Available from: http://www.ncbi.nlm.nih.gov/pubmed/12368666

- 6.
Schmidt KC, Turkheimer FE. Kinetic modeling in positron emission tomography. Q J Nucl Med. 2002;46(1):70–85.

- 7.
Hall AB, Ziadi MC, Leech JA, Chen S-Y, Burwash IG, Renaud J, et al. Effects of short-term continuous positive airway pressure on myocardial sympathetic nerve function and energetics in patients with heart failure and obstructive sleep apnea: a randomized study. Circulation. 2014;130(11):892–901. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24993098

- 8.
Rosenspire KC, Haka MS, Van Dort ME, Jewett DM, Gildersleeve DL, Schwaiger M, et al. Synthesis and preliminary evaluation of carbon-11-meta-hydroxyephedrine: a false transmitter agent for heart neuronal imaging. J Nucl Med. 1990;31(8):1328–34. Available from: http://www.ncbi.nlm.nih.gov/pubmed/2384800

- 9.
Yoshinaga K, Burwash IG, Leech JA, Haddad H, Johnson CB, deKemp RA, et al. The effects of continuous positive airway pressure on myocardial energetics in patients with heart failure and obstructive sleep apnea. J Am Coll Cardiol. 2007;49(4):450–8.

- 10.
Allman KC, Wieland DM, Muzik O, Degrado TR, Wolfe ER, Schwaiger M. Carbon-11 hydroxyephedrine with positron emission tomography for serial assessment of cardiac adrenergic neuronal function after acute myocardial infarction in humans. J Am Coll Cardiol. 1993;22(2):368–75.

- 11.
Cherry SR, Sorenson JA, Phelps ME. Physics in Nuclear Medicine; Chapter 18: Positron Emission Tomography. Philadelphia: Elsevier Saunders; 2012.

- 12.
Thackeray JT, Beanlands RS, Dasilva JN. Presence of specific 11C-meta-Hydroxy-ephedrine retention in heart, lung, pancreas, and brown adipose tissue. J Nucl Med. 2007;48(10):1733–40. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17873125

- 13.
Innis RB, Cunningham VJ, Delforge J, Fujita M, Gjedde A, Gunn RN, et al. Consensus nomenclature for in vivo imaging of reversibly binding radioligands. J Cereb Blood Flow Metab. 2007;27:1533–9.

- 14.
Wahl R, Buchanan J. Principles and Practice of positron emission tomography. Philidelphia: Lippincott Williams & Wilkins; 2002.

- 15.
Hutchins GD, Caraher JM, Raylman RR. A region of interest strategy for minimizing resolution distortions in quantitative myocardial PET studies. J NucI Med. 1992;33:1243–50.

- 16.
Iida H, Rhodes CG, de Silva R, Yamamoto Y, Araujo LI, Maseri A, et al. Myocardial tissue fraction--correction for partial volume effects and measure of tissue viability. J Nucl Med. 1991;32(11):2169–75. Available from: http://www.ncbi.nlm.nih.gov/pubmed/1941156

- 17.
Hany TF, Steinert HC, Goerres GW, Buck A, von Schulthess GK. PET diagnostic accuracy: improvement with in-line PET-CT system: initial results. Radiology. 2002;225(2):575–81.

- 18.
Patlak CS, Blasberg RG. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. Generalizations. J Cereb Blood Flow Metab. 1985;5(4):584–90. Available from: http://journals.sagepub.com/doi/pdf/10.1038/jcbfm.1985.87.

- 19.
Klein R, Renaud JM, Ziadi MC, Thorn SL, Adler A, Beanlands RS, deKemp RA. Intra- and inter-operator repeatability of myocardial blood flow and myocardial flow reserve measurements using rubidium-82 pet and a highly automated analysis program. J Nucl Cardiol. 2010;17(4):600–16. https://doi.org/10.1007/s12350-010-9225-3.

- 20.
Shrout PE, Fleiss JL. Intraclass correlations: uses in assessing rater reliability. Psychol Bull. 1979;86(2):420–8.

- 21.
Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;327:307–10. Available from: http://www.sciencedirect.com/science/article/pii/S0140673686908378

- 22.
Koo TK, Li MY. A guideline of selecting and reporting Intraclass correlation coefficients for reliability research. J Chiropr Med. 2016;15(2):155–63. Available from: http://linkinghub.elsevier.com/retrieve/pii/S1556370716000158

- 23.
Fleiss JL. Reliability of measurement. In: The design and analysis of clinical experiments [Internet]; 1999. p. 1–32. Available from: http://dx.doi.org/10.1002/9781118032923.ch1.

- 24.
Klein R, Ocneanu A, Renaud JM, Ziadi MC, Beanlands RSB, DeKemp RA. Consistent tracer administration profile improves test-retest repeatability of myocardial blood flow quantification with (82)Rb dynamic PET imaging. J Nucl Cardiol. 2016; Available from: http://www.ncbi.nlm.nih.gov/pubmed/27804067

- 25.
Rizzo G, Veronese M, Zanotti-Fregonara P, Bertoldo A. Voxelwise quantification of [11C](R)-rolipram PET data: a comparison between model-based and data-driven methods. J Cereb Blood Flow Metab. 2013;33(7):1032–40. Available from: http://journals.sagepub.com/doi/abs/10.1038/jcbfm.2013.43.

- 26.
Treyer V, Streffer J, Wyss MT, Bettio A, Ametamey SM, Fischer U, et al. Evaluation of the metabotropic glutamate receptor subtype 5 using PET and 11C-ABP688: assessment of methods. J Nucl Med. 2007;48(7):1207–15.

- 27.
Grafström J, Stone-Elander S. Comparison of methods for evaluating radiolabelled Annexin A5 uptake in pre-clinical PET oncological studies. Nucl Med Biol. 2014;41(10):793–800.

- 28.
Wong K, Kepe V, Small G, Satyamurthy N, Barrio J, Huang S. Comparison of simplified methods for quantitative analysis of [F-18]FDDNP PET data. IEEE Nucl Sci Symp Conf Rec. 2007;1–11:3146.

- 29.
Wilks MQ, Knowles SM, Wu AM, Huang S-C. Improved modeling of in vivo kinetics of slowly diffusing radiotracers for tumor imaging. J Nucl Med. 2014;55(9):1539–44. Available from: http://jnm.snmjournals.org/content/early/2014/07/02/jnumed.114.140038

- 30.
Boutagy N, Sinusas A. Recent Advances and Clinical Applications of PET Cardiac Autonomic Nervous System ImagingNo Title. Curr Cardiol Rep. 2017;19(4):33. https://doi.org/10.1007/s11886-017-0843-0.

- 31.
Magota K, Hattori N, Manabe O, Naya M, Oyama-Manabe N, Shiga T, et al. Electrocardiographically gated 11C-hydroxyephedrine PET for the simultaneous assessment of cardiac sympathetic and contractile functions. Ann Nucl Med. 2014;28(3):187–95.

- 32.
Capitanio S, Nanni C, Marini C, Bonfiglioli R, Martignani C, Dib B, et al. Heterogeneous response of cardiac sympathetic function to cardiac resynchronization therapy in heart failure documented by 11[C]-hydroxy-ephedrine and PET/CT. Nucl Med Biol. 2015;42(11):858–63.

- 33.
Harms HJ, Lubberink M, de Haan S, Knaapen P, Huisman MC, Schuit RC, et al. Use of a Single 11C-Meta-Hydroxyephedrine Scan for Assessing Flow-Innervation Mismatches in Patients with Ischemic Cardiomyopathy. J Nucl Med. 2015;56(11):1706–11. Available from: http://jnm.snmjournals.org/cgi/doi/10.2967/jnumed.115.154377

- 34.
Rijnierse MT, Allaart CP, De Haan S, Harms HJ, Huisman MC, Wu L, et al. Sympathetic denervation is associated with microvascular dysfunction in non-infarcted myocardium in patients with cardiomyopathy. Eur Heart J Cardiovasc Imaging. 2015;16(7):788–98.

- 35.
Thackeray JT, Bengel FM. PET imaging of the autonomic nervous system. Q J Nucl Med Mol Imaging. 2016;60(4):362–82. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27611712

- 36.
Werner RA, Rischpler C, Onthank D, Lapa C, Robinson S, Samnick S, et al. Retention kinetics of the 18F-labeled sympathetic nerve PET tracer LMI1195: comparison with 11C-hydroxyephedrine and 123I-MIBG. J Nucl Med. 2015;56(9):1429–33. Available from: http://www.ncbi.nlm.nih.gov/pubmed/26182969

### Funding

Networks of Centres of Excellence of Canada (NCE-15-P06-001), Ontario Research Foundation (ORF-RE07-021).

### Availability of data and materials

The data will not be shared because it will be used in other upcoming studies.

## Author information

### Affiliations

### Contributions

TW performed kinetic analysis, created figures, and wrote the manuscript. KYW performed statistical analysis, created figures and tables, and wrote the manuscript with TW. TW and RdK formulated partial volume correction. RCM processed clinical studies and assisted in creation of figures. JMR assisted in implementation of kinetic analysis tools. RSB and RdK supervised project development and analysis. All authors were involved in the editing process. All authors read and approved the final manuscript.

### Corresponding author

## Ethics declarations

### Ethics approval and consent to participate

All research subjects provided written informed consent, as approved by the Human Research Ethics Board at the University of Ottawa Heart Institute.

### Consent for publication

Consent has been obtained from participants to publish this work.

### Competing interests

RSB and RdK have received unrestricted university-industry grant funding from the Ontario Research Fund and Lantheus Medical Imaging.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Additional file

### Additional file 1:

**Figure S1.** Unchanged parent fraction in plasma (cyan) is calculated as the product of the plasma-to-whole blood fraction (green) times the unchanged parent fraction (blue) curves, derived from the human data presented in Harms et al. [3]. (DOCX 61 kb)

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## About this article

### Cite this article

Wang, T., Wu, K.Y., Miner, R.C. *et al.* Reproducible quantification of cardiac sympathetic innervation using graphical modeling of carbon-11-meta-hydroxyephedrine kinetics with dynamic PET-CT imaging.
*EJNMMI Res* **8, **63 (2018). https://doi.org/10.1186/s13550-018-0421-5

Received:

Accepted:

Published:

### Keywords

- Logan
- MA1
- One tissue compartment
- Sympathetic nervous system
- HED