Next Article in Journal
Latent TB Infection, Vitamin D Status and COVID-19 Severity in Mongolian Patients
Next Article in Special Issue
Cancer Metabolism as a Therapeutic Target and Review of Interventions
Previous Article in Journal
Implications of Dietary Guidance for Sport and Exercise
Previous Article in Special Issue
Inflammatory and Metabolic Biomarker Assessment in a Randomized Presurgical Trial of Curcumin and Anthocyanin Supplements in Patients with Colorectal Adenomas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pro-Vegetarian Food Patterns and Cancer Risk among Italians from the Moli-Sani Study Cohort

by
Claudia Francisca Martínez
1,2,
Augusto Di Castelnuovo
3,
Simona Costanzo
1,
Teresa Panzera
1,
Simona Esposito
1,
Chiara Cerletti
1,
Maria Benedetta Donati
1,
Giovanni de Gaetano
1,
Licia Iacoviello
1,4,*,
Marialaura Bonaccio
1 and
on behalf of the Moli-Sani Study Investigators
1
Department of Epidemiology and Prevention, IRCCS NEUROMED, 86077 Pozzilli, Italy
2
Population Health Research Center, National Institute of Public Health, Cuernavaca 62100, Mexico
3
Mediterranea Cardiocentro, 80122 Naples, Italy
4
Department of Medicine and Surgery, Research Center in Epidemiology and Preventive Medicine (EPIMED), University of Insubria, 21100 Varese, Italy
*
Author to whom correspondence should be addressed.
Moli-Sani Study Investigators are listed in the Supplementary Appendix in Supplementary Materials.
Nutrients 2023, 15(18), 3976; https://doi.org/10.3390/nu15183976
Submission received: 1 August 2023 / Revised: 6 September 2023 / Accepted: 11 September 2023 / Published: 14 September 2023

Abstract

:
Besides the Mediterranean diet, there is a paucity of studies examining plant-based diets in relation to cancer outcomes in Mediterranean populations. We analyzed 22,081 apparently cancer-free participants (mean age 55 ± 12 year) from the Moli-sani study (enrollment period 2005–2010; Italy). A general pro-vegetarian food pattern was computed by assigning positive or negative scores to plant- or animal-derived foods, respectively from a 188-item FFQ. A priori healthful or unhealthful pro-vegetarian food patterns distinguished between healthy plant foods (e.g., fruits, vegetables) and less-healthy plant foods (e.g., fruit juices, refined grains). Cancer incidence was defined as the earliest diagnosis of cancer from hospital discharge records over a median follow-up of 12.9 years. In multivariable-adjusted analyses, a general pro-vegetarian food pattern was associated with a lower rate of cancer incidence (HR = 0.85; 95%CI 0.75–0.97 for Q5 vs. Q1); no association was observed between the healthful or unhealthful pro-vegetarian food patterns and overall cancer incidence. A healthful pro-vegetarian pattern, however, was inversely associated with digestive cancer (HR = 0.76; 95%CI 0.58–0.99 for Q5 vs. Q1), while the unhealthful pro-vegetarian pattern was directly linked to respiratory cancer (HR = 1.68; 95%CI 1.06–2.68 for Q5 vs. Q1). Our findings in a Mediterranean population support the hypothesis that some, but not all pro-vegetarian diets, might prevent some cancers.

1. Introduction

Cancer represents a leading cause of premature death in Europe and is progressively surpassing cardiovascular diseases (CVD) as the principal cause of death [1]. However, it has been estimated that 30–50% of all cancer cases are preventable both by modifying or avoiding key risk factors (e.g., smoking, alcohol abuse) and through a healthy diet and lifestyle [2]. Regarding diet, the World Cancer Research Fund/American Institute for Cancer Research strongly recommends, as a key strategy for cancer prevention, adherence to a balanced nutritious diet rich in plant-derived foods (e.g., whole grain, vegetables, fruit, and beans), while limiting intakes of red and processed meat, fast foods, and other processed foods high in fats, starches, or sugars [3].
Several healthful dietary patterns are favorably associated with cancer outcomes. These are mostly plant-based diets characterized by a large amount of food (e.g., fruits, vegetables, nuts, and seeds), that are natural sources of bioactive compounds with potential anti-cancer activities [4].
To date, a number of studies have examined the potential advantages of a traditional Mediterranean diet against certain cancers, but evidence varies greatly according to different cancer sites, with moderate evidence being reached only for cancer mortality and colorectal cancer risk, while being rated as “low” or “very low” for other cancer subtypes (e.g., bladder, breast, gastric, prostate) [5]. Although mainly based on the consumption of plant-derived foods, the Mediterranean diet still includes moderate consumption of animal-derived foods, such as poultry, fish, and dairy products (mostly in the form of long-preservable cheese), as well as moderate amounts of wine with meals [6].
On the contrary, strictly vegetarian diets completely exclude meat, meat-derived foods, and other products obtained from animals [7]. They are reported to be favorably associated with several health outcomes [8], but their impact on cancer is still controversial.
Although the proportion of true vegetarians is relatively scarce in most populations [9], the concept of pro-vegetarian diets, emphasizing the consumption of plant-derived foods, but allowing limited amounts of animal-based products [10], has gained increasing attention over the last years.
To estimate adherence to a pro-vegetarian food pattern distinct from full vegetarianism, a general pro-vegetarian food pattern index was first developed, weighing both plant-derived and animal-derived foods [9]. Subsequently, two additional indices were developed to examine the quality of the plant-based foods by distinguishing a priori between healthful and unhealthful pro-vegetarian food patterns, based on existing literature data [11]. In fact, not all plant foods might be equally healthy due to their different nutritional composition and the extent of food processing. This led to the differentiation of the quality of plant-based diets, as a healthful plant-based diet prioritizes the consumption of, e.g., vegetables, fruit, and whole grains, and an unhealthful plant-based diet mainly comprises a large amount of fruit drinks and juices, potatoes (mostly in the form of French fries and chips), sugary beverages, and cereals made with refined grains [11].
In the last 10 years, several large cohort studies supported the health advantages of a plant-based/pro-vegetarian diet, mostly in relation to mortality [9,12,13,14,15,16] and to the incidence of major chronic diseases, including CVD [13,17] and type 2 diabetes [11]. Research examining the association of a priori-defined plant-based indices with risk of cancer (all or site specific) [18,19,20,21,22] or cancer death [12] is less robust and has yielded inconsistent findings. Also, most of the available evidence stems from US populations, with limited research in Mediterranean populations. Moreover, not all cohort studies distinguished the quality of plant foods in relation to cancer outcomes.
To fill this knowledge gap and contribute to clarifying the uncertain relationship between plant-based diets and cancer risk, we sought to examine the association of three previously defined pro-vegetarian food patterns, defined as general, healthful, and unhealthful, in relation to cancer incidence in a general population residing in a Southern Italian region.

2. Methods

2.1. Study Population

We analyzed data from the Moli-sani study, which is a prospective population-based cohort consisting of 24,325 men and women aged ≥35 years, who were randomly recruited from the general population of the Molise region in Southern Italy. The cohort was established in 2005–2010 with the main aim of investigating the contribution of environmental and genetic risk factors in the onset of major non-communicable diseases and mortality. Pregnancy at the time of recruitment, disturbances in mental or decision-making impairments, current poly-traumas or coma, or refusal to sign the informed consent were exclusion criteria. The study design was described previously [23]. From these analyses, we excluded those participants with a prevalent registered or self-reported malignant cancer, had missing data on diet, reported extreme energy intakes (<800 or >4000 kcal/d in men and <500 or >3500 kcal/d in women), their dietary or medical questionnaires were judged as unreliable by interviewers or were lost to follow-up. After exclusions, a total of 22,081 cancer-free participants were included in the analyses (Supplementary Figure S1). The Moli-sani study complies with the Declaration of Helsinki and was granted the approval of the Ethics Committee of the Catholic University in Rome, Italy, ID Prot. pdc. P.99 (A.931/03-138-04)/C.E./2004. Written informed consent was collected from all participants.

2.2. Dietary Data Collection and Computation of Three Different Pro-Vegetarian Food Patterns

At baseline visit, data on participants’ dietary intake in the preceding year were assessed by an interviewer-administered, semi-quantitative EPIC (European Prospective Investigation into Cancer and Nutrition) food frequency questionnaire (FFQ), which was validated and adapted to the Italian population [24]. The FFQ contained 14 sections (i.e., pasta/rice, soup, meat, fish, raw vegetables, cooked vegetables, eggs, sandwiches, salami, and other cured meats, cheese, fruit, bread/wine, milk/coffee/cakes, and herbs/spices) with 248 questions regarding 188 different foods [24]. Using specialized software, the frequency and quantity of each food item were linked to Italian Food Tables [25,26] to estimate daily energy intake and macro- and micro-nutrients.
We computed a general pro-vegetarian food pattern which assigned positive scores to plant-derived foods and reverse scores to animal-derived foods, as proposed by Martinez-Gonzalez et al. [9]. For the healthful and unhealthful pro-vegetarian food patterns, we used as reference the work by Satija et al. [11], and the modification proposed by Oncina et al., which evaluated only olive oil instead of all vegetable oils [21]. The three patterns were determined by using dietary information of 17 food groups (5 animal-derived and 12 plant-derived food groups), which are all described in Table 1. To compute the three pro-vegetarian food patterns, the consumption (g/d) of food groups was adjusted for total energy for men and women, respectively using the residual method. The energy-adjusted residuals were grouped according to quintiles. To obtain the general pro-vegetarian food pattern, quintile values of 7 plant-derived foods (i.e., fruits, legumes, whole grain products, refined grain products, potatoes, nuts and dried fruit, olive oil) and reverse quintile values of 5 animal-derived foods (i.e., meat and meat products, animal fats, eggs, fish and seafood, milk and dairy) were summed. The final score of the general pro-vegetarian food pattern potentially ranged from 12 (minimum adherence) to 60 (maximum adherence). To compute the healthful pro-vegetarian food pattern, positive scores were given to 7 healthy plant food groups (i.e., fruits, vegetables, legumes, whole grains, nuts and dried fruit, olive oil, tea and coffee), and reverse scores to 5 less-healthy plant food groups (i.e., refined grains, potatoes, fruit juices, sugar-sweetened beverages, sweets and desserts). For the unhealthful pro-vegetarian food pattern, less-healthy plant-derived foods were assigned positive scores, whereas reverse scores were given to the healthy plant-derived food groups and foods of animal origin.
The five animal-derived food groups (i.e., meat and meat products, animal fats, eggs, fish and seafood, and milk and dairy products) were scored negatively both in healthful and unhealthful pro-vegetarian food patterns, whose final scores potentially ranged from 17 to 85 points (Table 1).
Adherence to a traditional MD was measured through the a priori 9-item Mediterranean diet score (MDS) proposed by Trichopoulou et al. [27], which reflects the dietary habits of the olive tree growing throughout countries bordering the Mediterranean Sea basin.

2.3. Baseline Covariate Assessment

At study entry, participants completed an interview that included sociodemographic information, lifestyle, and clinical variables. Based on the highest qualification attained, education was categorized as (1) up to lower secondary school (approximately ≤8 years), (2) upper secondary school (>8–13 years), and (3) post-secondary education (>13 years). Personal history of cardiovascular disease (angina, myocardial infarction, revascularization procedures, peripheral artery diseases, and cerebrovascular events) was self-reported and confirmed by clinical records and therapy. Comorbidities, such as hypertension, hyperlipidemia, or diabetes, were defined if the participant reported having a treatment with disease-specific drugs. Regular use of low-dose aspirin for primary or secondary CVD prevention (no/yes) was recorded. Leisure-time physical activity for sport, walking, and gardening was assessed by a structured interviewer-administered questionnaire [28] and expressed as daily energy expenditure in metabolic equivalent task hours (MET-h/d) [29]. Height and weight were measured to estimate body mass index (BMI) that was classified as normal (≤25 kg/m2), overweight (>25 and <30 kg/m2), or obese (≥30 kg/m2). Participants’ smoking habit was categorized as never, current, or former smoker (reporting to have not smoked in the previous 12 months or more at the time of interview). Housing tenure was classified as rented, ownership of one dwelling, and ownership of more than one dwelling. Marital status was coded as married/live-in couple, separated/divorced, single, and widower. Additional covariates were considered for females only and included the menopausal status, use of oral contraceptives, and use of hormone replacement therapy.

2.4. Outcome Definition and Assessment

The outcome examined in this study was cancer incidence, which was defined as the earliest diagnosis of cancer from hospital discharge records. A hospitalization was defined as any length of stay of at least 24 h in a hospital or clinic. If a patient was transferred to another hospital or facility, it was considered a single hospitalization. Hospitalizations for rehabilitation and chemotherapy and/or radiotherapy (i.e., elective treatment in a hospital-based unit) were excluded. Primary and secondary diagnoses for hospitalization were coded using the 9th version of the International Classification of Diseases. A cancer hospitalization was assessed if the primary or the first secondary diagnoses of admission to the hospital were coded as ICD-9: 140–172, 174–208. The Moli-sani study cohort was followed-up for cancer hospitalization from study entry (2005–2010) through 31 December 2020.

3. Statistical Analysis

We quantified differences in the distribution of baseline covariates across fifths of a general pro-vegetarian food pattern using generalized linear models for both continuous and categorical variables adjusted for age, sex, and energy intake (kcal/d), (GENMOD procedure for categorical variables and GLM procedure for continuous variables using SAS software). Associations of three pro-vegetarian food patterns (both as continuous or categorical dependent variables) with the rate of cancer hospitalization were calculated as hazard ratios (HRs) and 95% confidence intervals (CIs) using Cox proportional hazard models. No violations of the proportional hazard assumption (log(−log) survival plots curves) were visually identified.
Potential confounders (measured at baseline) were defined a priori and identified based on existing literature, rather than deferring to statistical criteria [30], on the basis of their previously documented association both with diet and cancer risk [31,32,33]. The following multivariable models were ultimately fitted: Model 1 included age at baseline, sex, energy intake (kcal/d), and alcohol intake (g/d; continuous). Model 2, as in Model 1, further controlled for educational level, housing tenure, place of residence, smoking, body mass index (categorical), leisure-time physical activity, history of CVD, diabetes, hypertension, hyperlipidemia (no/yes), aspirin use, hormone replacement therapy, oral contraception, menopausal status, family history of cancer. Participants contributed person-time until the date of cancer hospitalization, death, date of emigration or loss to follow-up, or end of follow-up, whichever occurred first. We used cancer-free survival curves (adjusted for potential confounding variables as in Model 2) to describe cancer incidence over time across categories of adherence to pro-vegetarian food patterns. To test the robustness of the results, several sensitivity analyses were conducted by excluding (a) participants with a history of CVD, (b) participants with diabetes, and (c) cancer hospitalizations that occurred within the first 2 years of follow-up. Main analyses (Model 2) were re-run in men and women, and across age groups (<65 and ≥65 years), BMI categories, and by smoking status. Appropriate terms for testing interaction were included in the multivariable model to test for a difference of effect of the pro-vegetarian food patterns (1 SD increase) in relation to cancer risk within pre-specified population groups.
Missing data on covariates were handled using multiple imputation (SAS PROC MI, followed by PROCMIANALYZE), over different simulations (n = 10 imputed datasets) to avoid bias introduced by not-at-random missing data patterns. The data analysis was generated using SAS/STAT software, version 9.4 (SAS Institute Inc., Cary, NC, USA).

4. Results

Sociodemographic characteristics of the analyzed sample across fifths of adherence to a general pro-vegetarian food pattern are shown in Table 2. The median value of the general pro-vegetarian food pattern in this population was 36 (mean = 36; SD = 5.6; min-to-max = 14 to 57), median values for the healthful and unhealthful pro-vegetarian food patterns were 47 (mean = 47; SD = 6.0; min-to-max = 28 to 69) and 47 (mean = 47; SD = 6.5; min-to-max = 22 to 70), respectively. The general pro-vegetarian was directly associated with the healthful pro-vegetarian food pattern (Spearman correlation coefficient = 0.72; p-value <0.0001) and showed a weak inverse relationship with an unhealthful pro-vegetarian food pattern (Spearman correlation coefficient = −0.03; p-value <0.0001); the healthful and unhealthful pro-vegetarian food pattern were moderately correlated (Spearman correlation coefficient = −0.40; p-value <0.0001). The Spearman correlation coefficients of a general, healthful, and unhealthful pro-vegetarian food pattern with the Mediterranean diet score were 0.52, 0.49, and −0.30, respectively (all p-values <0.0001), suggesting weak to moderate strengths of these associations [34].
Subjects with higher adherence to the general pro-vegetarian food pattern tended to have higher adherence to the traditional Mediterranean diet, to be older, to have higher household income, and were more physically active. However, participants following a general pro-vegetarian food pattern were more likely to report a history of cardiovascular disease, hypertension, hyperlipidemia, and daily aspirin use than those with lower adherence to this food pattern (Table 2).
In terms of consumption, the general pro-vegetarian was directly associated with intake of plant-based food, and inversely with animal-based food consumption (Supplementary Table S1). Increasing adherence to the healthful pro-vegetarian food pattern was positively associated with total energy intake, dietary fiber, alcohol (g/d), and carbohydrate and fat (% of total energy intake; TEI); higher adherence to an unhealthful pro-vegetarian food pattern was linked to reduced total energy intake, dietary fiber, and fat intake (%TEI), and higher daily energy from carbohydrate and alcohol intake (g/d) (Supplementary Table S1). Over 12.9 years (median) of follow-up (interquartile ranges = 11.8 to 14.1 years; 271,891 person-years), we documented a total of 2306 hospital admissions for cancer (any site). Of these, 598 were digestive cancers (25.9% of all cancer hospitalizations; ICD-9 150–159); 183 (7.9%) respiratory and intrathoracic organs (ICD-9 160–165); 395 (17.1%) genitourinary organs (ICD-9 179–189); 285 (12.4%) breast cancer (ICD-9 174); 219 (9.5%) prostate (ICD-9 185); 178 (7.7%) lymphatic and hematopoietic tissue (ICD-9 200–209); 41 (1.8%) brain and nervous system (ICD-9 191–192); 407 (17.6%) were other cancer types. In multivariable-adjusted analyses controlled for known risk factors, higher adherence to the general pro-vegetarian food pattern was associated with 15% lower rate of cancer hospitalization (HR = 0.85; 95%CI 0.75–0.97 for Q5 (score > 42) vs. Q1 (score < 31); p-value for trend across fifths = 0.045). Similarly, an increment of 1 SD in the general pro-vegetarian food pattern was linked to a reduction of 5% in cancer rate (HR = 0.95; 95%CI 0.91–0.99) (Table 3). Further adjustment for the MDS did not alter these associations (HR = 0.84; 95%CI 0.72–0.97 for Q5 vs. Q1 of the general pro-vegetarian food pattern; HR = 0.94; 95%CI 0.90–0.99 for 1 SD increment in the pro-vegetarian food pattern). Multivariable-adjusted cancer-free survival curves across categories of the general pro-vegetarian food pattern were well separated with a tendency to diverge over time (Figure 1). Neither healthful nor unhealthful pro-vegetarian food patterns were associated with overall cancer risk (Table 3; Supplementary Figure S2).

Subgroup and Sensitivity Analyses

None of the baseline risk factors act as an effect modifier of the association between the pro-vegetarian food patterns (general, healthful, or unhealthful) with cancer risk (Table 4). Sensitivity analyses showed consistent inverse associations between a general pro-vegetarian food pattern with cancer risk in all the analyzed scenarios (excluding baseline CVD, diabetes, or early cancer hospitalizations) (Table 4). When hospitalizations were separately examined by specific cancer site, a general pro-vegetarian food pattern was associated with a lower rate of respiratory cancer (HR = 0.61; 95%CI 0.38–0.98 for Q5 vs. Q1), or digestive neoplasm (HR = 0.74; 95%CI 0.57–0.95 for Q5 vs. Q1) (Table 5). For the healthful pro-vegetarian food pattern, a lower rate of digestive cancer was also observed (HR = 0.76; 95%CI 0.58–0.99 for Q5 vs. Q1). Higher adherence to an unhealthful pro-vegetarian food pattern was otherwise directly associated with the rate of respiratory tract cancer (HR = 1.68; 95%CI 1.06–2.68 for Q5 vs. Q1) (Table 5).

5. Discussion

In this large prospective cohort of Italian men and women, who were apparently without cancer at baseline, a higher adherence to the general pro-vegetarian food pattern was associated with a lower rate of cancer incidence overall, specifically for digestive or respiratory cancers. Adherence to healthful or unhealthful pro-vegetarian food patterns was not associated with overall cancer risk. However, when specific cancer sites were considered, a healthful pro-vegetarian food pattern was linked to a reduced rate of digestive cancer, whereas an unhealthful pro-vegetarian food pattern was associated with a higher likelihood of respiratory cancers.
The relationships between plant-derived foods and diets with cancer risk have currently different levels of evidence, based on the food group and the cancer site. Probably, some key food groups that are positively scored in the general and healthful pro-vegetarian food patterns (e.g., dietary fiber and whole grains) could reduce the risk of colorectal cancer, while coffee consumption (that is positively scored in the healthful pro-vegetarian diet) has been associated with lower incidence of liver cancer [3]. Limited evidence exists on the potential protective effects of diets rich in fruits, vegetables, polyphenols, and vitamin C against respiratory cancer [3], while foods with high-glycemic index, such as sugary beverages, cereals made with refined grains, fruit juices and potatoes, that scored high in the unhealthful pro-vegetarian diet, could possibly increase the risk of lung cancer [35], although evidence is not conclusive [36].
Also, we acknowledge that although the total number of cancer hospitalizations was considerably large, the number of specific cancers was restricted; therefore, possibly limiting the statistical power of analyses for different cancer sites (e.g., brain cancer).
Of interest, in our study, the inverse association of a general pro-vegetarian food pattern with overall cancer risk was independent of adherence to the Mediterranean diet. This is possibly due to the fact that there are some important differences between the pro-vegetarian food pattern and its healthful version and the traditional Mediterranean diet, mostly in relation to the scoring criteria for three foods: potatoes, fish, and alcohol. This is further supported by the moderate correlation observed between pro-vegetarian and Mediterranean diet scores, indicating that these novel dietary indices might be able to capture novel unique relations of a pro-vegetarian diet with respect to cancer outcomes [37].
This is possibly the first study that has been conducted on Italians to examine three different pro-vegetarian food patterns in relation to cancer risk and acknowledge the quality of plant-derived foods.
Our findings on the general pro-vegetarian food pattern and cancer risk are in line with prior analyses in the large NutriNet-Santè cohort in France, showing a reduced risk of overall cancer for participants with strict adherence to a pro-vegetarian diet [18]. As shown in our study and in this large French cohort, a general pro-vegetarian food pattern was mostly effective in reducing the risk of digestive or respiratory cancers. This is possibly due to the high intake of fruits, vegetables, and fiber reported by participants with higher adherence to the overall pro-vegetarian food pattern, with a potential protective role against carcinogenesis [34,35]. In the UK Biobank study, a difference in overall cancer risk was observed between healthy and unhealthy plant-based diets, but not for individual cancers [21]. However, the potential effectiveness of a healthful plant-based diet against digestive cancers was documented in a previous hospital-based, multi-case control study [21], that showed a substantial reduction in the risk of oesophagus, stomach, and pancreatic cancers for participants with a very high adherence to this food pattern. Consistently, a large dietary share of unhealthful plant-based products was associated with a higher likelihood of oesophagus and stomach cancers, but not pancreatic neoplasm [21]. In our analyses, we did not observe any association of these three pro-vegetarian food patterns with breast or prostate cancers. In this study, our data are in contrast with analyses on a large sample of Spanish women from the SUN cohort, where a moderate adherence to a pro-vegetarian diet was associated with lower breast cancer incidence compared to lower consumption, and mostly amongst pre-menopausal women. In addition, a moderate conformity to an unhealthful pro-vegetarian food pattern was associated with a two-fold increased breast cancer risk, but only amongst post-menopausal women [19]. However, the lack of an association between a pro-vegetarian diet with breast cancer was reported amongst women participants in the NutriNet-Santè cohort from France [18].
Prospective analyses from the Health Professionals Follow-Up Study in US supported evidence that greater consumption of healthful plant-based foods was associated with a lower risk of prostate cancer, specifically aggressive forms among men aged <65 years [20]. However, analyses from the NutriNet-Santè cohort did not highlight a role of a general pro-vegetarian against prostatic neoplasm [18].
A relatively low number of cohort studies have analyzed plant-based diets in relation to cancer death, but the evidence is not conclusive. A prospective cohort study conducted in the US found an inverse association of an overall plant-based diet with cancer mortality, but the quality of plant foods was unrelated to the outcome [15]. Similarly, there were no meaningful associations between changes in plant-based diet indices and cancer mortality in a well-established large cohort of US men and women [12]. Finally, analyses on a large cohort of South Korean adults showed a higher risk of cancer death for participants with the highest adherence to a diet rich in unhealthful plant-derived foods, while overall or healthful plant-based diets seemed unrelated to the outcome [13].
There are several potential mechanisms through which greater adherence to a plant-based diet may influence cancer risk. High amounts of vegetables, fruits, and legumes increase the dietary intake of phytochemicals and vitamins with antioxidant and antiproliferative activity [12,38,39,40]. In addition, higher consumption of fiber and whole grains [41,42,43], and low consumption of sugar-sweetened beverages may play an important role in regulating mechanisms associated with glucose control and insulin growth factor, which are related to a higher risk of cancer [44,45,46,47,48]. In addition, fiber from a plant-based diet constitutes a beneficial substrate for intestinal microbiota and metabolites associated with certain cancers [49]. Also, both the general and healthful pro-vegetarian food patterns emphasize the consumption of olive oil, which has long been associated with lowering the risk of developing or dying from cancer, possibly through its anti-inflammatory, antioxidant, anti-atherogenic, and anti-thrombotic properties [50].
The health advantages of a pro-vegetarian diet could also be due to a low dietary share of animal-based foods; indeed, elevated intakes of some animal-derived foods, such as red and processed meat, have been reportedly associated with overall cancer risk [51,52] and certain cancers [52,53]. On the other hand, sugary beverages have been associated with increased incidence of and mortality from breast cancer among post-menopausal women [44,46,54]. Also, evidence suggests an increased risk of colon and gastric cancer associated with refined grains; however, findings remain unclear [52].
Other detrimental effects, aside from the poor quality of some unhealthful plant-based foods, could be due to non-nutritional factors, such as increased exposure to additives, toxic contaminants migrated from food packaging, and exposure to acrylamide, which is produced during cooking procedures at high-temperature and may increase cancer risk [55,56,57]. In fact, the unhealthful plant-based diets are also characterized by the presence of ultra-processed foods as defined by the Nova classification [58], and thus share all the critical features of these foods which go well beyond their poor nutritional composition [59]. Therefore, analyses (and interpretations) on the quality of plant-based diets in relation to health outcomes, including cancer, should not be limited to the evaluation of the nutritional content of plant-derived foods, but should consider the extent of food processing, as well.

6. Strengths and Limitations

The strengths of this study include the large sample size, a prospective population-based cohort design, a follow-up of more than 10 years, and the large number of covariates used to minimize, at least in part, confounding. Moreover, we computed three indices of pro-vegetarian food patterns that are easily reproducible for further comparisons. Also, we used sensitivity analyses to support the robustness of the results, by excluding participants who might have modified habitual dietary intakes due to illness.
Several limitations should be considered when interpreting our findings. The possibility of residual confounding remains, given the observational design, so that any causality relationship is only proposed. However, the main potential confounders were considered; therefore, it is unlikely that residual confounding entirely explains the observed results. Limitations in the outcome definition should be also acknowledged, since in this study, we only analyzed cancer hospitalizations, and consequently, cases of cancer that were managed in outpatient settings were missed.
Additionally, interpretation of risk estimations for specific cancer sites should consider the restricted number of events in subgroups, which may affect the statistical power and precision of CIs. Another weakness is the lack of repeated dietary assessments during follow-up, although dietary habits in adulthood tend to remain stable over time [60]. Moreover, dietary data were self-reported with an inherent risk of error and bias, possibly attenuated by exclusion of participants with implausible energy intakes. Furthermore, the FFQ used in the present study was extensively validated against diet records and biomarkers [61]. We also acknowledge that the association of each dietary component with cancer risk was not examined since diet was considered as a synergistic interaction of multiple components. Although based on prior evidence, the identification of healthful plant foods has an element of subjectivity [11]. Moreover, we could not differentiate potatoes by cooking method or processing level; therefore, all types of potatoes were scored as “unhealthy” in the healthful pro-vegetarian food pattern. The lack of correction for family-wise error is another limitation of this study. Finally, the generalizability of the results could be limited to an Italian population.

7. Conclusions

Our findings in a Mediterranean population suggest that, as already known for a traditional Mediterranean diet, a pro-vegetarian diet, based on healthful plant-derived foods and limited animal-derived products intakes, could be effective in reducing cancer incidence. Moreover, our data support prior evidence showing that a large dietary share of unhealthful, though plant-based foods, mostly highly processed, is associated with an increased risk of certain cancers [19,22] or cancer death [13,15,62].
In conclusion, a diet that contains high amounts of healthful plant products and minimal amounts of animal products should be encouraged, when a traditional Mediterranean diet cannot be followed. Future targeted analyses could be helpful to understand whether the observed health advantages associated with this pro-vegetarian dietary pattern are mostly due to the large amount of plant foods or to the lower consumption of animal foods or both. In addition, studies could attempt to establish the ideal ratios of high plant food consumption and low animal food consumption to maximize health benefits. In all instances, it is important in public health recommendations and dietary guidelines worldwide to discourage the consumption of unhealthful, especially ultra-processed, plant-derived products.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu15183976/s1, Figure S1: Flowchart for selection of study participants from the Moli-sani Study, 2005–2010.; Figure S2: Association between adherence to (a) healthful pro-vegetarian food pattern; and (b) unhealthful provegetarian food pattern with cancer risk in participants from the Moli-sani study cohort (n = 22,081).; Table S1: Scoring criteria for the pro-vegetarian (PVG) food patterns in the Moli-sani Study cohort, Italy 2005–2010; Table S2: Associations of pro-vegetarian (PVG) Food Patterns with nutritional factors in cancer-free participants from the Moli-sani Study cohort (n = 22,081).

Author Contributions

M.B. and C.F.M. contributed to the conception, design, and interpretation of data; S.C., T.P. and S.E. managed data collection; A.D.C. and M.B. analyzed the data; M.B. and C.F.M. wrote the manuscript; M.B.D., C.C., G.d.G. and L.I. originally inspired the Moli-sani study and critically reviewed the manuscript. Full approval of the manuscript by all authors was completed. All authors have read and agreed to the published version of the manuscript.

Funding

The present analyses were partially supported by the Italian Ministry of Health (Ricerca Corrente 2022–2024). Funders had no role in any phase of the study or in the analysis and interpretation of findings, nor in the preparation of the manuscript or in the publication process. All authors were and are independent of funders.

Institutional Review Board Statement

The Moli-sani study complies with the Declaration of Helsinki and was granted the approval of the Ethics Committee of the Catholic University in Rome, Italy, ID Prot. pdc. P.99 (A.931/03-138-04)/C.E./2004.

Informed Consent Statement

Written informed consent was obtained from all participants at study entry.

Data Availability Statement

The data underlying this article will be shared on reasonable request to the corresponding author. The data are stored in an institutional repository (https://repository.neuromed.it), and their access is restricted by ethical approvals and the legislation of the European Union.

Acknowledgments

We wish to thank the Pfizer Foundation (Rome, Italy), the Italian Ministry of University and Research (MIUR, Rome, Italy)– Decreto no.1588 Programma Triennale di Ricerca, and Instrumentation Laboratory, Milan, Italy for their support in the enrolment phase of the Moli-sani study. We are grateful to the Moli-sani study participants who enthusiastically joined the study and to the Associazione Cuore Sano ETS (Campobasso, Italy) for its cultural support of the research activities. A post-doctoral fellowship was granted to C.F.M. and the present study was performed in the context of the Fondazione Umberto Veronesi-IRCCS Neuromed framework agreement. E.R. was supported by the Fondazione Umberto Veronesi, which is gratefully acknowledged. C.F.M. gratefully recognizes the support allocated by the National Council of Humanities Science and Technology (México).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bray, F.; Laversanne, M.; Weiderpass, E.; Soerjomataram, I. The Ever-Increasing Importance of Cancer as a Leading Cause of Premature Death Worldwide. Cancer 2021, 127, 3029–3030. [Google Scholar] [CrossRef] [PubMed]
  2. World Health Organization. Cancer Prevention. Available online: https://www.who.int/health-topics/cancer#tab=tab_2 (accessed on 17 July 2023).
  3. World Cancer Research Fund. American Institute for Cancer Research. Diet, Nutrition, Physical Activity and Cancer: A Global Perspective A Summary of the Third Expert Report 2018. Available online: https://www.wcrf.org/wp-content/uploads/2021/02/Summary-of-Third-Expert-Report-2018.pdf (accessed on 17 July 2023).
  4. Subramaniam, S.; Selvaduray, K.R.; Radhakrishnan, A.K. Bioactive Compounds: Natural Defense Against Cancer? Biomolecules 2019, 9, 758. [Google Scholar] [CrossRef] [PubMed]
  5. Morze, J.; Danielewicz, A.; Przybyłowicz, K.; Zeng, H.; Hoffmann, G.; Schwingshackl, L. An Updated Systematic Review and Meta-Analysis on Adherence to Mediterranean Diet and Risk of Cancer. Eur. J. Nutr. 2021, 60, 1561–1586. [Google Scholar] [CrossRef] [PubMed]
  6. Trichopoulou, A.; Martínez-González, M.A.; Tong, T.Y.N.; Forouhi, N.G.; Khandelwal, S.; Prabhakaran, D.; Mozaffarian, D.; de Lorgeril, M. Definitions and Potential Health Benefits of the Mediterranean Diet: Views from Experts around the World. BMC Med. 2014, 12, 112. [Google Scholar] [CrossRef] [PubMed]
  7. Hargreaves, S.M.; Rosenfeld, D.L.; Moreira, A.V.B.; Zandonadi, R.P. Plant-Based and Vegetarian Diets: An Overview and Definition of These Dietary Patterns. Eur. J. Nutr. 2023, 62, 1109–1121. [Google Scholar] [CrossRef]
  8. Oussalah, A.; Levy, J.; Berthezène, C.; Alpers, D.H.; Guéant, J.L. Health Outcomes Associated with Vegetarian Diets: An Umbrella Review of Systematic Reviews and Meta-Analyses. Clin. Nutr. 2020, 39, 3283–3307. [Google Scholar] [CrossRef]
  9. 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]
  10. Hemler, E.C.; Hu, F.B. Plant-Based Diets for Personal, Population, and Planetary Health. Adv. Nutr. 2019, 10, S275–S283. [Google Scholar] [CrossRef]
  11. Satija, A.; Bhupathiraju, S.N.; Rimm, E.B.; Spiegelman, D.; Chiuve, S.E.; Borgi, L.; Willett, W.C.; Manson, J.A.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]
  12. Baden, M.Y.; Liu, G.; Satija, A.; Li, Y.; Sun, Q.; Fung, T.T.; Rimm, E.B.; Willett, W.C.; Hu, F.B.; Bhupathiraju, S.N. Changes in Plant-Based Diet Quality and Total and Cause-Specific Mortality. Circulation 2019, 140, 979–991. [Google Scholar] [CrossRef]
  13. Kim, J.; Kim, H.; Giovannucci, E.L. Plant-Based Diet Quality and the Risk of Total and Disease-Specific Mortality: A Population-Based Prospective Study. Clin. Nutr. 2021, 40, 5718–5725. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, H.; Shen, J.; Xuan, J.; Zhu, A.; Ji, J.S.; Liu, X.; Cao, Y.; Zong, G.; Zeng, Y.; Wang, X.; et al. Plant-Based Dietary Patterns in Relation to Mortality among Older Adults in China. Nat. Aging 2022, 2, 224–230. [Google Scholar] [CrossRef] [PubMed]
  15. Li, H.; Zeng, X.; Wang, Y.; Zhang, Z.; Zhu, Y.; Li, X.; Hu, A.; Zhao, Q.; Yang, W. A Prospective Study of Healthful and Unhealthful Plant-Based Diet and Risk of Overall and Cause-Specific Mortality. Eur. J. Nutr. 2022, 61, 387–398. [Google Scholar] [CrossRef] [PubMed]
  16. Shah, S.; Mahamat-Saleh, Y.; Ait-Hadad, W.; Koemel, N.A.; Varraso, R.; Boutron-Ruault, M.C.; Laouali, N. Long-Term Adherence to Healthful and Unhealthful Plant-Based Diets and Breast Cancer Risk Overall and by Hormone Receptor and Histologic Subtypes among Postmenopausal Females. Am. J. Clin. Nutr. 2023, 117, 467–476. [Google Scholar] [CrossRef]
  17. Weston, L.J.; Kim, H.; Talegawkar, S.A.; Tucker, K.L.; Correa, A.; Rebholz, C.M. Plant-Based Diets and Incident Cardiovascular Disease and All-Cause Mortality in African Americans: A Cohort Study. PLoS Med. 2022, 19, e1003863. [Google Scholar] [CrossRef]
  18. Kane-Diallo, A.; Srour, B.; Sellem, L.; Deschasaux, M.; Latino-Martel, P.; Hercberg, S.; Galan, P.; Fassier, P.; Guéraud, F.; Pierre, F.H.; et al. Association between a pro Plant-Based Dietary Score and Cancer Risk in the Prospective NutriNet-Santé Cohort. Int. J. Cancer 2018, 143, 2168–2176. [Google Scholar] [CrossRef] [PubMed]
  19. Romanos-Nanclares, A.; Toledo, E.; Sánchez-Bayona, R.; Sánchez-Quesada, C.; Martínez-González, M.Á.; Gea, A. Healthful and Unhealthful Provegetarian Food Patterns and the Incidence of Breast Cancer: Results from a Mediterranean Cohort. Nutrition 2020, 79–80, 110884. [Google Scholar] [CrossRef] [PubMed]
  20. Loeb, S.; Fu, B.C.; Bauer, S.R.; Pernar, C.H.; Chan, J.M.; Van Blarigan, E.L.; Giovannucci, E.L.; Kenfield, S.A.; Mucci, L.A. Association of Plant-Based Diet Index with Prostate Cancer Risk. Am. J. Clin. Nutr. 2022, 115, 662–670. [Google Scholar] [CrossRef]
  21. Oncina-Cánovas, A.; González-Palacios, S.; Notario-Barandiaran, L.; Torres-Collado, L.; Signes-Pastor, A.; de-Madaria, E.; Santibañez, M.; García-de la Hera, M.; Vioque, J. Adherence to Pro-Vegetarian Food Patterns and Risk of Oesophagus, Stomach, and Pancreas Cancers: A Multi Case–Control Study (The PANESOES Study). Nutrients 2022, 14, 5288. [Google Scholar] [CrossRef]
  22. Thompson, A.S.; Tresserra-Rimbau, A.; Karavasiloglou, N.; Jennings, A.; Cantwell, M.; Hill, C.; Perez-Cornago, A.; Bondonno, N.P.; Murphy, N.; Rohrmann, S.; et al. Association of Healthful Plant-Based Diet Adherence with Risk of Mortality and Major Chronic Diseases Among Adults in the UK. JAMA Netw. Open 2023, 6, e234714. [Google Scholar] [CrossRef]
  23. Iacoviello, L.; Bonanni, A.; Costanzo, S.; Curtis, A.; Castelnuovo, A.; Olivieri, M.; Zito, F.; Donati, M.B.; Gaetano, G. The Moli-Sani Project, a Randomized, Prospective Cohort Study in the Molise Region in Italy; Design, Rationale and Objectives. Ital. J. Public Health 2007, 4, 110–118. [Google Scholar]
  24. Pisani, P.; Faggiano, F.; Krogh, V.; Palli, D.; Vineis, P.; Berrino, F. Relative Validity and Reproducibility of a Food Frequency Dietary Questionnaire for Use in the Italian EPIC Centres. Int. J. Epidemiol. 1997, 26 (Suppl. S1), S152. [Google Scholar] [CrossRef] [PubMed]
  25. Salvini, S.; Parpinel, M.; Gnagnarella, P.; Maissoneuve, P.; Turrini, A. Banca Dati Tabella Di Composizione Degli Alimenti—Tabella Di Composizione Degli Alimenti—CREA. Available online: https://www.crea.gov.it/-/tabella-di-composizione-degli-alimenti (accessed on 24 November 2022).
  26. Pala, V.; Sieri, S.; Palli, D.; Salvini, S.; Berrino, F.; Bellegotti, M.; Frasca, G.; Tumino, R.; Sacerdote, C.; Fiorini, L.; et al. Diet in the Italian Epic Cohorts: Presentation of Data and Methodological Issues. Tumori J. 2018, 89, 594–607. [Google Scholar] [CrossRef] [PubMed]
  27. 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]
  28. Mannocci, A.; Thiene, D.D.; Cimmuto, A.D.; Masala, D.; Vito, E.; Torre, G. International Physical Activity Questionnaire: Validation and Assessment in an Italian Sample. Ital. J. Public. Health 2010, 7, 369–376. [Google Scholar] [CrossRef]
  29. Ainsworth, B.E.; Haskell, W.L.; Whitt, M.C.; Irwin, M.L.; Swartz, A.M.; Strath, S.J.; O’Brien, W.L.; Bassett, J.; Schmitz, K.H.; Emplaincourt, P.O.; et al. Compendium of Physical Activities: An Update of Activity Codes and MET Intensities. Med. Sci. Sports Exerc. 2000, 32, S498–S504. [Google Scholar] [CrossRef]
  30. Hernán, M.A.; Hernández-Diaz, S.; Werler, M.M.; Mitchell, A.A. Causal Knowledge as a Prerequisite for Confounding Evaluation: An Application to Birth Defects Epidemiology. Am. J. Epidemiol. 2002, 155, 176–184. [Google Scholar] [CrossRef]
  31. Kamińska, M.; Ciszewski, T.; Łopacka-Szatan, K.; Miotła, P.; Starosławska, E. Breast Cancer Risk Factors. Prz. Menopauzalny 2015, 14, 196. [Google Scholar] [CrossRef]
  32. Momenimovahed, Z.; Salehiniya, H. Epidemiological Characteristics of and Risk Factors for Breast Cancer in the World. Breast Cancer Targets Ther. 2019, 11, 151. [Google Scholar] [CrossRef]
  33. Risk Factors for Cancer—NCI. Available online: https://www.cancer.gov/about-cancer/causes-prevention/risk (accessed on 7 July 2023).
  34. Correlation and Regression. Available online: https://www.bmj.com/about-bmj/resources-readers/publications/statistics-square-one/11-correlation-and-regression (accessed on 27 July 2023).
  35. Chang, C.P.; Meyers, T.J.; Fu, A.; Zhang, M.Y.; Tashkin, D.P.; Rao, J.Y.; Cozen, W.; Mack, T.M.; Hashibe, M.; Morgenstern, H.; et al. Dietary Glycemic Index, Glycemic Load, and Lung Cancer Risk: A Case-Control Study in Los Angeles County. Cancer Epidemiol. 2020, 69, 101824. [Google Scholar] [CrossRef]
  36. Du, H.; Zhang, T.; Lu, X.; Chen, M.; Li, X.; Li, Z. Glycemic Index, Glycemic Load, and Lung Cancer Risk: A Meta-Analysis of Cohort and Case-Control Studies. PLoS ONE 2022, 17, e0273943. [Google Scholar] [CrossRef] [PubMed]
  37. 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. [Google Scholar] [CrossRef] [PubMed]
  38. Maru, G.B.; Hudlikar, R.R.; Kumar, G.; Gandhi, K.; Mahimkar, M.B. Understanding the Molecular Mechanisms of Cancer Prevention by Dietary Phytochemicals: From Experimental Models to Clinical Trials. World J. Biol. Chem. 2016, 7, 88. [Google Scholar] [CrossRef] [PubMed]
  39. Bertram, J.S. Dietary Carotenoids, Connexins and Cancer: What Is the Connection? Biochem. Soc. Trans. 2004, 32, 985–989. [Google Scholar] [CrossRef]
  40. Toledo, E.; Salas-Salvado, J.; Donat-Vargas, C.; Buil-Cosiales, P.; Estruch, R.; Ros, E.; Corella, D.; Fito, M.; Hu, F.B.; Aros, F.; et al. Mediterranean Diet and Invasive Breast Cancer Risk Among Women at High Cardiovascular Risk in the PREDIMED Trial: A Randomized Clinical Trial. JAMA Intern. Med. 2015, 175, 1752–1760. [Google Scholar] [CrossRef]
  41. Farvid, M.S.; Cho, E.; Eliassen, A.H.; Chen, W.Y.; Willett, W.C. Lifetime Grain Consumption and Breast Cancer Risk. Breast Cancer Res. Treat. 2016, 159, 335–345. [Google Scholar] [CrossRef]
  42. Kang, J.H.; Peng, C.; Rhee, J.J.; Farvid, M.S.; Willett, W.C.; Hu, F.B.; Rosner, B.A.; Tamimi, R.; Heather Eliassen, A. Prospective Study of a Diabetes Risk Reduction Diet and the Risk of Breast Cancer. Am. J. Clin. Nutr. 2020, 112, 1492–1503. [Google Scholar] [CrossRef]
  43. Chen, S.; Chen, Y.; Ma, S.; Zheng, R.; Zhao, P.; Zhang, L.; Liu, Y.; Yu, Q.; Deng, Q.; Zhang, K. Dietary Fibre Intake and Risk of Breast Cancer: A Systematic Review and Meta-Analysis of Epidemiological Studies. Oncotarget 2016, 7, 80980–80989. [Google Scholar] [CrossRef]
  44. Romanos-Nanclares, A.; Collins, L.C.; Hu, F.B.; Willett, W.C.; Rosner, B.A.; Toledo, E.; Eliassen, A.H. Sugar-Sweetened Beverages, Artificially Sweetened Beverages, and Breast Cancer Risk: Results from 2 Prospective US Cohorts. J. Nutr. 2021, 151, 2768–2779. [Google Scholar] [CrossRef]
  45. Schernhammer, E.S.; Bertrand, K.A.; Birmann, B.M.; Sampson, L.; Willett, W.C.; Feskanich, D. Consumption of Artificial Sweetener- and Sugar-Containing Soda and Risk of Lymphoma and Leukemia in Men and Women. Am. J. Clin. Nutr. 2012, 96, 1419–1428. [Google Scholar] [CrossRef]
  46. Romanos-Nanclares, A.; Toledo, E.; Gardeazabal, I.; Jiménez-Moleón, J.J.; Martínez-González, M.A.; Gea, A. Sugar-Sweetened Beverage Consumption and Incidence of Breast Cancer: The Seguimiento Universidad de Navarra (SUN) Project. Eur. J. Nutr. 2019, 58, 2875–2886. [Google Scholar] [CrossRef] [PubMed]
  47. Kaaks, R.; Johnson, T.; Tikk, K.; Sookthai, D.; Tjønneland, A.; Roswall, N.; Overvad, K.; Clavel-Chapelon, F.; Boutron-Ruault, M.C.; Dossus, L.; et al. Insulin-like Growth Factor I and Risk of Breast Cancer by Age and Hormone Receptor Status—A Prospective Study within the EPIC Cohort. Int. J. Cancer 2014, 134, 2683–2690. [Google Scholar] [CrossRef] [PubMed]
  48. Duchaine, C.S.; Dumas, I.; Diorio, C. Consumption of Sweet Foods and Mammographic Breast Density: A Cross-Sectional Study. BMC Public. Health 2014, 14, 554. [Google Scholar] [CrossRef]
  49. Louis, P.; Hold, G.L.; Flint, H.J. The Gut Microbiota, Bacterial Metabolites and Colorectal Cancer. Nat. Rev. Microbiol. 2014, 12, 661–672. [Google Scholar] [CrossRef]
  50. Martínez-González, M.A.; Sayón-Orea, C.; Bullón-Vela, V.; Bes-Rastrollo, M.; Rodríguez-Artalejo, F.; Yusta-Boyo, M.J.; García-Solano, M. Effect of Olive Oil Consumption on Cardiovascular Disease, Cancer, Type 2 Diabetes, and All-Cause Mortality: A Systematic Review and Meta-Analysis. Clin. Nutr. 2022, 41, 2659–2682. [Google Scholar] [CrossRef]
  51. Watling, C.Z.; Schmidt, J.A.; Dunneram, Y.; Tong, T.Y.N.; Kelly, R.K.; Knuppel, A.; Travis, R.C.; Key, T.J.; Perez-Cornago, A. Risk of Cancer in Regular and Low Meat-Eaters, Fish-Eaters, and Vegetarians: A Prospective Analysis of UK Biobank Participants. BMC Med. 2022, 20, 73. [Google Scholar] [CrossRef]
  52. Alexander, D.D.; Morimoto, L.M.; Mink, P.J.; Cushing, C.A. A Review and Meta-Analysis of Red and Processed Meat Consumption and Breast Cancer. Nutr. Res. Rev. 2010, 23, 349–365. [Google Scholar] [CrossRef] [PubMed]
  53. World Cancer Research Fund International. Meat, Fish, Dairy and Cancer Risk—WCRF International. Available online: https://www.wcrf.org/diet-activity-and-cancer/risk-factors/meat-fish-dairy-and-cancer-risk/ (accessed on 17 July 2023).
  54. Koyratty, N.; McCann, S.E.; Millen, A.E.; Nie, J.; Trevisan, M.; Freudenheim, J.L. Sugar-Sweetened Soda Consumption and Total and Breast Cancer Mortality: The Western New York Exposures and Breast Cancer (WEB) Study. Cancer Epidemiol. Biomark. Prev. 2021, 30, 945–952. [Google Scholar] [CrossRef] [PubMed]
  55. Debras, C.; Chazelas, E.; Srour, B.; Druesne-Pecollo, N.; Esseddik, Y.; de Edelenyi, F.S.; Agaësse, C.; De Sa, A.; Lutchia, R.; Gigandet, S.; et al. Artificial Sweeteners and Cancer Risk: Results from the NutriNet-Santé Population-Based Cohort Study. PLoS Med. 2022, 19, e1003950. [Google Scholar] [CrossRef]
  56. Chazelas, E.; Pierre, F.; Druesne-Pecollo, N.; Esseddik, Y.; Szabo De Edelenyi, F.; Agaesse, C.; De Sa, A.; Lutchia, R.; Gigandet, S.; Srour, B.; et al. Nitrites and Nitrates from Food Additives and Natural Sources and Cancer Risk: Results from the NutriNet-Santé Cohort. Int. J. Epidemiol. 2022, 51, 1106–1119. [Google Scholar] [CrossRef]
  57. Adani, G.; Filippini, T.; Wise, L.A.; Halldorsson, T.I.; Blaha, L.; Vinceti, M. Dietary Intake of Acrylamide and Risk of Breast, Endometrial, and Ovarian Cancers: A Systematic Review and Dose⇓response Meta-Analysis. Cancer Epidemiol. Biomark. Prev. 2020, 29, 1095–1106. [Google Scholar] [CrossRef] [PubMed]
  58. Monteiro, C.A.; Cannon, G.; Levy, R.; Moubarac, J.-C.; Jaime, P.; Paula Martins, A.; Canella, D.; Louzada, M.; Parra Also with Camila Ricardo, D.; Calixto, G.; et al. NOVA. The Star Shines Bright. World Nutr. 2016, 7, 28–38. [Google Scholar]
  59. Monteiro, C.A. Nutrition and Health. The Issue Is Not Food, nor Nutrients, so Much as Processing. Public. Health Nutr. 2009, 12, 729–731. [Google Scholar] [CrossRef] [PubMed]
  60. Edefonti, V.; De Vito, R.; Salvatori, A.; Bravi, F.; Patel, L.; Dalmartello, M.; Ferraroni, M. Reproducibility of A Posteriori Dietary Patterns across Time and Studies: A Scoping Review. Adv. Nutr. 2020, 11, 1255–1281. [Google Scholar] [CrossRef] [PubMed]
  61. Kaaks, R.; Slimani, N.; Riboli, E. Pilot Phase Studies on the Accuracy of Dietary Intake Measurements in the EPIC Project: Overall Evaluation of Results. European Prospective Investigation into Cancer and Nutrition. Int. J. Epidemiol. 1997, 26, S26–S36. [Google Scholar] [CrossRef]
  62. Pagliai, G.; Dinu, M.; Madarena, M.P.; Bonaccio, M.; Iacoviello, L.; Sofi, F. Consumption of Ultra-Processed Foods and Health Status: A Systematic Review and Meta-Analysis. Br. J. Nutr. 2021, 125, 308–318. [Google Scholar] [CrossRef]
Figure 1. Association between adherence to a general pro-vegetarian food pattern with cancer risk in participants from the Moli-sani study cohort (n = 22,081). Legend for Figure 1. Cancer-free survival curves were obtained from the multivariable model adjusted for age, sex, energy intake, alcohol intake, residence, educational level, housing tenure, occupational class, smoking status, body mass index, leisure-time physical activity, history of CVD, diabetes, hypertension, hyperlipidemia, aspirin use, oral contraception use, hormone replacement therapy, menopausal status, and family history of cancer. Cancer-free survival curves were generated using the first imputed dataset. The other imputed datasets were similar, and thus omitted.
Figure 1. Association between adherence to a general pro-vegetarian food pattern with cancer risk in participants from the Moli-sani study cohort (n = 22,081). Legend for Figure 1. Cancer-free survival curves were obtained from the multivariable model adjusted for age, sex, energy intake, alcohol intake, residence, educational level, housing tenure, occupational class, smoking status, body mass index, leisure-time physical activity, history of CVD, diabetes, hypertension, hyperlipidemia, aspirin use, oral contraception use, hormone replacement therapy, menopausal status, and family history of cancer. Cancer-free survival curves were generated using the first imputed dataset. The other imputed datasets were similar, and thus omitted.
Nutrients 15 03976 g001
Table 1. Scoring criteria for the pro-vegetarian food patterns in the Moli-sani study cohort, Italy, 2005–2010.
Table 1. Scoring criteria for the pro-vegetarian food patterns in the Moli-sani study cohort, Italy, 2005–2010.
ComponentIncluded FoodsGeneral
Pro-Vegetarian
Healthful
Pro-Vegetarian
Unhealthful
Pro-Vegetarian
Plant food groups (n = 12)
1. VegetablesSpinach, turnip greens, salad, green pepper, pumpkin, tomatoes, carrot, beet, broccoli, brussel sprouts, cauliflower, cabbage, kale, mushrooms, garlic, onions, zucchini, artichoke, fennel, olivesPositivePositiveReverse
2. FruitsCitrus, apple, pear, banana, kiwi, grape, peach, apricot, prune, strawberries, melon, fruit salad, figs, cherries, persimmonPositivePositive Reverse
3. LegumesBeans, chickpeas, lentils, peas, broad beansPositivePositive Reverse
4. Whole grainWhole grain breadPositivePositive Reverse
5. Refined grainsCrispbread/rusks, breakfast cereals, white bread, other bread, rice,
pasta, and other grains
PositiveReversePositive
6. PotatoesPotatoesPositiveReversePositive
7. Nuts and dried fruitWalnut, hazelnut, almond, peanut, dried fruitPositivePositiveReverse
8. Olive oilCommon olive oilPositivePositiveReverse
9. Tea and coffeeTea, caffeinated coffee, decaffeinated coffeeNot scoredPositiveReverse
Fruit juicesFruit juicesNot scoredReversePositive
Sugar-sweetened beveragesCarbonated/soft/isotonic drinks, diluted syrupsNot scoredReversePositive
Sweets and dessertsChocolate, nut spread, candies, cakes, pies, pastries, puddings (non-milk based), biscuits, dry cakes, honey, jam, and sugarNot scoredReversePositive
Animal food groups (n = 5)
Meat and meat productsChicken or turkey, rabbit, pork, beef, lamb, veal, offal, ham, cured meats, salami, mortadella, sausage, hamburgerReverseReverseReverse
Animal fats for cooking or as a spreadButter, another animal fatReverseReverseReverse
EggsEggsReverseReverseReverse
Fish and other seafoodHake, sole, sardines, trout, swordfish, shrimp, prawns, squid, cuttlefish, octopus, clams, stock fish, canned fishReverseReverseReverse
Milk and dairy productsWhole milk, partially-skimmed or skimmed milk, plain yogurt, low-fat yogurt, fruit yogurt, hard cheese, soft cheese, ice creamReverseReverseReverse
Positive scores indicate that a higher consumption of this food group receives higher scores. Reverse scores indicate that a higher consumption of this food group receives lower scores. In the general pro-vegetarian food pattern, consumption of whole grains and refined grains was aggregated as the “grains” food group.
Table 2. Baseline characteristics of study participants across levels of adherence to a general pro-vegetarian food pattern (FP) in the Moli-sani study cohort (n = 22,081) 1.
Table 2. Baseline characteristics of study participants across levels of adherence to a general pro-vegetarian food pattern (FP) in the Moli-sani study cohort (n = 22,081) 1.
General Pro-Vegetarian Food Pattern (Quintiles of)
All Q1Q2Q3Q4Q5
General pro-vegetarian FP (median, IQR)36 (32–40)29 (27–30)33 (32–34)36 (35–37)39 (38–40)43 (42–45)<0.0001
Healthful pro-vegetarian FP (median, IQR)47 (43–51)41 (38–44)44 (42–47)47 (44–50)49 (46–52)53 (50–56)<0.0001
Unhealthful pro-vegetarian FP (median, IQR)47 (43–51)47 (43–52)47.5 (43–52)47 (42–51)47 (42–51)47 (43–51)<0.0001
Mediterranean diet score (median, IQR)4.0 (1.6)3.0 (2.0–4.0)4.0 (3.0–5.0)4.0 (3.0–5.0)5.0 (4.0–6.0)6.0 (5.0–7.0)<0.0001
No. of subjects (%)22,0814720 (21.4)4042 (18.3)4670 (21.2)4037 (18.3)4612 (20.9)-
Age (years; mean, SD)55.2 (11.7)52.1 (11.2)54.6 (11.8)55.6 (11.7)56.5 (11.6)57.2 (11.3)<0.0001
Men48.048.546.848.347.748.50.47
Urban residence67.066.166.666.867.468.20.64
Educational level 0.073
Up to lower secondary52.150.151.951.753.453.6
Upper secondary35.036.635.335.334.233.2
Post-secondary12.913.212.712.912.213.1
Missing data0.10.020.10.10.20.1
Housing <0.0001
Rent8.910.19.38.78.38.0
One dwelling ownership82.283.083.282.182.580.2
More than one dwelling ownership8.96.77.48.99.111.6
Missing data0.20.20.10.30.10.2
Occupational class 0.0032
Professional/managerial20.719.919.620.721.121.9
Skilled non-manual occupations36.537.637.037.136.334.3
Skilled manual occupations18.019.517.817.217.417.9
Partly skilled/Unskilled18.717.019.519.019.319.1
Unemployed/unclassified6.16.06.16.05.96.8
Smoking status 0.10
Non-smokers49.547.550.848.150.551.2
Current23.326.824.023.821.120.4
Former27.125.725.128.028.328.3
Missing data0.10.00.10.10.10.1
Leisure-time PA, MET-h/day (mean, SD) 23.6 (4.0)3.3 (3.9)3.3 (3.7)3.5 (3.9)3.7 (4.1)4.0 (4.5)<0.0001
Body mass index 0.12
Normal weight27.730.228.126.126.827.1
Overweight42.941.443.243.443.243.3
Obese29.328.328.630.429.929.5
Missing data0.10.10.10.10.10.1
Cardiovascular disease 0.0043
Yes5.03.44.15.55.66.6
Missing data1.61.41.31.41.91.8
Diabetes 0.0001
Yes4.73.95.05.45.14.1
Missing data1.31.41.31.11.31.1
Hypertension 0.027
Yes27.821.625.629.131.231.9
Missing data0.70.80.80.80.60.6
Hyperlipidemia <0.0001
Yes7.64.76.58.08.610.3
Missing data0.90.70.80.80.81.3
Aspirin use 0.030
Yes4.62.84.14.75.75.7
Missing data1.91.61.61.81.82.5
Menopausal status 0.91
Yes56.446.454.457.360.964.1
Missing data0.10.00.10.20.10.1
Hormone replacement therapy 0.028
Yes5.64.74.66.65.66.4
Missing data0.010.00.050.00.00.0
Oral contraception use 0.72
Yes28.231.828.228.826.425.6
Missing data0.020.040.050.00.00.0
Family history of cancer40.439.140.240.640.741.50.91
1 Values are numbers (percentages) unless stated otherwise. IQR = interquartile range. Leisure-time physical activity (PA) levels are reported as means adjusted for age, sex, and energy intake. 2 Available for 21,889 participants. The p-values were obtained using generalized linear models both for continuous and categorical dependent variables adjusted for age, sex, and energy intake.
Table 3. Risk of cancer hospitalization in association with pro-vegetarian food patterns (FP) in cancer-free participants from the Moli-sani study cohort (n = 22,081).
Table 3. Risk of cancer hospitalization in association with pro-vegetarian food patterns (FP) in cancer-free participants from the Moli-sani study cohort (n = 22,081).
Quintiles of Dietary Scores
Q1Q2Q3Q4Q5P for Trend1 SD Increment
General pro-vegetarian FP
No. of cancer hospitalizations/no. of subjects470/4720407/4042485/4670448/4037496/4612--
Person-years, n58,22549,52357,28549,38557,472--
Event rates per 10,000 person-years80.782.284.790.786.3--
Model 1 (HR; 95%CI)-1-0.90 (0.79–1.03)0.89 (0.78–1.01)0.92 (0.80–1.04)0.84 (0.74–0.95)0.0200.94 (0.90–0.98)
Model 2 (HR; 95%CI)-1-0.91 (0.80–1.04)0.89 (0.78–1.01)0.93 (0.82–1.06)0.85 (0.75–0.97)0.0450.95 (0.91–0.99)
Healthful pro-vegetarian FP
No. of cancer hospitalizations/no. of subjects382/4017512/5121465/4263501/4690446/3990
Person-years, n48,74362,38052,23258,22250,313
Event rates per 10,000 person-years78.482.189.086.088.6
Model 1 (HR; 95%CI)-1-0.93 (0.82–1.07)1.01 (0.88–1.15)0.94 (0.83–1.08)0.98 (0.85–1.12)0.900.99 (0.95–1.03)
Model 2 (HR; 95%CI)-1-0.92 (0.80–1.05)0.98 (0.85–1.12)0.91 (0.80–1.05)0.93 (0.81–1.07)0.410.97 (0.93–1.01)
Unhealthful pro-vegetarian FP
No. of cancer hospitalizations/no. of subjects474/4427509/4639396/3901490/4752437/4362
Person-years, n55,34757,13847,92058,32853,157
Event rates per 10,000 person-years85.689.182.684.082.2
Model 1 (HR; 95%CI)-1-1.04 (0.91–1.17)0.97 (0.84–1.10)0.97 (0.86–1.10)0.97 (0.85–1.11)0.410.99 (0.95–1.03)
Model 2 (HR; 95%CI)-1-1.05 (0.93–1.19)1.00 (0.87–1.14)1.01 (0.89–1.15)1.04 (0.91–1.18)0.851.01 (0.97–1.06)
Hazard ratios (HR) with 95% confidence interval (95%CI) obtained from multivariable cause-specific Cox proportional hazards regression models. The multivariable-adjusted Model 1 was controlled for age, sex, energy intake, and alcohol intake (g/d; continuous). The multivariable-adjusted Model 2 was controlled for sex, age, energy intake, alcohol intake (g/d; continuous), residence, educational level, housing tenure, occupational class, smoking status, body mass index (categorical), leisure-time physical activity, history of CVD, diabetes, hypertension, hyperlipidemia, aspirin use, oral contraception use, hormone replacement therapy, menopausal status, and family history of cancer.
Table 4. Subgroup and sensitivity analyses for the association of pro-vegetarian food patterns (FP) (1 SD increase) with risk of cancer hospitalization among cancer-free participants from the Moli-sani study cohort (n = 22,081).
Table 4. Subgroup and sensitivity analyses for the association of pro-vegetarian food patterns (FP) (1 SD increase) with risk of cancer hospitalization among cancer-free participants from the Moli-sani study cohort (n = 22,081).
General Pro-Vegetarian FPHealthful Pro-Vegetarian FPUnhealthful Pro-Vegetarian FP
No. of Cases/
No. of Subjects
HR (95%CI)p-Value P for
Interaction
HR (95%CI)p-Value P for
Interaction
HR (95%CI)p-Value P for
Interaction
Subgroup analyses
Sex
Men1335/10,6000.93 (0.88–0.98)0.0110.550.96 (0.90–1.01)0.130.561.04 (0.98–1.10)0.190.32
Women971/11,4810.97 (0.91–1.04)0.450.98 (0.92–1.05)0.621.00 (0.94–1.07)0.92
Age
<65 y1434/17,2390.93 (0.88–0.98)0.0110.860.95 (0.90–1.00)0.0540.731.03 (0.97–1.08)0.320.96
≥65 y872/48420.94 (0.88–1.01)0.100.97 (0.91–1.05)0.491.02 (0.95–1.10)0.57
Body mass index
Normal weight (≤25 kg/m2)540/61190.93 (0.85–1.02)0.110.830.93 (0.86–1.02)0.110.591.10 (1.01–1.20)0.0350.12
Overweight (>25 ≤30 kg/m2)991/94730.95 (0.89–1.02)0.140.98 (0.91–1.04)0.471.01 (0.95–1.08)0.68
Obese (>30 kg/m2)775/64890.96 (0.89–1.04)0.310.99 (0.92–1.07)0.890.96 (0.90–1.04)0.33
Smoking status
Non-smokers919/10,9500.96 (0.90–1.03)0.270.940.99 (0.92–1.06)0.720.801.00 (0.94–1.07)0.890.28
Current smokers588/51400.93 (0.85–1.01)0.0760.97 (0.89–1.05)0.420.98 (0.90–1.07)0.69
Former smokers799/59910.95 (0.89–1.03)0.220.95 (0.88–1.02)0.181.06 (0.99–1.14)0.10
Sensitivity analyses
Excluding CVD2143/20,9430.95 (0.91–0.99)0.0190.97 (0.93–1.02)0.211.00 (0.96–1.05)0.84
Excluding diabetes2122/21,0100.94 (0.90–0.98)0.00460.96 (0.92–1.01)0.111.01 (0.97–1.06)0.52
Excluding early cancer hospitalizations2031/21,8060.96 (0.91–1.00)0.0600.99 (0.94–1.03)0.591.00 (0.96–1.05)0.85
Hazard ratios (HR) with 95% confidence interval (95%CI) obtained from multivariable cause-specific Cox proportional hazards regression models including age, sex, energy intake, alcohol intake, residence, educational level, housing tenure, occupational class, smoking status, body mass index (categorical), leisure-time physical activity, history of CVD, diabetes, hypertension, hyperlipidemia, aspirin use, oral contraception use, hormone replacement therapy, menopausal status, and family history of cancer.
Table 5. Pro-vegetarian food patterns (FP) and risk of hospitalization by cancer site in cancer-free participants from the Moli-sani study cohort (n = 22,081).
Table 5. Pro-vegetarian food patterns (FP) and risk of hospitalization by cancer site in cancer-free participants from the Moli-sani study cohort (n = 22,081).
General Pro-Vegetarian FPHealthful Pro-Vegetarian FPUnhealthful Pro-Vegetarian FP
Cancer site1 SD Increase
HR (95%CI)
Q5 vs. Q1
HR (95%CI)
1 SD Increase
HR (95%CI)
Q5 vs. Q1
HR (95%CI)
1 SD Increase
HR (95%CI)
Q5 vs. Q1
HR (95%CI)
Respiratory tract (n = 183)0.88 (0.76–1.03)0.61 (0.38–0.98)0.87 (0.75–1.02)0.70 (0.41–1.19)1.14 (0.98–1.33)1.68 (1.06–2.68)
Digestive (n = 598)0.88 (0.81–0.96)0.74 (0.57–0.95)0.92 (0.85–1.00)0.76 (0.58–0.99)1.00 (0.92–1.08)0.92 (0.70–1.21)
Genitourinary organs (n = 395)0.97 (0.87–1.08)0.86 (0.63–1.18)1.00 (0.90–1.11)1.07 (0.76–1.51)0.99 (0.90–1.10)0.94 (0.68–1.29)
Breast (n = 285)0.98 (0.87–1.10)0.94 (0.66–1.36)0.98 (0.87–1.10)0.85 (0.58–1.24)1.09 (0.96–1.22)1.40 (0.95–2.05)
Prostate (n = 219)1.08 (0.94–1.24)1.37 (0.88–2.13)1.07 (0.94–1.23)1.37 (0.84–2.23)1.03 (0.90–1.19)1.21 (0.80–1.82)
Lymphatic and Hematopoietic Tissue (n = 178)1.03 (0.88–1.20)1.03 (0.65–1.63)1.06 (0.90–1.23)1.03 (0.631.71)1.04 (0.90–1.22)1.24 (0.77–2.00)
Brain and nervous system (n = 41)0.93 (0.68–1.27)0.92 (0.37–2.29)1.15 (0.85–1.56)1.36 (0.54–3.44)1.04 (0.76–1.43)0.82 (0.31–2.18)
Other (n = 407)0.94 (0.86–1.05)0.80 (0.59–1.08)0.96 (0.87–1.06)1.03 (0.73–1.44)0.94 (0.85–1.04)0.78 (0.56–1.08)
Hazard ratios (HR) with 95% confidence interval (95%CI) obtained from multivariable cause-specific Cox proportional hazards regression models including age, sex, energy intake, alcohol intake, residence, educational level, housing tenure, occupational class, smoking status, body mass index (categorical), leisure-time physical activity, history of CVD, diabetes, hypertension, hyperlipidemia, aspirin use, oral contraception use, hormone replacement therapy, menopausal status, and family history of cancer.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Martínez, C.F.; Di Castelnuovo, A.; Costanzo, S.; Panzera, T.; Esposito, S.; Cerletti, C.; Donati, M.B.; de Gaetano, G.; Iacoviello, L.; Bonaccio, M.; et al. Pro-Vegetarian Food Patterns and Cancer Risk among Italians from the Moli-Sani Study Cohort. Nutrients 2023, 15, 3976. https://doi.org/10.3390/nu15183976

AMA Style

Martínez CF, Di Castelnuovo A, Costanzo S, Panzera T, Esposito S, Cerletti C, Donati MB, de Gaetano G, Iacoviello L, Bonaccio M, et al. Pro-Vegetarian Food Patterns and Cancer Risk among Italians from the Moli-Sani Study Cohort. Nutrients. 2023; 15(18):3976. https://doi.org/10.3390/nu15183976

Chicago/Turabian Style

Martínez, Claudia Francisca, Augusto Di Castelnuovo, Simona Costanzo, Teresa Panzera, Simona Esposito, Chiara Cerletti, Maria Benedetta Donati, Giovanni de Gaetano, Licia Iacoviello, Marialaura Bonaccio, and et al. 2023. "Pro-Vegetarian Food Patterns and Cancer Risk among Italians from the Moli-Sani Study Cohort" Nutrients 15, no. 18: 3976. https://doi.org/10.3390/nu15183976

APA Style

Martínez, C. F., Di Castelnuovo, A., Costanzo, S., Panzera, T., Esposito, S., Cerletti, C., Donati, M. B., de Gaetano, G., Iacoviello, L., Bonaccio, M., & on behalf of the Moli-Sani Study Investigators. (2023). Pro-Vegetarian Food Patterns and Cancer Risk among Italians from the Moli-Sani Study Cohort. Nutrients, 15(18), 3976. https://doi.org/10.3390/nu15183976

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