| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Circulation. 2007;115:2722-2730.)
© 2007 American Heart Association, Inc.
Epidemiology |
From the Department of Biostatistics, University of Washington, Seattle (R.A.K., R.L.M.); Division of Cardiology, HarborUCLA Medical Center, Los Angeles, Calif (R.D.); Departments of Medicine and Epidemiology, Columbia University, New York, NY (S.S.); Department of Cardiology, Johns Hopkins University, Baltimore, Md (J.A.L.); Department of Medicine, University of Vermont, Burlington (M.C.); National Institutes of Health, Bethesda, Md (D.E.B.); and Department of Public Health Sciences, Wake Forest University, Winston-Salem, NC (G.L.B.).
Correspondence to Robyn McClelland, Collaborative Health Studies Coordinating Center, Department of Biostatistics, University of Washington, Bldg 29, Suite 310, 6200 NE 74th St, Seattle, WA 98115. E-mail rmcclell{at}u.washington.edu
Received November 6, 2006; accepted March 26, 2007.
| Abstract |
|---|
|
|
|---|
Methods and Results Follow-up CAC measurements were available for 5756 participants with an average of 2.4 years between scans. The incidence of newly detectable CAC averaged 6.6% per year. Incidence increased steadily across age, ranging from <5% annually in those <50 years of age to >12% in those >80 years of age. Median annual change in CAC for those with existing calcification at baseline was 14 Agatston units for women and 21 Agatston units for men. Most traditional cardiovascular risk factors were associated with both the risk of developing new incident coronary calcium and increases in existing calcification. These included age, male gender, white race/ethnicity, hypertension, body mass index, diabetes mellitus, glucose, and family history of heart attack. Factors also existed that were related only to incident CAC risk, such as low- and high-density lipoprotein cholesterol and creatinine. Diabetes mellitus had the strongest association with CAC progression for blacks and the weakest for Hispanics, with intermediate associations for whites and Chinese.
Conclusions This is the first large multiethnic study reporting on the incidence and progression of CAC. Standard coronary risk factors were generally related to both CAC incidence and progression. Whites had more incident CAC and CAC progression than the other 3 racial/ethnic groups. Except for diabetes mellitus, risk factor relationships were similar across racial/ethnic groups.
Key Words: arteries calcium coronary disease epidemiology imaging risk factors
| Introduction |
|---|
|
|
|---|
Clinical Perspective p 2730
These previous studies on CAC progression have been limited by issues such as small sample sizes,1012,1519 retrospective designs,1012,17 and/or restricted populations such as high-risk hypertensives15,16 or familial hypercholesterolemics.18 Certain studies consisted predominantly of men,11,13,18 and many focused only on particular risk factors of interest such as cholesterol levels or statin use,12,13,18 obesity,14 or renal dysfunction.15 The Multi-Ethnic Study of Atherosclerosis (MESA) provides a unique opportunity to study the risk factors for progression of coronary artery calcification in a large, gender-balanced, multiethnic, nonreferred, asymptomatic cohort. We consider here traditional cardiovascular risk factors as they relate to the risk of developing incident detectable CAC and to the magnitude of CAC progression given that some calcium is present.
| Methods |
|---|
|
|
|---|
Measurement of CAC
CAC was measured with electron-beam computed tomography (EBCT) at 3 field centers and multidetector row helical computed tomography (MDCT) at 3 field centers. Each participant was scanned twice consecutively, and these scans were read independently at a centralized reading center. We collimated and reconstructed at least 40 scan slices using a 3.0-mm slice thickness for EBCT scanners and 2.5-mm slice thickness for MDCT scanners. The methodology for acquisition and interpretation of the scans has been reported previously.21 The results of the 2 scans were averaged to provide a more accurate point estimate of the amount of calcium present. The amount of calcium was quantified with the Agatston scoring method.22 Calcium scores were adjusted with a standard calcium phantom that was scanned along with the participant.23 The phantom contained 4 bars of known calcium density and provided a way to calibrate the x-ray attenuation level between measurements conducted on different machines. This was important because scanners were changed between baseline and follow-up at 3 of the 6 sites. Any detectable calcium was defined as a CAC score >0; a minimum focus of calcification was based on at least 4 contiguous voxels, which resulted in identification of calcium of 1.15 mm3 for the MDCT scanners and 1.38 mm3 for the EBCT scanners.21 The nominal section thickness was 3.0 mm for EBCT scanners and 2.5 mm for MDCT scanners. Interobserver agreement and intraobserver agreement were found to be very high (
=0.93 and 0.90, respectively). Follow-up CAC measurements were performed on half the cohort (randomly selected) at a second examination (September 2002 through January 2004) and the other half of the cohort at a third examination (March 2004 through July 2005) at an average of 1.6 and 3.2 years after the participants first examination, respectively. As for the baseline examination, the results from 2 consecutive scans were averaged. A full characterization of the distribution of CAC at baseline by age, gender, and race has been published previously.24
Measurement of Covariates
Information on demographics, smoking, medical conditions, and family history was collected by questionnaire at the initial examination. Height and weight were also measured at the baseline examination, and blood was drawn for measurements, including lipids, inflammation, fasting glucose, fibrinogen, and creatinine. Resting blood pressure was measured 3 times in the seated position, and the average of the last 2 measurements was used in the analysis. Medication use was determined by questionnaire. Additionally, the participant was asked to bring to the clinic containers for all medications used during the 2 weeks before the visit. The interviewer then recorded the name of each medication, the prescribed dose, and frequency of administration from the containers. Further details and additional references have been published previously.20
Statistical Analysis
All participants with both a baseline and a follow-up CAC measurement were included in the analysis. The presence of CAC was defined as an Agatston score >0. Sensitivity analyses revealed that using various alternative cut points (such as 10 Agatston units) did not materially alter the risk factor relationships (data not shown). Progression of CAC was defined in 2 ways: incident CAC defined as detectable CAC at the follow-up examination (either examination 2 or 3) in a participant free of detectable CAC at examination 1 and change in CAC score in participants who had detectable CAC at examination 1. These 2 end points were modeled separately. Risk factors considered include age, gender, race/ethnicity, education, income, systolic and diastolic blood pressures, use of antihypertensive medications, diabetes status (normal, impaired fasting glucose 100 to 125 mg/dL, untreated diabetes mellitus [>125 mg/dL fasting glucose], and treated diabetes mellitus), smoking (never, former, current), pack-years of smoking, body mass index, LDL and HDL cholesterol, triglycerides, use of lipid-lowering medication, fibrinogen, creatinine, and C-reactive protein (CRP). Triglycerides and CRP were highly skewed and were log-transformed. For creatinine, there was some evidence of nonlinearity; hence, creatinine was categorized with cut points selected empirically on the basis of observed risk of incident CAC. Three groups were defined as creatinine
0.9, 1.0, and
1.1 mg/dL.
Relative risk regression25 was used to model the probability of incident-detectable CAC among those free of CAC at examination 1. That is, the probability of incident CAC was modeled as a function of covariates using a generalized linear model with log link and binomial error distribution. In cases in which the model failed to converge with the binomial error (
5% of the models), we substituted gaussian error and used robust standard error estimates. We used relative risk regression rather than logistic regression because the incidence of new calcification is not rare (occurring in >16% of participants during the follow-up), so the odds ratio is an overestimate of the relative risk. Age- and gender-adjusted models for each risk factor were estimated, followed by a multivariable model constructed via a backward elimination variable selection process. The backward selection process started with all the candidate variables (regardless of statistical significance) and sequentially removed variables with P>0.05, beginning with the least significant variable. The time between scans was included as a covariate in all models, and interactions of each risk factor with gender and race/ethnicity were tested.
Among those with some detectable CAC at examination 1, we defined progression as the absolute difference between follow-up and examination 1 CAC, and this was treated as a continuous end point. The modeling strategy for this end point was analogous to the models for incident CAC; however, outliers in the calcium measurements posed an additional analytic challenge. To overcome this, we fit robust linear regression models. The robust regression algorithms in Stata26 were used, which iteratively downweight observations with large residuals.27 In these models,
300 participants (10.7%) were given extremely low or zero weight in the analysis. Generally, these were participants with very fast progression relative to their predicted progression based on risk factors. Scanner changes at some of the sites also influenced progression magnitude, and a term for scanner pair was included in all the models for progression. The multivariable model was again constructed via a backward elimination variable selection process as described above for the incident CAC analysis. Our primary analysis for absolute progression did not control for baseline CAC score. In a secondary analysis, we explored which risk factor associations persisted after adjustment for baseline score.
In addition to the raw change in calcium, it was also of interest to study the relative change in calcium (relative to baseline). An initial analysis considered the percentage change in CAC; however, this end point is highly influenced by very small baseline calcium scores, for which a very small absolute increase in CAC may result in a huge percentage increase. A common approach to this problem is to exclude those with small initial calcium, say CAC <10; however, this would have resulted in exclusion of >400 selectively low-risk participants. We expressed relative change in CAC as the change from baseline to follow-up in the log of CAC plus a constant: ie, ln(CACFU+25)ln(CACBL+25). The constant was chosen so that the resulting change was fairly symmetric and normally distributed. The addition of the constant prevents the very small baseline CAC scores from leading to unreasonably large relative increases, and use of the log-scale downweights the exceptionally large follow-up CAC scores.
The authors had full access to and take full responsibility for the integrity of the data. All authors have read and agree to the manuscript as written. All statistical analyses were performed with either Stata version 9.126 or R version 2.3.1.28
| Results |
|---|
|
|
|---|
9% had treated diabetes mellitus, and 40% had a family history of heart attack. Participants who were excluded because of a lack of follow-up calcium were slightly older (64 versus 62 years), were more likely to have CAC at baseline (56% versus 49%), had higher systolic blood pressure (130 versus 126 mm Hg), and were more likely to be diabetic (19% versus 13%) or current smokers (16.6% versus 12.4%). Among the 2948 participants (51.2%) without detectable CAC at baseline, 475 (16.1%) went on to develop incident CAC over an average of 2.4 years of follow-up, for an estimated incidence rate of 6.6% per year. Figure 1 illustrates the linearly increasing yearly rate of incident CAC as a function of baseline age. The incidence was <5% per year at younger ages and >12% per year at the highest ages. Age differences were similar across gender and race/ethnicity subgroups. Table 1 provides yearly rates of incident CAC by gender and race/ethnicity. In women, black participants had the highest yearly rate of incident CAC at 6.3%, whereas Chinese women had the lowest at 4.7%. White and Hispanic women were intermediate at 5.6% and 6.0%, respectively. After age adjustment, these racial/ethnic differences among women were nonsignificant. With the exception of Chinese, men had higher incidence rates than women (P<0.001). White men had the highest rate at 10.3%, followed by blacks at 8.2% (P=0.032 compared with white men after age adjustment) and Hispanics at 7.2% (age-adjusted P=0.039 compared with white men). Chinese men were again the lowest, with a rate of only 4.4% (age adjusted P<0.001 compared with white men).
|
|
Table 2 displays the associations of traditional cardiovascular risk factors with the risk of incident detectable CAC. Each additional 10 years of age was associated with a 39% higher risk (95% CI, 28 to 52), similar to comparisons of men relative to women (43% higher rate for men adjusted for age and length of follow-up; 95% CI, 22 to 67). After adjustment for age, gender, and length of follow-up, many traditional cardiovascular risk factors were associated with an increased risk of incident CAC, including white race, higher body mass index, high blood pressure, abnormal lipids (higher LDL or triglycerides, lower HDL, or use of lipid-lowering medication), glucose disorders (impaired fasting glucose, untreated or treated diabetes mellitus), family history of heart attack, higher fibrinogen, higher CRP, and low creatinine. Education, family income, and pack-years of smoking were not related to incident CAC risk. Backward selection through all the variables resulted in a model including age, gender, race/ethnicity, body mass index, use of antihypertensive medications, LDL cholesterol, triglycerides, use of lipid-lowering medications, diabetes status, family history of heart attack, and creatinine group. Current and former smokers had higher incident CAC rates than never smokers, but this difference was not statistically significant once other risk factors were considered. Fibrinogen and CRP were no longer significantly associated with CAC incidence after adjustment for body mass index. If triglycerides were not included as a covariate, HDL cholesterol entered instead with a protective association. The risk of incident CAC did not differ by type of scanner (MDCT versus EBCT) or by study site (data not shown).
|
Table 3 provides summary characteristics of annualized CAC progression. Three-hundred-fifty-eight (12.8%) participants revealed a negative calcium change on follow-up, and 691 (24.6%) had a very small amount of progression (<10 Agatston units per year). The majority of participants (1367 or 48.7%) had modest progression (between 10 and 99 Agatston units per year). Two-hundred-forty-nine (8.9%) had substantial progression (between 100 and 199 Agatston units annually) and 143 (5.1%) had extremely large progression (between 200 to >700 Agatston units per year). Median progression was 18 Agatston units overall (14 for women, 21 for men), with a mean of 46 Agatston units (36 for women, 54 for men). The distribution was heavily skewed because of the presence of several large progressors. Figure 2 is a histogram showing the distribution of annual change in CAC. Figure 3 illustrates a smoothed plot of annual change in calcium by age among the 2808 participants with detectable CAC at baseline. The increasing trend is linear across age and was similar for men and women, although men had
17 to 20 Agatston units more progression per year on average than women at any given age.
|
|
|
Table 4 displays the age-, gender-, scanner type, and follow-up timeadjusted associations of risk factors with CAC progression, followed by the multivariable results. Each additional year of follow-up was associated with an extra 17.5 Agatston units of progression on average (95% confidence interval, 14.6 to 20.4), and each additional 10 years of baseline age was associated with 8.8 more Agatston units of progression (95% confidence interval, 6.4 to 11.2). Men had an average of 10.9 more Agatston units of progression than women (95% confidence interval, 6.3 to 15.5). Other risk factors associated with more progression included white race/ethnicity, greater body mass index, higher blood pressure (particularly systolic), higher triglycerides or use of lipid-lowering medication, pack-years of smoking, diabetes mellitus, higher CRP, and a family history of heart attack. In the multivariable models, factors associated with increased progression included time between scans, age at baseline, male gender, white race/ethnicity, higher body mass index, high blood pressure, family history of heart attack, and diabetes mellitus. The amount of progression did not differ by type of scanner (MDCT versus EBCT) or by study site (data not shown). The R2 for the multivariable model is
12%; thus, much of the variation remains unexplained. There was also a significant interaction between race and diabetes mellitus (P=0.0002), indicating that the association of diabetes mellitus with progression was strongest for blacks than for the other races and weakest among Hispanics, with whites and Chinese intermediate. Treated diabetes mellitus was associated with an estimated 48 more Agatston units of progression in blacks, 28 more Agatston units in whites, 19 more Agatston units in Chinese, and only 7 more Agatston units in Hispanics.
|
Forcing in the same risk factors as the fully adjusted robust regression model (model 2 in Table 4), adjustment for baseline CAC removed the association of age, gender, race/ethnicity, smoking, and systolic blood pressure; body mass index, family history of heart attack, and diabetes mellitus remained significant. Within the subset of participants who had only a small amount of CAC (1 to 100 Agatston units) at baseline, factors that remained strong modifiable predictors of progression included blood pressure and glucose (data not shown). Factors that were associated with a higher relative increase in CAC (expressed relative to their baseline CAC) were longer time between scans and higher systolic blood pressure, diabetes mellitus, and CRP (even after adjustment for body mass index, which was not significant).
After exclusion of those taking medication for diabetes mellitus (335 who had detectable CAC at baseline and a follow-up CAC measurement), glucose level (log transformed) was significantly associated with progression. Among the 327 participants with treated diabetes mellitus, self-reported information on diabetes mellitus duration was available for 271. In this subset, diabetes duration was positively associated with calcium progression, with each additional 10 years of diabetes duration associated with 24 more Agatston units of progression (95% confidence interval, 7 to 416; P=0.005) after adjustment for length of follow-up, scanner, age, and gender. Aside from age and follow-up time, none of the other risk factors were significant in this subset.
| Discussion |
|---|
|
|
|---|
Diabetes mellitus was the strongest risk factor for CAC progression. In fact, glucose was associated with higher risk of incident CAC even at levels well below the standard thresholds for diabetes mellitus or impaired fasting glucose and among nondiabetic participants. Among treated diabetics, duration of diabetes was the only significant risk factor for CAC progression after adjustment for age, gender, scanner, and duration of follow-up. However, because of the small number of such patients (n=281 with both CAC change and diabetes duration), our power for detecting other associations was limited.
The associations that persisted after adjustment for baseline CAC were body mass index, family history of heart attack, diabetes mellitus, and glucose. This indicates that these factors have an association with CAC that is not yet completely reflected by their lifetime CAC accumulation. This is in contrast to the more chronic associations of age, gender, race/ethnicity, hypertension, and smoking, the influence of which is already manifest in the baseline CAC level. Although both fibrinogen and CRP were related to an increased risk of incident CAC, this became nonsignificant after adjustment for body mass index.
In the present study, we observed a modest U-shaped relationship of creatinine with the risk of incident CAC. Specifically, those with creatinine levels of 1 mg/dL had the lowest risk, those with creatinine <1 mg/dL had significantly higher risk, and those with creatinine >1 mg/dL had higher risk, but not significantly so. Renal dysfunction has been associated with a greater amount and progression of coronary calcium.15,29,30 Very few participants in MESA have renal dysfunction, which may explain why the association of higher creatinine levels with risk was nonsignificant. Low levels of creatinine can indicate the presence of severe liver disease or decreased muscle mass, but it seems unlikely that this explains the increased risk of incident CAC in this relatively healthy cohort. This finding was data driven and may simply be a false-positive. Further study is needed to confirm or refute this association.
Comparison of our results with previous studies on progression is challenging because of differing analytic strategies. Almost all previous studies have controlled for baseline calcium in the models, which may obscure important associations. These studies tend to conclude that the strongest (or even the only) predictor of calcium change is baseline calcium level and that factors such as age and gender are not associated.10,31 The rate of progression before baseline will be a strong determinant of the baseline CAC score. Furthermore, it is reasonable to assume that the rate of progression before baseline is predictive of future progression rates. In other words, baseline CAC is not a confounder but part of the causal pathway for the relationship between long-term risk factors and CAC progression. For example, despite the fact that men clearly have more progression than women, once the baseline CAC is controlled for, no gender difference exists. The reason is that they already have been progressing more and already have more CAC. They also will continue to progress more rapidly. It is difficult to interpret the results from such a model because the effects that the variables have had up until baseline are conditioned out. If baseline CAC is included in the model, only associations that are acute (eg, a new medication, recent infection, newly diagnosed condition) or possibly those that have a strong association that changes over time will be easily detectable.
The choice of scale for the analysis of continuous change in CAC score is a problematic analytic issue, one that makes comparison across various studies in the literature difficult. Different choices (eg, raw change versus percentage change versus change in log calcium plus a constant) yield different conclusions. We found certain standard analyses, eg, simply modeling raw change or percent change, to be heavily influenced by outliers. To overcome this issue, we used robust regression techniques when modeling absolute change and a logarithmic transformation of calcium plus a constant when modeling relative change. The addition of the constant reduces the impact of small baseline calcium, resulting in very large relative changes, and the log scale decreases the association of the extremely large follow-up calcium measurements.
Both antihypertensive medication and lipid-lowering medication use were associated with a significantly increased risk of incident CAC and significantly more CAC progression. This does not imply that these medications have a negative impact; more likely, their use indicates more severe underlying hypertension or dyslipidemia, which may have an impact beyond that of the concurrently measured blood pressure or cholesterol values. These may be normal as a result of treatment despite underlying disease. Several previous studies have found evidence that use of statins slows, or in some cases even reverses, the progression of calcium11,12,32 among those with high cholesterol. Those studies were based on relatively small sample sizes, and although the participants were free of clinical coronary artery disease, they were referred for assessment of their CAC. These participants tended to have more risk factors and more calcium at baseline than the participants in MESA. Additionally, the findings of statin benefit refer to those with hypercholesterolemia, and the untreated participants in MESA do not necessarily have high cholesterol.
A unique feature of our data is the ability to study racial/ethnic differences. Whites had far more incident CAC and CAC progression than the other 3 racial/ethnic groups. Risk factor relationships, however, were largely the same across racial/ethnic groups. For incident CAC, there was no evidence that the effect of risk factors differed by race. For CAC progression, although most risk factor relationships were similar across races, the effect of diabetes mellitus was strongest in blacks and weakest in Hispanics. This is consistent with at least 1 previous study33 reporting decreased congestive heart disease and cardiovascular disease mortality in Hispanic diabetics compared with non-Hispanic white diabetics.
Strengths of the MESA study include the large sample size, inclusion of 4 racial/ethnic groups, and the community-based (as opposed to referral-based) nature of the sample. Additionally, the prospective nature of the study allowed calcium measurement and risk factor assessment to proceed in a standardized manner. The sample did not include participants with clinical cardiovascular disease; hence, results on progression of calcium cannot be generalized to that population. Participants without a follow-up CAC tended to have a somewhat worse risk factor profile on average than included participants. This may have weakened certain risk factor associations. Some differences in scanning methodology exist between our study and previously published studies. The field of view for the MESA scans was 350 mm, which is larger than previous studies, resulting in larger pixels. We also required calcium to be present in 4 adjacent pixels to register as a calcified area, whereas other studies required between 2 and 4 pixels. We used backward selection to construct our multivariable models, and a risk exists that
1 confounding variables (ie, variables that might influence the association of other variables despite not being statistically significant themselves) were removed. Many statistical comparisons were performed in this study, and the some of the associations may be false-positives. A considerable amount of nonbiological variability exists in the measurement of CAC because of computed tomography technical factors. As a result, we may have been unable to detect certain more modest associations as a result of the presence of noise. In part, because of this variability, the R2 for our multivariable model of CAC change was relatively low at 12%. However, we were able to demonstrate several factors that are strongly associated with CAC progression.
| Conclusions |
|---|
|
|
|---|
| Acknowledgments |
|---|
Sources of Funding
This research was supported by contracts N01-HC-95159 through N01-HC-95165 and N01-HC-95169 from the National Heart, Lung, and Blood Institute.
Disclosures
None.
| References |
|---|
|
|
|---|
6.6% per year, ranging from <5% per year below age 50 to >12% per year at ages 80 and over. Median annual change in CAC for those with existing calcification at baseline was 14 Agatston units for women and 21 for men. In addition, although whites have more progression of CAC than do blacks, Hispanics, or Chinese, risk factor associations with progression are largely similar across different races/ethnicities.
Related Article:
Circulation 2007 115: 2683.
This article has been cited by other articles:
![]() |
M. M. Vigoda, L. I. Rodriguez, E. Wu, K. Perry, R. Duncan, D. J. Birnbach, and D. A. Lubarsky The Use of an Anesthesia Information System to Identify and Trend Gender Disparities in Outpatient Medical Management of Patients with Coronary Artery Disease Anesth. Analg., July 1, 2008; 107(1): 185 - 192. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Gottlieb and J. A.C. Lima Screening High-Risk Patients With Computed Tomography Angiography Circulation, March 11, 2008; 117(10): 1318 - 1332. [Full Text] [PDF] |
||||
![]() |
S. Schroeder, S. Achenbach, F. Bengel, C. Burgstahler, F. Cademartiri, P. de Feyter, R. George, P. Kaufmann, A. F. Kopp, J. Knuuti, et al. Cardiac computed tomography: indications, applications, limitations, and training requirements: Report of a Writing Group deployed by the Working Group Nuclear Cardiology and Cardiac CT of the European Society of Cardiology and the European Council of Nuclear Cardiology Eur. Heart J., February 2, 2008; 29(4): 531 - 556. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. V. Anand, E. Lim, D. Darko, P. Bassett, D. Hopkins, D. Lipkin, R. Corder, and A. Lahiri Determinants of Progression of Coronary Artery Calcification in Type 2 Diabetes: Role of Glycemic Control and Inflammatory/Vascular Calcification Markers J. Am. Coll. Cardiol., December 4, 2007; 50(23): 2218 - 2225. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||