| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2007;115:459-467.)
© 2007 American Heart Association, Inc.
Epidemiology |
From the Department of Neurology (M.W.L., M.S.), Johann Wolfgang Goethe University, Frankfurt am Main, Germany; the St. Georges University of London (H.S.M.), London, United Kingdom; the Julius Center for Health Sciences and Primary Care (M.L.B.), University Medical Center Utrecht, Utrecht, the Netherlands; the Department of Epidemiology and Biostatistics (M.L.B.), Erasmus Medical Center, Rotterdam, the Netherlands; and the Department of Community Medicine (M.R.), Lund University, Malmo University Hospital, Malmo, Sweden.
Correspondence to Matthias W. Lorenz MD, Johann Wolfgang Goethe-University, Department for Neurology, Schleusenweg 2-16, D-60528 Frankfurt/Main, Germany. E-mail matthias.lorenz{at}em.uni-frankfurt.de
Received March 23, 2006; accepted November 29, 2006.
| Abstract |
|---|
|
|
|---|
Methods and Results— Using a prespecified search strategy, we identified 8 relevant studies and compared study design, measurement protocols, and reported data. We identified sources of heterogeneity between studies. The assumption of a linear relationship between IMT and risk was challenged by use of a graphical technique. To obtain a pooled estimate of the relative risk per IMT difference, we performed a meta-analysis based on random effects models. The age- and sex-adjusted overall estimates of the relative risk of myocardial infarction were 1.26 (95% CI, 1.21 to 1.30) per 1–standard deviation common carotid artery IMT difference and 1.15 (95% CI, 1.12 to 1.17) per 0.10-mm common carotid artery IMT difference. The age- and sex-adjusted relative risks of stroke were 1.32 (95% CI, 1.27 to 1.38) per 1–standard deviation common carotid artery IMT difference and 1.18 (95% CI, 1.16 to 1.21) per 0.10-mm common carotid artery IMT difference. Major sources of heterogeneity were age distribution, carotid segment definition, and IMT measurement protocol. The relationship between IMT and risk was nonlinear, but the linear models fitted relatively well for moderate to high IMT values.
Conclusions— Carotid IMT is a strong predictor of future vascular events. The relative risk per IMT difference is slightly higher for the end point stroke than for myocardial infarction. In future IMT studies, ultrasound protocols should be aligned with published studies. Data for younger individuals are limited and more studies are required.
Key Words: atherosclerosis carotid arteries meta-analysis myocardial infarction stroke
| Introduction |
|---|
|
|
|---|
Clinical Perspective p 467
Different longitudinal studies have used different measurement methods and studied different populations. Furthermore, some studies have investigated relatively small populations, which resulted in large confidence intervals (CIs) for the risk estimates. We carried out a systematic review of the literature to provide an overview of the relevant studies, critically appraise the methods used, and where possible perform a meta-analysis to gain more robust estimates of the predictive value of increased IMT for clinical cardiovascular end points.
| Methods |
|---|
|
|
|---|
If the publications did not contain the full information necessary for meta-analysis, we obtained the missing information directly from the authors (see Acknowledgments). The National Heart, Lung, and Blood Institute in the United States provided the original data (limited-access dataset) of the Atherosclerosis Risk in Communities (ARIC) study, from which the hazard ratios (HRs) were newly calculated for uniformity of the estimates. Exclusion criteria and end point definitions were applied exactly as stated in the publications, but because the available dataset contained a later follow-up, the number of events is higher than in the publications (Table 1).
|
Statistical Analyses
The HRs and CIs of comparable studies were illustrated with forest plots. Pooled estimates were calculated with a random effects model as recommended by Blettner et al,13 with the method of moments approach from DerSimonian and Laird.14 Different standard deviations (SDs) of IMT distributions can be caused by population properties (eg, different age distributions12) or properties of the ultrasound technique (eg, reproducibility), or a combination of both aspects. If one assumes that the populations in the different studies share the same underlying IMT distribution and the same hazard function, then different SDs of IMT measurements are likely to be the result of measurement error, and the risk per IMT difference of 1 SD can be compared without conversion. On the other hand, if we compare populations with different baseline hazards that are influenced by a common IMT hazard function, we have to compare HRs per absolute IMT difference. We decided to implement both strategies: to compare HRs per 1 SD and those expressed as HRs per (absolute) 0.10-mm difference.
To ensure uniformity of estimates, we compared HRs that were adjusted by age and sex. Because the publications from the ARIC study4,6 were stratified by sex, we recalculated age- and sex-adjusted HRs from the original dataset (see Study Selection and Data Collection). To compare fully adjusted models, we tried to gain uniformly adjusted HRs by recalculation of our own datasets (M.R. for the Malmo Diet and Cancer Study [MDCS], M.B. for the Rotterdam Study, and M.W.L. for the Carotid Atherosclerosis Progression Study [CAPS]) and the ARIC data, and by personal communication with the other authors.
From every study, the HRs per IMT difference were log-transformed to obtain symmetric CIs. The standard error of this point estimator was recalculated for every study by log-transformation of the 95% CI limits. The difference between the log(lower limit) and log(upper limit) was divided by 2 x 1.96 (the 97.5% normal percentile) to obtain the standard error. The point estimator from every study was weighted with the inverse of the square of its standard error to calculate compound estimators. To investigate heterogeneity, we used the I2 statistics provided by Higgins et al.15 To identify study variables that contributed particularly to heterogeneity, we performed a meta-regression analysis with restricted maximum likelihood estimates that accounted for residual heterogeneity, as recommended by Thompson and Sharp.16 To assess possible publication bias, we used a nonparametric iterative trim-and-fill procedure suggested by Duval and Tweedy.17,18
The reporting of HRs per IMT difference is based on the assumption that there is a linear relationship between IMT and hazard (in the strict sense, a relationship implied by a linear model is log-linear but is mostly referred to as linear4,6). To reassess this assumption with reference to the publications of interest, we used a graphical technique introduced by Chambless et al.4,6 The relative hazards for vascular events are plotted against baseline IMT for a set of 3 to 6 IMT categories and compared with the assumed linear relationship. Because the IMT cut points must define large proportions of the samples, they depend on the distribution and cannot use the same reference category. The relative hazards for the IMT categories were recalculated or taken from the publications, and the mean IMT values for every IMT category were acquired from the corresponding author directly.
The recalculation of the ARIC, MDCS, and CAPS HRs and the meta-regression were performed with SAS (The SAS Institute, Cary, NC). Forest plots were produced by a Microsoft Excel macro program provided by Peter Rothwell and Enrico Flossmann and modified by M.W.L. Intraclass correlation coefficients were recalculated with SPSS 11.5 (SPSS Inc., Chicago, Ill). The trim-and-fill procedure was calculated manually with assistance of an Excel spreadsheet. M.W.L. had full access to and takes responsibility for the integrity of the data. All authors have read and agree to the manuscript as written.
| Results |
|---|
|
|
|---|
Table 1 provides an overview of the studies in terms of the study design, definition of clinical end points, observed population, sample size, follow-up period, carotid segments studied, definition of IMT, model construction (which accounts for potential bias factors), and report of results. The studies exhibited considerable variation in all of these aspects. In particular, there are marked differences in age distribution, with a lower limit that ranges from 1912 to 75 years.19 Furthermore, the publications can be divided into 2 categories: those that use mean carotid IMT, determined from a number of IMT measurements at specific positions or by automated software over a segment of the artery, and those that use maximal carotid IMT. The Rotterdam Study is the only one in which IMT was determined by both methods in the same population.7,8
For further analyses, the Kuopio Ischemic Heart Disease Risk Factor Study,2 the study reported by Kitamura et al,9 and the Longitudinal Investigation for the Longevity and Aging in Hokkaido County19 were excluded because they did not provide estimates that were comparable to the other studies.
Assessment of Potential Sources of Heterogeneity and Confounding
To determine the comparability of the populations, we calculated event rates (events per 1000 person-years) and absolute common carotid artery (CCA) IMT values for all the studies (Table 2). The event rates show distinct differences that may be largely attributable to the definition of end points. Maximal IMT was more or less comparable in terms of mean values, but the SDs were different. Mean IMT values were similar, but values from different studies exhibited different distributions.
|
The study design was longitudinal in all publications apart from those by Bots et al and Iglesias del Sol et al,3,7 who used a nested case-control design and a case-cohort design, respectively. The work by Iglesias del Sol et al7 was analyzed with the longitudinal studies, as the Cox regression model estimators of a case-cohort design are comparable with those of a longitudinal design.23 The publication by Bots et al3 was omitted because it refers to the same population as the other 2 publications on the Rotterdam Study and has a study design that provides results that are somewhat difficult to compare with those from longitudinal studies.
Figure 1 illustrates the definitions of the carotid segments for each study. This is obviously a major source of heterogeneity, because, between some studies, the "same" segments do not even overlap. Furthermore, the publications differed in terms of how the results were reported (Table 1). The majority of studies reported calculations from age- and sex-adjusted models. The most commonly reported estimator was the relative risk (provided as HR) per (linear) IMT difference of 1 SD. Another source of heterogeneity is the reproducibility of the IMT measurement procedures, as shown in Table 2. The reproducibility has clearly improved over time.
|
By means of a meta-regression analysis (see Methods), we determined the influence of sample properties (mean age, duration of follow-up, mean deviation, and SD of IMT) and measurement properties (maximal or mean IMT, near and far wall versus far wall IMT, 1 side versus both sides, position and length of the examined carotid segment, and reproducibility) on the hazard estimate between studies. The duration of follow-up was the strongest (and the only significant) predictor of the magnitude of the log-HR for stroke per IMT difference (P=0.0424). However, this result withstood neither multivariate modeling nor adjustment for multiple testing.
Assessment of the Relationship Between IMT and Hazard
To challenge the assumption of a linear IMT–hazard relationship, we used the graphical method described above (Figure 2). The confidence bands cross the dashed lines, which indicate that the linear model violates the assumption of linearity in the 2 largest studies4–6 (data not shown). For the majority of populations, the slope of the relationship between IMT and hazard is steeper for low to average IMT values (as can be seen in Figure 2), with the slope for the highest IMT categories being relatively moderate. In most of the studies, however, the linear model corresponds well to the categorized analysis for higher IMT values.
|
HR per IMT Difference Meta-Analysis
As stated in the Methods, different assumptions of the underlying distributions lead to different reference units of IMT. In Figure 3 we display the risk of MI and stroke in the 5 remaining studies, expressed as HR per CCA IMT difference of 1 SD of baseline IMT. Figure 4 shows the plots that correspond to HR per absolute CCA IMT difference of 0.1 mm. Figures 5 and 6
display the same calculations for the HRs of fully adjusted models. In comparison of the plots, we observed that HRs decrease when adjusted for a complete set of cardiovascular risk factors as compared with age- and sex-adjusted estimates. Furthermore, we can see that the HR per IMT difference is slightly higher for stroke than for MI. The heterogeneity is larger for the estimates of MI risk. In an attempt to reduce heterogeneity by strict selection of studies that used identical, or at least comparable, definitions of the carotid segments (Figure 1), we calculated the overall estimates shown in Table 3. The estimates or relative risk for stroke still remained higher than for MI, and the studies are more heterogeneous for MI.
|
|
|
|
|
The trim-and-fill procedure found no evidence of publication bias.
| Discussion |
|---|
|
|
|---|
Low Specificity of IMT Location for a Clinical End Point
Most studies that investigated the risk of both myocardial and cerebrovascular end points resulted in similar HRs for both end points. The ARIC study4,6 revealed the largest differences (eg, MI versus stroke risk, 1.13 versus 1.21, translated into risk ratio per 0.10-mm CCA IMT difference, adjusted for age and sex). The Cardiovascular Health Study5 shows a higher relative risk for stroke, too. The other 3 large studies7,8,10–12 resulted in an insignificantly higher relative risk of MI. The difference found in this meta-analysis is just statistically significant, as the CIs do not overlap for the models of HR per SD IMT difference. However, the difference is small. This lack of specificity for the target organ has been described by other authors,5–6 who drew the conclusion that early atherosclerotic changes are noncausal and nonspecific markers of atherosclerotic complications.6
Heterogeneity of Ultrasound Protocols
We determined multiple sources of heterogeneity between the studies. These include the details of the ultrasound protocol, namely the precise definitions of the carotid segments investigated, the use of mean or maximal IMT, the measurement of near and far wall or only far wall IMT, and whether IMT was measured only on 1 side or on both sides. Reproducibility has improved considerably over the past decades. The study populations differed with respect to age distribution, selection criteria, and length of follow-up. With the exception of the duration of follow-up, none of the sample and measurement properties significantly explained the heterogeneity between the studies. The fact that parameters such as the position of the carotid segment (Figure 1) do not significantly influence the hazard estimate puts the multiple differences between the ultrasound protocols into perspective. However, with only 5 studies to compare, the statistical power of the meta-regression is limited.
There are good arguments for and against each IMT measurement protocol used. Nevertheless, we strongly recommend the use of a standardized protocol for future studies. A recent consensus paper24 has given some guidance on this. It is noteworthy that this consensus paper recommends measurement of IMT in plaque-free areas. The ultrasound protocols of the epidemiological studies discussed in the present study do not explicitly state how to deal with plaques, but, given that all of them use clear anatomic landmarks, it is likely that, in the IMT measurements, plaque thickness would have been incorporated and measurements would not all have been in plaque-free regions. However, in a general population-based survey with low prevalence of plaque, this is unlikely to be a major confounding factor. For epidemiological studies of low-risk populations, we suggest adherence to exact anatomic definitions. Whenever risk populations are surveyed, the consensus paper24 recommendation applies. To ensure clinically meaningful results, we suggest adoption of the ultrasound protocol of one of the large studies mentioned in the present study (Table 1), preferably of a study with a target population similar to that of the intended study. If a single protocol is established in future, it should have demonstrated good reproducibility and a solid database that demonstrates the clinical relevance of IMT measured by this method. Both are true for the measurement and analysis protocol of CAPS.12
Nonlinear Relationship Between IMT and Risk
When IMT is used as a risk predictor, it is often assumed that the relationship between IMT and risk is linear. Our analysis demonstrates that the relationship between CCA IMT and the relative risk of vascular events is not strictly linear in most populations. This deviation from linearity was first observed by the ARIC investigators4,6 and found to be statistically significant for MI (P=0.002 for men, P=0.04 for women)4 and nearly significant for stroke (P=0.08 for men, P=0.06 for women).6 To understand this relationship, we have to be aware of the link with age: Young individuals with increased IMT are at considerably lower absolute risk but higher relative risk of vascular events.12 Furthermore, individuals with a high risk factor load (and consequently higher IMT) are more frequently under a physicians surveillance and pharmacological intervention, which may reduce their individual risk more than the corresponding IMT.
Vascular events are rare in young individuals, which makes IMT particularly attractive as an end point in epidemiological and treatment studies in young populations. The nonlinear relationship has major implications for the design of IMT studies in young populations, as it limits the transferability of the results from older samples. It is of the utmost importance to investigate the actual risk linked with increased IMT. More studies are needed to provide a larger database on the association between IMT and vascular risk in young individuals. As long as the data on the risk in young individuals are limited, we recommend that IMT studies on young populations should adopt the ultrasound protocol of a study that provides prospective data (at the time of submission there was only one such study12) to be able to express the IMT findings as a relative risk of clinical events.
Limitations
Our review and meta-analysis has a number of potential limitations. First, there is the issue of potential publication bias. We have tried to identify all those studies that calculated the relative risks of future cerebrovascular or cardiovascular events or death that depend on individual baseline carotid IMT. By means of an extensive Medline-based search, we aimed to retrieve all published data. Given that such observational studies require large samples and long follow-up periods to facilitate the detection of a considerable number of events, it seems very unlikely that a study of this dimension will have been published and not identified by our search. The trim-and-fill procedure, a statistical technique to detect imbalances within the data that possibly indicate the existence of unpublished studies, found no evidence of such an imbalance in the data discussed here.
Second, the graphical method used to investigate the linearity of the IMT-risk relationship relies on the assumption of section-wise linearity. The exact underlying function is determined by polynomial regression, where the maximal power of the function depends on the number of categories. With a minimum of 3 categories, the spline curves can be second-order, which is sufficient to detect the nonlinearity described. A more complex function may be overlooked.
A third potential limitation stems from the multiple sources of heterogeneity between the studies. To obtain compound estimates, we had to take into account the large number of heterogeneities on the one hand and the possible selection bias on the other. Selection bias could be introduced if the heterogeneities were informative, eg, if the proportion of advanced atherosclerosis and the frequency of plaques in a population could influence the authors decision to report mean or maximal IMT, depending on how well either measurement algorithm accounted for plaques. To account for both problems, we provided calculations and forest plots for both selected and unselected studies.
Fourth, a potential limitation is the use of HR per (linear) IMT difference. The use of categorized data has other problems, as the choice of the cutoff points must not be data dependent to avoid introduction of a new bias factor. For the higher age group, however, the linear model still fits surprisingly well. The results of this meta-analysis must be interpreted with caution if young populations with low IMT are of interest.
| Acknowledgments |
|---|
We thank Daniel H. OLeary (Department of Radiology, Tufts–New England Medical Center, Boston, Mass), Richard A. Kronmal (Department of Biostatistics, University of Washington, Seattle, Wash), and Lloyd Chambless (Department of Biostatistics, School of Public Health, University of North Carolina, Chapel Hill, NC) for generously sharing their data with us. We thank Peter M. Rothwell and Enrico Flossmann (both from the Department of Clinical Neurology, University of Oxford, Oxford, United Kingdom) for supplying the macro program for forest plots. One of our data sources, the ARIC study is conducted and supported by the National Heart, Lung, and Blood Institute in collaboration with the ARIC Investigators. This manuscript was prepared with a limited-access dataset obtained by the National Heart, Lung, and Blood Institute and does not necessarily reflect the opinions or views of the ARIC or the National Heart, Lung, and Blood Institute.
Disclosures
None.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. R. Folsom, R. A. Kronmal, R. C. Detrano, D. H. O'Leary, D. E. Bild, D. A. Bluemke, M. J. Budoff, K. Liu, S. Shea, M. Szklo, et al. Coronary Artery Calcification Compared With Carotid Intima-Media Thickness in the Prediction of Cardiovascular Disease Incidence: The Multi-Ethnic Study of Atherosclerosis (MESA) Arch Intern Med, June 23, 2008; 168(12): 1333 - 1339. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. J. O'Donnell, S. Demissie, M. Kimura, D. Levy, J. P. Gardner, C. White, R. B. D'Agostino, P. A. Wolf, J. Polak, L. A. Cupples, et al. Leukocyte Telomere Length and Carotid Artery Intimal Medial Thickness: The Framingham Heart Study Arterioscler. Thromb. Vasc. Biol., June 1, 2008; 28(6): 1165 - 1171. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. D. Solomon, J. Wittes, P. V. Finn, R. Fowler, J. Viner, M. M. Bertagnolli, N. Arber, B. Levin, C. L. Meinert, B. Martin, et al. Cardiovascular Risk of Celecoxib in 6 Randomized Placebo-Controlled Trials: The Cross Trial Safety Analysis Circulation, April 22, 2008; 117(16): 2104 - 2113. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Chonchol, H. Gnahn, and D. Sander Impact of subclinical carotid atherosclerosis on incident chronic kidney disease in the elderly Nephrol. Dial. Transplant., April 3, 2008; (2008) gfn021v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Sanz, P. R. Moreno, and V. Fuster The year in atherothrombosis. J. Am. Coll. Cardiol., March 4, 2008; 51(9): 944 - 955. [Full Text] [PDF] |
||||
![]() |
V. Stangl, V. Witzel, G. Baumann, and K. Stangl Current diagnostic concepts to detect coronary artery disease in women Eur. Heart J., March 2, 2008; 29(6): 707 - 717. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. C. Crawford, A. S. Nord, M. D. Badzioch, J. Ranchalis, L. A. McKinstry, M. Ahearn, C. Bertucci, C. Shephard, M. Wong, M. J. Rieder, et al. A common VLDLR polymorphism interacts with APOE genotype in the prediction of carotid artery disease risk J. Lipid Res., March 1, 2008; 49(3): 588 - 596. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. R. Feldman, G. J. Bosl, J. Sheinfeld, and R. J. Motzer Medical Treatment of Advanced Testicular Cancer JAMA, February 13, 2008; 299(6): 672 - 684. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Stern Are We Getting Nearer to Screening for Atherosclerosis? Circulation, January 1, 2008; 117(1): 122 - 126. [Full Text] [PDF] |
||||
![]() |
K. O. Stumpe, E. Agabiti-Rosei, T. Zielinski, D. Schremmer, J. Scholze, P. Laeis, P. Schwandt, and M. Ludwig Original Research: Carotid intima-media thickness and plaque volume changes following 2-year angiotensin II-receptor blockade. The Multicentre Olmesartan atherosclerosis Regression Evaluation (MORE) study Therapeutic Advances in Cardiovascular Disease, December 1, 2007; 1(2): 97 - 106. [Abstract] [PDF] |
||||
![]() |
M. Luque-Ramirez, C. Mendieta-Azcona, F. Alvarez-Blasco, and H. F. Escobar-Morreale Androgen excess is associated with the increased carotid intima-media thickness observed in young women with polycystic ovary syndrome Hum. Reprod., December 1, 2007; 22(12): 3197 - 3203. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Yumino and T. D. Bradley Pathogenesis of Atherosclerosis: Is Obstructive Sleep Apnea the New Kid on the Block? Am. J. Respir. Crit. Care Med., October 1, 2007; 176(7): 634 - 635. [Full Text] [PDF] |
||||
![]() |
L. F. Drager, L. A. Bortolotto, A. C. Figueiredo, E. M. Krieger, and G. Lorenzi-Filho Effects of Continuous Positive Airway Pressure on Early Signs of Atherosclerosis in Obstructive Sleep Apnea Am. J. Respir. Crit. Care Med., October 1, 2007; 176(7): 706 - 712. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Rundek, R. D. Brook, and J. D. Spence Letter by Rundek et al Regarding Article, "Prediction of Clinical Cardiovascular Events With Carotid Intima-Media Thickness: A Systematic Review and Meta-Analysis" Circulation, August 28, 2007; 116(9): e317 - e317. [Full Text] [PDF] |
||||
![]() |
W. Koenig Cardiovascular Biomarkers: Added Value With an Integrated Approach? Circulation, July 3, 2007; 116(1): 3 - 5. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||