Multivariable-adjusted relations of serum phosphorus levels to the incidence of cardiovascular disease (CVD) on follow-up. Estimated multivariable hazard ratios for CVD (y-axis) in relation to serum phosphorus levels (x-axis) as a function of penalized regression splines are shown. Vertical lines indicate quartiles of phosphorus levels. The dotted lines indicate confidence intervals. To convert serum phosphorus levels to millimoles per liter, multiply by 0.3229.
Dhingra R, Sullivan LM, Fox CS, Wang TJ, D’Agostino RB, Gaziano JM, Vasan RS. Relations of Serum Phosphorus and Calcium Levels to the Incidence of Cardiovascular Disease in the Community. Arch Intern Med. 2007;167(9):879-885. doi:10.1001/archinte.167.9.879
Copyright 2007 American Medical Association. All Rights Reserved. Applicable FARS/DFARS Restrictions Apply to Government Use.2007
Higher levels of serum phosphorus and the calcium-phosphorus product are associated with increased mortality from cardiovascular disease (CVD) in patients with chronic kidney disease (CKD) or prior CVD. However, it is unknown if serum phosphorus levels influence vascular risk in individuals without CKD or CVD.
We prospectively evaluated 3368 Framingham Offspring study participants (mean age, 44 years; 51% were women) free of CVD and CKD. We used multivariable Cox models to relate serum phosphorus and calcium levels to CVD incidence.
On follow-up (mean duration, 16.1 years), there were 524 incident CVD events (159 in women). In multivariable analyses and adjusting for established risk factors and additionally for glomerular filtration rate and for hemoglobin, serum albumin, proteinuria, and C-reactive protein levels, a higher level of serum phosphorus was associated with an increased CVD risk in a continuous fashion (adjusted hazard ratio per increment of milligrams per deciliter, 1.31; 95% confidence interval, 1.05-1.63; P = .02; P value for trend across quartiles = .004). Individuals in the highest serum phosphorus quartile experienced a multivariable-adjusted 1.55-fold CVD risk (95% confidence interval, 1.16%-2.07%; P = .004) compared with those in the lowest quartile. These findings remained robust in time-dependent models that updated CVD risk factors every 4 years and in analyses restricted to individuals without proteinuria and an estimated glomerular filtration rate greater than 90 mL/min per 1.73 m2. Serum calcium was not related to CVD risk.
Higher serum phosphorus levels are associated with an increased CVD risk in individuals free of CKD and CVD in the community. These observations emphasize the need for additional research to elucidate the potential link between phosphorus homeostasis and vascular risk.
Serum phosphorus levels are tightly regulated within a narrow range in healthy individuals because inorganic phosphorus plays a fundamental physiological role in energy production, membrane transport, and signal transduction.1 Recent experimental and clinical data both implicate higher serum phosphorus levels in the pathogenesis of vascular disease. In murine models, deletion of the fibroblast growth factor–23 gene2 (that regulates renal and intestinal phosphorus absorption) or the Klotho gene3 (that encodes a transmembrane protein in the renal distal convoluted tubule) are associated with a phenotype that is characterized by hyperphosphatemia, premature aging, arteriosclerosis, and vascular calcification. In addition, the phenotype in the Klotho null mice is reversed by a low-phosphorus diet.4 On a parallel note, several clinical studies have shown that higher serum phosphorus levels and the calcium-phosphorus product in patients with CKD are associated with increased carotid atherosclerosis and coronary calcification cross-sectionally5- 8 and with an elevated risk of all-cause mortality and cardiovascular disease (CVD) mortality longitudinally.9- 11 Recently, Tonelli et al12 reported an association of higher serum phosphorus levels within the reference range with an increased risk of CVD events during a 5-year follow-up period in patients with a previous myocardial infarction who have dyslipidemia but normal kidney function. Of note, polymorphisms of the Klotho gene have been associated with an adverse risk factor profile,13 and with atherosclerotic coronary disease in some reports14 but not in others.15
The aforementioned investigations12- 14 raise the possibility that serum phosphorus levels within the reference range may influence the risk of CVD in asymptomatic individuals free of prevalent CKD or CVD in the general population. Accordingly, we investigated the clinical correlates of serum phosphorus levels and related levels of serum phosphorus and calcium to the incidence of CVD in a large, community-based sample of individuals who were free of CVD and CKD at baseline.
The design and selection criteria of the Framingham Offspring Study have been described previously.16 All participants undergo routine physical examinations at the Framingham Heart Study approximately every 4 years. For the study described herein, we evaluated participants who had levels of serum phosphorus and calcium measured routinely at the second examination cycle (1979-1982) (n = 3676). We excluded participants with prevalent CVD (n = 167), defined as any history of coronary heart disease, cerebrovascular disease, peripheral vascular disease, or heart failure.17 The Modification of Diet in Renal Disease equation was used to calculate an estimated glomerular filtration rate (eGFR): 186.3 × (serum creatinine)−1.154 × age−0.203 × (0.742 for women).18 Individuals with an eGFR of less than 60 mL/min per 1.73 m2 at baseline were considered to have impaired kidney function18 and were excluded (n = 141). After exclusions, a total of 3368 individuals (51% of whom were women) were eligible for the present investigation. All participants provided written informed consent, and the study protocol was approved by the institutional review board of the Boston Medical Center, Boston, Mass.
At each examination visit, medical histories were taken, and participants underwent a physical examination (including blood pressure measurements), anthropometry, and laboratory assessment of risk factors. Hypertension was defined as a systolic blood pressure of 140 mm Hg or higher, diastolic blood pressure of 90 mm Hg or higher, or use of antihypertensive medications.19 Smoking 1 or more cigarettes daily on a regular basis in the year prior to the examination visit was considered current smoking. Total intake of alcohol was assessed by averaging the self-reported weekly consumption of alcoholic drinks. Fasting levels of total cholesterol, high-density lipoprotein cholesterol, serum albumin, blood glucose, and hemoglobin were measured using standardized assays and serum creatinine levels by using the modified Jaffe method. To be consistent with other laboratories, serum creatinine values were calibrated to the Cleveland Clinic Laboratory standard using a correction factor of 0.23 mg/dL (20.33 μmol/L).20 Diabetes mellitus was defined as a fasting blood glucose level of 126 mg/dL (7 mmol/L) or higher or the use of insulin or any hypoglycemic agents. Proteinuria was measured by a urine dipstick test (Ames Labstix; Ames Co Inc, Elkhart, Ind) and coded as none, trace, 1+, or 2+ or more. For analytical purpose, we dichotomized proteinuria as “none” vs “trace or more.” High-sensitivity C-reactive protein (hsCRP) was measured with a nephelometer (model BN100; Dade Behring; Deerfield, Ill) with an average intra-assay coefficient of variation of 3.8%.
Serum phosphorus and calcium levels were measured using a standard colorimetric method (Roche Diagnostics, Alameda, Calif). Serum phosphorus had an intra-assay coefficient of variation of 5.6%, and serum calcium had an intra-assay coefficient of variation of 2.5%.
All participants in the Framingham Heart Study are under continuous surveillance for CVD events and death. A committee of 3 experienced investigators (who were blinded to serum phosphorus and calcium levels) reviewed all suspected CVD events by examining hospitalization records, physician office visit notes, and pathology reports, if appropriate. A Framingham Heart Study neurologist evaluated participants with suspected cerebrovascular events, and a separate review committee that included a neurologist adjudicated these events. Criteria for CVD have been described previously.17 For the present study, incident CVD was defined as fatal or nonfatal myocardial infarction, angina pectoris (stable or unstable), cerebrovascular events (stroke or transient ischemic attacks), peripheral vascular disease, or congestive heart failure. This composite definition of CVD has been used consistently by Framingham investigators consistently over the years.17
To maximize our statistical power, we evaluated first-onset CVD (and not individual CVD outcomes separately). Primary analyses focused on serum phosphorus and calcium and the calcium-phosphorus product, each considered separately. Baseline characteristics of study participants were compared across quartiles of serum phosphorus using analysis of variance. Multivariable regression analysis was used to test for trends in clinical characteristics across phosphorus quartiles adjusted for age and sex. We evaluated the cross-sectional correlates of the continuous measure of serum phosphorus using sex-pooled multivariable linear regression analyses.
We calculated age- and sex-adjusted incidence of CVD on follow-up according to quartiles of serum phosphorus and calcium and the calcium-phosphorus product. Cox proportional hazard regression models21 were estimated to analyze the relations of serum phosphorus and calcium levels, and the calcium-phosphorus product to incident CVD events after confirming the assumption of proportionality of hazards for each variable.22 Serum phosphorus and calcium levels and the calcium-phosphorus product were analyzed as continuous variables and as quartiles. Because tests of interaction for sex with serum phosphorus and calcium levels and the calcium-phosphorus product were not statistically significant, all analyses were performed for pooled sexes to maximize statistical power.
Two sets of multivariable models were constructed to begin with: (1) initially adjusting for age and sex and (2) adjusting for other covariates (including established vascular risk factors), age, body mass index (BMI) (calculated as weight in kilograms divided by height in meters squared), diabetes mellitus, systolic blood pressure, treatment for hypertension, smoking, alcohol consumption, total/high-density lipoprotein (HDL) cholesterol ratio, levels of hemoglobin and serum albumin, eGFR, proteinuria, and hsCRP (multivariable model 1). Proteinuria and eGFR were used as covariates to account for subclinical renal disease. We adjusted for hsCRP because of its association with CVD incidence.23 We also adjusted for the total/HDL cholesterol ratio instead of low-density lipoprotein cholesterol levels because Framingham data indicate that the total/HDL cholesterol ratio is a better predictor of coronary disease risk relative to low-density lipoprotein cholesterol.24 In models considering the quartiles, we compared hazard ratios in the second, third, and fourth quartiles with those in the first quartile (referent) and tested for a linear trend in risk across quartiles. We also tested for interactions between serum phosphorus and calcium levels and the calcium-phosphorus product (modeled as continuous variables and using separate analyses for each) and age, hypertension status, and eGFR in the multivariable models.
In addition, the design of the Framingham Offspring Study includes assessment of CVD risk factors (BMI, diabetes mellitus, hypertension, smoking, alcohol consumption, current use of antihypertensive drug treatment, and total/HDL cholesterol ratio) at each quadrennial examination. We repeated the analyses described at the beginning of this subsection with time-dependent Cox models, updating established CVD risk factors every 4 years (multivariable model 2); variables that were not ascertained every 4 years (ie, serum phosphorus and calcium, serum albumin, and hemoglobin levels; eGFR; proteinuria; and hsCRP) were not updated. Because assessment of kidney disease with the Modification of Diet in Renal Disease equation may misclassify some people with CKD,25 we repeated the time-dependent multivariable analyses after excluding individuals with an eGFR of 90 mL/min per 1.73 m2 or less and participants with trace or more proteinuria (multivariable model 3).
To gain additional insights into the linearity of relations between serum phosphorus levels and incident CVD, we estimated generalized additive Cox models using penalized splines (covariates as in multivariable model 2).26 These models are useful for assessing the dose-response relations of serum phosphorus levels and the risk of CVD.27
In our sample, there were 9 participants with serum phosphorus levels exceeding 4.5 mg/dL (1.45 mmol/L), the upper limits of the reference range level for serum phosphorus level as described in the literature.1 Because we were interested in the relations of serum phosphorus levels within the reference range to CVD risk, analyses were repeated excluding these individuals.
In primary analyses, we used serum calcium levels and adjusted for serum albumin values. We performed secondary analyses using serum calcium levels corrected for serum albumin values using the following formula: corrected serum calcium (in milligrams per deciliter) = observed serum calcium + [0.8 × (4 − serum albumin)], if the serum albumin level is less than 4 g/dL.
Although the primary objective of our study was not to estimate CVD risk prediction models, we assessed if the addition of serum phosphorus levels improved the discrimination (assessed by C statistic) and/or the calibration (assessed by Hosmer-Lemeshow χ2 statistic) of a multivariable model that incorporated standard risk factors.
All analyses were performed with SAS statistical software (version 9.1; SAS Inc, Cary, NC), and the S-plus program (Insightful Corp, Basel, Switzerland) was used to plot the regression splines. A 2-sided P value of P<.05 was considered statistically significant.
Serum phosphorus and calcium levels and calcium-phosphorus product were normally distributed in both sexes. Quartile ranges and mean values of serum phosphorus and calcium were similar among both men and women. The baseline characteristics of participants are displayed in Table 1 accordingly to quartiles of serum phosphorus levels. The prevalence of several established CVD risk factors (except smoking) decreased with increasing levels of phosphorus (Table 1). In a multivariable linear regression model, serum phosphorus values were positively associated with age, hsCRP, eGFR, serum albumin level, and total/HDL cholesterol ratio but were inversely related to BMI, systolic blood pressure, and hemoglobin level (P<.02 for all). Multivariable-adjusted levels of serum phosphorus were lower in men compared with women (P<.001).
Serum phosphorus levels were weakly and positively correlated with serum calcium levels (r = 0.12 in men and 0.16 in women; P<.001 for both). The calcium-phosphorus product was highly correlated with serum phosphorus (r = 0.96 in both sexes; P<.001) but only modestly so with serum calcium levels (r = 0.39 in men and 0.42 in women; P<.001 for both). Consequently, the distribution of individuals in the quartiles of the calcium-phosphorus product mirrored their distribution in serum phosphorus level quartiles (data not shown).
During a follow-up period of 20 years (mean, 16.1 years), 524 participants experienced a first CVD event (159 in women); these included 138 myocardial infarctions, 173 angina pectoris events, 93 cerebrovascular events, 18 sudden cardiac deaths, 63 peripheral vascular disease, and 39 congestive heart failure events. The age- and sex-adjusted incidence of CVD rose across quartiles of serum phosphorus levels (and quartiles of calcium-phosphorus product; data not shown) but not across serum calcium level quartiles (Table 2).
In age- and sex-adjusted models, the hazard ratio for incident CVD relative to a 1-mg/dL (0.3229-mmol/L) increase in serum phosphorus levels was 1.17 (95% confidence interval [CI], 0.96-1.44; P = .12). However, after adjusting for all other covariates, there was a statistically significant increase in the risk of CVD with increasing serum phosphorus levels (hazard ratio for a 1-mg/dL (0.3229-mmol/L) increase of 1.31; 95% CI, 1.05-1.63; P = .02). In multivariable models, individuals in the highest quartile were 1.55 times more likely to experience CVD than those in first quartile (Table 3, multivariable model 1). We observed a statistically significant linear trend for increasing risk of CVD across quartiles of serum phosphorus levels (Table 3) that were supported by regression splines (Figure). The graded association of serum phosphorus with CVD risk remained robust in multivariable models updating established CVD risk factors as time-varying covariates every 4 years (Table 3, multivariable model 2). These results remained robust in analyses restricted to individuals with an eGFR greater than 90 mL/min per 1.73 m2 and without proteinuria (Table 3, multivariable model 3). We did not observe effect modification by age, sex, hypertension status, serum calcium level, or eGFR (interaction terms were not statistically significant).
The analyses of the calcium-phosphorus product to the risk of CVD yielded results similar to those for serum phosphorus levels. The calcium-phosphorus product also was positively related to CVD risk (hazard ratio = 1.12 increase per SD increment; 95% CI, 1.02-1.24; P = .02). There was also a significant trend of increasing CVD risk across quartiles of calcium-phosphorus product (P = .01). These relations were maintained in time-dependent Cox models.
Serum calcium levels were not related to increased CVD events in age-adjusted or in multivariable models (Table 3; a 1 mg/dL (0.25 mmol/L) increment was associated with a hazards ratio for CVD of 1.10; 95% CI, 0.81-1.49; P = .56).
The association of serum phosphorus levels with increased CVD risk remained robust when analyses were repeated excluding the 9 participants with values exceeding 4.5 mg/dL (1.45 mmol/L). We repeated analyses using corrected serum calcium levels (adjusted for serum albumin levels), and observed results identical to models using uncorrected serum calcium levels.
The C statistic for the multivariable model incorporating standard CVD risk factors was 0.79 (95% CI, 0.77-0.81) with good calibration (χ² = 14.1). On addition of serum phosphorus values (modeled as quartiles), the C statistic increased minimally to 0.80 (95% CI, 0.78-0.82), but the calibration improved modestly (χ² = 6.5).
Our principal findings are 3-fold. First, cross-sectionally, serum phosphorus levels demonstrated an interesting relation to CVD risk factors, being positively related to age and total/HDL cholesterol ratio but inversely related to BMI and systolic blood pressure. Other investigators have reported similar cross-sectional correlates of serum phosphorus levels in patients with CKD28 and in people without CKD.29 Second, serum phosphorus levels (and the highly correlated calcium-phosphorus product) were associated with an increased risk of incident CVD in a continuous fashion in multivariable models adjusting for hsCRP and in a model incorporating established CVD risk factors as time-varying covariates. Values of serum phosphorus higher than 3.5 mg/dL (.382 mmol/L) (defining the lower limit of the highest quartile and well within the reference range) were associated with a 55% increased CVD risk. Furthermore, these associations remained robust in a subgroup of individuals with an eGFR greater than 90 mL/min per 1.73 m2 and without proteinuria. Of note, relations of serum phosphorus levels to CVD incidence strengthened in multivariable models compared with age- and sex-adjusted models, possibly because of confounding of the latter models by the inverse relations of serum phosphorus levels to select CVD risk factors. To our knowledge, the present investigation is the first to demonstrate a graded independent relation of serum phosphorus levels (within the reference range) to CVD risk in a community-based sample of men and women without either CKD or CVD at baseline. Third, serum calcium levels were not related to incident CVD in any of the models. Given the very high correlation of serum phosphorus levels and the calcium-phosphorus product, and the lack of any association of serum calcium with CVD risk, the association of the calcium-phosphorus product with CVD risk was likely driven by the association of serum phosphorus levels with CVD risk.
There are several potential mechanisms that may explain the association of serum phosphorus levels with greater CVD risk. First, high phosphorus levels are known to inhibit 1,25-dihydroxyvitamin D synthesis in some studies.30,31 Lower levels of 1,25-dihydroxyvitamin D are hypothesized to decrease cardiac contractility32 and to increase coronary calcification.33 Additional research studies may be required to evaluate if lower 1,25-dihydroxyvitamin D levels in individuals with higher serum phosphorus levels contribute to a greater vascular risk.
Second, higher serum phosphorus levels may directly promote vascular injury. Investigators have reported that higher phosphorus levels increase the propensity of mineral deposition in vascular smooth muscle cells in vitro.34 This increased predisposition can be partly explained by increased osteopontin expression35,36 Serum phosphorus may also directly increase vascular calcification, especially when levels of calcium-phosphorus product are high, as observed in patients with CKD8 or in individuals without CKD (referred to as dystrophic calcification).37 Most observational studies of patients with CKD have noted significant associations of vascular calcification with increased levels of serum phosphorus and calcium-phosphorus product. It is not clear, however, if higher serum phosphorus levels within the reference range are associated with increased vascular calcification.
Third, higher serum phosphorus levels may indicate subclinical renal dysfunction, with associated cardiovascular sequelae. Although our data lacked information on parathyroid hormone (PTH), we performed additional analyses excluding individuals with an eGFR of 90 mL/min per 1.73 m2 or less and/or proteinuria; the relations of serum phosphorus levels to CVD risk remained robust in these analyses. It is also important to note that none of the participants in our sample had a calcium-phosphorus product level of more than 55 mg2/dL2 (a threshold above which adverse events have been reported to escalate in patients with CKD9). Therefore, we do not believe that secondary hyperparathyroidism caused by renal dysfunction is a likely mechanism that would explain the observed relations in our study.
Fourth, higher serum phosphorus levels increase circulating PTH levels even in healthy individuals.38 Higher PTH levels may be proinflammatory; PTH has been reported to induce production of IL-6 (interleukin 6) via bone resorption39 or by increased hepatic synthesis.40 Higher levels of serum IL-6 and hsCRP have been associated with increased CVD risk. In our study, we adjusted for hsCRP and still observed a significant independent association of serum phosphorus levels and CVD risk
Surveys conducted by the US Department of Agriculture demonstrate an increase in the dietary intake of phosphorus among Americans, accompanied by a gradual decline in intake of calcium over the past 2 decades.41 Although the significance of increased consumption of phosphorus is unclear, limited data suggest that greater dietary intake of phosphorus may be associated with higher serum phosphorus levels.30,38,42,43
The Institute of Medicine has formulated dietary reference guidelines for phosphorus intake for Americans that suggest higher phosphorus levels may not be associated with adverse health effects in the general population.44 These dietary references have been questioned because the rising dietary phosphorus intake is paralleled by a greater risk of fractures and lower bone density among children.45,46 In the present study, we did not assess dietary phosphorus intake. Further research is warranted, therefore, to evaluate if a higher dietary intake of phosphorus translates into a greater risk for CVD in the general population by contributing to higher serum phosphorus levels.
The large, community-based sample of men and women; exclusion of individuals with CKD at baseline; adjustment for eGFR, proteinuria, and hsCRP; and the extended follow-up of 20 years strengthen our study.
Several limitations of our study deserve comment. We did not measure circulating PTH or 1,25-dihydroxyvitamin D levels in our sample to gain insights into mechanisms by which serum phosphorus levels may increase vascular risk. We attempted to control for secondary hyperparathyroidism by excluding individuals with CKD at baseline and by adjusting for eGFR and proteinuria in multivariable models. We also did not adjust for physical activity or other dietary determinants of CVD. In addition, we performed sex-pooled analyses and analyzed the composite outcome of first CVD events; we had limited power to evaluate the risk of individual CVD outcomes. Furthermore, we did not adjust for multiple statistical testing; however, all analyses were defined a priori. Finally, all of our participants were white, which may limit the generalizability of our results to nonwhite individuals.
In conclusion, in a large community-based sample of individuals free of baseline CKD and CVD, higher levels of serum phosphorus are associated with an increased risk of incident CVD in a graded fashion. Additional research is warranted to confirm or refute our findings and to elucidate the mechanisms underlying the potential link between phosphorus homeostasis and vascular risk.
Correspondence: Ramachandran S. Vasan, MD, Framingham Heart Study, 73 Mount Wayte Ave, Suite 2, Framingham, MA 01702-5803 (firstname.lastname@example.org).
Accepted for Publication: November 28, 2006.
Author Contributions: Dr Vasan had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: Dhingra and Vasan. Acquisition of data: Dhingra, Sullivan, and Vasan. Analysis and interpretation of data: Dhingra, Sullivan, Fox, Wang, Gaziano, and Vasan. Drafting of the manuscript: Dhingra, Sullivan, Fox, and Vasan. Critical revision of the manuscript for important intellectual content: Dhingra, Sullivan, Fox, Wang, D’Agostino, Gaziano, and Vasan. Statistical expertise: Sullivan, D'Agostino, and Vasan. Administrative, technical, or material support: D'Agostino and Vasan. Study supervision: Dhingra, Fox, Sullivan, D'Agostino, Gaziano, and Vasan.
Financial Disclosure: Dr Gaziano has received funding from McNeil Consumer Products Co and Pliva; has received research support in the form of pills or packaging from BASF Group, DMS Pharmaceuticals Group Inc, and Wyeth Pharmaceuticals; has received honoria from Bayer and Pfizer Pharmaceutical Co for speaking engagements; has served as a consultant for McNeil Consumer Products Co and Wyeth Pharmaceuticals; and has served as an expert witness for Merck & Co Inc, NutraQuest Inc, and GlaxoSmithKline.
Funding/Support: This work was supported through the National Institute of Health/National Heart, Lung, and Blood Institute (NHLBI) contracts N01-HC-25195, 1R01HL67288, and 2K24HL04334 (Dr Vasan) and K23HL74077 (Dr Wang). Dr Gaziano has received investigator-initiated federal funding from the National Institutes of Health (National Cancer Institute, NHLBI, National Institute on Aging, National Eye Institute), and the Veterans Administration (Cooperative Study Program).