Next Article in Journal
The Impact of Pectin Supplementation on Intestinal Barrier Function in Healthy Young Adults and Healthy Elderly
Next Article in Special Issue
The Willingness to Modify Portion Sizes or Eat New Protein Foods Largely Depends on the Dietary Pattern of Protein Intake
Previous Article in Journal
Plasma Vitamin C Levels: Risk Factors for Deficiency and Association with Self-Reported Functional Health in the European Prospective Investigation into Cancer-Norfolk
Previous Article in Special Issue
Brazilian Food Truck Consumers’ Profile, Choices, Preferences, and Food Safety Importance Perception
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Provegetarian Food Pattern Emphasizing Preference for Healthy Plant-Derived Foods Reduces the Risk of Overweight/Obesity in the SUN Cohort

1
Department of Preventive Medicine and Public Health, School of Medicine, University of Navarra, 31008 Pamplona, Spain
2
Biomedical Research Centre Network of Physiopathology of Obesity and Nutrition (CIBERobn), Institute of Health Carlos III, 28029 Madrid, Spain
3
Navarra Institute for Health Research (IdiSNA), 31008 Pamplona, Spain
4
Department of Nutrition, Harvard TH Chan School of Public Health, Boston, MA 02115, USA
5
Department of Nutrition, Food Sciences and Physiology, School of Pharmacy and Nutrition, University of Navarra, 31008 Pamplona, Spain
6
Department of Pharmacology and Toxicology, School of Pharmacy and Nutrition, University of Navarra, 31008 Pamplona, Spain
7
Food Design and Consumer Behavior Section, Department of Food Science, University of Copenhagen, Rolighedsvej 26, 1958 Frederiksberg C, Denmark
*
Author to whom correspondence should be addressed.
Nutrients 2019, 11(7), 1553; https://doi.org/10.3390/nu11071553
Submission received: 21 May 2019 / Revised: 24 June 2019 / Accepted: 5 July 2019 / Published: 9 July 2019
(This article belongs to the Special Issue Food Choice and Nutrition)

Abstract

:
Provegetarian diets (i.e., preference for plant-derived foods but not exclusion of animal foods) have been associated with a reduced risk of long-term weight gain and could be more easily embraced than strict vegetarian diets. However, not all plant-derived foods are equally healthy. In the “Seguimiento Universidad de Navarra” (SUN) cohort, we prospectively evaluated the association between different provegetarian food patterns and the incidence of overweight/obesity in 11,554 participants with initial body mass index <25 kg/m2. A provegetarian food pattern (FP) was built by assigning positive scores to plant foods and reverse scores to animal foods. A healthful and an unhealthful provegetarian FP, which distinguished between healthy (fruits/vegetables/whole grains/nuts/legumes/olive oil/coffee) and less-healthy plant foods (fruit juices/potatoes/refined grains/pastries/sugary beverages), were also built. A total of 2320 new cases of overweight or obesity were identified after a median follow-up of 10.3 years. Higher baseline conformity with the overall provegetarian FP was inversely associated with overweight/obesity (HR comparing extreme quintiles: 0.85; 95% CI: 0.75 to 0.96; p-trend: 0.014). This association was stronger for the healthful FP (HR: 0.78; 95% CI: 0.67 to 0.90; p-trend: <0.001) and was not apparent for the unhealthful FP (HR: 1.07; 95% CI: 0.92 to 1.23; p-trend: 0.551). In a large prospective cohort of relatively young adults, better conformity with a healthy provegetarian diet was associated with a reduced long-term risk of overweight/obesity, whereas no consistent trend was found for a FP that emphasized less-healthy plant foods.

1. Introduction

It is well-established that diets higher in plant-based foods and with fewer animal source foods confer both health and environmental benefits [1,2,3]. Plant-based diets have been associated with a lower risk of various diseases [4,5,6,7,8,9,10,11], including a lower risk of obesity [12,13,14,15,16,17,18,19], which is among the greatest contributors to premature mortality [20,21,22]. Moreover, whereas strict vegan or vegetarian diets might not be easily embraced within the general population, provegetarian diets (i.e., preferentially but not exclusively consuming plant-derived foods) based on gradual dietary changes may be easier to adopt. Although this type of plant-based diet is also commonly known as flexitarian, we used the provegetarian operational score because it has been defined in a quantitative and easily reproducible way [10]. However, with the exception of several recent studies [23,24,25,26,27], prior studies of plant-based diets do not take into account the fact that not all plant-derived foods are equally healthy.
Since the potential mechanisms responsible for the beneficial effects of plant-based diets on body weight [28,29,30,31,32,33,34,35,36,37] are mainly supported by dietary constituents like fiber and other biologically active micronutrients which are characteristic of healthy plant foods (i.e., fruits and vegetables), it is warranted to evaluate whether the healthfulness of plant-derived foods influences the risk of long-term weight gain. Additionally, it is important to provide evidence towards the advantages and disadvantages of strategies that consumers may put in practice when they adopt plant-based diets following the societal trends, and to provide insights for tomorrow’s food innovation. Therefore, we aimed to prospectively examine the association between different versions of a provegetarian food pattern and the incidence of overweigh/obesity in the SUN (Seguimiento Universidad de Navarra/University of Navarra Follow-up) cohort.

2. Materials and Methods

2.1. Study Population

The “Seguimiento Universidad de Navarra” (SUN) prospective cohort study started on December 1999 following the models of the Nurses’ Health Study and the Health Professionals Follow-Up Study. It is a multipurpose and dynamic cohort (i.e., recruitment is permanently open) composed of young middle-aged university graduates from different Spanish universities, including the University of Navarra. Baseline assessment and follow-up information is gathered through postal or web-based questionnaires completed by participants every 2 years. Self-administered questionnaires include information on sociodemographic, lifestyle and dietary variables as well as the use of medication and the prevalence or incidence of diseases during follow-up. For the main outcomes, medical records are requested to participants or their next of kin in order to confirm the diagnoses by a panel of medical doctors of the study, blinded to the dietary habits of participants and other exposures of interest. For participants lost to follow-up, the National Death Index is periodically checked to determine their vital status and cause of death. The overall retention in the cohort is higher than 90%. All procedures were approved by the Institutional Review Board of the University of Navarra (2017.151) according to Declaration of Helsinki guidelines. Voluntary completion of the baseline self-administrated questionnaire was considered to imply informed consent, as participants received detailed information about the study. Further details about the methodology and characteristics of participants have been described in more detail elsewhere [38].
Up to July 2018, there were 22,791 participants who had completed the baseline questionnaire. In order to ensure that all participants had completed at least one follow-up questionnaire, only those who were recruited before October 2015 were considered eligible. Out of the 22,467 eligible participants, we excluded 6655 with prevalent overweight or obesity (body mass index (BMI) ≥ 25 kg/m2) at baseline, 1628 who reported implausible values for total energy intake according to predefined limits (>4000 kcal/d in men and >3500 kcal/d in women or <800 kcal/d in men and <500 kcal/d in women), 95 women who were pregnant at baseline, 714 with prevalent chronic disease (diabetes, cardiovascular disease and/or cancer), and 393 with a weight change > 10 kg in the 5 years before entering the study to reduce the potential source of confounding by other causes of weight changes. Among the remaining 12,982 participants, 1160 were lost to follow-up leaving a total of 11,822 (retention = 91%). Additionally, we excluded 268 participants who had missing values in the variables of interest regarding the outcome. After exclusions, a sample of 11,554 participants was included in this study (Figure 1).

2.2. Dietary Assessment

A 136-item semi-quantitative food frequency questionnaire (FFQ) repeatedly validated in Spanish participants [39,40] was completed at baseline. Frequencies of consumption were measured in 9 categories (ranging from never/almost never to >6 servings/day) for each food item, in order to assess typical food intake over the previous year. Daily intake (grams/day) was calculated by multiplying the specified portion size for each food item by the frequency of consumption. Nutrient intake was updated by a trained dietician using the latest available information on Spanish food composition tables [41,42]. The scoring criteria for the overall, healthful and unhealthful provegetarian food patterns (FPs) are showed in Table 1, as described in previous studies [10,23]. Information on the food items from the FFQ included in each food group is shown in Table S1. Differences in food consumption between the lowest and highest quintiles of the provegetarian FPs are graphically represented in Figures S1–S3, showing differences (in percentage) between medians of extreme quintiles (quintile 1 and 5) and the median of the whole sample.
Adherence to an overall provegetarian FP was calculated according to the score proposed by Martínez-González et al. [10], which quantifies the habit of preferentially consuming plant-derived foods instead of animal-derived foods without the need to follow a strict vegetarian diet. This score quantifies the energy-adjusted consumption (g/d) of seven plant food groups (fruits, vegetables, potatoes, nuts, legumes, cereal grains and olive oil) and five animal food groups (dairy, eggs, meat, fish and seafood and animal fat). The residual method was used separately for men and women to adjust for total energy intake and the energy-adjusted estimates (residuals) were ranked according to their sex-specific quintiles. The quintile values for animal food groups were reversed (a value of 5 was assigned for the first quintile, 4 for the second quintile, and successively until the value of 1 was assigned to the fifth quintile). To obtain the overall provegetarian FP score, quintile values of plant foods and reverse quintile values of animal foods were summed. Thus, the final scores assigned to participants could range from 12 (lowest adherence) to 60 (highest adherence).
Likewise, adherence to a healthful and unhealthful provegetarian FP was calculated according to the scoring criteria suggested by Satija et al. [23]. Healthy (fruits, vegetables, whole grains, nuts, legumes, olive oil, coffee) and less-healthy plant foods (fruit juices, potatoes, refined grains, pastries, sugary beverages) were distinguished. As shown in Table 1, to build the healthful provegetarian FP, positive scores were assigned to healthy plant foods and reverse scores to less-healthy plant foods as well as to animal foods. For the unhealthful provegetarian FP, positive scores were assigned to less-healthy plant foods and reverse scores to healthy plant foods and animal foods. Quintiles and reverse quintiles were summed to obtain scores of the healthful and unhealthful versions of a provegetarian FP. Participants were categorized into energy-adjusted quintiles of the healthful and unhealthful provegetarian FPs. Thus, their final scores could range from 18 (lowest adherence) to 90 (highest adherence).
Adherence to the Mediterranean diet was assessed using the a priori 9-item Mediterranean diet score (MDS) proposed by Trichopoulou et al. [43].

2.3. Other Covariates Assessment

The baseline questionnaire gathered information about sociodemographic and anthropometric characteristics, health-related habits, and clinical data including personal and family history of disease. Physical activity was assessed using a previously validated questionnaire, which included 17 different activities during leisure time [44]. Metabolic equivalents (METs) were calculated to yield MET-hours per week scores for each participant.
Baseline dietary intake data derived from food composition tables for Spain [41,42] were used to identify the fat quality index (FQI) and the carbohydrate quality index (CQI). The FQI was calculated using the ratio (monounsaturated + polyunsaturated fatty acids)/(saturated + trans fatty acids) as a continuous variable, ranging from 1.01 to 4.93. The CQI was determined by summing up quintiles of the following four criteria: Dietary fiber intake (g/day), glycemic index, the ratio whole grains/total grains and the ratio solid carbohydrates/solid + liquid carbohydrates. All criteria had the same weighting, and the CQI ranged from 4 to 20. The FQI and CQI have been previously used in this cohort [45] and also in the Prevención con Dieta Mediterránea (PREDIMED) study [46].

2.4. Outcome Ascertainment

Weight of participants was self-reported at baseline and every 2 years in the follow-up questionnaires. BMI was calculated as the self-reported weight in kilograms divided by the square of height in meters. Weight was self-reported by participants at baseline and every 2 years of follow-up. BMI was calculated as weight in kilograms divided by the square of height in meters, which was ascertained at baseline. Anthropometric data (i.e., self-reported weight and BMI) have been previously validated in a subsample of this cohort finding strong correlation results [47]. The outcomes were: (1) Incidence of overweight/obesity (BMI ≥ 25 kg/m2); and (2) average change in body weight (grams) per year (g/y).

2.5. Statistical Analysis

Inverse probability weighting was used to determine the age- and sex-adjusted baseline characteristics of participants according to quintiles of adherence to the different provegetarian FPs in order to remove differences that were only explained by the different age and sex distribution across baseline quintiles of provegetarian FPs.
Cox proportional hazards analyses were conducted to assess the relationship between provegetarian FPs with varying degree of healthfulness and the incidence of overweight/obesity. All models included age as underlying time variable and were stratified by age groups and year of recruitment. After crude analyses, a model adjusted for sex and age was fitted. Multivariable models were additionally adjusted for known risk factors of weight gain and potential confounders by using baseline values of the following covariates: Baseline BMI (kg/m2, continuous), physical activity (METs-h/week, quartiles), hours of TV watching (quartiles), smoking status (never, current, or former), marital status (single, married, other), years of university education (continuous), total energy intake (kcal/day, continuous), snacking between meals (yes, no), following a special diet at baseline (yes, no), parental family history of obesity (yes, no) and hours of siesta sleep (0, >0 to ≤0.5, >0.5). Robust standard errors were used. In order to reduce the effects of potential measurement errors due to diet variation during follow-up, dietary data were updated after 10 years of follow-up if the participant was followed for longer than 10 years and had completed the 10-year follow-up FFQ. Tests of linear trends were conducted using the medians of each provegetarian score for each quintile and treating the variable as continuous in Cox models. Moreover, to minimize the potential distortions produced by different measurement units in each provegetarian score, z scores were used. Each z score was calculated for each participant as the score value of the participant minus the mean value for the overall sample divided by the standard deviation of the score in the overall sample.
Multiple linear regression models were used to assess the association between the provegetarian FPs and the average yearly weight change during follow-up, considering the first quintile of adherence as the reference group. We estimated β regression coefficients (and their 95% confidence interval (CI)) for the other four quintiles, which should be interpreted as the difference in average yearly weight change (g/y) for each of the upper four quintiles versus the lowest quintile. To take into account repeated measurements of diet, generalized estimating equations with an unstructured correlation matrix were used to assess the relationship between updated scores of the provegetarian FPs and average yearly weight change. Analyses were adjusted for the same confounding factors as the main Cox regression analysis.
Additional analyses were conducted after excluding participants with a baseline BMI ≥ 30 kg/m2 instead of BMI ≥ 25 kg/m2 (as in the main analyses, which only included normal weight participants) in order to assess the incidence of obesity among overweight participants (i.e., change of BMI category from overweight to obese).
Sensitivity analyses were conducted to explore the robustness of our findings by rerunning the analyses under the following scenarios: Excluding participants with energy intake between the 5th and 95th percentiles; excluding participants with no answer in > 12 items in the 136-item baseline FFQ; including participants with weight change > 10 kg over the previous 5 years before entering the study; excluding participants with BMI > 24.5 kg/m2 at baseline (even though self-reported BMI has been previously validated finding good correlation results, some degree of misclassification might still exist); excluding early incident cases of overweight (participants who became overweight or obese only during the first 2 y of follow-up); additionally adjusting for weight gain ≥ 3 kg over the previous 5 years before entering the cohort; considering as the outcome only obesity (BMI ≥ 30 kg/m2) instead of overweight (BMI ≥ 25 kg/m2) and truncating follow-up at 10 years.
All analyses were performed using STATA/SE version 12.0 (StataCorp, College Station, TX, USA).

3. Results

A total of 11,554 participants (8419 women and 3135 men; 34.7 ± 10.8 years [mean age ± SD]) were followed for a median of 10.3 years. Age-and sex-adjusted baseline characteristics of participants according to quintiles of adherence to the different provegetarian FPs are presented in Tables S2–S4. Compared with participants in the first quintile, those in the fifth quintile with higher scores on the overall and healthful provegetarian FPs were more likely to live alone, follow a special diet, be non-smokers, be more physically active, have a family history of obesity, and have a higher total energy and dietary fiber intake. Conversely, those with high scores on the unhealthful provegetarian FP were less physically active, less likely to live alone, more likely to smoke, have snacks, have a higher total energy intake and lower total fiber intake than those with low scores.
Over 113,213 person-years of follow-up, a total of 2320 new cases of overweight/obesity (BMI ≥ 25 kg/m2) were identified. Hazard ratios (HR) and their 95% confidence intervals (CI) for the risk of overweight/obesity in the age-and sex-adjusted and multivariable-adjusted models are shown in Table 2. In the fully adjusted models, the overall provegetarian FP was inversely associated with overweight/obesity incidence (HR comparing extreme quintiles: 0.85; 95% CI: 0.75 to 0.96; HR for each additional 2 z-score units: 0.89; 95% CI: 0.82 to 0.97; p-trend = 0.014). When we distinguished among healthful and unhealthful provegetarian FPs, we found a stronger inverse association between the healthful provegetarian FP and overweight/obesity incidence (HR comparing extreme quintiles: 0.78; 95% CI: 0.67 to 0.90; HR for each additional 2 z-score units: 0.84; 95% CI: 0.77 to 0.92; p-trend < 0.001), whereas null findings were observed for the unhealthful provegetarian FP (HR comparing extreme quintiles: 1.07; 95% CI: 0.92 to 1.23; HR for each additional 2 z-score units: 1.00; 95% CI: 0.91 to 1.10; p-trend = 0.551). Results from the multivariable-adjusted models are also graphically represented in Figure 2, showing HRs and 95% CIs of overweight/obesity according to baseline quintiles of the 3 provegetarian FPs.
Consistently, absolute average yearly weight change (g/y) modestly decreased across quintiles of overall and healthful provegetarian FPs (i.e., higher adherence to both overall and healthful provegetarian FPs was associated with less weight gain over time) and remained fairly constant across quintiles of the unhealthful provegetarian FP (Table 3). The average difference in yearly weight gain between extreme quintiles of adherence to the overall provegetarian FP was −167 g (95% CI: −249 to −84) and for the healthful provegetarian FP it was −202 g (95% CI: −294 to −110).
In an additional analysis conducted to assess the effect of several provegetarian FPs on incidence of obesity among overweight participants, a total of 669 new cases of obesity (BMI ≥ 30 kg/m2) were identified. HRs and their 95% CIs for the risk of obesity in the age-and sex-adjusted and multivariable-adjusted models are shown in Table 4. Fully adjusted estimates suggested that the inverse association between the overall and healthful provegetarian food patterns and overweight/obesity among normal weight participants is not as evident among initially overweight subjects (HR comparing extreme quintiles of the overall provegetarian food pattern: 0.79; 95% CI: 0.61 to 1.02; p-trend = 0.056. HR comparing extreme quintiles of the healthful provegetarian food pattern: 0.79; 95% CI: 0.59 to 1.05; p-trend = 0.066).
In ancillary analyses (data not shown), we investigated the association between adherence to an unhealthful provegetarian FP exclusively composed of plant-based foods (less-healthy plant foods positively weighted and healthy plant foods negatively weighted) and risk of overweight/obesity. We also explored the association between an unhealthful provegetarian FP composed of less-healthy plant foods (positively weighted) and animal foods (negatively weighted) and risk of overweight/obesity. The purpose of running these ancillary analyses was to study the independent effects of consuming more less-healthy plant foods compared to healthier plant foods (HR comparing extreme quintiles: 1.21; 95% CI: 1.05 to 1.40; p-trend: (0.021) and compared to more animal foods (HR comparing extreme quintiles: 0.91; 95% CI: 0.80 to 1.05; p-trend: (0.130).
Sensitivity analyses were conducted to assess the robustness of our results under a variety of scenarios (Figure 3 and Table S5). The associations of overall and healthful provegetarian FPs with the risk of overweight/obesity were attenuated when participants with no answer in >12 out of the 136 items included in the baseline FFQ were excluded and when obesity (BMI ≥ 30 kg/m2) was the main outcome considered. Slightly stronger associations were found when participants with baseline BMI > 24.5 kg/m2 were excluded. The adverse association with the unhealthful provegetarian FP became stronger when the follow-up period was truncated at 10 years.

4. Discussion

In a Mediterranean dynamic prospective cohort study composed of university graduates with an initial BMI < 25 kg/m2 (normal weight), we found that a general preference for plant-derived foods (provegetarian FP) was modestly associated with lower risk of subsequently developing overweight/obesity during follow-up. This inverse association was considerably stronger when we specifically evaluated the preference for healthy plant-derived foods (healthful provegetarian FP), whereas no such association was found for the preference for less healthy plant-derived foods (unhealthful provegetarian FP). Our findings were robust in sensitivity analyses and remained consistent and significant when updating diet measurements with the use of repeated nutritional data obtained from the FFQ administered after 10 years of follow-up. Additional analyses including normal weight and overweight participants with an initial BMI < 30 kg/m2 revealed that we cannot be equally certain about the inverse association between these provegetarian food patterns and obesity among overweight subjects. These results are in line with recent findings from three ongoing prospective cohort studies in the United States, which were published while our paper was under review [25].
A-priori defined provegetarian indices were used to assess the effect of preference for plant-derived foods compared with preferential selection of foods from animal sources on the risk of new-onset overweight/obesity. This approach builds on previous research and includes well-defined food groups with varying nutritional content. Thus, results are in close agreement with previous findings, which could be expected, as most healthy plant foods positively weighed in the healthful provegetarian FP (vegetables, fruits, legumes, whole grains, nuts, olive oil, coffee) have been consistently associated with lower risk of weight gain. Likewise, less-healthy plant foods positively weighed in the unhealthful provegetarian FP (fruit juices, potatoes, refined grains, sugary beverages, pastries) have been consistently associated with higher risk of weight gain [48]. Noteworthy, most of the cereal grains consumed by our participants were refined, especially white bread, which accounted for 49% of the cereal group. In addition, according to the NOVA (a name, not an acronym) food classification system that categorizes foods based on the extent and purpose of food processing, plant-based foods such as sugary beverages and pastries can be considered ultra-processed foods. Their convenience (imperishable, ready-to-consume), hyper-palatability, branding and aggressive marketing give ultra-processed foods enormous market advantages over minimally processed healthy food, and their consumption has been systematically associated with the deterioration of overall nutritional quality of diets and with several chronic non-communicable diseases including obesity [49]. Therefore, from a public health perspective, focusing solely on plant or animal origin of foods may not be the most effective message in terms of obesity prevention.
The associations observed between a general preference for plant-derived foods and lower risk of overweight/obesity are consistent with previous research [12,13,14,15,16,17,18,19,20]. In particular, preference for less-healthy plant foods was associated with an increased risk of overweight/obesity compared with preferential selection of healthier plant foods, as could be expected. However, preference for less-healthy plant foods did not affect the risk of overweight/obesity compared with preference for foods from animal sources. These results suggest that plant-based foods should be prioritized over animal foods, and reinforce the importance of a preference for healthier plant foods in terms of weight gain prevention. Similarly, a recent study found that exclusion of less healthy plant-based foods from a plant-based diet index did not substantially change its beneficial effect, while exclusion of the healthy plant-based foods combined from the plant-based diet index moderately attenuated the inverse association. Therefore, it seems that healthy plant foods could apparently be driving the observed beneficial effect. The associations were also moderately attenuated by giving healthier animal foods positive scores [12]. Additionally, another study found that a greater consumption of sugar-sweetened beverages and ultra-processed food was associated with a higher risk of being overweight among vegetarians [50], which further supports our results. Hence, an additional implication of this study is to call for evidence regarding the effects on health and sustainability of innovatively processed plant-based products that are becoming increasingly common among consumers.
The potential mechanisms underlying our findings may be related to the preventive effect against overweight/obesity of a high intake of foods containing dietary fiber and bioactive compounds, which are naturally found in fruits and vegetables. Minimally processed plant-based foods are typically low in energy density and contain high concentrations of micronutrients and phytochemicals like flavonoids that may have a beneficial influence on energy homeostatic pathways [51]. Other components like fiber may also increase satiation and positively modulate gut microbiota composition [30]. In fact, we found evidence that fiber may mediate the beneficial effect of a provegetarian diet on overweight/obesity (i.e., the inverse association was attenuated after additionally adjusting for total fiber intake [g/day]). Another study in the SUN cohort found a significant inverse association between the CQI and the incidence of overweight/obesity, which highlights the importance of preferentially consuming high-quality carbohydrates in terms of obesity prevention [52]. Moreover, vegetable-and fruit-based dietary patterns have been suggested to modulate inflammation status [53], and a previous study in our cohort found that a higher proinflammatory diet was associated with a higher risk of overweight and obesity [54].
Beyond the health of individuals, and taking into consideration the health of our planet, it is also known that the adoption of provegetarian food patterns could significantly reduce greenhouse gas emissions [1] as well as other environmental impacts, and it is not as expensive as adopting other healthy and environmentally friendly dietary patterns such as the Mediterranean Diet (MedDiet) [55]. This focus may be a useful strategy to change dietary behaviors in the general public towards a healthier diet, as it is a way to connect with deeper needs and desires of people while improving their health as a side effect [56]. Although the MedDiet is also associated with a significantly lower risk of overweight/obesity in our cohort (HR [95% CI] for the highest category of adherence to the MDS [range 6–9] vs the lowest [range 0–3]: 0.78 [0.69 to 0.87]) and represents a food pattern rich in plant-derived foods, there are some important differences between the provegetarian food pattern and the MedDiet in relation to the scoring criteria for fish, potatoes and alcohol. In fact, the provegetarian FP and its healthful version were only moderately correlated (Pearson correlation coefficient: 0.54) with the traditional MedDiet (as defined by the MDS score proposed by Trichopoulou et al. [43]), reflecting that these are novel diet indices that capture unique aspects of a pro-vegetarian diet. Importantly, linking obesity with other major global challenges like climate change and health of the planet focuses attention on the scale and urgency of addressing these challenges and may increase demand for policy action by the public [57]. Our study reinforces current evidence and recommendations that suggest a shift towards diets rich in plant foods with lower intake of animal foods, and emphasizes the importance of prioritizing healthy plant foods. Finally, the present study also provides insights for the development of innovative foods of plant origin by food scientists and the food industry. The design of new foods should be performed within boundaries of key nutrients (e.g., added sugar and salt), but also boundaries in the number of ingredients and their level of processing.
The strengths of the present study include the prospective longitudinal design with a long follow-up period and a relatively large sample size; the restriction to participants with baseline BMI < 25 kg/m2, which minimized the possibility of reverse causation bias; and the use of repeated measurements of the diet, which reduced potential measurement errors and increase validity of well-validated self-reported comprehensive food-frequency methods such as the SUN cohort FFQ that has been repeatedly validated in Spanish participants [36,37].
The present study also has some limitations. Although we adjusted for a wide array of potential confounders, we cannot exclude the possibility of residual confounding given the observational nature of the study design and the fact that adherence to the overall and healthful provegetarian was higher among participants with healthier characteristics (non-smokers, higher levels of physical activity). However, the restriction to highly educated participants was applied to control for confounding by socioeconomic status, which is an approach considered in epidemiology as an excellent technique to reduce confounding [58]. This also implies that our sample from the SUN cohort is not representative of the general Spanish population. However, most cohorts are usually non-representative, and “it is not representativeness of the study subjects that enhances the generalization, it is knowledge of specific conditions and an understanding of mechanism that makes for a proper generalization” [59]. Therefore, the dose-response trend found in our results, as well as the consistency with previous studies and the existence of biological plausibility confer high potential of causality for the benefits of a plant-based dietary pattern, particularly a healthy provegetarian pattern, in the prevention of long-term weight gain. Another potential limitation is that weight and dietary information were self-reported so, although they were validated as previously mentioned, we cannot rule out some degree of nondifferential misclassification. Moreover, whole grain consumption has only recently started to be more common in Spain, so the only whole grain product included in the FFQ was brown bread. Therefore, we could not fully consider the effect of a higher consumption of other whole grains on the risk of incident overweight and obesity.
In conclusion, higher adherence to a provegetarian FP emphasizing preference for healthy plant-derived foods was associated with a lower risk of developing overweight and obesity in a cohort of Spanish university graduates with initial low body mass index. Given that excess weight is also a critical risk factor for developing other increasingly frequent non-communicable chronic diseases, a generalized adoption of healthful provegetarian diets is an important target that should be considered in future public health policies.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6643/11/7/1553/s1, Table S1. Food groups included in the provegetarian food patterns and food items from the SUN Project baseline FFQ, Table S2. Age-and sex-adjusted baseline characteristics (means [SDs] for continuous variables and percentages for dichotomous variables) by quintiles of the provegetarian food pattern in the SUN cohort, Table S3. Age-and sex-adjusted baseline characteristics (means [SDs] for continuous variables and percentages for dichotomous variables) by quintiles of the healthful provegetarian food pattern in the SUN cohort, Table S4. Age-and sex-adjusted baseline characteristics (means [SDs] for continuous variables and percentages for dichotomous variables) by quintiles of the unhealthful provegetarian food pattern in the SUN cohort, Table S5: Sensitivity analyses. Hazard Ratios (95% confidence intervals) of incident overweight/obesity for extreme quintiles (Q5 vs. Q1) of adherence to the different provegetarian food patterns. Figure S1. Percentage differences between extreme quintiles (quintile 1 and 5) and median scores (quintile 3) for the consumption of each food category of the provegetarian food pattern: Overall study population (The SUN Project, 1999–2015). Figure S2. Percentage differences between extreme quintiles (quintile 1 and 5) and median scores (quintile 3) for the consumption of each food category of the healthful provegetarian food pattern: Overall study population. Figure S3. Percentage differences between extreme quintiles (quintile 1 and 5) and median scores (quintile 3) for the consumption of each food category of the unhealthful provegetarian food pattern: overall study population.

Author Contributions

Conceptualization, M.Á.M.G., M.B.R. and F.J.A.P.C.; methodology, M.Á.M.G. and M.B.R.; validation, M.Á.M.G. and M.B.R.; formal analysis, C.G.D., A.G. and J.S.S.; investigation, M.Á.M.G. and M.B.R.; resources, M.Á.M.G. and M.B.R.; data curation, M. Á.M.G.; writing—original draft preparation, C.G.D.; writing—review and editing, all authors; supervision, M.B.R. and J.A.M.; project administration, M.Á.M.G. and M.B.R.; funding acquisition, M. Á.M.G. and M.B.R.

Funding

This research was funded by the Spanish Government-Instituto de Salud Carlos III, the 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, and G03/140), the Navarra Regional Government (45/2011, 122/2014), and the University of Navarra. C.G.D. was supported by a predoctoral contract for training in health research (PFIS) (FI18/00073) of the Instituto de Salud Carlos III.

Acknowledgments

The authors are indebted to the participants of the SUN cohort for their continued cooperation and participation. We are also grateful to the members of the Department of Nutrition of the Harvard School of Public Health (Willett W.C., Hu F.B., and Ascherio A.) who helped us to design the SUN study. We also thank the other members of the SUN Group.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Springmann, M.; Clark, M.; Mason-D’Croz, D.; Wiebe, K.; Bodirsky, B.; Lassaletta, L.; de Vries, W.; Vermeulen, S.; Herrero, M.; Carlson, K.; et al. Options for keeping the food system within environmental limits. Nature 2018, 562, 519–525. [Google Scholar] [CrossRef] [PubMed]
  2. Godfray, H.C.J.; Aveyard, P.; Garnett, T.; Hall, J.W.; Key, T.J.; Lorimer, J.; Pierrehumber, R.T.; Scarborough, P.; Springmann, M.; Jebb, S.A. Meat consumption, health, and the environment. Science 2018, 361, 5324. [Google Scholar] [CrossRef] [PubMed]
  3. Willett, W.; Rockström, J.; Loken, B.; Springmann, M.; Lang, T.; Vermeulen, S.; Garnett, T.; Tilman, D.; Declerck, F.; Wood, A.; et al. Food in the Anthropocene: The EAT–Lancet Commission on healthy diets from sustainable food systems. Lancet 2019, 393, 447–492. [Google Scholar] [CrossRef]
  4. Key, T.; Fraser, G.; Thorogood, M.; Appleby, P.; Beral, V.; Reeves, G.; Burr, M.; Chang-Claude, J.; Frentzel-Beyme, R.; Kuzma, J.; et al. Mortality in Vegetarians and Nonvegetarians: Detailed Findings from A Collaborative Analysis Of 5 Prospective Studies. Am. J. Clin. Nutr. 1999, 70, 516s–524s. [Google Scholar] [CrossRef] [PubMed]
  5. Fraser, G.E. Vegetarian diets: What do we know of their effects on common chronic diseases? Am. J. Clin. Nutr. 2009, 89, 1607S–1612S. [Google Scholar] [CrossRef] [PubMed]
  6. McEvoy, C.T.; Temple, N.; Woodside, J.V. Vegetarian diets, low-meat diets and health: A review. Public Health Nutr. 2012, 15, 2287–2294. [Google Scholar] [CrossRef] [PubMed]
  7. Huang, T.; Yang, B.; Zheng, J.; Li, G.; Wahlqvist, M.; Li, D. Cardiovascular disease mortality and cancer incidence in vegetarians: A meta-analysis and systematic review. Ann. Nutr. Metab. 2012, 60, 233–240. [Google Scholar] [CrossRef] [PubMed]
  8. Orlich, M.; Singh, P.; Sabaté, J.; Jaceldo-Siegl, K.; Fan, J.; Knutsen, S.; Beeson, W.; Fraser, G. Vegetarian Dietary Patterns and Mortality in Adventist Health Study 2. JAMA Intern. Med. 2013, 173, 1230–1238. [Google Scholar] [CrossRef]
  9. Kwok, C.; Umar, S.; Myint, P.; Mamas, M.; Loke, Y. Vegetarian Diet, Seventh Day Adventists and Risk of Cardiovascular Mortality: A Systematic Review and Meta-Analysis. Int. J. Cardiol. 2014, 176, 680–686. [Google Scholar] [CrossRef]
  10. Martínez-González, M.A.; Sanchez-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 Prevencion con Dieta Mediterranea (PREDIMED) study. Am. J. Clin. Nutr. 2014, 100, 320S–328S. [Google Scholar] [CrossRef]
  11. Chen, Z.; Zuurmond, M.; van der Schaft, N.; Nano, J.; Wijnhoven, H.; Ikram, M.; Franco, O.; Voortman, T. Plant Versus Animal Based Diets and Insulin Resistance, Prediabetes and Type 2 Diabetes: The Rotterdam Study. Eur. J. Epidemiol. 2018, 33, 883–893. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, Z.; Schoufour, J.; Rivadeneira, F.; Lamballais, S.; Ikram, M.; Franco, O.; Voortman, T. Plant-Based Diet and Adiposity Over Time in A Middle-Aged and Elderly Population. Epidemiology 2018, 30, 303–310. [Google Scholar] [CrossRef] [PubMed]
  13. Turner-McGrievy, G.M.; Mandes, T.; Crimarco, A. A plant-based diet for overweight and obesity prevention and treatment. J. Geriatr. Cardiol. 2017, 14, 369–374. [Google Scholar] [PubMed]
  14. Huang, R.Y.; Huang, C.C.; Hu, F.B.; Chavarro, J.E. Vegetarian diets and weight reduction: A meta-analysis of randomized controlled trials. J. Gen. Intern. Med. 2016, 31, 109–116. [Google Scholar] [CrossRef] [PubMed]
  15. Turner-McGrievy, G.M.; Davidson, C.R.; Wingard, E.E.; Wilcox, S.; Frongillo, E.A. Comparative effectiveness of plant-based diets for weight loss: A randomized controlled trial of five different diets. Nutrition 2015, 31, 350–358. [Google Scholar] [CrossRef] [PubMed]
  16. Tonstad, S.; Butler, T.; Yan, R.; Fraser, G.E. Type of vegetarian diet, body weight, and prevalence of type 2 diabetes. Diabetes Care 2009, 32, 791–796. [Google Scholar] [CrossRef] [PubMed]
  17. Rosell, M.; Appleby, P.; Spencer, E.; Key, T. Weight gain over 5 years in 21 966 meat-eating, fish-eating, vegetarian, and vegan men and women in EPIC-Oxford. Int. J. Obes. 2006, 30, 1389–1396. [Google Scholar] [CrossRef] [Green Version]
  18. Newby, P.; Tucker, K.L.; Wolk, A. Risk of overweight and obesity among semivegetarian, lactovegetarian, and vegan women. Am. J. Clin. Nutr. 2005, 81, 1267–1274. [Google Scholar] [CrossRef]
  19. Spencer, E.A.; Appleby, P.N.; Davey, G.K.; Key, T.J. Diet and Body Mass Index in 38 000 EPIC-Oxford Meat-Eaters, Fish-Eaters, Vegetarians and Vegans. Int. J. Obes. 2003, 27, 728–734. [Google Scholar] [CrossRef]
  20. Schulz, M.; Nöthlings, U.; Hoffmann, K.; Bergmann, M.M.; Boeing, H. Identification of a food pattern characterized by high-fiber and low-fat food choices associated with low prospective weight change in the EPIC-Potsdam cohort. J. Nutr. 2005, 135, 1183–1189. [Google Scholar] [CrossRef]
  21. Hruby, A.; Manson, J.E.; Qi, L.; Malik, V.S.; Rimm, E.B.; Sun, Q.; Willett, W.C.; Hu, F.B. Determinants and consequences of obesity. Am. J. Public Health 2016, 106, 1656–1662. [Google Scholar] [CrossRef]
  22. Foreman, K.; Marquez, N.; Dolgert, A.; Fukutaki, K.; Fullman, N.; McGaughey, M.; Pletcher, M.; Smith, A.; Tang, K.; Yuan, C.; et al. Forecasting Life Expectancy, Years of Life Lost, And All-Cause and Cause-Specific Mortality For 250 Causes of Death: Reference and Alternative Scenarios For 2016–40 For 195 Countries and Territories. Lancet 2018, 392, 2052–2090. [Google Scholar] [CrossRef]
  23. Satija, A.; Bhupathiraju, S.N.; Rimm, E.B.; Spiegelman, D.; Chiuve, S.E.; Borgi, L.; Willett, W.C.; Manson, J.E.; Sun, Q.; Hu, F.B. Plant-Based Dietary Patterns and Incidence of Type 2 Diabetes in US Men and Women: Results from Three Prospective Cohort Studies. PLoS Med. 2016, 13, e1002039. [Google Scholar] [CrossRef] [PubMed]
  24. Satija, A.; Bhupathiraju, S.; Spiegelman, D.; Chiuve, S.; Manson, J.; Willett, W.; Rexrode, K.; Rimm, E.; Hu, F.B. Healthful and Unhealthful Plant-Based Diets and The Risk of Coronary Heart Disease In, U.S. Adults. J. Am. Coll. Cardiol. 2017, 70, 411–422. [Google Scholar] [CrossRef] [PubMed]
  25. Satija, A.; Malik, V.; Rimm, E.B.; Sacks, F.; Willett, W.; Hu, F.B. Changes in intake of plant-based diets and weight change: Results from 3 prospective cohort studies. Am. J. Clin. Nutr. 2019. [Google Scholar] [CrossRef] [PubMed]
  26. Kim, H.; Caulfield, L.E.; Rebholz, C.M. Healthy Plant-Based Diets Are Associated with Lower Risk of All-Cause Mortality in US Adults. J. Nutr. 2018, 148, 624–631. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kim, H.; Caulfield, L.E.; Garcia-Larsen, V.; Steffen, L.M.; Grams, M.E.; Coresh, J.; Rebholz, C.M. Plant-Based Diets and Incident CKD and Kidney Function. Clin. J. Am. Soc. Nephrol. 2019, 14, 682–691. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. McKeown, N.; Jacobs, D. In Defence Of Phytochemical-Rich Dietary Patterns. Br. J. Nutr. 2010, 104, 1–3. [Google Scholar] [CrossRef]
  29. Liang, N.; Kitts, D.D. Role of Chlorogenic Acids in Controlling Oxidative and Inflammatory Stress Conditions. Nutrients 2015, 8, 16. [Google Scholar] [CrossRef]
  30. Cuevas-Sierra, A.; Ramos-Lopez, O.; Riezu-Boj, J.I.; Milagro, F.I.; Martinez, J.A. Diet, Gut Microbiota, and Obesity: Links with Host Genetics and Epigenetics and Potential Applications. Adv. Nutr. 2019, 10, S17–S30. [Google Scholar] [Green Version]
  31. Song, M.; Fung, T.; Hu, F.; Willett, W.; Longo, V.; Chan, A.; Giovannucci, E. Association of Animal and Plant Protein Intake with All-Cause and Cause-Specific Mortality. JAMA Intern. Med. 2016, 176, 1453–1463. [Google Scholar] [CrossRef] [PubMed]
  32. Westerterp-Plantenga, M.; Nieuwenhuizen, A.; Tomé, D.; Soenen, S.; Westerterp, K. Dietary Protein, Weight Loss, and Weight Maintenance. Annu. Rev. Nutr. 2009, 29, 21–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Lopez-Garcia, E.; Schulze, M.; Fung, T.; Meigs, J.; Rifai, N.; Manson, J.; Hu, F. Major Dietary Patterns Are Related to Plasma Concentrations of Markers of Inflammation and Endothelial Dysfunction. Am. J. Clin. Nutr. 2004, 80, 1029–1035. [Google Scholar] [CrossRef] [PubMed]
  34. Watzl, B.; Kulling, S.; Möseneder, J.; Barth, S.; Bub, A. A 4-Wk Intervention with High Intake of Carotenoid-Rich Vegetables and Fruit Reduces Plasma C-Reactive Protein in Healthy, Nonsmoking Men. Am. J. Clin. Nutr. 2005, 82, 1052–1058. [Google Scholar] [CrossRef] [PubMed]
  35. Ma, Y.; Griffith, J.A.; Chasan-Taber, L.; Olendzki, B.C.; Jackson, E.; Stanek, E.J.; Li, W.; Pagoto, S.L.; Hafner, A.R.; Ockene, I.S. Association between dietary fiber and serum C-reactive protein. Am. J. Clin. Nutr. 2006, 83, 760–766. [Google Scholar] [CrossRef]
  36. Hermsdorff, H.H.; Zulet, M.; Puchau, B.; Martínez, J.A. Fruit and vegetable consumption and proinflammatory gene expression from peripheral blood mononuclear cells in young adults: A translational study. Nutr. Metab. 2010, 7, 42. [Google Scholar] [CrossRef] [PubMed]
  37. Lopez-Legarrea, P.; Iglesia, R.D.L.; Abete, I.; Navas-Carretero, S.; Martinez, J.A.; Zulet, M.A. The protein type within a hypocaloric diet affects obesity-related inflammation: The RESMENA project. Nutrition 2014, 30, 424–429. [Google Scholar] [CrossRef]
  38. Carlos, S.; 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]
  39. De la Fuente-Arrillaga, C.; Ruiz, Z.V.; 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]
  40. 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]
  41. Mataix, J. Tabla de Composición de Alimentos, 5th ed.; Universidad de Granada: Granada, Spain, 2009. [Google Scholar]
  42. Moreiras, O.; Carbajal, A.; Cabrera, L.; Cuadrado, C. Tablas de Composición de Los Alimentos; de Prácticas, G., Ed.; Pirámide: Madrid, Spain, 2015. [Google Scholar]
  43. 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, 2595–2596. [Google Scholar] [CrossRef] [PubMed]
  44. 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]
  45. Zazpe, I.; Sánchez-Taínta, A.; Santiago, S.; Fuente-Arrillaga, C.D.L.; Bes-Rastrollo, M.; Martínez, J.A.; Martínez-González, M.A. 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]
  46. Sanchez-Tainta, A.; Zazpe, I.; Bes-Rastrollo, M.; Salas-Salvado, J.; Bullo, M.; Sorli, 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]
  47. Bes-Rastrollo, M.; Pérez, J.R.; Sánchez-Villegas, A.; Alonso, A.; Martínez-González, M.A. Validation of self-reported weight and body mass index in a cohort of university graduates in Spain. Rev. Esp. Obes. 2005, 3, 352–358. [Google Scholar]
  48. Mozaffarian, D. Dietary and Policy Priorities for Cardiovascular Disease, Diabetes, and Obesity: A Comprehensive Review. Circulation 2016, 133, 187–225. [Google Scholar] [CrossRef] [PubMed]
  49. Monteiro, C.A.; Cannon, G.; Moubarac, J.C.; Levy, R.B.; Louzada, M.L.C.; Jaime, P.C. The UN Decade of Nutrition, the NOVA food classification and the trouble with ultra-processing. Public Health Nutr. 2018, 21, 5–17. [Google Scholar] [CrossRef]
  50. Silveira, J.A.C.D.; Meneses, S.S.; Quintana, P.T.; Santos, V.D.S. Association between Overweight and Consumption of Ultra-Processed Food and Sugar-Sweetened Beverages among Vegetarians. Rev. Nutr. 2017, 30, 431–441. [Google Scholar] [CrossRef]
  51. Bertoia, M.L.; Rimm, E.B.; Mukamal, K.J.; Hu, F.B.; Willett, W.C.; Cassidy, A. Dietary Flavonoid Intake and Weight Maintenance: Three Prospective Cohorts of 124 086 US Men and Women Followed for up to 24 Years. BMJ 2016, 352, i17. [Google Scholar] [CrossRef]
  52. Santiago, S.; Zazpe, I.; Bes-Rastrollo, M.; Sánchez-Tainta, A.; Sayón-Orea, C.; Fuente-Arrillaga, C.D.L.; 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. 2014, 69, 297–302. [Google Scholar] [CrossRef]
  53. Barbaresko, J.; Koch, M.; Schulze, M.B.; Nöthlings, U. Dietary Pattern Analysis and Biomarkers of Low-Grade Inflammation: A Systematic Literature Review. Nutr. Rev. 2013, 71, 511–527. [Google Scholar] [CrossRef] [PubMed]
  54. Ramallal, R.; Toledo, E.; Martínez, J.A.; Shivappa, N.; Hébert, J.R.; Martínez-González, M.A.; Ruiz-Canela, M. Inflammatory Potential of Diet, Weight Gain, and Incidence of Overweight/Obesity: The SUN Cohort. Obesity 2017, 25, 997–1005. [Google Scholar] [CrossRef] [PubMed]
  55. Fresán, U.; Martínez-González, M.A.; Sabaté, J.; Bes-Rastrollo, M. Global Sustainability (Health, Environment and Monetary Costs) of Three Dietary Patterns: Results from a Spanish Cohort (the SUN Project). BMJ Open 2019, 9, e021541. [Google Scholar] [CrossRef] [PubMed]
  56. Hekler, E.B.; Gardner, C.D.; Robinson, T.N. Effects of a College Course About Food and Society on Students Eating Behaviors. Am. J. Prev. Med. 2010, 38, 543–547. [Google Scholar] [CrossRef] [PubMed]
  57. Swinburn, B.A.; Kraak, V.I.; Allender, S.; Atkins, V.J.; Baker, P.I.; Bogard, J.R.; Brinsden, H.; Calvillo, A.; Schutter, O.D.; Devarajan, R.; et al. The Global Syndemic of Obesity, Undernutrition, and Climate Change: The Lancet Commission Report. Lancet 2019, 393, 791–846. [Google Scholar] [CrossRef]
  58. Rothman, K.J. Epidemiology: An Introduction, 2nd ed.; Oxford University Press: New York, NY, USA, 2012. [Google Scholar]
  59. Rothman, K.J.; Gallacher, J.E.; Hatch, E.E. Why representativeness should be avoided. Int. J. Epidemiol. 2013, 42, 1012–1014. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flow chart depicting the selection process among participants of the Seguimiento Universidad de Navarra (SUN) cohort.
Figure 1. Flow chart depicting the selection process among participants of the Seguimiento Universidad de Navarra (SUN) cohort.
Nutrients 11 01553 g001
Figure 2. Hazard ratios (HR) and 95% confidence intervals (CI) of overweight/obesity according to baseline quintiles of the provegetarian food patterns.
Figure 2. Hazard ratios (HR) and 95% confidence intervals (CI) of overweight/obesity according to baseline quintiles of the provegetarian food patterns.
Nutrients 11 01553 g002
Figure 3. Sensitivity Analyses (results under different scenarios). Hazard Ratios (95% CI) of incident overweight/obesity for extreme quintiles (Q5 vs. Q1) of adherence to the overall, healthful and unhealthful provegetarian food patterns.
Figure 3. Sensitivity Analyses (results under different scenarios). Hazard Ratios (95% CI) of incident overweight/obesity for extreme quintiles (Q5 vs. Q1) of adherence to the overall, healthful and unhealthful provegetarian food patterns.
Nutrients 11 01553 g003
Table 1. Scoring criteria for the provegetarian food patterns.
Table 1. Scoring criteria for the provegetarian food patterns.
Provegetarian Food Pattern (Potential Range of 12–60)Healthful/Unhealthful Provegetarian Food Patterns (Potential Range of 18–90)
ComponentCriteriaComponentCriteria
Plant food GroupsEnergy-Adjusted QuintilesPlant Food GroupsEnergy-Adjusted Quintiles
HealthyHealthfulUnhealthful
1. VegetablesPositive1. VegetablesPositiveReverse
2. FruitsPositive2. FruitsPositiveReverse
3. LegumesPositive3. LegumesPositiveReverse
4. Cereal grainsPositive4. Whole grainsPositiveReverse
5. PotatoesPositive5. NutsPositiveReverse
6. NutsPositive6. Olive oilPositiveReverse
7. Olive oilPositive7. CoffeePositiveReverse
Less-healthy
Animal Food Groups 8. Fruit juicesReversePositive
8. DairyReverse9. PotatoesReversePositive
9. EggsReverse10. Refined grainsReversePositive
10. MeatReverse11. Sugary beveragesReversePositive
11. Fish & seafoodReverse12. PastriesReversePositive
12. Animal fatReverseAnimal Food Groups
13. DairyReverseReverse
14. EggsReverseReverse
15. MeatReverseReverse
16. Fish & seafoodReverseReverse
17. Miscellaneous foodReverseReverse
18. Animal fatReverseReverse
Table 2. Risk of overweight/obesity (HR and 95% CI) 1 according to quintiles (Q) of different provegetarian food patterns (FP).
Table 2. Risk of overweight/obesity (HR and 95% CI) 1 according to quintiles (Q) of different provegetarian food patterns (FP).
Provegetarian FPQ1Q2Q3Q4Q5p-Trendz-Continuous 2
Cases/Person-years622/29,366494/24,015358/17,697429/20,395417/21,740 2320/113,213
Age-and sex-adjusted1 (ref)0.97 (0.86 to 1.10)0.93 (0.82 to 1.06)0.94 (0.83 to 1.06)0.77 (0.68 to 0.88)<0.0010.83 (0.77 to 0.90)
Repeated measures1 (ref)0.98 (0.87 to 1.10)0.92 (0.81 to 1.05)0.94 (0.83 to 1.07)0.76 (0.67 to 0.86)<0.0010.82 (0.75 to 0.88)
Multivariable-adjusted 11 (ref)0.97 (0.85 to 1.09)0.95 (0.83 to 1.09)0.93 (0. 82 to1.06)0.85 (0.75 to 0.96)0.0140.89 (0.82 to 0.97)
Repeated measures1 (ref)0.98 (0.87 to 1.11)0.94 (0.83 to 1.08)0.94 (0. 83 to1.06)0.83 (0.73 to 0.95)0.0050.88 (0.81 to 0.95)
Healthful Provegetarian FP
Cases/Person-years490/24,57652
8/24,878
413/19,903516/25,293373/18,563 2320/113,213
Age-and sex-adjusted1 (ref)0.99 (0.87 to 1.12)0.95 (0.84 to 1.09)0.89 (0.79 to 1.01)0.81 (0.71 to 0.93)0.0010.87 (0.80 to 0.94)
Repeated measures1 (ref)0.99 (0.88 to 1.12)0.94 (0.82 to 1.07)0.92 (0.81 to 1.04)0.80 (0.70 to 0.92)0.0010.87 (0.80 to 0.94)
Multivariable-adjusted 11 (ref)0.98 (0.86 to 1.12)0.87 (0.76 to 1.00)0.87 (0.76 to 0.99)0.78 (0.67 to 0.90)<0.0010.84 (0.77 to 0.92)
Repeated measures1 (ref)0.98 (0.86 to 1.11)0.85 (0.74 to 0.97)0.89 (0.78 to 1.01)0.76 (0.66 to 0.88)<0.0010.84 (0.77 to 0.92)
Unhealthful Provegetarian FP
Cases/Person-years544/23,243516/24,135514/25,004388/20,549358/20,282 2320/113,213
Age-and sex-adjusted1 (ref)0.94 (0.83 to 1.06)0.93 (0.82 to 1.05)0.88 (0.77 to 1.00)0.84 (0.73 to 0.96)0.0060.88 (0.81 to 0.96)
Repeated measures1 (ref)0.95 (0.85 to 1.08)0.92 (0.81 to 1.04)0.87 (0.76 to 0.99)0.86 (0.75 to 0.99)0.0120.88 (0.81 to 0.96)
Multivariable-adjusted 11 (ref)1.03 (0.91 to 1.17)1.04 (0.92 to 1.18)1.00 (0.87 to 1.15)1.07 (0.92 to 1.23)0.5511.00 (0.91 to 1.10)
Repeated measures1 (ref)1.04 (0.92 to 1.17)1.03 (0.91 to 1.16)0.99 (0.86 to 1.13)1.09 (0.94 to 1.26)0.4161.01 (0.92 to 1.10)
Age was used as the underlying time variable in all the models. All the models were stratified by age groups (10-year periods) and year of recruitment (4-year periods). 1 Additionally adjusted for baseline body mass index (BMI, kg/m2, continuous), physical activity (metabolic equivalents (METs)-h/week, quartiles), hours of TV watching (quartiles), smoking status (current, never, former), marital status (single, married, other), years of university education (years, continuous), total energy intake (kcal/day, continuous), snacking between meals (yes, no), following a special diet at baseline (yes, no), parental family history of obesity (yes, no), hours of siesta (0, >0 to ≤0.5, >0.5). Robust standard errors were used. 2 Standardized z-scores divided by 2 (i.e., HR and 95% CI for each additional 2 z-score units).
Table 3. Estimates (differences and 95% confidence intervals) for average yearly weight change (g/y) according to updated quintiles (Q) of adherence to a Provegetarian Food Pattern (FP), a Healthful Provegetarian FP and an Unhealthful Provegetarian FP in the SUN Project.
Table 3. Estimates (differences and 95% confidence intervals) for average yearly weight change (g/y) according to updated quintiles (Q) of adherence to a Provegetarian Food Pattern (FP), a Healthful Provegetarian FP and an Unhealthful Provegetarian FP in the SUN Project.
Provegetarian FP
Q1Q2Q3Q4Q5p-Trend
Absolute Yearly Weight Change (g), Adjusted Mean 1 401 (343 to 460)413 (344 to 481)343 (280 to 405)274 (212 to 336)
Age-and Sex-Adjusted Differences Versus Q10 (ref)−52 (−132 to 27)-43 (-131 to 44)−112 (−195 to −30)−164 (−246 to −81)<0.001
Multivariable-Adjusted Differences Versus Q1 10 (ref)−39 (−118 to 39)−28 (−115 to 58)−98 (−181 to −16)−167 (−249 to −84)<0.001
Updated Provegetarian FP (repeated measures)
Absolute Yearly Weight Change (g), Adjusted Mean 1367 (321 to 413)361 (310 to 411)342 (282 to 402)302 (251 to 352)221 (168 to 274)
Multivariable-Adjusted Differences Versus Q1 (GEE) 20 (ref)−6 (−70 to 57)−26 (−96 to 45)−65 (−133 to 2)−146 (−214 to -78)<0.001
Healthful Provegetarian FP
Absolute Yearly Weight Change (g), Adjusted Mean 1472 (413 to 531)416 (359 to 473)347 (283 to 411)355 (297 to 411)270 (203 to 338)
Age-and Sex-Adjusted Differences Versus Q10 (ref)−69 (−151 to 13)−134 (−221 to −46)−123 (−206 to −40)−198 (−289 to −107)<0.001
Multivariable-adjusted differences versus Q1 10 (ref)−56 (−138 to 26)−125 (−212 to −38)−117 (−200 to −35)−202 (−294 to −110)<0.001
Updated Healthful Provegetarian FP (repeated measures)
Absolute Yearly Weight Change (g), Adjusted Mean 1416 (367 to 465)342 (295 to 390)285 (233 to 337)321 (273 to 370)216 (160 to 271)
Multivariable-Adjusted Differences versus Q1 (GEE) 20 (ref)−73 (−139 to−7)−130 (−201 to−61)−94 (−163 to −26)−200 (−275 to −125)<0.001
Unhealthful Provegetarian FP
Absolute Yearly Weight Change (g), Adjusted Mean 1365 (306 to 425)346 (288 to 404)401 (343 to 458)393 (330 to 457)387 (320 to 455)
Age-and Sex-Adjusted Differences Versus Q10 (ref)−11 (−93 to 71)50 (−32 to 133) 38 (−49 to 124)51 (−37 to 139)0.149
Multivariable-Adjusted Differences Versus Q1 10 (ref)−19 (−101 to 62)35 (−47 to 118)28 (−60 to 116)22 (−70 to 114)0.434
Updated Unhealthful Provegetarian FP (repeated measures)
Absolute Yearly Weight Change (g), Adjusted Mean 1306 (256 to 355)303 (254 to 351)334 (287 to 381)330 (278 to 382)338 (283 to 393)
Multivariable-Adjusted Differences Versus Q1 (GEE) 20 (ref)−3 (−70 to 64)28 (−39 to 96)24 (−47 to 96)32 (−43 to 108)0.324
1 Adjusted for sex, age, baseline BMI, physical activity, hours of TV watching, smoking status, marital status, years of university education, total energy intake, snacking between meals, following a special diet at baseline, parental family history of obesity, hours of siesta sleep and year of recruitment. 2 GEE, generalized estimating equations for multivariable-adjusted differences (95% CIs) of yearly weight change (g/y) using repeated measures of diet after 10 years of follow-up.
Table 4. Risk of obesity (HR and 95% CI) 1 in overweight participants (n = 3935) according to quintiles of different provegetarian food patterns (FP).
Table 4. Risk of obesity (HR and 95% CI) 1 in overweight participants (n = 3935) according to quintiles of different provegetarian food patterns (FP).
Provegetarian FPQ1Q2Q3Q4Q5p-Trendz-Continuous 2
 Cases/Person-years178/9650162/817896/5655126/7693107/8288 2320/113,213
 Age-and sex-adjusted1 (ref)1.06 (0.86 to 1.32)0.93 (0.72 to 1.20)0.90 (0.72 to 1.14)0.71 (0.55 to 0.90)0.0030.78 (0.67 to 0.91)
  Repeated measures1 (ref)1.06 (0.86 to 1.32)0.86 (0.67 to 1.12)0.92 (0.73 to 1.16)0.74 (0.58 to 0.94)0.0070.80 (0.68 to 0.93)
 Multivariable-adjusted 11 (ref)1.09 (0.87 to 1.36)0.91 (0.70 to 1.18)0.96 (0.75 to 1.23)0.79 (0.61 to 1.02)0.0560.82 (0.69 to 0.98)
  Repeated measures1 (ref)1.03 (0.83 to 1.29)0.83 (0.64 to 1.08)0.97 (0.77 to 1.24)0.83 (0.65 to 1.07)0.1230.85 (0.72 to 1.01)
Healthful Provegetarian FP
 Cases/Person-years141/7423139/7876135/6808149/9846105/7510 2320/113,213
 Age-and sex-adjusted1 (ref)0.95 (0.75 to 1.20)1.04 (0.82 to 1.32)0.80 (0.63 to 1.01)0.73 (0.56 to 0.95)0.0060.79 (0.68 to 0.93)
  Repeated measures1 (ref)0.98 (0.77 to 1.23)1.05 (0.83 to 1.33)0.78 (0.62 to 0.99)0.74 (0.57 to 0.96)0.0040.79 (0.68 to 0.92)
 Multivariable-adjusted 11 (ref)0.93 (0.73 to 1.20)1.03 (0.80 to 1.33)0.86 (0.67 to 1.10)0.79 (0.59 to 1.05)0.0660.87 (0.73 to 1.04)
  Repeated measures1 (ref)0.94 (0.74 to 1.21)1.04 (0.81 to 1.34)0.82 (0.64 to 1.05)0.79 (0.60 to 1.03)0.0360.86 (0.73 to 1.02)
Unhealthful Provegetarian FP
 Cases/Person-years179/9930143/8761144/8338108/644395/5991 2320/113,213
 Age-and sex-adjusted1 (ref)0.93 (0.74 to 1.16)0.99 (0.79 to 1.24)0.93 (0.73 to 1.19)0.87 (0.68 to 1.11)0.3050.92 (0.78 to 1.07)
  Repeated measures1 (ref)0.92 (0.74 to 1.14)0.97 (0.78 to 1.21)0.90 (0.70 to 1.15)0.86 (0.67 to 1.11)0.2560.92 (0.78 to 1.08)
 Multivariable-adjusted 11 (ref)1.12 (0.90 to 1.40)1.05 (0.83 to 1.32)1.03 (0.80 to 1.33)0.96 (0.73 to 1.26)0.7630.96 (0.81 to 1.14)
  Repeated measures1 (ref)1.12 (0.90 to 1.39)1.04 (0.83 to 1.31)1.02 (0.79 to 1.30)0.95 (0.73 to 1.25)0.6980.97 (0.82 to 1.15)
Age was used as the underlying time variable in all the models. All the models were stratified by age groups (10-year periods) and year of recruitment (4-year periods). 1 Additionally adjusted for baseline BMI (kg/m2, continuous), physical activity (METs-h/week, quartiles), hours of TV watching (quartiles), smoking status (current, never, former), marital status (single, married, other), years of university education (years, continuous), total energy intake (kcal/day, continuous), snacking between meals (yes, no), following a special diet at baseline (yes, no), parental family history of obesity (yes, no), hours of siesta (0, >0 to ≤0.5, >0.5). Robust standard errors were used. 2 Standardized z-scores divided by 2 (i.e., HR and 95% CI for each additional 2 z-score units).

Share and Cite

MDPI and ACS Style

Gómez-Donoso, C.; Martínez-González, M.Á.; Martínez, J.A.; Gea, A.; Sanz-Serrano, J.; Perez-Cueto, F.J.A.; Bes-Rastrollo, M. A Provegetarian Food Pattern Emphasizing Preference for Healthy Plant-Derived Foods Reduces the Risk of Overweight/Obesity in the SUN Cohort. Nutrients 2019, 11, 1553. https://doi.org/10.3390/nu11071553

AMA Style

Gómez-Donoso C, Martínez-González MÁ, Martínez JA, Gea A, Sanz-Serrano J, Perez-Cueto FJA, Bes-Rastrollo M. A Provegetarian Food Pattern Emphasizing Preference for Healthy Plant-Derived Foods Reduces the Risk of Overweight/Obesity in the SUN Cohort. Nutrients. 2019; 11(7):1553. https://doi.org/10.3390/nu11071553

Chicago/Turabian Style

Gómez-Donoso, Clara, Miguel Ángel Martínez-González, J. Alfredo Martínez, Alfredo Gea, Julen Sanz-Serrano, Federico J. A. Perez-Cueto, and Maira Bes-Rastrollo. 2019. "A Provegetarian Food Pattern Emphasizing Preference for Healthy Plant-Derived Foods Reduces the Risk of Overweight/Obesity in the SUN Cohort" Nutrients 11, no. 7: 1553. https://doi.org/10.3390/nu11071553

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