# Degree of food processing and breast cancer risk: a prospective study in 9 European countries

Manon Cairat, Jessica Blanco Lopez, Aline Al Nahas, Carine Biessy, Vivian Viallon, Corinne Casagrande, Genevieve Nicolas, Marc J. Gunter, Inge Huybrechts, Thibault Violet, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Manon Cairat, Jessica Blanco Lopez, Aline Al Nahas, Carine Biessy, Vivian Viallon, et al. Degree of food processing and breast cancer risk: a prospective study in 9 European countries. Food Production, Processing and Nutrition, 2024, 6 (1), pp.89. 10.1186/s43014-024-00264-2. hal-04831673

HAL Id: hal-04831673

https://hal.science/hal-04831673v1

Submitted on 21 Jan 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# Degree of food processing and breast cancer risk: a prospective study in 9 European countries

Manon Cairat¹,²,³,⁴*†, Sahar Yammine³†, Thibault Fiolet², Agnès Fournier², Marie-Christine Boutron-Ruault², Nasser Laouali², Francesca Romana Mancini², Gianluca Severi², Fernanda Morales Berstein⁴,⁵, Fernanda Rauber⁶, Renata Bertazzi Levy⁶, Guri Skeie⁷, Kristin Benjaminsen Borch⁷, Anne Tjønneland⁸,⁹, Lene Mellemkjær⁸, Yan Borné¹⁰, Ann H. Rosendahl¹¹, Giovanna Masala¹², Maria Teresa Giraudo¹³, Maria Santucci de Magistris¹⁴, Verena Katzke¹⁵, Rashmita Bajracharya¹⁵, Carmen Santiuste¹⁶,¹⁷, Pilar Amiano¹⁷,¹⁸,¹⁹, Stina Bodén²⁰, Carlota Castro-Espin²¹,²², Maria-Jose Sánchez¹⁷,²³,²⁴,²⁵, Mathilde Touvier³,²⁶, Mélanie Deschasaux-Tanguy³,²⁶, Bernard Srour³,²⁶, Matthias B. Schulze²⁷,²⁸, Marcela Guevara¹⁷,²⁹,³⁰, Nathalie Kliemann³¹, Jessica Blanco Lopez¹, Aline Al Nahas¹, Kiara Chang³², Eszter P. Vamos³², Christopher Millett³², Elio Riboli³³, Alicia K. Heath³³, Carine Biessy¹, Vivian Viallon¹, Corinne Casagrande¹, Genevieve Nicolas¹, Marc J. Gunter¹,³³ and Inge Huybrechts¹

## Abstract

Recent epidemiological studies have suggested a positive association between ultra-processed food consumption and breast cancer risk, although some studies also reported no association. Furthermore, the evidence regarding the associations between intake of food with lower degrees of processing and breast cancer risk is limited.

Thus, we investigated the associations between dietary intake by degree of food processing and breast cancer risk, overall and by breast cancer subtypes in the European Prospective Investigation into Cancer and Nutrition (EPIC) study.

Dietary intake of EPIC participants was assessed via questionnaires at baseline. More than 11,000 food ingredients were classified into four groups of food processing levels using the NOVA classification system: unprocessed/minimally processed (NOVA 1), culinary ingredients (NOVA 2), processed (NOVA 3) and ultra-processed (NOVA 4). Cox proportional hazards models were used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) of breast cancer per standard deviation increase in daily consumption (grams) of foods from each NOVA group.

The current analysis included 14,933 breast cancer cases, diagnosed among the 318,686 EPIC female participants, (median follow-up of 14.9 years). No associations were found between breast cancer risk and the level of dietary intake from NOVA 1 [HR per 1 SD =0.99 (95% CI 0.97 – 1.01)], NOVA 2 [HR per 1 SD =1.01 (95% CI 0.98 – 1.03)] and NOVA 4 [HR per 1 SD =1.01 (95% CI 0.99 – 1.03)] foods. However, a positive association was found between NOVA 3 and breast cancer risk [HR per 1 SD =1.05 (95% CI 1.03 – 1.07)] which became non-significant after adjustment for alcohol intake [HR per 1 SD =1.01 (95% CI 0.98 – 1.05)] or when beer and wine were excluded from this group [HR per 1 SD =0.99 (95% CI 0.97 – 1.01)]. The associations did not differ by breast cancer subtype, menopausal status or body mass index.

<footnote>†Manon Cairat and Sahar Yammine equally contributed.</footnote>

<footnote>*Correspondence:
Manon Cairat
manon.cairat@gustaveroussy.fr
Full list of author information is available at the end of the article</footnote>
Findings from this large-scale prospective study suggest that the positive association between processed food intake and breast cancer risk was likely driven by alcoholic beverage consumption.

**Keywords** Epidemiology, Prospective study, Breast cancer, NOVA classification, Food processing

## Graphical Abstract

### Degree of food processing and breast cancer risk: a prospective study in 9 European countries

Aim: To investigate the associations between dietary intakes based on the degree of food processing and risk of breast cancer, overall and by breast cancer subtypes in the EPIC cohort

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

**Methods:**

The European Prospective Investigation into Cancer and Nutrition (EPIC) study includes 521,323 participants recruited between 1992 and 2000 across 10 European countries.

Dietary intake of EPIC participants was assessed via questionnaires at baseline. food ingredients were classified into the four groups of NOVA.

The current analysis included 14,933 breast cancer cases.

Cox proportional hazards models were used to estimate hazard ratios of breast cancer per standard deviation increase in daily consumption of foods from each NOVA group (g/day).

<table><caption>Table 1: Associations between NOVA group 3 and breast cancer risk</caption><thead><tr><th>Breast cancer subtypes</th><th>N Cases</th><th>NOVA 3 HR per 1 SD (95%CI)</th><th>P homogeneity</th></tr></thead><tbody><tr><td>Overall</td><td>14,933</td><td>1.05 (1.03-1.07)</td><td></td></tr><tr><td>In situ</td><td>1603</td><td>1.01 (0.96-1.07)</td><td>0.18</td></tr><tr><td>Invasive</td><td>13,320</td><td>1.05 (1.03-1.07)</td><td></td></tr><tr><td>Invasive ER+</td><td>7789</td><td>1.06 (1.04-1.08)</td><td>0.86</td></tr><tr><td>Invasive ER-</td><td>1736</td><td>1.06 (1.01-1.11)</td><td></td></tr><tr><td>Invasive PR+</td><td>5268</td><td>1.05 (1.03-1.08)</td><td>0.51</td></tr><tr><td>Invasive PR-</td><td>2726</td><td>1.07 (1.03-1.11)</td><td></td></tr><tr><td>Invasive HER2+</td><td>901</td><td>1.04 (0.98-1.12)</td><td>0.57</td></tr><tr><td>Invasive HER2-</td><td>3676</td><td>1.07 (1.03-1.10)</td><td></td></tr></tbody></table>

Models were stratified by age and center and adjusted for education, height, physical activity, age at menarche, oral contraceptive use, age at first full-term pregnancy, parity, breastfeeding, menopausal status and menopausal hormone therapy use. Each NOVA group was mutually adjusted for the other NOVA groups.

**Results:**

No associations were found between breast cancer risk and the level of dietary intake from NOVA 1, NOVA 2 and NOVA 4. However, a positive association was found between NOVA 3 and breast cancer risk which became non-significant after adjustment for alcohol intake [HR per 1 SD =1.01 (95% CI 0.98 - 1.05)].

Conclusion: In this large-scale prospective cohort, a positive association between processed food intake and breast cancer risk was observed, but the association was likely driven by alcohol intake.

## Introduction

Diets have transitioned from fresh, unprocessed, and minimally processed foods toward a rise in the consumption of ultra-processed foods. These now contribute roughly 25 to 60% of the total daily energy intake of individuals across countries (Adams & White, 2015; Juul et al., 2022; Latasa et al., 2018; Levy et al., 2022; Madruga et al., 2022; Marrón-Ponce et al., 2019; Moubarac et al., 2014a; Wang et al., 2021). In the past decades, several food processing frameworks have been developed, such as the NOVA classification, the Food Standards Australia New Zealand, the International Food and Information Council, and the International Food Policy and Research Institute (Bleiweiss-Sande et al., 2019; Crino et al., 2017; Moubarac et al., 2014b). Although most of these similarly classify basic foods as processed or unprocessed, the NOVA system is the most widely applied to scientific studies and may be more useful for monitoring changes in the food supply and evaluating associations with health outcomes (Crino et al., 2017; Monteiro

et al., 2019). The NOVA classification categorizes foods into four groups based on the extent and purpose of food processing they undergo: unprocessed/minimally processed foods (NOVA 1), processed culinary ingredients (NOVA 2), processed foods (NOVA 3) and ultra-processed foods (NOVA 4) (Monteiro et al., 2019). Foods classified as NOVA 4, which are foods that undergo multiple physical, biological, and/or chemical processes, have been suggested to have detrimental health effects due to their poorer nutritional qualities on average (e.g. they are often energy dense and/or rich in saturated and trans-fatty acids) or the presence of a wide range of additives and contaminants formed during processing (Lane et al., 2021; Pagliai et al., 2021; Srour et al., 2022). Epidemiological studies investigating the association between the consumption of foods classified by the NOVA scale and the risk of breast cancer are sparse. Breast cancer risk was inversely associated with NOVA 1 (Fiolet et al., 2018; Jacobs et al., 2022), positively associated with NOVA 3 (Kliemann et al., 2023), and either not associated (Jacobs
et al., 2022; Romaguera et al., 2021) or positively (Chang et al., 2023; Fiolet et al., 2018; Romieu et al., 2022) associated with NOVA 4. A recent analysis in the European Prospective Investigation into Cancer and Nutrition (EPIC), performed by our group, suggested that higher consumption of NOVA 1 was associated with lower breast cancer risk while higher consumption of NOVA 3 was associated with higher breast cancer risk (Kliemann et al., 2023). In that study, we did not stratify our analyses by breast cancer subtype although such analyses could shed light on the potential mechanisms underlying the associations observed between different degrees of food processing and breast cancer, since the etiologies behind breast cancer subtypes are different. Furthermore, few previous studies were able to stratify their analyses by breast cancer subtypes or by alcohol intake and body mass index (BMI). Therefore, the aim of this study was to investigate the associations between diet according to the degree of food processing and breast cancer risk, overall and by breast cancer subtype, menopausal status, alcohol intake and BMI, within the EPIC cohort.

# METHODS

## EPIC Cohort

Between 1992 and 2000, a total of 521,323 EPIC participants were recruited from 23 centers across 10 European countries (Riboli et al., 2002). At recruitment, socio-demographic, dietary, lifestyle, anthropometric and medical data were collected for all participants by administration of validated country-specific questionnaires. Ethical approval for the study was obtained from the relevant ethical review boards of the participating centers of EPIC as well as from the ethics committee of the International Agency for Research on Cancer (IARC).

## Study population and follow-up

The selection of the study population is shown in Supplementary Figure S1. We used data from all participating countries apart from Greece due to a lack of data access. Participants were further excluded if they (i) had any cancer diagnosis before recruitment, (ii) had no follow-up, (iii) had no lifestyle or dietary information, (iv) had an energy intake-to-requirement ratio within the extreme ranking (top and bottom 1%, which are implausible dietary exposure values), or (v) were men. Women were followed from study inclusion until the date of their latest known contact, cancer diagnosis, death, emigration, or the end of the follow-up period (between 2008 and 2014 depending on the center), whichever occurred first. The analytical sample included 318,686 women who were free of cancer at recruitment.

## Identification of incident breast cancer cases

In Italy, Spain, the United Kingdom, the Netherlands, Denmark, Norway, and Sweden, population-based cancer registries were used to identify breast cancer cases. In France and Germany, a combination of methods was used, including health insurance records, contacts with cancer and pathology registries, and active follow-up of participants and their next of kin. Almost all centers (except Malmö, Granada, and Murcia) had information on tumor characteristics, including invasiveness status (in situ/invasive/unknown), estrogen receptor (ER) status (ER-positive/ER-negative/unknown), progesterone receptor (PR) status (PR-positive/PR-negative/unknown), and human epidermal growth factor receptor 2 (HER2) status (HER2-positive/HER2-negative/unknown). The diagnosis of breast cancer cases was based on the 2nd or the 3rd revision (depending on the year of diagnosis) of the International Classification of Diseases for Oncology (ICD-O-2 or ICD-O-3) (International Statistical Classification of Diseases and Related Health Problems 10th Revision, n.d.). In the present work, the first diagnosis of breast cancer was identified as primary incident breast tumors. Vital status was collected from regional or national mortality registries.

## Dietary data and NOVA classification

Country-specific dietary questionnaires were used in EPIC and validated at the center level (Huybrechts et al., 2022). Semi-quantitative food frequency questionnaires, extensive quantitative dietary questionnaires, and combined methods (i.e. a 7-day record on hot meals was combined with quantitative food-frequency questionnaires in Malmö, Sweden) were used to collect dietary data at baseline. These were center specific to account for local dietary habits and were either self-administered or administered in-person by trained interviewers. These data were then harmonized to obtain a standardized food list with comparable detail across countries. The dietary questionnaires and their mode of administration were described in detail in previous publications (Riboli et al., 2002, Huybrechts I et al., 2022). Then, more than 11,000 foods/ingredients/beverages were categorized into one of the four NOVA groups based on their degrees of food processing. The different NOVA groups are defined in the Additional File (Text S1). The classification of EPIC foods into NOVA groups has been described in depth elsewhere (Huybrechts et al., 2022). To account for potential changes in industrialization over time, lower, middle, and upper bound scenarios were created. The "middle bound" scenario, deemed most likely in the past 25 years, was used for the primary analysis. In the "lower bound" scenario, foods with potential for less processing were
assigned to a lower processed NOVA group, while in the “upper bound” scenario, foods with potential for more processing were assigned to a higher processed.

For each study participant, we calculated the dietary intake from each NOVA group as expressed by (1) the total absolute intake in grams/day (g/day) and (2) the total absolute intake in kcal/day (kcal/d). We also calculated the relative contribution of each NOVA food group to the total daily dietary intake in grams (%g/day) and kcal (%kcal/day). The g/d unit was considered the primary exposure because it better captures industrial foods with zero calorie content (e.g., artificially sweetened drinks) and food processing factors (e.g., neoformed contaminants or food additives).

## Covariates at recruitment

Information on lifestyle, reproductive/hormonal factors, and medical history was gathered using baseline questionnaires. All EPIC centers collected information on educational level, age at menarche, age at first full term pregnancy and parity, breastfeeding, and use of oral contraceptives and menopausal hormone therapy (MHT). Menopausal status was determined by combining different baseline information. Women who reported to have menstrual cycles, had at least nine menstrual periods over the previous 12 months, or were younger than 42 years were considered as premenopausal women. Women who reported fewer than four menses in the past year, a bilateral ovariectomy, or were older than 55 years were considered as postmenopausal women. Otherwise, women were considered perimenopausal.

Body weight and height were either measured by a health care professional or self-reported in each center. Weight and height were used to calculate BMI defined as weight in kilograms divided by height in meters squared (kg/m²). Physical activity levels were estimated using a questionnaire focused on past-year physical activity in occupational, leisure, and household domains. The Cambridge physical activity index was then created by combining occupational physical activity with time spent in physical exercise (such as cycling, swimming, and jogging) (Wareham et al., 2003). Alcohol intake in grams per day was based on the number of standard glasses of wine, beer, cider, sweet liquor, distilled spirits, or fortified wines consumed daily or weekly during the 12 months before recruitment. The Mediterranean diet score was calculated using a methodology previously described (Couto et al., 2011).

## Statistical analysis

In the main analyses, the middle-bound scenario for the NOVA classification and the absolute g/d of the four NOVA food groups were used.

Baseline characteristics were examined by quartiles of each NOVA food group intake. Multivariable Cox proportional hazards regression models were performed to estimate hazard ratios (HRs) and their corresponding 95% confidence intervals (CIs) for the associations between the intake of each NOVA food group [1 standard deviation (SD) increment] and breast cancer incidence, overall and by breast cancer subtypes. Age served as the primary time scale. All four NOVA groups were simultaneously included in the Cox model.

All models were stratified by age at recruitment in 1-year categories and study center and adjusted for potential confounding factors including educational level (none, primary school, technical/professional school, higher education), physical activity (inactive, moderately inactive, moderately active, active), height in cm (continuous), age at menarche in years (≤ 13, > 13), oral contraceptive use (never, ever, unknown), pregnancies (nulliparous, 1 or 2 children, >3 children), age at first full-term pregnancies (continuous), breastfeeding (never, ever, unknown), menopausal status (pre, peri, post-menopause), and menopausal hormone therapy (MHT) use (never, ever, unknown). We investigated whether adding different dietary-related factors/components (BMI, total energy intake, total fat, sodium intake, carbohydrate intake, Mediterranean diet, and alcohol intake) to the model changed the HR associated with NOVA groups. Only alcohol consumption modified the HRs and was therefore included in an additional model (Supplementary Table S1). Ever use of oral contraceptives and MHT and breastfeeding had >5% missing values, which were accommodated by using a “missing” category in the models. All other covariates had <5% missing values, which were replaced with the mode for categorical variables, or the median for continuous variable values observed among the subjects with complete data.

Heterogeneities according to the invasiveness status or hormonal receptor status were evaluated with competing risk analyses. In these analyses, cases with missing information on the studied subtype were excluded from the corresponding analysis and those who developed the competing breast cancer subtypes were censored at the time of occurrence (Lunn & McNeil, 1995). Heterogeneities were calculated as the deviations of logistic beta-coefficients observed in each of the subgroup relative to the overall beta-coefficient.

As subgroup and sensitivity analyses, we repeated the analyses (1) by using lower and upper bound scenarios for the NOVA classification, (2) by using consumption of NOVA groups measured as %g/d, kcal/d and %kcal/d instead of g/d, and (3) by removing alcoholic drinks (present in NOVA groups 3 and 4). When we analyzed the intake of NOVA food groups measured as the proportion
of overall daily food intake (%g/d) or daily kcal intake (%kcal/d), Cox regression analyses were performed separately for each NOVA group. Finally, we explored whether associations between NOVA groups and breast cancer risk varied by alcohol intake, BMI categories, menopausal status and country. Effect modification was evaluated by using likelihood ratio tests to compare models with and without cross-product interaction terms. Statistical analyses were conducted using SAS software (version 9.4, Copyright © 2017, SAS Institute Inc.).

## Results

During a median follow-up time of 14.9 years (13.5-16.4), 14,933 breast cancer cases were diagnosed (1,603 in situ, 13,320 invasive, and 10 of unknown invasiveness status) among the 318,686 participants. Among the invasive breast cancer cases, 9,525 had information on ER status (7,789 ER-positive and 1,736 ER-negative), 7,994 on PR status (5,268 PR-positive and 2,726 PR-negative), and 4,577 on HER2 status (901 HER2-positive and 3,676 HER2-negative). There were 573 ER+PR±HER2+, 3023 ER+PR±HER2-, 264 ER-PR-HER2-, and 419 ER-PR-HER2+ breast cancer cases.

Using the middle-bound scenario expressed as g/d, consumption of food classified as NOVA 1 contributed 74% of the total diet (Table 1). The main foods contributing to this group were coffee/tea (31%), water (19%), fruits (12%) and milk/plain yogurt (12%) (Table 2). The contribution of NOVA 2 to the total diet was 1%, with plant oils being the highest contributor to the group (37%) followed by table sugar (29%), and animal fats (28%). The contribution of processed foods (NOVA 3) to the total diet was 11%, with an important contribution of beer/wine (35%) and processed bread (26%). Overall, ultra-processed foods (NOVA 4) contributed 13% to the total diet with dairy desserts and drinks among the top group contributors (14%), followed by soft drinks (13%), ultra-processed breads (12%) and sweetened beverages (11%). The relative intake of food classified as NOVA 1 was highest in France (80%) and Denmark (79%). The relative intake of food classified as NOVA 2 was highest in Italy (3%) and Spain (2%). NOVA 3 foods were mostly highly consumed in Italy (23%), Spain (14%) and Germany (14%), while NOVA 4 foods were highly consumed in Norway (23%) and the United Kingdom (19%).

The main baseline characteristics of participants by quartiles of intake of NOVA 1, NOVA 2, NOVA 3, and NOVA 4 are presented in supplementary Tables S2, S3, S4 and S5, respectively.

Table 3 shows the associations between the middle-bound scenario of each NOVA group intake (in g/d) and breast cancer, overall and by breast cancer subtypes. Overall, intake of NOVA 1 [HR<sub>per 1 SD</sub>=0.99 (95% CI 0.97

<table><caption>**Table 1** NOVA group intake and relative and absolute contributions to total diet overall and by country.</caption><thead><tr><th rowspan="2">NOVA group<br>Country</th><th rowspan="2">N</th><th colspan="2">g/d</th><th colspan="2">%g/d</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td colspan="6"><strong>NOVA 1</strong></td></tr><tr><td>All</td><td>318,686</td><td>1979.9</td><td>853.5</td><td>74.2</td><td>10.7</td></tr><tr><td>France</td><td>67,403</td><td>2492.1</td><td>796.2</td><td>79.6</td><td>8.0</td></tr><tr><td>Italy</td><td>30,513</td><td>1113.9</td><td>361.9</td><td>63.7</td><td>10.1</td></tr><tr><td>Spain</td><td>24,850</td><td>1344.5</td><td>379.8</td><td>75.4</td><td>10.3</td></tr><tr><td>The United Kingdom</td><td>52,566</td><td>2068.3</td><td>621.4</td><td>74.3</td><td>10.3</td></tr><tr><td>The Netherlands</td><td>26,912</td><td>2216.5</td><td>619.3</td><td>76.3</td><td>8.3</td></tr><tr><td>Germany</td><td>27,379</td><td>1921.3</td><td>787.9</td><td>69.3</td><td>11.0</td></tr><tr><td>Sweden</td><td>26,368</td><td>1984.4</td><td>762.8</td><td>77.6</td><td>8.1</td></tr><tr><td>Denmark</td><td>28,720</td><td>2849.4</td><td>802.3</td><td>78.7</td><td>8.9</td></tr><tr><td>Norway</td><td>33,975</td><td>1190.6</td><td>381.1</td><td>68.1</td><td>9.7</td></tr><tr><td colspan="6"><strong>NOVA 2</strong></td></tr><tr><td>All</td><td>318,686</td><td>27.0</td><td>21.6</td><td>1.2</td><td>1.0</td></tr><tr><td>France</td><td>67,403</td><td>42.8</td><td>19.4</td><td>1.4</td><td>0.7</td></tr><tr><td>Italy</td><td>30,513</td><td>47.4</td><td>19.7</td><td>2.8</td><td>1.0</td></tr><tr><td>Spain</td><td>24,850</td><td>37.5</td><td>18.5</td><td>2.1</td><td>1.0</td></tr><tr><td>The United Kingdom</td><td>52,566</td><td>13.9</td><td>14.8</td><td>0.5</td><td>0.6</td></tr><tr><td>The Netherlands</td><td>26,912</td><td>18.5</td><td>17.9</td><td>0.7</td><td>0.7</td></tr><tr><td>Germany</td><td>27,379</td><td>25.7</td><td>20.8</td><td>1.0</td><td>0.8</td></tr><tr><td>Sweden</td><td>26,368</td><td>20.9</td><td>18.1</td><td>0.9</td><td>0.7</td></tr><tr><td>Denmark</td><td>28,720</td><td>15.1</td><td>12.3</td><td>0.4</td><td>0.4</td></tr><tr><td>Norway</td><td>33,975</td><td>12.5</td><td>8.7</td><td>0.7</td><td>0.5</td></tr><tr><td colspan="6"><strong>NOVA 3</strong></td></tr><tr><td>All</td><td>318,686</td><td>278.0</td><td>190.5</td><td>11.3</td><td>7.7</td></tr><tr><td>France</td><td>67,403</td><td>356.5</td><td>196.7</td><td>11.9</td><td>6.4</td></tr><tr><td>Italy</td><td>30,513</td><td>399.2</td><td>189.8</td><td>23.0</td><td>9.2</td></tr><tr><td>Spain</td><td>24,850</td><td>254.9</td><td>169.0</td><td>14.2</td><td>8.2</td></tr><tr><td>The United Kingdom</td><td>52,566</td><td>182.3</td><td>146.9</td><td>6.7</td><td>4.8</td></tr><tr><td>The Netherlands</td><td>26,912</td><td>248.4</td><td>137.1</td><td>8.8</td><td>4.6</td></tr><tr><td>Germany</td><td>27,379</td><td>366.5</td><td>202.5</td><td>13.9</td><td>6.8</td></tr><tr><td>Sweden</td><td>26,368</td><td>238.4</td><td>139.3</td><td>9.5</td><td>4.7</td></tr><tr><td>Denmark</td><td>28,720</td><td>302.2</td><td>220.0</td><td>8.5</td><td>5.9</td></tr><tr><td>Norway</td><td>33,975</td><td>140.9</td><td>74.5</td><td>8.3</td><td>4.2</td></tr><tr><td colspan="6"><strong>NOVA 4</strong></td></tr><tr><td>All</td><td>318,686</td><td>333.2</td><td>240.8</td><td>13.3</td><td>8.9</td></tr><tr><td>France</td><td>67,403</td><td>215.5</td><td>135.6</td><td>7.2</td><td>4.5</td></tr><tr><td>Italy</td><td>30,513</td><td>184.6</td><td>138.9</td><td>10.5</td><td>6.6</td></tr><tr><td>Spain</td><td>24,850</td><td>146.4</td><td>125.5</td><td>8.2</td><td>6.5</td></tr><tr><td>The United Kingdom</td><td>52,566</td><td>508.5</td><td>296.0</td><td>18.5</td><td>9.3</td></tr><tr><td>The Netherlands</td><td>26,912</td><td>403.6</td><td>200.9</td><td>14.2</td><td>6.7</td></tr><tr><td>Germany</td><td>27,379</td><td>417.9</td><td>254.3</td><td>15.9</td><td>8.8</td></tr><tr><td>Sweden</td><td>26,368</td><td>294.1</td><td>172.8</td><td>12.0</td><td>6.2</td></tr><tr><td>Denmark</td><td>28,720</td><td>434.4</td><td>280.0</td><td>12.3</td><td>7.3</td></tr><tr><td>Norway</td><td>33,975</td><td>386.9</td><td>163.2</td><td>22.8</td><td>8.8</td></tr></tbody></table>

NOVA 1 Unprocessed/minimally processed foods, NOVA 2 Processed culinary ingredients, NOVA 3 Processed foods, NOVA 4 Ultra-processed foods, SD Standard Deviation
– 1.01)], NOVA 2 [HR<sub>per 1 SD</sub>=1.01 (0.98 – 1.03)], and NOVA 4 [HR<sub>per 1 SD</sub>=1.01 (0.99 – 1.03)] were not associated with breast cancer risk. However, intake of processed foods (NOVA 3) was associated with a higher risk of breast cancer [HR<sub>per 1 SD</sub>=1.05 (1.03 – 1.07)]. Estimates did not differ by invasiveness or hormone receptor status (Table 3, P<sub>homogeneity</sub>≥ 0.11). When the model was further adjusted for alcohol intake (Table 4), the positive association between NOVA 3 and breast cancer risk was attenuated and no longer statistically significant [HR<sub>per 1 SD</sub>=1.01 (0.98 – 1.03)]. Furthermore, when alcoholic drinks were excluded from NOVA 3 the association with breast cancer risk was also null [HR<sub>per 1 SD</sub>=0.99 (0.97 – 1.01), Table 5]. Associations were similar when models were stratified by alcohol intake, BMI at recruitment (P<sub>interaction</sub> ≥ 0.17, Table 6) or menopausal status (Table S6).

In secondary analyses using %g/d, kcal/d or %kcal/d as the exposure, the results were consistent with those obtained in the main analyses (Supplementary Table S7). However, when using %g/d as an exposure variable, a higher intake of NOVA 1 was associated with a slightly lower risk of breast cancer [HR<sub>per 1 SD</sub>=0.96 (0.94 - 0.98)]. Nevertheless, this association was no longer statistically significant when the model was further adjusted for alcohol intake [HR<sub>per 1 SD</sub>=0.98 (0.94-1.00), Supplementary Table S7]. In addition, the results were similar when we used lower and upper bound scenarios (data not shown). Finally, no heterogeneity was reported by country (P<sub>homogeneity</sub> ≥ 0.09, Supplementary Table S8).

## Discussion

In this large-scale prospective analysis, we found a positive association between the consumption of processed foods and breast cancer, which was likely driven by alcohol – an already established risk factor for breast cancer. The association between the degree of food processing and breast cancer risk did not differ by breast cancer subtype, menopausal status, alcohol intake or BMI.

In this study, no associations were found between the consumption of food included in the NOVA 1 group and breast cancer risk, overall or by breast cancer subtypes when the absolute values of intake were evaluated (g/day or kcal/day). Although a slight inverse association was reported when the %g/day values were used, this association disappeared when models were further adjusted for alcohol intake. Furthermore, because we only observed an inverse association when the variable was expressed as %g/d, these results might be because individuals who consumed more food from NOVA group 1 also consumed less food from NOVA group 3. Only one cohort study (Fiolet et al., 2018) and one case-control study (Jacobs et al., 2022) investigated associations between

NOVA 1 and breast cancer risk and reported an inverse association (using %g/d and %kcal/d, respectively). Although NOVA 1 foods have low energy density and are rich in phytochemicals (carotenoids, flavonoids, dietary fiber), vitamins and minerals, known to be anticancerogenic (Bakker et al., 2016), our study does not support the hypothesis of a lower breast cancer risk with higher intake of NOVA 1 foods. In addition, we found no evidence of an association between NOVA 2 and breast cancer risk, as also reported in a South African case-control study (Jacobs et al., 2022).

We observed a positive association between NOVA 3 intake and breast cancer risk. To our knowledge, no previous population study has reported a positive association between NOVA 3 and breast cancer risk (Fiolet et al., 2018; Jacobs et al., 2022). Interestingly, in our study population, the positive association disappeared when models were adjusted for alcohol consumption or when alcoholic drinks were excluded from the NOVA groups. Indeed, alcohol is an established risk factor for breast cancer and, as such, could drive the positive association between processed foods and breast cancer risk. Of note, in this study we observed that, on average, beer and wine made up 35% of NOVA 3 g/day intake.

Furthermore, we found no evidence of an association between the consumption of NOVA 4 and breast cancer risk. Other studies have also reported no association between NOVA 4 intake and breast cancer risk (Jacobs et al., 2022; Romaguera et al., 2021). However, our results differ from those from the NutriNet-Santé French cohort and two case-control studies, which reported a positive association between the consumption of NOVA 4 and breast cancer risk (Fiolet et al., 2018; Queiroz et al., 2018; Romieu et al., 2022). It has been suggested that NOVA 4 foods may increase breast cancer risk through several factors such as their high energy density due to added sugars and fats, the presence of a variety of additives, preservatives and processing contaminants (e.g. acrylamide, *trans*-fatty acids, endocrine disrupters, etc.) or lack of fiber, proteins, and other components that are associated with fullness and satisfaction, leading individuals to eat more in an attempt to feel satisfied/saturated (Friedman, 2015; Luiten et al., 2016; Moubarac et al., 2013; Pouzou et al., 2018). In addition, we might have expected to observe a positive association between NOVA 4 and breast cancer risk due to alcoholic distilled drinks, however, in this study population alcoholic distilled drinks made up 2.2% of NOVA 4 g/day intake. The lack of association between NOVA 4 and breast cancer risk in the current study might be explained by the fact that the consumption of NOVA 4 in EPIC was quite low as this was based on dietary intakes at recruitment (during the nineties); since then, the consumption of NOVA
<table><caption>**Table 2** Absolute intakes and relative contributions of food to total diet and to each NOVA group</caption><thead><tr><th rowspan="2">NOVA groups</th><th rowspan="2">Food groups</th><th colspan="2">g/d</th><th colspan="2">% in total diet</th><th>% within the food group</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>Mean</th></tr></thead><tbody><tr><td rowspan="15">1</td><td>Water</td><td>382.23</td><td>531.38</td><td>11.81</td><td>14.65</td><td>19.3</td></tr><tr><td>Fruit</td><td>232.79</td><td>176.68</td><td>9.68</td><td>7.68</td><td>11.8</td></tr><tr><td>Milk and plain yoghurt</td><td>234.29</td><td>205.38</td><td>9.22</td><td>7.63</td><td>11.8</td></tr><tr><td>Cereal, grains and flour made from these foods</td><td>40.89</td><td>50.14</td><td>1.67</td><td>1.95</td><td>2.1</td></tr><tr><td>Potatoes</td><td>75.8</td><td>58.63</td><td>3.15</td><td>2.65</td><td>3.8</td></tr><tr><td>Fresh pasta</td><td>34.52</td><td>40.19</td><td>1.57</td><td>2.15</td><td>1.7</td></tr><tr><td>Beans, lentils and chickpeas</td><td>18.47</td><td>24.06</td><td>0.79</td><td>1.14</td><td>0.9</td></tr><tr><td>Vegetables</td><td>186.24</td><td>117.89</td><td>7.54</td><td>4.74</td><td>9.4</td></tr><tr><td>Nuts and Seeds</td><td>2.28</td><td>5.22</td><td>0.09</td><td>0.21</td><td>0.1</td></tr><tr><td>Eggs</td><td>17.55</td><td>16.37</td><td>0.72</td><td>0.66</td><td>0.9</td></tr><tr><td>Poultry</td><td>17.3</td><td>17.63</td><td>0.75</td><td>0.86</td><td>0.9</td></tr><tr><td>Red meat</td><td>41.65</td><td>34.35</td><td>1.7</td><td>1.44</td><td>2.1</td></tr><tr><td>Fish</td><td>23.56</td><td>27.11</td><td>1.07</td><td>1.47</td><td>1.2</td></tr><tr><td>Seafood</td><td>3.01</td><td>5.56</td><td>0.14</td><td>0.26</td><td>0.2</td></tr><tr><td>Fungi</td><td>5.93</td><td>8.81</td><td>0.22</td><td>0.32</td><td>0.3</td></tr><tr><td>Coffee/tea</td><td>606.68</td><td>462.33</td><td>21.97</td><td>13.6</td><td>30.6</td></tr><tr><td>Fruit juice fresh and smoothies</td><td>10.58</td><td>25.65</td><td>0.4</td><td>1.02</td><td>0.5</td></tr><tr><td>Fruit juice UHT or pasteurised</td><td>42.78</td><td>70.3</td><td>1.66</td><td>2.75</td><td>2.2</td></tr><tr><td>Homemade broth</td><td>3.31</td><td>7.93</td><td>0.11</td><td>0.28</td><td>0.2</td></tr><tr><td rowspan="7">2</td><td>Table sugar</td><td>7.8</td><td>11.62</td><td>0.33</td><td>0.51</td><td>28.9</td></tr><tr><td>Plant oil</td><td>9.97</td><td>11.95</td><td>0.46</td><td>0.62</td><td>36.9</td></tr><tr><td>Animal fats</td><td>7.52</td><td>9.76</td><td>0.3</td><td>0.39</td><td>27.9</td></tr><tr><td>Other processed culinary ingredients</td><td>1.3</td><td>2.71</td><td>0.04</td><td>0.09</td><td>4.8</td></tr><tr><td>Table salt</td><td>0.4</td><td>0.65</td><td>0.01</td><td>0.02</td><td>1.5</td></tr><tr><td>Cheese</td><td>37.35</td><td>33.57</td><td>1.54</td><td>1.45</td><td>13.4</td></tr><tr><td>Salted, smoked or canned meat, without additives</td><td>5.83</td><td>12.64</td><td>0.28</td><td>0.7</td><td>2.1</td></tr><tr><td rowspan="12">3</td><td>Salted, smoked or canned fish</td><td>6.66</td><td>9.69</td><td>0.28</td><td>0.42</td><td>2.4</td></tr><tr><td>Processed bread</td><td>72.83</td><td>80.59</td><td>3.14</td><td>3.76</td><td>26.2</td></tr><tr><td>Vegetables and other plant foods preserved</td><td>16.49</td><td>26.1</td><td>0.69</td><td>1.2</td><td>5.9</td></tr><tr><td>Legumes preserved</td><td>4.41</td><td>11.32</td><td>0.19</td><td>0.56</td><td>1.6</td></tr><tr><td>Fruit preserved</td><td>17.8</td><td>25.28</td><td>0.68</td><td>0.94</td><td>6.4</td></tr><tr><td>Nuts salted and nut spreads</td><td>1.57</td><td>4.28</td><td>0.06</td><td>0.18</td><td>0.6</td></tr><tr><td>Beer and Wine</td><td>96.93</td><td>150.13</td><td>3.68</td><td>5.31</td><td>34.9</td></tr><tr><td>Condensed milk, yogurt plain sweetened</td><td>6.15</td><td>17.53</td><td>0.24</td><td>0.66</td><td>2.2</td></tr><tr><td>Bread crumbs</td><td>0.21</td><td>0.59</td><td>0.01</td><td>0.02</td><td>0.1</td></tr><tr><td>Meringue, non ultra-processed bakeries</td><td>8.74</td><td>17.42</td><td>0.37</td><td>0.7</td><td>3.1</td></tr><tr><td>Cheese</td><td>3.06</td><td>6.1</td><td>0.14</td><td>0.3</td><td>1.1</td></tr></tbody></table>
<table><caption>**Table 2** (continued)</caption><thead><tr><th rowspan="2">NOVA groups</th><th rowspan="2">Food groups</th><th colspan="2">g/d</th><th colspan="2">% in total diet</th><th>% within the food group</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th><th>Mean</th></tr></thead><tbody><tr><td rowspan="20">4</td><td>Ultra-processed breads</td><td>40.57</td><td>62.77</td><td>1.72</td><td>2.95</td><td>12.1</td></tr><tr><td>Pastries, buns, and cakes</td><td>21.09</td><td>26.95</td><td>0.87</td><td>1.2</td><td>6.3</td></tr><tr><td>Biscuits</td><td>14.48</td><td>18.15</td><td>0.63</td><td>0.89</td><td>4.3</td></tr><tr><td>Breakfast cereals</td><td>4.62</td><td>9.76</td><td>0.17</td><td>0.37</td><td>1.4</td></tr><tr><td>Ice cream, ice pops and frozen yogurts</td><td>7.11</td><td>11.53</td><td>0.31</td><td>0.56</td><td>2.1</td></tr><tr><td>Industrial desserts</td><td>1.26</td><td>8.96</td><td>0.05</td><td>0.33</td><td>0.4</td></tr><tr><td>Packaged salty snacks</td><td>1.75</td><td>5.14</td><td>0.07</td><td>0.21</td><td>0.5</td></tr><tr><td>Potato products</td><td>6.95</td><td>15.08</td><td>0.26</td><td>0.58</td><td>2.1</td></tr><tr><td>Pizza and focaccia</td><td>6.55</td><td>10.41</td><td>0.29</td><td>0.52</td><td>1.9</td></tr><tr><td>Filled Pasta (industrial pasta)</td><td>3.09</td><td>7.96</td><td>0.14</td><td>0.4</td><td>0.9</td></tr><tr><td>Instant and canned soups</td><td>9</td><td>19.67</td><td>0.35</td><td>0.78</td><td>2.7</td></tr><tr><td>Dairy substitute products</td><td>3.2</td><td>32.73</td><td>0.13</td><td>1.31</td><td>1</td></tr><tr><td>Processed cheese</td><td>3.05</td><td>6.68</td><td>0.14</td><td>0.31</td><td>0.9</td></tr><tr><td>Sauces, dressing and gravies</td><td>9.39</td><td>12.12</td><td>0.34</td><td>0.43</td><td>2.8</td></tr><tr><td>Vegetable spread and products</td><td>0.19</td><td>1.18</td><td>0.01</td><td>0.05</td><td>0.1</td></tr><tr><td>Soft drinks</td><td>43.32</td><td>109.71</td><td>1.8</td><td>4.35</td><td>12.9</td></tr><tr><td>Dairy desserts and drinks</td><td>46.27</td><td>69.3</td><td>1.77</td><td>2.5</td><td>13.8</td></tr><tr><td>Sweetened beverages</td><td>37.88</td><td>113.82</td><td>1.31</td><td>3.53</td><td>11.3</td></tr><tr><td>Beverages (dry weight) (ex. Coffee powder, chocolate powder and milk powder)</td><td>0.86</td><td>3.71</td><td>0.03</td><td>0.14</td><td>0.3</td></tr><tr><td>Alcoholic distilled drinks</td><td>7.5</td><td>20.23</td><td>0.27</td><td>0.71</td><td>2.2</td></tr><tr><td>Artificial sweeteners</td><td>0.48</td><td>2.55</td><td>0.02</td><td>0.08</td><td>0.1</td></tr><tr><td>Sweet snacks</td><td>11.32</td><td>18.44</td><td>0.43</td><td>0.7</td><td>3.4</td></tr><tr><td>Processed meat</td><td>31.76</td><td>29.04</td><td>1.38</td><td>1.47</td><td>9.4</td></tr><tr><td>Meat alternatives</td><td>0.85</td><td>3.93</td><td>0.03</td><td>0.15</td><td>0.3</td></tr><tr><td>Nutrition powders and drinks</td><td>0.01</td><td>0.47</td><td>0</td><td>0.02</td><td>0</td></tr><tr><td>Margarine</td><td>10.6</td><td>13.64</td><td>0.43</td><td>0.59</td><td>3.2</td></tr><tr><td>Ready meals</td><td>5.12</td><td>10.77</td><td>0.17</td><td>0.36</td><td>1.5</td></tr><tr><td>Alcohol-free versions of alcoholic beverages</td><td>1.62</td><td>19.3</td><td>0.06</td><td>0.7</td><td>0.5</td></tr><tr><td>Ultra-processed vegetables and legumes</td><td>3.31</td><td>12.18</td><td>0.12</td><td>0.45</td><td>1.0</td></tr><tr><td>Rice-based dishes</td><td>0</td><td>0.13</td><td>0</td><td>0.01</td><td>0</td></tr></tbody></table>

NOVA 1 Unprocessed/minimally processed foods, NOVA 2 Processed culinary ingredients, NOVA 3 Processed foods, NOVA 4 Ultra-processed foods

4 has replaced the consumption of other NOVA groups (e.g. recipes that were made at home in the 1990s may currently be industrially processed), which may bias the associations with breast cancer risk. Indeed, NOVA 4 has been recently suggested to represent up to 60% of total

daily energy intake in some countries of the European area such as the UK (Rauber et al., 2018) while in the current EPIC study population, NOVA 4 contributed to 31% of total energy intake (data not shown). However, when we classified food products based on the modern/current
<table><caption>**Table 3** Associations between NOVA groups (in g/d) and breast cancer risk, overall and by breast cancer subtypes</caption><thead><tr><th rowspan="2">Breast cancer subtypes</th><th rowspan="2">N Cases</th><th colspan="2">NOVA 1</th><th colspan="2">NOVA 2</th><th colspan="2">NOVA 3</th><th colspan="2">NOVA 4</th></tr><tr><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th></tr></thead><tbody><tr><td>Overall</td><td>14,933</td><td>0.99 (0.97-1.01)</td><td></td><td>1.01 (0.98-1.03)</td><td></td><td>1.05 (1.03-1.07)</td><td></td><td>1.01 (0.99-1.03)</td><td></td></tr><tr><td>In situ</td><td>1603</td><td>0.97 (0.91-1.04)</td><td>0.63</td><td>1.00 (0.94-1.07)</td><td>0.87</td><td>1.01 (0.96-1.07)</td><td>0.18</td><td>1.01 (0.99-1.03)</td><td>0.87</td></tr><tr><td>Invasive</td><td>13,320</td><td>0.99 (0.97-1.01)</td><td></td><td>1.01 (0.98-1.03)</td><td></td><td>1.05 (1.03-1.07)</td><td></td><td>1.01 (0.95-1.07)</td><td></td></tr><tr><td>Invasive ER+</td><td>7789</td><td>0.99 (0.96-1.02)</td><td>0.54</td><td>1.00 (0.97-1.03)</td><td>0.29</td><td>1.06 (1.04-1.08)</td><td>0.86</td><td>1.00 (0.98-1.03)</td><td>0.11</td></tr><tr><td>Invasive ER-</td><td>1736</td><td>1.01 (0.95-1.07)</td><td></td><td>1.04 (0.98-1.10)</td><td></td><td>1.06 (1.01-1.11)</td><td></td><td>1.06 (1.00-1.11)</td><td></td></tr><tr><td>Invasive PR+</td><td>5268</td><td>0.99 (0.96-1.03)</td><td>0.56</td><td>1.01 (0.98-1.04)</td><td>0.34</td><td>1.05 (1.03-1.08)</td><td>0.51</td><td>0.99 (0.95-1.02)</td><td>0.15</td></tr><tr><td>Invasive PR-</td><td>2726</td><td>0.98 (0.93-1.03)</td><td></td><td>0.98 (0.93-1.03)</td><td></td><td>1.07 (1.03-1.11)</td><td></td><td>1.03 (0.98-1.09)</td><td></td></tr><tr><td>Invasive HER2+</td><td>901</td><td>0.99 (0.91-1.08)</td><td>0.97</td><td>1.00 (0.92-1.09)</td><td>0.70</td><td>1.04 (0.98-1.12)</td><td>0.57</td><td>1.04 (0.97-1.13)</td><td>0.39</td></tr><tr><td>Invasive HER2-</td><td>3676</td><td>0.99 (0.95-1.03)</td><td></td><td>1.02 (0.98-1.06)</td><td></td><td>1.07 (1.03-1.10)</td><td></td><td>1.01 (0.97-1.03)</td><td></td></tr></tbody></table>

CI Confidence Interval, ER Estrogen Receptor, HER2 Human epidermal growth factor receptor 2, HR Hazard ratio, NOVA 1 unprocessed/minimally processed foods, NOVA 2 Processed culinary ingredients, NOVA 3 Processed foods, NOVA 4 Ultra-processed foods, PR Progesterone receptor, SD Standard deviation

<table><caption>**Table 4** Associations between NOVA groups (in g/d) and breast cancer risk with further adjustment for alcohol consumption</caption><thead><tr><th rowspan="2">Breast cancer subtypes</th><th rowspan="2">N Cases</th><th colspan="2">NOVA 1</th><th colspan="2">NOVA 2</th><th colspan="2">NOVA 3</th><th colspan="2">NOVA 4</th></tr><tr><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th></tr></thead><tbody><tr><td>Overall</td><td>14,933</td><td>0.99 (0.97-1.01)</td><td></td><td>1.01 (0.99-1.03)</td><td></td><td>1.01 (0.98-1.03)</td><td></td><td>1.00 (0.98-1.03)</td><td></td></tr><tr><td>In situ</td><td>1603</td><td>0.98 (0.91-1.04)</td><td>0.61</td><td>1.00 (0.94-1.07)</td><td>0.84</td><td>0.98 (0.90-1.07)</td><td>0.51</td><td>1.01 (0.95-1.08)</td><td>0.89</td></tr><tr><td>Invasive</td><td>13,320</td><td>0.99 (0.97-1.02)</td><td></td><td>1.01 (0.99-1.03)</td><td></td><td>1.01 (0.98-1.04)</td><td></td><td>1.01 (0.99-1.03)</td><td></td></tr><tr><td>Invasive ER+</td><td>7789</td><td>0.99 (0.96-1.02)</td><td>0.58</td><td>1.00 (0.97-1.04)</td><td>0.33</td><td>1.00 (0.97-1.04)</td><td>0.50</td><td>1.00 (0.98-1.03)</td><td>0.11</td></tr><tr><td>Invasive ER-</td><td>1736</td><td>1.01 (0.95-1.08)</td><td></td><td>1.04 (0.98-1.11)</td><td></td><td>1.03 (0.96-1.12)</td><td></td><td>1.06 (1.00-1.11)</td><td></td></tr><tr><td>Invasive PR+</td><td>5268</td><td>1.00 (0.96-1.04)</td><td>0.51</td><td>1.01 (0.98-1.05)</td><td>0.30</td><td>1.00 (0.95-1.06)</td><td>0.30</td><td>0.99 (0.95-1.03)</td><td>0.15</td></tr><tr><td>Invasive PR-</td><td>2726</td><td>0.98 (0.93-1.03)</td><td></td><td>0.98 (0.94-1.03)</td><td></td><td>1.04 (0.98-1.11)</td><td></td><td>1.03 (0.98-1.09)</td><td></td></tr><tr><td>Invasive HER2+</td><td>901</td><td>0.99 (0.91-1.08)</td><td>0.91</td><td>1.00 (0.92-1.08)</td><td>0.56</td><td>1.07 (0.96-1.19)</td><td>0.29</td><td>1.04 (0.97-1.12)</td><td>0.44</td></tr><tr><td>Invasive HER2-</td><td>3676</td><td>1.00 (0.95-1.04)</td><td></td><td>1.02 (0.98-1.07)</td><td></td><td>1.00 (0.95-1.05)</td><td></td><td>1.01 (0.97-1.05)</td><td></td></tr></tbody></table>

CI Confidence Interval, ER Estrogen Receptor, HER2 Human epidermal growth factor receptor 2, HR Hazard ratio, NOVA 1 Unprocessed/minimally processed foods, NOVA 2 Processed culinary ingredients, NOVA 3 Processed foods, NOVA 4 Ultra-processed foods; PR, progesterone receptor, SD Standard deviation

food environment (upper bound scenario), the association between NOVA 4 and breast cancer risk was still null.

Finally, in our study, the associations between the degree of food processing and breast cancer risk did not differ by invasiveness and hormone receptor status. Our results are consistent with a previous study that investigated the association between NOVA 4 and breast cancer risk by molecular status suggesting no differences in risk estimates [ER+ or PR+: OR<sub>10% increase</sub>=1.04 (0.96 – 1.13); HER2+: OR<sub>10% increase</sub>=0.96 (0.84 – 1.10); ER-PR-HER2-: OR<sub>10% increase</sub>=0.93 (0.75 – 1.15)] (Romaguera et al., 2021). However, another study reported a positive association

between NOVA 4 and ER+ breast cancer [OR<sub>T3vsT1</sub>=2.44 (1.01 - 5.90)] and no clear association with ER- breast cancer [OR<sub>T3vsT1</sub>=1.87 (0.43 - 8.13)] (Romieu et al., 2022). No other studies reported results for other NOVA groups by invasiveness and molecular status of breast cancer, which makes any comparison with our results challenging.

Our study has several strengths including the prospective design, the multicenter aspect, the long-term follow-up, and the availability of a comprehensive assessment of participant characteristics, as well as the large number of incident breast cancer cases. We had self-reported data on lifestyle, reproductive, and medical factors and were

<footnote>a Models were stratified by age and center and adjusted for education, height, physical activity, age at menarche, oral contraceptive use, age at first full-term pregnancy, parity, breastfeeding, menopausal status and menopausal hormone therapy use. Each NOVA group was mutually adjusted for the other NOVA groups</footnote>

<footnote>a Models were stratified by age and center and adjusted for education, height, physical activity, age at menarche, oral contraceptive use, age at first full-term pregnancy, parity, breastfeeding, menopausal status, menopausal hormone therapy use and alcohol consumption. Each NOVA group was mutually adjusted for the other NOVA groups</footnote>
<table><caption>**Table 5** Associations between NOVA 3 intake after excluding alcoholic drinks (in g/d) and breast cancer risk</caption><thead><tr><th>Breast cancer subtypes</th><th>HR per 1 SD (95%CI)<sup>a</sup></th><th>P<sub>homogeneity</sub></th></tr></thead><tbody><tr><td>Overall</td><td>0.99 (0.97-1.01)</td><td></td></tr><tr><td>Invasive</td><td>0.99 (0.97-1.01)</td><td>0.78</td></tr><tr><td>In situ</td><td>0.98 (0.92-1.04)</td><td></td></tr><tr><td>Invasive ER+</td><td>0.98 (0.95-1.00)</td><td>0.08</td></tr><tr><td>Invasive ER-</td><td>1.03 (0.98-1.10)</td><td></td></tr><tr><td>Invasive PR+</td><td>0.99 (0.96-1.02)</td><td>0.91</td></tr><tr><td>Invasive PR-</td><td>0.99 (0.95-1.04)</td><td></td></tr><tr><td>Invasive HER2+</td><td>1.04 (0.96-1.12)</td><td>0.15</td></tr><tr><td>Invasive HER2-</td><td>0.97 (0.93-1.01)</td><td></td></tr></tbody></table>

CI Confidence Interval, ER Estrogen Receptor, HER2 Human epidermal growth factor receptor 2, HR Hazard ratio, NOVA 3 Processed foods, PR Progesterone receptor, SD Standard deviation

ª Models were stratified by age and center and adjusted for education, height, physical activity, age at menarche, oral contraceptive use, age at first full-term pregnancy, parity, breastfeeding, menopausal status, and menopausal hormone therapy use. Each NOVA group was mutually adjusted for the other NOVA groups

therefore able to consider a wide range of potential founders. We also stratified the analyses according to menopausal status, alcohol intake, obesity breast cancer subtype. However, it should be noted that the statistical power for some of these stratified analyses was rather low (e.g. for some of the breast cancer subtypes), therefore, these results should be interpreted with caution. The major limitation of the study is that dietary data were only collected at baseline (in the 1990s) while the food environment changed in the intervening years, exposing the EPIC participants to potentially different degrees of

food processing over the course of their follow-up. However, three different scenarios were created considering that the food environment may have changed over time compared to the baseline. The lower and upper bound scenarios were used in sensitivity analyses to explore the potential impact of further industrialization of food products and of changes in consumer habits to convenience foods over time, and results were virtually unchanged. The fact that dietary data were collected only once may cause random measurement error and may fail to reflect long-term habits; any such bias would likely lead to an underestimation of true associations (i.e. regression dilution bias) (Clarke et al., 1999; Hutcheon et al., 2010). However, it is noteworthy that recent analyses comparing dietary follow-up data in some of the EPIC countries demonstrate only minor changes in dietary intakes among the EPIC participants, potentially due to the relatively older age of the participants included in the cohort (unpublished data). Finally, the dietary questionnaires used in EPIC were not designed to identify different food processing categories. Therefore, several assumptions had to be made when insufficient information about the processing of the food item was available, potentially contributing to measurement error. However, byproducts of processing (e.g. trans fatty acids or syringol metabolites) have been positively associated with NOVA group 4 in EPIC which is a sign of a good measurement of ultra-processed foods in EPIC (Huybrechts et al., 2022).

<table><caption>**Table 6** Associations between NOVA intake (in g/d) and breast cancer risk, stratified by alcohol intake and BMI at recruitment</caption><thead><tr><th>Characteristics at recruitment</th><th>N cases</th><th>NOVA 1<br>HR per 1 SD (95%CI)<sup>a</sup></th><th>NOVA 2<br>HR per 1 SD (95%CI)<sup>a</sup></th><th>NOVA 3<br>HR per 1 SD (95%CI)<sup>a</sup></th><th>NOVA 4<br>HR per 1 SD (95%CI)<sup>a</sup></th></tr></thead><tbody><tr><td><strong>Alcohol intake (g/d)</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Non-drinkers</td><td>2014</td><td>0.99 (0.93-1.05)</td><td>1.00 (0.94-1.05)</td><td>1.00 (0.91-1.10)</td><td>1.01 (0.95-1.07)</td></tr><tr><td>&gt; 0 - ≤3</td><td>4311</td><td>0.97 (0.93-1.01)</td><td>1.03 (0.99-1.07)</td><td>1.00 (0.93-1.07)</td><td>1.01 (0.98-1.05)</td></tr><tr><td>&gt; 3 - ≤12</td><td>4594</td><td>1.02 (0.98-1.06)</td><td>1.02 (0.98-1.06)</td><td>0.99 (0.93-1.05)</td><td>1.00 (0.97-1.04)</td></tr><tr><td>&gt; 12 - ≤24</td><td>2298</td><td>0.98 (0.93-1.03)</td><td>0.97 (0.92-1.03)</td><td>1.05 (0.99-1.12)</td><td>1.01 (0.96-1.07)</td></tr><tr><td>&gt;24</td><td>1716</td><td>0.99 (0.93-1.05)</td><td>0.99 (0.93-1.05)</td><td>1.02 (0.99-1.06)</td><td>1.01 (0.95-1.07)</td></tr><tr><td>P<sub>interaction</sub></td><td></td><td>0.79</td><td>0.17</td><td>0.68</td><td>0.52</td></tr><tr><td><strong>Body mass index</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>&lt;25</td><td>8871</td><td>0.99 (0.97-1.02)</td><td>1.01 (0.98-1.04)</td><td>1.05 (1.03-1.08)</td><td>1.01 (0.98-1.04)</td></tr><tr><td>≥25</td><td>6062</td><td>0.97 (0.94-1.01)</td><td>1.01 (0.97-1.04)</td><td>1.05 (1.02-1.08)</td><td>1.00 (0.97-1.03)</td></tr><tr><td>P<sub>interaction</sub></td><td></td><td>0.57</td><td>0.31</td><td>0.27</td><td>0.88</td></tr></tbody></table>

CI Confidence Interval, NOVA 1 Unprocessed/minimally processed foods, NOVA 2 Processed culinary ingredients, NOVA 3 Processed foods, NOVA 4 Ultra-processed foods, SD Standard deviation

ª Models were stratified by age and center and adjusted for education, height, physical activity, age at menarche, oral contraceptive use, age at first full-term pregnancy, parity, breastfeeding, menopausal status and menopausal hormone therapy use. Each NOVA group was mutually adjusted for the other NOVA groups
## Conclusion

This large-scale prospective analysis among European women suggests that the positive association between processed food intake and breast cancer risk was likely driven by alcoholic beverage consumption. Other degrees of food processing were not associated with breast cancer risk.

### Abbreviations

BMI Body Mass Index

CIs Confidence Intervals

EPIC European Prospective Investigation into Cancer and Nutrition

ER Estrogen Receptor

HER2 Human Epidermal Growth Factor Receptor 2

HRS Hazard Ratios

IARC International Agency for Research on Cancer

ICD-O International Classification of Diseases for Oncology

MHT Menopausal Hormone Therapy

OR Odd Ratio

PR Progesterone Receptor

SD Standard Deviation

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s43014-024-00264-2.

Supplementary Material 1.

### Acknowledgements

The authors thank all participants in the EPIC cohort for their invaluable contribution to the study. The authors also thank Bertrand Hémon (IARC) for his previous help with the EPIC database and the National Institute for Public Health and the Environment (RIVM), Bilthoven, the Netherlands, for their contribution and ongoing support to the EPIC Study. We acknowledge the use of data from the EPIC-Oxford cohort.

### Disclaimer

Where authors are identified as personnel of the International Agency for Research on Cancer/World Health Organization, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy, or views of the International Agency for Research on Cancer/World Health Organization.

### Authors' contributions

The food processing indicators were generated by FR, RBL, IH, GN and CC. The statistical analyses were performed by MC. MC, SY and IH were responsible for drafting the manuscript. All authors had the opportunity to comment on the analysis and interpretation of the findings and approved the final version for publication.

### Funding

FMB was supported by a Wellcome Trust PhD studentship in Molecular, Genetic and Lifecourse Epidemiology (224982/Z/22/Z). The coordination of EPIC is financially supported by IARC and by the Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London which has additional infrastructure support provided by the NIHR Imperial Biomedical Research Centre (BRC). The national cohorts are supported by: Danish Cancer Society (Denmark); Ligue Contre le Cancer, Institut Gustave Roussy, Mutuelle Générale de l'Éducation Nationale, Institut National de la Santé et de la Recherche Médicale (Inserm) (France); German Cancer Aid, German Cancer Research Center (DKFZ), German Institute of Human Nutrition Potsdam-Rehbruecke (DIfE), Federal Ministry of Education and Research (BMBF) (Germany); Associazione Italiana per la Ricerca sul Cancro-AIRC-Italy, Compagnia di San Paolo and National Research Council (Italy); Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek

Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands (The Netherlands); Health Research Fund (FIS) - Instituto de Salud Carlos III (ISCIII), Regional Governments of Andalucía, Asturias, Basque Country, Murcia and Navarra, and the Catalan Institute of Oncology - ICO (Spain); Swedish Cancer Society, Swedish Research Council and County Councils of Skåne and Västerbotten (Sweden); Cancer Research UK (14136 to EPIC-Norfolk; C8221/A29017 to EPIC-Oxford), Medical Research Council (1000143 to EPIC-Norfolk; MR/M012190/1 to EPIC-Oxford). (United Kingdom).
The work reported in this paper was performed during Agnès Fournier's term as a Visiting Scientist at the IARC. This work was supported by Cancer Research UK C33493/A29678; l'Institut National du Cancer (« INCa ») APPEL À PROJET SHS 2021 Grant number: INCa N°2021-138 and APPEL À PROJET 2021 Grant number: INCa N°SHSESP 21-074 and the World Cancer Research Fund (WCRF) Grant reference number: IIG_FULL_2020_033. Researchers were independent from the funders. Funders had no role in the collection, analysis, and interpretation of data, the writing of the report, and the decision to submit the article for publication.

### Availability of data and materials

For information on how to apply for gaining access to EPIC data and/or bio-specimens, please follow the instructions at https://epic.iarc.fr/access/submit_appl_access.php.

## Declarations

### Ethics approval and consent to participate

The EPIC cohort was conducted in accordance with the Declaration of Helsinki and received ethical approval from the IARC ethics committee and all study centers. Participants provided written informed consent.

### Consent for publication

Not applicable.

### Competing interests

The authors declare no potential conflicts of interest.

### Author details

¹Nutrition and Metabolism Branch, International Agency for Research on Cancer, Lyon, France. ²Université Paris-Saclay, UVSQ, Inserm, Gustave Roussy, CESP, 94805 Villejuif, France. ³Université Sorbonne Paris Nord and Université Paris Cité, INSERM, INRAE, CNAM, Center of Research in Epidemiology and Statistics (CRESS), Nutritional Epidemiology Research Team (EREN), F-93017 Bobigny, France. ⁴Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, United Kingdom. ⁵MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom. ⁶Department of Preventive Medicine, School of Medicine, University of São Paulo, São Paulo, Brazil. ⁷Department of Community Medicine, Faculty of Health Sciences,UiT - the Arctic University of Norway, Tromsø, Norway. ⁸Danish Cancer Society Research Center, Copenhagen, Denmark. ⁹Department of Public Health, University of Copenhagen, Copenhagen, Denmark. ¹⁰Nutrition Epidemiology, Department of Clinical Sciences Malmö, Faculty of Medicine, Lund University, Lund, Sweden. ¹¹Department of Clinical Sciences Lund, Oncology, Lund University and Skåne University Hospital, Lund, Sweden. ¹²Institute for Cancer Research, Prevention and Clinical Network (ISPRO), Florence, Italy. ¹³Department of Clinical and Biological Sciences and Centre for Biostatistics, Epidemiology, and Public Health (C-BEPH), Turin, Italy. ¹⁴Dipartimento di medicina clinica e chirurgia Federico ii university, Naples, Italy. ¹⁵Department of Cancer Epidemiology, German Cancer research Center (DKFZ), Heidelberg, Germany. ¹⁶Department of Epidemiology, Murcia Regional Health Council, IMIB-Arrixaca, Murcia, Spain. ¹⁷Spanish Consortium for Research on Epidemiology and Public Health (CIBER-ESP), Instituto de Salud Carlos III, Madrid, Spain. ¹⁸Epidemiology of Chronic and Communicable Diseases Group, Biodonostia Health Research Institute, San Sebastián, Spain. ¹⁹Sub Directorate for Public Health and Addictions of Gipuzkoa, Ministry of Health of the Basque Government, San Sebastian, Spain. ²⁰Department of Clinical Sciences, Pediatrics, Umeå University, Umeå, Sweden. ²¹Unit of Nutrition and Cancer, Cancer Epidemiology Research Program, Catalan Institute of Oncology (ICO), Av. Granvia de L'Hospitalet 199-203, L'Hospitalet de Llobregat, Barcelona, Spain. ²²Unit of Nutrition and Cancer, Epidemiology, Public Health, Cancer Prevention and Palliative Program, Bellvitge Biomedical Research Institute (IDIBELL), L'Hospitalet de Llobregat, Av.
Granvia de L'Hospitalet, 199-203 Barcelona, Spain. ²³Escuela Andaluza de Salud Pública (EASP), 18011 Granada, Spain. ²⁴Instituto de Investigación Biosanitaria ibs.GRANADA, 18012 Granada, Spain. ²⁵Department of Preventive Medicine and Public Health, University of Granada, 18071 Granada, Spain. ²⁶Nutrition And Cancer Research Network (NACRE Network), Paris, France. ²⁷Department of Molecular Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, Nuthetal, Germany. ²⁸Institute of Nutritional Science, University of Potsdam, Nuthetal, Germany. ²⁹Instituto de Salud Pública y Laboral de Navarra, 31003 Pamplona, Spain. ³⁰Navarra Institute for Health Research (IdiSNA), 31008 Pamplona, Spain. ³¹Cancer Hospital and Research Centre of Santa Catarina (CEPON), Florianópolis, Brazil. ³²Public Health Policy Evaluation Unit, School of Public Health, Imperial College London, London, United Kingdom. ³³Department of Epidemiology and Biostatistics, School of Public Health, Imperial College London, London, United Kingdom. ³⁴Inserm U1018, "Exposome, Heredity, Cancer, and Health" Team, Hôpital Paul Brousse, Bâtiment 15/16, 16 avenue Paul Vaillant Couturier, 94807 VILLEJUIF CEDEX, France.

Received: 19 September 2023 Accepted: 25 April 2024
Published online: 09 October 2024

## References

Adams, J., & White, M. (2015). Characterisation of UK diets according to degree of food processing and associations with socio-demographics and obesity: cross-sectional analysis of UK National Diet and Nutrition Survey (2008-12). *International Journal of Behavioral Nutrition and Physical Activity, 12*(1), 160.

Bakker, M. F., Peeters, P. H., Klaasen, V. M., Bueno-de-Mesquita, H. B., Jansen, E. H., Ros, M. M., et al. (2016). Plasma carotenoids, vitamin C, tocopherols, and retinol and the risk of breast cancer in the European Prospective Investigation into Cancer and Nutrition cohort. *The American Journal of Clinical Nutrition, 103*(2), 454-64.

Bleiweiss-Sande, R., Chui, K., Evans, E. W., Goldberg, J., Amin, S., & Sacheck, J. (2019). Robustness of Food Processing Classification Systems. *Nutrients, 11*(6), 1344.

Chang, K., Gunter, M. J., Rauber, F., Levy, R. B., Huybrechts, I., Kliemann, N., et al. (2023). Ultra-processed food consumption, cancer risk and cancer mortality: a large-scale prospective analysis within the UK Biobank. *EClinicalMedicine, 56*, 101840.

Clarke R, Shipley M, Lewington S, Youngman L, Collins R, Marmot M, et al. (1999) Underestimation of risk associations due to regression dilution in long-term follow-up of prospective studies. *American Journal Of Epidemiology, 150*(4):341-53.

Couto, E., Boffetta, P., Lagiou, P., Ferrari, P., Buckland, G., Overvad, K., et al. (2011). Mediterranean dietary pattern and cancer risk in the EPIC cohort. *British Journal of Cancer, 104*(9), 1493-9.

Crino M, Barakat T, Trevena H, Neal B (2017) Systematic Review and Comparison of Classification Frameworks Describing the Degree of Food Processing. *Nutrition Food Technology, 3*(1).

Fiolet, T., Srour, B., Sellem, L., Kesse-Guyot, E., Allès, B., Méjean, C., et al. (2018). Consumption of ultra-processed foods and cancer risk: results from NutriNet-Santé prospective cohort. *BMJ: British medical journal / British Medical Association, 360*, k322.

Friedman, M. (2015). Acrylamide: inhibition of formation in processed food and mitigation of toxicity in cells, animals, and humans. *Food & Function, 6*(6), 1752-72.

Hutcheon, J. A., Chiolero, A., & Hanley, J. A. (2010). Random measurement error and regression dilution bias. *BMJ, 340*, c2289. British Medical Journal Publishing Group.

Huybrechts I, Rauber F, Nicolas G, Casagrande C, Kliemann N, Wedekind R, et al (2022) Characterization of the degree of food processing in the European Prospective Investigation into Cancer and Nutrition: Application of the Nova classification and validation using selected biomarkers of food processing. *Frontiers in Nutrition*.

International Statistical Classification of Diseases and Related Health Problems 10th Revision. https://icd.who.int/browse10/2010/en. Accessed 11 Nov 2023.

Jacobs I, Taljaard-Krugell C, Wicks M, Cubasch H, Joffe M, Laubscher R, et al (2022) Degree of food processing and breast cancer risk in black urban

women from Soweto, South African: the South African Breast Cancer study. *British Journal of Nutrition*: 1-12.

Juul, F., Parekh, N., Martinez-Steele, E., Monteiro, C. A., & Chang, V. W. (2022). Ultra-processed food consumption among US adults from 2001 to 2018. *The American Journal of Clinical Nutrition, 115*(1), 211-21.

Kliemann, N., Rauber, F., Bertazzi Levy, R., Viallon, V., Vamos, E. P., Cordova, R., et al. (2023). Food processing and cancer risk in Europe: results from the prospective EPIC cohort study. *Lancet Planet Health, 7*(3), e219-32.

Lane, M. M., Davis, J. A., Beattie, S., Gómez-Donoso, C., Loughman, A., O'Neil, A., et al. (2021). Ultraprocessed food and chronic noncommunicable diseases: a systematic review and meta-analysis of 43 observational studies. *Obesity Reviews, 22*(3), e13146.

Latasa, P., Louzada, M. L. D. C., Martinez Steele, E., & Monteiro, C. A. (2018). Added sugars and ultra-processed foods in Spanish households (1990-2010). *European Journal of Clinical Nutrition, 72*(10), 1404-12. Nature Publishing Group.

Levy RB, Andrade GC, Cruz GL da, Rauber F, Louzada ML da C, Claro RM, et al. (2022) Três décadas da disponibilidade domiciliar de alimentos segundo a NOVA - Brasil, 1987-2018. 56:75.

Luiten, C. M., Steenhuis, I. H., Eyles, H., Ni Mhurchu, C., & Waterlander, W. E. (2016). Ultra-processed foods have the worst nutrient profile, yet they are the most available packaged products in a sample of New Zealand supermarkets. *Public Health Nutrition, 19*(3), 530-8.

Lunn, M., & McNeil, D. (1995). Applying Cox regression to competing risks. *Biometrics, 51*(2), 524-32.

Madruga M, Martínez Steele E, Reynolds C, Levy RB, Rauber F (2022) Trends in food consumption according to the degree of food processing among the UK population over 11 years. *British Journal of Nutrition*: 1-24.

Marrón-Ponce, J. A., Flores, M., Cediel, G., Monteiro, C. A., & Batis, C. (2019). Associations between Consumption of Ultra-Processed Foods and Intake of Nutrients Related to Chronic Non-Communicable Diseases in Mexico. *Journal of the Academy of Nutrition and Dietetics, 119*(11), 1852-65.

Monteiro, C. A., Cannon, G., Levy, R. B., Moubarac, J.-C., Louzada, M. L., Rauber, F., et al. (2019). Ultra-processed foods: what they are and how to identify them. *Public Health Nutrition, 22*(5), 936-41.

Moubarac, J.-C., Martins, A. P. B., Claro, R. M., Levy, R. B., Cannon, G., & Monteiro, C. A. (2013). Consumption of ultra-processed foods and likely impact on human health. Evidence from Canada. *Public Health Nutrition, 16*(12), 2240-8.

Moubarac, J.-C., Batal, M., Martins, A. P. B., Claro, R., Levy, R. B., Cannon, G., et al. (2014). Processed and ultra-processed food products: consumption trends in Canada from 1938 to 2011. *Canadian Journal of Dietetic Practice and Research, 75*(1), 15-21.

Moubarac, J.-C., Parra, D. C., Cannon, G., & Monteiro, C. A. (2014). Food classification systems based on food processing: significance and implications for policies and actions: a systematic literature review and assessment. *Current Obesity Reports, 3*(2), 256-72.

Pagliai, G., Dinu, M., Madarena, M. P., Bonaccio, M., Iacoviello, L., & Sofi, F. (2021). Consumption of ultra-processed foods and health status: a systematic review and meta-analysis. *British Journal of Nutrition, 125*(3), 308-18.

Pouzou, J. G., Costard, S., & Zagmutt, F. J. (2018). Probabilistic assessment of dietary exposure to heterocyclic amines and polycyclic aromatic hydrocarbons from consumption of meats and breads in the United States. *Food and Chemical Toxicology, 114*, 361-74.

Queiroz SA, Sousa IM de, Silva FR de M, Lyra C de O, Fayh APT (2018) Nutritional and environmental risk factors for breast cancer: a case-control study. *Scientia Medica*. 28(2):ID28723-ID28723.

Rauber F, da Costa Louzada ML, Steele EM, Millett C, Monteiro CA, Levy RB. (2018). Ultra-Processed Food Consumption and Chronic Non-Communicable Diseases-Related Dietary Nutrient Profile in the UK (2008-2014). *Nutrients, 10*(5), 587

Riboli, E., Hunt, K. J., Slimani, N., Ferrari, P., Norat, T., Fahey, M., et al. (2002). European Prospective Investigation into Cancer and Nutrition (EPIC): study populations and data collection. *Public Health Nutrition, 5*(6B), 1113-24.

Romaguera, D., Fernández-Barrés, S., Gracia-Lavedán, E., Vendrell, E., Azpiri, M., Ruiz-Moreno, E., et al. (2021). Consumption of ultra-processed foods and drinks and colorectal, breast, and prostate cancer. *Clinical Nutrition, 40*(4), 1537-45.

Romieu I, Khandpur N, Katsikari A, Biessy C, Torres-Mejía G, Ángeles-Llerenas A, Alvarado-Cabrero I, Sánchez GI, Maldonado ME, Porras C, Rodriguez AC, Garmendia ML, Chajés V, Aglago EK, Porter PL, Lin M, His M, Gunter MJ,
Huybrechts I, Rinaldi S; PRECAMA team (2022) Consumption of industrial processed foods and risk of premenopausal breast cancer among Latin American women: the PRECAMA study. *BMJ Nutr Prev Health, 5*(1), 1–9

Srour, B., Kordahi, M. C., Bonazzi, E., Deschasaux-Tanguy, M., Touvier, M., & Chassaing, B. (2022). Ultra-processed foods and human health: from epidemiological evidence to mechanistic insights. *The Lancet Gastroenterology & Hepatology, 5*(2), 12468–1253(22), 00169–8.

Wang, L., Martínez Steele, E., Du, M., Pomeranz, J. L., O'Connor, L. E., Herrick, K. A., et al. (2021). Trends in Consumption of Ultraprocessed Foods Among US Youths Aged 2–19 Years, 1999–2018. *JAMA, 326*(6), 519–30.

Wareham, N. J., Jakes, R. W., Rennie, K. L., Schuit, J., Mitchell, J., Hennings, S., et al. (2003). Validity and repeatability of a simple index derived from the short physical activity questionnaire used in the European Prospective Investigation into Cancer and Nutrition (EPIC) study. *Public Health Nutrition, 6*(4), 407–13.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.