Genetic Evidence for a Causal Role of Serum Phosphate in Coronary Artery Calcification: The Rotterdam Study.

Background Hyperphosphatemia has been associated with coronary artery calcification (CAC) mostly in chronic kidney disease, but the association between phosphate levels within the normal phosphate range and CAC is unclear. Our objectives were to evaluate associations between phosphate levels and CAC among men and women from the general population and assess causality through Mendelian randomization. Methods and Results CAC, measured by electron-beam computed tomography, and serum phosphate levels were assessed in 1889 individuals from the RS (Rotterdam Study). Phenotypic associations were tested through linear models adjusted for age, body mass index, blood pressure, smoking, prevalent cardiovascular disease and diabetes, 25-hydroxyvitamin D, total calcium, C-reactive protein, glucose, and total cholesterol : high-density lipoprotein cholesterol ratio. Mendelian randomization was implemented through an allele score including 8 phosphate-related single-nucleotide polymorphisms. In phenotypic analyses, serum phosphate (per 1 SD) was associated with CAC with evidence for sex interaction (Pinteraction=0.003) (men β, 0.44 [95% CI, 0.30-0.59]; P=3×10-9; n=878; women β, 0.24 [95% CI, 0.08-0.40]; P=0.003; n=1011). Exclusion of hyperphosphatemia, chronic kidney disease (estimated glomerular filtration rate <60 mL/min per 1.73 m2) and prevalent cardiovascular disease yielded similar results. In Mendelian randomization analyses, instrumented phosphate was associated with CAC (total population β, 0.93 [95% CI: 0.07-1.79]; P=0.034; n=1693), even after exclusion of hyperphosphatemia, chronic kidney disease and prevalent cardiovascular disease (total population β, 1.23 [95% CI, 0.17-2.28]; P=0.023; n=1224). Conclusions Serum phosphate was associated with CAC in the general population with stronger effects in men. Mendelian randomization findings support a causal relation, also for serum phosphate and CAC in subjects without hyperphosphatemia, chronic kidney disease, and cardiovascular disease. Further research into underlying mechanisms of this association and sex differences is needed.

1,25‐dihydroxyvitamin D

25‐hydroxyvitamin D

alkaline phosphatase

fibroblast growth factor 23

genome wide association study

Mendelian randomization

prevalence ratio

parathyroid hormone

Rotterdam Study

What Is New?

Higher serum phosphate levels are associated with coronary artery calcifications in the general population, with stronger effects in men.

Mendelian randomization findings support that this association is causal, also for subjects with normal serum phosphate levels and without chronic kidney disease.

What Are the Clinical Implications?

Serum phosphate and coronary artery calcifications are associated in the absence of chronic kidney disease, hyperphosphatemia, and prevalent cardiovascular disease, challenging the concept that only severe hyperphosphatemia, in the setting of chronic kidney disease, is associated with coronary artery calcifications.

Our findings support a causal role for phosphate in the emerging epidemiological findings that higher serum phosphate levels are associated with increased mortality and cardiovascular events in the general population.

Arterial calcification is defined as the deposition of calcium and phosphate in the wall of arteries. , It was considered a passive consequence of aging until identified as a risk factor for cardiovascular events. Current evidence supports that a complex cellular‐mediated process underlies arterial calcification. At least 2 different layers can calcify: the intima, characteristic of atherosclerosis, and the media, typical of chronic kidney disease (CKD). The arterial mineralization process can result in hydroxyapatite [Ca10(PO4)6(OH)2] formation, as is found in bone. Coronary artery calcification (CAC) is one of the most studied calcification processes because of its specificity for atherosclerosis, its correlation to plaque burden and its ability to predict cardiovascular events. , Currently, it is widely accepted that CAC occurs mainly in the intima. Calcification mechanisms are multifactorial, but ectopic bone formation is considered the basis of CAC.

Serum phosphate has been related to arterial calcification in several human and animal disorders, , where genetically induced severe hyperphosphatemia leads to extensive calcification. Hyperphosphatemia‐induced calcification was described as the main mechanism of increased mortality in CKD. Similarly, the role of phosphate in CAC has been restricted mainly to hyperphosphatemia. However, Park et al reported an association of serum phosphate and CAC even within the normal phosphate range in Korean subjects with normal renal function. Interestingly, it has been reported that increasing yet normal serum phosphate is also a risk factor for cardiovascular morbidity and mortality in the general population. , These findings have been consistent, although a clear sex difference, with stronger associations in men, became evident. , The mechanisms underlying the associations between serum phosphate and cardiovascular morbidity and mortality and the reported sex differences remain unexplained.

Two non–mutually exclusive mechanisms have been described for serum phosphate in CAC: (1) passive deposition of calcium phosphate, inhibited by pyrophosphate; and (2) active induction of osteoblastic differentiation of vascular cells. Additionally, serum phosphate is regulated by several factors including parathyroid hormone (PTH), 1,25‐dihydroxyvitamin D3 (1,25(OH)2D3) and α‐klotho and fibroblast growth factor 23 (FGF23). PTH and 1,25(OH)2D3 exert inductive effects on arterial calcification. α‐klotho and FGF23 exert protective effects on arterial calcification. , Nevertheless, reports on the association between serum phosphate and CAC in the general population are not consistent as Gronhoj et al did not find an association between serum phosphate and prevalent CAC, but the same group reported an association of serum phosphate with progression of CAC over time. ,

As literature on serum phosphate and CAC in the general population is scarce and contradictive, , we aimed to analyze this association in the population‐based RS (Rotterdam Study) and test for potential sex differences. Because results from epidemiologic associations can be affected by reverse causation and unmeasured confounding, we also aimed to test for causality applying Mendelian randomization (MR), a statistical technique whereby genetic variants are used as instrumental variables for the exposure with the purpose of avoiding these sources of bias. If the relevant assumptions are valid, significant MR results can be interpreted as evidence of causality of the exposure on the outcome.

METHODS

Study Population

The RS is a prospective cohort study of men and women in the district of Ommoord, Rotterdam, designed to investigate the incidence and determinants of chronic disabling diseases. Rationale and design have been described elsewhere. The RS is now composed of 4 cohorts, named RS‐I, RS‐II, RS‐III, and RS‐IV (initiated in 1989, 2000, 2006, and 2017; total n ≈18 000 subjects). Subjects were assessed at baseline and through follow‐up visits. For the current study, we included participants from RS‐I (n=7983). Fasting serum phosphate and CAC were assessed during the second follow‐up visit of RS‐I. A total of 1889 subjects with both measurements available were included in the analyses. The RS was approved by the Medical Ethics Committee of the Erasmus MC (registration number MEC 02.1015) and by the Dutch Ministry of Health, Welfare and Sport (Population Screening Act WBO, license number 1071272‐159521‐PG). The Rotterdam Study was entered into the Netherlands National Trial Register ( www.trialregister.nl ) and into the World Health Organization International Clinical Trials Registry Platform (ICTRP; https://apps.who.int/trialsearch/ ) under shared catalog number NTR6831. All participants provided written informed consent to participate in the study and to have their information obtained from treating physicians. Data can be obtained upon request. Requests should be directed to the management team of the RS (datamanagement.ergo@erasmusmc.nl), which has a protocol for approving data requests. Because of restrictions based on privacy regulations and informed consent of the participants, data cannot be made freely available in a public repository.

Coronary Calcification Assessment

Coronary artery calcification was visualized using electron‐beam computed tomography (C‐150 Imatron Scanner, GE Healthcase, South San Francisco, CA). From the level of the aortic root through the heart, 38 images were obtained with a 100‐ms scan time and a 3‐mm slice thickness. During 1 breath hold, images were acquired at 80% of the cardiac cycle by using echocardiographic triggering. Quantification of coronary calcification was performed with AccuImage software (AccuImage Diagnostics Corporation, South San Francisco, CA). The presence of calcification was defined as a minimum of 2 adjacent pixels (area=0.65 mm2) with a density >130 HU. Following Agatston’s method, calcium scores were calculated by multiplying the area in square meters of individual calcified lesions with a factor based on the peak density of the lesion. The total score for the entire epicardial coronary vascular system comprised the sum of the scores for all individual lesions. Scores were log transformed (+1) to reduce the sensitivity to observations with extremely high CAC values.

Laboratory Measurements

Fasting blood samples were obtained in the second follow‐up visit, and serum phosphate was measured with a method based on the formation of ammonium phosphomolybdate that corresponds to the inorganic fraction of total phosphorus. Total calcium concentrations were assessed through a colorimetric o‐cresolphthalein complex one method (Merck Diagnostica, Amsterdam, the Netherlands). Levels of 25‐hydroxyvitamin D [25(OH)D] were determined through an electrochemiluminescense immunoassay, adjusting for seasonality through cosinor models. Creatinine concentrations were determined through a sarcosine‐based colorimetric assay and standardized. Subsequently, the Chronic Kidney Disease Epidemiology Collaboration equations were applied to estimate glomerular filtration rate. C‐reactive protein, glucose, cholesterol, and alkaline phosphatase (ALP) levels were measured through standard methods. , Ionized calcium was measured through a colorimetric detection assay using a Hitachi 917 (Roche, Mannheim, Germany). Assessments for ALP and ionized calcium were done at baseline visit and are therefore not simultaneous with serum phosphate.

Genotyping

Participants were genotyped in the Illumina HumanHap550 BeadChip SNP array. Variants were filtered on call rate <98%, minor allele frequency <0.01 and Hardy Weinberg equilibrium P<10−6, and imputed to the Haplotype Reference consortium panel, release 1.1. Genetic instruments for serum phosphate were selected from genome‐wide association analysis (GWAS) significant independent single‐nucleotide polymorphisms (SNPs) identified in the European GWAS by Kestenbaum et al and in a Japanese GWAS by Kanai et al. , Variants selected for analyses were checked for Hardy Weinberg equilibrium P>0.05 for genotyped SNPs, imputation quality for imputed SNPs (>0.8) and allele frequencies for palindromic SNPs to decrease the possibility of strand inconsistencies.

Other Covariates

Prevalent cardiovascular disease was defined as prevalent myocardial infarction, revascularization, stroke, and heart failure. Prevalent cardiovascular disease was assessed using general practitioners’ records and hospital discharge letters. Smoking status was assessed during home interviews. Blood pressure, height, and weight were measured during visits. Body mass index was calculated as weight in kilograms divided by height in square meters. Prevalent diabetes was assessed using general practitioners’ records, information on antidiabetic medication use, and fasting blood glucose levels.

Statistical Analysis

Phenotypic Associations

The association between serum phosphate and CAC as a continuous variable was assessed through linear regression models. The analysis was stratified according to estimated glomerular filtration rate (eGFR; <60 mL/min per 1.73 m2: CKD). Serum phosphate was assessed continuously and in quintiles. In addition, we tested the linearity of the association between serum phosphate and CAC using restricted cubic spline functions with 5 default knots (5th, 27.5th, 50th, 72.5th, and 95th percentile). , We tested if the calcium×phosphate product was associated with CAC, as previously reported in CKD. In addition, we assessed the association between serum phosphate and CAC as a categorical variable through prevalence ratios (PRs), for CAC scores >100, 300, 400, and 1000. We included interaction terms of phosphate with sex in age‐adjusted models to explore potential sex differences in the association between phosphate and CAC. When a sex difference in serum phosphate and CAC was confirmed, we performed sex‐stratified analysis. Because calcium is synergistic to phosphate in arterial calcification, we assessed the relation between serum phosphate and ionized calcium.

Model I included adjustments for age, body mass index, and smoking. Model II included also blood pressure, 25(OH)D, total serum calcium, C‐reactive protein, total cholesterol:high‐density lipoprotein cholesterol ratio and glucose levels, prevalent CVD and prevalent diabetes.

Sensitivity analyses

We restricted the analyses to subjects without hyperphosphatemia (hyperphosphatemia: serum phosphate>1.45 mmol/L, >4.5 mg/dL), without CKD and without prevalent CVD. To explore for interaction, we created 25(OH)D categories splitting at 48 nmol/L, which is a threshold related to CVD in the Dutch population. We assessed whether serum tissue nonspecific ALP was associated with CAC. ALP generates phosphate through hydrolysis of inorganic pyrophosphate, a potent arterial calcification inhibitor, and has been independently related to CAC in the general population, in CVD and in CKD.

Mendelian Randomization

To test for causality, we assessed whether phosphate was associated with continuous CAC. We selected 8 SNPs, assumed an additive model, and built a phosphate‐increasing allele score aligning the alleles: a higher score predicted a higher serum phosphate. Scores (as a single instrument) are associated with lower risk of weak instrument bias than the simultaneous use of multiple SNPs. The SNPs included in the score were derived from a GWAS meta‐analysis within the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium composed of individuals of European ancestry, and a recent GWAS within the Biobank Japan Project composed of individuals of East Asian ancestry. ,

In the former, the population under study was part of the discovery sample of the GWAS, which could result in bias from winner’s curse. However, all SNPs reported in the above‐mentioned GWA studies have been replicated in a UK biobank GWAS. Summary statistics of this GWAS are publicly available.

To properly account for uncertainty in imputed SNPs, genotypes were extracted from dosage files and therefore its values span between 0 and 2, reflecting the probability of getting up to 2 risk (phosphate‐increasing) alleles. If the SNP was genotyped, we report Hardy Weinberg equilibrium test; otherwise, the imputation quality is displayed.

Palindromic SNPs were checked for allele frequency concordance between RS and the GWAS catalog. In addition, one SNP mapping to the nonpseudoautosomal region of the X chromosome was coded as 0/2 in men (for 0/1 risk allele) and as 0/1/2 in women, following recent guidelines when assuming a pattern of X‐chromosome inactivation.

Inclusion of correlated SNPs, in linkage disequilibrium, is a potential source for bias in standard MR analyses without covariance matrix, leading to increased type I error rates. Therefore, we included only independent SNPs in the score, applying a threshold of r 2<0.01.

We used the score to genetically predict serum phosphate and tested for causality applying a 2‐stage least square regression, where the first stage regresses the exposure on the instruments and the second stage regresses the outcome on the fitted values of the exposure estimated in first stage. We applied robust standard errors. We assessed MR assumptions as follows:

Assumption N°1: the instrument must be associated with the exposure, relevance condition. We regressed serum phosphate levels on phosphate SNPs/score in the population with serum phosphate levels available and with serum phosphate and CAC levels available. We considered β, P, and F‐statistic. In general, an instrument with an F‐statistic <10 is considered a weak instrument and a higher F‐statistic reflects a better instrument. However, no cutoffs should be made, as bias is a continuous phenomenon. ,

Assumption N°2: The instrument must not be associated with potential confounders, independence condition. We regressed potential confounders on phosphate scores. Evidently, it is not possible to assess the association of instruments with unmeasured confounders.

Assumption N°3: The instrument must be related to the outcome only through the exposure, meaning the absence of horizontal pleiotropy, exclusion‐restriction condition. To test this assumption, both a frequentist and a Bayesian approach were applied, the latter as sensitivity analysis. MR‐Egger regression was not applied because of our 1‐sample setting and its low statistical power when the SNP‐exposure associations are homogenous. We implemented instead an adaptive lasso regression (sivreg, Stata) that provides estimates while allowing less than 50% of instruments to be invalid by horizontal pleiotropy.

Analyses were performed with unweighted scores. Although weighted scores may increase power, use of internal weights derived from the data should be avoided because of the severe bias that this approach induces and the RS represented ≈40% of the Cohorts for Heart and Aging Research in Genomic Epidemiology GWAS. ,

To test if results from the score were driven by only 1 SNP, we applied the leave‐1‐out approach, excluding 1 SNP at a time from the score and testing for each reduced score whether genetically predicted phosphate was still associated with CAC. This penalization technique is considered a robust method : if results are driven by only 1 SNP, a high index of pleiotropy should be suspected and properly assessed. To obtain results for the Sargan test, which is an indirect measure of heterogeneity among instruments, we did not apply robust standard errors to the leave‐1‐out analysis.

Similar to phenotypic associations, we performed subgroup analyses excluding participants with hyperphosphatemia, CKD, and prevalent CVD.

Sensitivity analysis

In addition of applying a single score, we genetically predicted phosphate through the combination of all SNPs simultaneously (joint instruments analyses). This approach might have more power but might suffer from weak instrument bias. We applied in this setting the Sargan test, an overidentification test, to assess whether all instruments included in the regression are valid in linear combination and not correlated with error terms. Additionally, Sargan provides a test of heterogeneity among instruments. We also applied the leave‐1‐out approach to the joint instruments analysis.

Consistent with ancestry of participants, we built a score derived only from SNPs from the Cohorts for Heart and Aging Research in Genomic Epidemiology meta‐analysis (EUR‐score) and tested for reproducibility of results.

Finally, we applied a Bayesian approach designed for a 1‐sample setting that allows several SNPs to exert pleiotropic effects and that does not rely in the assumption of no correlation between SNP strength and pleiotropic effects. The method is described in Data S1. This approach allows the implementation of pleiotropic effects for a subset (49%) of SNPs, incorporated in their prior distribution, and applies variational Bayes through a modified Markov chain Monte Carlo to estimate the posterior mean and its 95% CI. Genotypes were centered to improve convergence. Analyses were performed after imputing missing values through multiple imputation with chained equations. We used SPSS version 21.0 (IBM Corp, Armonk, NY), Stata version 16 (Stata Corp LP, College Station, TX) and R version 3.5.0 (R Foundation for Statistical Computing, Vienna, Austria). A 2‐sided P<0.05 was considered significant.

RESULTS

The general characteristics of the study population are shown in Table 1. Mean serum phosphate levels were higher in women than in men (P women: 1.18 mmol/L [3.65 mg/dL]; P men: 1.02 mmol/L [3.16 mg/dL]; P t‐test<0.001 [data not shown]) (reference range of serum phosphate, 0.8–1.4 mmol/L=2.5–4.5 mg/dL). Median CAC score levels were higher in men than in women (CACmen: 310.2; CACwomen: 54.7; P M Whitney<0.001). Total calcium and calcium×phosphate product was positively associated with serum phosphate in both sexes while ionized calcium was negatively associated with serum phosphate in women only. All results are expressed per 1‐SD increase of serum phosphate (0.16 mmol/L=0.49 mg/dL), unless otherwise stated.

Table 1

General Characteristics of Study Population, per Quintiles of Fasting Phosphate Levels

	Men					Women
	Phosphate in quintiles					Phosphate in quintiles
	1	2	3	4	5	1	2	3	4	5
N (phosphate in mmol/L)	175 (0.83)	176 (0.95)	175 (1.02)	176 (1.09)	176 (1.21)	202 (0.98)	202 (1.11)	202 (1.17)	202 (1.24)	203 (1.37)
Age, y	70.8	71.3	70.6	70.6	70.7	70.7	70.2	71.2	70.8	70.6
BMI, kg/m2	26.8	26.7	26.3	26.8	26.1	29.1	27.6	27.3	26.9	26.3
Ever smoke, %	93.1	89.1	93.1	94.3	94.3	48.8	52.0	55.9	58.2	53.0
Systolic BP, mm Hg	145.1	143.4	148.0	146.1	144.0	145.8	143.9	147.9	139.3	142.5
Diastolic BP, mm Hg	76.9	76.5	82.5	78.1	76.5	78.4	75.8	79.3	73.9	74.4
Ionized Ca++, mmol/L*	1.29	1.29	1.29	1.29	1.28	1.31	1.30	1.29	1.30	1.28
Calcium, mmol/L	2.39	2.41	2.40	2.40	2.42	2.43	2.43	2.44	2.44	2.46
Ca×P product, mmol2/L2	1.98	2.28	2.44	2.61	2.94	2.38	2.69	2.87	3.04	3.38
ALP, U/L*	79.6	76.5	75.6	74.4	75.0	80.5	76.7	77.7	80.3	76.5
25(OH)D, nmol/L	66.7	63.1	65.6	60.8	60.5	49.6	51.5	48.4	49.2	50.8
CAC score	446.7	668.5	696.9	840.3	1059.8	234.2	193.4	319.0	342.6	279.0
CRP, mg/L	3.95	3.37	3.74	3.64	5.20	4.78	3.77	3.57	3.36	3.05
Glucose, mmol/L	6.08	6.00	6.14	6.09	6.02	6.12	5.77	5.80	5.77	5.64
eGFR, mL/min per 1.73 m2	72.7	72.5	73.9	73.1	74.7	71.8	72.6	71.0	71.8	72.8
Chol to HDL ratio	4.78	4.86	4.82	4.65	4.45	4.39	4.40	4.36	4.21	4.10
Prevalent CVD, %	12.0	17.6	16.6	17.0	19.9	6.0	7.4	4.5	8.4	6.4
Prevalent diabetes, %	13.1	13.6	13.7	17.0	15.9	16.8	11.9	12.4	10.4	9.9

Continuous values are displayed as means; categorical variables are displayed in percentages. 25(OH)D indicates 25‐hydroxyvitamin D levels; ALP, alkaline phosphatase levels; BMI, body mass index; BP, blood pressure; Ca×P product, total calcium×phosphate levels; calcium, total calcium levels; Chol to HDL ratio, total cholesterol to HDL cholesterol ratio; CRP, C‐reactive protein; eGFR, estimated glomerular filtration rate; ionized Ca++, ionized calcium levels; and prevalent CVD, prevalent cardiovascular disease.

Ionized calcium and alkaline phosphatase (ALP) levels were not assessed simultaneously with serum phosphate levels.

Results From Phenotypic Associations

CAC as Continuous Trait

Table 2 shows that serum phosphate was associated with CAC in the total population after adjustments for model I (β, 0.37 [95% CI, 0.26–0.48]; P=3×10−11; n=1889). We found a significant interaction between serum phosphate and CAC across sexes (P interaction=0.003). Sex‐stratified analyses showed that the association between serum phosphate and CAC was stronger in men (β, 0.52 [95 CI, 0.38–0.67]; P=5×10−12; n=878) than women (β, 0.22 [95 CI, 0.06–0.38]; P=0.006; n=1011). Further adjustments (model II) induced a slight attenuation of the association in men (men β, 0.44 [95 CI, 0.30–0.59]; P=3×10−9; women β, 0.24 [95 CI, 0.08–0.40]; P=0.003).

Table 2

Association Between Serum Phosphate Levels and Coronary Artery Calcification Scores

	Model I			Model II
	n	β (95% CI)*	P value	n	β (95% CI)*	P value
Men	878	0.52 (0.38–0.67)	<0.001	878	0.44 (0.30–0.59)	<0.001
Women	1011	0.22 (0.06–0.38)	0.006	1011	0.24 (0.08–0.40)	0.003
Total	1889	0.37 (0.26–0.48)	<0.001	1889	0.34 (0.23–0.45)	<0.001

BMI indicates body mass index; and HDL, high‐density lipoprotein.

βs were obtained from linear regression models and expressed per 1‐SD increase in phosphate (0.16 mmol/L=0.49 mg/dL). Model I: adjusted for age, BMI, smoking. Model II: adjusted for age; BMI; blood pressure; smoking; prevalent cardiovascular disease; prevalent diabetes; and serum levels of 25‐hydroxyvitamin D, total calcium, C‐reactive protein, total cholesterol to HDL cholesterol ratio, and glucose.

The stratified analyses (Table 3) showed that serum phosphate was associated with CAC across the spectrum of kidney function in men (eGFR ≥60 mL/min per 1.73 m2 β, 0.53 [95 CI, 0.35–0.70]; P<0.001; n=736; eGFR <60 mL/min per 1.73 m2 β, 0.53 [95 CI, 0.31–0.75]; P<0.001; n=142), while in women this association was constrained to normal eGFR (eGFR ≥60 mL/min per 1.73 m2 β, 0.22 [95 CI, 0.04–0.39]; P=0.016; n=839; eGFR <60 mL/min per 1.73 m2 β, 0.25 [95 CI, −0.17 to 0.66]; P=0.238; n=172). Adjustments for model II induced a slight attenuation in men (eGFR ≥60 mL/min per 1.73 m2 β, 0.44 [95 CI, 0.27–0.62]; P<0.001; eGFR<60 mL/min per 1.73 m2 β, 0.45 [95 CI, 0.21–0.68]; P=0.001).

Table 3

Association Between Serum Phosphate Levels and Coronary Artery Calcification Scores, Stratified by eGFR

	eGFR ≥60 mL/min per 1.73 m2 *			eGFR <60 mL/min per 1.73 m2 *
	n	β (95% CI)†	P value	n	β (95% CI)†	P value
Model I
Men	736	0.53 (0.35 to 0.70)	<0.001	142	0.53 (0.31 to 0.75)	<0.001
Women	839	0.22 (0.04 to 0.39)	0.016	172	0.25 (−0.17 to 0.66)	0.238
Total	1575	0.36 (0.24 to 0.49)	<0.001	314	0.42 (0.20 to 0.64)	<0.001
Model II
Men	736	0.44 (0.27 to 0.62)	<0.001	142	0.45 (0.21 to 0.68)	<0.001
Women	839	0.22 (0.05 to 0.40)	0.011	172	0.30 (−0.12 to 0.72)	0.154
Total	1575	0.33 (0.21 to 0.46)	<0.001	314	0.36 (0.14 to 0.58)	0.002

BMI indicates body mass index; eGFR, estimated glomerular filtration rate; and HDL, high‐density lipoprotein.

eGFR estimated from creatinine‐based Chronic Kidney Disease Epidemiology Collaboration equations.

βs were obtained from linear regression models and expressed per 1‐SD increase in phosphate (0.16 mmol/L=0.49 mg/dL). Model I: adjusted for age, BMI, smoking. Model II: adjusted for age; BMI; blood pressure; smoking; prevalent cardiovascular disease; prevalent diabetes; and serum levels of 25‐hydroxyvitamin D, total calcium, C‐reactive protein, total cholesterol to HDL cholesterol ratio, and glucose.

The analyses in quintiles suggested a threshold for the association of serum phosphate and CAC (Table 4): setting the first quintile as reference, men with serum phosphate above 1.09 mmol/L displayed a significant trend for higher CAC (β for fourth quartile, 0.87 [95 CI, 0.46–1.28]; P<0.001; β for fifth quintile, 1.18 [95 CI, 0.77–1.59]; P<0.001; P trend=3×10−10; n=878). For women the threshold was >1.37 mmol/L (β for fifth quintile, 0.67 [95 CI, 0.22–1.11]; P=0.003; P trend=0.002; n=1011). A 5‐knot restricted cubic spline function did not find evidence for a nonlinear association between serum phosphate and CAC in the total population and in men and women separately (Figure S1).

Table 4

Association Between Serum Phosphate Levels and Coronary Artery Calcification Scores, per Quintiles of Phosphate Levels

Men				Women
n	Phosphate levels mean (range)*	β (95% CI)†	P value	n	Phosphate levels mean (range)*	β (95% CI)†	P value
175	0.83 (0.63–0.91)	1 (Ref)		202	0.98 (0.74 to 1.06)	1 (Ref)
176	0.95 (0.91 to 0.98)	0.28 (−0.13 to 0.69)	0.178	202	1.11 (1.06 to 1.14)	0.04 (−0.40 to 0.47)	0.869
175	1.02 (0.98 to 1.05)	0.37 (−0.04 to 0.78)	0.078	202	1.17 (1.14 to 1.20)	0.03 (−0.41 to 0.47)	0.889
176	1.09 (1.05 to 1.13)	0.87 (0.46 to 1.28)	<0.001	202	1.24 (1.20 to 1.28)	0.26 (−0.18 to 0.69)	0.247
176	1.21 (1.13 to 2.47)	1.18 (0.77 to 1.59)	<0.001	203	1.37 (1.28 to 1.70)	0.67 (0.22 to 1.11)	0.003
P trend		<0.001				0.002

BMI indicates body mass index.

Phosphate quintiles are expressed in mmol/L.

Betas were obtained from linear regression models. First quintile of phosphate was set as reference. Analyses were adjusted for age, BMI, and smoking.

Similarly, the results of the calcium×phosphate product with CAC (Table 5) suggested threshold values: setting the first quintile as reference, men with a product above 2.44 mmol2/L2 displayed a significant trend for higher CAC (β for third quintile, 0.49 [95 CI, 0.09–0.90]; P=0.017; β for fourth quintile, 0.77 [95 CI, 0.37–1.18; P<0.001; β for fifth quintile, 1.17 [95 CI, 0.77–1.58]; P<0.001; P trend=5×10−11; n=878). In women the threshold was above 3.05 mmol2/L2 (β for fourth quintile, 0.47 [95 CI, 0.04–0.91]; P=0.034; β for fifth quintile, 0.64 [95 CI, 0.19–1.08]; P=0.005; P trend=0.001; n=1011). All coefficients are per unit increase in calcium×phosphate product. ALP was not associated with CAC (men β, −0.001 [95 CI, −0.009 to 0.007]; P=0.873; n=574; women β, −0.005 [95 CI, −0.012 to 0.001]; P=0.094; n=747) (data not shown).

Table 5

Association Between Serum Calcium×Phosphate Product Levels and Coronary Artery Calcification Scores, per Quintiles Of Calcium×Phosphate Product Levels

Men				Women
N	Product mean (range)*	β (95% CI)†	P value	n	Product mean (range)*	β (95% CI)†	P value
175	1.97 (1.50 to 2.16)	1 (Ref)		202	2.35 (1.67 to 2.57)	1 (Ref)
176	2.27 (2.16 to 2.36)	0.03 (−0.37 to 0.44)	0.868	202	2.68 (2.58 to 2.77)	0.001 (−0.43 to 0.44)	0.995
175	2.44 (2.36 to 2.52)	0.49 (0.09 to 0.90)	0.017	202	2.86 (2.77 to 2.96)	0.29 (−0.14 to 0.73)	0.185
176	2.62 (2.52 to 2.71)	0.77 (0.37 to 1.18)	<0.001	202	3.05 (2.96 to 3.15)	0.47 (0.04 to 0.91)	0.034
176	2.96 (2.71 to 6.57)	1.17 (0.77 to 1.58)	<0.001	203	3.40 (3.16 to 4.20)	0.64 (0.19 to 1.08)	0.005
P trend		<0.001				0.001

BMI indicates body mass index; and HDL, high‐density lipoprotein.

Calcium×phosphate product levels are expressed in mmol2/L2.

Betas were obtained from linear regression models. First quintile of calcium×phosphate product level was set as reference. Analyses were adjusted for age, BMI, and smoking.

CAC as a Categorical Trait

After adjustments for model II (Table S1), each 1‐SD increase in serum phosphate (0.16 mmol/L) was associated with an increased PR for CAC >100 of 8% in men only (95 CI, 4%–12%]; P=6×10−5); a PR for CAC >300 of 10% in men only (95 CI, 5%–15%; P=2×10−4); a PR for CAC >400 of 10% in men only (95 CI, 5%–16%; P=2×10−4), and a PR for CAC >1000 of 20% in men (95 CI, 12%–28%; P=4×10−7) and 36% in women (95 CI, 18%–55%; P<0.001).

Sensitivity Analysis

After exclusion of subjects with hyperphosphatemia, serum phosphate was associated with CAC in both men (β, 0.53 [95 CI, 0.37–0.69]; P=1×10−10; n=873) and women (β, 0.21 [95 CI, 0.03–0.40]; P=0.020; n=974) (Table S2). In men, results from model II showed a slight attenuation (men β, 0.46 [95 CI, 0.31–0.62]; P=1×10−8). Moreover, exclusion of CKD and hyperphosphatemia did not change the association between phosphate and CAC (men β, 0.53 [95 CI, 0.35–0.71]; P=1×10−8; n=733; women β, 0.22 [95 CI, 0.03–0.42]; P=0.026; n=808). These results were slightly attenuated in men after adjustments for model II (men β, 0.45 [95 CI, 0.27–0.63]; P=1×10−6). Finally, after exclusion of CKD, hyperphosphatemia and prevalent CVD, the association between serum phosphate and CAC remained (men β, 0.55 [95 CI, 0.35–0.74]; P=4×10−8; n=627; women β, 0.20 [95 CI, 0.0003–0.40]; P=0.050; n=765). Adjustments for model II yielded similar results (men β, 0.50 [95 CI, 0.31–0.69]; P=5×10−7; women β, 0.22 [95 CI, 0.02–0.43]; P=0.029). We found no evidence of effect modification by 25(OH)D (data not shown).

Results From MR Analyses

No evidence of departure from Hardy Weinberg equilibrium was observed in genotyped SNPs; for imputed SNPs, the imputation quality was above 96%. No frequency/strand inconsistency between the original GWAS and RS was detected for the only palindromic SNP (rs2970818; Adenine/Thymine); therefore, it was included in the score (Table S3).

Concerning MR first assumption, the strength of the instruments was tested regressing serum phosphate levels (expressed in entire units of mmol/L) on the SNPs in all subjects with serum phosphate levels available and further restricted to those with serum phosphate and CAC (Table S4). The allele score was significantly associated with serum phosphate in both the whole cohort with available phosphate levels and in the subset of the population with available phosphate and CAC levels. The F‐statistic for the allele score was >25, meaning results from MR analyses should not be affected by weak instrument bias. Three SNPs from Biobank Japan showed a significant association with serum phosphate in the whole cohort but further restriction to subjects with serum phosphate and CAC decreased power, and the associations were no longer significant. Therefore, we first built a score including all 8 SNPs and subsequently applied sensitivity analyses.

Concerning MR second assumption, we regressed potential confounders on the phosphate scores: the 8‐SNP score and the 5‐SNP EUR‐score (Table S5). We found no association of the scores with total calcium, 25(OH)D, C‐reactive protein, total cholesterol : high‐density lipoprotein cholesterol ratio and glucose; nor with body mass index, smoking, prevalent CVD, or DM.

Allele Score Analyses

Figure 1 shows the results of the 2‐stage least square regression regression of phosphate, genetically predicted through the unweighted score, and CAC adjusted for age, sex, and 10 principal components in the whole cohort. The allele score has been scaled to be associated to 1 SD of serum phosphate (0.16 mmol/L=0.49 mg/dL). A significant relation was found between the unweighted phosphate score and CAC (β, 0.93 [95 CI, 0.05–1.82]; P=0.039).

Figure 1

Mendelian randomization results for serum phosphate and coronary artery calcification: allelic score method and leave‐1‐SNP‐out approach applied to the whole cohort.

Betas were derived from 2‐stage least square for the score as a single instrument and adjusted for age, sex, and 10 principal components. Results are expressed as change in outcome per 1‐SD increase in phosphate (0.16 mmol/L=0.49 mg/dL). Leave‐1‐SNP‐out approach: allelic score analyses with the subtraction of 1 SNP at‐a‐time. Closest annotated gene is displayed if known to be associated with (or possible related to) phosphate homeostasis.

Sex‐stratified analysis suggested that the association between genetically predicted phosphate and CAC was more consistent in men than in women (men β, 1.31 [95 CI, −0.02 to 2.64]; P=0.053; n=782; women β, 0.56 [95 CI, −0.61 to 1.75]; P=0.347; n=911).

When we applied the leave‐1‐SNP‐out approach (from the score), we found that results lost significance after extracting rs1697421 (β, 0.48 [95 CI, −0.51 to 1.46]; P=0.344); rs2970818 (β, 0.96 [95 CI, −0.02 to 1.93]; P=0.054) and rs35186465 (β, 0.85 [95 CI, −0.06 to 1.77]; P=0.068), one at a time.

Subgroup Analysis According to Serum Phosphate Levels, Kidney Function, and CVD

The unweighted phosphate score remained associated to CAC after exclusion of hyperphosphatemia (β, 1.10 [95 CI, 0.10–2.10]; P=0.031; n=1659); after exclusion of hyperphosphatemia and CKD (β, 1.30 [95 CI, 0.33–2.27]; P=0.009; n=1377); and after exclusion of hyperphosphatemia, CKD, and prevalent CVD (β, 1.23 [95 CI, 0.17–2.28]; P=0.023; n=1244) (Figure 2, Table S6).

Figure 2

Mendelian randomization results for serum phosphate and coronary artery calcification: allelic score method applied in subgroup analyses according to serum phosphate levels, kidney function, and prevalent cardiovascular disease.

Betas were derived from 2‐stage least square for the score as a single instrument and adjusted for age, sex, and 10 principal components. Results are expressed as change in outcome per 1‐SD increase in phosphate (0.16 mmol/L=0.49 mg/dL). CKD indicates chronic kidney disease, defined as a glomerular filtration rate <60 mL/min per 1.73 m2. Prevalent CVD, prevalent cardiovascular disease, defined as prevalent myocardial infarction, revascularization, stroke, and heart failure; HyperP, hyperphosphatemia, defined as a phosphate level >1.45 mmol/L (=4.5 mg/dL).

Joint Instruments Analysis and Sargan Statistics

Phosphate genetically predicted through all 8 SNPs simultaneously (Table S7) was associated with CAC (β, 0.81 [95 CI, 0.04–1.58]; P=0.038). The exclusion of 1 SNP at a time yielded similar results as their exclusion from the score: phosphate is not associated with CAC if the following SNPs are excluded, 1 at a time: rs1697421, rs2970818, and rs35186465.

The Sargan tests could not reject the null hypothesis in any case, providing validity of the instruments and an indirect evidence of low heterogeneity among them.

EUR‐Score

Table S6 shows that restriction to EUR‐score did not attenuate the results of genetically predicted phosphate and CAC (β, 1.04 [95 CI, 0.07–2.01]; P=0.036; n=1693). Results remained significant after excluding hyperphosphatemia (β, 1.06 [95 CI, 0.04–2.08]; P=0.041; n=1659), after excluding hyperphosphatemia and CKD (β, 1.30 [95 CI, 0.29–2.30]; P=0.011; n=1377), and after excluding hyperphosphatemia, CKD, and prevalent CVD (β, 1.05 [95 CI, 0.08–2.03]; P=0.035; n=1244).

Assessment of Potential Horizontal Pleiotropy

Concerning MR’s third assumption (Table S8), from a frequentist approach the adaptive lasso regression found no evidence of invalidity of instruments; therefore, the inference was similar as that obtained from 2‐stage least square regression. There was no rejection of the Hansen test, meaning that all SNPs are valid and uncorrelated with error terms.

From a Bayesian approach, we applied a method that incorporates pleiotropic effects into a fraction of SNPs and tests whether the association between the genetically predicted exposure and outcome remains. Allowing 49% of SNPs to exert pleiotropic effects and assuring model convergence, we found that genetically predicted phosphate was still associated with CAC (posterior mean, 0.94 [95% CI,] 0.13–1.45).

DISCUSSION

Our analyses showed that serum phosphate was strongly associated with CAC, even after excluding subjects with hyperphosphatemia, CKD, and prevalent CVD, in a Dutch population‐based cohort study. The implementation of MR methods, where genetically predicted phosphate is consistently associated with CAC, strengthens the inference of our observational findings and supports causality. It has been previously shown that hyperphosphatemia or a supraphysiological high‐phosphate medium is able to induce an osteoblastic transformation of vascular smooth muscle cells with subsequent medial calcification, , characteristic of CKD. Nevertheless, whether normal serum phosphate is associated with intimal calcification has been much less explored. , Our data demonstrated important sex differences in this association, adding to sex differences previously found concerning serum phosphate and all‐cause mortality, CVD mortality, and atherosclerosis, where consistently a stronger (or unique) relation has been described in men.

When analyzed as a categorical trait, serum phosphate increases the PR in men for CAC >100, which is considered of moderate risk. In men, serum phosphate increases PR for CAC >300 and >400, which is considered high risk. Remarkably, serum phosphate increases PR for CAC >1000 in both sexes, a category that confers a high mortality risk. To the best of our knowledge, this is the first study to assess the relation between serum phosphate and CAC through MR, by definition less prone to be affected by confounding and reverse causation. An important source of bias in MR might be horizontal pleiotropy. Nevertheless, we applied elaborated regression models that assess or allow pleiotropy. , The persistence of similar results to standard 2‐stage least square regression, through lasso regression, and the obtainment of significant results despite allowing almost half of the SNPs to exert pleiotropic effects, through Bayesian modeling, confirm the robustness of our findings.

The association between serum phosphate and CAC after restriction of MR to subjects without hyperphosphatemia, CKD and prevalent CVD, supports that increasing (but within normal range) serum phosphate in the general population without clinical CVD is a pathogenic factor for increasing the CAC burden. This finding challenges the concept that only severe hyperphosphatemia, in the uremic context of CKD, is associated with CAC. More importantly, it might provide an explanation for the emerging epidemiologic associations of serum phosphate and increased mortality and cardiovascular events in cohorts with mostly normal serum phosphate; and to CVD mortality and atherosclerosis in men with strict normal serum phosphate. , If these associations are causal, there must be an underlying mechanism. CAC induction by an increasing, yet normal, serum phosphate might be one of these mechanisms.

The approach of leaving 1 SNP out has recently been acknowledged as a robust penalization method to test validity in MR. We found that results were not significant when specific SNPs were omitted from the score, 1 at a time:

rs1697421: Its omission results in the nullification of the association. This SNP is intergenic, but its positional candidate gene is ALPL, which encodes for tissue‐nonspecific ALP. ALP was not associated with CAC in our population study; but it hydrolyzes pyrophosphate into phosphate. Pyrophosphate is one of the most potent calcification inhibitors. The condition where an SNP affects the outcome through a pathway affected by the risk factor of interest is termed vertical pleiotropy and does not invalidate MR findings. If this SNP influences ALP and its downstream activity/levels, pyrophosphate and phosphate, it will correspond to mediation of the effect.

rs2970818: This SNP is also intergenic but one of the positional candidate genes is FGF23, which encodes for a key hormone in phosphate homeostasis through increased renal phosphate excretion. In contrast to observational studies linking higher FGF23 levels to arterial calcification, research at the cellular level has shown that FGF23 inhibits osteoblastic differentiation of vascular cells, partially through α‐klotho actions. , Therefore, horizontal pleiotropy is unlikely.

rs35186465: There are not any known phosphate‐related genes annotated to this SNP.

Therefore, the association of genetically predicted phosphate with CAC is explained mostly by the contribution of 3 SNPs from the allele score located in chromosomes 1, 12, and 17. Though a role from several SNPs located throughout the genome improves the validity from MR, it seems that rs1697421 near ALPL plays a key role.

Besides FGF23 and α‐klotho, serum phosphate is regulated by 1,25(OH)2D3 and PTH levels. Both 1,25(OH)2D3 and PTH are positively associated with CAC. The significant results from MR analyses decrease their likelihood as confounders. Nevertheless, as PTH has been related to arterial calcification even at normal levels and PTH increases with increased serum phosphate, our data cannot rule out a role of PTH on CAC.

Two main pathways of phosphate‐induced calcification have been described in the coronary bed: (1) a passive deposition of calcium and phosphate, strongly regulated by ALP‐pyrophosphate‐P, and (2) an active process of osteoblastic differentiation of vascular pericytes and calcifying vascular cells, able to synthesize matrix vesicles, which start the mineralization process. Current evidence has shown that ALP, pyrophosphate and phosphate are present in matrix vesicles surfaces of atherosclerotic plaques, linking closely both mechanisms of calcification in CAC and potentially providing a biological explanation for our results.

Although a restricted cubic spline model did not find evidence for a nonlinear association, we found an apparent dose‐effect relation in phosphate and CAC, with normal serum phosphate thresholds of 1.09 mmol/L (3.35 mg/dL) and 1.37 mmol/L (4.24 mg/dL) in men and women. Interestingly, Dhingra et al (Framingham study) described a close cutoff for serum phosphate of 3.5 mg/dL (1.13 mmol/L) above which CVD mortality and morbidity increased. The authors stated that it was not clear whether increasing serum phosphate within normal range was associated with CVD risk. Our data suggests that this question can be answered in a confirmatory way.

We also found that normal levels of the calcium×phosphate product were associated with CAC. Recent literature highlights that circulating calciprotein particles, composed of calcium, phosphate, and calcification inhibitors such as fetuin A, are crucial in calcification and that its composition dictates whether pathologic mineralization is inhibited or not. The calcium×phosphate product within calciprotein particles has been identified as the culprit in this process. Similar to serum phosphate, it might be that a normal product in serum does not reflect a safe product at the cellular level.

The stronger associations observed in men are consistent to previous research on serum phosphate and CVD (atherosclerosis, cardiovascular event rates and CVD mortality , ). These results are unexpected, especially because women have higher serum phosphate and because the protective effect of 17β‐estradiol in arterial calcification is predominant in premenopausal women. We can only speculate whether the association between calcium and phosphate levels plays a role in the sex difference as we (and others) have found an inverse relation between serum phosphate and ionized calcium in women but not in men, and an inverse relation between them seems necessary to keep a constant calcium×phosphate product in serum. , It is important to add that phosphate intake has also been related to arterial calcification, as an abrupt postprandial phosphate increase suffices to initiate mineralization within seconds and to decrease α‐klotho expression.

This study has several limitations. We had a small sample size and no measurements of 1,25(OH)2D3, FGF23, and PTH levels, nor information on phosphate intake. CAC and ALP were not determined simultaneously. Several tests were not suitable because of our 1‐sample MR. Moreover, the RS is a Dutch population‐based cohort study, precluding inference to other populations or ethnic groups. Finally, we did not have prospective data on CAC and coronary events, which prevents us from drawing conclusions about the clinical implications of our findings. But there are several strengths, especially concerning the results from MR, that provides a formal test of causality provided the assumptions are fulfilled. Results from F‐statistics strongly suggest that our results are not affected by weak‐instrument bias. We were able to perform important stratified and subgroup analyses and to test instruments’ validity.

To conclude, we hereby provide both frequentist and Bayesian evidence from the MR approach that normal phosphate is a causative factor in CAC in the general population without hyperphosphatemia, without prevalent CVD and with normal kidney function. We add more evidence to support the concept of phosphotoxicity, and our results call for a review of the current normal serum phosphate range. We agree with the European Food and Safety Agency that more research is needed to study the relationships between dietary intake of phosphate and serum phosphate levels and adverse health outcomes. Public health policies might be needed to decrease phosphate intake because of the growing evidence of phosphate as a continuous risk factor for adverse outcomes such as atherosclerosis, CVD mortality, and, now, CAC. Further research should focus on unveiling the underlying mechanisms of the detrimental effects of phosphate in human health and to establish a threshold above which phosphate must be considered harmful for men and women, especially because a large fraction of the population appears to be exposed to nonsafe P levels.

Sources of Funding

The generation and management of GWAS genotype data for the Rotterdam Study (RS‐I, RS‐II, RS‐III) was executed by the Human Genotyping Facility of the Genetic Laboratory of the Department of Internal Medicine, Erasmus Medical Center, Rotterdam, the Netherlands. The GWAS data sets are supported by the Netherlands Organization of Scientific Research Investments (No. 175.010.2005.011, 911‐03‐012); the Genetic Laboratory of the Department of Internal Medicine, Erasmus Medical Center; the Research Institute for Diseases in the Elderly (014‐93‐015); and the Netherlands Genomics Initiative/Netherlands Organization for Scientific Research Netherlands Consortium for Healthy Aging, Project No. 050‐060‐810. The Rotterdam Study is funded by Erasmus Medical Center and Erasmus University, Rotterdam; Netherlands Organization for the Health Research and Development; the Research Institute for Diseases in the Elderly; the Ministry of Education, Culture and Science; the Ministry for Health, Welfare and Sports; the European Commission (DG XII); and the Municipality of Rotterdam. Bosman is supported by a grant from Health~Holland (PhosphoNorm; LSHM18029).

Disclosures

None.

Data S1

Tables S1–S8

Figure S1

Click here for additional data file.