|
|
|||||||||
Clinical Investigations |
1 Department of Medicine, Feinberg School of Medicine, Northwestern University, Chicago, Illinois
2 Division of Cardiology, Department of Medicine, The Weatherhead PET Center for Preventing and Reversing Atherosclerosis, University of Texas Medical School at Houston and the Memorial Hermann Hospital, Houston, Texas
| ABSTRACT |
|---|
|
|
|---|
2 analysis (P < 0.001). The relative odds ratios of having stress-induced myocardial perfusion abnormalities for a resting homogeneity index outside 1 SD of healthy reference subjects was highly predictive and substantially greater than for standard risk factors. Conclusion: Patchy heterogeneous resting myocardial perfusion by noninvasive cardiac PET quantified objectively using Markovian homogeneity analysis, and its improvement after dipyridamole, are powerful independent predictors of even mild stress-induced perfusion abnormalities, more than standard risk factors, consistent with coronary microvascular dysfunction as an early marker of preclinical CAD for potential preventive treatment.
Key Words: myocardial perfusion PET coronary artery disease heterogeneity homogeneity microvascular function risk factors
| INTRODUCTION |
|---|
|
|
|---|
The hallmark of coronary endothelial dysfunction is mild heterogeneous vasoconstriction of coronary arteries or coronary microvasculature under a wide spectrum of different conditions and vasomotor stimuli involving many different mechanisms, including inhibition of vasodilator mechanisms or activation of vasoconstrictor mechanisms by many different interacting vasoactive mediators. Heterogeneity of coronary endothelial function has been well documented in humans (1012) with associated altered coronary blood flow or perfusion reflecting coronary arteriolar as well as epicardial arterial endothelial dysfunction (1317). Resting coronary flow falls by approximately 20% after inhibition of coronary endothelial nitric oxide production without significant reduction in maximum coronary flow or coronary flow reserve (1822), thereby reflecting altered resting microvascular function.
Therefore, we have hypothesized that the visually apparent heterogeneity of resting myocardial perfusion or its improvement after dipyridamole stress on high-quality, noninvasive PET images outside the limits of healthy control subjects is one manifestation of coronary microvascular dysfunction associated with endothelial dysfunction (23). Testing this hypothesis is important because coronary endothelial dysfunction is associated with early preclinical coronary atherosclerosis, increased coronary events (3,4), and subsequent clinically manifest CAD many years later (6), thereby providing a basis for intense, lifelong, pharmacologic preventive treatment.
Coronary flow reserve and myocardial perfusion imaging after pharmacologic arteriolar vasodilation for identifying flow-limiting coronary artery stenosis as first reported from this laboratory (2427) is now widespread as a routine clinical diagnostic procedure. In this paradigm, the resting perfusion image serves as a baseline for comparison with the stress perfusion image for identifying discrete regional perfusion abnormalities due to flow-limiting coronary artery stenosis, myocardial scar, or hibernating myocardium. However, the current report analyzes and quantifies the distinctly different diffuse patchy heterogeneity of resting myocardial perfusion as a marker of coronary endothelial dysfunction associated with coronary atherosclerosis, independently from and around these traditional discrete regional myocardial perfusion defects caused by flow-limiting stenosis or myocardial scar.
In this study, we use a mathematic technique from Markovian homogeneity analysis (28) to provide precise, objective, automated quantification of resting perfusion heterogeneity in 1,034 subjects, its normal limits in 50 healthy reference subjects, and its close association with documented CAD, thereby demonstrating a basic new observation in myocardial perfusion imaging with important clinical implications.
| MATERIALS AND METHODS |
|---|
|
|
|---|
PET
Patients were instructed to fast for 4 h and abstain from caffeine, theophylline, and cigarettes for 24 h before study. As previously described (2933), PET was performed using the University of Texas designed, Positron Posicam Auricle, bismuth germanate, 2-dimensional (2D) multislice tomograph with a reconstructed resolution of 10-mm full width at half maximum. Using a rotating rod source containing 148185 MBq (45 mCi) of 68Ge, transmission images to correct for photon attenuation contained approximately 4060 million counts. Emission images obtained after intravenous injection of 9251,850 MBq (2550 mCi) of generator-produced 82Rb contained 2050 million counts depending on the age of the generator and size of the patient. After resting 82Rb data acquisition, dipyridamole (0.142 mg/kg/min) was infused for 4 min. At 4 min after completion of the dipyridamole infusion, the same dose of 82Rb was given intravenously.
Automated Quantitative Analysis of PET Images
Completely automated analysis of the severity and size of PET abnormalities was performed by previously described software (2933). A 3-dimensional (3D) restructuring algorithm generates true short- and long-axis views from PET transaxial cardiac images acquired in 2D tomographic mode to minimize scatter. From circumferential profiles, 3D topographic views of the left ventricle are reconstructed showing relative regional activity distribution divided into lateral, inferior, septal, and anterior quadrant views of the 3D topographic display corresponding to the coronary arteries illustrated in Figure 1.
|
Activity is normalized to the maximum 2% of pixels in the whole heart dataset. Regions of each quadrant are identified as outside 97.5% confidence intervals (CI) or 2.5 SD of 50 healthy control subjects with no risk factors by complete medical history. The percentage of circumferential profile units outside 97.5% CI is calculated automatically after correcting for any misregistration of attenuation and emission images that commonly cause artifactual defects (34).
Markovian Homogeneity Analysis.
Markovian texture or homogeneity analysis characterizes an image by examining the probability that a pixel with a given intensity will have a neighbor with a different intensity (28), where Pd(m) is the probability that 2 adjacent pixels have intensity values that differ by m. For this study, the homogeneity index, H, is given by the following equation:
![]() | (Eq. 1) |
The homogeneity index thus quantifies mathematically the intuitive notion of homogeneity. A perfusion image that is inhomogeneous or diffusely patchy has a small homogeneity index near 0, whereas a uniform image has a large index near 1. Conversely, an image with a small homogeneity index can be considered to be heterogeneous and vice versa. Each pixel can have a maximum of 8 neighbors: above, below, left, right, above left, above right, below left, and below right, where a complete topographic map is like the surface of a cylinder and "wraps" along the radial dimension. There are 5,184 unique pixel pairs for a 21 x 64 matrix, assuming the 64-axis wraps. Intensity of the image matrix is normalized to 1,000. More generally, for an N-by-M matrix, assuming the M-axis wraps, there are (4NM 3M) unique pixel pairs. Radial pixel size and, hence, heart size does not impact the calculation of the heterogeneity index.
Equation 1 shows that as the differences among neighboring pixel units become larger for severe defects, the coefficient 1/(1 + m)2 decreases rapidly with increasing m, corresponding to increasing differences of intensity between neighboring pixels. Expanding the summation in Equation 1 for differences of m = 0, 1, 2, 3, etc., yields the following:
![]() |
Application of Homogeneity Analysis to PET Perfusion Images.
The rest and stress scans are displayed as topographic displays in 4 quadrant views (lateral, inferior, septal, and anterior).
For applying Equation 1 to this topographic map, 3 modifications were made as follows: (i) The basal 4 slices are discarded to avoid count variability in the membranous septum and the apical 2 slices are discarded to minimize partial-volume effects and variability in locating the last apical slice; (ii) pixels with intensity values below 500 are reset to 500, and pixels with intensity values above 850 are reset to 850 to eliminate any effect on the homogeneity index, H, of very low activity levels of myocardium scar and to eliminate effects of the highest activity levels of normal myocardium.
In effect, these limits further confine the homogeneity analysis to relative activity values ranging from 50% to 85% of maximum on each PET image, thereby excluding extreme values as from severe defects or hot spots that would bias the value of H for quantifying more subtle differences among pixel units; (iii) these modified intensity values, 500850 inclusive, are scaled into an integer range of 35 levels, so that each new intensity level represents 1% of the range, thereby mathematically further restricting the analysis to myocardium that is not scarred, severely ischemic, or maximally perfused. Consequently, the homogeneity index is not greatly influenced by severe perfusion defects. The degree of small-scale diffuse heterogeneous "patchiness" is objectively quantified on resting and stress images and the rest-to-stress change independently of, separately from, or around severe discrete regional perfusion defects due to myocardial scar, reduced coronary flow reserve of flow-limiting stenosis, or the base-to-apex longitudinal perfusion gradient due to diffuse disease (30). The scaling impacts the heterogeneity index and determines the "coarseness" or the degree of patchiness being quantified.
Mean and SD values for the homogeneity index, H, were computed for the 50 healthy control subjects just as for the other automated quantitative measurements on the PET images for comparison with the patients in this study.
Statistical Analysis
All statistical analyses were performed using SPSS version 11.5 (SPSS Inc.). Data are reported as mean ± 1 SD or SEM as appropriate.
Analysis 1.
Multivariate logistic regression analysis was performed with the independent variables being the continuous values of the resting homogeneity index (rH), the rest-to-stress change in H (rsH
), and all discrete risk factors of age, sex, history of diabetes, hypertension, high cholesterol, family history of vascular disease, excess weight, smoking, menopausal status, and lack of exercise. The dependent variables were an abnormal PET after dipyridamole (stress PET) defined as either the lowest mean quadrant activity on the stress PET image, Q, outside 1 SD of healthy reference subjects (Q < 1 SD), indicating flow-limiting stenosis, or the base-to-apex longitudinal perfusion gradient (L) after dipyridamole outside 1 SD of healthy reference subjects (L < 1 SD), indicating diffuse coronary artery narrowing as previously demonstrated (30). An abnormal stress PET therefore includes all cases with any abnormality of either Q or L and excludes those with completely normal Q and Lthat is, both Q > 1 SD and L > 1 SD. It indicates a not-normal PET perfusion scan after dipyridamole stress attributed to either a localized regional defect or an abnormal base-to-apex longitudinal perfusion abnormality outside 1 SD of healthy reference subjects, thereby objectively documenting even mild CAD.
Analysis 2.
Multivariate linear regression analysis was performed with the same independent variables as abovethat is, the continuous values of the resting homogeneity index (rH), the rest-to-stress change in H (rsH
), and all discrete risk factors. The dependent variable is the continuous value of the lowest mean quadrant activity of the stress PET image (Q), indicating the quantitative severity of regional perfusion defects after dipyridamole caused by flow-limiting stenosis.
A Pearson
2 analysis was performed for the discrete variables as follows: abnormal homogeneity (rH < 2 SD or rsH
< 2 SD), borderline homogeneity (rH and rsH
within 12 SD), normal homogeneity (rH > 1 SD and rsH
> 1 SD), abnormal stress scans (Q < 2 SD or L < 2 SD), borderline stress scans (Q and L within 12 SD), and normal stress scans (Q > 1 SD and L > 1 SD). A 2-tailed P value < 0.05 was considered statistically significant.
| RESULTS |
|---|
|
|
|---|
|
|
Figure 4 illustrates the base-to-apex longitudinal perfusion gradient of these same 3 pairs of reststress PET scans, being within normal limits for the first pair (Fig. 4A) and markedly abnormal for the other 2 examples (Figs. 4B and 4C), indicating diffuse CAD plus a severe localized flow-limiting stenosis (Fig. 4B) or without a severe localized flow-limiting stenosis (Fig. 4C).
|
), and the risk factors as the independent variables. The discrete dependent variable is any abnormality of the stress perfusion scan, either the minimum quadrant average activity outside, or greater than, Q > 1 SD of healthy reference subjects or the base-to-apex longitudinal perfusion gradient outside, or greater than, L > 1 SD of healthy reference subjects on stress PET images. As expected, standard risk factors are predictive of abnormal stress perfusion images. A family history of vascular disease and smoking were not significantly predictive because of the brevity of details of the history recorded in the database options. The family history did not differentiate among parents, siblings, or remote relations. Smoking did not differentiate among remote brief smoking, active current smoking, or amount of smoking.
|
The negative values of B for the homogeneity index (rH) and its rest-to-stress change (rsH
) indicate that high values of rH or rsH
are associated with very high probability of normal stress perfusion images (the minimum quadrant average activity Q less than or within 1 SD of healthy reference subjects) and a low probability of abnormal stress perfusion images (Q greater than or outside 1 SD of healthy reference subjects). Similarly, low values of rH or its rest-to-stress change, rsH
, are associated with a very high probability of abnormal stress perfusion defects (Q > 1 SD) and a low probability of normal stress perfusion images (Q < 1 SD).
For logistic regression, which uses a sigmoid/logistic model instead of a linear one, B is the logarithm of the odds ratio, and EXP(B) is the ratio of the oddsthat is, the probability of something occurring divided by the probability of something not occurring, quantified as EXP(B). For example, in this analysis for diabetes, B is 1.64, the EXP(B) is e1.64 or 5.15,that is, the odds ratio for diabetes. Therefore, a person with a history of diabetes carries a 5.15 times greater odds of an abnormal stress PET scan than the odds for a person without a history of diabetes.
For the homogeneity index, B is 23.05, the EXP(B) is e23.05 or 1.03 x 1010that is, the odds ratio for the homogeneity index. Therefore, the odds of a patient with a normal homogeneity index (HI = 1.0) having an abnormal stress PET scan is only 0.000000000103 of the odds for a patient with an abnormal homogeneity index (HI = 0.0+) of having an abnormal stress PET. Similarly, the odds for a patient with an abnormal homogeneity index of having a normal stress image is an equally small percent of the odds for a patient with a normal homogeneity index of having a normal stress image. From the other viewpoint, the odds for a person with an abnormal resting homogeneity index (HI = 0.0+) of having an abnormal stress PET is 1/0.0000000103 or 9.7 x 109 times the odds of a person with a normal homogeneity index (HI + 1.0) of having an abnormal stress PET. Separately and independently of the resting homogeneity index, rH, the rest-to-stress improvement in the homogeneity index, rsH
, after dipyridamole stress is comparably predictive of CAD with a comparable odds ratio.
Table 2 summarizes the multivariate linear regression analysis with the independent variables being the resting perfusion homogeneity index (rH), the reststress change in homogeneity index (rsH
), and all risk factors together in the first 5 rows and for the risk factor alone without the PET data in rows 7 through 17. The single dependent variable is the continuous quantitative severity of stress perfusion defectsthat is, the minimum average quadrant activity on the stress PET images (Q). This analysis shows that the resting homogeneity index and its rest-to-stress improvement are closely correlated with stress-induced regional perfusion defects separately from and independently of other risk factors (P < 0.001). For linear regression analysis of continuous variables, B is the "slope" of the regression equation (coefficient of the independent variable in the fitted equation) that is substantially greater for the resting homogeneity index and its rest-to-stress change than for any of the standard risk factors.
|
Table 3 shows the
2 analysis with numbers of subjects in each category where homogeneity was defined as "abnormal" if either the resting homogeneity index or its rest-to-stress improvement were outside or < 2 SD of healthy reference subjects, "normal" if both were >2 SD, and "borderline" for all other combinations of the rest and rest-to-stress change as mixed 12 SD, < 2 SD, and > 2 SD. Similarly, stress images were defined as "abnormal" if either the lowest mean quadrant average activity (Q) caused by flow-limiting stenosis or the longitudinal base-to-apex perfusion gradient (L) due to diffuse CAD were < 2 SD of healthy reference subjects, normal if both were > 2 SD, and borderline for all other combinations. In Table 3, the distribution of subjects in the binary discrete categories of homogeneity and stress perfusion categories is significant with P < 0.001. Table 4 shows the relative or percentage distribution of the
2 distribution of raw numbers in Table 3 expressed in Table 4 as the percentage of the patients in each homogeneity category in rows from left to right across the table. Figure 5 is a bar graph of the relative percentage distribution of patients in each of the homogeneity categories derived from the
2 analysis in Table 3 and the percentage distributions in Table 4.
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Endothelial dysfunction as a cause of coronary arterial and microvascular vasoconstriction is well documented in experimental studies and in humans by coronary arteriography or Doppler flow-velocity wires or catheters. Vascular mediators derived from coronary endothelium include prostacyclin, nitric oxide, thromboxane, endothelin, bradykinin, angiotensin, serotonin, substance P, C-type naturetic peptide ([CNP], an endothelium-derived hyperpolarizing factor), and others. The mechanisms may be inhibition of normal vasodilatory mediators such as nitric oxide or activation of vasoconstrictor mediators such as endothelin.
The stimuli, mediators, and the vascular responses of epicardial coronary arteries and the coronary microvasculature are quite different, even divergent. For example, in the epicardial coronary arteries, acetylcholine-induced vasodilation is mediated by nitric oxide (20,35,36). However, in the coronary microcirculation, acetylcholine-induced arteriolar vasodilation and increased coronary flow are not mediated by nitric oxide (18,37,38). With epicardial artery endothelial dysfunction, acetylcholine causes arterial vasoconstriction while arteriolar vasodilation with increased flow remains intact as an example of divergent pathophysiologic behavior of the macrovasculature and microvasculature of the heart (1822).
As a further example, endothelial nitric oxide production mediates epicardial coronary artery vasodilation during exercise (19) but is not involved in arteriolar vasodilation and increased coronary flow during exercise (19) unless there is a flow-limiting stenosis in which nitric oxide helps maintain perfusion during exercise (39). In opposition to these vasodilator mechanisms, endothelin is a powerful coronary arteriolar vasoconstrictor that is activated in coronary atherosclerosis in parallel with inhibition of nitric oxide production.
Thus, there is no single specific vasomotor abnormality, gold standard, diagnostic test, or even definition that identifies or defines coronary endothelial dysfunction. Our data indicate that resting myocardial perfusion heterogeneity is one manifestation of this wide spectrum of coronary vascular behavior that is a powerful independent predictor of preclinical CAD, more than standard risk factors. In view of the different, sometimes divergent, arterial and arteriolar behaviors in response to the wide variety of vasoactive mediators, resting perfusion heterogeneity would not necessarily be expected to parallel the effects of intracoronary acetylcholine or cold pressor testing, just as the arteriolar response to acetylcholine with increased blood flow does not parallel its vasoconstrictive effect on epicardial coronary arteries in CAD.
The limitations of this study deserve comment. A criticism may be that the limited resolution of PET cannot resolve the small regions of heterogeneous perfusion previously described in experimental animals (40) or the subendocardial underperfusion that is an effect of flow-limiting stenosis. The heterogeneity that is visually apparent and objectively quantified in this study involves regions of myocardium greater than the 1-cm3 scanner resolution, consistent with the arterial distribution of coronary arteries and their secondary or tertiary branches demonstrated to have heterogeneous endothelial function by coronary arteriography and intracoronary Doppler flow-velocity measurements. Therefore, the heterogeneity that we observe by PET perfusion imaging is separate and unrelated to the dispersion of perfusion in small 1-mm myocardial samples for microsphere measurements of perfusion reported for experimental animals (40).
The heterogeneous resting perfusion in this study was quantified separately from, independently of, and around significant regional perfusion defects caused by flow-limiting stenosis and, therefore, does not involve subendocardial hypoperfusion due to reduced perfusion pressure or reduced early diastolic subendocardial filling caused by flow-limiting stenosis.
The limits of heterogeneity were determined from 50 healthy control subjects imaged on the same scanner and software as the patients so that the technical limitations of PET or any potential effects of microscopic dispersion apply equally to both sets of subjects with the significant differences reported here. The application of homogeneity analysis to PET perfusion images requires careful attention to the technical details of cardiac PET with 82Rb that are different than those required for cancer PET, including the necessity of lower spatial resolution in favor of a high count density, 2D imaging to reduce scattered radiation, high-count, low-noise, filtered backprojection reconstruction, and compulsive correction of emissiontransmission image coregistration (34).
Coronary arteriography was not performed or used as a comparative gold standard in all of these patients. The percentage diameter stenosis as a measure of the severity of CAD is notoriously inadequate because of diffuse disease. The base-to-apex longitudinal perfusion gradient by PET perfusion imaging identifies early diffuse CAD better than regional stress-induced perfusion defects of flow-limiting stenosis as we have demonstrated (30). Because coronary atherosclerosis is a continuous spectrum from early mild stages to severe stenosis, the conventional categorization of arteriograms or perfusion images into "normal" or "abnormal" for determination of sensitivity or specificity is artificial and incorrect, particularly when defined as outside 2 SD of normal. Our multivariate regression analysis using continuous quantitative variables confirms the continuous spectrum of these endpoints. Accordingly, for added certainty, we performed logistic regression analysis using as thresholds of our endpoints a cutoff of < 1 SD as "not normal"that is, a greater probability of being "abnormal" than "normal" to include the great extent of mild preclinical CAD with potential for plaque rupture and coronary events.
Although resting perfusion heterogeneity or its improvement is associated with "not normal" stress perfusion PET scans, some patients with resting perfusion heterogeneity had normal stress perfusion PET scans. By association with otherwise comparable patients with stress-induced perfusion changes, our findings suggest that such patients with heterogeneity or its improvement without stress-induced defects are at risk for vascular disease. However, our data do not prove that point. Proof would require many years follow-up of such patients off lipid medications to determine the outcomes, a difficult study that is not ethically appropriate.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: K. Lance Gould, MD, The Weatherhead PET Center, University of Texas Medical School, Room 4.256MSB, 6431 Fannin St., Houston, TX 77030.
E-mail: gould{at}pet.med.uth.tmc.edu
| REFERENCES |
|---|
|
|
|---|
Related articles in JNM:
This article has been cited by other articles:
![]() |
R. Djaberi, J. o. t Roodt, J. D. Schuijf, T. J. Rabelink, E. J. de Koning, A. M. Pereira, M. P. Stokkel, J. W. Smit, J. J. Bax, and J. W. Jukema Endothelial Dysfunction in Diabetic Patients with Abnormal Myocardial Perfusion in the Absence of Epicardial Obstructive Coronary Artery Disease J. Nucl. Med., December 1, 2009; 50(12): 1980 - 1986. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Gupta, J. Tongers, and D. W. Losordo Human Studies of Angiogenic Gene Therapy Circ. Res., October 9, 2009; 105(8): 724 - 736. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Gould Does Coronary Flow Trump Coronary Anatomy? J. Am. Coll. Cardiol. Img., August 1, 2009; 2(8): 1009 - 1023. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Fricke, S. Eckert, A. Dongas, H. Fricke, R. Preuss, O. Lindner, D. Horstkotte, and W. Burchert Myocardial Sympathetic Innervation in Patients with Symptomatic Coronary Artery Disease: Follow-up After 1 Year with Neurostimulation J. Nucl. Med., September 1, 2008; 49(9): 1458 - 1464. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Gould, T. Pan, C. Loghin, N. P. Johnson, and S. Sdringola Reducing Radiation Dose in Rest-Stress Cardiac PET/CT by Single Poststress Cine CT for Attenuation Correction: Quantitative Validation J. Nucl. Med., May 1, 2008; 49(5): 738 - 745. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Gould, T. Pan, C. Loghin, N. P. Johnson, A. Guha, and S. Sdringola Frequent Diagnostic Errors in Cardiac PET/CT Due to Misregistration of CT Attenuation and Emission PET Images: A Definitive Analysis of Causes, Consequences, and Corrections J. Nucl. Med., July 1, 2007; 48(7): 1112 - 1121. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Loghin, S. Sdringola, and K. L. Gould Does coronary vasodilation after adenosine override endothelin-1-induced coronary vasoconstriction? Am J Physiol Heart Circ Physiol, January 1, 2007; 292(1): H496 - H502. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. J. Gibbons, P. A. Araoz, and E. E. Williamson The Year in Cardiac Imaging J. Am. Coll. Cardiol., December 5, 2006; 48(11): 2324 - 2339. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | RSS | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |