Next Article in Journal
Hydrogen Breath Tests: Are They Really Useful in the Nutritional Management of Digestive Disease?
Next Article in Special Issue
High-Protein, Low-Glycaemic Meal Replacement Decreases Fasting Insulin and Inflammation Markers—A 12-Month Subanalysis of the ACOORH Trial
Previous Article in Journal
Cloudy Apple Juice Fermented by Lactobacillus Prevents Obesity via Modulating Gut Microbiota and Protecting Intestinal Tract Health
Previous Article in Special Issue
The Mechanisms of the Anti-Inflammatory and Anti-Apoptotic Effects of Omega-3 Polyunsaturated Fatty Acids during Methotrexate-Induced Intestinal Damage in Cell Line and in a Rat Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Macronutrient Quality and All-Cause Mortality in the SUN Cohort

by
Susana Santiago
1,
Itziar Zazpe
1,
Cesar I. Fernandez-Lazaro
2,
Víctor de la O
2,
Maira Bes-Rastrollo
2 and
Miguel Ángel Martínez-González
2,*
1
Department of Nutrition and Food Sciences and Physiology, Campus Universitario, School of Pharmacy and Nutrition, University of Navarra, 31008 Pamplona, Spain
2
Department of Preventive Medicine and Public Health, Campus Universitario, School of Medicine, University of Navarra, 31008 Pamplona, Spain
*
Author to whom correspondence should be addressed.
Submission received: 2 February 2021 / Revised: 12 March 2021 / Accepted: 15 March 2021 / Published: 17 March 2021

Abstract

:
No previous study has assessed the relationship between overall macronutrient quality and all-cause mortality. We aimed to prospectively examine the association between a multidimensional macronutrient quality index (MQI) and all-cause mortality in the SUN (Seguimiento Universidad de Navarra) (University of Navarra Follow-Up) study, a Mediterranean cohort of middle-aged adults. Dietary intake information was obtained from a validated 136-item semi-quantitative food-frequency questionnaire. We calculated the MQI (categorized in quartiles) based on three quality indexes: the carbohydrate quality index (CQI), the fat quality index (FQI), and the healthy plate protein source quality index (HPPQI). Among 19,083 participants (mean age 38.4, 59.9% female), 440 deaths from all causes were observed during a median follow-up of 12.2 years (IQR, 8.3–14.9). No significant association was found between the MQI and mortality risk with multivariable-adjusted hazard ratio (HR) for the highest vs. the lowest quartile of 0.79 (95% CI, 0.59–1.06; Ptrend = 0.199). The CQI was the only component of the MQI associated with mortality showing a significant inverse relationship, with HR between extreme quartiles of 0.64 (95% CI, 0.45–0.90; Ptrend = 0.021). In this Mediterranean cohort, a new and multidimensional MQI defined a priori was not associated with all-cause mortality. Among its three sub-indexes, only the CQI showed a significant inverse relationship with the risk of all-cause mortality.

1. Introduction

Healthy eating results essential to reduce the risk of developing non-communicable diseases and increases the number of years lived with good health [1]. However, the potential impact of diet quality on chronic diseases and mortality is frequently underestimated. Findings of the Global Burden of Disease Study suggested that dietary risks were responsible for approximately 22% of all adult deaths (11 million deaths), which calls for an action to improve overall dietary quality intake [2].
Although diet quality is a complex concept that is difficult to measure, general recommendations for a healthy eating include fostering bioactive-rich foods (vegetables, fruits, seeds, whole grains, nuts, plant oils, fish, and yogurt), selecting non-processed or minimally-processed foods and avoiding ultra-processed products, which are important sources of sugars, refined starch, and industrial additives such as trans fatty acids and sodium [3]. In this context, multiple dietary indices have gained attention to measure diet quality. These indices are used in nutritional epidemiology to assess compliance with national nutrition guidelines, adherence to predefined high-quality or healthy dietary patterns, dietary changes and risk of developing chronic diseases [4]. Diet quality indices usually include items related to nutrients, foods/food groups or, most frequently, a combination of both. Diet quality indices differ in the number of foods and nutrients included, the cut-off values, and the scoring algorithms used to build them [5].
In general, an adequate balance of macronutrients and micronutrients is a common dimension of diet quality, which has been recently included in several emerging indices. Specifically, the proportionality of the energy-yielding macronutrients (proteins, carbohydrates, and fats) and/or fatty acids (monounsaturated, polyunsaturated, saturated, and trans fatty acids) (MUFA, PUFA, SFA and TFA) has been mainly addressed by preferred intake ranges and intake ratios for fatty acids [6]. Additionally, nutritional dietary guidelines across countries have proposed ranges of energy macronutrient distribution to improve overall health and prevent chronic diseases [7]. However, optimal macronutrient composition for supporting longevity remains uncertain [8,9]. In fact, some authors [10] have investigated global associations between macronutrient supply and age-specific mortality, concluding that the macronutrient composition of energy supply that minimizes mortality varies with age.
Regarding energy from carbohydrates, both high and low percentages, have been associated with increased mortality, observing the lowest risk of mortality between 50–55% of calories from carbohydrate intake [11]. In the PURE Study, the authors found that energy intake from carbohydrates higher than 60% was associated with total mortality risk, whereas no association was found for total fat and types of fat [12]. Growing evidence on dietary fats, particularly saturated fats, in relation to health and mortality is also conflicting [13,14].
Evidence from prospective cohort studies suggests that animal protein intake is associated with increased all-cause mortality, as opposed to plant protein, whereas evidence of the effect of total protein intake is more inconsistent [15,16,17].
Beyond intake or percentage of total energy intake, the potential association of macronutrients with mortality ultimately relies on overall dietary quality and specific food sources [9,18]. However, quality of macronutrients, primarily carbohydrate quality intake, has been investigated in relation to mortality [19], chronic diseases [20,21], and inflammatory markers [22]. To the best of our knowledge, to date, there is no study that has specifically investigated the association between global macronutrient quality intake and all-cause mortality. Hence, we aimed to assess the association of global macronutrient quality, through a multidimensional macronutrient quality index (MQI) with all-cause mortality in the “Seguimiento Universidad de Navarra (SUN)” (University of Navarra Follow-Up) Project.

2. Materials and Methods

2.1. Study Design

The SUN Project is a Spanish dynamic, multipurpose, prospective, and permanently open cohort of university graduates that began the recruitment of participants in 1999. Its design was based on other American large cohort studies. More detailed information about the Sun cohort has been previously described elsewhere [23]. Its primary aim was to assess the impact of lifestyle and diet on non-communicable diseases. Self-reported questionnaires at baseline and every 2 years are used to collect information mainly through mailed or electronically. The baseline questionnaire collects information related to socio-demographics variables, lifestyle, medical history, and anthropometric variables.
By December 2019, a total of 22,894 participants were recruited. For this study, individuals recruited after March 2017 (n = 341) were excluded to ensure a minimum follow-up of 2 years. We additionally excluded 1478 individuals without follow-up (retention rate 90.6%), and 1992 participants with energy intake outside of predefined limits (<500 or >3500 kcal/d for women, and <800 or >4000 kcal/d for men) [24] (Figure 1). Therefore, our analysis included a total of 19,083 subjects.
At entry to the SUN project, participants receive detailed written information and give their permission to participate in the study before any follow-up on their medical history. The voluntary completion of questionnaire at baseline implies informed consent. According to the Declaration of Helsinki, potential candidates are informed of their right to drop out of the SUN project or withdraw their consent to continue in the study at any time without reprisal. The Institutional Review Board of the University of Navarra approved the SUN project, which is registered at clinicaltrials.gov (accessed on 2 February 2021) as NCT02669602.

2.2. Exposure Assessment

Participants provided their dietary intake information at the beginning of the study (baseline) and after 10 years of follow-up. A validated self-administrated semi-quantitative 136-item food-frequency questionnaire (FFQ) was used to assess participants ‘dietary intake [25,26,27]. The FFQ classifies food items into general food groups: dairy, eggs, meat and fish, vegetables, fruits, legumes and cereals, oils and fats, pastries, beverages, and miscellaneous. FFQ items included a standard portion size and 9 consumption frequency options, from “never/almost never” to “more than 6 times/day”. Daily consumption of every food item was estimated by multiplying its typical portion size with its frequency of consumption using an ad hoc computer program specifically developed for this aim. The calculation of dietary intakes was performed by a team of trained nutritionist using the daily intake of every food item and the corresponding nutrient composition based on Spanish Food Composition Tables [28,29]. Missing information on any dietary item was considered as having no consumption for such items.
The MQI was calculated from dietary intake collected in the FFQ and was defined summing up three sub-indexes: the carbohydrate quality index (CQI), the fat quality index (FQI), and the healthy plate protein source quality index (HPPQI), as follows, MQI = CQI + FQI + HPPQI.
In brief, the CQI is based on 4 carbohydrate quality domains: glycemic index (GI), total dietary fiber intake (g/d), ratio of carbohydrates from solid/total (solid + liquid) carbohydrates, and ratio of carbohydrates from whole grains/carbohydrates from total grains (whole grains + refined grains or their products). Liquid carbohydrates were calculated summing up alcoholic and sugar-sweetened beverages and fruit juices consumption (4 items in the FFQ), while solid carbohydrate corresponded to the carbohydrate content of other foods with any carbohydrate content. On the other hand, refined grains were derived from the 13 items in the validated FFQ related to the consumption of white rice, refined bread, refined pasta, and different bakery products made with refined grains, and finally, whole grain consumption corresponded to the item “whole bread consumption” (serving size 60 g, 9 categories in the validated FFQ). The CQI was has been used in previous studies in the SUN Project, in the PREDIMED (Prevención con Dieta Mediterránea) and in the PREDIMED-Plus trials [30,31,32,33,34]. For each of its four components, we categorized participants into quintiles and received a value (ranging from 1 to 5) according to each quintile (only for GI the values were in reverse). Finally, we constructed the CQI summing for each participant the four values. The metabolomic signature of this dietary quality index has been specifically evaluated in the PREDIMED study [35].
To calculate the FQI as a continuous variable, we used the following ratio:
FQI = (Monounsaturated + Polyunsaturated)/(Saturated + Trans Fatty Acids) (MUFA + PUFA)/(SFA + TFA).
These four components were given the same weighting in the FQI. This dietary quality index has been previously used to investigate the association with nutritional adequacy in this cohort [31] and also in the PREDIMED study [33] and also with the risk of CVD [36] and with telomere length [37], among other outcomes.
Finally, the novel HPPQI has been created based on the nutritional protein quality recommendations found in the last international dietary guidelines [38,39,40,41]. This new dietary quality index was calculated with the following ratio:
HPPQI = (seafood + poultry + pulses + nuts)/(red and processed meats + cheese).
Thus, the foods located in the numerator and denominator reflect healthy and unhealthy source of proteins, respectively, according to the Harvard’s Healthy Eating Plate [41].
Lastly, to calculate each criterion of the MQI, the weighting of each criteria, CQI, FQI, and HPPQI was equal. We categorized participants into quintiles for each sub-index and we assigned values ranging from 1 point (first quintile) to 5 points (fifth quintile). Lastly, we calculated the MQI score adding the values of each sub-index, ranging from 3 (the worst quality) to 15 points (the best quality) (Table 1).

2.3. Outcome Assessment: All-Cause Mortality

The primary outcome was all-cause mortality. Information on mortality and its causes was identified through active and permanent contact with participants. Participants who did not respond to any of the 5 follow-up mailings were contacted by email or phone. We were able to identify approximately 75% of deaths through reports next-to skin, work associates, professional associations, or the authorities’ postal system. Additionally, we checked, at least once a year, both the National Statistical Institute and the Spanish National Death Index—a central index that collects information of any death registered in Spain, to confirm the status vital of all participants and to identify their cause of death if unknown. The last time when mortality and its causes were checked was December 2018. Death certificates and medical records of its causes were obtained.

2.4. Ascertainment of Covariates

When participants enter in the cohort, they provided additional information about socio-demographic characteristics (marital status and level of education, among others), anthropometric measures (weight, height, and weight change during the past 5 years), lifestyle and health-related habits (alcohol intake, smoking status, physical activity, leisure-time watching TV, lifestyle behavior, supplement use, and special diet), and family and personal medical history. Two subsamples of the SUN cohort were used to assess accuracy of physical activity [42] and of self-reported weight and height [43]. A 9-item Mediterranean Diet Score [44], frequently used in nutritional epidemiology [45] was used to assess the adherence to this dietary pattern. The prevalence and history of several non-communicable diseases was ascertained at baseline and updated until the exit of the cohort or until death was reported. Cardiovascular disease included any reported medical diagnosis of myocardial infarction, stroke, atrial fibrillation, coronary artery bypass grafting or other re-vascularization procedures, paroxysmal tachycardia, peripheral venous thrombosis, aortic aneurism, pulmonary embolism or heart failure. Self-reported cardiovascular events, cancer, and diabetes have always been confirmed by medical records [43,46,47]. These self-reports can be considered as reliable given than all participants of the SUN project are graduates, and more than 50% were health professionals themselves. Moreover, the validity of self-reported obesity, dyslipidemia, and hypertension diagnoses was evaluated in different subsamples of the SUN project.

2.5. Statistical Analyses

We described the baseline characteristics of participants adjusted for age and sex using the inverse probability weighting method, using proportions and means according to quartiles of the MQI. We use crude and multivariable Cox regressions models to assess the association of the MQI and each of its sub-index (CQI, FQI, and HPPQI) with all-cause mortality. We calculated hazard ratios (HRs) with their 95 % confidence interval (CIs) and we considered always the lowest quartile as the reference category.
All models included age (underlying time variable), and were stratified by recruitment period (5 categories) and deciles of age. Time at entry was considered as the date of completion of first questionnaire and exit time was date of death or date when participants completed the last follow-up questionnaire. We determined follow-up from the date when the questionnaire completed at baseline to the last follow-up questionnaire was returned, for whom death was not reported.
Participants completed a full-length FFQ after 10 years of follow-up allowed. We used repeated measurements using both cumulative average information and updated data of the MQI and its sub-indexes to reduce any effect of variation in dietary pattern and to obtain a more realistic measurement of long-term diet during follow-up.
After crude analyses, we fitted three multivariate Cox regression models with potential confounders that may affect the effect of CQI on mortality risk: model 1 was further adjusted for sex and total energy intake (kcal/d, continuous); model 2 was additionally adjusted for marital status (single, married, widowed, and others), educational level (years of higher education, continuous), smoking (never, current, and former smoker), cumulative smoking habit (packs/year, continuous), alcohol intake (never, <5 women or <10 men g/d, 5–25 women or 10–50 men g/d, and >25 women or >50 men g/d), leisure-time physical activity (metabolic equivalent-h/week, continuous), body mass index (BMI [kg/m2, linear and quadratic terms, continuous]), time spent sitting (hours/week, continuous), weight gain in the previous 5 years before entering the cohort (<3 kg and ≥3 kg), and following a special diet at baseline (yes/no); and model 3 was further adjusted for family history of CVD (yes/no), and the following prevalent diseases (yes/no): diabetes, CVD, cancer, depression, hypertension, hypercholesterolemia, and hypertriglyceridemia.
The individual contribution of each specific sub-index was also evaluated. Moreover, when we calculated HR for each component of MQI, we fitted an additional analysis adjusting for the intake of the rest of macronutrients in the fully adjusted model.
Potential confounders were selected a priori based on prior knowledge and previous findings of the SUN cohort on all-cause mortality [19,48]. We assessed linear trends across increasing categories by assigning medians to each quartile, and this variable was treated as continuous.
We assessed the interaction between age (continuous) and the MQI (continuous and categorized in quartiles) by testing an interaction product-term with the maximum likelihood ratio test.
The following sensitivity analyses were performed: (a) selection of only male or only female participants, (b) selection of participants aged <45 or ≥45 years, (c) censoring the analysis in ≥50 years, (d) selection of only health or non-health professionals, (d) excluding subjects with prevalent hypertension or hypercholesterolemia, (e) definition of a different predefined limits of energy intake (percentiles 5 and 95), (f) excluding subjects with prevalent CVD or cancer, (g) excluding subjects who were following a special diet at baseline, (h) excluding subjects with ≥30 missing items in FFQ, (i) redefining the MQI used the following ratio to calculate PQIa = (fish + poultry and lean meats + eggs + pulses + nuts + reduced-fat dairy + whole grains)/(red and processed meats + full-fat dairy + refined cereals), (j) redefining the MQI used the following ratio to calculate PQIb = (seafood + poultry + pulses + nuts + dairy products)/(red and processed meats) and finally (k) changing the weights of the 3 components of the MQI as follows, CQI (55%), FQI (30%) and HPPQI (15%) according to an acceptable macronutrient distribution [40].
We used in all analyses STATA version 16 (STATA Corp., College Station, TX, USA) with the SUN database updated in December 2019. All values presented are two-tailed; p-value < 0.05 was deemed as statistically significant.

3. Results

The present analysis included 19,083 Spanish adults (mean [SD] age 38.4 [12.4] years; 7654 men [40.5%]) that were followed up for a median time of 12.2 years (interquartile range 8.3 to 14.9). During this period, we identified 440 deaths from all causes, including 84 deaths due to CVD, 226 due to cancer, and 130 due to other causes. Among deceased participants, the mean age (SD) at death was 57.8 (15.5) years.
Table 2 shows age- and sex-adjusted baseline characteristics of participants according to quartiles of MQI. The median value of the MQI was 9 (range 3 to 15). Participants with higher MQI scores (range 12 to 15) were more likely to be single, never smokers, physically active, follow a special diet, and had higher baseline prevalence of chronic conditions. Individuals with better MQI scores had lower energy intake from fat. Conversely, participants with lower adherence to the MQI (range 3 to 7) were more likely to be married, and to snack be between meals and had lower adherence to a Mediterranean diet.
Table 3 presents the associations between all-cause mortality associated with quartiles of the MQI. No significant association was observed between mortality and the MQI when comparing the highest vs. the lowest quartile (reference category). In model 3, the adjusted HRs across successive quartiles compared with Q1 were: 1.03 (95% CI, 0.79–1.34) for Q2; 1.15 (95 % CI, 0.89–1.50) for Q3; and 0.79 (95% CI, 0.59–1.06) for Q4, with a non-significant linear trend (Ptrend = 0.199). We did not find any interaction between age and MQI (continuous, p-interaction = 0.640 and categorized in quartiles, p-interaction = 0.451).
The relationship between quartiles of each individual dimension of the MQI and mortality was additionally assessed (Table 4). The multivariate-adjusted HR between extreme quartiles of the CQI was 0.64 (95% CI, 0.45–0.90), with a significant inverse dose-response relation (Ptrend = 0.021). However, neither FQI nor HPPQI were associated with lower all-cause mortality risk. For the FQI, the multivariable-adjusted HRs were 0.85 (95% CI, 0.64–1.13), 0.86 (95% CI, 0.65–1.14) and 0.75 (95% CI, 0.56–1.00) in successive quartiles in comparison with the lowest quartile (Ptrend = 0.070). On the other hand, for the HPPQI the HRs in the most adjusted model were: 0.80 (95% CI, 0.60–1.07), 0.74 (95% CI, 0.55–0.99) and 0.93 (95% CI, 0.69–1.25) for the Q2, Q3, and Q4 respectively (Ptrend = 0.844).
We conducted Cox regression models to calculate the risk of mortality, comparing quartiles of MQI (Table 5) and CQI (Table 6) for update and cumulative average dietary information after 10 years of follow-up. After using repeated nutritional measurements in our analyses, we found that the MQI was not associated with lower mortality when comparing extreme quartiles. Thus, the MQI was not associated with significantly different rates of all-cause mortality with HRs of 0.90 (95% CI, 0.66–1.23) when using an Updated Diet and 0.84 (95% CI, 0.63–1.12) when using Cumulative Diet Average with the information repeated after 10-year follow-up. However, the individual analysis of CQI revealed an inverse association with all-cause of mortality when we compared extreme quartiles using both approaches: multivariable-adjusted HRs 0.68 (95% CI, 0.49–0.94; Ptrend = 0.031) for updated, and 0.69 (95% CI, 0.50–0.93; Ptrend = 0.016) for cumulative average information, respectively.
We conducted multiple sensitivity analyses to assess the robustness of our findings. Overall, these analyses were consistent with our main analyses and estimates remained similar (Figure 2).

4. Discussion

In this Mediterranean population, we found no statistically significant association between a new multidimensional score of macronutrient quality, MQI, and all-cause mortality. However, higher quality of carbohydrate showed an association with lower mortality.
The role of carbohydrate intake on health continues to be a matter of debate. Concerning quantity, a previous investigation found a U-shaped association between carbohydrate intake and mortality [8]. Nevertheless, other authors have suggested that carbohydrate quality and food sources seem to play a more important role in population health than quantity of carbohydrate intake [9,20]. In the majority of the studies, carbohydrate quality has been operationally measured using unidimensional indicators, mainly dietary GI or glycemic load, dietary fiber, or whole grains intake. In fact, both higher intakes of whole grains and dietary fiber are associated with reduced incidence of mortality, whereas less useful markers of carbohydrate quality on mortality are GI/load [21,49,50]. Furthermore, higher consumption of starchy carbohydrates and sugar are associated with a higher mortality risk [51]. Possible biological explanation for this association is that higher fiber intake could reduce concentrations of serum inflammatory biomarkers, and inflammation could be associated with mortality risk [22,52]. Previous investigation in this cohort, conclude that CQI seems to reflect the combined effect and potential synergies of different dimensions of carbohydrate quality [19].
Regarding fat quality, we used the ratio MUFA + PUFA /SFA + TFA. This criterion was based on current dietary guidelines and previous evidence on the association of diets high in SFA and TFA and with higher mortality, whereas in the case of diets high in PUFA and MUFA, there was an inverse association [14]. In particular, the anti-inflammatory effect of omega-3 fatty acids could explain, at least partly, their beneficial effect on mortality. Furthermore, the pathophysiology of neurodegenerative or CVD, cancers, among other age-related morbidities, has been associated with inflammation [53]. However, in our study, FQI was not independently associated with all-cause mortality. A previous large research work including 18 countries found that higher intake of MUFA, PUFA, and SFA was associated with lower mortality risk [12]. In this sense, emerging evidence is questioning current dietary fat recommendations, especially in case of SFA and food sources of specific fatty acids, arguing that SFA can have different biologic effects, which could be further modified by the carbohydrate content of diet food matrix [13,51]. Some evidence suggests that SFA (and especially dairy-derived SFA) [54], may not predict overall mortality [13,55] and that MUFA may also be not relevant, since, just as an example, a PREDIMED cohort derived paper, comparing extra-virgin and normal olive oil association with CVD, yields opposite results, suggesting that MUFA may have no direct effect on this outcome, be just be a marker of the intake of some phytochemical (polyphenols) [56,57]. However, isocaloric substitution models [58,59] and the 2015 meta-analysis by the Cochrane Collaboration [60] support small, but potentially important, reductions in CVD and total mortality by reducing SFA intake, and replacing it with PUFA or MUFA. Replacement with carbohydrate appears less useful.
In relation to protein quality, higher HPPQI was not associated with lower mortality risk in this study. Previous investigations on the influence of protein in all-cause mortality have centered on animal or plant proteins effect, or replacement of animal sources by plant sources [15,17]. In general, protein-rich diets, particularly including red and processed meats, and also animal sources increase health risks and mortality [17,61]. On the other hand, among elderly people, higher protein intake is crucial to increase muscle mass and strength, which are both independent predictors of mortality [53,62]. Regarding protein quality intake, food content of essential amino acids and their digestibility, have been widely used as traditional indicators, but in this context, food sources of traditional high-quality protein, do not imply higher diet quality or health outcomes. For instance, certain animal sources of protein have been related to increased morbidity, whereas consumption of plant proteins and plant-based dietary patterns are associated with more health benefits [63]. As mentioned before, in our investigation, protein quality was established, taking into account food protein sources, classified into healthy or unhealthy, independently of its animal or plant origin. Thus, the healthy protein group included both animal (seafood and poultry). In relation to meat and dairy, which are sources of SFA, TFA, and protein, but also micronutrients, Astrup et al. [13], found these food groups are not associated with increased risk of CVD, and suggested that limiting their intake is not supported by recent evidence. In this sense, the complexity and diversity of the relation between isolated or global macronutrient intake, sources, and health outcomes, should be highlighted. Thus, health effects of foods may not be predicted by their content in any nutrient group without considering the overall macronutrient distribution and food matrix [51].
Several limitations should be acknowledged. First, dietary variables were self-reported through a FFQ, which could have a certain degree of measurement bias. However, this method is considered the most appropriate approach to assess food and nutrient intake in large cohorts and moreover, the FFQ used in this adult cohort, has been repeatedly validated [25,26]. Besides, participants with energy intakes outside predefined limits were excluded [24]. Second, habitual diet was assessed at baseline, and changes in diet during follow-up are possible. For this reason, we repeated the main analysis after 10-years of follow-up using two different approaches to better capture the exposure. Third, the dietary quality indexes were not formally validated as such, but we have calculated them through a FFQ repeatedly validated in Spain [25,26,27]. Fourth, our sample was not representative of the general population. Thus, the generalization of our findings should be based on plausible biological mechanisms rather than on statistical representativeness. However, this homogeneity among participants reduces the likelihood of misclassification bias, reduces potential confounding and increase internal validity. Fifth, this is cohort of university graduates, which may influence diet and lifestyles. However, both high level of education and homogeneity of the participants add validity to the self-reported data collection and prevent confounding by education and other socioeconomic factors. Sixth, residual confounding cannot be totally eliminated, even though analyses were adjusted for major potential confounders. Finally, this is a young cohort with low rate of deaths that may have affected the main analyses. Thus, our estimates for all-cause mortality could be underpowered. Seventh, we acknowledge that the MQI and its components have not been formally validated, and their building might raise some concerns about arbitrary decisions. However, we feel that such concerns are unfounded because we defined all indices a priori, meaning that we used previously peer-reviewed published literature to compute them. Therefore, the arbitrariness of indices should not be claimed. In fact, we created our indices based on the best scientific evidence available at the moment of conducting our analyses. In mainstream nutritional science, most a priori-defined dietary indexes assign equal weights to each item. For example, three a priori dietary quality indexes have used our same approach to equally weight each component: the Provegetarian Food Pattern [64], the Mediterranean Diet Score [44], and the Dietary approaches to Stop Hypertension [65].
To date, overall quality of macronutrient intake has not been prospectively assessed at the same time in a large adult cohort. Furthermore, this study includes the use of the HPPQI, a novel index to measure protein quality, and the MQI have been calculated using dietary quality indexes, the FQI and the CQI, repeatedly used in previous research. Other strengths are the relatively large sample size, its prospective design with along follow-up period, the high retention proportion, and the use of multivariate analyses, adjusted for a wide number of potential confounders. Besides, additional assets of this investigation are the use of validated questionnaires and methods, the utilization of repeated dietary measurements and the inclusion of multiple sensitivity analyses, which adds robustness to findings. It has to be also mentioned the verification of deaths by medical records or consultation of the National Death Index. Finally, to our knowledge, there are not previous studies that have examined the association of a global MQI with all-cause mortality.

5. Conclusions

In conclusion, in this Mediterranean cohort, we did not find a significant inverse relationship between a new multidimensional MQI and lower risk of all-cause mortality, but the association was significant for the CQI. These results emphasize the difficulty in establishing associations between global quality of macronutrients and mortality risk among middle-aged adults.

Author Contributions

I.Z. and M.Á.M.-G. were involved with study conception and design; I.Z., S.S., C.I.F.-L. and M.Á.M.-G. performed the data analysis and interpretation; S.S., I.Z., C.I.F.-L. drafted the manuscript; V.d.l.O., M.B.-R. and M.Á.M.-G. provided critical edits to the manuscript. All authors have revised the manuscript for important intellectual content and approved the final version to be published. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the Instituto de Salud Carlos III and European Regional Development Fund (FEDER) (RD 06/0045, CIBER-OBN, Grants PI10/02658, PI10/02293, PI13/00615, PI14/01668, PI14/01798, PI14/01764, PI17/01795, PI20/00564 and G03/140), the Navarra Regional Government (45/2011, 122/2014, 41/2016), and the University of Navarra.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki. The present study was approved by the Institutional Review Board of the University of Navarra and this cohort is registered at clinicaltrials.gov (accessed on 2 February 2021) as NCT02669602.

Informed Consent Statement

Detailed written information was provided to candidates before participation. The voluntary completion of the baseline questionnaire was considered to imply informed consent. All participants gave their permission before any follow-up on their medical history. Potential candidates were informed of their right to refuse to participate in the SUN project or to withdraw their consent to participate at any time without reprisal, according to the principles of the Declaration of Helsinki.

Data Availability Statement

The data that support the findings of this study are available from the SUN Project at [email protected], upon reasonable request.

Acknowledgments

We would like to thank all participants of the SUN project for their continued cooperation and participation. Special thanks go to all the members of the SUN project for administrative, technical, and material support. We also would like to thank the members of the Department of Nutrition of the Harvard School of Public Health (Willet WC., Hu FB., and Ascherio A.) who helped us to design the SUN study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chatelan, A.; Muka, T.; Franco, O.H. Eating to add years of life and life to these years: What should be in the menu? Am. J. Clin. Nutr. 2020, 111, 733–734. [Google Scholar] [CrossRef]
  2. Afshin, A.; Sur, P.J.; Fay, K.A.; Cornaby, L.; Ferrara, G.; Salama, J.S.; Mullany, E.C.; Abate, K.H.; Abbafati, C.; Abebe, Z.; et al. Health effects of dietary risks in 195 countries, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 2019, 393, 1958–1972. [Google Scholar] [CrossRef] [Green Version]
  3. Mozaffarian, D.; Forouhi, N.G. Dietary guidelines and health—Is nutrition science up to the task? BMJ 2018, 360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Burggraf, C.; Teuber, R.; Brosig, S.; Meier, T. Review of a priori dietary quality indices in relation to their construction criteria. Nutr. Rev. 2018, 76, 747–764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Waijers, P.M.C.M.; Feskens, E.J.M.; Ocké, M.C. A critical review of predefined diet quality scores. Br. J. Nutr. 2007, 97, 219–231. [Google Scholar] [CrossRef] [Green Version]
  6. Miller, V.; Webb, P.; Micha, R.; Mozaffarian, D. Defining diet quality: A synthesis of dietary quality metrics and their validity for the double burden of malnutrition. Lancet Planet. Health 2020, 4, e352–e370. [Google Scholar] [CrossRef]
  7. Venn, B.J. Macronutrients and human health for the 21st century. Nutrients 2020, 12, 2363. [Google Scholar] [CrossRef]
  8. Kwon, Y.-J.; Lee, H.S.; Park, J.-Y.; Lee, J.-W. Associating Intake Proportion of Carbohydrate, Fat, and Protein with All-Cause Mortality in Korean Adults. Nutrients 2020, 12, 3208. [Google Scholar] [CrossRef]
  9. Shan, Z.; Guo, Y.; Hu, F.B.; Liu, L.; Qi, Q. Association of Low-Carbohydrate and Low-Fat Diets with Mortality among US Adults. JAMA Intern. Med. 2020, 180, 513–523. [Google Scholar] [CrossRef]
  10. Senior, A.M.; Nakagawa, S.; Raubenheimer, D.; Simpson, S.J. Global associations between macronutrient supply and age-specific mortality. Proc. Natl. Acad. Sci. USA 2020, 117, 30824–30835. [Google Scholar] [CrossRef]
  11. Seidelmann, S.B.; Claggett, B.; Cheng, S.; Henglin, M.; Shah, A.; Steffen, L.M.; Folsom, A.R.; Rimm, E.B.; Willett, W.C.; Solomon, S.D. Dietary carbohydrate intake and mortality: A prospective cohort study and meta-analysis. Lancet Public Health 2018, 3, e419–e428. [Google Scholar] [CrossRef] [Green Version]
  12. Dehghan, M.; Mente, A.; Zhang, X.; Swaminathan, S.; Li, W.; Mohan, V.; Iqbal, R.; Kumar, R.; Wentzel-Viljoen, E.; Rosengren, A.; et al. Associations of fats and carbohydrate intake with cardiovascular disease and mortality in 18 countries from five continents (PURE): A prospective cohort study. Lancet 2017, 390, 2050–2062. [Google Scholar] [CrossRef] [Green Version]
  13. Astrup, A.; Magkos, F.; Bier, D.M.; Brenna, J.T.; de Oliveira Otto, M.C.; Hill, J.O.; King, J.C.; Mente, A.; Ordovas, J.M.; Volek, J.S.; et al. Saturated Fats and Health: A Reassessment and Proposal for Food-Based Recommendations: JACC State-of-the-Art Review. J. Am. Coll. Cardiol. 2020, 76, 844–857. [Google Scholar] [CrossRef]
  14. Kim, Y.; Je, Y.; Giovannucci, E.L. Association between dietary fat intake and mortality from all-causes, cardiovascular disease, and cancer: A systematic review and meta-analysis of prospective cohort studies. Clin. Nutr. 2020. [Google Scholar] [CrossRef]
  15. Qi, X.X.; Shen, P. Associations of dietary protein intake with all-cause, cardiovascular disease, and cancer mortality: A systematic review and meta-analysis of cohort studies. Nutr. Metab. Cardiovasc. Dis. 2020, 30, 1094–1105. [Google Scholar] [CrossRef] [PubMed]
  16. Naghshi, S.; Sadeghi, O.; Willett, W.C.; Esmaillzadeh, A. Dietary intake of total, animal, and plant proteins and risk of all cause, cardiovascular, and cancer mortality: Systematic review and dose-response meta-analysis of prospective cohort studies. BMJ 2020, 370. [Google Scholar] [CrossRef]
  17. Chen, Z.; Glisic, M.; Song, M.; Aliahmad, H.A.; Zhang, X.; Moumdjian, A.C.; Gonzalez-Jaramillo, V.; van der Schaft, N.; Bramer, W.M.; Ikram, M.A.; et al. Dietary protein intake and all-cause and cause-specific mortality: Results from the Rotterdam Study and a meta-analysis of prospective cohort studies. Eur. J. Epidemiol. 2020, 35, 411–429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Shan, Z.; Rehm, C.D.; Rogers, G.; Ruan, M.; Wang, D.D.; Hu, F.B.; Mozaffarian, D.; Zhang, F.F.; Bhupathiraju, S.N. Trends in Dietary Carbohydrate, Protein, and Fat Intake and Diet Quality among US Adults, 1999–2016. J. Am. Med. Assoc. 2019, 322, 1178–1187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Fernandez-Lazaro, C.I.; Zazpe, I.; Santiago, S.; Toledo, E.; Barbería-Latasa, M.; Martínez González, M.Á. Association of carbohydrate quality and all-cause mortality in the Sun Project: A prospective cohort study. Clin. Nutr. 2020. [Google Scholar] [CrossRef] [PubMed]
  20. Ludwig, D.S.; Hu, F.B.; Tappy, L.; Brand-Miller, J. Dietary carbohydrates: Role of quality and quantity in chronic disease. BMJ 2018, 361. [Google Scholar] [CrossRef] [Green Version]
  21. Reynolds, A.; Mann, J.; Cummings, J.; Winter, N.; Mete, E.; Te Morenga, L. Carbohydrate quality and human health: A series of systematic reviews and meta-analyses. Lancet 2019, 393, 434–445. [Google Scholar] [CrossRef] [Green Version]
  22. Buyken, A.E.; Goletzke, J.; Joslowski, G.; Felbick, A.; Cheng, G.; Herder, C.; Brand-Miller, J.C. Association between carbohydrate quality and inflammatory markers: Systematic review of observational and interventional studies. Am. J. Clin. Nutr. 2014, 99, 813–833. [Google Scholar] [CrossRef] [Green Version]
  23. Carlos, S.; De La Fuente-Arrillaga, C.; Bes-Rastrollo, M.; Razquin, C.; Rico-Campà, A.; Martínez-González, M.A.; Ruiz-Canela, M. Mediterranean diet and health outcomes in the SUN cohort. Nutrients 2018, 10, 439. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Willet, W.C. Nutritional Epidemiology, 3rd ed.; Oxford University Press: New York, NY, USA, 2013. [Google Scholar]
  25. De La Fuente-Arrillaga, C.; Vzquez Ruiz, Z.; Bes-Rastrollo, M.; Sampson, L.; Martinez-González, M.A. Reproducibility of an FFQ validated in Spain. Public Health Nutr. 2010, 13, 1364–1372. [Google Scholar] [CrossRef]
  26. Fernández-Ballart, J.D.; Piñol, J.L.; Zazpe, I.; Corella, D.; Carrasco, P.; Toledo, E.; Perez-Bauer, M.; Martínez-González, M.Á.; Salas-Salvadó, J.; Martn-Moreno, J.M. Relative validity of a semi-quantitative food-frequency questionnaire in an elderly Mediterranean population of Spain. Br. J. Nutr. 2010, 103, 1808–1816. [Google Scholar] [CrossRef] [Green Version]
  27. Martin-moreno, J.M.; Boyle, P.; Gorgojo, L.; Maisonneuve, P.; Fernandez-rodriguez, J.C.; Salvini, S.; Willett, W.C. Development and validation of a food frequency questionnaire in Spain. Int. J. Epidemiol. 1993, 22, 512–519. [Google Scholar] [CrossRef] [PubMed]
  28. Mataix Verdu, J. Tabla de Composición de Alimentos Españoles (Spanish Food Composition Tables), 4th ed.; Universidad de Granada: Granada, Spain, 2003. [Google Scholar]
  29. Moreiras, O.; Carbajal, Á.; Cabrera, L.; Cuadrado, C. Tablas de Composición de Alimentos (Food Composition Tables), 9th ed.; Pirámide: Madrid, Spain, 2005. [Google Scholar]
  30. Zazpe, I.; Santiago, S.; Gea, A.; Ruiz-Canela, M.; Carlos, S.; Bes-Rastrollo, M.; Martínez-González, M.A. Association between a dietary carbohydrate index and cardiovascular disease in the SUN (Seguimiento Universidad de Navarra) Project. Nutr. Metab. Cardiovasc. Dis. 2016, 26, 1048–1056. [Google Scholar] [CrossRef]
  31. Zazpe, I.; Sánchez-Taínta, A.; Santiago, S.; De La Fuente-Arrillaga, C.; Bes-Rastrollo, M.; Martínez, J.A.; Martínez-González, M.Á. Association between dietary carbohydrate intake quality and micronutrient intake adequacy in a Mediterranean cohort: The SUN (Seguimiento Universidad de Navarra) project. Br. J. Nutr. 2014, 111, 2000–2009. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Santiago, S.; Zazpe, I.; Bes-Rastrollo, M.; Sánchez-Tainta, A.; Sayón-Orea, C.; De La Fuente-Arrillaga, C.; Benito, S.; Martínez, J.A.; Martínez-González, M. Carbohydrate quality, weight change and incident obesity in a Mediterranean cohort: The SUN Project. Eur. J. Clin. Nutr. 2015, 69, 297–302. [Google Scholar] [CrossRef] [Green Version]
  33. Sánchez-Tainta, A.; Zazpe, I.; Bes-Rastrollo, M.; Salas-Salvadó, J.; Bullo, M.; Sorlí, J.V.; Corella, D.; Covas, M.I.; Arós, F.; Gutierrez-Bedmar, M.; et al. Nutritional adequacy according to carbohydrates and fat quality. Eur. J. Nutr. 2016, 55, 93–106. [Google Scholar] [CrossRef] [PubMed]
  34. Martínez-González, M.A.; Fernandez-Lazaro, C.I.; Toledo, E.; Diáz-López, A.; Corella, D.; Goday, A.; Romaguera, D.; Vioque, J.; Alonso-Gómez, Á.M.; Wärnberg, J.; et al. Carbohydrate quality changes and concurrent changes in cardiovascular risk factors: A longitudinal analysis in the PREDIMED-Plus randomized trial. Am. J. Clin. Nutr. 2020, 111, 291–306. [Google Scholar] [CrossRef] [PubMed]
  35. Bulló, M.; Papandreou, C.; Ruiz-Canela, M.; Guasch-Ferré, M.; Li, J.; Hernández-Alonso, P.; Toledo, E.; Liang, L.; Razquin, C.; Corella, D.; et al. Plasma Metabolomic Profiles of Glycemic Index, Glycemic Load, and Carbohydrate Quality Index in the PREDIMED Study. J. Nutr. 2021, 151, 50–58. [Google Scholar] [CrossRef]
  36. Santiago, S.; Zazpe, I.; Gea, A.; Nuñez-Córdoba, J.M.; Carlos, S.; Bes-Rastrollo, M.; Martínez-González, M.A. Fat Quality Index and Risk of Cardiovascular Disease in the Sun Project. J. Nutr. Health Aging 2018, 22, 526–533. [Google Scholar] [CrossRef]
  37. Ojeda-Rodríguez, A.; Zazpe, I.; Alonso-Pedrero, L.; Zalba, G.; Guillen-Grima, F.; Martinez-Gonzalez, M.A.; Marti, A. Association between diet quality indexes and the risk of short telomeres in an elderly population of the SUN project. Clin. Nutr. 2020, 39, 2487–2494. [Google Scholar] [CrossRef] [PubMed]
  38. U.S. Department of Health and Human Services; U.S. Department of Agriculture. 2015–2020 Dietary Guidelines for Americans. Available online: http://health.gov/dietaryguidelines/2015/guidelines/ (accessed on 15 December 2020).
  39. Harvard T.H. Chan School of Public. The Nutrition Source. What Should I Eat? Available online: https://www.hsph.harvard.edu/nutritionsource/what-should-you-eat/protein/ (accessed on 15 December 2020).
  40. Institute of Medicine. Dietary Reference Intakes: Dietary Reference Intakes for Energy, Carbohydrate, Fiber, Fat, Fatty Acids, Cholesterol, Protein, and Amino Acids (Macronutrients); National Academy Press: Washington, DC, USA, 2005. [Google Scholar] [CrossRef]
  41. Harvard T.H. Chan School of Public. Healthy Eating Plate | The Nutrition Source. Available online: https://www.hsph.harvard.edu/nutritionsource/healthy-eating-plat (accessed on 15 December 2020).
  42. Martínez-González, M.A.; López-Fontana, C.; Varo, J.J.; Sánchez-Villegas, A.; Martinez, J.A. Validation of the Spanish version of the physical activity questionnaire used in the Nurses’ Health Study and the Health Professionals’ Follow-Up Study. Public Health Nutr. 2005, 8, 920–927. [Google Scholar] [CrossRef] [PubMed]
  43. Bes-Rastrollo, M.; Pérez Valdivieso, J.R.; Sánchez-Villegas, A.; Alonso, Á.; Martínez-González, M. Validación del peso e índice de masa corporal auto-declarados de los participantes de una cohorte de graduados universitarios [Validation of self-reported weight and body mass index in a cohort of university graduates]. Rev. Esp. Obes. 2005, 3, 352–358. [Google Scholar]
  44. Trichopoulou, A.; Costacou, T.; Bamia, C.; Trichopoulos, D. Adherence to a Mediterranean Diet and Survival in a Greek Population. N. Engl. J. Med. 2003, 348, 2599–2608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Martínez-González, M.A.; Gea, A.; Ruiz-Canela, M. The Mediterranean Diet and Cardiovascular Health: A Critical Review. Circ. Res. 2019, 124, 779–798. [Google Scholar] [CrossRef]
  46. Alonso, A.; Beunza, J.J.; Delgado-Rodríguez, M.; Martínez-González, M.A. Validation of self reported diagnosis of hypertension in a cohort of university graduates in Spain. BMC Public Health 2005, 5, 1–7. [Google Scholar] [CrossRef] [Green Version]
  47. Fernández-Montero, A.; Beunza, J.J.; Bes-Rastrollo, M.; Barrio, M.T.; de la Fuente-Arrillaga, C.; Moreno-Galarraga, L.; Martínez-González, M.A. Validación de los componentes del síndrome metabólico autodeclarados en un estudio de cohortes. Gac. Sanit. 2011, 25, 303–307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Morales, G.; Martínez-González, M.A.; Barbería-Latasa, M.; Bes-Rastrollo, M.; Gea, A. Mediterranean diet, alcohol-drinking pattern and their combined effect on all-cause mortality: The Seguimiento Universidad de Navarra (SUN) cohort. Eur. J. Nutr. 2020. [Google Scholar] [CrossRef]
  49. Yang, Y.; Zhao, L.G.; Wu, Q.J.; Ma, X.; Xiang, Y.B. Association between dietary fiber and lower risk of all-cause mortality: A meta-analysis of cohort studies. Am. J. Epidemiol. 2015, 181, 83–91. [Google Scholar] [CrossRef] [PubMed]
  50. Clar, C.; Al-Khudairy, L.; Loveman, E.; Kelly, S.A.M.; Hartley, L.; Flowers, N.; Germanò, R.; Frost, G.; Rees, K. Low glycaemic index diets for the prevention of cardiovascular disease. Cochrane Database Syst. Rev. 2017, 7, CD004467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Ho, F.K.; Gray, S.R.; Welsh, P.; Petermann-Rocha, F.; Foster, H.; Waddell, H.; Anderson, J.; Lyall, D.; Sattar, N.; Gill, J.M.R.; et al. Associations of fat and carbohydrate intake with cardiovascular disease and mortality: Prospective cohort study of UK Biobank participants. BMJ 2020, 368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Hardy, D.S.; Garvin, J.T.; Xu, H. Carbohydrate quality, glycemic index, glycemic load and cardiometabolic risks in the US, Europe and Asia: A dose–response meta-analysis. Nutr. Metab. Cardiovasc. Dis. 2020, 30, 853–871. [Google Scholar] [CrossRef]
  53. Ekmekcioglu, C. Nutrition and longevity–From mechanisms to uncertainties. Crit. Rev. Food Sci. Nutr. 2020, 60, 3063–3082. [Google Scholar] [CrossRef]
  54. Basterra-Gortari, F.J.; Bes-Rastrollo, M.; Forga, L.; Martínez, J.A.; Martínez-González, M.A. Validity of self-reported body mass index in the National Health Survey. Anales Sist. Sanit. Navar. 2007, 30, 373–381. [Google Scholar]
  55. Harcombe, Z. US dietary guidelines: Is saturated fat a nutrient of concern? Br. J. Sports Med. 2019, 53, 1393–1396. [Google Scholar] [CrossRef]
  56. Billingsley, H.E.; Carbone, S. The antioxidant potential of the Mediterranean diet in patients at high cardiovascular risk: An in-depth review of the PREDIMED. Nutr. Diabetes 2018, 8, 13. [Google Scholar] [CrossRef] [Green Version]
  57. Estruch, R.; Ros, E.; Salas-Salvadó, J.; Covas, M.-I.; Corella, D.; Arós, F.; Gómez-Gracia, E.; Ruiz-Gutiérrez, V.; Fiol, M.; Lapetra, J.; et al. Primary Prevention of Cardiovascular Disease with a Mediterranean Diet Supplemented with Extra-Virgin Olive Oil or Nuts. N. Engl. J. Med. 2018, 378, e34. [Google Scholar] [CrossRef] [PubMed]
  58. Guasch-Ferré, M.; Babio, N.; Martínez-González, M.A.; Corella, D.; Ros, E.; Martín-Peláez, S.; Estruch, R.; Arós, F.; Gómez-Gracia, E.; Fiol, M.; et al. Dietary fat intake and risk of cardiovascular disease and all-cause mortality in a population at high risk of cardiovascular disease. Am. J. Clin. Nutr. 2015, 102, 1563–1573. [Google Scholar] [CrossRef] [Green Version]
  59. Li, Y.; Hruby, A.; Bernstein, A.M.; Ley, S.H.; Wang, D.D.; Chiuve, S.E.; Sampson, L.; Rexrode, K.M.; Rimm, E.B.; Willett, W.C.; et al. Saturated Fats Compared with Unsaturated Fats and Sources of Carbohydrates in Relation to Risk of Coronary Heart Disease A Prospective Cohort Study. J. Am. Coll. Cardiol. 2015, 66, 1538–1548. [Google Scholar] [CrossRef] [Green Version]
  60. Hooper, L.; Martin, N.; Abdelhamid, A.; Davey Smith, G. Reduction in saturated fat intake for cardiovascular disease. Cochrane Database Syst. Rev. 2015, 6, CD011737. [Google Scholar] [CrossRef] [PubMed]
  61. Wang, X.; Lin, X.; Ouyang, Y.Y.; Liu, J.; Zhao, G.; Pan, A.; Hu, F.B. Red and processed meat consumption and mortality: Dose-response meta-analysis of prospective cohort studies. Public Health Nutr. 2016, 19, 893–905. [Google Scholar] [CrossRef] [Green Version]
  62. Burd, N.A.; McKenna, C.F.; Salvador, A.F.; Paulussen, K.J.M.; Moore, D.R. Dietary protein quantity, quality, and exercise are key to healthy living: A muscle-centric perspective across the lifespan. Front. Nutr. 2019, 6, 83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Katz, D.L.; Doughty, K.N.; Geagan, K.; Jenkins, D.A.; Gardner, C.D. Perspective: The Public Health Case for Modernizing the Definition of Protein Quality. Adv. Nutr. 2019, 10, 755–764. [Google Scholar] [CrossRef]
  64. Martínez-González, M.A.; Sánchez-Tainta, A.; Corella, D.; Salas-Salvadó, J.; Ros, E.; Arós, F.; Gómez-Gracia, E.; Fiol, M.; Lamuela-Raventós, R.M.; Schröder, H.; et al. A provegetarian food pattern and reduction in total mortality in the Prevención con Dieta Mediterránea (PREDIMED) study. Am. J. Clin. Nutr. 2014, 100 (Suppl. S1), 320S–328S. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Fung, T.T.; Chiuve, S.E.; McCullough, M.L.; Rexrode, K.M.; Logroscino, G.; Hu, F.B. Adherence to a DASH-style diet and risk of coronary heart disease and stroke in women. Arch. Intern. Med. 2008, 168, 713–720. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart of participants included in the study.
Figure 1. Flowchart of participants included in the study.
Nutrients 13 00972 g001
Figure 2. Hazard ratios (HR) and 95% confidence intervals (CI) between extreme quartiles for the association between Macronutrient Quality Index (MQI) and mortality in the SUN cohort. Abbreviations: FFQ, food frequency questionnaire, MQI, Macronutrient Quality Index; HPPQI, Healthy Plate Protein Source Quality Index. PQIa was defined as the following ratio: (fish + poultry and lean meats + eggs + pulses + nuts + reduced-fat dairy + whole grains)/(red and processed meats + full-fat dairy + refined cereals); PQIb was defined as the following ratio: (seafood + poultry + pulses + nuts + dairy products)/(red and processed meats).
Figure 2. Hazard ratios (HR) and 95% confidence intervals (CI) between extreme quartiles for the association between Macronutrient Quality Index (MQI) and mortality in the SUN cohort. Abbreviations: FFQ, food frequency questionnaire, MQI, Macronutrient Quality Index; HPPQI, Healthy Plate Protein Source Quality Index. PQIa was defined as the following ratio: (fish + poultry and lean meats + eggs + pulses + nuts + reduced-fat dairy + whole grains)/(red and processed meats + full-fat dairy + refined cereals); PQIb was defined as the following ratio: (seafood + poultry + pulses + nuts + dairy products)/(red and processed meats).
Nutrients 13 00972 g002
Table 1. Components of the Macronutrient Quality Index.
Table 1. Components of the Macronutrient Quality Index.
Components of Macronutrient Quality IndexIndex Range (Points)Criteria for Minimum IndexCriteria for Maximum Index
Carbohydrate quality index (CQI)1–5Minimum CQI (first quintile) Maximum CQI (fifth quintile)
Fat quality index (FQI)1–5Minimum FQI (first quintile)Maximum FQI (fifth quintile)
Healthy Plate Protein quality index (HPPQI)1–5Minimum HPPQI (first quintile)Maximum HPPQI (fifth quintile)
Total index (range)3–15
Table 2. Age and sex-adjusted baseline characteristics of participants according to quartiles of the Macronutrient Quality Index (MQI): the SUN (Seguimiento Universidad de Navarra) cohort 1,2.
Table 2. Age and sex-adjusted baseline characteristics of participants according to quartiles of the Macronutrient Quality Index (MQI): the SUN (Seguimiento Universidad de Navarra) cohort 1,2.
CharacteristicsQ1Q2Q3Q4
n (frequency)6646450939064022
MQI range<88–910–11>11
Marital status
Single42.843.443.445.9
Married52.150.850.947.5
Widowed1.01.21.00.8
Others4.24.74.75.7
Years at university 5.1 (1.5)5.1 (1.6)5.0 (1.5)5.0 (1.5)
Smoking (%)
Never smoker47.046.049.649.4
Current smoker 28.628.729.732.0
Former smoker24.425.320.718.6
Cumulative smoking habit (packs-years)6.1 (10.1)6.1 (10.1)5.7 (9.6)5.5 (9.1)
Alcohol intake (g/d)
Never18.517.317.318.0
<5 women/<10 men48.448.448.349.4
5–25 women/10–50 men31.332.232.331.2
>25 women/>50 men1.92.12.11.4
Physical activity (METs-h/week)19.1 (20.7)20.7 (21.8)23.0 (24.5)26.3 (25.9)
Body mass index (kg/m2)23.5 (3.6)23.6 (3.5)23.7 (3.6)23.5 (3.5)
Time spent sitting (h/d)5.3 (2.1)5.3 (2.0)5.2 (2.0)5.2 (2.0)
Snacking between meals (% Yes)35.434.331.529.5
Special diet (% Yes)5.27.19.314.6
Medically-diagnosed conditions at baseline
Prevalent diabetes (%)1.62.02.52.2
Prevalent hypertension (%)10.310.811.512.7
Prevalent dyslipemia (%)6.16.97.77.8
Prevalent cardiovascular disease (%) 34.54.24.95.0
Prevalent cancer (%)2.62.42.42.8
Dietary variables
Trichopoulou’s MedDiet score2.6 (1.3)3.6 (1.3)4.4 (1.3)5.3 (1.3)
Total energy intake (kcal/d)2304 (613)2355 (618)2369 (633)2364 (620)
Carbohydrate intake, % E41.7 (7.3)43.0 (7.1)43.9 (7.2)45.6 (7.8)
Fat intake, % E37.8 (6.1)36.7 (6.3)35.9 (6.6)34.3 (7.0)
Protein intake, % E18.4 (3.3)18.2 (3.1)18.1 (3.5)18.2 (3.5)
Abbreviations: MQI, Macronutrient Quality Index; Q, quartile; %E, percentage of energy intake. 1 Adjusted through inverse probability weighting. 2 Values are means ± SDs or numbers of participants (percentages) unless otherwise indicated. 3 Prevalent CVD was considered as any of the following medical diagnosis occurring prior to entering in the study: myocardial infarction, stroke, atrial fibrillation, coronary artery bypass grafting or other re-vascularization procedures, paroxysmal tachycardia, peripheral venous thrombosis, aortic aneurism, pulmonary embolism or heart failure.
Table 3. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Macronutrient Quality Index (MQI) and all-cause mortality in the SUN cohort (n = 19,083).
Table 3. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Macronutrient Quality Index (MQI) and all-cause mortality in the SUN cohort (n = 19,083).
Q1Q2Q3Q4p for Trend
n (frequency)6646450939064022
MQI range<88–910–11>11
Deaths13810211387
Person-years79,52652,39644,30443,661
Mortality rate/1000 person years1.741.952.551.99
Crude model1.00 (Ref.)1.10 (0.85–1.42)1.43 (1.12–1.84)1.10 (0.84–1.44)0.214
Model 11.00 (Ref.)0.98 (0.75–1.27)1.10 (0.85–1.42)0.71 (0.54–0.94)0.035
Model 21.00 (Ref.)1.07 (0.83–1.39)1.20 (0.93–1.56)0.84 (0.63–1.12)0.408
Model 31.00 (Ref.)1.03 (0.79–1.34)1.15 (0.89–1.50)0.79 (0.59–1.06)0.199
Abbreviations: MQI, Macronutrient Quality Index; Q, quartile; Ref., referent value. Model 1: adjusted for sex (dichotomous), age (underlying variable), total energy intake (continuous), and stratified by recruitment period (5 categories), and deciles of age. Model 2: additionally adjusted for marital status (four categories), educational level (years of higher education, continuous), smoking (three categories), cumulative smoking habit (packs/year, continuous), alcohol intake (four categories), leisure-time physical activity (metabolic equivalent-h/week, continuous), BMI (kg/m2, linear, continuous), time spent sitting (hours/week, continuous), weight gain in the previous 5 years before entering the cohort (dichotomous) and following a special diet at baseline (dichotomous). Model 3: additionally adjusted for family history of CVD (dichotomous), and prevalent diagnosis of diabetes (dichotomous), hypertension (dichotomous), hypercholesterolemia (dichotomous), hypertriglyceridemia (dichotomous), CVD (dichotomous), depression (dichotomous), and cancer (dichotomous).
Table 4. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Carbohydrate Quality Index (CQI), Fat Quality Index (FQI), Protein Quality Index (PQI), and all-cause mortality in the SUN cohort (n = 19,083).
Table 4. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Carbohydrate Quality Index (CQI), Fat Quality Index (FQI), Protein Quality Index (PQI), and all-cause mortality in the SUN cohort (n = 19,083).
CQIQ1Q2Q3Q4p for Trend
n (frequency)6038428153663398
CQI range<1010–1112–14>14
Deaths13910912963
Person-years71,60750,01460,96237,303
Mortality rate/1000 person years1.942.182.121.69
Crude model1.00 (Ref.)1.12 (0.87–1.44)1.06 (0.83–1.35)0.84 (0.62–1.13)0.337
Model 11.00 (Ref.)0.87 (0.67–1.12)0.86 (0.67–1.10)0.57 (0.42–0.78)0.001
Model 21.00 (Ref.)0.98 (0.75–1.27)0.97 (0.75–1.25)0.66 (0.48–0.91)0.024
Model 31.00 (Ref.)0.99 (0.76–1.28)0.99 (0.77–1.28)0.65 (0.47–0.90)0.023
Model 41.00 (Ref.)0.97 (0.74–1.27)0.96 (0.73–1.27)0.64 (0.45–0.90)0.021
FQIQ1Q2Q3Q4p for Trend
n (frequency)4771477147714770
FQI range<1.401.40–1.611.62–1.90>1.90
Deaths11693108123
Person-years55,42456,03954,35554,069
Mortality rate/1000 person years2.101.661.992.27
Crude model1.00 (Ref.)0.78 (0.59–1.03)0.94 (0.72–1.22)1.07 (0.83–1.37)0.276
Model 11.00 (Ref.)0.90 (0.69–1.19)0.90 (0.69–1.18)0.91 (0.70–1.18)0.547
Model 21.00 (Ref.)0.93 (0.70–1.24)0.99 (0.75–1.30)0.94 (0.72–1.22)0.723
Model 31.00 (Ref.)0.89 (0.67–1.18)0.92 (0.70–1.21)0.86 (0.66–1.12)0.337
Model 51.00 (Ref.)0.85 (0.64–1.13)0.86 (0.65–1.14)0.75 (0.56–1.00)0.070
HPPQIQ1Q2Q3Q4p for Trend
n (frequency)4771477147714770
PQI range<0.650.65–0.930.93–1.38>1.38
Deaths11587103135
Person-years58,04755,98853,81952,032
Mortality rate/1000 person years1.981.551.912.59
Crude model1.00 (Ref.)0.77 (0.59–1.02)0.93 (0.71–1.21)1.26 (0.98–1.61)0.007
Model 11.00 (Ref.)0.71 (0.54–0.94)0.72 (0.55–0.94)0.82 (0.63–1.05)0.470
Model 21.00 (Ref.)0.81 (0.61–1.08)0.79 (0.60–1.04)0.97 (0.74–1.25)0.781
Model 31.00 (Ref.)0.79 (0.60–1.05)0.73 (0.55–0.97)0.92 (0.70–1.19)0.945
Model 61.00 (Ref.)0.80 (0.60–1.07)0.74 (0.55–0.99)0.93 (0.69–1.25)0.844
Abbreviations: Ref., referent value; Q, quartile. Model 1: adjusted for sex (dichotomous), age (underlying variable), total energy intake (continuous) and stratified by recruitment period (5 categories) and deciles of age. Model 2: additionally adjusted for marital status (four categories), educational level (years of higher education, continuous), smoking (three categories), cumulative smoking habit (package/year, continuous), alcohol intake (four categories), leisure-time physical activity (metabolic equivalent-h/week, continuous), BMI (kg/m2, linear, continuous), time spent sitting (hours/week, continuous), weight gain in the previous 5 years before entering the cohort (dichotomous) and following a special diet at baseline (dichotomous). Model 3: additionally adjusted for family history of CVD (dichotomous), and prevalent diagnosis of diabetes (dichotomous), hypertension (dichotomous), hypercholesterolemia (dichotomous), hypertriglyceridemia (dichotomous), CVD (dichotomous), depression (dichotomous), and cancer (dichotomous). Model 4: additionally adjusted for snacking (dichotomous), animal protein (g/d, continuous), plant protein (g/d, continuous), MUFA (%E, continuous), PUFA (%E, continuous), SFA (%E, continuous), and TFA (%E, continuous). Model 5: additionally adjusted for snacking (dichotomous), animal protein (g/d, continuous), plant protein (g/d, continuous), and carbohydrate intake g/d, continuous). Model 6: additionally adjusted for snacking (dichotomous), MUFA (%E, continuous), PUFA (%E, continuous), SFA (%E, continuous), TFA (%E, continuous), and carbohydrate intake (g/d, continuous).
Table 5. Updated 1 and cumulative average 2 dietary measurements after 10 years of follow-up. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Macronutrient Quality Index (MQI) and mortality in the SUN cohort.
Table 5. Updated 1 and cumulative average 2 dietary measurements after 10 years of follow-up. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Macronutrient Quality Index (MQI) and mortality in the SUN cohort.
Updated Diet 1Q1Q2Q3Q4p for Trend
MQI range<77–910–11>11
Deaths9314910791
Person-years56,29777,42842,97243,190
Mortality rate/1000 person years1.651.922.492.11
Crude model1.00 (Ref.)1.18 (0.91–1.53)1.53 (1.16–2.02)1.25 (0.94–1.67)0.044
Model 11.00 (Ref.)1.04 (0.80–1.36)1.14 (0.86–1.52)0.81 (0.60–1.10)0.211
Model 21.00 (Ref.)1.11 (0.85–1.44)1.25 (0.93–1.67)0.97 (0.71–1.32)0.953
Model 31.00 (Ref.)1.06 (0.81–1.38)1.19 (0.89–1.60)0.90 (0.66–1.23)0.610
Cumulative Diet Average 2Q1Q2Q3Q4p for Trend
MQI range>7.57.5–99.5–11>11
Deaths13710411089
Person-years79,57853,32644,22842,754
Mortality rate/1000 person years1.721.952.492.08
Crude model1.00 (Ref.)1.10 (0.85–1.42)1.41 (1.10–1.81)1.17 (0.89–1.52)0.095
Model 11.00 (Ref.)0.97 (0.75–1.25)1.09 (0.84–1.41)0.75 (0.57–0.99)0.079
Model 21.00 (Ref.)1.06 (0.82–1.38)1.19 (0.92–1.55)0.90 (0.67–1.20)0.630
Model 31.00 (Ref.)1.02 (0.78–1.33)1.16 (0.89–1.51)0.84 (0.63–1.12)0.367
Abbreviations: Ref., referent value; Q, quartile. 1 Repeated measures: update information of MQI after 10 years of follow-up. 2 Repeated measures: cumulative average information of MQI (at baseline and after 10 years of follow-up). Model 1: adjusted for sex (dichotomous), age (underlying variable), total energy intake (continuous) and sex (dichotomous) stratified by recruitment period (5 categories) and deciles of age and recruitment period (5 categories). Model 2: additionally adjusted for marital status (four categories), educational level (years of higher education, continuous), smoking (three categories), cumulative smoking habit (package/year, continuous), alcohol intake (four categories), smoking (three categories), smoking habit (package/year, continuous), educational level (years of higher education, continuous), leisure-time physical activity (metabolic equivalent-h/week, continuous), BMI (kg/m2, linear, continuous), marital status (four categories), time spent sitting (hours/week, continuous), weight gain in the previous 5 years before entering the cohort (dichotomous) and following a special diet at baseline (dichotomous). Model 3: additionally adjusted for family history of CVD (dichotomous), and prevalent diagnosis of diabetes (dichotomous), hypertension (dichotomous), hypercholesterolemia (dichotomous), hypertriglyceridemia (dichotomous), CVD (dichotomous), depression (dichotomous), and cancer (dichotomous).
Table 6. Updated 1 and cumulative average 2 dietary measurements after 10 years of follow-up. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Carbohydrate Quality Index (CQI) and mortality in the SUN cohort.
Table 6. Updated 1 and cumulative average 2 dietary measurements after 10 years of follow-up. Hazard ratios (HR) and 95% confidence intervals (CI) for the association between quartiles of the Carbohydrate Quality Index (CQI) and mortality in the SUN cohort.
Updated Diet 1Q1Q2Q3Q4p for Trend
CQI range<1010–1112–14>14
Deaths13610812967
Person-years69,85448,97162,04339,018
Mortality rate/1000 person years1.952.212.081.71
Crude model1.00 (Ref.)1.13 (0.87–1.45)1.02 (0.80–1.29)0.82 (0.61–1.10)0.194
Model 11.00 (Ref.)0.92 (0.69–1.16)0.87 (0.68–1.12)0.60 (0.44–0.81)0.002
Model 21.00 (Ref.)1.00 (0.77–1.30)0.97 (0.75–1.25)0.69 (0.51–0.94)0.031
Model 31.00 (Ref.)1.01 (0.78–1.31)1.00 (0.77–1.29)0.68 (0.50–0.93)0.032
Model 41.00 (Ref.)1.00 (0.76–1.31)0.97 (0.74–1.28)0.68 (0.49–0.94)0.031
Cumulative Diet Average 2Q1Q2Q3Q4p for Trend
CQI range<9.59.5–1111.5–13.5>13.5
Deaths13311110393
Person-years68,63551,15246,44343,657
Mortality rate/1000 person years1.942.172.211.73
Crude model1.00 (Ref.)0.88 (0.68–1.14)0.93 (0.72–1.21)0.59 (0.46–0.78)<0.001
Model 11.00 (Ref.)0.89 (0.69–1.16)0.94 (0.72–1.23)0.61 (0.46–0.80)0.001
Model 21.00 (Ref.)1.01 (0.78–1.32)1.04 (0.79–1.37)0.71 (0.53–0.94)0.017
Model 31.00 (Ref.)1.02 (0.79–1.33)1.07 (0.82–1.41)0.70 (0.53–0.94)0.017
Model 41.00 (Ref.)1.01 (0.77–1.32)1.05 (0.79–1.39)0.69 (0.50–0.93)0.016
Abbreviations: Ref., referent value; Q, quartile. 1 Repeated measures: update information of MQI after 10 years of follow-up. 2 Repeated measures: cumulative average information of CQI (at baseline and after 10 years of follow-up). Model 1: adjusted for sex (dichotomous), age (underlying variable), total energy intake (continuous) and stratified by recruitment period (5 categories) and deciles of age. Model 2: additionally adjusted for marital status (four categories), educational level (years of higher education, continuous), smoking (three categories), cumulative smoking habit (package/year, continuous), alcohol intake (four categories), leisure-time physical activity (metabolic equivalent-h/week, continuous), BMI (kg/m2, linear, continuous), time spent sitting (hours/week, continuous), weight gain in the previous 5 years before entering the cohort (dichotomous) and following a special diet at baseline (dichotomous). Model 3: additionally adjusted for family history of CVD (dichotomous), and prevalent diagnosis of diabetes (dichotomous), hypertension (dichotomous), hypercholesterolemia (dichotomous), hypertriglyceridemia (dichotomous), CVD (dichotomous), depression (dichotomous), and cancer (dichotomous). Model 4: additionally adjusted for snacking (dichotomous), animal protein (g/d, continuous), plant protein (g/d, continuous), MUFA (%E, continuous), PUFA (%E, continuous), SFA (%E, continuous), and TFA (%E, continuous).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Santiago, S.; Zazpe, I.; Fernandez-Lazaro, C.I.; de la O, V.; Bes-Rastrollo, M.; Martínez-González, M.Á. Macronutrient Quality and All-Cause Mortality in the SUN Cohort. Nutrients 2021, 13, 972. https://0-doi-org.brum.beds.ac.uk/10.3390/nu13030972

AMA Style

Santiago S, Zazpe I, Fernandez-Lazaro CI, de la O V, Bes-Rastrollo M, Martínez-González MÁ. Macronutrient Quality and All-Cause Mortality in the SUN Cohort. Nutrients. 2021; 13(3):972. https://0-doi-org.brum.beds.ac.uk/10.3390/nu13030972

Chicago/Turabian Style

Santiago, Susana, Itziar Zazpe, Cesar I. Fernandez-Lazaro, Víctor de la O, Maira Bes-Rastrollo, and Miguel Ángel Martínez-González. 2021. "Macronutrient Quality and All-Cause Mortality in the SUN Cohort" Nutrients 13, no. 3: 972. https://0-doi-org.brum.beds.ac.uk/10.3390/nu13030972

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop