Dynamic Quantitative Iodine Myocardial Perfusion Imaging with Dual-Layer CT using a Porcine Model

Ischemic heart disease is the globally leading cause of death. When using coronary CT angiography, the functional hemodynamics within the myocardium remain uncertain. In this study myocardial CT perfusion imaging using iodine contrast agent demonstrated to strongly improve the assessment of myocardial disorders. However, a retrieval of such dynamics using Hounsfield units from conventional CT poses concerns with respect to beam-hardening effects and low contrast-to-noise ratio (CNR). Dual-energy CT offers novel approaches to overcome aforementioned limitations. Quantitative peak enhancement, perfusion, time to peak and iodine volume measurements inside the myocardium were determined resulting in 0.92 mg/ml, 0.085 mg/ml/s 17.12 s and 29.89 mg/ml*s, respectively. We report on the first extensive quantitative and iodine-based analysis of myocardial dynamics in a healthy porcine model using a dual-layer spectral CT. We further elucidate on the potential of reducing the radiation dose from 135 to 18 mGy and the contrast agent volume from 60 to 30 mL by presenting a two-shot acquisition approach and measuring iodine concentrations in the myocardium in-vivo down to 1 mg/ml, respectively. We believe that dynamic quantitative iodine perfusion imaging may be a highly sensitive tool for the precise functional assessment and monitoring of early myocardial ischemia.

helps to avoid movement artifacts. In general, a main advantage of CT is its high resolution. Considering the relatively small thickness of the myocardium, this has a high potential to improve the detection of even small regions of subendocardial ischemia. Multiple studies showed evidence that myocardial perfusion imaging can significantly increase the accuracy of assessing flow limiting coronary artery diseases or areas of prior infarction [6][7][8] . Especially in the case of moderate coronary artery disease, where an anatomical evaluation of the stenosis is ambiguous, perfusion imaging is considered helpful for diagnosis.
Among CT myocardial perfusion (CTMP) imaging two approaches have been developed. With the static technique a single scan is performed when the enrichment of iodine within the myocardium is considered maximal, so that an optimal contrast between healthy, under-perfused and infarcted myocardial tissue is achieved. To ensure correct timing, the (test) bolus has to be tracked and the enrichment within the aorta determined 9 . The major advantage of static perfusion imaging is that myocardial blood-flow can be rapidly assessed, and that the underlying radiation dose is comparatively low. Using rest/stress protocols, especially in combination with CTA, a good assessment of suspected coronary artery disease has been reported 10,11 . However, static CTMP only provides a (semi-)qualitative analysis, while also being highly susceptible to temporal perfusion artifacts and thereby heavily reliant on the correct timing. In comparison, dynamic CTMP is performed by constantly acquiring scans within a 20-40 seconds timeframe and hence imaging the complete first pass of contrast agent through the heart 12 . While this technique is rigid and allows for a quantitative analysis of hemodynamics, the underlying Hounsfield-based assessment is susceptible to beam-hardening effects induced by the iodine contrast agent 13 . Further physiological models must be considered which relate Hounsfield units to the actual blood or iodine flow 14 . Finally, the continuous acquisition mode poses concerns with respect to an increased radiation dose.
Here the introduction of dual-energy CT (DECT) in clinical routine offers potential to overcome the aforementioned limitations. Different technological approaches have been established for acquiring dual-energy data. Mainly three technical implementations are used in clinical CT scanners: Dual-source CT (DSCT), rapid kVp switching CT (KVSCT) and Dual-layer CT (DLCT). A DSCT system consists of two X-ray tubes and two associated detectors, whereas the KVSCT machines use one detector and a source with two alternating kVp settings during the scan. To acquire spectral information in DSCT and KVSCT systems a scan must be designed accordingly before its start. As the photon energy separation takes place in the detector, in DLCT conventional and spectral data are obtained simultaneously without the necessity of a spectral scan pre-setting. Spectral data provides strongly increased contrast and sensitivity (via the virtual monoenergetic image (MonoE) and iodine density image channels) towards soft-tissue and iodine-based imaging, thereby opening up new perspectives for a reduction of radiation dose and volume of injected contrast agent. In 2012 So. et al. showed that by using rapid kVp switching DECT a higher reproducibility of dynamic myocardial perfusion measurements in comparison to a Hounsfield units based evaluation can be assured 15 . Further Sánchez-Gracián et al. demonstrated, that iodine quantification in myocardial perfusion stress DECT is beneficial for the differentiation of healthy and ischemic or necrotic myocardium 16 . Finally, Fahmi et al. reported on quantitative myocardial perfusion images and thereby demonstrated a superior assessment of myocardial infarcts in comparison to a conventional measurement 17 . Sellerer et al. showed an improved root mean square deviation (RMSD) of observed iodine concentrations (with respect to true values) for different measurement configurations using a modern clinical DLCT system 18 . This is of particular importance, as the quantification accuracy of iodine densities in regions of low contrast agent accumulation like in the myocardium plays a crucial role in the generation of our results.
In this study, we report on the very first extensive quantitative and iodine-based analysis of myocardial dynamics in a healthy porcine model in rest, using a dual-layer spectral CT. We include the presentation of iodine perfusion, time to peak, peak enhancement and volume maps. In a first step, we show in congruence with previous studies, that low-energy virtual MonoE-keV images and iodine density maps are superior in comparison to conventional images in the depiction of temporal opacification of the myocardium with respect to relative signal enhancement and contrast-to-noise ratio. Afterwards the temporal uptake of iodine is modelled using a gamma-variate function and corresponding dynamics deduced from the model in a pixel-wise manner. Further, we elucidate on the possibility and performance of quantitatively assessing the myocardial peak enhancement map using a semi-static two-shot approach. Based on the accuracy and consistency of determining low iodine densities with modern dual-energy systems, we suggest that dynamic quantitative iodine perfusion imaging may be a highly sensitive tool for the precise functional assessment and monitoring of early myocardial ischemia, if a valid database for stress and/or rest can be established 18 .

Results
A healthy pig (31.75 cm in mean diameter, 78 kg) was scanned using a state of the art dual-layer spectral CT system from Philips with a detector coverage of 4 cm and a peak voltage of 120 kVp. An overall volume of iodine contrast agent of 40 ml (Ultravist 300, Bayer, Bayer AG, Leverkusen, Germany) was injected into the ear vein, which corresponds to a concentration of 154 mg/kg-bodyweight. Within this animal study, scans were immediately acquired after injection to track the bolus and enable hemodynamics analysis throughout the whole thorax. Electrocardiographically (ECG)-triggered dynamic axial scans were adjusted to the cardiac cycle with an increment of zero, with the scanning area covering the mid part of the left ventricular myocardium. A total of 36 consecutive scans were acquired. Finally, no (commercial) automatic motion-correction algorithm for the heart was available as the entire heart volume was imaged only fragmentary. Therefore, data points with extensive motion in heart or other thorax structures were manually abandoned, resulting in a non-equidistant sampling of the curve with 15 scan points. Figure 1a shows the corresponding Hounsfield units for conventional and mono-energetic reconstruction as well as the iodine density of the porcine myocardium at the initial stage and at high contrast agent enrichment, i.e. 23.3 seconds after the contrast agent injection. In congruence with previous dual-energy studies on other indications, the temporal relative signal increase during iodine uptake with respect to image noise is highest in case of the low mono-energetic channels as well as the iodine density map itself 19,20 . This is also well reflected in the temporal contrast-to-noise ratio ( Fig. 1b) and relative signal enhancement (Fig. 1c), where the virtual MonoE-40 keV and iodine density images outperform the conventional HU image by far.
Since the temporal opacification of the myocardium is quantitatively and most sensitively depicted (as an iodine density offset arising from soft tissue is nearly absent) within the iodine density map, this channel was used for further iodine-dynamics analysis. Figure 2a displays a transversal porcine thorax slice for different scan times after the initial injection of contrast agent. Notice a subsequent bolus pass/enrichment of iodine within the right and left heart ventricle as well as myocardium for the exemplary scan times of 7.5 seconds, 10.9 seconds and 23.3 seconds, respectively. In a first step, imaging data of 15 scan points was temporally modelled in a pixel-wise manner using gamma variate fit functions, as proposed in literature 21,22 . Using a two-step fitting routine, only the first pass of contrast agent was considered so that artifacts arising from recirculation were excluded. Figure 2b displays data-points and the corresponding fit for an exemplary pixel in various regions of interest (ROI) showing a high degree of correlation of model and measurements (as quantified by the coefficient of determination R 2 being close to 1), even in case of a non-equidistant temporal sampling. The corresponding transversal thorax slices obtained from the model (Fig. 2c) are of smoother appearance than the original data, contingent on the fact that outliners and noise are damped via the fitting routine. This however is an essential premise for the meaningful retrieval of perfusion and peak enhancement map. Please note distinctive discrepancies between model and measured data for late scan times arising from a correction of recirculation, as indicated by low iodine density areas, among others, in the right ventricle.
In a next step the fitted data was used to derive the absolute, quantitative iodine peak enhancement, perfusion, time to peak and volume maps as depicted in Fig. 3, using the maximal slope method 23 . A peak enhancement value ( Fig. 3a), which is given as the maximal absolute increase in iodine density within the measured time frame, of 0.92 mg/ml (11,86% of the value found in the descending aorta) was determined within the myocardium. The corresponding perfusion value (Fig. 3b), i.e. the maximal slope of the temporal iodine density curve, was determined to 0.085 mg/ml/s. This value amounts to only 6,7% of the value found within the descending aorta, contingent on a limited influx of contrast agent through the coronary arteries in combination with a slow uptake/diffusion within the dense muscular tissue. Interestingly similar values were found in the ventral region of the lungs. As expected for a healthy animal, the time to peak map (Fig. 3c) indicates a simultaneous and consistent perfusion throughout the entire myocardium, approximately 4-5 seconds after the bolus passes through the descending aorta. Only a small ventral region with a delayed maximal opacification was found. However as depicted in the very right line plot in Fig. 2b the iodine density itself only yields minimal changes within the timeframe of 16-24 seconds. Finally, the overall volume of iodine perfusing the porcine thorax was calculated by integrating the temporal iodine density curves. As the myocardial enrichment remains high for late scan times and a flush out of iodine from the muscle is relatively slow in comparison to vascular and ventricular heart structures, the fitted curves were integrated within a timeframe of 50 seconds (so that the iodine density drops below a threshold value of 0.25 mg). Here a value of 29.89 mg/ml*s was calculated within the myocardium which accounts to 29,43% of the volume dwelling in the aorta. An overview of all perfusion parameters for various exemplary locations in the thorax is given in Table 1.
The simultaneous analysis of all perfusion parameters allows for a specific diagnosis of myocardial disorders: a reduced peak enhancement with a normal time-to-peak value indicates scarred/fibrotic myocardium, limiting the absolute capacity of the myocardium, while in presence of a delayed peak value an ischemic disorder is considered 24 . Further, the perfusion pattern, which is related to the temporal change of blood flow in both volume and velocity, is thereby a good indicator for a blocked or limited supply of blood through the coronary arteries. Therefore, quantitative CT perfusion of the myocardium assessed by absolute iodine concentration measurements over time has the potential to significantly increase the performance of CT for detection of myocardial www.nature.com/scientificreports www.nature.com/scientificreports/ ischemia in patients with ischemic heart disease. The blood volume (or in this case the iodine volume map), gives additional information on the overall blood intake capacity, its storing as well as the flush out behaviour.
A major drawback of a precise flow dynamics analysis is the necessity of a sufficient temporal sampling of the iodine density curve. This implies acquisitions before and during iodine contrast enrichment until myocardial saturation is reached, which poses concerns with respect to radiation dose. For an initial assessment of the myocardial blood flow (for instance in case of a screening modality or when myocardial ischemia is uncertain), we propose to derive the quantitative, iodine peak enhancement map from a semi-static two shot approach in stress and/or rest without the need of modelling/fitting the data. For this, two scans are obtained in the myocardial saturation phase (in this case at 23.3 and 24.9 seconds). To obtain a homogeneous pattern within the myocardium, the corresponding two iodine densities are averaged and smoothed by a minor Gaussian-blur (to mimic the fitting routine). Finally, to directly calculate the peak enhancement, an offset value from literature or a scout scan/control measurement is subtracted. Note that in a conventional HU-based assessment an arterial curve must be obtained in order to normalize data, which however is obsolete in case of quantitative, absolute iodine density values. Iodine density values vs. scan time as obtained from the dynamic dual-energy perfusion CT for an exemplary pixel for various region of interest as indicated by yellow boxes in (a), which were selected by an experienced radiologist. A high accordance of data and fit was obtained using a gamma variate fit model. Note that only the first pass of contrast agent was considered, and a recirculation of blood was correspondingly neglected for the fit. (c) Transversal quantitative iodine density slices of the porcine thorax obtained from the model being in good agreement with the measured data as shown in (a). Note that distinctive discrepancies for late scan times arise from the neglection of recirculation. Figure 4a,b compare the ground truth peak enhancement map of the fully sampled scan with the one obtained by the two-shot acquisition method, respectively. Obviously, the peak enhancement values within the ventricles and arteries is underestimated by far. However, in case of the saturated myocardium the two-shot approach yields very similar values to the model. The relative error within the peak enhancement (Fig. 4c) accounts to values in the range of 5-15% and is considerably homogeneous within the overall myocardium. Note that the remaining dynamics parameters as derived by the maximum slope method remain unassessed as here the precise onset as well as saturation of iodine throughout the myocardium must be tracked in time.
perfusion, (c) time to peak and (d) volume maps of the porcine thorax, derived from the gamma variate fit model using the slope method. The peak enhancement indicates the maximal increase in iodine density, the perfusion relates to the highest temporally gradient in influx of iodine, the time to peak indicates the point in time when maximal peak enhancement is reached, and the iodine volume is related to uptake/storing/flush out behaviour within the respective structures. Note that areas, where the gamma variate fit failed (for instance in the cava inferior) or areas where the iodine density remained below 0.55 mg/ml during the scan are color-coded black.

Discussion
Within this study, we showed that by using a state-of-the-art dual-energy perfusion CT imaging technique, it is possible to deduce absolute iodine-based flow dynamics of the heart, arteries and myocardium. This allows for an assessment of functional properties of the latter in a quantitative manner. Thereby, major advantages arise in comparison to conventional HU-based diagnosis: first of all, the increase in image contrast during perfusion is most pronounced in the iodine density map, helping with a better visual assessment of the myocardium and potential ischemic or infarcted regions. This allows for a more rigid fit of the imaging data as well as the potential to reduce the amount of contrast agent and radiation dose. Secondly, the image values examined by the radiologist are not affected by beam hardening artifacts and hence are independent from factors such as patient size and variations between different CT platforms, facilitating a value-based assessment. Finally, the obtained image data is quantitative, which renders an acquisition and consequent normalization with an arterial curve unnecessary, provided that the concentration of contrast agent with respect to patient characteristics is reliably controlled. Given the establishment of an extensive clinical database comparing iodine-based dynamical indicators for various myocardial dysfunctions, in both stress and rest, we believe that DECT may allow for a fast, rigid, highly sensitive and specific assessment of heart diseases.
However, an essential requirement for the establishment of such a database and using data without normalization is that the quantitative iodine densities present in the human myocardium can be measured precisely and in a reproducible manner. In the here presented animal study, a volume of 40 ml of Ultravist 300 was injected and an offset adjusted iodine density of 0.92 mg/ml found within the myocardium. Note that the corresponding iodine in-vivo concentration of 154 mg/kg-bodyweight, is well below volumes typically applied in clinical routine. In a very recent study comparing various dual-energy scanners Sellerer et al. demonstrated, by using an abdominal phantom and mimicking various patient sizes, that even low iodine concentrations below 1 mg/ml (dependent on underlying dual-energy technique and patient size) can be reliably measured with a relative error around 10% only 18 . Older studies reported on similar values, i.e. a constant absolute error of 0.1 mg/ml in case of iodine concentrations below 1 mg/ml 25 . Given a relative increase of iodine density within the healthy myocardium of 575% (0.16 mg/ml to 1.08 mg/ml) during contrast uptake within our study, the accuracy of modern systems can be considered as sufficiently high.
We further demonstrated that a semi-static two-shot approach yields a deviation of 5-15% in the iodine peak enhancement value in comparison to one derived from the full scan. Ensuing from previous studies showing that ischemic and infarcted regions of the myocardium exhibit a by 23% and 47% decreased mean iodine density, we believe that the quantitative semi-static approach may not be highly specific, however may be a good first indicator for myocardial malfunction (e.g. in case of a reduced peak enhancement or delayed time to peak value) in an initial assessment 16 . Further, previous studies on quantitative iodine perfusion implied the usage of 60 mL of Ultravist 370 (approximate iodine concentration of 317 mg/kg-bodyweight) resulting in an iodine density of 2.56 mg/ml within the healthy myocardium 16 . As we found a value of approximately 1 mg/ml in the myocardium (iodine concentration of 154 mg/kg-bodyweight) which can be reliably measured with modern DECT systems and allows for a meaningful modelling of the data, we believe that a reduction of applied iodine volume/concentration by a factor of 2 is possible.
Finally, in this animal study 15 scan points were used for the analysis, in order to model early enrichment of arteries and ventricles, accounting to a computed tomography (CT) dose index (CTDIvol) of 15 × 9 mGy equalling 135 mGy. In a clinical setting when only the myocardium is of interest early scan times after the injection can be neglected and an equidistant sampling of myocardial perfusion slope and saturation with 8-12 points is generally sufficient. Correspondingly a CTDIvol of 72-108 mGy can be expected, being well below a critically considered value of 250 mGy.
The purpose of this experiment is considered a proof of principle investigation. As the DECT technique offers several advantages regarding quantitative image results using iodine contrast agent, the main point of our manuscript focuses on a initial investigation of quantitative perfusion measurements. The added value of this study is to show the advantage of quantitative iodine density maps in comparison to conventional CT images. In terms of the improved contrast-to-noise ratio and the relative signal enhancement using contrast agent density maps in DECT perfusion imaging, our results are relevant to several clinical perfusion use cases with animals or even humans. Indeed, the improvement of clinical diagnosis should be considered in further statistically meaningful studies with the benefit of various patient groups and myocardial disorders in rest and stress condition. This will permit the establishment of a valid database of quantitative iodine-based hemodynamical parameters. Hereby it will be of essential importance to clarify which of the quantitative hemodynamical indicators are of major diagnostic importance. Further, within this study a dual-layer IQon Scanner from Philips was used, providing a detector coverage of 4 cm only. A spatially limited examination of the heart (especially the myocardium) however bears the risk of overlooking locally confined flow defects. Under these circumstances, vendors are obligated to increase detector coverage in order to screen the entire myocardium at once. were acquired with a 64-slice single source dual-layer spectral CT IQon scanner with a detector coverage of 4 cm and a rotation time of 0.27 seconds (Philips Healthcare, The Netherlands). The spatial resolution is 2.47 pixels per mm resulting in a pixel size of 0.40 × 0.40 mm 2 . The ECG triggered temporal resolution aggregates to 1.08 scans per second or a full acquisition of the specified range every 0,92 second. Firmware version on the scanner was 4.1.0.0. The perfusion scan was conducted in stationary mode with 120 kVp and a mean tube current of 100 mAs. Scans were reconstructed with a Philips B-kernel and a slice thickness of 3 mm. Within 26 seconds 36 scans were recorded, so that the full first pass of contrast bolus was imaged. Hereby the scans were electrocardiographically (ECG)-triggered and adjusted to the cardiac cycle with an increment of zero. Each scan was conducted with a CTDIvol of 9 mGy, resulting in a CTDIvol of 135 mGy and 18 mGy for the extensive analysis of myocardial dynamics (15 scan points) and the two-shot approach, respectively. The perfusion scan was conducted using 40 ml of contrast agent (Ultravist 300, Bayer, Bayer AG, Leverkusen, Germany, iodine content 300 mg/ml) at a flow rate of 4 ml/s. The contrast agent was injected into the ear vein via an 18-gauge catheter using a dual syringe injection system (Stellant, MEDRAD, Inc., Indianola, PA, USA).

Materials and Methods
Contrast-to-Noise and relative signal enhancement. Conventional, virtual Mono-keV and Iodine density reconstructions alongside with their standard deviations (within respective ROIs of 50 mm 2 ) were directly obtained from the Philips IntelliSpace Portal. The relative increase in signal (RIS, Fig. 1c) and temporal contrast-to-noise ratio (CNR, Fig. 1b Retrieval of myocardial dynamics. A certain slice within the thorax was determined for further processing and 15 data points in time, with minor structural movement, were selected from the iodine density images. In order to model the first pass of contrast bolus through the heart, arteries and myocardium, gamma variate functions were used, as proposed in literature for both magnetic resonance imaging (MRI) and CT 21,26 : t t 0 ( ) 0 whereas I(t) is the time dependent iodine density and A, t 0 , α, β are fit parameters. Thereby the data was fitted in a two-step sequence: in a first fit all 15 data points were used and a preliminary peak time determined. In a second fit data exceeding the initial peak by more than 12 seconds, i.e. accounting to the second pass of contrast agent, was correspondingly neglected. Afterwards quantitative iodine perfusion, peak enhancement, time to peak and iodine volume (in imitation of the conventional blood volume) were calculated based on the maximal slope method 23 : Finally, aforementioned maps were filtered with a median-filter with a 3 × 3-kernel to reduce noise arising from the pixel-wise assessment. Note that data was not normalized or scaled with respect to arterial input functions, hemocrit-scale, cerebral blood volume (CBV) factor or other physiological models for multiple reasons: first of all, within this animal study corresponding values and models are not available. Secondly, we wanted to provide absolute and quantitative, i.e. unscaled, values to the readers, which are independent from assumptions made in underlying physiological models. Correspondingly the iodine volume map presented here does not represent iodine volumes in the sense of millilitres, but rather gives quantitative measure on the area under the curve of the temporal iodine density in the unit of mg/ml*s. The latter is thereby related to the tissues overall blood intake capacity, its storing as well as the flush out behaviour.

Data availability
The data generated in the current study is available from the corresponding author on reasonable request.