Advertisement
Research paper| Volume 17, ISSUE 3, P185-191, May 2023

Download started.

Ok

Quantification of myocardial blood flow using dynamic myocardial CT perfusion compared with 82Rb PET

Open AccessPublished:April 04, 2023DOI:https://doi.org/10.1016/j.jcct.2023.03.007

      Abstract

      Purpose

      Absolute measures of myocardial blood flow (MBF) obtained with dynamic myocardial CT perfusion (DM-CTP) are underestimated when compared with reference standards. This is to some extent explained by incomplete extraction of iodinated contrast agent (iCA) to the myocardial tissue. We aimed to establish an extraction function for iCA, use the function to calculate MBFCT and to compare this with MBF measured with 82Rb positron emission tomography (PET).

      Materials and methods

      Healthy individuals without coronary artery disease (CAD) were examined with 82Rb PET and DM-CTP. The factors a and β of the generalized Renkin-Crone model were estimated using a non-linear least squares model. The factors providing the best fit for the data were subsequently used to calculate MBFCT.

      Results

      Of consecutive 91 individuals examined, 79 were eligible for analysis. The factors a and β providing the best fit of the nonlinear least-squares model to the data were a ​= ​0.614 and β ​= ​0.218 (R-squared ​= ​0.81). Conversion of the CT inflow parameter (K1) values using the derived extraction function resulted in a significant correlation between MBF measured during stress using CT and PET (P ​= ​0.039).

      Conclusion

      In healthy individuals, flow estimates obtained with dynamic myocardial CT perfusion during stress were, after conversion to MBF using the extraction of iodinated CT contrast agent, correlated with absolute MBF quantified with 82Rb PET.

      Keywords

      Abbreviations:

      CAD (coronary artery disease), CCTA (coronary computed tomography angiography), DM-CTP (dynamic myocardial CT perfusion), iCA (iodinated contrast agent), MBF (myocardial blood flow), PET (positron emission tomography)

      1. Introduction

      Coronary computed tomography angiography (CCTA) has become a mainstay of coronary artery disease (CAD) evaluation and is now recommended as the first-line test for chronic coronary syndromes in clinical guidelines.
      • Gulati M.
      • Levy P.D.
      • Mukherjee D.
      • et al.
      AHA/ACC/ASE/CHEST/SAEM/SCCT/SCMR guideline for the evaluation and diagnosis of chest pain: executive summary: a report of the American college of cardiology/American heart association joint committee on clinical practice guidelines.
      • Knuuti J.
      • Wijns W.
      • Saraste A.
      • et al.
      ESC guidelines for the diagnosis and management of chronic coronary syndromes.
      • Moss A.J.
      • Williams M.C.
      • Newby D.E.
      • Nicol E.D.
      The updated NICE guidelines: cardiac CT as the first-line test for coronary artery disease.
      However, symptomatic myocardial ischemia leading to ischemic heart disease (IHD) can also exist in absence of obstructive CAD and can be caused by other structural and functional alterations such as microvascular dysfunction. Therefore, methods integrating evaluation of myocardial perfusion are valuable when selecting patients with suspected IHD who may benefit from treatment.
      • Tonino P.A.
      • De Bruyne B.
      • Pijls N.H.
      • et al.
      Fractional flow reserve versus angiography for guiding percutaneous coronary intervention.
      ,
      • Dewey M.
      • Siebes M.
      • Kachelriess M.
      • et al.
      Clinical quantitative cardiac imaging for the assessment of myocardial ischaemia.
      Impaired myocardial perfusion resulting in ischemia can be quantified objectively as myocardial blood flow (MBF). The gold standard for animal experimental estimation of MBF is microspheres and for non-experimental use, positron emission tomography (PET) is established as the clinical gold standard.
      Nevertheless, other imaging modalities such as echocardiography, single photon emission computed tomography, magnetic resonance imaging and cardiac CT are also capable of estimating both relative distribution of flow and absolute MBF.
      • Dewey M.
      • Siebes M.
      • Kachelriess M.
      • et al.
      Clinical quantitative cardiac imaging for the assessment of myocardial ischaemia.
      Among these, dynamic myocardial CT perfusion (DM-CTP) has emerged as a method that in combination with CCTA can provide anatomic and functional evaluation using a single modality.
      • Patel A.R.
      • Bamberg F.
      • Branch K.
      • et al.
      Society of cardiovascular computed tomography expert consensus document on myocardial computed tomography perfusion imaging.
      The extent of ischemia determined by DM-CTP can be defined using a relative cutoff value, which in general perform better compared with a cutoff based on absolute results.
      • Kono A.K.
      • Coenen A.
      • Lubbers M.
      • et al.
      Relative myocardial blood flow by dynamic computed tomographic perfusion imaging predicts hemodynamic significance of coronary stenosis better than absolute blood flow.
      Nonetheless, precise quantification of absolute MBF may be important, especially in case of potential three-vessel disease or suspected microvascular dysfunction.
      • Camici P.G.
      • Crea F.
      Coronary microvascular dysfunction.
      Quantification of MBF using DM-CTP is accomplished by the acquisition of multiple images as blood mixed with contrast moves through the myocardium.
      • Ho K.T.
      • Chua K.C.
      • Klotz E.
      • Panknin C.
      Stress and rest dynamic myocardial perfusion imaging by evaluation of complete time-attenuation curves with dual-source CT.
      Some of the most recent larger DM-CTP accuracy studies have reported absolute MBF values which were substantially lower than the results reported using microspheres and 15O PET.
      • Cobb F.R.
      • Bache R.J.
      • Greenfield Jr., J.C.
      Regional myocardial blood flow in awake dogs.
      • Danad I.
      • Uusitalo V.
      • Kero T.
      • et al.
      Quantitative assessment of myocardial perfusion in the detection of significant coronary artery disease: cutoff values and diagnostic accuracy of quantitative [(15)O]H2O PET imaging.
      • Coenen A.
      • Rossi A.
      • Lubbers M.M.
      • et al.
      Integrating CT myocardial perfusion and CT-FFR in the work-up of coronary artery disease.
      • Nous F.M.A.
      • Geisler T.
      • Kruk M.B.P.
      • et al.
      Dynamic myocardial perfusion CT for the detection of hemodynamically significant coronary artery disease.
      It has been suggested that “MBF” calculated as the ratio between the maximum slope of the fit myocardial curve and the peak arterial input function is likely to be computing the inflow rate constant K1 and not true MBF.
      • Ishida M.
      • Kitagawa K.
      • Ichihara T.
      • et al.
      Underestimation of myocardial blood flow by dynamic perfusion CT: explanations by two-compartment model analysis and limited temporal sampling of dynamic CT.
      K1 describes the inflow of tracer from capillaries to tissue and in order to convert the inflow rate to MBF, the extraction of iodinated contrast agent (iCA) is required. The extraction rate describes the proportion of tracer in the capillary blood which is extracted from the blood to the myocardial tissue. Iodinated contrast agent, as well as other small molecules used as tracers, display a distinct roll-off phenomenon where extraction decrease with increasing flow. To account for the incomplete and nonlinear extraction of tracer, an extraction function E(MBF) is needed to convert K1 to MBF.
      The present study aimed to establish a relation between flow and extraction of CT iCA in healthy individuals without cardiovascular risk factors and with CCTA verified absence of any CAD. In addition, we aimed to use the derived extraction function to calculate MBFCT and to compare this with MBF measured with 82Rb PET.

      2. Methods

      2.1 Study population and procedures

      This study (The DYNAMITE-PILOT study) was a prospective, single-center observational study performed from October 2020 to December 2021. The study was conducted according to the Declaration of Helsinki II and was approved by The Danish Committee on Health Research Ethics (protocol number H-19043153). Oral and written consent was obtained from all participants. Participants in the DYNAMITE-PILOT study were recruited from the ongoing Copenhagen General Population CT study, in which individuals from the Greater Copenhagen Area are undergoing a research protocol consisting of a non-contrast and a contrast Cardiac CT.
      • Sigvardsen P.E.
      • Fuchs A.
      • Kuhl J.T.
      • et al.
      Left ventricular trabeculation and major adverse cardiovascular events: the Copenhagen General Population Study.
      Individuals with a calcium score of zero, a normal CCTA and age >45 years were eligible for inclusion in the DYNAMITE-PILOT study.
      Specific exclusion criteria of the DYNAMITE-PILOT study included: 1) BMI >35, 2) diabetes, 3) hypertension, 4) hypercholesterolemia, 5) moderate or heavy smoking (>20 pack years), 6) inability to receive study information, 7) bronchospastic lung disease (COPD or asthma), and 8) pregnancy or active breastfeeding.
      Participants were assigned for two separate study procedures according to radiation safety constraints. Participants either underwent 82Rb PET (rest and stress) and DM-CTP examination during adenosine stress (PET and DM-CTP group), or DM-CTP at rest and during stress (DM-CTP only group). Paired flow measures were thus obtained for stress 82Rb PET and stress DM-CTP whereas no paired rest flow measures were acquired.
      Both paired results during stress and the unpaired results from examinations at rest were included in the derivation of the extraction function for iCA. The extraction as a function of MBF, E(MBF) was applied to the DM-CTP flow measures to obtain MBFCT, which were then compared with MBFPET.

      2.2 DM-CTP image acquisition

      Image acquisition was performed using a 320 detector row CT scanner (Aquilion ONE/PRISM Edition, Canon Medical, Otawara, Japan). Participants were instructed to strictly avoid caffeine 16 ​h prior to examination. Two intravenous lines were inserted, 12-lead ECG was obtained, and oral Beta blockers were administered to participants with a pre-examination heart rate above 65. An individualized contrast volume based on estimated body water was injected at a flow rate of 6 ​mL/s. Iodinated contrast agent with an iodine content of 350 ​mg I/ml and a molecular size of 1.2 ​nm was used (Omnipaque (Iohexol), GE Healthcare, Oslo, Norway).
      After a delay of 6 ​s, low-dose prospective acquisitions were performed at every other heartbeat for 30 consecutive seconds covering the base of the heart to the apex.
      • Moller M.B.
      • Sorgaard M.H.
      • Linde J.J.
      • Kober L.V.
      • Kofoed K.F.
      Optimization of image sampling rate to lower the radiation dose of dynamic myocardial CT perfusion.
      Tube voltage was 80 ​kV and tube current was fixed at 300 ​mA. The scan range was 120 ​mm covering the entire left ventricle for each rotation, the gantry rotation time was 275 ​ms and detector collimation was 0.5 ​mm ​× ​320. Dynamic myocardial CT perfusion was performed as an electrocardiography-triggered prospective target scan acquired at 77–87% of the RR interval. When heart rates above 70 were observed, the scan time was reduced to 25 ​s to reduce radiation dose.
      Adenosine (0.14 ​mg/kg/min) was used to induce myocardial hyperemia and the stress acquisitions were initiated 3 ​min after starting the infusion pump. Blood pressure and heart rate was recorded before and 2 ​min after initiating the infusion of adenosine. The rate pressure product (RPP) was calculated as the product of heart rate and systolic blood pressure. The RPP was calculated for all DM-CTP and PET studies at rest and during stress.
      In the DM-CTP only group, the stress scan was preceded by an examination at rest performed according to the same scan protocol as the stress. In this group, a wait of 15 ​min between the first and the second injection of iCA was ensured.

      2.3 DM-CTP evaluation and quantification of myocardial blood flow

      For the evaluation of rest and stress DM-CTP a low motion phase was determined by an automatic raw data analysis tool (phase eXact, Canon Medical Systems Corporation). Images were reconstructed at 1-mm slice thickness using Iterative dose reduction (Adaptive Iterative Dose Reduction (AIDR) 3D) and a kernel (FC03) with beam hardening correction. Before evaluation, all acquisitions were corrected for breathing related displacement using a non-rigid body registration algorithm. Also, a four-dimensional noise reduction filter (4D similarity filter) was applied.
      Interpretation of the image-data and calculation of MBF was performed on a workstation using dedicated software (Vitrea research, version 7.11, Vital Images, a Canon group company, Minnetonka, MN). The semi-automated software set the axis of the left ventricle, identified the contours of the myocardium, and placed a region of interest in the ascending aorta used to obtain the arterial input function. Manual augmentation of the contours was applied when necessary. The left ventricular myocardial tissue was subdivided into 16 segments (excluding the apical cap) according to the American Heart Association myocardial segmentation model. The 16 segments were divided into a total of 864 subsegments. The MBFCT was calculated according to the method described below for every subsegment.
      • Cerqueira M.D.
      • Weissman N.J.
      • Dilsizian V.
      • et al.
      Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association.

      2.4 Quantification of myocardial blood flow

      2.4.1 Signal and concentration of contrast

      The CT signal is measured in Hounsfield units (HU) and represents the combination of attenuation due to baseline attenuation without contrast and attenuation due to iCA.
      HUmeasured(t)=HUbaseline+HUcontrast(t)


      Baseline attenuation was measured at the first acquisition before any new contrast had entered and was subtracted from all subsequent measurements.
      HUcontrast(t)=HUmeasuredHUbaseline(t)


      Using CT, there is a specific linear relationship between signal (HU) and iodine concentration (mg I/ml) at a certain tube voltage (kV).
      Kikuchi et al. performed a phantom study using 80 ​kV and found the following relationship that was adapted in this study.
      • Kikuchi Y.
      • Oyama-Manabe N.
      • Naya M.
      • et al.
      Quantification of myocardial blood flow using dynamic 320-row multi-detector CT as compared with (1)(5)O-H(2)O PET.
      Ccontrast(t)=HUcontrast(t)47.9


      2.4.2 Modeling and quantifying flow rates

      The flow kinetics within the myocardium were described using a single tissue (two compartment) model developed by Kety and Schmidt together with its assumptions.
      • Kety S.S.
      • Schmidt C.F.
      The nitrous oxide method for the quantitative determination of cerebral blood flow in man: theory, procedure and normal values.
      This is given by the differential equation:
      dCt(t)dt=K1Ca(t)k2Ct(t)


      The derivative of the tissue concentration Ct(t) describes the change in the tissue concentration over time. This is determined by the flow into the tissue (the inflow rate constant K1 multiplied by the concentration of contrast in the arterial blood at the site of extraction (capillaries) at time t) subtracted by the outflow from the tissue (the outflow rate constant k2 multiplied by the concentration of contrast in the tissue at time t).
      K1 and k2 were estimated by the non-linear least-squares method where values of K1 and k2 were initially assumed. Then optimizations were performed to minimize the least squares error between the result on the right side of the equation and the slope of the tissue curve for all instances of t (left side).

      2.4.3 Calculation of myocardial blood flow

      Only a fraction of the contrast flowing through the capillaries is extracted to the extravascular compartment.
      • Canty Jr., J.M.
      • Judd R.M.
      • Brody A.S.
      • Klocke F.J.
      First-pass entry of nonionic contrast agent into the myocardial extravascular space. Effects on radiographic estimates of transit time and blood volume.
      In order to quantitate perfusion, the estimated K1 must therefore be divided by the extraction fraction (E) in order to quantitate perfusion.
      MBF=K1E
      (1)


      The extraction fraction is defined by an extraction function that is unique for a particular tracer. It is a function of MBF and is dependent on the product of permeability (P) and surface area (S). The permeability-surface area product (PS) describes the flow of molecules through the capillary membranes in a certain volume of tissue. The relation is described by the Renkin-Crone model.
      • Crone C.
      The permeability of capillaries in various organs as determined by use of the 'indicator diffusion' method.
      ,
      • Renkin E.M.
      Transport of potassium-42 from blood to tissue in isolated mammalian skeletal muscles.
      E(MBF)=1ePSMBF
      (2)


      This model is consistent with the observation that tracer extraction typically decreases with flow, despite the PS product increasing due to capillary recruitment.
      • Klein R.
      • Beanlands R.S.
      • deKemp R.A.
      Quantification of myocardial blood flow and flow reserve: technical aspects.
      The PS function is typically presented on the following form
      PS=αMBF+β
      (3)


      By combining equations ((2), (3)), the extraction function can be expressed as
      E(MBF)=1eαMBF+βMBF=1eαβMBF=1eαeβMBF
      (4)


      Substituting eα with a
      E(MBF)=1aeβMBF
      (5)


      This is inserted instead of E in equation (1) and K1 can be defined by the following expression
      K1=(1aeβMBF)MBF
      (6)


      The two extraction fraction parameters a and β were estimated for iCA as described in section 2.6. Then, Equations ((1), (5)) were combined and MBFCT values were quantified.
      MBFCT=K11aeβMBF
      (7)


      The myocardial flow reserve (MFR) was calculated as the ratio between the MBF during stress and at rest
      MFR=MBFStressMBFRest


      2.5 PET image acquisition, evaluation and quantification of MBF

      Perfusion imaging with DM-CTP and 82Rb PET was separated by a minimum of 48 ​h. All patients were asked to abstain from caffeine 12 ​h prior to imaging. PET myocardial perfusion imaging was performed during rest and stress conditions in a single session. For each condition, participants received 1110 MBq (±10%) 82Rb supplied from a CardioGen-82 Sr-82/Rb-82 generator manufactured for Bracco (Bracco Diagnostics Inc., Princeton, NJ).
      Rest and stress images were acquired electrocardiography-triggered for 6 ​min from the start of the 82Rb infusion on a Siemens Biograph mCT/PET 128-slice scanner (Siemens Healthcare, Knoxville, TN, USA). Adenosine (0.14 ​mg/kg/min) was used to induce myocardial hyperemia for 6 ​min, and the 82Rb infusion was initiated 2.5 ​min after starting the infusion of adenosine.
      • von Scholten B.J.
      • Hasbak P.
      • Christensen T.E.
      • et al.
      Cardiac (82)Rb PET/CT for fast and non-invasive assessment of microvascular function and structure in asymptomatic patients with type 2 diabetes.
      The 6 ​min long PET acquisitions were reconstructed into dynamic image-series consisting of 18 frames (1 ​× ​10s, 8 ​× ​5s, 3 ​× ​10s, 2 ​× ​20s, 4 ​× ​60s) using an iterative 3D ordered subset expectation maximization model (2 iterations and 21 subsets) with corrections for time-of-flight and point-spread function. All data were smoothened using a 6.5-mm Gaussian post-filtering.
      The left ventricular myocardial tissue was subdivided into 17 segments.
      • Cerqueira M.D.
      • Weissman N.J.
      • Dilsizian V.
      • et al.
      Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association.
      No manual adjustments were made. Pre-study evaluations showed inclusion of extramyocardial structures in the basal segments and segments 1–6 were therefore excluded from analysis. In addition, segment 17 was also excluded from the analysis to comply with the DM-CTP evaluation protocol. The 82Rb PET MBF values (MBFPET) were calculated according to the Lortie model using dedicated software (QPET, Cedars-Sinai Medical Center, Los Angeles, California, USA).
      • Lortie M.
      • Beanlands R.S.
      • Yoshinaga K.
      • Klein R.
      • Dasilva J.N.
      • DeKemp R.A.
      Quantification of myocardial blood flow with 82Rb dynamic PET imaging.

      2.6 Derivation of extraction function and calculation of MBF

      Both flow measures obtained at rest and during stress were used for derivation of the extraction function. Since the rest observations were unpaired, the rest observations were divided in three groups of ten DM-CTP observations in each based on the size of the RPP, which was from the literature expected to be related to MBF at rest.
      • Czernin J.
      • Muller P.
      • Chan S.
      • et al.
      Influence of age and hemodynamics on myocardial blood flow and flow reserve.
      The three groups were matched with the same RPP interval measured at the 82Rb PET rest examination. The mean of the three RPP group values given the weight equivalent to the number of CT observations were utilized as data points for the derivation of the extraction function in addition to the stress results.
      The parameters a and β of the generalized Renkin-Crone model in Equation (6) were estimated using a non-linear least squares model. Values of a and β providing the best fit for the data were used to calculate MBFCT by solving Equation (7).
      Then, the derived extraction function was applied to the DM-CTP flow measures (K1 at rest and during stress) in order to obtain MBFCT (Equation (7)), which were finally compared with absolute flow quantitated with 82Rb PET (MBFPET).

      2.7 Comparison of extraction with existing literature

      In order to assess the derived extraction function, the extraction fraction (E) of iCA was calculated for the mean MBFCT at rest and during stress of the present study as well as for theoretical MBF values at rest, MBF ​= ​1 ​mL/min/g and during stress, 4 ​= ​mL/min/g. The current literature was reviewed for research describing the extraction of iodinated contrast from blood to the myocardial tissue. Two clinical and one animal experimental study were identified (Table 1).
      • Kikuchi Y.
      • Oyama-Manabe N.
      • Naya M.
      • et al.
      Quantification of myocardial blood flow using dynamic 320-row multi-detector CT as compared with (1)(5)O-H(2)O PET.
      ,
      • Canty Jr., J.M.
      • Judd R.M.
      • Brody A.S.
      • Klocke F.J.
      First-pass entry of nonionic contrast agent into the myocardial extravascular space. Effects on radiographic estimates of transit time and blood volume.
      ,
      • Takafuji M.
      • Kitagawa K.
      • Ishida M.
      • et al.
      Clinical validation of the accuracy of absolute myocardial blood flow quantification with dual-source CT using (15)O-water PET.
      Since no experimental data describing the extraction at rest were identified, two studies on the fractional extraction of comparable molecules were included.
      • Svendsen J.H.
      • Bjerrum P.J.
      • Haunso S.
      Myocardial capillary permeability after regional ischemia and reperfusion in the in vivo canine heart. Effect of superoxide dismutase.
      ,
      • Svendsen J.H.
      • Efsen F.
      • Haunso S.
      Capillary permeability of 99 mTc-DTPA and blood flow rate in the human myocardium determined by intracoronary bolus injection and residue detection.
      Estimates for parameters a and β acquired from the literature were used to generate extraction functions, on the generalized Renkin Crone formula, for graphical comparison.
      Table 1Extraction (E) of iodinated CT contrast agent and other small molecules.
      TracerE at restE at MBF ​= ​1 ​mL/min/gE during stressE at MBF ​= ​4 ​mL/min/g
      Experimental estimation
       Canty et al., 1991CT CA (Ioversol)0.33 ​± ​0.02
       Svendsen et al., 1991Small molecules
      Small hydrophilic molecules with comparable size and properties as iodinated CT contrast agent. MBF at rest and stress are presented ± 1SD. For MBF ​= ​1 and 4, no SD's are available since the value is derived from the extraction function.
      0.46 ​± ​0.04
       Svendsen et al., 199299mTc-DTPA0.55 ​± ​0.09
      Reference standard estimation
       Kikuchi et al., 2014CT CA (Iohexol)0.730.33
       Takafuji et al., 2021CT CA (iodinated)0.810.39
      All studies are included in the reference list. E ​= ​1 and 4 are approximate reference flows at rest and during stress for healthy individuals. Studies investigating extraction of CT CA are highlighted in grey.
      CA: contrast agent; DTPA: diethylenetriamine pentaacetate; EDTA: ethylenediaminetetraacetic acid; MBF: myocardial blood flow.
      a Small hydrophilic molecules with comparable size and properties as iodinated CT contrast agent. MBF at rest and stress are presented ± 1SD. For MBF ​= ​1 and 4, no SD's are available since the value is derived from the extraction function.

      2.8 Statistical analysis

      Using descriptive statistics, the continuous variables are expressed as mean ​± ​SD and the categorical data are expressed as frequency (percentage). The Mann-Whitney U test (Wilcoxon rank-sum test) was used for comparison of non-parametric data and the students t-test was used to asses differences in parametric data across groups. Sex and family history of cardiovascular disease are expressed as frequency (percentage) and the Fisher's exact test was used to assess differences in proportion. The relationship between continuous variables was assessed using the Pearson correlation coefficient. P-values of <0.05 were considered significant. All statistical analyses were performed using R statistical package, version 3.6.1 (R Foundation for Statistical Computing).

      3. Results

      Due to the COVID-19 pandemic, the planned enrollment of 100 participants were hindered and a total of 91 healthy individuals were enrolled and underwent either rest and stress 82Rb PET and stress DM-CTP (PET and DM-CTP group, n ​= ​57) or DM-CTP only including both rest and stress (DM-CTP only, n ​= ​34). Of the consecutive 91 participants examined, 79 were eligible for analysis (Fig. 1). Twelve participants were excluded due to technical failure (n ​= ​4), structural heart disease discovered after inclusion (n ​= ​3), PET examination not performed (n ​= ​3), CA extravasation (n ​= ​1) and discomfort (severe dizziness) during adenosine (n ​= ​1)). There were no significant differences in sex, age, BMI or prior smoking between the PET and DM-CTP and DM-CTP only groups (Table 2). Participants in the PET and DM-CTP group had a median time between examinations of 14 days (range; 2–91 days).
      Fig. 1
      Fig. 1Participant flowchart. DM-CTP ​= ​dynamic myocardial CT perfusion, CAD ​= ​coronary artery disease, CCTA ​= ​coronary CT angiography, PET ​= ​positron emission tomography.
      Table 2Characteristics of the study population.
      PET and DM-CTPDM-CTP onlyP-value
      Baseline characteristics
       N4930
       Female27 (55)16 (53)1.00
       Age (years)56.29 ​± ​8.5657.27 ​± ​8.560.60
       BMI (kg/m2)24.44 ​± ​3.5123.86 ​± ​2.630.40
       Smoking (pack-years)2.55 ​± ​4.723.70 ​± ​4.700.76
       Family history of CVD6 (12)1 (3)0.24
      Data are presented as counts (percentage) or mean ​± ​standard deviation.
      CCTA: coronary CT angiography; CVD: Cardiovascular disease.

      3.1 Hemodynamics

      A significant increase in heart rate as well as RPP during stress conditions in the PET and DM-CTP group was observed for both PET (heart rate; 64 vs. 80 bpm, P ​< ​0.01 and RPP; 7485 vs. 9405 bpm · mmHg, P ​< ​0.01) and DM-CTP (heart rate; 58 vs. 75 bpm, P ​< ​0.01 and RPP; 7624 vs. 9591 bpm · mmHg, P ​< ​0.01). No significant differences in the RPP were observed between the CT (7624 ​± ​1363 ​at rest and 9591 ​± ​2446 during stress) and PET examination (7485 ​± ​1558 bpm · mmHg at rest and 9405 ​± ​1989 bpm · mmHg during stress). Between procedure reproducibility of the RPP during stress for all subjects of the PET and DM-CTP group was high (r ​= ​0.98) (Fig. 2).
      Fig. 2
      Fig. 2Rate pressure product at 82Rb PET and DM-CTP. PET ​= ​positron emission tomography, RPP ​= ​rate pressure product.

      3.2 Extraction function for iodinated contrast

      In Fig. 3, the K1 estimates obtained with DM-CTP are plotted against the absolute MBF obtained with 82Rb PET. K1 was 0.48 ​± ​0.17 ​mL/min/g at rest and 1.50 ​± ​0.45 ​mL/min/g during stress. The factors a and β providing the best fit of the nonlinear least-squares model (see equation (6)) to the data was a ​= ​0.61 and β ​= ​0.22 (R-squared ​= ​0.81).
      Fig. 3
      Fig. 3K1 obtained with DM-CTP and MBF obtained with 82Rb PET. Solid black line represents the nonlinear least squares model with equation K1 = (1–0.61 × e−0.22/MBF) × MBF with best fit to the data (a ​= ​0.614 and β ​= ​0.218).
      The resulting extraction function E ​= ​1–0.61 × e−0.22/MBF lead to a mean extraction fraction at rest of 0.52 ​± ​0.05 and 0.43 ​± ​0.03 during stress. The extraction function derived in the present study is plotted in Fig. 4 along with the extraction functions derived by the two clinical studies by Kikuchi et al. and Takafuji et al.
      • Kikuchi Y.
      • Oyama-Manabe N.
      • Naya M.
      • et al.
      Quantification of myocardial blood flow using dynamic 320-row multi-detector CT as compared with (1)(5)O-H(2)O PET.
      ,
      • Takafuji M.
      • Kitagawa K.
      • Ishida M.
      • et al.
      Clinical validation of the accuracy of absolute myocardial blood flow quantification with dual-source CT using (15)O-water PET.
      Fig. 4
      Fig. 4Two representations of the extraction function. A, Extraction fraction as a function of MBF. MBF ​= ​0.5 to 5.0 and extraction fraction of 0.31–0.64 represent physiological range. Blue dots represent extraction fraction estimates obtained by animal experiments at rest and during stress which are presented in (Rest ​= ​0.5 and stress ​= ​0.33). B, Relation between K1 and MBF. Grey zones represent DM-CTP K1 ​± ​one standard deviation obtained in the current study at rest and during stress. MBF ​= ​myocardial blood flow.

      3.3 Estimates of myocardial blood flow using DM-CTP

      A comparison of MBF estimated with 82Rb PET and MBFCT obtained by conversion of DM-CTP K1 values to MBFCT, using the derived extraction function and equation (7), are presented in Fig. 5. A significant correlation between MBF measured during stress using CT and PET was found (P ​= ​0.039).
      Fig. 5
      Fig. 5Myocardial blood flow during stress evaluated with dynamic myocardial CT perfusion and 82Rb PET. A significant correlation between MBF measured during stress using CT and PET was found (P ​= ​0.039). Dotted line represents identity line and solid line represents best fit to the data. Since no paired data w obtained at rest, the regression line was forced through the origin, which is physiologically rational. Slope ​= ​1.01 and R2 ​= ​0.90. R-squared should be interpreted with caution due to the exclusion of the intercept in the best fit. MBF ​= ​myocardial blood flow.
      A group comparison of MBF measured with the two modalities is included in the supplementary material.

      4. Discussion

      In the present study, an expression determining extraction as a function of MBF based on the generalized Renkin Crone model was established by using 82Rb PET as the reference standard; E ​= ​1–0.61 × e−0.22/MBF.
      The extraction of a tracer can be estimated indirectly, as shown in the present study, by using a reference representing “true” MBF but can also be estimated by direct measurement of arterial and venous concentration, although this is difficult to achieve in humans.
      Direct experimental estimation of extraction of iCA has been undertaken by Canty et al. in a very elegant study.
      • Canty Jr., J.M.
      • Judd R.M.
      • Brody A.S.
      • Klocke F.J.
      First-pass entry of nonionic contrast agent into the myocardial extravascular space. Effects on radiographic estimates of transit time and blood volume.
      They compared the CT myocardial signal of a conventional nonionic CA (ioversol) to a particulate emulsion (ethiodol) that remained in the vascular space in closed-chest dogs during pharmacological vasodilation and found an extraction fraction of 0.33 ​± ​0.02. The study was not performed at rest. However, they also reported extraction of 50 ​± ​3% in areas supplied by vessels with induced stenosis, which might be somewhat comparable to the resting condition since the hyperemic response is blunted or absent in areas with ischemia during stress.
      The indirect method has been utilized to establish the extraction of iCA by Kikuchi et al. and Takafuji et al. who determined the relation between extraction of iCA and MBF using the current non-invasive gold standard 15O PET and found E ​= ​1–0.90 × e−1.20/MBF and E ​= ​1–0.90 × e−1.54/MBF respectively (Fig. 4-A).
      • Kikuchi Y.
      • Oyama-Manabe N.
      • Naya M.
      • et al.
      Quantification of myocardial blood flow using dynamic 320-row multi-detector CT as compared with (1)(5)O-H(2)O PET.
      ,
      • Takafuji M.
      • Kitagawa K.
      • Ishida M.
      • et al.
      Clinical validation of the accuracy of absolute myocardial blood flow quantification with dual-source CT using (15)O-water PET.
      Characteristics separating these studies from the present are sample size (N ​= ​12 and N ​= ​17 vs. 49 matched and 30 unmatched individuals), CAD status (absence of significant CAD and known or suspected CAD vs. CCTA verified absence of any CAD in our study) and finally the reference method. These differences complicate the direct comparison. However, inspection of the extraction functions of the two studies reveals supraphysiological extraction values at MBF <1.5 ​mL/min/g when compared with experimental data. On the other hand, it seems that the extraction functions derived by Kikuchi et al. and Takafuji et al. reflect the experimental data better at MBF >3.0 ​mL/min/g compared with the extraction function of the present study.
      The consensus of the present and the two other clinical studies is that extraction of iCA to the myocardial tissue during first pass is incomplete. Assuming complete extraction results in MBFCT equal to K1 which would cause significant underestimation of MBFCT. In addition, the extraction fraction of iCA decrease with increasing flow and using a constant for conversion is therefore not optimal and would potentially result in overestimation of resting and underestimation of hyperemic flow.
      In the present study, we found a significant correlation between MBF measured during stress using CT and PET. This adds to the growing body of evidence regarding the accuracy of myocardial blood flow estimation from DM-CTP compared with PET.
      • Kikuchi Y.
      • Oyama-Manabe N.
      • Naya M.
      • et al.
      Quantification of myocardial blood flow using dynamic 320-row multi-detector CT as compared with (1)(5)O-H(2)O PET.
      ,
      • Takafuji M.
      • Kitagawa K.
      • Ishida M.
      • et al.
      Clinical validation of the accuracy of absolute myocardial blood flow quantification with dual-source CT using (15)O-water PET.
      ,
      • Alessio A.M.
      • Bindschadler M.
      • Busey J.M.
      • Shuman W.P.
      • Caldwell J.H.
      • Branch K.R.
      Accuracy of myocardial blood flow estimation from dynamic contrast-enhanced cardiac CT compared with PET.
      However, when inspecting the individual observations, results are scattered. A contributing factor could be the discrepancy between the blood flow kinetics visualized using the two modalities. Because of a short injection time of the CT tracer (5–10 ​s), a high temporal resolution and an acquisition time of 30 ​s, only the first pass is visualized. In comparison 82Rb is injected much slower, the temporal resolution is lower and the total acquisition time is 6 ​min. The kinetics visualized using PET are therefore different compared with DM-CTP.
      In addition, scattered results, represented by large distance to the best fit to the data, are also seen when comparisons within the same modality are undertaken. One example of similar scattered results is the study by Lortie et al. which compares 82Rb PET with 13N-ammonia PET.
      • Lortie M.
      • Beanlands R.S.
      • Yoshinaga K.
      • Klein R.
      • Dasilva J.N.
      • DeKemp R.A.
      Quantification of myocardial blood flow with 82Rb dynamic PET imaging.
      Nevertheless, the general impression is that the results obtained with PET is reproduced using DM-CTP, which is reflected by the significant correlation.

      4.1 Limitations

      The primary limitation of this study is the single center, single scanner, single perfusion CT software and single CA injection protocol setup. Apart from the injection protocol, the extraction is not dependent on these factors. Nevertheless, when using the indirect extraction function derivation method utilized in the present study, these factors will impact the result. In addition, the estimated extraction function is based on MBF obtained with the chosen reference standard which in this study was 82Rb PET. Despite the limited clinical implementation, 15O PET is still considered the gold standard for quantification. Another important limitation of this study is the lack of paired 82Rb PET and DM-CTP examinations at rest that would have improved the total number of paired observations and most likely, robustness of the results.

      4.2 Clinical implications of this study

      This study aims to increase the awareness that the inflow rate (K1) is significantly different from MBFCT and that an extraction function is necessary for absolute quantification of MBFCT if comparison of absolute numbers with other modalities is intended (Fig. 4-B).
      Since quantification depends on parameters mentioned in the limitations, the extraction function is ideally evaluated before undertaking clinical DM-CTP using site-specific set-up. However, this approach is not feasible for most sites and as an alternative, the extraction function presented in this study can be adopted and tedious work bypassed. This especially applies when similar set-up is utilized.
      Because a sampling rate of one acquisition for every second heartbeat was selected in the present study, the scan protocol is also applicable to dual-source CT scanners (that must alternate table position due to limited z-axis range) which has been employed in most of the previous DM-CTP studies. Before implementation, use of the presented extraction function should undergo clinical validation.

      5. Conclusions

      In this study, the first pass extraction of iodinated CT contrast agent to the myocardial tissue as a function of flow was defined in a cohort with coronary CT angiography verified absence of coronary artery disease. Flow estimates obtained with dynamic myocardial CT perfusion during stress were, after conversion to myocardial blood flow using the derived extraction function, correlated with absolute myocardial blood flow quantified with 82Rb positron emission tomography. Experimental estimation of the extraction of iodinated contrast agent in humans over a range of myocardial blood flows, which is unaffected by scanner, protocol and software, is warranted but difficult to achieve.

      Declaration of competing interest

      This study was supported by Sygeforsikring “danmark”; AP Møller og hustru Chastine McKinney Møllers Fond; The Danish Agency for Science Technology and Innovation; The Danish Council for Strategic Research, Copenhagen; The Lundbeck Foundation; The Research Council of Rigshopitalet; The Research Council of the Heart.
      Center, Rigshospitalet, Copenhagen; The Danish Research Foundation, Canon Medical Systems Corporation; The Danish Heart Foundation; The John and Birthe Meyer Foundation; and The University of Copenhagen.
      The funding sources did not have any role in the study design; conduct of the study; data analysis; data interpretation; or writing of this report.

      Appendix A. Supplementary data

      The following is/are the supplementary data to this article.

      References

        • Gulati M.
        • Levy P.D.
        • Mukherjee D.
        • et al.
        AHA/ACC/ASE/CHEST/SAEM/SCCT/SCMR guideline for the evaluation and diagnosis of chest pain: executive summary: a report of the American college of cardiology/American heart association joint committee on clinical practice guidelines.
        Circulation. 2021 Nov 30; 144: e368-e454https://doi.org/10.1161/CIR.0000000000001029
        • Knuuti J.
        • Wijns W.
        • Saraste A.
        • et al.
        ESC guidelines for the diagnosis and management of chronic coronary syndromes.
        Eur Heart J. 2019; (2019)
        • Moss A.J.
        • Williams M.C.
        • Newby D.E.
        • Nicol E.D.
        The updated NICE guidelines: cardiac CT as the first-line test for coronary artery disease.
        Curr Cardiovasc Imaging Rep. 2017; 10: 15
        • Tonino P.A.
        • De Bruyne B.
        • Pijls N.H.
        • et al.
        Fractional flow reserve versus angiography for guiding percutaneous coronary intervention.
        N Engl J Med. 2009; 360: 213-224
        • Dewey M.
        • Siebes M.
        • Kachelriess M.
        • et al.
        Clinical quantitative cardiac imaging for the assessment of myocardial ischaemia.
        Nat Rev Cardiol. 2020; 17: 427-450
        • Patel A.R.
        • Bamberg F.
        • Branch K.
        • et al.
        Society of cardiovascular computed tomography expert consensus document on myocardial computed tomography perfusion imaging.
        J Cardiovasc Comput Tomogr. 2020; 14: 87-100
        • Kono A.K.
        • Coenen A.
        • Lubbers M.
        • et al.
        Relative myocardial blood flow by dynamic computed tomographic perfusion imaging predicts hemodynamic significance of coronary stenosis better than absolute blood flow.
        Invest Radiol. 2014; 49: 801-807
        • Camici P.G.
        • Crea F.
        Coronary microvascular dysfunction.
        N Engl J Med. 2007; 356: 830-840
        • Ho K.T.
        • Chua K.C.
        • Klotz E.
        • Panknin C.
        Stress and rest dynamic myocardial perfusion imaging by evaluation of complete time-attenuation curves with dual-source CT.
        JACC Cardiovasc Imaging. 2010; 3: 811-820
        • Cobb F.R.
        • Bache R.J.
        • Greenfield Jr., J.C.
        Regional myocardial blood flow in awake dogs.
        J Clin Invest. 1974; 53: 1618-1625
        • Danad I.
        • Uusitalo V.
        • Kero T.
        • et al.
        Quantitative assessment of myocardial perfusion in the detection of significant coronary artery disease: cutoff values and diagnostic accuracy of quantitative [(15)O]H2O PET imaging.
        J Am Coll Cardiol. 2014; 64: 1464-1475
        • Coenen A.
        • Rossi A.
        • Lubbers M.M.
        • et al.
        Integrating CT myocardial perfusion and CT-FFR in the work-up of coronary artery disease.
        JACC Cardiovasc Imaging. 2017; 10: 760-770
        • Nous F.M.A.
        • Geisler T.
        • Kruk M.B.P.
        • et al.
        Dynamic myocardial perfusion CT for the detection of hemodynamically significant coronary artery disease.
        JACC Cardiovasc Imaging. 2021 Jan; 15: 75-87https://doi.org/10.1016/j.jcmg.2021.07.021
        • Ishida M.
        • Kitagawa K.
        • Ichihara T.
        • et al.
        Underestimation of myocardial blood flow by dynamic perfusion CT: explanations by two-compartment model analysis and limited temporal sampling of dynamic CT.
        J Cardiovasc Comput Tomogr. 2016; 10: 207-214
        • Sigvardsen P.E.
        • Fuchs A.
        • Kuhl J.T.
        • et al.
        Left ventricular trabeculation and major adverse cardiovascular events: the Copenhagen General Population Study.
        Eur Heart J Cardiovasc Imaging. 2021; 22: 67-74
        • Moller M.B.
        • Sorgaard M.H.
        • Linde J.J.
        • Kober L.V.
        • Kofoed K.F.
        Optimization of image sampling rate to lower the radiation dose of dynamic myocardial CT perfusion.
        J Cardiovasc Comput Tomogr. 2021 Sep-Oct; 15: 457-460https://doi.org/10.1016/j.jcct.2021.04.001
        • Cerqueira M.D.
        • Weissman N.J.
        • Dilsizian V.
        • et al.
        Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A statement for healthcare professionals from the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association.
        Circulation. 2002; 105: 539-542
        • Kikuchi Y.
        • Oyama-Manabe N.
        • Naya M.
        • et al.
        Quantification of myocardial blood flow using dynamic 320-row multi-detector CT as compared with (1)(5)O-H(2)O PET.
        Eur Radiol. 2014; 24: 1547-1556
        • Kety S.S.
        • Schmidt C.F.
        The nitrous oxide method for the quantitative determination of cerebral blood flow in man: theory, procedure and normal values.
        J Clin Invest. 1948; 27: 476-483
        • Canty Jr., J.M.
        • Judd R.M.
        • Brody A.S.
        • Klocke F.J.
        First-pass entry of nonionic contrast agent into the myocardial extravascular space. Effects on radiographic estimates of transit time and blood volume.
        Circulation. 1991; 84: 2071-2078
        • Crone C.
        The permeability of capillaries in various organs as determined by use of the 'indicator diffusion' method.
        Acta Physiol Scand. 1963; 58: 292-305
        • Renkin E.M.
        Transport of potassium-42 from blood to tissue in isolated mammalian skeletal muscles.
        Am J Physiol. 1959; 197: 1205-1210
        • Klein R.
        • Beanlands R.S.
        • deKemp R.A.
        Quantification of myocardial blood flow and flow reserve: technical aspects.
        J Nucl Cardiol. 2010; 17: 555-570
        • von Scholten B.J.
        • Hasbak P.
        • Christensen T.E.
        • et al.
        Cardiac (82)Rb PET/CT for fast and non-invasive assessment of microvascular function and structure in asymptomatic patients with type 2 diabetes.
        Diabetologia. 2016; 59: 371-378
        • Lortie M.
        • Beanlands R.S.
        • Yoshinaga K.
        • Klein R.
        • Dasilva J.N.
        • DeKemp R.A.
        Quantification of myocardial blood flow with 82Rb dynamic PET imaging.
        Eur J Nucl Med Mol Imag. 2007; 34: 1765-1774
        • Czernin J.
        • Muller P.
        • Chan S.
        • et al.
        Influence of age and hemodynamics on myocardial blood flow and flow reserve.
        Circulation. 1993; 88: 62-69
        • Takafuji M.
        • Kitagawa K.
        • Ishida M.
        • et al.
        Clinical validation of the accuracy of absolute myocardial blood flow quantification with dual-source CT using (15)O-water PET.
        Radiol Cardiothorac Imaging. 2021; 3e210060
        • Svendsen J.H.
        • Bjerrum P.J.
        • Haunso S.
        Myocardial capillary permeability after regional ischemia and reperfusion in the in vivo canine heart. Effect of superoxide dismutase.
        Circ Res. 1991; 68: 174-184
        • Svendsen J.H.
        • Efsen F.
        • Haunso S.
        Capillary permeability of 99 mTc-DTPA and blood flow rate in the human myocardium determined by intracoronary bolus injection and residue detection.
        Cardiology. 1992; 80: 18-27
        • Alessio A.M.
        • Bindschadler M.
        • Busey J.M.
        • Shuman W.P.
        • Caldwell J.H.
        • Branch K.R.
        Accuracy of myocardial blood flow estimation from dynamic contrast-enhanced cardiac CT compared with PET.
        Circ Cardiovasc Imaging. 2019; 12e008323