Ability of a biomarker-based score to predict death from circulatory disease and cancer in NHANES III

Abstract

Background

A score based on serum concentrations of C-reactive protein (CRP), albumin, gamma-glutamyl transferase (GGT), and HDL cholesterol was positively associated with death from cancer, circulatory disease, and all-cause mortality. We replicated this in the third National Health and Nutrition Examination Survey (NHANES III), a US nationally representative survey conducted between 1988–1994.

Methods

Baseline measurements of CRP, albumin, GGT, and HDL were available for participants with mortality follow-up (n=13,056). A biomarker score, ranging 0–4, was created by adding number of markers with abnormal values (cut-off: CRP>10mg/L, albumin<35mg/L, GGT>36U/L, HDL<1.04mmol/L). Its association with mortality was analyzed with multivariate Cox proportional hazards models.

Results

The score was positively associated with death from all causes, cancer and circulatory disease [e.g. HR all-cause mortality: 1.21 (95% CI: 1.09, 1.35), 1.92 (1.67, 2.20), 3.38 (2.62, 4.36), and 7.93 (5.77, 10.89), for score 1, 2, 3, 4 vs.0]. These patterns were found across the Charlson Comorbidity Index (CCI). Where CCI =3, risk of cancer death was 1.09 (0.93, 1.28), 1.81 (1.43, 2.29), 4.67 (3.05, 7.14), and 6.97 (5.32, 9.14) for score 1, 2, 3, 4 vs. 0. No effect-modification by sex or race/ethnicity was observed.

Conclusions

These findings correlate with results from a Swedish study. This biomarker-based score could help clinicians make decisions in prevention and disease management.

Background

In a large prospective cohort study from Sweden, the combination of easily and inexpensively measurable biomarkers clearly predicted mortality in persons aged 50+ years [1]. This mortality score included serum levels of C-reactive protein (CRP), albumin, gamma-glutamyltransferase (GGT), and high-density lipoprotein (HDL) cholesterol, all reflecting possible mechanisms contributing to early death: inflammation [2, 3], hepatic function [4], and lipid metabolism [5, 6]. The score was positively and progressively associated with all-cause mortality as well as cancer and circulatory disease-specific death. Even among cancer patients with no other co-morbidities, according to the Charlson Comorbidity Index (CCI), the score was predictive for mortality [1]. The mortality score’s ability to predict death is important, especially given the need for an assessment of functional rather than chronological age when treating cancer patients, as recently expressed by the European Organization for Research and Treatment of Cancer (EORTC) Elderly Task Force [7]. Following the rising life expectancy in many countries, it is of importance to accurately predict cancer burden among aging cancer populations suffering from concomitant disease at time of cancer diagnosis.

We examined the association of this biomarker-based score with mortality and with respect to CCI and other potential modifying factors among 13,056 participants in the Third National Health and Nutrition Examination Survey (NHANES III), a large study designed to be nationally representative of the United States population.

Methods

Study population and data collection

This study is based on the NHANES III Mortality linkage which provides follow-up data from the date of NHANES III survey participation (1988–1994) through December 31, 2006.

NHANES III is a cross-sectional study conducted by the National Center for Health Statistics (NCHS) between 1988 and 1994 [8]. It was designed as a multistage stratified, clustered probability sample of the US civilian non-institutionalized population at least two months old. Subjects participated in an interview conducted at home and an extensive physical examination, which included a blood sample, conducted in a mobile examination center [8]. NHANES III was conducted in two phases (1988–1991 and 1991–1994). Unbiased national estimates of health and nutrition characteristics can be independently produced for each phase. Within each phase, subjects were randomly assigned to participate in either the morning or afternoon/evening examination session. In total, 30,818 people were interviewed in NHANES III and had a physical examination and a blood sample taken. The NCHS has updated the mortality linkage of NHANES III to death certificate data found in the National Death Index (NDI) so that mortality ascertainment is based upon the results from a probabilistic match between NHANES III and NDI death certificate records [9]. From this linkage, we selected all men and women eligible for mortality follow-up (all participants 17 years of older) with baseline measurements of CRP, GGT, HDL, and albumin (n=13,056).

Serum CRP was quantified using latex-enhanced nephelometry, standardized according to the World Health Organization's international reference preparation of CRP with the coefficient of variation ranging from 3.2–16.1% (median, 6.3%) through the study period [10]. Serum GGT concentration was assayed with a Hitachi 737 Analyzer (Boehringer–Mannheim Diagnostics, Indianapolis, IN, USA) with coefficients of variation ranging from 1.9 to 3.5% [10]. HDL cholesterol was measured using standard enzymatic methods and serum albumin was assessed with the bromcresol purple method. Details about all laboratory procedures have been published elsewhere [10].

Other covariates in the analysis included age, race/ethnicity, poverty-to-income-ratio, comorbidity index, smoking, alcohol consumption, and physical activity. Race and ethnicity were combined into four race–ethnic groups: Non-Hispanic white, Non-Hispanic black, Mexican American and Other. The poverty to income ratio is an index of poverty status which is calculated by dividing family income by a poverty threshold specific to family size.

Comorbidity was evaluated with a comorbidity coefficient similar to the CCI, as used in other NHANES III-based analyses [11]. Each of the comorbidities available in the dataset contributed one point to the composite index with additional points given for older age. Cigarette smoking, alcohol consumption, and physical activity were assessed using questionnaires. The protocols for the conduct of NHANES III were approved by the institutional review board of the National Center for Health Statistics, Centers for Disease Control and Prevention. Written informed consent was obtained from all participants.

Data analysis

To address the combined association between the four biomarkers of interest and all-cause mortality we started with a multivariate Cox proportional hazards regression using continuous variables of CRP (mg/L), albumin (g/L), GGT (U/L), and HDL (mmol/L). A mortality score, ranging from 0 to 4, was then calculated as the number of biomarkers with abnormal values according to their clinical cut-offs (CRP>10mg/L, albumin<35g/L, GGT>35U/L, and HDL<1.03mmol/L) [3, 12, 13]. Multivariate Cox proportional hazards regression was used to investigate this mortality score in relation to all-cause mortality as well as cancer-specific (ICD10:C00-C99) and circulatory-specific death (ICD10: I00-I99). Sampling weights for NHANES III were used to account for sampling variability and to adjust for differential probability of selection of persons [8]. All models were adjusted for age (continuous), race/ethnicity, poverty to income ratio, comorbidity index, smoking, alcohol consumption, physical activity, and BMI. A test for trend was conducted by using assignment to quartiles as an ordinal scale, meaning a continuous variable ranging from 1 to 4 indicating the different quartiles. Stratified analyses were conducted for men and women [14], age-groups (<65 and 65+), different categories of the CCI, and ethnic groups. The latter offers the opportunity to evaluate whether the score holds across different ethnic groups and is in concordance with many other studies performed in NHANES [15–17]. We also assessed the association between the mortality score and the CCI by calculating the correlation coefficient and kappa’s coefficient of agreement between both measurements. All tests were two-sided; p-values <0.05 were considered to be statistically significant. All analyses were conducted with Statistical Analysis Systems (SAS) release 9.1.3 (SAS Institute, Cary, NC).

Results

A total of 3,081 persons died during follow-up time, of whom 653 (21.19 %) died of cancer and 1,388 (45.05%) of circulatory disease. Population characteristics are shown in Table 1. A higher CCI was observed for those who died during follow-up than for those who were alive at the end of follow-up (32.91% versus 8.11% with CCI=4+) (Table 1).

Table 1 Descriptive statistics of study population by vital status

A multivariate Cox proportional hazards model including continuous variables of CRP (mg/L), albumin (g/L), GGT (U/L), and HDL (mmol/L) showed that CRP and GGT were statistically significantly associated with all-cause mortality: HR per one unit increase: 1.19 (95%CI: 1.11, 1.27), 0.95 (0.93, 0.96), 1.31 (1.25, 1.39), and 0.98 (95%CI: 0.91, 1.05), respectively. A log transformation was performed for CRP, GGT, and HDL due to their skewed distributions.

When using the mortality score, there was a statistically significant positive association between the score and all-cause mortality as well as cancer and circulatory disease-specific death (e.g. HR for all-cause mortality: 1.21 (95%CI: 1.09, 1.35), 1.92 (95%CI: 1.67, 2.20), 3.38 (95%CI: 2.62, 4.36), and 7.93 (95%CI: 5.77, 10.89) for score=1, 2, 3 and 4 compared to score=0). Similar patterns were found when stratifying by sex and age groups (Table 2 and 3). Figure 1 also shows how each marker of the mortality score contributes to the risk of death. There is no clear pattern by specific biomarkers, but GGT and albumin seem to contribute more than CRP and HDL cholesterol.

Table 2 Hazard Ratio (HR) and 95% Confidence Intervals (CI) for risk of all-cause, cancer-specific, and circulatory disease death

Table 3 Age-group specific analysis: Hazard Ratio (HR) and 95% Confidence Intervals (CI) for risk of all cause, cancer-specific, and circulatory disease death

Figure 1

Hazard ratio for risk of all-cause death by different values of the mortality score. The model was adjusted for age, gender, poverty to income ratio, race/ethnicity, smoking behavior, alcohol consumption, vigorous physical activity, Charlson Comorbidity Index, and BMI.

To compare the mortality score with the more commonly used CCI, we estimated the association between the CCI and all-cause mortality, which showed analogous risks (HR: 1.06 (95%CI: 0.45, 2.49), 1.55 (95%CI: 0.65, 3.71), and 2.41 (95%CI: 1.00, 5.76) for CCI=2, 3, and 4+ compared to CCI=0; P for trend<0.001). Effect modification by CCI for the association between the mortality score and risk of dying was then assessed with stratified analyses by categories of CCI (Table 4). The patterns observed in Table 2 were seen in each category of CCI, even among those with a CCI < 3. For example, among those with CCI =3, the risk of cancer-specific death was 1.10 (95%CI: 0.93, 1.28), 1.81 (95%CI: 1.43, 2.29), 4.67 (95%CI: 3.05, 7.14), and 6.97 (95%CI: 5.32, 9.14) for score=1, 2, 3 and 4 compared to score=0 (Table 4). The univariate association between the mortality score and the CCI did not show a strong correlation between both measurements (correlation coefficient: 0.15; P<0.001 and kappa’s coefficient of agreement: 0.01; P<0.001).

Table 4 Hazard Ratio (HR) and 95% Confidence Intervals (CI) for risk of all-cause, cancer-specific, and circulatory death, stratified by Charlson Comorbidity index

A stratified analysis by race/ethnicity showed that the score predicted mortality in a similar way for non-Hispanic white, non-Hispanic blacks, and Mexican American (Table 5). The results for cancer and circulatory death were also comparable with the findings in Table 2.

Table 5 Hazard Ratio (HR) and 95% Confidence Intervals (CI) for risk of all-cause, cancer-specific, and circulatory death, stratified by race/ethnicity

Excluding those with follow-up of <1 year, as part of a sensitivity analysis, showed similar patterns to those observed in Table 2 (results not shown). For example, the HRs for all-cause mortality increased with values of the mortality score: 1.20 (95%CI: 1. 07, 1.34), 1.81 (95%CI: 1.58, 2.08), 2.99 (95%CI: 2.33, 3.84), and 5.55 (95%CI: 4.32, 7.14), for score=1, 2, 3 and 4 compared to score=0. Similar observations were made when excluding those with follow-up <3 or <5 years (results not shown).

Discussion

Using NHANES-III data we could replicate the results of the Swedish AMORIS study [1]. The combination of serum levels of CRP, albumin, GGT, and HDL was positively associated with mortality and was even predictive in patients with low comorbidity based on CCI. Furthermore, the association between the score and mortality was found to be similar among non-Hispanic whites, non-Hispanic blacks, and Mexican Americans.

The three suggested mechanisms, inflammation, liver dysfunction, and lipid metabolism, have been associated with early death and comorbidities in several different research settings. The inflammation-based Glasgow Prognostic Score (GPS), based on serum levels of CRP and albumin, has repeatedly been shown to be a predictor of survival, independent of tumor stage, performance status and treatment [3, 18–20]. In a study of 540 cancer patients, increasing GPS correlated with more aggressive tumor biology in terms of tumor size, presence of lymph node metastasis, and higher tumor recurrence rate [20]. In addition to being a marker of liver dysfunction, high levels of GGT (>28 U/L) have been found to be positively associated with risk of all cause mortality and risk of developing cancer [21, 22]. It is also thought that persistent production of reactive oxygen stress following increased GGT expression in tumor cells contributes to genetic instability and tumor progression [19]. Finally, HDL cholesterol, as a component of the metabolic syndrome, has been extensively studied in relation to early death and comorbidities [12, 23]. Apart from its link with the lipid metabolism, HDL is also associated with inflammation [6, 23, 24]. More specifically, HDL has been linked with pro-inflammatory cytokines such as tumor necrosis factor-α (TNF-α) [25]. It has been shown that HDL reduces free TNF-α resulting in reduced tissue damage, reduced infiltration of macrophages and neutrophils, and potential attenuated tumor formation [26]. A prospective study of 989 persons aged 65 and over, showed that high levels of HDL were significantly associated with low total risk of death in men [27].

As in the Swedish AMORIS study, the HRs found using data from NHANES III were also of a magnitude that may be of clinical relevance. Even though Figure 1 did not show a clear pattern by specific markers of the mortality score, albumin and GGT seemed to contribute slightly more to the prediction of death than the other markers. It is important to note that these results are slightly different from our analysis based on continuous levels of the four studied biomarkers, which is likely explained by the more medically relevant cut-offs used in our mortality score. The findings in Figure 1 suggest that markers of inflammation and reactive oxygen species are more strongly associated with risk of death. Because of limited sample size, it was not possible to perform identical age-stratified analysis so that we could not identify the predictive strength of the score in those aged 50–65. In our analysis, the score was found to be equally predictive of death among those younger and older than 65. Interestingly, we also replicated the findings stratified by categories of the CCI. Among those with CCI<3, there was a positive trend between the mortality score and all-cause mortality as well as circulatory disease-specific death. The association was not statistically significant with cancer-specific death; this may be due to the comparably small number of cancer deaths (24 versus 42%). Nevertheless, the results by categories of CCI indicate that the mortality score is predictive for mortality over and above the prediction by CCI.

In our study we could also evaluate the mortality score in different ethnic groups. A study evaluating the clinical utility of the triglycerides-to-HDL ratio showed that this ratio may identify insulin-resistance in Aboriginals, Chinese, and Europeans, but not South Asians, thus reflecting possible ethnic differences in normal lipid levels [28]. Also, levels of albumin also differed by ethnicity. Among UK resident patients of white European or south Asian ethnicity with type 2 diabetes mellitus there were significant differences in microalbuminuria [29]. In addition, the association between CRP and disease outcome appeared to be influenced by ethnicity. In a multi-ethnic cohort including 2362 Caucasians, 1601 African Americans, 1353 Hispanics, and 751 Chinese, CRP was only associated with the risk of circulatory disease in Caucasians [30]. This variation by ethnicity was smaller for GGT. In another study based on NHANES III, the correlation between GGT and risk of hypertension was consistent by race-ethnicity [31]. In our study there was also no effect-modification by race-ethnicity; the mortality score seemed to predict all cause, cancer-specific and circulatory death in a consistent way across ethnic groups. This may suggest that a combination of these four biomarkers is more general and less susceptible to subgroup differences.

This study has several strengths including its generalizibility following the use of nationally representative data. Therefore it was also possible in our analysis to perform a stratified analysis by race/ethnicity. We were able to adjust for many potential confounding factors and examine interactions by age and the CCI. In a sensitivity analysis, we used different medical cut-offs by sex for HDL and GGT, but this did not materially alter our findings (results not shown). Another limitation of this study and the mortality score is that it relies on one single measurement so that it may be prone to measurement error and within-person variation. Repeated measurements may strengthen the accuracy of the score. In addition, it is important to note that the current mortality score only represents a selection of mechanisms and markers associated with early death. It is possible that other markers such as hypertension or the ratio of total to HDL cholesterol may also be good components of a mortality score, however we believe that the current mortality score already represents three important mechanisms associated with early death. A next step in the process of potentially taking this score in clinic is the performance of sensitivity and specific tests in order to determine the most predictive cut-offs for the different biomarkers used.

Conclusion

The combination of four biomarkers predicted all-cause mortality as well as cancer and circulatory disease-specific death in NHANES III. This is in agreement with a Swedish study and was consistent across different ethnic groups. This trans-ethnic validity and the positive association for all categories of the CCI suggests that this score may help clinicians making decisions in prevention and disease management.