Association of dietary vitamin K and risk of coronary heart disease in middle-age adults: the Hordaland Health Study Cohort
Haugsgjerd, T.R., Egeland, G.M., Nygard, O.T., et al.
Introduction
Vitamin K is a fat-soluble vitamin including vitamin K1 (K1; phylloquinone) from green leafy vegetables and vegetable oils as the main dietary sources, and vitamin K2 (K2; menaquinones) from dairy products, meat and egg yolk as the main dietary sources in Europe.1–3 K2 has a longer half-life in the circulation than K1.4 Both are absorbed from the small ileum and jejunum. K1 and K2 are incorporated into chylomicrons and delivered to the liver. K2 is also transported via low-density lipoprotein and high-density lipoprotein (HDL) particles to extrahepatic tissue.4 5
Vitamin K functions as a cofactor for the enzyme gamma-glutamyl carboxylase which converts protein-bound glutamate residues into gammacarboxyglutamate (Gla).6 7 Gla-containing proteins are involved in, for example, the coagulation of blood,8 inhibition of arterial calcification (Matrix Gla Protein) and vascular smooth muscle cell apoptosis and movement that is considered protective against vascular injury (Gas-6).9 Matrix Gla Protein is involved in both medial and intimal calcification, and low vitamin K status has been associated with both types of calcification.10–14 In addition, a study that examined the effect of warfarin (a vitamin K antagonist) on medial and intimal plaque calcification in apoE−/− mice concluded that warfarin accelerates both medial and intimal calcification of atherosclerotic plaque.15 Patients with both medial and intimal calcification have a higher cardiovascular risk when compared with similar patients without calcification.16 17 Therefore, an inverse association between vitamin K intake and coronary heart disease (CHD) could be expected. Results from observational studies on the association between intake of vitamin K and CHD are inconsistent.18–24 Among the identified studies, three found reduced risk of CHD in multivariable adjusted analyses at higher dietary K219 20 or K1.24
Nordic Nutrition Recommendations include a provisional recommended intake of vitamin K of 1 µg/kg body weight per day,3 while adequate intake is 90 µg/day for women and 120 µg/day for men.25 However, these recommendations may not be sufficient to attain complete carboxylation of extrahepatic vitamin K-dependent proteins.26 27
Given the limited number of epidemiological studies,18–24 and the fact that dietary vitamin K sources and content differ between countries,28–31 further research is warranted. The purpose of the current study was to evaluate the association between intake of both K1 and K2 and subsequent CHD events among community-living middle-age adults in Norway.
Subjects and methods
Study population
The current study is a prospective, community-based cohort study of participants living in Hordaland County, Norway (known as The Hordaland Health Study (HUSK); https://husk.w.uib.no/). The recruitment was based on a cohort from 1992 - 1993 (The Hordaland Homocysteine Study), where eligible subjects (Hordaland County residents born 1950–1951) were identified from the National Population Register on 31 December 1992.32 33 In 1997–1999, all living Homocysteine Study cohort members born 1950–1951 and residing in the city of Bergen or the neighbouring suburban municipalities were invited to participate in HUSK. The baseline examinations were conducted during 1997–1999 as a collaboration between the National Health Screening Service (now The Norwegian Institute of Public Health), The University of Bergen and local health services. Participation rate was 77%. Participants underwent a brief health examination and provided a non-fasting blood sample. Information on lifestyle was collected via self-administered questionnaires. A semiquantitative Food Frequency Questionnaire (FFQ) was completed by 87% of the participants yielding 3107 men and women age 47–49 years eligible for the current study.
We excluded from the analyses 27 men and 35 women who reported extreme energy intakes (below the first percentile: <1125 kcal for men and <705 kcal for women; or above the 99th percentile: >4519 kcal for men and >3571 kcal for women). Further, we excluded 27 participants (22 men and 5 women) who had prior CHD based on self-reported information and/or prior CHD hospitalisations during 1994–1999. Additionally, those with missing information on self-reported myocardial infarction from the Homocysteine Study (1992–1993; 4 men and 19 women) were excluded. Further, we excluded two participants (one man and one woman) who reported use of warfarin and six participants (two men and four women) with missing measurement on dietary vitamin K intake. The final study population thus included 1279 men and 1708 women.
Patient and public involvement
Participants were not involved in designing the research question, conducting the study, or in the interpretation, or writing of the results. There are no plans to involve participants or relevant patient communities in dissemination of results. Results are disseminated to study participants via website (https://husk.w.uib.no).
Dietary assessment
Information on food intake was obtained at baseline (1997–1999) using a slightly modified version of a previously described34 past-year 169-item semiquantitative FFQ. The FFQ was handed out on the health examination day, filled out at home and returned by mail to the HUSK project centre. The questionnaire included frequency alternatives (from once a month to several times per day), the number of units consumed and portion sizes (eg, slices, glasses, spoons) to capture the habitual diet during the past year. The dietary information presented includes individual food or beverage items, food groups and nutrient intakes. Daily nutrient intakes were computed from a database and software system developed at the Department of Nutrition, University of Oslo (KBS, V.3.2). The nutrient database is primarily based on the official Norwegian food composition table35 and available literature.28 Data for K1 are mostly developed by public authorities in Finland,36 Sweden37 and USA.38 For some Norwegian food products, analyses were performed using high performance liquid chromatography of fermented foods.39 K2 was evaluated as one entity with no distinction between the different menaquinones. At time of study, the most commonly used dietary supplements in Norway did not contain vitamin K. Thus, vitamin K intake reflects only dietary sources. Measurements used as independent variables in this study are the total dietary amount of K1 and K2, expressed as energy-adjusted residuals.40
Health examination and health habits
Baseline examinations included measurements of height, weight, waist circumference, resting blood pressure (Dinamap 845 XT equipment (Criticon)) and non-fasting venous blood samples for evaluation of serum lipids and glucose. Serum samples of total cholesterol, HDL cholesterol, triglycerides and glucose were analysed within 7 days at the department of Clinical Chemistry, Ullevål University Hospital, Oslo, using enzymatic methods with reagents from Boehringer Mannheim (Roche, Basel, Switzerland).
Information on educational level and medication use was collected through self-administered questionnaires.
Hypertension was considered present if the mean of at least two consecutive measurements of systolic blood pressure was ≥140 mm Hg or of diastolic blood pressure ≥90 mm Hg or if use of medication for hypertension was reported.
Diabetes was diagnosed according to diagnostic criteria at the time of the screening/survey. Participants taking diabetic medications or who reported a diagnosis of diabetes were defined as having diabetes mellitus. Also, participants with a serum glucose level >7 mmol/L who had not eaten a meal during the last 8 hours, or with glucose level >11.1 mmol/L and less than 8 hours since their last meal, were defined as having diabetes. Pre-diabetes was defined as having glucose levels between 5.6 and 7 mmol/L at least 8 hours after their last meal or between 7.8 and 11 mmol/L less than 8 hours after their last meal.
Participants answered one categorical question on past-year vigorous physical activity resulting in sweating or breathlessness (none, <1 h/week, 1–2 h/week, or ≥3 h/week). This variable was treated as a categorical variable with none as the reference.
Participants were classified as non-smokers, former smokers or current smokers and this variable was treated as a categorical variable with non-smokers as reference.
Outcome
The study endpoints were incident (first time) hospitalisation with CHD (International classification of diseases (ICD)9 codes 410–414, ICD10 codes I20–I25) as primary or secondary diagnosis or death with CHD as the underlying cause of death. Participants were followed from baseline through 31 December 2009 for CHD events through the Cardiovascular Disease in Norway project database (CVDNOR, www.cvdnor.no)41 42 and The Cause of Death Registry. Follow-up time represented time from baseline (1997–1999) until CHD, death from other causes, emigration or 31 December 2009, whichever came first. During follow-up, there were 107 non-fatal and 5 fatal events of interest while 60 participants died due to other causes and were censored at date of death.
Statistical analyses
Energy-adjusted residuals were obtained from linear regression models with total energy intake as independent variable and K1 or K2 as dependent variables. The residuals measure the difference between actual intake and expected intake predicted by total energy intake40 and thereby provides an assessment of K1 and K2 intake relative to energy consumed. Residuals were then categorised into sex-specific quartiles.
Descriptive characteristics included counts with percents and medians (interquartile range) for categorical and continuous variables, respectively. Trends in dichotomous, categorical and continuous baseline characteristics across energy-adjusted quartiles of K1 and K2 were evaluated using logistic, ordinal logistic and linear regression analyses, respectively. The median residuals for each quartile group was specified as a continuous independent variable in the regression models.
Cox proportional hazards models were used to calculate adjusted HR and 95% CIs for CHD associated with sex-specific energy-adjusted quartiles and per 10 µg increments of K1 or K2 intake. The covariates included were either those associated with intake of K1 or K2 and with CHD, or those that modified the association of either K1 or K2 with CHD when included in the multivariable model. Analyses included adjustments for sex, age (years) and total energy intake (kcal/day) (model 1), with additional adjustment for categories of vigorous physical activity (none vs <1 hour/week, 1–2 hours/week and ≥3 hours/week), smoking habits (previous smokers and current smokers, respectively vs non-smokers) and education (high school or vocational school, and any college or university, respectively vs primary school(≤10 years)) (model 2). In a third model, K1 was adjusted additionally for energy-adjusted intake of fibre (g/day) and folate (mg/day); K2 was additionally adjusted for energy-adjusted intake of saturated fatty acids (SFA)(g/day) and calcium (mg/day) (model 3). The following additional potential confounders were also evaluated but not included in the tables as they did not noticeably alter the vitamin K1 or K2 coefficients for CHD: family history of myocardial infarction and energy-adjusted alcohol intake (g/day). Further, adjusting for the following intermediate factors: body mass index (BMI, kg/m2), diabetes mellitus (pre-diabetes and diabetes, respectively, vs no diabetes), hypertension, serum total cholesterol (mmol/L) and statin use, only attenuated the association to a small degree.
To test for linear trends across energy-adjusted quartiles of K1 and K2 intakes, the median value of the residuals within each quartile group was entered as a continuous independent variable. Supplementary analyses re-evaluated K1 and K2 intake as sex-specific quartiles of absolute intake rather than energy-adjusted residuals.
Missing data on physical activity (3.8%), education (0.8%) and smoking habits (2.1%) were handled with listwise deletion in all analyses included in the main manuscript. In supplementary analyses, missing values for physical activity, smoking and education were imputed using ordinal logistic regression as the imputation model in MICE (multiple imputation using chained equation) with 20 imputations. All variables in the Cox regression models were included as imputation variables together with total cholesterol, HDL cholesterol, triglycerides and BMI as auxiliary variables due to their correlation with physical activity, smoking and education.
The proportional hazards assumption was evaluated using Schoenfeld’s test and log–log test.
In Cox regression with penalised splines, the functional form of the association between absolute K2 intake (not residuals) and risk of CHD was estimated by smoothing splines, in which the estimated smooth functions were used to plot the relative hazards of CHD.43 Intakes above the 95th percentile and below the 5th percentile are excluded in the figure.
To test for sex interactions between K1 and K2, we compared models with and without an interaction term using likelihood-ratio test.
Statistical analyses were performed using Stata V.15 (Stata Corp LP) and R V.3.4.0 (https://www.r-project.org/, The R Foundation for Statistical Computing, Vienna, Austria). P<0.05 were considered statistically significant.
Consent to participate
All subjects gave their written consent to participate in the study.
Results
Reported intake of energy-adjusted K1 ranged from 8 to 1063 µg/day/1000 kcal (median 48 µg/day/1000 kcal) and were higher for women compared with men. Intake of energy-adjusted K2 ranged from 1 to 31 µg/day/1000 kcal (median 7 µg/day/1000 kcal) and were slightly higher for women compared with men.
The major dietary sources of K1 were vegetables (64%), fruits and berries (6%) and milk and milk products (6%), while sources of K2 were cheese (40%), other dairy products (14%), meat (24%) and eggs (13%).
In the evaluation of baseline characteristics associated with K1 intake, the concentration of HDL cholesterol and the proportion of participants who were highly educated and reported at least 1 hour of vigorous physical activity per week were higher with higher quartiles of K1 intake (table 1). Further, intakes of energy-adjusted total vitamin K, folate and fibre were higher with higher K1 intake quartiles. In contrast, a lower proportion with a family history of CHD and lower energy-adjusted K2, SFA and carbohydrate intakes were noted with higher K1 intake. In addition, intake of fruit and berries and vegetables were higher with higher K1 intake quartiles, while intake of cheese, milk and milk products and soft drinks with sugar were lower with higher intake quartiles of energy-adjusted K1 intake.
Evaluation of baseline characteristics by quartiles of energy-adjusted K2 intake identified that the proportion of participants highly educated, and the concentration of HDL cholesterol were higher with higher quartiles of K2 intake, while the concentration of triglycerides was lower (table 2). Further, intake of energy-adjusted total fat, SFA and calcium was higher with higher K2 intake quartiles. In contrast, lower intake of energy-adjusted K1 and carbohydrates was noted with higher K2 intake. In addition, intake of butter, eggs, cheese, meat and minced meat were higher with higher quartiles, while intake of soft drinks with sugar and fruit and berries were lower with higher quartiles of energy-adjusted K2 intake.
Association between dietary vitamin K1 and CHD
During a mean 10.8 (SD 1.3) years follow-up, representing 32 362 person years among 2987 participants, we documented 112 incident CHD events. Due to listwise deletion of missing values (2.1% for smoking habits, 0.8% for education and 3.8% for physical activity), multivariable-adjusted analyses included 6.5% fewer participants compared with model 1 analyses (ie, 2792 (1213 men and 1579 women) participants and 100 CHD events).
When adjusting for age, sex and total energy intake, there was no association between intake of energy-adjusted K1 and CHD comparing the fourth to the first quartile and there was no trend (table 3, model 1). The results were similar when further adjusting for physical activity, smoking habits and education (HRQ4vs.Q1=0.92 (0.54 to 1.57), p for trend 0.64; table 3, model 2). In analyses of energy-adjusted K1 intake as a continuous variable (per 10 µg increase), there was no association between K1 and CHD in the adjusted analysis (table 3, model 2). Additional adjustments for energy-adjusted fibre and folate did not materially change the results (table 3, model 3). In supplementary analyses, where missing data were handled with multiple imputation, results were similar to those presented in table 3 (online supplementary table 1, models 2 and 3).
Results were consistent with the above analyses in the supplemental analyses evaluating sex-specific quartiles of absolute K1 intake rather than energy-adjusted residuals (online supplementary table 2).
Association between dietary vitamin K2 and CHD
When adjusting for age, sex and total energy intake, there was a lower risk of CHD with energy-adjusted K2 in the fourth compared with the first quartile (HRQ4vs.Q1=0.50 (0.28 to 0.88), p for trend 0.02; table 3, model 1). Results were consistent when further adjusting for physical activity, smoking habits and education (HRQ4vs.Q1=0.52 (0.29 to 0.94), p for trend 0.03; table 3, model 2). Consistency in results was observed in analyses of K2 intake as a continuous variable (per 10 µg increase; HR=0.74 (0.52 to 1.05), p=0.09). Additional adjustments for energy-adjusted SFA and calcium slightly attenuated the risk estimates for the association between K2 intake and CHD (HRQ4vs.Q1=0.58 (0.28 to 1.19), p for trend 0.16; table 3, model 3). Similar results were found in supplementary analyses where missing data were handled with multiple imputation (model 2: HR per 10µg increase=0.70 (0.50 to 0.98), p=0.04) (online supplementary table 1, models 2 and 3).
When evaluating sex-specific absolute K2 intake rather than energy-adjusted residuals, HRs were similar to those observed in the primary analyses (HRQ4vs.Q1=0.72 (0.36 to 1.45), p for trend 0.25; online supplementary table 2, model 2). Similarly, the penalised spline figure for absolute K2 intake and its association with CHD adjusting for model 2 covariates showed a tendency towards lower risk of CHD with higher K2 intake (figure 1).
Discussion
Among community-dwelling middle-age adults in Western Norway, a higher energy-adjusted reported intake of K2 was associated with a lower risk of subsequent CHD events, whereas intake of K1 was not associated with incident CHD. Similar direction of associations was observed when further adjusting for potential dietary confounders.
Strengths and weaknesses
Strengths of our study include a long follow-up time with minimal competing risk from other causes of death in this relatively young study population. Linkage to the CVDNOR project database assured complete follow-up. In addition, we had information on several possible confounders including baseline health status, medication use, health habits and history of CHD which enabled us to evaluate incident CHD.
Weaknesses include self-reported information on dietary intake, health habits and medication use which may lead to misclassification in covariates used in the analyses. Furthermore, we lack information on changes in diet including intake of vitamin K, medications and risk factors over time. Inherent problems with FFQs are systematic under-reporting and over-reporting. However, the FFQ is well suited to rank individuals when adjusting for total energy intake.44 This FFQ is not validated specifically for intake of vitamin K. Vitamin K content of foods differs according to production conditions, and the bioavailability is dependent on preparation, fat content of meals, the food matrix and subtypes of vitamin K. Therefore, although the rank ordering of participants may be valid, the absolute vitamin K intake based on FFQ is likely inaccurate.5 29 30 45 46 In addition, K2 can be produced by intestinal gut microbiota, but little is known about its contribution to vitamin K status since the majority is located in bacterial membranes in the colon and is probably not available for absorption.26
Further, vitamin K2 intake may be underestimated in this study population since much of the information on vitamin K content of food used in this study comes from a Dutch study,28 and recent research have shown that Norwegian cheeses are especially rich in vitamin K2.29 47 In addition, we could not differentiate the subtypes of vitamin K219 in our study.
Although we performed multivariable analyses, residual confounding may still be present.
Results in relation to other studies
This study did not find an association between intake of K1 and CHD, in line with results from most previous studies.19 20 22 23 However, Erkkilä et al concluded that high K1 intake may be a marker of low CHD risk, but that dietary patterns associated with K1 intakes, rather than intake of K1 itself might account for this association.21 Similarly, Juanola-Falgarona et al studied the association between dietary K1 and K2 with mortality in a cohort with high cardiovascular disease (CVD) risk, and found that an increase in dietary intakes of both K1 and K2 were associated with a reduced risk of all-cause mortality, while only K1 was associated with a reduced risk of CVD mortality.24 However, since participants came from a Mediterranean country in which the consumption of fruit, vegetables and vegetable oils was quite high, K1 could be regarded as a marker of adherence to a healthy diet.24
Regarding vitamin K2, Geleijnse et al found a 41% lower risk of CHD comparing the highest with the lowest tertile of K2 intake in the Rotterdam study,20 similar to our results. Gast et al found a dose–response relationship with a 9% lower CHD risk with each 10 µg higher K2 intake, with the strongest association shown for long-chain K2.19 Opposite from this study where the association attenuated when further adjusting for SFA and calcium, the association became stronger in both these studies in their multivariable models.19 20 Adjusting for SFA and calcium, however, may be an overadjustment as dairy products rich in SFA and calcium are also major sources of K2.24 Zwakenberg et al found that only intake of long-chain K2 (per10µg) was borderline significantly inversely associated with CHD mortality (p for trend 0.06).23 This discrepancy compared with our study may be because lifestyle factors such as diet may have a larger impact on total CHD events (fatal and non-fatal) given the underlying mechanisms of reducing calcification. Treatment, however, is probably of larger importance for CHD mortality.48
In the western diet, cheese is the most important source of long-chain K2, and hard cheese is, in general, richer in K2 than soft cheese.28 29 Fu et al showed that vitamin K concentrations in cheese ranged from 40 µg to 850 µg per 100 g, and that reduced-fat products contained 5%–22% of the vitamin K found in full-fat equivalents.30 K2 in cheese originates from bacterial processes present at the start of the cheese-making process.26 As different lactic acid bacteria are used in cheese making, a large variability in K2 content is found. Cheeses from Norway are among those with the highest long-chain K2 content.29 47
A meta-analysis including only two published studies of sufficient quality did not conclude that there is a lower risk of cardiovascular events with higher intake of K2.49 However, a systematic review and meta-analysis found that supplementation with vitamin K (K1 and K2) significantly reduced vascular calcification, but not vascular stiffness, compared with controls.50 However, Shea et al studied supplementation with K1 on coronary artery calcification progression in older men and women and found no difference between the control and treatment groups in the main analyses. Less progression of coronary artery calcification was found in participants who were ≥85% adherent to supplementation and in those with pre-existing coronary artery calcification.51
Potential mechanisms
The lower risk of CHD with a high intake of K2 may have different explanations. Jakobsen et al showed that intake of SFA seems to be preferable compared with intake of carbohydrates with high glycaemic index in order to reduce risk of myocardial infarction.52 Further, intake of K2 correlates positively with intake of SFA, especially dairy sources as cheese, and negatively with intake of carbohydrates, especially sugar-rich sources as soft drinks with sugar and fruit and berries (both fresh and canned). Further, cheese has been associated with lower risk of CHD,53 and the median intake of cheese more than triples between the first and fourth quartile of K2 intake. This may also partly explain why adjusting for SFA and calcium attenuated the association between intake of K2 and CHD (table 3, model 3).
The lower risk of CHD with a high intake of K2 may further be explained by carboxylation of vascular Matrix Gla Protein and consequently less arterial calcification.54 Intimal calcification starts in the inner layer of large arteries, is associated with dyslipidaemia and may cause ischemia and arterial infarction, while calcification of the medial layer occurs even in small arteries and may lead to arterial stiffness, hypertension and left ventricular hypertrophy that further increases risk of CHD.55
The observed association for K2 only may be due to the fact that in addition to being cleared by the liver, it is also transported to extrahepatic tissues.5 However, extrahepatic vitamin K-dependent proteins seem to be of lower priority compared with those in the liver.26 27 Thus, one hypothesis is that intake of vitamin K has to be of a certain magnitude in order to have an effect on CHD. The different results on K1 and K2 may also be due to biological differences between K1 and K2 or to lower ability of the FFQ to estimate K1.4 26 56 Due to different bioavailability, the contribution of K2 to vitamin K status is at least equal to that of K1 even though dietary K1 contributes to the majority of the total vitamin K intake.26 46 57
Alternatively, our findings may reflect that K2 may be a marker of another nutrient or food constituent that has heart-healthy properties.
Implications and future research
Our findings contribute to the sparse literature relating dietary vitamin K to future CHD risk. Current dietary guidelines are based on insufficient knowledge with regard to vitamin K metabolism and the different characteristics of K1 and K2. Therefore, our results indicate a need for more studies on the association between K2 and CHD. In addition, more knowledge about the absorption, transport and bioactivity of K2 is warranted.
Conclusion
In this Norwegian community-based study population, we observed that intake of K2 was associated with lower risk of CHD, while there was no association between intake of K1 and CHD. These results are considered generalisable to other middle-aged Western populations in which dairy products are the primary source of K2.