Association between ultra-processed food consumption and gut microbiota in senior subjects with overweight/obesity and metabolic syndrome

The production and consumption of ultra-processed foods (UPF) has increased considerably during the last years worldwide. Collective evidence shows the association between UPF consumption and adverse health outcomes, including inflammatory gastro-intestinal disorders and obesity. The gut microbiota has been suggested as potential mediator of the effects of UPF consumption on metabolism and health. However, few studies have been conducted in order to elucidate these aspects. Therefore, the aim of the present study was to assess the cross-sectional associations between UPF consumption and gut microbiota in a population of senior subjects (n = 645) within the frame of the PREDIMED-Plus trial. Eligible participants were men and women (aged 55–75 years), without documented history of cardiovascular disease at enrollment, with overweight/obesity (body mass index ≤ 27 and <40 kg/m²) and metabolic syndrome. Using the information of food frequency questionnaires, the consumption of UPF, expressed as a percentage of total dietary energy intake in kcal/day, was calculated considering those food items classified in group 4 of NOVA system. Population was categorized according to tertiles of UPF consumption. Taxonomic fecal microbiota information, along with blood biochemical parameters, anthropometric measurements and clinical data were obtained. Bioinformatics analysis was performed to study the association between fecal microbiota composition and UPF consumption. We observed that subjects allocated in the highest tertile of UPF consumption (21.4 ± 5.0 % kcal/day) presented lower adherence to MedDiet (p < 0.001) and higher total energy intake (p < 0.001). The taxonomic analysis of the fecal microbiota revealed a significant (Benjamini-Hochberg adjusted p < 0.2) positive association between specific taxa and tertiles (T) of UPF consumption: Alloprevotella (p = 0.041 vs. T2; p = 0.065 vs. T3), Negativibacillus (p = 0.096 vs. T3), Prevotella (p = 0.116 vs. T3), and Sutterella (p = 0.116 vs. T2). UPF consumption was positively associated with lower adherence to MedDiet and higher total energy intake in senior subjects with overweight obesity and metabolic syndrome. In addition, positive association with specific fecal microbiota taxa related to inflammatory gastro-intestinal diseases and low consumption of fruits and vegetables, was observed.

Introduction

Production and consumption of ultra-processed food (UPF) has increased considerably during the last years worldwide, displacing foods that are usually consumed in healthy eating patterns, abundant in unprocessed vegetables, fruits, legumes and nuts among others (1, 2).

UPF are defined as multi-ingredient industrial formulations manufactured from substances derived from foods and additives or ingredients not normally used in home food preparation (such as hydrogenated or unesterified oils, protein isolate, maltodextrin, casein, and gluten), retaining little or none of their original food characteristics. UPF include food substances not used in culinary preparations, in particular flavorings, colorings, sweeteners, emulsifiers and other additives used to mimic the sensory qualities of unprocessed or to mask undesirable qualities of the final product (3, 4). They are frequently highly palatable and easy to eat, as they usually can be consumed at the time of purchase or with very little preparation. As a result, the UPF are commonly rich in energy density, added sugars, salt, and saturated and trans fatty acids, and low in fiber, protein, and micronutrients. According to the NOVA food classification system (4, 5) foods and food products can be classified into four categories according to the degree of processing: low or unprocessed foods (i.e., fresh fruits, vegetables or legumes), culinary ingredients (i.e., honey or sugar), processed foods (fruit juices or canned vegetables with salt), and ultra-processed foods (such as pizza or biscuits).

A growing body of evidence shows an association between UPF consumption and an increased risk of chronic diseases such as obesity (5–7), metabolic syndrome (MetS) (5), diabetes (8, 9), dyslipidemia (10), cardiovascular disease (CVD) (5, 11), cancer (12), and all-cause mortality (13, 14). Most of the studies relating UPF consumption and health outcomes are cross-sectional or prospective epidemiologic studies with a relatively low level of scientific evidence of cause-effect relationship. However, some recent clinical studies indicate that the decrease in the amount of UPF consumed should help improve strategies for cardiovascular risk factors in patients in secondary care (14, 15). These findings should serve as an incentive for limiting consumption of UPF, and promote natural or minimally processed foods, as several national nutritional policies recommend (8, 16).

Considering the global increase in consumption of UPFs and the growing evidence linking its excessive consumption to disease outcomes, it is important to outline the possible mechanisms by which this type of food may have detrimental health effects. One of the suggested mechanisms by which UPF consumption may have potential effects on metabolism and health is through the modulation of gut microbiota composition and function (7, 17). In fact, there is growing evidence demonstrating that diet may have an effect on cardiovascular risk factors and metabolism through gut microbiota modulation. Food processing alters the provision of macronutrient substrate to the colon due to differing digestibility, and thereby may impact the microbiota and its metabolic activity (17). Other additives and nutrient constituents have also demonstrated to modulate the host gut microbiota after consumption (7, 18). However, few studies have been conducted in order to explore the relationship between UPF and gut microbiota, and therefore further investigation is required in this sense.

The aim of the present study was to assess the baseline cross-sectional associations between UPF consumption and gut microbiota in a population of senior subjects with overweight/obesity and MetS within the frame of the PREDIMED-Plus trial.

Materials and methods

Study design and participants

The present cross-sectional study was conducted within the frame of the PREDIMED (PREvención con DIeta MEDiterránea)-Plus clinical trial. This trial aims to assess the long-term effects of an intensive weight loss lifestyle intervention based on an energy-restricted Mediterranean diet (MedDiet), physical activity (PA) promotion, and behavioral support (intervention group) on CVD events and mortality, vs. a control group following an energy-unrestricted MedDiet without any advice to increase PA (19). The trial was registered at the International Standard Randomized Controlled Trial (Number: ISRCTN89898870 -data of registration: 2014-) and approved by the institutional review board of all participating institutions. A detailed protocol is available at the web page http://predimedplus.com/. Eligible participants were men and women (aged 55–75 years), without documented history of CVD at enrollment, with overweight/obesity (body mass index (BMI) ≤ 27 and <40 kg/m²) and who met at baseline at least three components of the MetS.

The present study encompasses a subsample of participants from the PREDIMED-Plus recruiting centers of Reus, Barcelona (IMIM), Alicante, and Valencia, with available fecal microbiota data at baseline (n = 656).

Dietary assessments

Dietary intake was assessed using a validated 143 item semiquantitative food frequency questionnaires (FFQs) (20) administered by trained dietitians. Participants reported their average frequency and quantity of food consumed during the previous year. The intake of each item was calculated by multiplying a typical portion size by frequency of consumption (9 possible responses ranging from never to > 6 times/day). Spanish food composition tables (21) were used to derive nutrient (sodium, saturated and trans fatty acids), fiber, alcohol (g/day), and total energy intake (kcal/day), as well as to determine consumption of specific food groups, such as fruits and vegetables (g/day).

Food and beverage items of FFQ were classified in accordance with the NOVA system proposed by Monteiro et al. (22). Foods and drinks were allocated in four different groups according to the nature, extent and purpose of their industrial processing, considering a variety of physical, biological, and chemical methods, as well as use of additives. For the purposes of the current study, we considered food items classified in the NOVA group 4 (G4), including “ultra-processed food and drink products.” These items were further classified into subgroups: dairy products; processed meat; pre-prepared dishes, snacks and fast-foods; sweets; non-alcoholic beverages; alcoholic beverages (Table 1). Using the information of FFQs we calculated the consumption of UPF G4 and expressed as a percentage of total energy intake (in kcal/day). The proportion of UPF in the total diet was calculated for each participant and categorized in tertiles, corresponding to low (first tertile), medium (second tertile), and high (third tertile) consumption.

TABLE 1

Table 1. List of food items from the PREDIMED-Plus food frequency questionnaires allocated in Group 4 (ultra-processed food and drink products) according to NOVA classification system.

General assessments, anthropometric and biochemical measurements

Information about sociodemographic and lifestyle aspects, education level, individual and family medical history, and current medication use was collected. Leisure-time PA was measured using the validated REGICOR questionnaire (23) which included questions to collect information about the type of activity, frequency (number of days) and duration (min/day) performed during a representative month.

Waist circumference was measured midway between the lowest rib and the iliac crest using an anthropometric tape, body weight was measured twice using high-quality electronic calibrated scales and height was measured twice using a wall-mounted stadiometer. Systolic and diastolic blood pressure was measured at rest 3 times using a validated semiautomatic oscillometer (Omron HEM-705CP, Kyoto, Japan) and the mean of repeated measures was used.

Blood samples were collected after an overnight fast. Tubes of serum and plasma were divided into aliquots, coded and stored at −80°C in a central laboratory until analyses. Plasma levels of glucose, insulin, total cholesterol, high density lipoprotein (HDL) cholesterol and triglycerides were measured using standard enzymatic methods, low density lipoprotein (LDL) cholesterol was calculated with the Friedewald formula (whenever triglycerides were <300 mg/dL).

Stool samples collection, DNA extraction and 16S amplicon sequencing

Stool samples were collected at home by participants and kept frozen till the delivery to the laboratory. Stool samples were then separated into 250 mg aliquots and stored at −80°C, until analysis. Patients exposed to antibiotics treatment (N = 6) were excluded from the present analysis.

Microbial DNA was extracted using the QIAmp PowerFecal DNA kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. Previous extraction samples were disrupted in a 5 min lysis step using FastPrep-24™ 5G Homogenizer (MP Biomedicals, Santa Ana, CA, USA). After extraction, the DNA was stored at −20°C until further processing. DNA concentration and purity were assessed with the Qubit 2.0 Fluorometer-dsDNA Broad Range (BR) Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA).

The V4 region (515F-806R) of the prokaryotic 16S ribosomal RNA (rRNA) gene was amplified in triplicate PCR reactions. PCR was performed in a total reaction volume of 35 μl, and the master mix contained 0.7 μl of a unique barcoded primer, 515F-n and 806R-n (10 μM each per reaction), 0.7 μl dNTPs mixture, 0.35 μl Phusion Green Hot Start II High-Fidelity DNA Polymerase (2 U/μl; Thermo Scientific, Landsmeer, The Netherlands), 7 μl 5× Phusion Green HF Buffer, and 25.5 μl DNAse- RNAse-free water (Promega, Madison, WI, USA). The amplification program included 30 s of initial denaturation step at 98°C, followed by 25 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 10 s, elongation at 72°C for 10 s, and a final extension step at 72°C for 7 min. The PCR products were visualized in 1% agarose gel (~290 bp) and purified with CleanPCR kit (CleanNA, Alphen Aan den Rijn, The Netherlands). The DNA concentration of the purified PCR products was measured with Qubit dsDNA BR Assay Kit, and 200 ng from each purified PCR product was used for the construction of the sequencing libraries, which were sequenced on an Illumina Novaseq platform. To ensure the quality of the sequencing data, artificial mock communities with known composition (24) were included in each library as positive controls, and 7 fecal samples were sequenced in duplicates. Finally, to control for potential contaminant sequences, negative control samples were included in each library, using as DNA template, nuclease free-water or material from DNA extraction blanks.

Raw sequence data were processed using the NG-Tax pipeline (24), with a read length of 100 nt. Paired-end libraries were demultiplexed and only read pairs with perfectly matching barcodes were used for downstream steps. Amplicon sequence variants (ASVs), were determined with the default settings and taxonomy was assigned to each ASV, using the USEARCH algorithm (25) and the Silva database (v138.1) (26).

Statistical analysis

The clinical characteristics of the study participants were described using the software IBM SPSS Statistics version 23 (SPSS Inc., Chicago, Illinois, USA). Numerical variables were considered normally distributed according to the central limit theorem and described as means and standard deviations, whereas categorical variables were described as numbers and percentages. Differences across tertiles of UPF consumption (in % of kcal/day from total energy intake) were tested with Pearson's chi square test or one-way ANOVA as appropriate and p < 0.05 deemed as significant.

Fecal microbiota data was analyzed using R (version 4.0.5) and R Studio (version 1.2.5033). Prior filtering steps, a total of 192,011,642 sequence reads for 656 samples was obtained, clustering into 6285 ASVs. These ASVs represented 342 bacterial taxa at genus level. A cut off value of 10% of prevalence at genus level on the absolute abundances of the ASV counts was used to remove ASVs with a prevalence ≤10% between samples. After filtering steps, a total of 179,946,951 sequence reads for 645 samples and 91 genera represented were obtained. The minimum number of sequence reads detected in the dataset was 33, the maximum number of sequence reads was 988,905, the average number of sequence reads was 278,986.

The total relative abundance and percentage of each taxon per tertile was assessed at genus level. Principal component analysis (PCA) was calculated over centered log-ratio (CLR) transformed counts to evaluate the fecal microbiota distribution of the study population among the tertiles of UPF consumption. Fecal microbial richness and diversity was assessed by calculating Chao1, Shannon and Simpson indices on absolute abundance counts and differences across tertiles of UPF consumption tested with one-way ANOVA. Beta diversity was calculated in terms of Euclidean distance over CLR-transformed counts and permutational multivariate analysis of variance (PERMANOVA) performed to test differences in microbiota dissimilarity across tertiles of UPF consumption. The log-normalized Firmicutes-to-Bacteroidetes ratio was computed based on the relative abundance between the phylum Firmicutes and Bacteroidetes and differences across tertiles of UPF consumption tested using one-way ANOVA.

The association between tertiles of UPF consumption and fecal microbiota was assessed using the R package MaAsLin2 (27) (version 1.2.0). Specifically, counts were CLR-transformed, and model adjusted for different potential confounders (age, sex, recruiting center, smoking habits, diabetes prevalence, BMI categories (overweight, with BMI 27–29.9 kg/m²; obese, with BMI 30–39.9 kg/m²; severe obese, with BMI ≥40 kg/m²), and tertiles of PA. Multiple testing correction was performed using Benjamini-Hochberg procedure, and p < 0.2 deemed as significant. Spearman's correlation was used to test the relationship between taxa significantly associated with tertiles of UPF consumption, different G4 UPF items categories and cardiovascular risk factors: body weight, waist circumference, and body mass index (BMI), systolic blood pressure (SBP), diastolic blood pressure (DBP), fasting plasma glucose (FPG), insulin, total cholesterol, HDL cholesterol, LDL cholesterol, triglycerides, and glycated hemoglobin (HbA1c).

Results

General characteristics of the study population

A total of 656 participants in the framework of the PREDIMED-Plus clinical trial with available fecal microbiota data were included in this study. After applying quality filtering cutoffs and excluding those participants exposed to antibiotic treatment, the study population was reduced to 645 individuals stratified according to tertiles of UPF consumption. The characteristics of the study populations are described in Table 2. A total of 215 individuals were allocated in tertile 1 (T1) representing those individuals with lower consumption of UPF (7.2 ± 2.3%); 217 individuals in T2 (13.1 ± 1.6%) and 213 individuals T3 (21.4 ± 5.0%), representing those individuals with higher consumption of UPF. Participants in the highest tertiles of UPF (T2 and T3) consumption tended to be younger, with a higher energy consumption and lower adherence to MedDiet.

TABLE 2

Table 2. Characteristics of the study population according to tertiles of UPF consumption.

From the total energy of UPF consumed, the main contributors (Figure 1) were represented by non-alcoholic beverages (26%), sweets (26%), processed meats (21%), preprepared dishes, snacks and fast-foods (14%), dairy products (11%) and alcoholic beverages (2%).

FIGURE 1

Figure 1. Pie chart showing the proportion (in % kcal/day) of UPF consumption (group 4 of the NOVA classification system) in the total diet of the study population.

Fecal microbiota diversity and dissimilarity and main phyla distribution across study groups

No statistically significant differences in alpha diversity indices calculated, Chao1, Shannon, and Simpson, across tertiles of UPF consumption (Supplementary Figure 1 and Supplementary Table 1) were observed. PCA over the CLR-transformed taxonomic feature counts shown that PC1 and PC2 account ~10.15% and 8% to the total variation, respectively. PCA plot showed that the fecal microbiota samples did not cluster based on tertiles of UPF consumption (Supplementary Figure 2). Furthermore, the result of the PERMANOVA test based on Euclidean distance did not show statistically significant differences across tertiles of percentage in UPF energy intake (Supplementary Table 2). There were no differences in Firmicutes-to-Bacteroidetes ratio across tertiles UPF consumption (Supplementary Table 3).

Association between specific fecal microbiota taxa and UPF consumption

The analysis performed with MaAsLin2 shown a positive association between fecal microbiota genera Alloprevotella and Sutterella, representing, respectively, 1.2% and 0.6% of the total relative abundance (Supplementary Table 4), and T2 of percentage in UPF energy intake, and a positive association between fecal microbiota genera Alloprevotella, Negativibacillus, and Prevotella, representing respectively 0.3%, 0.03% and 0.3% of the total relative abundance (Supplementary Table 4), and T3 of percentage in UPF energy intake (Table 3).

TABLE 3

Table 3. Results of the association analysis performed with MaAsLin2 between taxonomic feature counts at genus level and tertiles of percentage in UPF energy intake.

In addition, the results of the Spearman correlation did not show any relevant association between taxa significantly associated with UPF consumption and UPF items categories (Supplementary Figure 3 and Supplementary Table 5) nor with cardiovascular risk factors (Supplementary Figure 4 and Supplementary Table 6).

Discussion

The current cross-sectional study aimed to assess the association between UPF consumption and fecal microbiota composition in senior subjects with overweight/obesity and MetS within the frame of the PREDIMED-Plus trial. A previous study conducted in the whole PREDIMED-Plus cohort shown that higher consumption of UPF was associated with higher visceral and overall adiposity accumulation (28). Similar results have also been reported some years before by Mendoça et al., in the context of a prospective cohort study in healthy middle-aged Spanish individuals, in which those participants in the highest quartile of UPF consumption showed an increased risk of developing overweight or obesity (29). A recent comprehensive systematic review explored extant relevant nutrition literature reporting increasing evidence of UPF exposure associated with adverse health outcomes such as overweight, obesity and cardio-metabolic risk (30), suggesting that a high frequency of UPF consumption could be related to obesity and several cardiometabolic alterations.

In the present study, the population was categorized according to tertiles of UPF consumption, and we observed that subjects allocated in the highest tertile presented lower adherence to MedDiet and higher total energy intake. These results were in line with the aforementioned articles, and very similar to the ones obtained previously by Cuevas-Sierra et al. where higher values of total energy intake and a lower adherence to MedDiet were observed in subjects consuming more than 5 UPF servings/day (18).

According to our findings, we did not observe any association between alpha and beta diversity nor with Firmicutes-to-Bacteroidetes ratio in our study population. In relation to that, it has been reported that food processing may disturb gut microbiota composition also reducing bacterial diversity. However, the mechanisms by which UPF affects gut microbial homeostasis remain unclear (31). In addition, there are studies linking a varied and omnivorous diet, and consequently higher adherence to MedDiet, with an increase in the ratio Firmicutes-to-Bacteroidetes (32, 33). In line with these studies, Partula et al. aimed to explore the association between usual diet and gut microbiota composition, established through 16S rRNA gene sequencing in stool samples from healthy French adults (34). They observed that food items for which a limited consumption is generally recommended (fried products, sodas or sugary drinks, fatty sweet products, processed meats, ready-cooked meals, desserts, and cheese) were negatively associated with alpha-diversity and contributed to shifts within microbiota composition (beta-diversity). These contradictory results in relation to the aforementioned studies can be explained by a lower consumption of UPF in our senior study population.

In the present study we observed a positive association between 4 different taxa abundances and higher UPF consumption. Specifically, we found that abundance of genera Alloprevotella, Negativibacillus, Prevotella and Sutterella, was positively associated with the highest consumption of UPF. These positive associations are in accordance with other recently published studies. The association between specific food items and relative abundances of taxa tested by Partula et al. shown that Alloprevotella was negatively associated with the consumption of fruits (34). Alloprevotella was found enriched in the gut microbial community determined by absolute quantification 16S rRNA gene amplicon sequencing of healthy individuals exposed to 4 days dietary fiber intervention (35). However, these findings should be interpreted with caution as obtained in a small sample size.

Negativibacillus has been recently discovered in the microbial community of human stool samples and it has been found enriched in patients with ulcerative colitis (UC) (36). There is not much information on this genus in the bibliography, however, in a previous published study it has been found that the abundance of Negativibacillus was positively associated with obesity in mice (37).

Prevotella species, as Negativibacillus, have been associated with patients with UC (36). This genus has also been associated with plant-based diet, rich in fiber and carbohydrates, and has been considered as a biomarker of health. Nevertheless, there are important confounding factors such as microbial exposure, lifestyle and geography that may affect this associative evidence (38). In addition, the potential pathogenic role of certain Prevotella strains has been also highlighted (39). The strain-level diversity may explain the controversial role of Prevotella species in health and disease (40). This could partially explain the association that we found between this taxon and higher UPF consumption, as the response to different dietary patterns may be strain dependent and undetectable due to the next generation sequencing techniques utilized that limited us to genus level resolution, as also emphasized by De Filippis et al. (40).

Finally, Sutterella, has been recently linked to inflammatory gastrointestinal diseases (41–44). The relationship between diet and development of inflammatory gastrointestinal tract disorders has been indicated in several epidemiological studies (45). Previous meta-analysis studies showed that soft drink consumption and sucrose intake were associated with 69% and 10% increased risk of UC development, respectively (46, 47). In addition, consumption of fruits and vegetables was related to decreased UC development (48). This could in part explain the positive association between Sutterella and UPF consumption in our cohort.

Despite these findings, in the current study we did not observe any correlation between the above-mentioned UPF-associated taxa, UPF items categories and CVD risk factors, respectively. Nevertheless, lower adherence to MedDiet was observed in those participants in the highest tertile of UPF energy consumption. These results were in consonance with previous studies (49–51).

Our study has some limitations that deserve to be discussed. First, the observational design of the study did not make possible to conclude on causality or directionality as these results indicate just associations. However, observational design is the most indicated approach to explore these type of associations, as longitudinal studies based on interventions including UPF intake may be detrimental for health and could cause ethical issues. Secondly, the results cannot be extrapolated to other populations since participants included in the analysis were senior Mediterranean subjects with overweight/obesity and MetS. Under the methodological point of view, the nature of 16S sequencing limits the taxonomic profiling to genus-level resolution. Finally, the assessment of food intake through a FFQ is sensitive to possible measurements errors. Nevertheless, food-based FFQs have been widely used as a tool in epidemiological studies since the 1990's (52).

On the other hand, the current study also has some strengths, such as the large cohort sample size, the extensive data collection, and the use of several potential confounding factors which minimizes the possibility of reverse causality bias. Finally, it is also important to highlight that the study setup provides an ideal scenario to explore specific gut microbial profiles related to dietary patterns.

Conclusions

In conclusion, the results of this present observational study conducted within the frame of the PREDIMED-Plus cohort revealed that the consumption of UPF (rich in energy density, added sugars, salt, and saturated and trans fatty acids, and low in fiber, protein, and micronutrients) was positively associated with lower adherence to MedDiet and higher total energy intake. In addition, we observed a positive association with specific taxa (Alloprevotella, Negativibacillus, Prevotella and Sutterella) related to inflammatory gastro-intestinal diseases and low consumption of fruits and vegetables, among others, across tertiles of UPF consumption.

Diet and nutritional status are important determinants in human health. The role of diet in shaping the gut microbiota could be key to improve host's health. Consequently, further future studies including whole metagenome shotgun sequencing would be crucial to obtain more detailed taxonomic identification. Furthermore, information generated from metabolomics analysis would be useful to improve and strengthen the results. In addition, more in-depth studies regarding the association between taxa and UPF consumption are needed in order to further reduce the controversies between the authors. The detection of unhealthy dietary patterns and association with gut microbial profiles would be essential for preventing related chronic diseases and to improve public health strategies.

Data availability statement

The datasets generated and analyzed during the current study are not publicly available due to data regulations and for ethical reasons, considering that this information might compromise research participants' consent because our participants only gave their consent for the use of their data by the original team of investigators. However, collaboration for data analyses can be requested by sending a letter to the PREDIMED-Plus steering Committee (cHJlZGltZWRfcGx1c19zY29tbWl0dGVlJiN4MDAwNDA7Z29vZ2xlZ3JvdXBzLmNvbQ==). The request will then be passed to all the members of the PREDIMED-Plus Steering Committee for deliberation.

Ethics statement

The studies involving human participants were reviewed and approved by the study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Hospital Universitari Sant Joan de Reus Ethics Committee (protocol code 13-07-25/7proj2, date of approval: 25/07/2013), and the Comité de Ética de la Investigación Provincialde Málaga (protocol code Predimed+DM/01, date of approval: 27/11/ 2014). The patients/participants provided their written informed consent to participate in this study. The patients/participants provided their written informed consent to participate in this study.

Author contributions

All the principal PREDIMED-Plus investigators from the respective centers contributed to study concept and design and to data extraction from the participants (FT, JVio, DC, MF, JVid, and JS-S). AA and MM performed the statistical analyses. AA, MM, and JS-S drafted the manuscript. All authors reviewed the manuscript for important intellectual content and approved the final version to be published.

Funding

This work was supported by the official Spanish Institutions for funding scientific biomedical research, CIBER Fisiopatología de la Obesidad y Nutrición (CIBEROBN) and Instituto de Salud Carlos III (ISCIII), through the Fondo de Investigación para la Salud (FIS), which is co-funded by the European Regional Development Fund (three coordinated FIS projects leaded by JS-S, including the following projects: PI13/00462, PI16/00501 and PI19/00576); The Generalitat Valenciana PROMETEO 17/2017, PROMETEO 21/2021 and APOSTD/2020/164 (to AÁ-S). This research was also partially funded by EU-H2020 Grant Eat2beNICE/H2020-SFS-2016-2 and by NIH grant R01DK127601. None of the funding sources took part in the design, collection, analysis, interpretation of the data, or writing the report, or in the decision to submit the manuscript for publication. Food companies Hojiblanca (Lucena, Spain) and Patrimonio Comunal Olivarero (Madrid, Spain) donated extra virgin olive oil; and the Almond Board of California (Modesto, CA), American Pistachio Growers (Fresno, CA), and Paramount Farms (Wonderful Company, LLC, Los Angeles, CA) donated nuts for the PREDIMED-Pilot study. AA was supported by the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 713679 and from the Universitat Rovira i Virgili (URV). MM was supported by Sara Borrell -CD21/00045- Instituto de Salud Carlos III (ISCIII). IM-I was supported by a Miguel Servet II contract -CPII21-00013- Instituto de Salud Carlos III (ISCIII).

Acknowledgments

We thank all the volunteers for the participation and personnel for the contribution in the PREDIMED-Plus trial. We also thank all the investigators of the PREDIMED-Plus study. CIBEROBN (Centros de Investigación Biomédica en Red: Obesidad y Nutrición) is an initiative of ISCIII, Madrid, Spain. The authors also thank the PREDIMED-Plus Biobank Network as a part of the National Biobank Platform of the ISCIII for storing and managing the PREDIMED-Plus biological samples. JS-S senior author of this study gratefully acknowledges the financial support by ICREA under the ICREA Academia programme. The research groups want to thank Jasper Koehorst and Bart Nijsse from the Unlock platform for their support in processing the data with the NG-tax pipeline. We also want to thank Athanasia Ioannou for her work with the stool samples preparation.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.976547/full#supplementary-material

References

1. Salas-Salvadó J, Martinez-González MÁ, Bulló M, Ros E. The role of diet in the prevention of type 2 diabetes. Nutr Metab Cardiovasc Dis. (2011) 21:09. doi: 10.1016/j.numecd.2011.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Martínez Leo EE, Segura Campos MR. Effect of ultra-processed diet on gut microbiota and thus its role in neurodegenerative diseases. Nutrition. (2020) 71:110609. doi: 10.1016/j.nut.2019.110609

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Steele EM, Baraldi LG, Da Costa Louzada ML, Moubarac JC, Mozaffarian D, Monteiro CA. Ultra-processed foods and added sugars in the US diet: Evidence from a nationally representative cross-sectional study. BMJ Open. (2016) 6:1–8. doi: 10.1136/bmjopen-2015-009892

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Monteiro CA, Cannon G, Moubarac J-C, Levy RB, Louzada MLC, Jaime PC. The UN decade of nutrition, the NOVA food classification and the trouble with ultra-processing. Public Health Nutr. (2018) 21:5–17. doi: 10.1017/S1368980017000234

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Martinez-Perez C, San-Cristobal R, Guallar-Castillon P, Martínez-González MÁ, Salas-Salvadó J, Corella D, et al. Use of different food classification systems to assess the association between ultra-processed food consumption and cardiometabolic health in an elderly population with metabolic syndrome (PREDIMED-Plus Cohort). Nutrients. (2021) 13:2471. doi: 10.3390/nu13072471

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Nardocci M, Leclerc B-S, Louzada M-L, Monteiro CA, Batal M, Moubarac J-C. Consumption of ultra-processed foods and obesity in Canada. Can J Public Heal. (2019) 110:4–14. doi: 10.17269/s41997-018-0130-x

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Cuevas-Sierra A, Ramos-Lopez O, Riezu-Boj JI, Milagro FI, Martinez JA. Diet, gut microbiota, and obesity: links with host genetics and epigenetics and potential applications. Adv Nutr. (2019) 10:S17–30. doi: 10.1093/advances/nmy078

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Srour B, Fezeu LK, Kesse-Guyot E, Allès B, Debras C, Druesne-Pecollo N, et al. Ultraprocessed food consumption and risk of type 2 diabetes among participants of the nutrinet-santé prospective cohort. JAMA Intern Med. (2020) 180:283. doi: 10.1001/jamainternmed.2019.5942

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Llavero-Valero M, Escalada-San Martín J, Martínez-González MA, Basterra-Gortari FJ, de la Fuente-Arrillaga C, Bes-Rastrollo M. Ultra-processed foods and type-2 diabetes risk in the SUN project: A prospective cohort study. Clin Nutr. (2021) 40:2817–24. doi: 10.1016/j.clnu.2021.03.039

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Donat-Vargas C, Sandoval-Insausti H, Rey-García J, Moreno-Franco B, Åkesson A, Banegas JR, et al. High consumption of ultra-processed food is associated with incident dyslipidemia: a prospective study of older adults. J Nutr. (2021) 151:2390–8. doi: 10.1093/jn/nxab118

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Srour B, Fezeu LK, Kesse-Guyot E, Allès B, Méjean C, Andrianasolo RM, et al. Ultra-processed food intake and risk of cardiovascular disease: prospective cohort study (NutriNet-Santé). BMJ. (2019) 2019:1451. doi: 10.1136/bmj.l1451

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Fiolet T, Srour B, Sellem L, Kesse-Guyot E, Allès B, Méjean C, et al. Consumption of ultra-processed foods and cancer risk: results from NutriNet-Santé prospective cohort. BMJ. (2018) 2018:k322. doi: 10.1136/bmj.k322

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Hall KD, Ayuketah A, Brychta R, Cai H, Cassimatis T, Chen KY, et al. Ultra-processed diets cause excess calorie intake and weight gain: an inpatient randomized controlled trial of ad libitum food intake. Cell Metab. (2019) 30:67–77.e3. doi: 10.1016/j.cmet.2019.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Schnabel L, Kesse-Guyot E, Allès B, Touvier M, Srour B, Hercberg S, et al. Association between ultraprocessed food consumption and risk of mortality among middle-aged adults in France. JAMA Intern Med. (2019) 179:490. doi: 10.1001/jamainternmed.2018.7289

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Beslay M, Srour B, Méjean C, Allès B, Fiolet T, Debras C, et al. Ultra-processed food intake in association with BMI change and risk of overweight and obesity: A prospective analysis of the French NutriNet-Santé cohort. Souza Lopes AC, editor. PLOS Med. (2020) 17:e1003256. doi: 10.1371/journal.pmed.1003256

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Bonaccio M, Di Castelnuovo A, Costanzo S, De Curtis A, Persichillo M, Sofi F, et al. Ultra-processed food consumption is associated with increased risk of all-cause and cardiovascular mortality in the Moli-sani Study. Am J Clin Nutr. (2021) 113:446–55. doi: 10.1093/ajcn/nqaa299

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Muralidharan J, Moreno-Indias I, Bulló M, Lopez JV, Corella D, Castañer O, et al. Effect on gut microbiota of a 1-y lifestyle intervention with Mediterranean diet compared with energy-reduced Mediterranean diet and physical activity promotion: PREDIMED-Plus Study. Am J Clin Nutr. (2021) 114:1148–58. doi: 10.1093/ajcn/nqab150

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Cuevas-Sierra A, Milagro FI, Aranaz P, Martínez JA, Riezu-Boj JI. Gut microbiota differences according to ultra-processed food consumption in a Spanish population. Nutrients. (2021) 13:1–20. doi: 10.3390/nu13082710

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Salas-Salvadó J, Díaz-López A, Ruiz-Canela M, Basora J, Fitó M, Corella D, et al. Effect of a lifestyle intervention program with energy-restricted mediterranean diet and exercise on weight loss and cardiovascular risk factors: one-year results of the PREDIMED-Plus trial. Diabetes Care. (2018) 42:dc180836. doi: 10.2337/dc18-0836

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Fernández-Ballart JD, Piñol JL, Zazpe I, Corella D, Carrasco P, Toledo E, et al. Relative validity of a semi-quantitative food-frequency questionnaire in an elderly Mediterranean population of Spain. Br J Nutr. (2010) 103:1808–16. doi: 10.1017/S0007114509993837

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Moreiras O, Carbajal Á, Cabrera L, Cuadrado C. Tablas de composición de alimentos. 19a edició. Pirámide (2018).

Google Scholar

22. Monteiro CA, Cannon G, Jaime PC, Canella DS, Louzada ML, Calixto G, et al. Food Classification. Public Health NOVA. The Star Shines Bright. (2016).

Google Scholar

23. Molina L, Sarmiento M, Peñafiel J, Donaire D, Garcia-Aymerich J, Gomez M, et al. Validation of the regicor short physical activity questionnaire for the adult population. Lucía A, editor. PLoS ONE. (2017) 12:e0168148. doi: 10.1371/journal.pone.0168148

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ramiro-Garcia J, Hermes GDA, Giatsis C, Sipkema D, Zoetendal EG, Schaap PJ, et al. NG-Tax, a highly accurate and validated pipeline for analysis of 16S rRNA amplicons from complex biomes. F1000Research. (2018) 5:1791. doi: 10.12688/f1000research.9227.2

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. (2010) 26:2460–1. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. (2012) 41:D590–6. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Mallick H, Rahnavard A, McIver LJ, Ma S, Zhang Y, Nguyen LH, et al. Multivariable association discovery in population-scale meta-omics studies. Coelho LP, editor. PLoS Comput Biol. (2021) 17:e1009442. doi: 10.1371/journal.pcbi.1009442

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Konieczna J, Morey M, Abete I, Bes-Rastrollo M, Ruiz-Canela M, Vioque J, et al. Contribution of ultra-processed foods in visceral fat deposition and other adiposity indicators: Prospective analysis nested in the PREDIMED-Plus trial. Clin Nutr. (2021) 40:4290–300. doi: 10.1016/j.clnu.2021.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Mendonça R de D, Pimenta AM, Gea A, de la Fuente-Arrillaga C, Martinez-Gonzalez MA, Lopes ACS, et al. Ultraprocessed food consumption and risk of overweight and obesity: the University of Navarra Follow-Up (SUN) cohort study. Am J Clin Nutr. (2016) 104:1433–40. doi: 10.3945/ajcn.116.135004

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Elizabeth L, Machado P, Zinöcker M, Baker P, Lawrence M. Ultra-processed foods and health outcomes: A narrative review. Nutrients. (2020) 12:1–36. doi: 10.3390/nu12071955

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Miclotte L, Van de Wiele T. Food processing, gut microbiota and the globesity problem. Crit Rev Food Sci Nutr. (2020) 60:1769–82. doi: 10.1080/10408398.2019.1596878

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Cust AE, Skilton MR, van Bakel MME, Halkjær J, Olsen A, Agnoli C, et al. Total dietary carbohydrate, sugar, starch and fibre intakes in the European Prospective Investigation into Cancer and Nutrition. Eur J Clin Nutr. (2009) 63:S37–60. doi: 10.1038/ejcn.2009.74

PubMed Abstract | CrossRef Full Text | Google Scholar

33. De Filippis F, Pellegrini N, Vannini L, Jeffery IB, La Storia A, Laghi L, et al. High-level adherence to a Mediterranean diet beneficially impacts the gut microbiota and associated metabolome. Gut. (2016) 65:1812–21. doi: 10.1136/gutjnl-2015-309957

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Partula V, Mondot S, Torres MJ, Kesse-Guyot E, Deschasaux M, Assmann K, et al. Associations between usual diet and gut microbiota composition: results from the Milieu Intérieur cross-sectional study. Am J Clin Nutr. (2019) 109:1472–83. doi: 10.1093/ajcn/nqz029

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Tian T, Zhang X, Luo T, Wang D, Sun Y, Dai J. Effects of short-term dietary fiber intervention on gut microbiota in young healthy people. Diabetes, Metab Syndr Obes Targets Ther. (2021) 14:3507–16. doi: 10.2147/DMSO.S313385

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Gryaznova MV, Solodskikh SA, Panevina AV, Syromyatnikov MY, Dvoretskaya YD, Sviridova TN, et al. Study of microbiome changes in patients with ulcerative colitis in the Central European part of Russia. Heliyon. (2021) 7:e06432. doi: 10.1016/j.heliyon.2021.e06432

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Wang W, Pan Y, Wang L, Zhou H, Song G, Wang Y, et al. Optimal dietary ferulic acid for suppressing the obesity-related disorders in leptin-deficient obese C57BL/6J -ob/ob mice. J Agric Food Chem. (2019) 67:4250–8. doi: 10.1021/acs.jafc.8b06760

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Dahl WJ, Rivero Mendoza D, Lambert JM. Diet, Nutrients and the Microbiome. Progress in Molecular Biology and Translational Science. Elsevier Inc. (2020). doi: 10.1016/bs.pmbts.2020.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Precup G, Vodnar DC. Gut Prevotella as a possible biomarker of diet and its eubiotic versus dysbiotic roles: A comprehensive literature review. Br J Nutr. (2019) 122:131–40. doi: 10.1017/S0007114519000680

PubMed Abstract | CrossRef Full Text | Google Scholar

40. De Filippis F, Pasolli E, Tett A, Tarallo S, Naccarati A, De Angelis M, et al. Distinct genetic and functional traits of human intestinal prevotella copri strains are associated with different habitual diets. Cell Host Microbe. (2019) 25:444–53.e3. doi: 10.1016/j.chom.2019.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Mukhopadhya I, Hansen R, Nicholl CE, Alhaidan YA, Thomson JM, Berry SH, et al. A comprehensive evaluation of colonic mucosal isolates of sutterella wadsworthensis from inflammatory bowel disease. Bereswill S, editor. PLoS One. (2011) 6:e27076. doi: 10.1371/journal.pone.0027076

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Hiippala K, Kainulainen V, Kalliomäki M, Arkkila P, Satokari R. Mucosal prevalence and interactions with the epithelium indicate commensalism of Sutterella spp. Front Microbiol. (2016) 7:1–13. doi: 10.3389/fmicb.2016.01706

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Kaakoush NO. Sutterella species, IgA-degrading bacteria in ulcerative colitis. Trends Microbiol. (2020) 28:519–22. doi: 10.1016/j.tim.2020.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Derqaoui S, Oukessou M, Attrassi K, Elftouhy FZ, Nassik S. Detection of Sutterella spp. in Broiler Liver and Breast. Front Vet Sci. (2022) 9:e859902. doi: 10.3389/fvets.2022.859902

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Hou JK, Abraham B, El-Serag H. Dietary intake and risk of developing inflammatory bowel disease: a systematic review of the literature. Am J Gastroenterol. (2011) 106:563–73. doi: 10.1038/ajg.2011.44

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Wang F, Feng J, Gao Q, Ma M, Lin X, Liu J, et al. Carbohydrate and protein intake and risk of ulcerative colitis: Systematic review and dose-response meta-analysis of epidemiological studies. Clin Nutr. (2017) 36:1259–65. doi: 10.1016/j.clnu.2016.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Nie J-Y, Zhao Q. Beverage consumption and risk of ulcerative colitis. Medicine. (2017) 96:e9070. doi: 10.1097/MD.0000000000009070

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Li F, Liu X, Wang W, Zhang D. Consumption of vegetables and fruit and the risk of inflammatory bowel disease. Eur J Gastroenterol Hepatol. (2015) 27:623–30. doi: 10.1097/MEG.0000000000000330

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Gomez-Arango LF, Barrett HL, McIntyre HD, Callaway LK, Morrison M, Dekker Nitert M. Increased systolic and diastolic blood pressure is associated with altered gut microbiota composition and butyrate production in early pregnancy. Hypertension. (2016) 68:974–81. doi: 10.1161/HYPERTENSIONAHA.116.07910

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Mendonça R de D, Lopes ACS, Pimenta AM, Gea A, Martinez-Gonzalez MA, Bes-Rastrollo M. Ultra-processed food consumption and the incidence of hypertension in a mediterranean cohort: the seguimiento universidad de navarra project. Am J Hypertens. (2016) 2016:hpw137. doi: 10.1093/ajh/hpw137

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Ascher S, Reinhardt C. The gut microbiota: An emerging risk factor for cardiovascular and cerebrovascular disease. Eur J Immunol. (2018) 48:564–75. doi: 10.1002/eji.201646879

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Shim J-S, Oh K, Kim HC. Dietary assessment methods in epidemiologic studies. Epidemiol Health. (2014) 2014:e2014009. doi: 10.4178/epih/e2014009

PubMed Abstract | CrossRef Full Text | Google Scholar