Learn more: PMC Disclaimer | PMC Copyright Notice
Risk of cancer in regular and low meat-eaters, fish-eaters, and vegetarians: a prospective analysis of UK Biobank participants
Abstract
Background
Following a vegetarian diet has become increasingly popular and some evidence suggests that being vegetarian may be associated with a lower risk of cancer overall. However, for specific cancer sites, the evidence is limited. Our aim was to assess the associations of vegetarian and non-vegetarian diets with risks of all cancer, colorectal cancer, postmenopausal breast cancer, and prostate cancer and to explore the role of potential mediators between these associations.
Methods
We conducted a prospective analysis of 472,377 UK Biobank participants who were free from cancer at recruitment. Participants were categorised into regular meat-eaters (n = 247,571), low meat-eaters (n = 205,385), fish-eaters (n = 10,696), and vegetarians (n = 8685) based on dietary questions completed at recruitment. Multivariable-adjusted Cox regressions were used to estimate hazard ratios (HR) and 95% confidence intervals (CI) for all cancer incidence and separate cancer sites across diet groups.
Results
After an average follow-up of 11.4 years, 54,961 incident cancers were identified, including 5882 colorectal, 7537 postmenopausal breast, and 9501 prostate cancers. Compared with regular meat-eaters, being a low meat-eater, fish-eater, or vegetarian were all associated with a lower risk of all cancer (HR: 0.98, 95% CI: 0.96–1.00; 0.90, 0.84–0.96; 0.86, 0.80–0.93, respectively). Being a low meat-eater was associated with a lower risk of colorectal cancer in comparison to regular meat-eaters (0.91, 0.86–0.96); however, there was heterogeneity in this association by sex (p = 0.007), with an inverse association across diet groups in men, but not in women. Vegetarian postmenopausal women had a lower risk of breast cancer (0.82, 0.68–0.99), which was attenuated and non-significant after adjusting for body mass index (BMI; 0.87, 0.72–1.05); in mediation analyses, BMI was found to possibly mediate the observed association. In men, being a fish-eater or a vegetarian was associated with a lower risk of prostate cancer (0.80, 0.65–0.99 and 0.69, 0.54–0.89, respectively).
Conclusion
The lower risk of colorectal cancer in low meat-eaters is consistent with previous evidence suggesting an adverse impact of meat intake. The lower risk of postmenopausal breast cancer in vegetarian women may be explained by their lower BMI. It is not clear whether the other differences observed for all cancers and for prostate cancer reflect any causal relationships or are due to other factors such as residual confounding or differences in cancer detection.
Supplementary Information
The online version contains supplementary material available at 10.1186/s12916-022-02256-w.
Background
Cancer is a leading cause of death worldwide [1], and in the United Kingdom (UK), 28% of all deaths in 2017 were attributable to cancer [2–4]. Colorectal, breast, and prostate cancer collectively account for 39% of all new cancer diagnoses in the UK [5], and it has been estimated that nearly 40% of cancer cases may be preventable through modifiable factors [6, 7]. Although several dietary factors have been suggested to influence cancer risk, it remains unclear whether dietary patterns are related to the risk of developing cancer [8, 9].
It has been hypothesised that vegetarian diets, which exclude the consumption of all meat and fish, may be associated with a lower cancer risk. In addition to excluding red and processed meat, which are associated with an increased risk of colorectal cancer [8], vegetarians also generally consume higher amounts of plant foods such as fruits, vegetables, and whole grains compared to meat-eaters [10, 11], which might also contribute to lowering the risk of some site-specific cancers [8]. Evidence from two large cohorts which include a large proportion of vegetarians, the European Prospective Investigation into Cancer and Nutrition-Oxford (EPIC-Oxford) and the Adventist Health Study-2 (AHS-2), has suggested that vegetarians may have a lower risk of developing cancer (all types combined) in comparison to meat-eaters [12, 13], but the evidence remains unclear for individual cancer sites [12, 14–17]. Moreover, the risk of cancer in those who do not consume meat but do eat fish (fish-eaters or pescatarians) may differ from that of meat-eaters; some evidence has suggested that fish-eaters may have a lower overall risk of cancer [12], and a lower risk of colorectal cancer [12, 16] than meat-eaters, but no differences have been reported for breast [14] or prostate cancer risk [12, 15]. Despite the substantial number of vegetarians and fish-eaters in these cohorts (8000–25,000 participants), power to detect an association for specific cancer sites may be limited due to relatively small numbers of cancer cases (~ 5000 total cases) in these individual studies [12, 14–17].
Any difference in cancer risk between diet groups may be due to differences in physiological characteristics, including adiposity. In western populations, vegetarians and fish-eaters have been shown to have lower body mass indices in comparison with the body mass index (BMI) of meat-eaters [18–20] which is important for cancer risk because obesity is a known risk factor for several cancer sites [8]. Another hypothesised explanation for the lower risk of cancer observed amongst vegetarians and fish-eaters is the possible differences in hormone levels [21], such as insulin-like growth factor-I (IGF-I) and testosterone, which may be related to their dietary intakes [21–23]. Hormone difference may be important as higher levels of IGF-I have been associated with higher risks of colorectal, breast, and prostate cancer [24] and higher levels of free testosterone have been associated with prostate cancer [25] and postmenopausal breast cancer [26].
To further understand these relationships, we assessed the associations of diet groups with risks of all, colorectal, postmenopausal breast, and prostate cancer in the UK Biobank, which includes approximately 10,000 fish-eaters, 8000 vegetarians, and nearly 55,000 total incident cancer cases. We additionally aimed to assess the roles of BMI, circulating IGF-I, and calculated free testosterone as potential mediators of the observed associations between diet groups and cancer risk.
Methods
Study design and participants
Potential participants were first identified for the UK Biobank study using National Health Service (NHS) records, and 9.2 million eligible individuals, aged 40–70 and living within 25 miles of one of the assessment centres in the UK, were invited to participate in the study. Over 500,000 participants (5.5% response rate) consented to participate between 2006 and 2010 [27] and visited one of 22 assessments centres across England, Wales, and Scotland. A full description of the study protocol can be found on the UK Biobank website [27].
The UK Biobank was approved by the NHS North West Multicentre Research Ethics Committee (21/NW/0157). All participants provided informed consent at recruitment, allowing for follow-up using data-linkage to health records. The study was performed in accordance to the Declaration of Helsinki.
Exclusions
Participants were excluded from this analysis if they withdrew consent over the study period (n = 871), had a prevalent cancer diagnosis at recruitment (excluding non-melanoma skin cancer International Statistical Classification of Disease (ICD-10) code: C44; n = 29,504), their genetic sex was different from their reported sex (n = 321), or they did not contribute any follow-up time (n = 2; Additional File 1 Fig. S1). Participants who responded as ‘do not know’ or ‘prefer not to say’ for all dietary questions regarding meat intake were also excluded from the analyses (n = 282). This left a total of 472,337 participants, of whom 217,937 were males and 254,400 were females. For prostate cancer analyses, women were excluded, and for postmenopausal breast cancer analyses, women who were premenopausal at recruitment and did not reach the age of 55 over the follow-up time (n = 16,222), and men, were excluded.
Diet group classification
Diet groups were categorised using the touchscreen questionnaire completed at recruitment which asked participants about their frequency of consumption of processed meat, beef, lamb or mutton, pork, chicken, turkey or other poultry, and oily and non-oily fish. Participants chose a frequency of intake ranging from “Never” to “Once or more daily”. From these responses, participants were categorised into four diet groups (regular meat-eaters; low meat-eaters; fish-eaters; and vegetarians). Regular meat-eaters were participants who said they consumed processed, red meat (beef, pork, lamb), or poultry > 5 times a week. Low meat-eaters were participants who reported consuming processed, red meat, or poultry ≤ 5 times a week. Fish-eaters were participants who reported that they never consumed red meat, processed meat, or poultry but ate oily and/or non-oily fish. Vegetarians were defined as participants who reported that they never consumed any meat or fish. The vegetarian group also included vegans who reported not consuming any meat, fish, dairy, or eggs (n = 446).
Covariates and biomarkers
The baseline touchscreen questionnaire also asked participants about sociodemographic, reproductive, and lifestyle factors. In addition, all participants had their blood drawn and anthropometric measurements, including height and weight, taken by a trained professional. Further information on covariate data collection and classification can be found in the Additional File 1 Supplementary Methods.
Non-fasting blood samples were provided by 99.7% of participants at recruitment and were shipped to the central processing laboratory at 4 °C prior to serum preparation, aliquoting, and cryopreservation in the central working archive. Biochemistry markers were measured including insulin-like growth factor-I (IGF-I) and testosterone, as well as sex hormone-binding globulin which we used to calculate an estimate of free testosterone [28]. Further description of the UK Biobank biomarker measurements can be found online [29].
Follow-up and outcome ascertainment
Data on cancer diagnosis were ascertained using a combination of records from the NHS Digital (cancer registry) and Public Health England for participants from England and Wales, NHS Central Register for participants from Scotland [30] as well as the Hospital Episodes Statistics (HES) data for English participants and Scottish Morbidity Records (SMR) for Scottish participants (please see details in the Additional File 1 Supplementary Methods). Using the World Health Organization’s ICD-10 codes, participants were classified as having an event if they had an incident diagnosis of cancer recorded as: all cancer (C00-97 excluding non-melanoma skin cancer: C44), colorectal cancer (C18-C20), breast cancer (C50), or prostate cancer (C61), or if no prior incident diagnosis was reported their primary underlying cause of death was the respective cancer. Participants contributed follow-up time from the date of recruitment until the date of the first cancer registration or cancer first recorded on death certificate, date of death, or last day of follow-up available from HES and SMR data (28 February 2021 for England and Scotland). Cancer registry data were available until 31 July 2019 for England and Wales, and 31 October 2015 for Scotland; after this time, only HES and SMR data were used for the follow-up of participants. For Welsh participants, hospital episode data did not extend past the cancer registry censoring date and therefore were not used. For breast cancer, analyses were restricted to postmenopausal breast cancer and women contributed follow-up time beginning when they turned 55 years of age or their date at recruitment if they were categorised as being postmenopausal from questions asked at baseline (see Additional File 1 Supplementary Methods for further details) [31].
Statistical analyses
Baseline characteristics of UK Biobank participants were summarised across diet groups for all participants, and separately for men and women.
Cox proportional hazards regressions were used, with age as the underlying time variable, to estimate hazard ratios (HR) and 95% confidence intervals (CI). Minimally adjusted models were stratified by sex (for all cancer and colorectal cancer analyses only) and age at recruitment (< 45, 45–49, 50–54, 55–59, 60–64, ≥ 65 years) and adjusted for region at recruitment (North-West England, North-Eastern England, Yorkshire & the Humber, West Midlands, East Midlands, South-East England, South-West England, London, Wales, and Scotland).
Multivariable-adjusted Cox regression models for all analyses were further adjusted for height (eight sex-specific categories increasing by 5 cm, and unknown/missing (0.51%)), physical activity (low: 0–9.99, medium: 10–49.99, high: ≥ 50 metabolic equivalent of task-hours /week, and unknown/missing (4.04%)), Townsend deprivation index (quintiles from most deprived to least deprived, and unknown/missing (0.13%)), education (completion of national exam at age 16, completion of national exam at age 17–18, college or university degree, or other/unknown/missing (18.7%)), employment status (employed, retired, not in paid employment, or unknown (1.15%)), smoking status (never, former, light smoker: ≤ 15 cigarettes/day, medium smoker: 16–29 cigarettes/day, heavy smoker: ≥ 30 cigarettes/day, or missing/unknown (0.65%)), alcohol consumption (none drinkers, < 1, 1–9.99, 10–19.99, ≥ 20 g/day, or unknown/missing (0.73%)), ethnicity (White, Mixed race or other, Asian or British Asian, and Black or Black British, or missing/unknown (0.56%)), and diabetes status (no, yes, or unknown (0.53%)).
For colorectal cancer and for all cancer sites, multivariable models were further adjusted for female specific covariates: menopausal hormone therapy (MHT) use (no, former, current, or unknown (0.58%)) and menopausal status at recruitment (premenopausal, postmenopausal, or unknown (9.0%)). Moreover, for colorectal cancer, multivariable models were adjusted for non-steroidal anti-inflammatory drug use (NSAID; no reported use, irregular use, regular use of aspirin/ibuprofen). For prostate cancer, models were additionally adjusted for marital status (not living with a partner, living with a partner) [32]. For postmenopausal breast cancer, models were additionally adjusted for MHT use (same as above), age at menarche (≤ 12 years, 13 years old, ≥ 14 years, or unknown (22.5%)), parity and age at first birth (nulliparous, 1–2 children < 25 years old, 3+ children < 25 years old, 1–2 children 25–29.9 years old, 3+ children 25–29.9 years old, 1–2 children 30+ years old, 3+ children 30+ years old, or missing (0.3%)). Further information on covariate classification can be found in Additional File 1 Supplementary Methods. In all models, the proportional hazards assumption was evaluated using Schoenfeld residuals, and no violations were observed.
We considered BMI as a potential confounder as well as a mediator. When BMI was considered as a potential confounder, BMI measured at recruitment was added to multivariable models (multivariable adjusted + BMI; < 20, 20–22.49, 22.5–24.9, 25.0–27.49, 27.5–29.9, 30–32.49, 32.5–34.9, ≥ 35 kg/m2, or unknown/missing (0.57%)). Models assessing BMI as a mediator are explained below in the mediation analyses section.
To determine if there was heterogeneity in the associations of diet groups with cancer risk, and to assess the influence of confounder adjustments [33, 34], χ2 statistics and p-values for including the diet group in the model were estimated using likelihood ratio tests (LRT) comparing a model without the diet groups variable to the model with the diet groups variable.
Subgroup and sensitivity analyses
For all analyses, we assessed heterogeneity by subgroups of BMI (median: < 27.5 and ≥ 27.5 kg/m2) and smoking status (ever and never) by using a LRT comparing the main model to a model including an interaction term between diet groups and the subgroup variable (BMI and smoking status). For colorectal cancer, we further assessed heterogeneity by sex. For all cancer sites combined, we additionally explored heterogeneity by smoking status, censoring participants at baseline who were diagnosed with lung cancer.
In sensitivity analyses, we excluded cases and participants who had less than 2 years of follow-up and all participants with missing data on covariates. We also examined associations separately in white participants because a large proportion of the vegetarians in this cohort are of South Asian ethnicity (~ 17.5%). Furthermore, we additionally adjusted for fruit and vegetable intake in the multivariable adjusted model (< 3 servings/day, 3–3.99 servings/day, 4–5.99 servings/day, ≥ 6 servings/day, unknown) to control for this component of dietary intake as a proxy for a healthy diet. For prostate cancer analyses, we included in the multivariable adjusted model prostate-specific antigen (PSA) testing (no PSA testing, had PSA test, or unknown) reported at baseline in all men and during follow-up from general practice records in a subsample (n = 99,412 males; records available for participants until 31 May 2016 for England, 31 March 2017 for Scotland, and 31 August 2017 for Wales).
Mediation analyses
If a significant association was observed between a diet group and a cancer outcome in the main analyses, we then further explored potential mediators that have been shown to be associated or possibly associated with diet groups [19, 21] and were previously related to the cancer site of interest (BMI, IGF-I, and free testosterone) [25, 26]. To determine if differences in mediators were observed by diet group, we used multivariable linear regression to compare the selected biomarker measurements (IGF-I and free testosterone [28]) and BMI across dietary groups, adjusting for potential confounders (see Additional File 1 Supplementary Methods). We did not explore mediation if there was no significant difference in cancer risk between each diet group and regular meat-eaters or if the biomarker concentrations were not significantly different between diet groups. We explored mediation via BMI for all cancer, colorectal cancer, and postmenopausal breast cancer risk [8], but not for prostate cancer due to its heterogeneous association with risk by stage and grade [35] and as these data are not available in this cohort. For prostate cancer and postmenopausal breast cancer, we also explored potential mediation via circulating concentrations of IGF-I and calculated free testosterone [25, 26, 28]. We did not explore biomarker mediation for the all cancer–diet group associations as these biomarkers have not been associated with all cancer risk.
To assess for mediation, we used the inverse odds ratio weighting (IORW) method [36, 37]. This method aims to decompose associations between diet group mediated by the potential mediator (natural indirect effect [NIE]) and the estimated association of diet group with cancer risk not mediated by baseline BMI or biomarkers (natural direct effect [NDE]). The term “effect” is used here in concordance with the causal mediation literature but should not be interpreted as implying causality. To determine the proportion of the association between diet groups and cancer outcome mediated by the mediator of interest (e.g. BMI), we took the log of the indirect effect HR and divided it by the log of the total effect HR. Further details of the mediation analyses can be found in the Additional File 1 Supplementary Methods [38, 39].
All analyses were conducted using Stata version 17.0 (Stata Corp LP, College Station, TX). P-values were two-sided with p < 0.05 being considered statistically significant.
Results
Of the participants included in the analysis, 247,571 (52.4%) were classified as regular meat-eaters, 205,385 (43.5%) were low meat-eaters, 10,696 (2.3%) were fish-eaters and 8685 (1.8%) were vegetarians. After an average of 11.4 years of follow-up, 54,961 incident cases of any type of cancer were diagnosed; 5882 participants were diagnosed with colorectal cancer, 7537 women were diagnosed with postmenopausal breast cancer, and 9501 men were diagnosed with prostate cancer.
Table 1 presents participants’ baseline characteristics across diet groups. Vegetarians and fish-eaters had a lower BMI, were younger, more likely to be never smokers, have a university/college degree, and report consuming less alcohol at recruitment compared to regular meat-eaters. Vegetarian men were also less likely to have had a PSA test in comparison to meat-eaters (Table 1). Additional File 1 Table S1 presents the baseline characteristics across diet groups stratified by sex. Both men and women fish-eaters and vegetarians had lower BMIs and were younger at recruitment in comparison to regular meat-eaters.
Table 1
Baseline characteristics across diet groups in UK Biobank
Diet groups | ||||
---|---|---|---|---|
Regular meat-eaters | Low meat-eaters | Fish-eaters | Vegetarians | |
Number of participants | 247,571 | 205,385 | 10,696 | 8685 |
Age at recruitment—years, mean (SD) | 56.0 (8.2) | 56.9 (8.0) | 54.0 (8.0) | 53.0 (7.9) |
Sex—female | 114,849 (46.4%) | 126,165 (61.4%) | 7664 (71.7%) | 5722 (65.9%) |
BMI—kg/m2, mean (SD) | 27.9 (4.9) | 27.0 (4.7) | 25.3 (4.3) | 25.7 (4.7) |
Male height—cm, mean (SD) | 175.7 (6.8) | 175.4 (6.9) | 176.4 (6.9) | 175.5 (7.2) |
Female height—cm, mean (SD) | 162.4 (6.3) | 162.5 (6.3) | 163.5 (6.4) | 162.1 (6.8) |
Physical activity | ||||
Low | 72,811 (29.4%) | 59,752 (29.1%) | 2430 (22.7%) | 2371 (27.3%) |
Moderate | 116,591 (47.1%) | 98,692 (48.1%) | 5710 (53.4%) | 4346 (50.0%) |
High | 48,012 (19.4%) | 38,704 (18.8%) | 2273 (21.3%) | 1690 (19.5%) |
Townsend deprivation index | ||||
Q1—Most affluent | 51,117 (20.6%) | 40,433 (19.7%) | 1777 (16.6%) | 1258 (14.5%) |
Q5—Most deprived | 48,227 (19.5%) | 41,627 (20.3%) | 2385 (22.3%) | 2216 (25.5%) |
Ethnicity | ||||
White | 233,959 (94.5%) | 193,033 (94.0%) | 9922 (92.8%) | 6903 (79.5%) |
Mixed other | 3576 (1.4%) | 3294 (1.6%) | 172 (1.6%) | 152 (1.8%) |
Asian or British Asian | 4114 (1.7%) | 5054 (2.5%) | 369 (3.4%) | 1524 (17.5%) |
Black or Black British | 4218 (1.7%) | 3295 (1.6%) | 167 (1.6%) | 48 (0.6%) |
Education | ||||
National exam at 16 years of age | 41,764 (16.9%) | 34,271 (16.7%) | 1180 (11.0%) | 1099 (12.7%) |
National exam at 17-18 years of age | 13,750 (5.6%) | 10,805 (5.3%) | 578 (5.4%) | 551 (6.3%) |
Degree or college | 146,214 (59.1%) | 119,791 (58.3%) | 8015 (74.9%) | 6109 (70.3%) |
Employment | ||||
In paid employment | 146,078 (59.0%) | 115,579 (56.3%) | 7338 (68.6%) | 6065 (69.8%) |
Retired | 77,483 (31.3%) | 70,640 (34.4%) | 2341 (21.9%) | 1582 (18.2%) |
Not in paid employment | 21,068 (8.5%) | 17,028 (8.3%) | 877 (8.2%) | 921 (10.6%) |
Living with a partner—yes | 187,545 (75.8%) | 141,711 (69.0%) | 6930 (64.8%) | 5771 (66.4%) |
Smoking status | ||||
Never | 132,294 (53.4%) | 114,385 (55.7%) | 6075 (56.8%) | 5561 (64.0%) |
Previous | 85,319 (34.5%) | 69,642 (33.9%) | 3800 (35.5%) | 2480 (28.6%) |
Light smoker < 15 cig/day | 7594 (3.1%) | 6299 (3.1%) | 290 (2.7%) | 210 (2.4%) |
Medium smoker 15-29 cig/day | 10,101 (4.1%) | 6644 (3.2%) | 157 (1.5%) | 139 (1.6%) |
Heavy smoker 30+ cig/day | 10,418 (4.2%) | 7432 (3.6%) | 333 (3.1%) | 252 (2.9%) |
Alcohol intake g/day, mean (SD) | 19.9 (21.3) | 14.9 (16.6) | 13.6 (14.5) | 13.0 (16.1) |
Diabetic—yes | 15,603 (6.3%) | 10,748 (5.2%) | 290 (2.7%) | 465 (5.4%) |
Prostate-specific antigen test reported at baseline or in follow-up—yes, male only | 51,555 (38.8%) | 33,394 (42.2%) | 1125 (37.1%) | 929 (31.4%) |
Female specific covariates | ||||
Age at menarche—years, mean (SD) | 12.6 (2.8) | 12.5 (2.9) | 12.5 (2.9) | 12.4 (3.3) |
Menopausal status | ||||
Premenopausal | 24,939 (21.7%) | 23,360 (18.5%) | 2232 (29.1%) | 1843 (32.2%) |
Postmenopausal | 78,626 (68.5%) | 92,413 (73.2%) | 4717 (61.5%) | 3342 (58.4%) |
Menopausal hormone therapy use | ||||
Never | 70,830 (61.7%) | 76,747 (60.8%) | 5468 (71.3%) | 4386 (76.7%) |
Former | 34,149 (29.7%) | 38,960 (30.9%) | 1590 (20.7%) | 9,66 (16.9%) |
Current | 8993 (7.8%) | 9988 (7.9%) | 584 (7.6%) | 320 (5.6%) |
Parity | ||||
Nulliparous | 17,671 (15.4%) | 25,569 (20.3%) | 2330 (30.4%) | 1736 (30.3%) |
1-2 children | 67,306 (58.6%) | 70,915 (56.2%) | 3910 (51.0%) | 2770 (48.4%) |
3+ children | 29,343 (25.5%) | 29,526 (23.4%) | 1415 (18.5%) | 1203 (21.0%) |
Age at first birth—years, mean (SD) | 25.4 (4.6) | 25.2 (4.6) | 26.4 (5.1) | 26.0 (4.9) |
Percentages include missing values and therefore may not add up to 100%
Values are N (%) unless otherwise indicated
Abbreviations: BMI, body mass index, cig cigarette, Q quintile, SD standard deviation
The minimally adjusted models and sequential adjustments for the associations between diet groups and cancer risks are presented in Additional file 1 Table S2, and Fig. 1 shows the multivariable-adjusted models. In the multivariable-adjusted models (not including BMI), a vegetarian diet was associated with a lower risk of all cancer (HR: 0.86, 95% CI: 0.80–0.93), postmenopausal breast cancer (0.82, 0.68–0.99), and prostate cancer (0.69, 0.54–0.89; Fig. 1A) in comparison to regular meat-eaters. Furthermore, compared to being a regular meat-eater, fish-eaters had a lower risk of all cancers (0.90, 0.84–0.96) and prostate cancer (0.80, 0.65–0.99), and low meat-eaters had a lower risk of colorectal cancer (0.91, 0.86–0.96; Fig. 1A). When including BMI as a potential confounder, associations were slightly attenuated apart from prostate cancer which did not change (Fig. 1B). For postmenopausal breast cancer, after adjustment for BMI, the risk for vegetarians compared to regular meat-eaters was no longer statistically significant (0.87, 0.72–1.05; Fig. 1B).
Multivariable adjusted hazard ratios (95% CI) for diet groups and risk of all cancer, prostate cancer, postmenopausal breast cancer, and colorectal cancer not adjusting for BMI (A) and adjusting for BMI (B). Regular meat-eaters: consume red or processed meat or poultry > 5 times a week. Low meat-eaters: consume red and processed meat or poultry ≤5 times per week. Fish-eaters: do not consume red, processed meat, or poultry but consumed fish. Vegetarians (including vegans): do not consume any meat or fish. All models used age as the underlying time variable and are stratified by sex (for only all cancer and colorectal cancer), age groups, and adjusted for region of recruitment, height, physical activity, Townsend deprivation index, education, employment status, smoking status, alcohol consumption, ethnicity, and diabetes status. For all cancer and colorectal cancer analyses, models were further adjusted for menopausal hormone therapy use and menopausal status and colorectal cancer models are adjusted for non-steroid anti-inflammatory drug use. For prostate cancer models are further adjusted for marital status. For breast cancer models are further adjusted for menopausal hormone therapy use, age at menarche, and age at first birth/parity. Full details for each covariate are provided in the statistical analysis section in the main text. Multivariable + BMI models further adjusts for BMI. P-value for heterogeneity from likelihood ratio tests for model fit comparing a model without diet groups, to a model including diet group. BMI, body mass index; CI, confidence intervals; HR, hazard ratio; N, number of participants
Subgroup analyses
No evidence of heterogeneity was observed across BMI subgroups in the associations between diet groups and risk of all cancer, colorectal, postmenopausal breast, and prostate cancer (Additional File 1 Tables S3-S6). For colorectal cancer, there was evidence of heterogeneity by sex (Phet = 0.007), with male low meat-eaters, fish-eaters, and vegetarians having a lower risk of colorectal cancer (0.89, 0.83–0.95; 0.69, 0.47–1.01; 0.57, 0.36–0.91, respectively) in comparison to regular meat-eaters, whereas no significant association was observed across diet groups for females (Additional File 1 Table S4). For smoking status, some evidence of heterogeneity was observed in the association between diet groups and all cancer risk (phet = 0.056); amongst ever smokers, the low meat-eaters, fish-eaters, and vegetarians had lower risks of all cancer sites than regular meat-eaters (0.97, 0.94–0.99; 0.86, 0.78–0.95; 0.79, 0.70–0.90, respectively), whereas these associations were non-significant for non-smokers (Additional File 1 Table S3). However, when censoring participants who developed lung cancer during follow-up, the test for heterogeneity by smoking status became non-significant between diet groups (phet = 0.223; Additional File 1 Table S3) although the association with diet group for all cancer sites was still only significant amongst ever smokers.
Sensitivity analyses
Associations remained largely the same when analyses were restricted to participants of white European ancestry and when participants with missing data were excluded (Additional File 1 Fig. S2). When the participants who had an event or were censored in the first 2 years of follow-up were excluded, results remained mostly the same except that being a fish-eater was more strongly associated with a lower risk of prostate cancer (HR 0.69, 0.55–0.88) in comparison to regular meat-eaters (Additional File 1 Fig. S3). In analyses, additionally adjusted for intake of fruit and vegetables in the multivariable models, no changes in associations were observed (Additional File 1 Fig. S3). For prostate cancer risk, when PSA testing was added to multivariable models, the associations were not materially changed (Additional File 1 Table S2).
Mediation analyses
Adjusted and relative means of BMI, IGF-I, and free testosterone across diet groups are shown in Additional File 1 Table S7. Explorations of potential mediators for significant diet-cancer associations are shown in Table 2. When we considered the potential of mediation via BMI in the associations of diet groups and risk of all cancer, this was not found to substantially mediate the observed associations (Table 2). For colorectal cancer risk, BMI was not found to mediate the observed lower risk in low meat-eaters compared with regular meat-eaters (Table 2); hormonal biomarkers were not explored because no differences in concentrations were observed between regular and low meat-eaters (Additional File 1 Table S7). For postmenopausal breast cancer risk, BMI was found to be a potential mediator for the observed difference in the risk between vegetarians and regular meat-eaters, with a decomposed HRNIE of 0.83 (95% CI: 0.63–1.08) implying that BMI may explain nearly 93% of the lower risk observed in vegetarian women although this was not statistically significant (Table 2). When IGF-I was explored independently and after adjusting for BMI, a HRNIE of 0.91 (95% CI: 0.73–1.15) was observed. For prostate cancer risk, IGF-I and free testosterone concentrations did not seem to mediate the observed difference in risk between vegetarians and regular meat-eaters, and free testosterone was not found to mediate the difference in risk between fish-eaters and regular meat-eaters (Table 2).
Table 2
Summary of estimated direct effect, indirect effect, and total effect using potential mediators for the association of diet groups in comparison to regular meat-eaters and risk of all cancer, colorectal cancer, postmenopausal breast cancer, and prostate cancer risk
Potential mediators (hazard ratio; 95% CI) | |||
---|---|---|---|
All cancer | Mediation through BMI | Mediation through IGF-Ia | Mediation through free testosteronea |
Low meat-eaters versus regular meat-eaters | Mediation through BMI (n = 450,412) | ||
Total effect | 0.99 (0.96–1.00) | ||
Natural indirect effect | 0.99 (0.98–1.00) | ||
Natural direct effect | 0.99 (0.95–1.00) | ||
Mediation through BMI (n = 256,727) | |||
Fish-eaters versus regular meat-eaters | |||
Total effect | 0.90 (0.83–0.97) | ||
Natural indirect effect | 0.99 (0.90–1.01) | ||
Natural direct effect | 0.90 (0.83–0.98) | ||
Mediation through BMI (n = 254,709) | |||
Vegetarians versus regular meat-eaters | |||
Total effect | 0.86 (0.78–0.96) | ||
Natural indirect effect | 0.94 (0.81–1.08) | ||
Natural direct effect | 0.92 (0.77–1.09) | ||
Colorectal cancer | Mediation through BMI (n = 450,412) | ||
Low meat-eaters versus regular meat-eaters | |||
Total effect | 0.91 (0.85–0.97) | ||
Natural indirect effect | 1.00 (0.98–1.02) | ||
Natural direct effect | 0.91 (0.86–0.97) | ||
Postmenopausal breast cancer | Mediation through BMI (n = 111,574) | Mediation through IGF-I (n = 103,853) | Mediation through free testosterone (n = 93,662) |
Vegetarians versus regular meat-eaters | |||
Total effect | 0.82 (0.68–0.99) | 0.86 (0.71–1.05) | 0.86 (0.71–1.05) |
Natural indirect effect | 0.83 (0.63–1.08) | 0.91 (0.73–1.15) | 1.06 (0.76–1.38) |
Natural direct effect | 0.99 (0.79–1.23) | 0.94 (0.70–1.21) | 0.81 (0.62–1.05) |
Prostate cancerb | Mediation through IGF-I (n = 126,538) | Mediation through free testosterone (n = 116,087) | |
Vegetarians versus regular meat-eaters | |||
Total effect | 0.71 (0.56–0.92) | 0.71 (0.56–0.92) | |
Natural indirect effect | 1.10 (0.77–1.56) | 0.99 (0.67–1.48) | |
Natural direct effect | 0.64 (0.50–1.01) | 0.71 (0.51- 1.01) | |
Fish-eaters versus regular meat-eaters | Mediation through free testosterone (n = 116,186) | ||
Total effect | 0.80 (0.65–0.99) | ||
Natural indirect effect | 0.95 (0.70–1.29) | ||
Natural direct effect | 0.86 (0.56–1.32) |
All models used age as the underlying time variable and are stratified by sex (for only all cancer and colorectal cancer) and age groups at recruitment, and adjusted for region of recruitment, height, physical activity, Townsend deprivation index, education, employment status, smoking status, alcohol consumption, ethnicity, diabetes status, and body mass index (except when it was considered a potential mediator). For all cancer and colorectal cancer, models are further adjusted for menopausal hormone therapy use and menopausal status. Colorectal cancer models are adjusted for non-steroid anti-inflammatory drug use. For prostate cancer, models are further adjusted for marital status. For breast cancer, models are further adjusted for menopausal hormone therapy use, age at menarche, and age at first birth/ parity. Full details for each covariate are provided in the statistical analysis section in the main text.
Mediation analyses restricted to significant associations between diet-cancer in the main analyses (Fig. (Fig.1)1) and if there was a significant difference in biomarker concentrations between diet group (Additional File 1 Table S7).
Natural indirect effect represents the estimated association of diet group and cancer outcome through the potential mediator
Natural direct effect represents the estimated association of diet group and cancer outcome not through the potential mediator
Models exclude participants with missing values for mediator(s)
aModels are adjusted for BMI
bBMI not assessed as a mediator with total prostate cancer risk. Association of IGF-I and free testosterone presented as both hormones have been associated with prostate cancer risk. IGF-I concentrations not assessed for fish-eaters as no difference in concentrations in comparison to regular meat-eaters was observed.
Abbreviations: BMI body mass index, CI confidence intervals, IGF-I insulin like growth factor-I
Discussion
In this large British cohort, being a low meat-eater, fish-eater, or vegetarian was associated with a lower risk of all cancer sites when compared to regular meat-eaters. We also found a lower risk of colorectal cancer amongst low meat-eaters, a lower risk of postmenopausal breast cancer risk in vegetarian women, and a lower risk of prostate cancer amongst vegetarian men. The lower risk of postmenopausal breast cancer in vegetarians may be largely a result of vegetarians having a lower BMI than regular meat-eaters, with possibly some further impact due to vegetarian women in this population having slightly lower circulating IGF-I concentrations.
All cancer
In this study, vegetarians, fish-eaters, and low meat-eaters all had a lower risk of developing all cancer in comparison to regular meat-eaters. It is important to consider that although some cancers may have similar aetiologies, some cancer sites may not be associated with dietary or nutritional factors and that using all cancer incidence as an outcome may crudely capture other lifestyle factors, outside of diet, that may be associated with cancer risk and may confound the associations observed; therefore, these results should be interpreted carefully. In the two largest previous prospective studies following vegetarians, EPIC-Oxford and AHS-2 found that being a vegetarian was associated with a 10% and 8% lower risk of all cancer than being a meat-eater, respectively, after adjusting for lifestyle risk factors and BMI [12, 13]. Fish-eaters in EPIC-Oxford had a lower risk of developing all cancer [12], but no association with risk for all cancer was observed for fish-eaters in comparison to meat-eaters in AHS-2 [13]. In the current analysis, we observed some evidence of heterogeneity by smoking status, and when we removed lung cancer from all cancer cases, significant associations were only observed across diet groups within the ever smoker subgroup. Therefore, the differences observed between diet groups for all cancer outcomes combined may not be due to diet and might be due to residual confounding by differences in other lifestyle factors, such as smoking.
Colorectal cancer
The risk of colorectal cancer was lower in low meat-eaters in comparison to regular meat-eaters whereas there was no significant difference for fish-eaters and vegetarians, potentially due to lack of power as the point estimates suggested lower risks in both these non-meat-eating diet groups. In both EPIC-Oxford and AHS-2, being a fish-eater was associated with a lower risk of colorectal cancer in comparison with meat-eaters, whereas no association was observed for being vegetarian and risk of colorectal cancer compared to regular meat-eaters [12, 16]. We also observed heterogeneity by sex, in that significant inverse associations were observed with risk across diet groups in men, when compared to regular meat-eaters, but not for women. This may in part be due to dietary differences between sexes; however, the number of colorectal cancer cases in some diet groups was too small to draw a clear conclusion. The intake of processed meat has been evaluated by the World Health Organization and World Cancer Research Fund to be a definite cause of colorectal cancer [40] and red meat as a probable cause of colorectal cancer [40, 41]. This is likely to at least in part explain the lower risk of colorectal cancer in low meat-eaters, and mechanisms suggested include chemicals in meat such as nitrosamines [40, 42]. Overweight and obesity also increase the risk for colorectal cancer [43, 44], but in mediation analyses, BMI did not appear to mediate the difference observed between low meat-eaters and regular meat-eaters.
Postmenopausal breast cancer
A borderline significantly lower risk for postmenopausal breast cancer was observed for vegetarian women, which appeared to be largely due to their lower BMI as evidenced in mediation analyses and the attenuation of estimates when analyses were adjusted for BMI. We also observed a small potential effect for mediation for lower risk of postmenopausal breast cancer for vegetarians through lower IGF-I concentrations, perhaps influenced by the inclusion of vegans in this group [23]. To date, studies have reported a non-significantly lower risk of breast cancer for women following a vegetarian or pescatarian diet with or without adjustment for BMI [12, 14, 17, 45], which may be due to lack of power to detect modest associations in individual studies. Breast cancer is a heterogeneous disease, with differing risk factors by menopausal status and hormone receptor status [46]. BMI is robustly associated with higher postmenopausal breast cancer risk, probably due to higher circulating oestrogen derived from aromatisation of androgens in the adipose tissue [46]. As such, being vegetarian would be expected to confer a lower risk of postmenopausal breast cancer in comparison to meat-eaters because vegetarians generally have a lower BMI, but whether BMI is a confounder or a mediator for this association is not clear; if vegetarians have a lower BMI because of their diet then BMI may be considered a mediator, but if vegetarians have a lower BMI that is not due to their dietary intake but rather due to other non-dietary lifestyle factors (e.g. physical activity), then BMI would be considered to be a confounder.
Previous research has also suggested that vegetarian women are less likely to use MHT or to attend breast cancer screening [47]. In this analysis, we adjusted for MHT use at baseline, but residual confounding due to differences in use of MHT during follow-up is still possible. Data on breast cancer screening during follow-up were not available in this cohort; therefore, full adjustment for screening attendance was not possible and differences between diet groups in screening may have influenced our findings.
Prostate cancer
The risk of prostate cancer was lower in men who were vegetarians or fish-eaters in comparison to regular meat-eaters, but no difference in risk was observed for low meat-eaters. Previous analyses in the EPIC-Oxford cohort found a non-significantly lower risk of prostate cancer for British vegetarians and fish-eaters in comparison to meat-eaters [12]. In the AHS-2 study, no difference was found for vegetarians or fish-eaters, whereas being vegan was associated with a 35% lower risk of prostate cancer (based on 1079 cases in the cohort of which only 59 were in vegans) [15]. To date, no established dietary risk factor has been found in relation to prostate cancer risk, although there is some evidence which suggests that higher intake of dairy products, and possibly milk specifically, may increase the risk of prostate cancer [48]. This association has been proposed to be possibly mediated through IGF-I [22, 49], a hormone shown to be positively associated with both milk intake and prostate cancer risk [25, 50]. In this cohort, slightly lower IGF-I concentrations have been observed in vegetarians compared to regular meat-eaters [21], and IGF-I has also been associated with prostate cancer risk [25]; however, the difference in IGF-I concentrations between these diet groups is small and may not confer a substantial difference in prostate cancer risk. As might be expected, in mediation analyses, the estimates were imprecise and there was no evidence that the difference in IGF-I concentrations between diet groups mediates the observed associations with cancer risk.
In this cohort, vegetarian men were less likely than meat-eaters to have had a PSA screening test at recruitment; therefore, vegetarians may have a lower risk of having prostate cancer diagnosed following a PSA test. Similarly, two other cohorts have also reported that vegetarian men were less likely to have had a PSA test [47, 51]. When data at recruitment and available general practice records during follow-up were assessed for PSA testing in the UK Biobank, there was only a small difference with 40% of regular meat-eaters and 37% of vegetarians reporting having had a PSA test (although general practice records were only available for half of the participants) after adjusting for age differences. Adding PSA screening in the multivariable-adjusted model did not attenuate the estimates, suggesting the differences in PSA screening in vegetarians or fish-eaters in comparison to regular meat-eaters does not explain the observed associations, but other differences in attendance for medical examinations could possibly also contribute. Due to unavailable data in UK Biobank, we were also unable to assess associations by tumour subtypes which may be aetiologically different [35]. Considering the substantial difference in risk we observed for vegetarian men, differences in detection and residual confounding, as well as chance, may contribute to this observed difference.
Interpretation of results: role of confounding and mediation
The role of residual and unmeasured confounding must be considered when interpreting the findings from this study. Vegetarians and fish-eaters differ from meat-eaters in many non-dietary lifestyle factors such as lower smoking and alcohol consumption, and higher physical activity [52]. Although relevant potential confounders were added to the multivariable models to adjust for these differences, imperfect measurements and/or changes in these confounders over time may result in incomplete adjustment for these variables. For example, the evidence of heterogeneity by smoking status when looking at all cancer as an outcome suggested that residual confounding by smoking may be present.
Differences in BMI between diet groups have also been suggested to explain the lower cancer incidence observed amongst vegetarians [18]; however, when BMI was considered as a potential confounder and mediator, the difference between BMI by diet groups only slightly attenuated the estimates, with the exception of postmenopausal breast cancer. Whether differences in BMI by diet group is due solely to their diet or other lifestyle factors remains unclear, making it difficult to tease out whether BMI mediates or confounds the associations between diet group and cancer risk.
Strengths and limitations
Strengths of this study include the prospective nature and moderately long follow-up time of participants. Data-linkage to health records was used to determine cancer diagnoses which minimises misclassification and loss to follow-up of participants. The UK Biobank study also gathered data on an array of potential confounders and biochemical biomarkers amongst participants; thus, we were able to adjust the models for potential confounding as well as conduct mediation analysis exploring potential mediators between diet groups and cancer risk. When analyses excluded the first 2 years of follow-up, the results remained largely the same, reducing the chance that these associations are due to reverse causality.
There are some limitations to consider in these analyses. Although there were many cancer cases accrued during the follow-up period, these analyses may still be underpowered to detect moderate associations due to the relatively low numbers of cancer cases amongst vegetarians and fish-eaters in this cohort. We also used hospital admissions data to follow-up participants after 2015 in Scotland and 2019 in England because cancer registry data were not available in UK Biobank after this date, which may result in some missing cancer cases and relatively later dates of diagnosis. We also were unable to adjust for total energy intake as this could not be calculated due to the limited number of dietary questions asked at recruitment. As detailed above, the results may be influenced by unmeasured and residual confounding, as well as chance with numerous comparisons, and causality cannot be confirmed. Misclassification of diet may also be possible, as participants may have underreported their intake or changed their diet over the follow-up period, possibly resulting in attenuation of the risk estimates. Vegetarian diets are characterised by not consuming meat, however, this does not necessarily mean that all vegetarians follow a healthy diet, which may influence their risk of cancer and these results. The mediation analyses only explored three potential mediators, and other possible mediating factors such as other biomarkers relevant for cancer, such as oestradiol, were not available. Moreover, baseline BMI and biomarkers were used to assess mediation, and therefore, these measures may not represent BMI during the follow-up and long-term biomarker concentrations, although correlations with repeat measures of BMI and biomarkers showed high agreement [53]. As well, we used the IORW method to explore mediation with bootstrapping CIs, which may make estimates less statistically efficient in comparison to parametric methods, but the IORW has the advantage that it can be applied in survival analysis and provides estimates of the proportion mediated for the mediators of interest. The UK Biobank has been shown to have a healthier risk profile than the UK population [54] and only included British participants most of whom are of white European ancestry (94%); this may limit generalizability to other populations. However, the risks estimated may still be valid to estimate relative differences for risk-factor disease associations [55].
Conclusions
In conclusion, this study found that being a low meat-eater, fish-eater, or vegetarian was associated with a lower risk of all cancer, which may be a result of dietary factors and/or non-dietary differences in lifestyle such as smoking. Low meat-eaters had a lower risk of colorectal cancer, vegetarian women had a lower risk of postmenopausal breast cancer, and men who were vegetarians or fish-eaters had a lower risk of prostate cancer. BMI was found to potentially mediate or confound the association between vegetarian diets and postmenopausal breast cancer. It is not clear if the other associations are causal or a result of differences in detection between diet groups or unmeasured and residual confounding. Future research assessing cancer risk in cohorts with large number of vegetarians is needed to provide more precise estimates of the associations and to explore other possible mechanisms or explanations for the observed differences.
Acknowledgements
This work has been conducted using the UK Biobank Resource under Application Number 24494, and we wish to express our gratitude to the participants and those involved in building the resource. UK Biobank is an open access resource.
Abbreviations
AHS-2 | Adventist Health Study-2 |
BMI | Body mass index |
CI | Confidence intervals |
EPIC-Oxford | European Prospective Investigation into Cancer and Nutrition – Oxford |
HES | Hospital episode statistics |
HR | Hazard ratio |
ICD | International Statistical Classification of Diseases |
IGF-I | Insulin-like growth factor-I |
IORW | Inverse odds ratio weighting |
LRT | Likelihood ratio test |
MHT | Menopausal hormone therapy |
NDE | Natural direct effect |
NIE | Natural indirect effect |
NHS | National Health Service |
NSAID | Non-steroidal anti-inflammatory drugs |
PSA | Prostate-specific antigen |
SMR | Scottish Morbidity Records |
UK | United Kingdom |
Authors’ contributions
CZW, TJK, and AP-C conceived and designed the research question; AP-C managed the project and was responsible for acquiring the data; CZW analysed the data, prepared the figures and tables, and wrote the first draft of the manuscript; CZW prepared the data for analysis; and JAS, YD, TYNT, RKK, AK, RCT, TJK, AP-C provided input on data analysis and interpretation of results. All authors revised the manuscript critically for important intellectual content and read and approved the final manuscript.
Funding
This work is supported by the Cancer Research UK grant (C8221/A29017) and by the World Cancer Research Fund (WCRF UK), as part of the Word Cancer Research Fund International grant programme (2019/1953). CZW is supported by the Nuffield Department of Population Health Doctor of Philosophy student scholarship. YD is supported by the World Cancer Research Fund UK grant (2019/1953). TYNT is supported by the Nuffield Department of Population Health Intermediate Fellowship and the UK Medical Research Council (MR/M012190/1). RKK is supported by the Clarendon Scholarship from the University of Oxford. AK is supported by the Wellcome Trust, Our Planet Our Health (Livestock, Environment and People – LEAP; 205,212/Z/16/Z). APC is supported by a Cancer Research UK Population Research Fellowship (/ C60192). A28516
Availability of data and materials
Bona fide researchers can apply to use the UK Biobank dataset by registering and applying at http://ukbiobank.ac.uk/register-apply/.
Declarations
The UK Biobank protocol was approved by the NHS North West Multicentre Research Ethics Committee (21/NW/0157). All participants provided informed consent at recruitment, allowing for follow-up using data-linkage to health records.
Not applicable.
The authors declare no competing interests.
Footnotes
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.