Research Abstract Objective To jointly analyse two food dimensions, the Food Standards Agency Nutrient Profiling System (FSAm-NPS), used to derive the Nutri-Score front-of-pack label, and the NOVA classification in relation to mortality. Design Prospective cohort study. Setting Moli-sani Study, Italy 2005-10. Participants 22 895 participants (mean age 55 (SD […]
Click here to view original web page at www.bmj.com
Research
Abstract
Objective To jointly analyse two food dimensions, the Food Standards Agency Nutrient Profiling System (FSAm-NPS), used to derive the Nutri-Score front-of-pack label, and the NOVA classification in relation to mortality.
Design Prospective cohort study.
Setting Moli-sani Study, Italy 2005-10.
Participants 22 895 participants (mean age 55 (SD 12) years; 48% men).
Main outcomes measures Associations between dietary exposures and mortality risk, assessed using multivariable cause specific Cox proportional hazard models controlled for known risk factors.
Results A total of 2205 deaths occurred during 272 960 person years of follow-up. In the highest quarter of the FSAm-NPS index compared with the lowest quarter, multivariable adjusted hazard ratios for all cause and cardiovascular mortality were 1.19 (95% confidence interval 1.04 to 1.35; absolute risk difference 4.3%, 95% confidence interval 1.4% to 7.2%) and 1.32 (1.06 to 1.64; 2.6%, 0.3% to 4.9%), respectively. The hazard ratios were 1.19 (1.05 to 1.36; absolute risk difference 9.7%, 5.0% to 14.3%) and 1.27 (1.02 to 1.58; 5.0%, 1.2% to 8.8%), respectively, for all cause and cardiovascular mortality when the two extreme categories of ultra-processed food intake were compared. When these two indices were analysed jointly, the magnitude of the association of the FSAm-NPS dietary index with all cause and cardiovascular mortality was attenuated by 22.3% and 15.4%, respectively, whereas mortality risks associated with high ultra-processed food intake were not altered.
Conclusions Adults with the lowest quality diet, as measured using the FSAm-NPS dietary index (underpinning the Nutri-Score), and the highest ultra-processed food consumption (NOVA classification) were at the highest risk for all cause and cardiovascular mortality. A significant proportion of the higher mortality risk associated with an elevated intake of nutrient poor foods was explained by a high degree of food processing. In contrast, the relation between a high ultra-processed food intake and mortality was not explained by the poor quality of these foods.
Introduction
Poor diets are responsible for more deaths than any other risk factor globally and are the leading cause of obesity and non-communicable diseases.12 Traditional approaches to improve nutrition and health have focused on nutrients as the key determinants of a healthful diet. This nutrient based approach is the basis of the vast majority of food based dietary guidelines and educational programmes worldwide and emphasises the consumption of foods that are natural sources of fibre, vitamins, and minerals and low in saturated fat, cholesterol, and sodium.3 By contrast, the NOVA classification has been proposed as a novel way to rate foods on the basis of the degree of processing rather than their nutrient content; processing of foods might play a role in health beyond their nutritional composition, through a variety of mechanisms triggered by non-nutritional components, such as cosmetic additives, food contact materials, neoformed compounds, and degradation of the food matrix.45 Despite large and well conducted cohort studies indicating that increased amounts of ultra-processed food in the food supply are linked to a higher risk of non-communicable diseases, including obesity and mortality,678 the role of food processing is still largely ignored or minimised.59
Governmental strategies to improve diet quality for reducing the burden of non-communicable diseases include the implementation of labelling systems for printing on the front of packages of food products. A variety of front-of-pack labelling systems exist, including both nutrient specific warning labels (for example, to alert consumers when a given food contains an excessive amount of nutrients that should have lower consumption in the diet) and summary indicators that rate the overall nutritional quality of a product, based on algorithms that take both positive and negative attributes of a food’s nutritional composition into account.1011 Whatever the labelling system used, the main aim of front-of-pack labels is to inform consumers of the overall nutritional value of foods to help them to make healthier dietary choices and, in some cases, to incentivise food companies to reformulate the nutritional composition of their products.1011 Front-of-pack labels are agreed to be useful tools to improve people’s diet quality, as their use was found to encourage healthier food purchasing.1213
In May 2020 the European Commission announced the intended adoption of a harmonised and mandatory front-of-pack nutrition labelling scheme at the EU level by the end of 2022, as part of its Farm-to Fork Strategy.1014 The Nutri-Score,15 developed by French academic researchers and endorsed in seven European countries, is possibly the most well studied nutritional labelling system; several investigations have supported its effectiveness in improving diet quality and associating with health outcomes.1016171819 The Nutri-Score is a five colour coded scheme ranging from dark green to dark orange, associated with letters from A to E, to optimise logo accessibility and understanding by consumers of food products; its main goal is to provide the consumer with the information needed to compare the nutritional quality of products in the same category or consumed for the same purpose.10
Although they have gained much attention from the EU authorities, concerns exist about the effectiveness of the front-of-pack nutrition labelling schemes in improving people’s diet quality. Some of the arguments include the fact that the Nutri-Score system, and theoretically all systems focusing only on nutritional composition, does not cover all the dimensions of the food, which cannot be limited to its nutrient balance.2021 Some foods can be considered nutritionally adequate although being highly processed and containing many additives, as well as other substances that are not typically used in domestic culinary recipes (for example, modified starches, hydrogenated oils, syrups obtained from cereals, flavours).45 Moreover, a large scale meta-analysis of nationally representative samples recently found that the classification by degree of food processing correlates with nutritional quality, further complicating the overall understanding of the effect of these two dimensions on human health.22
To contribute to a better understanding of which dimension of foods plays a major role in defining health risk at the population level, we analysed the individual and joint associations of two food classification systems in relation to all cause and cause specific mortality in a large sample of Italian adults from the Moli-sani Study cohort. As a secondary aim, we analysed some biological mechanisms potentially linking these two dietary indicators to the health outcomes investigated. For the purpose of this study, we analysed the Food Standards Agency Nutrient Profiling System (FSAm-NPS), used to compute the Nutri-Score ranking foods according to their nutritional value,23 and the NOVA classification,4 which evaluates the degree of food processing.
Methods
We analysed data from the Moli-sani Study, a population based cohort study established in 2005-10 with an enrolment of 24 325 men and women (aged ≥35 years) randomly recruited from the general population of Molise, a southern Italian region. The main purpose was to investigate genetic and environmental risk factors in the onset and progression of cardiovascular, cerebrovascular, and cancer diseases. Exclusion criteria were pregnancy at the time of recruitment, disturbances in mental capacity or decision making impairments, current poly-traumas or coma, or refusal to give informed consent. Details of the study are available elsewhere.24
For the purpose of this study, we excluded participants with missing data on diet, implausible energy intakes (<800 or >4000 kcal/d in men and <500 or >3500 kcal/d in women), incomplete dietary or medical questionnaires, or missing data on cause specific death. We finally analysed 22 895 participants. Supplementary figure A shows the flowchart for selection of study participants.
Dietary assessment and calculation of FSAm-NPS dietary index
We assessed food intake during the year before enrolment by an interviewer administered semiquantitative European Prospective Investigation into Cancer and Nutrition (EPIC) food frequency questionnaire (FFQ) validated and adapted to the Italian population.25 The FFQ contains 14 sections (pasta/rice, soup, meat (excluding salami and other cured meats), fish, raw vegetables, cooked vegetables, eggs, sandwiches, salami and other cured meats, cheese, fruit, bread/wine, milk/coffee/cakes, and herbs/spices), with 248 questions concerning 188 different food items. We asked participants to indicate the number of times they consumed a certain food item (per day, week, month, or year), from which we calculated the absolute frequency of consumption of each item. We assessed the quantity of the food consumed by asking the participant to select an image of a food portion or a predefined standard portion when no image was available. We used 17 sets of pictures, each showing a small, medium, and large portion size, with additional quantifiers (for example, “smaller than the small portion” or “between the small and medium portion”).26
We linked frequencies and quantities of each food to Italian food tables by using specifically designed software to obtain estimates of daily intake of macro-nutrients and micro-nutrients plus energy.26 We evaluated adherence to the traditional Mediterranean diet by using the Mediterranean Diet Score developed by Trichopoulou and colleagues.27
The FSAm-NPS is a modified version of the nutrient profiling system (FSA-NPS), which was initially developed in the UK to regulate television food advertising to children.1015 We calculated the FSAm-NPS score as was previously done in other population cohorts.181923 For all foods and beverages included in the FFQ, on the basis of composition for each 100 g of content, we allocated 0 to 40 points for nutrients that should be consumed in limited amounts (A points; that is, total sugars (g), saturated fats (g), sodium (mg), and energy (kJ)) and 0 to 15 points for nutrients or components that should be promoted (C points; that is, dietary fibre (g), protein (g), and fruit, vegetables, legumes, and nuts (%)). We calculated the total score of the product by subtracting the sum of C points from the sum of A points. Thus, the final FSAm-NPS score for each food/beverage was based on a scale that could theoretically range from −15 (healthiest food) to 40 (least healthy food). On the basis of this overall FSAm-NPS score, the Nutri-Score labelling system categorises food products into five colours, each associated to letters from A (dark green) to E (dark orange) reflecting their nutritional quality (supplementary table A).15
We calculated the FSAm-NPS dietary index at the individual level as an energy weighted mean of the FSAm-NPS scores of all foods and beverages consumed by each participant by using the following equation, in which FSi represents the score of food/beverage i, Ei the energy intake from food/beverage i specific for each participant, and n the total number of foods/beverages consumed. Increasing values of the FSAm-NPS dietary index therefore reflect decreasing diet quality overall.
NOVA classification
To estimate ultra-processed food intake, we used the NOVA classification that groups foods into four categories reflecting different levels of processing: (a) fresh or minimally processed foods (for example, fruit, meat, milk)—that is, foods altered only by processes such as removal of inedible or unwanted parts, drying, crushing, grinding, fractioning, roasting, boiling, pasteurisation, refrigeration, freezing, placing in containers, vacuum packaging, or non-alcoholic fermentation without the addition of salt, sugar, oils or fats, or other food substances to the original food; (b) processed culinary ingredients (for example, oils, butter, sugar); (c) processed foods (for example, canned fish, unpackaged freshly made breads) that are manufactured by adding salt, sugar, oil, or other processed culinary ingredients to minimally processed foods; or (d) ultra-processed foods defined as industrial formulations made mostly or entirely from substances extracted from foods or derived from food constituents often containing added flavours, colours, emulsifiers, and other cosmetic additives and little or no whole foods (for example, carbonated drinks, processed meat, sweet or savoury packaged snacks).428 For the purpose of these analyses, we used the fourth ultra-processed foods category. We summed up the amount consumed (g/d) of each food group from the fourth category of NOVA (a total of 22 foods and beverages) and calculated the proportion (%) of ultra-processed food in the total weight of food and beverages consumed (g/d) by creating a weight ratio. Such an approach is more appropriate than an energy ratio, as it better accounts for non-nutritional factors pertaining to food processing (for example, neo-formed contaminants, additives, and alterations to the structure of raw foods).29
We then divided participants into quarters based on the proportion of ultra-processed food consumed over the total food intake. The full list of foods categorised according to the NOVA classification is available as supplementary table B.
Follow-up for vital status
We followed up the Moli-sani Study cohort for mortality from March 2005 to 31 December 2019. We assessed cause specific mortality by using the Italian mortality registry, validated by Italian death certificates (ISTAT form) and coded according to ICD-9 (international classification of diseases, revision 9). Cardiovascular mortality included deaths from diseases of the circulatory system, when the underlying cause of death included ICD-9 codes 390-459. We used ICD-9 codes 430-438 to define specific cause of death for cerebrovascular disease and ICD-9 codes 410-414 and 429 for ischaemic heart disease. Cancer death was when the underlying cause of death included ICD-9 codes 140-208. We included non-cardiovascular/non-cancer causes of death in an “other cause mortality” group.
Baseline covariates assessment
Personal history of cardiovascular disease (angina, myocardial infarction, revascularisation procedures, peripheral artery diseases, and cerebrovascular events) was self-reported and confirmed by medical records and therapy. Personal history of cancer was self-reported and confirmed by medical records. We considered participants to have diabetes, hypertension, or hyperlipidaemia at baseline if they were taking disease specific drugs.
We expressed leisure time physical activity as daily energy expenditure in metabolic equivalent task hours (MET-h/d) for sport, walking, and gardening. Height and weight were measured, and body mass index was calculated as kg/m2. Blood pressure was measured by an automatic device (OMRON-HEM-705CP) three times on the non-dominant arm, with the average of the last two values being taken as the blood pressure. Measurements were made in a quiet room with comfortable temperature, with the participants lying down for at least five minutes.
We classified participants as never, current, or former smokers (reported not having smoked at all over the previous 12 months or more). Educational level was based on the highest qualification attained and categorised as up to lower secondary (approximately ≤8 years of study), upper secondary school (>8 and ≤13 years), and postsecondary education (>13 years). We classified housing tenure as rented, ownership of one dwelling, and ownership of more than one dwelling.
Selection of cardiovascular risk factors
Key biological mechanisms through which nutrient poor diets and ultra-processed foods may adversely affect health include, among others, altered serum lipid concentrations, inflammation, oxidative stress, dysglycaemia, insulin resistance, and hypertension.30 We selected biomarkers reflecting different underlying pathways to incidence and progression of cardiovascular disease,3132 as potential mediators of an association between diet and mortality, by subject area knowledge according to the following criteria: previously studied for their relevance in pathways predisposing to cardiovascular disease, shown in epidemiological studies to be related to cardiovascular disease or mortality, and already investigated in the Moli-sani Study cohort. Assessment of blood biomarkers of cardiovascular disease risk in the Moli-sani Study cohort is described in the supplementary materials.
Statistical analysis
We report baseline characteristics of the analytical sample across quarters of the FSAm-NPS dietary index (sex specific) and ultra-processed food intake as means (standard deviations) or percentages for continuous and categorical traits, respectively. We calculated differences in the distribution of baseline covariates across quarters of the dietary exposure by using generalised linear models adjusted for age, sex, and energy intake (GENMOD procedure for categorical variables and GLM procedure for continuous variables in SAS software) (table 1 and table 2).
We examined associations of FSAm-NPS dietary index and ultra-processed food intake (quarters and per 1 standard deviation increment) with all cause and cause specific mortality through multivariable cause specific Cox proportional hazards with time on study as the timescale and adjusting for baseline age as a covariate in the model. We visually assessed the proportional hazards assumption (log(−log) plots of survival curves) and identified no violation. We calculated multivariable adjusted hazard ratios across quarters of both dietary exposures (quarter 1 as reference) and included sex, age (continuous), energy intake (continuous), educational level (up to lower secondary, upper secondary, post-secondary), housing tenure (rented, one dwelling ownership, more than one dwelling ownership), smoking (never, current, former smokers), body mass index (continuous), leisure time physical activity (continuous), history of cancer (no/yes), history of cardiovascular disease (no/yes), diabetes (no/yes), hypertension (no/yes), hyperlipidaemia (no/yes), and residence (urban, rural). Participants contributed person time until their date of death, date of emigration, or loss to follow-up or until end of follow-up, whichever occurred first. We included participants who died from a cause other than the one under study and censored them at the date of the competing death event. We defined potential confounders a priori and identified them on the basis of existing literature, rather than deferring to statistical criteria.33
We calculated the absolute risk difference for hazard ratios derived from the previously described models.34 We used the mean value of the included covariates and calculated the difference in absolute risk at the maximum follow-up time. We applied the bootstrap method to derive 95% confidence intervals for absolute risk difference based on 500 bootstrap samples. To maximise data availability, we handled missing data on covariates (see flowchart in supplementary figure A) by using multiple imputation (SAS PROC MI, followed by PROC MIANALYZE; n=10 imputed datasets).
We considered a biomarker as potentially mediating the association of the FSAm-NPS dietary index or ultra-processed food intake with all cause and cause specific mortality if it was on the causal pathway of these associations and if it was associated with both the exposure and the outcome, in accordance with predefined mediation principles.35 We tested these criteria in distinct multivariable regression models for each potential mediator individually (supplementary tables C and D) and through cause specific Cox models including ultra-processed food consumption (continuous) or the FSAm-NPS dietary index (continuous) as a covariate (supplementary tables E and F).
To quantify how much of the association between the FSAm-NPS dietary index (or ultra-processed food intake) and health outcomes was explained by ultra-processed food intake (or by the FSAm-NPS dietary index), we relied on a traditional change-in-estimate method—that is, the quantification of the percentage reduction in the β coefficient for the FSAm-NPS dietary index (or ultra-processed food intake; β0) after inclusion of ultra-processed food intake (or the FSAm-NPS dietary index; β1) to the multivariable model. We calculated the attenuation (%) according to the following equation: 100×(β0–β1)/(β0), where β0=natural log (HR0) and β1=natural log (HR1). HR0 is the multivariable hazard ratio for one dietary exposure (for example, FSAm-NPS dietary index) not adjusted for the other (for example, ultra-processed food intake), and HR1 is the multivariable hazard ratio for the same dietary exposure resulting from the multivariable model further including the other dietary exposure. A negative attenuation indicates no mediation effect (that is, HR1>HR0). We calculated a 95% confidence interval around the percentage attenuation by using a bootstrap method with 1000 re-samplings for each imputed dataset.
The multivariable model also served as the reference for the analysis used to estimate the extent to which selected cardiovascular risk factors explained the association of the FSAm-NPS dietary index or ultra-processed food intake with all cause and cause specific mortality; for this, each marker was alternately, and at the end simultaneously, included into the multivariable adjusted model. For the mediation analysis, we used the publicly available %MEDIATE macro in SAS,36 which calculates the point and interval estimates of the percentage of exposure effect explained by one or more intermediate variables, with 95% confidence intervals and P values.
To test the robustness of the associations, we did sensitivity analyses by excluding participants with a history of cancer, cardiovascular disease, and diabetes; excluding participants with baseline diabetes, hypertension, and hyperlipidaemia (therefore assessing a potential bias resulting from modified habitual dietary intakes due to illness, such as indications to switch to a healthier diet); and using ultra-processed foods as an energy ratio (percentage of calories from ultra-processed foods on the total calories consumed daily), both as exposure and explanatory factors. We used SAS/STAT software, version 9.4, for data analysis.
Patient and public involvement
No patients were involved in setting the research question or the outcome measures; nor were they involved in developing plans for recruitment, design, or implementation of the study. No patients were asked to advise on the interpretation or writing up of results. Experts in communication and scientific dissemination were originally involved in the design of the study in order to develop effective strategies for public dissemination. All participants are periodically informed about the advancements of the study via an annual calendar, which contains information on the research activity related to the study, along with newsletters that are periodically sent to participants or published on the website of the Moli-sani Study (http://www.moli-sani.org/). Regular meetings are also organised with participants from each of the 30 towns/villages from the Molise region, randomly recruited in the cohort, to share the results of the study.
Results
The analytical sample consisted of 11 973 (52.3%) women and 10 922 (47.7%) men, with a mean age at enrolment of 55.4 (SD 11.7) years, a mean FSAm-NPS dietary index of 7.4 (1.4), and a mean ultra-processed food weight ratio of 10.8% (6.7%); the average energy from ultra-processed foods was 18.3% (8.4%; range 8.4-70.0%) of the total calories consumed daily. Although the means of both dietary indices are broadly the same (although in different units), the standard deviations are very different, with the standard deviation for ultra-processed foods being more than four times larger than that for the FSAm-NPS dietary index. As a result, the means for the four quarters are much more different in table 2 than in table 1. The correlation of FSAm-NPS dietary index with ultra-processed food intake was low to moderate (Spearman correlation coefficient=0.34).
Compared with participants in the bottom quarter, those scoring higher on the FSAm-NPS dietary index (quarter 4) were younger and more educated, lived prevalently in urban areas, had lower body mass index, practised less physical activity, and had lower prevalence of chronic conditions at baseline, with the exception of cancer (table 1). Findings were similar for the comparison of the lowest and highest quarters of ultra-processed food intake, with some differences including sex (men tended to consume less ultra-processed food than women), smoking status (current smokers were more represented in the highest quarter of ultra-processed food intake than the lowest), and distribution of baseline cardiovascular disease (observed across FSAm-NPS dietary index quarters but not across ultra-processed food intake quarters) (table 2).
Regarding nutritional factors, the two scores shared many similarities. Higher levels of FSAm-NPS dietary index or ultra-processed food intake were inversely associated with adherence to a Mediterranean diet, monounsaturated-to-saturated fat ratio, and consumption of fruits and nuts, vegetables, cereals, legumes, fish, alcohol, starch, fibre, and protein. Increases in both scores positively correlated with energy intake, total fat, saturated fat, polyunsaturated fat, and dietary cholesterol (table 3).
We observed divergent association for meats, which were inversely associated with ultra-processed food intake but not with the FSAm-NPS dietary index, and for milk and dairy products, which were directly correlated only with ultra-processed food intake. Energy from carbohydrate decreased according to the FSAm-NPS dietary index but increased with ultra-processed food intake, and sugar was directly associated with ultra-processed food intake but not with the FSAm-NPS dietary index. Monounsaturated fatty acids increased with ultra-processed food intake, and sodium was positively associated with the FSAm-NPS dietary index and inversely with ultra-processed food intake (table 3).
An increase in the FSAm-NPS dietary index was directly associated with higher concentrations of C reactive protein, white blood cell count, granulocyte-to-lymphocyte ratio, insulin, C-peptide, apolipoprotein B100, cystatin C, and diastolic blood pressure, while being inversely linked to blood glucose concentrations, high density lipoprotein cholesterol, triglycerides, lipoprotein (a), and serum vitamin D (supplementary table G). A higher proportion of ultra-processed foods in the diet positively correlated with granulocyte-to-lymphocyte ratio, insulin, cystatin C, creatinine, and heart rate and inversely with blood glucose, blood cholesterol, high density lipoprotein cholesterol, triglycerides, apolipoprotein A, apolipoprotein B100, lipoprotein (a), and systolic and diastolic blood pressure (supplementary table G). Most of these associations withstood multivariable adjusted regression analyses, with some exceptions, such as total blood cholesterol, which became inversely associated with the FSAm-NPS dietary index (supplementary tables C and D).
Of 10 risk factors (insulin, C-peptide, total blood cholesterol, triglycerides, lipoprotein (a), C reactive protein, white blood cell count, granulocyte-to-lymphocyte ratio, diastolic blood pressure, and serum vitamin D) associated with the FSAm-NPS dietary index in cross sectional analysis, eight (insulin, C-peptide, total blood cholesterol, C reactive protein, white blood cell count, granulocyte-to-lymphocyte ratio, diastolic blood pressure, and serum vitamin D) and five (C reactive protein, white blood cell count, granulocyte-to-lymphocyte ratio, diastolic blood pressure, and serum vitamin D) were also associated with all cause and cardiovascular disease mortality, respectively, in a multivariable model also adjusted for the FSAm-NPS dietary index (supplementary table E); we therefore included them in the mediation analysis. For 17 of 18 risk factors associated with ultra-processed food intake (all markers with the exception of serum vitamin D), 12 were associated with all cause mortality (cystatin C, creatinine, blood glucose, insulin, C-peptide, total blood cholesterol, C reactive protein, white blood cell count, granulocyte-to-lymphocyte ratio, diastolic and systolic blood pressure, and heart rate) and eight with cardiovascular disease mortality (cystatin C, creatinine, C reactive protein, white blood cell count, granulocyte-to-lymphocyte ratio, diastolic and systolic blood pressure, and heart rate) in a multivariable adjusted model also including ultra-processed food intake (supplementary table F).
Association with mortality
Among 2205 deaths from any cause that occurred over a median follow-up of 12.2 (interquartile range 11.2-13.2) years (272 960 person years), 792 were attributed to cardiovascular disease, of which 426 were due to ischaemic heart disease/cerebrovascular disease, 820 to cancer, and 593 to other causes. In the multivariable adjusted model controlled for sociodemographic and clinical factors, the hazard ratios associated with the highest quarter of the FSAm-NPS dietary index were 1.19 (95% confidence interval 1.04 to 1.35; absolute risk difference 4.3%, 95% confidence interval 1.4% to 7.2%) for all cause mortality and 1.32 (1.06 to 1.64; 2.6%, 0.3% to 4.9%) for cardiovascular disease mortality (fig 1; supplementary table H). We found no associations with the other causes of death (fig 1).
Fig 1
All cause and cause specific mortality estimates across quarters of Food Standards Agency nutrient profiling system (FSAm-NPS) dietary index in Moli-sani Study cohort (n=22 895), using data obtained from multiple imputation. Hazard ratios with 95% confidence intervals obtained from multivariable cause specific Cox proportional hazards regression models including sex, age (continuous), energy intake (continuous), educational level (up to lower secondary, upper secondary, post-secondary), housing tenure (rent, 1 dwelling ownership, >1 dwelling ownership), smoking (never, current, former smokers), body mass index (continuous), leisure time physical activity (continuous), history of cancer (no/yes), history of cardiovascular disease (no/yes), diabetes (no/yes), hypertension (no/yes), hyperlipidaemia (no/yes), and residence (urban, rural). IHD=ischaemic heart disease
The inclusion of ultra-processed food intake (weight ratio; continuous) into the multivariable model attenuated the association of the FSAm-NPS dietary index with all cause mortality by 22.3% (95% confidence interval 16.4% to 30.2%) and mitigated that with cardiovascular disease mortality by 15.4% (10.5% to 22.6%) (table 4). Consistently, absolute risk differences were also reduced (supplementary table H).
Compared with the lowest quarter, an elevated intake of ultra-processed foods was associated with a higher hazard of all cause mortality (hazard ratio 1.19, 1.05 to 1.36 for quarter 4 versus quarter 1; absolute risk difference 9.7%, 5.0% to 14.3%), cardiovascular disease mortality (1.27, 1.02 to 1.58; 5.0%, 1.2% to 8.8%) (fig 2; supplementary table H), ischaemic heart disease/cerebrovascular disease mortality (1.39, 1.03 to 1.88; 3.0%, 0.1% to 5.9%), and other cause mortality (1.37, 1.07 to 1.70; 7.3%, −2.7% to 17.4%) (fig 2). By contrast with the FSAm-NPS dietary index, for which excess risk was restricted to the highest quarter, these hazard ratios increased close to monotonically across quarters.
Fig 2
All cause and cause specific mortality estimates across quarters of ultra-processed food intake (weight ratio) as defined by NOVA classification in Moli-sani Study cohort (n=22 895), using data obtained from multiple imputation. Hazard ratios with 95% confidence intervals obtained from multivariable cause specific Cox proportional hazards regression models including sex, age (continuous), energy intake (continuous), educational level (up to lower secondary, upper secondary, post-secondary), housing tenure (rent, 1 dwelling ownership, >1 dwelling ownership), smoking (never, current, former smokers), body mass index (continuous), leisure time physical activity (continuous), history of cancer (no/yes), history of cardiovascular disease (no/yes), diabetes (no/yes), hypertension (no/yes), hyperlipidaemia (no/yes), and residence (urban, rural). IHD=ischaemic heart disease
The inclusion of the FSAm-NPS dietary index (continuous) into the multivariable model did not substantially alter the strength of associations between ultra-processed food intake and mortality risks either as hazard ratios (table 4) or as absolute risk differences (supplementary table H). Analyses using ultra-processed food intake as an energy ratio, both as an exposure and as an explanatory factor, yielded similar results (supplementary table I).
Results did not substantially change in sensitivity analyses, showing that the excess risk of all cause and cardiovascular disease mortality associated with a higher FSAm-NPS dietary index was explained, at least in part, by an elevated degree of food processing, whereas this was not the case for the ultra-processed food-mortality associations (supplementary tables J and K).
Analyses of biological pathways
The excess risk of all cause mortality associated with higher FSAm-NPS dietary index was partly accounted for by altered markers of glucose metabolism and serum cholesterol, which explained 8.9% and 4.2% of this association; all markers explained 18.6% (P<0.001) (table 5). For ultra-processed food intake, the excess of all cause mortality in the highest quarter was largely explained by altered levels of biomarkers reflecting renal function (26.0%; P<0.001) followed by inflammatory markers (10.0%; P=0.007), total blood cholesterol (8.5%; P<0.001), and markers of glucose metabolism (5.8%; P=0.02). Altogether, these factors explained up to 32.7% (P<0.001) of the association of heavy ultra-processed food intake with all cause mortality. The association of high ultra-processed food consumption with cardiovascular disease mortality was mediated by biomarkers of renal function and inflammation (table 5).
Discussion
In a large prospective cohort of 22 895 Italian adults, both diets mainly composed of food products with higher FSAm-NPS, which reflects poor nutritional profiles, and those with a large dietary share of ultra-processed foods were associated with higher hazards of all cause and cardiovascular disease mortality. Increased ultra-processed food consumption, but not the FSAm-NPS dietary index, was also an independent risk factor for mortality due to ischaemic heart disease/cerebrovascular disease and other causes. When these two food dimensions (that is, nutrient balance and food processing) were analysed simultaneously, we observed that the associations of the FSAm-NPS dietary index with all cause and cardiovascular disease mortality were significantly attenuated, whereas estimations for ultra-processed food intake remained almost unchanged for all the outcomes under study. These findings suggest that highly processed foods are associated with poor health outcomes independently of their low nutritional composition, but not the other way around.
Comparison with other studies
This is the first study providing a targeted analysis of how nutritional and non-nutritional food dimensions are associated with mortality risk. Although no previous study has specifically investigated this, other investigations have indirectly considered both nutritional quality and food processing in association with a variety of health outcomes. A review of the literature analysing the relative effect of adjusting for diet quality/patterns on the reported associations between ultra-processed food intake and health related outcomes in more than 20 prospective cohort studies found that these adjustments did not explain the association between ultra-processed food intake and health related outcomes, with estimates remaining highly significant.37 On the contrary, none of the large cohort studies evaluating the health impact of the FSAm-NPS dietary index considered the degree of food processing as a potential covariate,1838 thereby introducing potential sources of bias in these studies.
One exception is the longitudinal analysis from the SUN cohort of Spanish graduates showing that the association between a higher FSAm-NPS dietary index and all cause mortality was stronger among participants with a high intake of ultra-processed foods.19 This may therefore be an underrated common factor in most of the large scale prospective cohort studies that have evaluated the impact of the Nutri-Score on health; much of the research in this field did not consider that food processing and nutritional quality are partially correlated, making disentangling their independent effects on human health difficult. Interestingly, a study comparing the nutritional quality (as assessed by the Nutri-Score) and the degree of processing (as assessed by the NOVA classification) of foods in the Open Food Facts database reported that ultra-processed foods were represented to different extents across all Nutri-Score categories, ranging from 26.1% in nutritional category A (highest nutritional quality) up to 83.7% in nutritional category E (lowest nutritional quality).20 As a result, the authors suggested that front-of-pack labels could benefit from more details besides the Nutri-Score and information on the nutrient content, including a warning on the level of food processing, to really improve people’s diets.20 This position is also endorsed by other experts who call for incorporation of a warning label for ultra-processed foods as an additional measure to guide people towards healthier eating.39
In our study, ultra-processed food consumption was more evidently associated with biomarkers potentially reflecting the biological mechanisms behind its possible effect on mortality risk than the FSAm-NPS dietary index. Moreover, we observed that, although a nutritionally poor diet, as reflected by consumption of foods with higher FSAm-NPS scores, may affect mortality risk through an unfavourable modulation of pathways that are known to be affected by diet (for example, markers of glucose metabolism), the higher mortality hazards associated with a diet rich in ultra-processed foods were ascribed to altered renal function and to a lesser extent to increased markers of inflammation.
Evidence on the potential effect of the overall diet quality on biomarkers of renal function, a well established risk factor for cardiovascular disease, is not robust, with only a few cohort studies supporting an association.4041 However, diets high in ultra-processed foods have been increasingly shown to be associated with altered renal function and higher inflammation,7424344 possibly through mechanisms that are triggered by non-nutritional components of the diet, such as food additives and contaminants present in highly processed foods, as well as food processing itself that affects both nutritional composition and the food matrix.45 For instance, the packaging of ultra-processed foods is a major source of synthetic chemicals, such as phthalates and bisphenols, that are among the so-called endocrine disrupting chemicals,46 which could have adverse effects on renal function and might contribute to progressive cumulative renal injury over a lifetime and are also associated with altered concentrations of inflammatory biomarkers.474849 Consistently, evidence shows that acrylamide, which is one of the most relevant contaminants produced when foods are heated to high temperatures (for example, French fries and potato crisps, cereal products, and roasted coffee), has documented nephrotoxic effects and has been linked to increased oxidative stress and inflammation,5051 as well as some food additives largely used in the food industry.52
Although access to edible, safe, and healthy food is essential, the usefulness of high levels of food processing has been strongly questioned,21 and it has several implications for human health. Modifications to the food matrix during processing can alter nutrient bio-accessibility and absorption kinetics, which may promote an inflammatory gut microbiota that in turn is associated with several cardiometabolic conditions.30 Food processing can also lead to the loss of some protective micronutrients and phytochemicals naturally present in plant foods.21 The differential effect of nutritional quality and food processing on biological pathways represent an additional valuable reason to consider both these two food dimensions as equally important for human health.
Strengths and limitations of study
This is the first study examining the possible health impact of the joint exposure to nutrient poor foods, according to the nutrient profile system underpinning the Nutri-Score front-of-pack label, and to foods characterised by an elevated degree of processing as described by the NOVA classification. Strengths of this analysis include the prospective design, long follow-up, and use of a large dataset, with a careful account for a large number of covariates to minimise confounding. However, the study also has several limitations. Firstly, owing to the observational design, we cannot fully rule out the potential role of residual confounding by unmeasured factors. Dietary data were self-reported, and this may lead to recall bias. Also, potential exists for social desirability bias that might lead to the underreporting of ultra-processed food consumption, which could bias studied associations towards the null. Moreover, the FFQ used in this study, like most FFQs used in large scale prospective cohort studies,1819 was not originally developed to assess the degree of food processing, so many food items were not included (for example, pre-prepared dishes, energy bars, slimming products).
Although no one gold standard for applying the NOVA categorisation exists,53 we recognise that FFQs may not cover the full spectrum of foods consumed, including ultra-processed foods, owing to the limited number of predefined food lists and the lack of supporting information on cooking methods, ingredients, eating place, and the brand names of the packaged foods,5354 which would be extremely useful in identifying ultra-processed foods. However, most existing large cohort studies conducted to explore the relation between ultra-processed food consumption and mortality risks have so far used FFQs.555657 Moreover, imprecision in the identification of ultra-processed foods may also pertain to 24 hour diet recalls or diet records.53 Finally, evidence suggests that classification of ultra-processed foods with an FFQ may be valid for the purpose of qualitative comparisons, although less appropriate in absolute intake estimations.54 However, we expect misclassification to be non-systematic, and this would likely lead to non-differential measurement error possibly resulting in an underestimation of the studied associations.
Another weakness is that diet and health data were measured at baseline only, so potential changes occurring over the life course might have modified the strength of the findings; nevertheless, evidence shows that diet in adulthood tends to remain stable over time,58 and most of the biomarkers we tested were not found to vary substantially over time.59 Finally, caution is needed in generalising these findings to other populations.
Conclusions and policy implications
The diet-health relation has been traditionally explored and explained almost exclusively by nutritional composition, leading to recommendations of limiting nutrients to be avoided or reduced (for example, sugar, salt, and fat) while favouring others (for example, fibre) to prevent major nutrition related non-communicable diseases.21 Findings from this large cohort of an Italian general population suggest that part of the mortality risk associated with a nutritionally unbalanced diet is due to an elevated degree of processing that usually characterises nutrient poor foods. Our findings suggest that food processing and the nutritional quality of food cover different but complementary dimensions, which should both be considered when analysing the diet-disease relation. This is in line with the assumption that a food’s health potential is not exclusively associated with its nutritional composition.21
From a public health perspective, this study reinforces the opportunity to reformulate dietary guidelines worldwide, by paying more attention to the degree of processing of foods along with nutrient based recommendations. We acknowledge that some progress has been made in this field; the latest version of the dietary guidelines to improve cardiovascular health released by the American Heart Association recommends choosing minimally processed foods instead of ultra-processed foods,60 in accordance with what has been already done in some countries.61 Finally, our findings will hopefully contribute to the ongoing discussions about the potential implementation of a nutrition labelling system at the European Union level.1014
What this study adds
In a large Italian population cohort, the Nutri-Score and the NOVA classification were independently associated with all cause and cardiovascular mortality
Part of the excess mortality risk associated with a nutrient poor diet, as reflected by increased values of the Nutri-Score, was significantly explained by a higher degree of food processing
Ultra-processed food intake, by contrast, remained associated with mortality even after the poor nutritional quality of the diet was accounted for
Data availability statement
The data underlying this article will be shared on reasonable request to the corresponding author. The data are stored in an institutional repository (https://repository.neuromed.it), and access is restricted by the ethical approvals and the legislation of the European Union.
Acknowledgments
We are grateful to the Moli-sani Study participants who enthusiastically joined the study and to the Associazione Cuore Sano ONLUS (Campobasso, Italy) for its support to our research communication activities. ER was supported by Fondazione Umberto Veronesi, which is gratefully acknowledged.
Footnotes
Contributors: MB, LI, ADiC, GdG, and GG conceived and designed the study. ADeC was in charge of the Moli-sani bio-bank and did laboratory tests. SC stored and managed the data. MB, ER, and ADiC analysed the data. MB and GG drafted the manuscript. CC, MBD, GdG, and LI originally promoted the Moli-sani study and critically revised this manuscript. All authors gave final approval and agree to be accountable for all aspects of the work ensuring integrity and accuracy. MB and LI are the guarantors. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.
Funding: The enrolment phase of the Moli-sani Study was supported by research grants from Pfizer Foundation (Rome, Italy), the Italian Ministry of University and Research (MIUR, Rome, Italy)–Programma Triennale di Ricerca, decreto No1588, and Instrumentation Laboratory, Milan, Italy. The analyses reported here were partially supported by a grant to LI (AIRC individual grant - project code: 25942) and by the Italian Ministry of Health (Ricerca Corrente 2022-2024). The funders had no role in considering the study design or in the collection, analysis, or interpretation of data, the writing of the report, or the decision to submit the article for publication. All authors were and are independent from the funders.
Competing interests: All authors have completed the ICMJE uniform disclosure form at https://www.icmje.org/disclosure-of-interest/ and declare: support for the submitted work as detailed above; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
The lead authors (the manuscript’s guarantors) affirm that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.
Dissemination to participants and related patient and public communities: Results from analyses conducted in this paper and their possible impact on real world life will be disseminated to both Moli-sani Study cohort participants and public communities through press releases and social media, as well as by newsletters and regular in-presence meetings.
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
