Next Article in Journal
Extraction and Quantitation of Phytosterols from Edible Brown Seaweeds: Optimization, Validation, and Application
Next Article in Special Issue
Red Meat Derived Glycan, N-acetylneuraminic Acid (Neu5Ac) Is a Major Sialic Acid in Different Skeletal Muscles and Organs of Nine Animal Species—A Guideline for Human Consumers
Previous Article in Journal
Hyperspectral Imaging with Machine Learning Approaches for Assessing Soluble Solids Content of Tribute Citru
Previous Article in Special Issue
Sustainable Future Protein Foods: The Challenges and the Future of Cultivated Meat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Red and White Meat Intake in Relation to Gut Flora in Obese and Non-Obese Arab Females

1
Department of Community Health Sciences, College of Applied Medical Sciences, King Saud University, Turki Alawwal street, Riyadh 11451, Saudi Arabia
2
Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, King Saud University, Riyadh 11433, Saudi Arabia
3
Chair for Biomarkers of Chronic Diseases, Biochemistry Department, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Foods 2023, 12(2), 245; https://doi.org/10.3390/foods12020245
Submission received: 30 October 2022 / Revised: 20 December 2022 / Accepted: 28 December 2022 / Published: 5 January 2023
(This article belongs to the Special Issue Animal-Based Food Consumption - Trends and Perspectives)

Abstract

:
Background: high meat intake may contribute to several chronic diseases including obesity. However, evidence is insufficient on the relation between red/white meat intake and gut flora among individuals with varying degrees of adiposity. Objective: investigate the association of red/white meat intake with gut flora in Saudi Arabian females with/without obesity. Methods: this observational study involved 92 females with and without obesity (n = 44, 48, respectively) aged 19–25 years. The whole-genome shotgun technique was used to analyze the gut flora. Shannon alpha and Bray–Curtis beta diversity as well as correlation coefficients were used. Results: in the total sample, there were positive correlations between Actinobacteria, Bacteroides (p ≤ 0.05), Flavonifractor plautii (p ≤ 0.0001), and total red meat intake. There were also positive correlations between total white meat intake, Bacteroides, and Faecalibacterium prausnitzii (p ≤ 0.05) in the total sample. In the group without obesity, there was a positive correlation between low white meat intake and Actinobacteria (p = 0.05). In the group with obesity, there was a positive correlation between high white meat intake and Bacteroides (p ≤ 0.001). Conclusion: our findings suggest that meat intake had an impact on the gut flora of Arab adult females, independent of adiposity. Specific strains identified in this study need further investigation to determine their relation to meat intake and obesity.

1. Introduction

In the past decade, global dietary patterns have shifted toward Western-style diets, characterized by higher consumption of refined carbohydrates, added sugars, and animal food sources, including eggs, red meat, processed meat, and white meat [1,2,3]. In 2018, the mean global intake of red meat was ≥50 g/day, having increased by about 88% since 1990 [4]. In Saudi Arabia, the average red meat consumption was high, amounting to ≈50 kg per capita of meat per year in 2017 [5], with an annual growth rate of about 0.24% since 1999 [6]. This increase in meat intake is concurrent with the progressive increase in obesity from 22% in 1990–1993 [7] to 40% in 2017 [8]. High intake of meat has been shown to increase weight gain due to its high energy density and fat content [9]. Compared with other food groups, meat intake was most highly correlated with the prevalence of obesity [10] and other adverse health conditions, such as diabetes [11], cardiovascular disease (CVD) [12], and increased proinflammatory blood markers such as C-reactive protein (CRP), a predictor for CVD risk [13], especially in women [13,14].
Recent studies suggest that diet is one of the most important factors that influence the composition and diversity of the gut flora, eventually modulating the risk of several chronic diseases [15]. The composition of gut bacteria is affected partly by the microorganisms in the ingested food, and partly by the host’s dietary behavior and lifestyle. For example, rats consuming a high-fat diet have a decreased relative abundance of Bacteroidetes and Bifidobacteria and an increased relative abundance of Firmicutes and Proteobacteria [16]. Randomized control trials (RCTs) have reported that meat protein intake correlates with overall microbial diversity [17,18]. However, very few RCTs have investigated the impact of different types of meat on gut flora [17,18,19,20,21], with some revealing an increase in Bacteroides enterotype and a decrease in Firmicutes [17,20] and Clostridium [18] after a high red meat meal. Other RCTs showed that the effect of meat on gut flora may be dependent on the type of meat, i.e., red meat or processed meat [19,21]. It has been suggested that microbial metabolites originating from meat protein, such as short-chain fatty acids, could affect gut flora indirectly by promoting the release of 5-serotonin, leading to raised intestinal motility and ion transport as well as alterations in the gut microbiota composition [22].
Nonetheless, the human gut flora has been observed to differ according to race, ethnicity, and gender, in addition to lifestyle and dietary factors [23,24]. Given the rapid shift in dietary patterns in Saudi Arabia, from a nutrient-rich traditional diet to an energy-dense Western diet [3], and the rising obesity rates, especially in females, it is important to identify the impact of meat intake on gut microbiota composition. To our knowledge, limited evidence is available on the association between the intake of different types of meat and gut flora in the Middle Eastern region. We aim to identify the impact of red/white meat intake on gut flora and the inflammatory marker hs-CRP in Saudi Arabian females with and without obesity. The current study is part of a previous study that aimed to identify the relation between gut microbiota and adiposity. In this study, we focus mainly on dietary factors related to gut flora, specifically meat intake, using the unique whole-genome shotgun (WGS) sequencing technique.

2. Materials and Methods

2.1. Study Design

The current study is part of a case-control study conducted between January 2019 and March 2020 at the clinic of the College of Applied Medical Sciences in King Saud University, Riyadh, Saudi Arabia, and aiming to identify the relation between gut microbiota and obesity markers. Details on the study design were published previously [25]. Briefly, the sample size was 92, based on the gut microbiota composition and its differences among Saudi females attending university at a 5% significance level and 80% power. The sample was calculated in a similar manner to a previous study [26]. We observed a Firmicutes: Bacteroidetes of 0.9 ± 0.4 in women with normal weight and 1.7 ± 1.7 in women with obesity, with a 95% confidence interval (CI) and 80% power. The participants were randomly recruited via flyers, faculty member assistance, oral presentations, and social media networks. We excluded those who were ≤18 years old, overweight (body mass index (BMI) 25.0–29.9 kg/m2), pregnant, or following specific diets (e.g., calorie-restricted diets), as well as those who reported the presence of gastrointestinal diseases in the past eight weeks, endocrine or oncological disease history, psychiatric disorders, anorexia, other medical conditions, and usage of multi-vitamins, vitamin B12, or antibiotics in the past 6 months (n = 193).
After submitting a signed consent form, each participant was given an appointment at the nutrition clinic in the same college and supplied with containers for collecting stool samples that were to be returned on a scheduled visit day. During the visit, demographic, dietary, and anthropometric data were collected, as well as fasting blood samples. Figure S1 represents the number of individuals assessed for participation in the study. The final sample included 92 Saudi female students aged 18 to 25 and classified into obese (BMI ≥ 30 kg/m2, n = 44) and non-obese (BMI = 18.50–24.99 kg/m2, n = 48). The study protocol was approved by the Institutional Review Board Committee of King Khalid University Hospital at King Saud University (IRB #E-19-3625). Clinical trial registration: URL: https://www.clinicaltrials.gov (accessed on 15 October 2022), unique identifier: NCT05664321.

2.2. Anthropometric Measurements

Anthropometric measurements were collected by trained staff using standardized methods. Measurements were recorded twice, and the average results were used for the final analysis. Body weight was measured and recorded to the nearest 0.10 kg and height to the nearest 0.50 cm using an international standard scale (Digital Pearson Scale, ADAM Equipment Inc., Oxford, CT, USA). BMI was calculated by dividing the weight in kilograms by the squared height in meters. BMI was classified into two groups: normal weight (18.50–24.99 kg/m2) and obese (≥30 kg/m2), using the World Health Organization (WHO) criteria [27].
The waist and hip circumferences were measured using a non-stretchable tape. The waist circumference was measured at the narrowest point between the lowest rib and the umbilicus, and hip circumference at the point of the great trochanter, with the two measurements recorded to the nearest 0.50 cm. The waist–hip ratio (WHR) was calculated by dividing the mean waist circumference by the mean hip circumference [28] and categorized as follows: (i) normal WHR (<0.83) and (ii) high WHR (≥0.83) [28]. Body composition, including fat percentage (BF%) and muscle mass, was measured using a bioelectrical impedance analysis method (770 BIA, Inbody, Seoul, South Korea) [29]. The categories are as follows: (i) normal BF% (≤35%) and (ii) high BF% (>35%) [30].

2.3. Dietary Data

Dietary data were collected by trained dietitians during a structured interview. The validated Saudi Food and Drug Authority Food Frequency Questionnaire (FFQ) was used to assess the participants’ dietary intake [31]. The participants were asked to report their frequency of intake of each food item during the past year. The questionnaire was provided in the Arabic language and included 133 food items [31]. Food modules were used to help the participants to determine portion sizes of meat and other foods. ‘Red meat’ referred to beef, camel, goat, or lamb, while ‘white meat’ indicated poultry (chicken, duck, and turkey) and fish [32]. ‘Total meat’ was the total of these two categories. Two 24-h recalls on consecutive days were completed by 20% of the participants to assess the accuracy of the FFQ. Nutrient intake was assessed by the FFQ, and the 24-h recalls correlated well, with a correlation coefficient ranging between 0.50–0.60 for macronutrients and 0.40–0.70 for micronutrients. ESHA Food Processor Software version 11.1 (ESHA Research, Salem, OR, USA) was used for nutrition analysis.

2.4. Lipid Profile and Hs-CRP Test

Blood samples were collected after ≥10 h of fasting from the cubital vein into two 5-mL tubes, an ethylenediaminetetraacetic acid tube for the whole blood sample and a gel tube for the serum sample. Serum samples were transferred on the day of collection to the study laboratory and stored at −80 °C to facilitate their availability for further analysis. Lipid profiles were measured using a biochemical analyzer (Konelab, Espoo, Finland). Friedewald’s equation was used to calculate low-density lipoprotein cholesterol (LDL-C) [33]. Serum hs-CRP was measured using commercial enzyme-linked immunosorbent assay kits.

2.5. Stool Analysis

Each fecal sample was collected in a clean, dry screw-top container and stored directly at −80 °C. All the samples were then transferred to the study laboratory and kept at −80 °C for further analysis. Later, the DNA was extracted from 0.25 g frozen stool aliquots using the QIAGEN PowerFecal DNA Kit (Catalogue: 12830-50). The purity of the isolated DNA (260/280 ratio) and its concentration was measured using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The DNA concentration was ≥1.60. DNA libraries were prepared using the Nextera XT DNA Library Preparation Kit (Illumina) and Nextera Index Kit (Illumina), with a total DNA input of 1 ng. Genomic DNA was fragmented using a proportional amount of Illumina Nextera XT fragmentation enzyme. Combinatory dual indices were added to each sample, followed by 12 cycles of PCR to construct the libraries. The DNA libraries were purified using AMpure magnetic beads (Beckman Coulter) and eluted in QIAGEN EB buffer. The samples were sequenced on an Illumina HiSeq 4000, 2 × 150 bp. DNA extractions were sent to CosmosID (Rockville, MD, USA) to identify the gut microbiota composition at the level of the main microbial phyla by identifying the total bacterial DNA and the Firmicutes, Bacteroidetes, Actinobacteria, Proteobacteria, Verrucomicrobia, and Fusobacteria DNA using the WGS metagenomic sequencing technique.

2.5.1. Bioinformatics Analysis Methods

Unassembled sequencing reads were directly analyzed using the CosmosID bioinformatics platform (CosmosID Inc., Rockville, MD, USA), described elsewhere [34,35,36], for multi-kingdom microbiome analysis, profiling of antibiotic resistance and virulence genes, and quantification of the relative abundance of organisms. The system utilizes a high-performance data-mining k-mer algorithm that rapidly disambiguates millions of short-sequence reads into the discrete genomes engendering the particular sequences. The pipeline has two separable comparators; the first consists of a pre-computation phase for reference databases and the second is a per-sample computation. The input for the pre-computation phase consists of databases of reference genomes, virulence markers, and antimicrobial resistance markers that are continuously curated by CosmosID scientists. The output of the pre-computational phase is a phylogeny tree of microbes together with sets of variable length k-mer fingerprints (biomarkers) uniquely associated with distinct branches and leaves of the tree.
The second per-sample computational phase searches the hundreds of millions of short-sequence reads, or alternatively contigs from draft de novo assemblies, against the fingerprint sets. This query enables the sensitive yet highly precise detection and taxonomic classification of microbial NGS reads. The resulting statistics are analyzed to return the fine-grain taxonomic and relative abundance estimates for the microbial NGS datasets. To exclude false positive identifications, the results are filtered using a filtering threshold based on internal statistical scores that are determined by analyzing a large number of diverse metagenomes. The same approach is applied to enable the sensitive and accurate detection of genetic markers for virulence and for resistance to antibiotics.

2.5.2. Alpha Diversity Boxplots (with Wilcoxon Rank-Sum)

Alpha diversity boxplots were calculated from the phylum, genus, species, and strain-level abundance score matrices from the CosmosID-HUB analysis. Chao, Simpson, and Shannon alpha diversity metrics were calculated in R using the R package vegan [37,38]. Wilcoxon Rank-Sum tests were performed between groups using the R package ggsignif [39]. Boxplots with overlaid significance in p-value format were generated using the R package ggpubr [40].

2.5.3. Beta Diversity PCoA (with PERMANOVA)

Beta diversity principal coordinate analyses were calculated from phylum, genus, species, and strain-level matrices for bacteria from CosmosID-HUB. Bray–Curtis dissimilarity was calculated in R using the vegan package with the function vegdist, and PCoA tables were generated using ape’s function pcoa [41]. PERMANOVA tests for each distance matrix were generated using vegan’s [38] function adonis2, and beta dispersion was calculated and compared using the anova method for the betadisper function from vegan [38]. Plots were visualized using the R package ggpubr [40].

2.6. Statistical Analysis

Quantitative variables were tested for normality prior to the analysis. The normality of the variables was assessed visually using a histogram and Q-Q plot and/or by evaluating skewness and kurtosis. The independent samples t-test was used for continuous variables and outcomes. Nonparametric tests were used for variables that were not normally distributed and that skewed toward one side. To identify the impact of meat by type, we divided the participants into distinct categories, i.e., high and low intake, using the median of each type of meat as the cut-off point. For the intake of white meat, the non-obese group was considered to have high white meat (HWM, n = 32) if the intake was >34 g/1000 kcal and low white meat (LWM, n = 16) if the intake was ≤34 g/1000 kcal. The obese group was stratified into high (HWM, n = 28) if the intake was >25 g/1000 kcal and low (LWM, n = 16) if the intake was ≤25 g/1000 kcal.
The non-obese group was also stratified according to red meat intake into high red meat intake (HRM, n = 21) if intake was >12 g /1000 kcal and low red meat intake (LRM, n = 27) if intake was ≤12 g/1000 kcal. The obese group was considered to have high red meat intake (HRM, n = 14) if intake was >16 g/1000 kcal and low red meat intake (LRM, n = 30) if the red meat intake was ≤16 g/1000 kcal.
A Pearson correlation coefficient was used to identify the correlations between meat (total and by type) and gut flora. All non-normal variables were transformed prior to parametric testing. A p-value of <0.05 and a 95% CI were used to report estimated statistical significance. Furthermore, a Benjamini-Hochberg critical value for a false discovery rate of 0.25 was computed for all correlations. The analysis was carried out using IBM SPSS Statistics for Windows (version 24, IBM Corp., Armonk, NY, USA).
To identify the distribution of species abundances and similarities between groups, the CosmosID application was used to examine the relationship between gut flora and different levels and types of meat intake in the obese and non-obese groups, and by body fat percentage and waist to hip ratio stratification. The Shannon test was used for alpha diversity, which represents the distribution of species abundances in a given sample as a number that depends on species evenness and richness, concentrating on community variation within a single sample [42]. For beta diversity, Bray–Curtis was used to measure the similarity or dissimilarity between samples [42].

3. Results

3.1. Characteristics of Participants

Blood biochemical data, including total cholesterol (TC), HDL-C, LDL-C, total cholesterol/HDL ratio, triglyceride (TG), and hs-CRP, were significantly higher in the obese group compared with the non-obese group (p < 0.05) (Table 1). Additionally, total red meat intake (g/1000 kcal) was significantly higher in the obese group compared with the non-obese group (p < 0.05) (Table 1).

3.2. Correlation between Gut Flora and Meat Intake

Among the total participants, there were positive correlations between Actinobacteria, Bacteroides (unidentified species) (p ≤ 0.05), Flavonifractor plautii (p ≤ 0.0001), and total red meat intake (Table 2). For total white meat intake, there were positive correlations between Bacteroides (unidentified species) and Faecalibacterium Prausnitzii (p ≤ 0.05).

3.2.1. White Meat Intake and Gut Flora

In the non-obese group, there was a positive correlation between LWM intake, Actinobacteria (p = 0.05), and Bifidobacterium longum (p = 0.02) (Table 3). There was also a positive correlation between HWM intake and Bacteroidetes (p = 0.03), Flavonifractor plautii (p ≤ 0.001), and Clostridium Bolteae (p ≤ 0.001), and an inverse correlation with Firmicutes (p = 0.04).
In the obese group, there was a positive correlation between HWM intake and Bacteroides (unidentified species) (p ≤ 0.001), Faecalibacterium Prausnitzii (p = 0.04), and Clostridium difficile (p = 0.03) (Table 3).

3.2.2. Red Meat Intake and Gut Flora

In the non-obese group, HRM intake correlated positively with Flavonifractor plautii (p ≤ 0.001) and inversely with Akkermansia muciniphila (p = 0.05) (Table 4).
In the obese group, there were no significant correlations between HRM intake and gut flora (p > 0.05) (Table 4).

3.2.3. Correlation between Gut Flora and Hs-CRP and Lipid Profile

In the total sample, hs-CRP was inversely correlated with Bifidobacterium adolescentis (p = 0.04) (Table S1). Lipid profile, including TC, HDL-C, LDL-C, TC/HDL ratio, and TG did not show significant correlations with gut bacteria.

3.2.4. Meat Intake and Gut Flora (Alpha and Beta Diversity)

Gut flora diversity results showed no significant differences in alpha diversity (p = 0.09) between total meat intake in the obese and non-obese groups (Figure S1). However, there was a significant difference in beta diversity (p = 0.05) between total meat intake in the obese and non-obese groups (Figure 1). Moreover, there were no significant differences in Alpha diversity (p > 0.05) between high and low white meat intake in the obese and non-obese groups (Figure S2). There were no significant differences in Alpha diversity (p > 0.05) between high and low red meat intake in the obese and non-obese groups (Figure S3).

4. Discussion

4.1. Main Findings

In this study, we found significant correlations between types of meat intake and certain gut bacteria among the total participants, independent of adiposity. We also found significant correlations between types of meat intake and specific gut bacteria which differed among obese and non-obese groups. There was a significant difference in beta diversity between total meat intake in the obese and non-obese groups, but not in alpha diversity. Further, we found that the hs-CRP level correlated negatively with Bifidobacterium adolescentis. No significant relationship was observed between lipid profile and gut flora.

4.2. Comparison with Previous Studies

Our findings show that Actinobacteria, Bacteroides (unidentified species), and Flavonifractor plautii correlated positively with total red meat intake, while both Bacteroides (unidentified species) and Faecalibacterium Prausnitzii correlated positively with total white meat intake among the total participants, indicating that these associations are independent of adiposity.
A previous cross-sectional study that included 98 healthy volunteers aged between 2 and 50 years showed that Bacteroides correlated directly with animal protein, implying that meat consumption is characterized by this enterotype. Furthermore, Bacteroidetes and Actinobacteria phyla were positively associated with fat intake, whereas Firmicutes and Proteobacteria showed inverse associations [20].
To get a clearer understanding of the relation between gut flora and different types of meat relative to obesity, we stratified the participants by both adiposity level and type of meat. In the non-obese group, Bifidobacterium longum was positively correlated with LWM. These findings are in agreement with a study by Reyeset et al., which reported that an abundance of Bifidobacterium longum was positively correlated with unsaturated fat intake, and that it was more abundant in the lean group than in the overweight and obese groups [43]. We also found that in the non-obese group, Flavonifractor plautii was positively correlated with both HWM and HRM intakes. In contrast, a previous RCT that determined the influence of fried meat intake on gut microbiota showed that fried meat intake reduced microbial richness and Flavonifractor abundance [44]. However, the cooking method, such as frying or grilling, may itself have influenced the level of this bacteria [45,46]. We also observed that for the non-obese group there was an inverse correlation between HRM intake and Akkermansia muciniphila. A previous animal study reported differences in the responses of Akkermansia muciniphila to the intake of two different types of protein (chicken protein and soy protein) [47], with the chicken protein-based diet maintaining the Akkermansia muciniphila level in the gut and the soy protein-based diet decreasing the abundance of this bacteria [47]. Nevertheless, several human and animal studies have reported that many types of food and medications may affect Akkermansia muciniphila [48]. Thus, more studies are needed to identify the mechanism behind the effect of red meat on Akkermansia muciniphila while considering other factors that might influence this bacterium. Furthermore, in the non-obese group, there was a positive correlation between LWM intake and Actinobacteria, in contrast to an animal study in which rats fed a high amount of chicken protein had a greater abundance of Actinobacteria than other protein groups [49].
In the obese group, we found that HWM intake was positively correlated with Clostridioides difficile, a finding consistent with a study by Heise et al. (2021) in which 364 different fresh poultry products were screened and 15.80% of samples were shown to have positive levels of Clostridioides difficile [50]. Further, we observed that in the obese group HWM intake was directly correlated with Faecalibacterium prausnitzii. However, a previous prospective observational study showed that following a ketogenic diet that included protein from animal sources such as meat and poultry for three months did not alter Faecalibacterium prausnitzii levels [51]. Faecalibacterium prausnitzii is one of the main butyrate producers in the intestine and it plays a crucial role in gut physiology and host wellbeing [52]. It is the main energy source for colonocytes and confers protective properties against colorectal cancer and inflammatory bowel diseases. However, Faecalibacterium prausnitzii levels have depleted considerably over the past few years [52]. Further investigations into which gut factors modulate its presence are warranted.
From a gut diversity perspective, our study showed that significant differences in gut flora diversity were found only in beta diversity between varying amounts and types of meats. Alpha and beta diversity were measured to observe the richness and variability of gut microbiota composition, respectively [42]. The results of the current study are consistent with a previous study by David et al. (2014), which revealed that no significant differences in alpha diversity were detected between plant- or animal-based diets. However, there was a significant difference in beta diversity that was specific to animal-based diets [17]. Nevertheless, several other factors may influence gut microbes more than diet effects [53]. For example, two large-scale observational population studies recognized 69 and 126 factors related to inter-individual and health traits associated with the gut microbiota [53,54]. Obesity was determined to be one of the important influential factors that affect gut microbiota [55], and low gut microbiota diversity and richness were correlated with an increased risk of obesity [56].
Regarding hs-CRP and gut flora, our results showed a significant inverse correlation between the hs-CRP level and Bifidobacterium adolescentis. A study by Rajkumar et al. (2014) showed that participants with a high hs-CPR level (>3 mg/L) had significantly lower bifidobacteria and lactobacilli compared with those who had a low hs-CRP level (<3 mg/L) [57]. Moreover, a systematic review of 14 observational studies that investigated the association between gut microbiota and hs-CPR level showed that inflammatory markers such as hs-CRP and IL-6 were inversely correlated with Bifidobacterium, Faecalibacterium, Ruminococcus, and Prevotella [58]. The inverse correlation between the hs-CRP level and Bifidobacterium adolescentis may be explained by the association of some Bifidobacterium taxa with reduced systemic inflammation [59,60]. Further, our study showed no significant correlation between lipid profile and gut flora in contrast with other studies that showed significant correlations between different gut flora and lipid profiles [61,62], which may be because specific operational taxonomic units reported to be associated with lipid profiles did not exist in the participants in our study. Moreover, other factors may influence the lipid profile besides the gut flora, such as age, sex, and genetics [61].
The main limitation of this study is the generalizability of these results, as the study population included only females with a specific age range. Moreover, this study is based on observations, and therefore causality is not possible. Furthermore, even with the extensive measures used in the data collection process, there was still a possibility of recall bias, instrument error, and confirmation bias. One of the study’s strengths is its use of the WGS technique, the gold standard method for identifying the gut flora composition, showing the highest sensitivity compared with other techniques [63]. Moreover, the design of the study allowed us to identify the differences in the intake of various amounts and types of meat in relation to gut flora. Additionally, to our knowledge, our work is the first to examine associations between types of meat intake, gut flora, lipid profile, and hs-CRP in females in the Middle Eastern region. Furthermore, precautions such as repeated measurements, trained staff, and interview-based multi-pass dietary assessment methods were used to increase the accuracy and precision of the data collection process.

5. Conclusions

In conclusion, our findings suggest that certain bacteria are correlated with meat intake (by type of meat), independent of adiposity. We also found significant correlations between types of meat intake and specific gut bacteria which differed among obese and non-obese groups. Further investigations are needed to explore the correlations we reported, which are not explained in the literature, between meat intake and different gut bacteria, including Flavonifractorplautii, Actinobacteria, Bifidobacteriumlongum, Clostridiumbolteae, Firmicutes, Faecalibacterium prausnitzii, Clostridioides difficile, Akkermansia muciniphila, and Bacteroides with unidentified species. Moreover, the relationship between gut flora and meat intake can differ in beta diversity between obese and non-obese groups, but not in alpha diversity. Future studies with a larger number of participants of different ages, from both genders, and from various regions are recommended.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods12020245/s1, Figure S1: Total meat intake across obese and non-obese groups (alpha diversity-Shanon); Figure S2: High and low white meat intake across cases and controls (Alpha diversity- Shannon); Figure S3: High and low red meat intake across cases and controls (Alpha diversity- Shannon). Table S1: Correlation between gut flora and lipid profile.

Author Contributions

Conceptualization, S.A.-M., S.S., G.S.A. and E.A.A.; methodology, M.A. and J.A.; software, M.A.; validation, S.S., M.A., G.S.A. and E.A.A.; formal analysis, E.A.A. and J.A.; investigation, J.A.; resources, G.S.A.; data curation, G.S.A.; writing—original draft preparation, E.A.A., S.A.-M. and G.S.A.; writing—review and editing, G.S.A., S.A.-M. and S.S.; visualization, S.A.-M.; supervision, G.S.A.; funding acquisition, S.A.-M. and G.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project No. IFKSURG-2-926.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of King Khalid University Hospital at King Saud University (IRB #E-19-3625) on 19 February 2019.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Popkin, B.M. Global nutrition dynamics: The world is shifting rapidly toward a diet linked with noncommunicable diseases. Am. J. Clin. Nutr. 2006, 84, 289–298. [Google Scholar] [CrossRef] [PubMed]
  2. Anand, S.S.; Hawkes, C.; De Souza, R.J.; Mente, A.; Dehghan, M.; Nugent, R.; Zulyniak, M.A.; Weis, T.; Bernstein, A.M.; Krauss, R.M. Food consumption and its impact on cardiovascular disease: Importance of solutions focused on the globalized food system: A report from the workshop convened by the World Heart Federation. J. Am. Coll. Cardiol. 2015, 66, 1590–1614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. DeNicola, E.; Aburizaiza, O.S.; Siddique, A.; Khwaja, H.; Carpenter, D.O. Obesity and public health in the Kingdom of Saudi Arabia. Rev. Environ. Health 2015, 30, 191–205. [Google Scholar] [CrossRef] [PubMed]
  4. Miller, V.; Reedy, J.; Cudhea, F.; Zhang, J.; Shi, P.; Erndt-Marino, J.; Coates, J.; Micha, R.; Webb, P.; Mozaffarian, D. Global, regional, and national consumption of animal-source foods between 1990 and 2018: Findings from the Global Dietary Database. Lancet Planet. Health 2022, 6, e243–e256. [Google Scholar] [CrossRef] [PubMed]
  5. Intelligence, G.M. Saudi Arabia Red Meat Market Analysis: Market Segmentation By Meat Type, By Storage, By Distributional Channel & By Region With Forecast 2017–2030. Available online: https://www.goldsteinresearch.com/report/saudi-arabia-red-meat-market (accessed on 18 September 2022).
  6. Ministry of Commerce and Industry. Annual Report 2009; European Central Bank: Frankfurt, Germany, 2009. [Google Scholar]
  7. Al-Nuaim, A.A.; Bamgboye, E.A.; Al-Rubeaan, K.A.; Al-Mazrou, Y. Overweight and obesity in Saudi Arabian adult population, role of sociodemographic variables. J. Community Health 1997, 22, 211–223. [Google Scholar] [CrossRef]
  8. Centers for Disease Control and Prevention. Adult Obesity Facts. Available online: https://www.cdc.gov/obesity/data/adult.html#print (accessed on 18 September 2022).
  9. Vergnaud, A.C.; Norat, T.; Romaguera, D.; Mouw, T.; May, A.M.; Travier, N.; Luan, J.; Wareham, N.; Slimani, N.; Rinaldi, S.; et al. Meat consumption and prospective weight change in participants of the EPIC-PANACEA study. Am. J. Clin. Nutr. 2010, 92, 398–407. [Google Scholar] [CrossRef] [Green Version]
  10. You, W.; Henneberg, M. Meat consumption providing a surplus energy in modern diet contributes to obesity prevalence: An ecological analysis. BMC Nutr. 2016, 2, 1–11. [Google Scholar]
  11. Pan, A.; Sun, Q.; Bernstein, A.M.; Schulze, M.B.; Manson, J.E.; Willett, W.C.; Hu, F.B. Red meat consumption and risk of type 2 diabetes: 3 cohorts of US adults and an updated meta-analysis. Am. J. Clin. Nutr. 2011, 94, 1088–1096. [Google Scholar] [CrossRef] [Green Version]
  12. Micha, R.; Michas, G.; Mozaffarian, D. Unprocessed red and processed meats and risk of coronary artery disease and type 2 diabetes–an updated review of the evidence. Curr. Atheroscler. Rep. 2012, 14, 515–524. [Google Scholar] [CrossRef] [Green Version]
  13. Papier, K.; Hartman, L.; Tong, T.Y.N.; Key, T.J.; Knuppel, A. Higher Meat Intake Is Associated with Higher Inflammatory Markers, Mostly Due to Adiposity: Results from UK Biobank. J. Nutr. 2022, 152, 183–189. [Google Scholar] [CrossRef]
  14. Choi, J.; Joseph, L.; Pilote, L. Obesity and C-reactive protein in various populations: A systematic review and meta-analysis. Obes. Rev. 2013, 14, 232–244. [Google Scholar] [CrossRef] [PubMed]
  15. Singh, R.K.; Chang, H.W.; Yan, D.; Lee, K.M.; Ucmak, D.; Wong, K.; Abrouk, M.; Farahnik, B.; Nakamura, M.; Zhu, T.H.; et al. Influence of diet on the gut microbiome and implications for human health. J. Transl. Med. 2017, 15, 73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Zhang, C.; Zhang, M.; Pang, X.; Zhao, Y.; Wang, L.; Zhao, L. Structural resilience of the gut microbiota in adult mice under high-fat dietary perturbations. ISME J. 2012, 6, 1848–1857. [Google Scholar] [CrossRef]
  17. David, L.A.; Maurice, C.F.; Carmody, R.N.; Gootenberg, D.B.; Button, J.E.; Wolfe, B.E.; Ling, A.V.; Devlin, A.S.; Varma, Y.; Fischbach, M.A. Diet rapidly and reproducibly alters the human gut microbiome. Nature 2014, 505, 559–563. [Google Scholar] [CrossRef] [Green Version]
  18. Foerster, J.; Maskarinec, G.; Reichardt, N.; Tett, A.; Narbad, A.; Blaut, M.; Boeing, H. The influence of whole grain products and red meat on intestinal microbiota composition in normal weight adults: A randomized crossover intervention trial. PLoS ONE 2014, 9, e109606. [Google Scholar] [CrossRef] [Green Version]
  19. Lunn, J.; Kuhnle, G.; Mai, V.; Frankenfeld, C.; Shuker, D.; Glen, R.; Goodman, J.; Pollock, J.; Bingham, S. The effect of haem in red and processed meat on the endogenous formation of N-nitroso compounds in the upper gastrointestinal tract. Carcinogenesis 2007, 28, 685–690. [Google Scholar] [CrossRef]
  20. Wu, G.D.; Chen, J.; Hoffmann, C.; Bittinger, K.; Chen, Y.-Y.; Keilbaugh, S.A.; Bewtra, M.; Knights, D.; Walters, W.A.; Knight, R. Linking long-term dietary patterns with gut microbial enterotypes. Science 2011, 334, 105–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Le Leu, R.K.; Winter, J.M.; Christophersen, C.T.; Young, G.P.; Humphreys, K.J.; Hu, Y.; Gratz, S.W.; Miller, R.B.; Topping, D.L.; Bird, A.R. Butyrylated starch intake can prevent red meat-induced O6-methyl-2-deoxyguanosine adducts in human rectal tissue: A randomised clinical trial. Br. J. Nutr. 2015, 114, 220–230. [Google Scholar] [CrossRef] [Green Version]
  22. Yano, J.M.; Yu, K.; Donaldson, G.P.; Shastri, G.G.; Ann, P.; Ma, L.; Nagler, C.R.; Ismagilov, R.F.; Mazmanian, S.K.; Hsiao, E.Y. Indigenous bacteria from the gut microbiota regulate host serotonin biosynthesis. Cell 2015, 161, 264–276. [Google Scholar] [CrossRef] [Green Version]
  23. Borrello, K.; Lim, U.; Park, S.Y.; Monroe, K.R.; Maskarinec, G.; Boushey, C.J.; Wilkens, L.R.; Randolph, T.W.; Le Marchand, L.; Hullar, M.A.; et al. Dietary Intake Mediates Ethnic Differences in Gut Microbial Composition. Nutrients 2022, 14, 660. [Google Scholar] [CrossRef]
  24. Dwiyanto, J.; Hussain, M.H.; Reidpath, D.; Ong, K.S.; Qasim, A.; Lee, S.W.H.; Lee, S.M.; Foo, S.C.; Chong, C.W.; Rahman, S. Ethnicity influences the gut microbiota of individuals sharing a geographical location: A cross-sectional study from a middle-income country. Sci. Rep. 2021, 11, 2618. [Google Scholar] [CrossRef] [PubMed]
  25. Aljazairy, E.A.; Al-Musharaf, S.; Abudawood, M.; Almaarik, B.; Hussain, S.D.; Alnaami, A.M.; Sabico, S.; Al-Daghri, N.M.; Clerici, M.; Aljuraiban, G.S. Influence of Adiposity on the Gut Microbiota Composition of Arab Women: A Case-Control Study. Biology 2022, 11, 1586. [Google Scholar] [CrossRef] [PubMed]
  26. Kasai, C.; Sugimoto, K.; Moritani, I.; Tanaka, J.; Oya, Y.; Inoue, H.; Tameda, M.; Shiraki, K.; Ito, M.; Takei, Y.; et al. Comparison of the gut microbiota composition between obese and non-obese individuals in a Japanese population, as analyzed by terminal restriction fragment length polymorphism and next-generation sequencing. BMC Gastroenterol. 2015, 15, 100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. World Health Organization. Obesity and Overweight. Available online: https://www.who.int/news-room/fact-sheets/detail/obesity-and-overweight (accessed on 19 September 2022).
  28. Nishida, C.; Ko, G.; Kumanyika, S. Body fat distribution and noncommunicable diseases in populations: Overview of the 2008 WHO Expert Consultation on Waist Circumference and Waist–Hip Ratio. Eur. J. Clin. Nutr. 2010, 64, 2–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. So, W.-Y.; Swearingin, B.; Crooms, B.; Lee, R.; Choi, Y.; Dail, T.K.; Ha, C. Body composition measurements determined by air displacement plethysmography and eight-polar bioelectrical impedance analysis are equivalent in African American college students. HealthMED 2012, 6, 1896–1899. [Google Scholar]
  30. Romero-Corral, A.; Somers, V.K.; Sierra-Johnson, J.; Thomas, R.J.; Collazo-Clavell, M.; Korinek, J.; Allison, T.G.; Batsis, J.; Sert-Kuniyoshi, F.; Lopez-Jimenez, F. Accuracy of body mass index in diagnosing obesity in the adult general population. Int. J. Obes. 2008, 32, 959–966. [Google Scholar] [CrossRef] [Green Version]
  31. Alkhalaf, M.; Edwards, C.; Combet, E. Validation of a food frequency questionnaire specific for salt intake in Saudi Arabian adults using urinary biomarker and repeated multiple pass 24-hour dietary recall. Proc. Nutr. Soc. 2015, 74, E337. [Google Scholar] [CrossRef] [Green Version]
  32. Ahmad, S.; Badpa, A.G. Meat products and Byproducts for value Addition. In Food Processing: Strategies for Quality Assessment; Springer: Berlin/Heidelberg, Germany, 2014; pp. 125–153. [Google Scholar]
  33. Friedewald, W.T.; Levy, R.I.; Fredrickson, D.S. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin. Chem. 1972, 18, 499–502. [Google Scholar] [CrossRef]
  34. Ottesen, A.; Ramachandran, P.; Reed, E.; White, J.R.; Hasan, N.; Subramanian, P.; Ryan, G.; Jarvis, K.; Grim, C.; Daquiqan, N.; et al. Enrichment dynamics of Listeria monocytogenes and the associated microbiome from naturally contaminated ice cream linked to a listeriosis outbreak. BMC Microbiol. 2016, 16, 275. [Google Scholar] [CrossRef] [Green Version]
  35. Ponnusamy, D.; Kozlova, E.V.; Sha, J.; Erova, T.E.; Azar, S.R.; Fitts, E.C.; Kirtley, M.L.; Tiner, B.L.; Andersson, J.A.; Grim, C.J. Cross-talk among flesh-eating Aeromonas hydrophila strains in mixed infection leading to necrotizing fasciitis. Proc. Natl. Acad. Sci. USA 2016, 113, 722–727. [Google Scholar] [CrossRef] [Green Version]
  36. Lax, S.; Smith, D.P.; Hampton-Marcell, J.; Owens, S.M.; Handley, K.M.; Scott, N.M.; Gibbons, S.M.; Larsen, P.; Shogan, B.D.; Weiss, S.; et al. Longitudinal analysis of microbial interaction between humans and the indoor environment. Science 2014, 345, 1048–1052. [Google Scholar] [CrossRef] [Green Version]
  37. Jost, L. Entropy and diversity. Oikos 2006, 113, 363–375. [Google Scholar] [CrossRef]
  38. Oksanen, J.; Simpson, G.L.; Guillaume Blanche, F.; Kindt, R.; Legendre, P.; Minchin, P.R.; O'Hara, R.B.; Solymos, P.; Stevens, M.H.H.; Szoecs, E.; et al. Package Vegan. 2022. Available online: https://github.com/vegandevs/vegan (accessed on 18 September 2022).
  39. Ahlmann-Eltze, C. Ggsignif: Significance Brackets for 'ggplot2'. Available online: https://cran.r-project.org/web/packages/ggsignif/index.html (accessed on 18 September 2022).
  40. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  41. Paradis, E.; Schliep, K. ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 2019, 35, 526–528. [Google Scholar] [CrossRef] [PubMed]
  42. Walters, K.; Martiny, J. Alpha-, beta-, and gamma-diversity of bacteria varies across habitats. PLoS ONE 2020, 15, e0233872. [Google Scholar] [CrossRef]
  43. Mayorga Reyes, L.; Gonzalez Vazquez, R.; Cruz Arroyo, S.M.; Melendez Avalos, A.; Reyes Castillo, P.A.; Chavaro Perez, D.A.; Ramos Terrones, I.; Ramos Ibanez, N.; Rodriguez Magallanes, M.M.; Langella, P. Correlation between diet and gut bacteria in a population of young adults. Int. J. Food Sci. Nutr. 2016, 67, 470–478. [Google Scholar] [CrossRef] [PubMed]
  44. Gao, J.; Guo, X.; Wei, W.; Li, R.; Hu, K.; Liu, X.; Jiang, W.; Liu, S.; Wang, W.; Sun, H. The association of fried meat consumption with the gut microbiota and fecal metabolites and its impact on glucose homoeostasis, intestinal endotoxin levels, and systemic inflammation: A randomized controlled-feeding trial. Diabetes Care 2021, 44, 1970–1979. [Google Scholar] [CrossRef]
  45. Pérez-Burillo, S.; Pastoriza, S.; Jiménez-Hernandez, N.; D’Auria, G.; Francino, M.P.; Rufián-Henares, J.A. Effect of food thermal processing on the composition of the gut microbiota. J. Agric. Food Chem. 2018, 66, 11500–11509. [Google Scholar] [CrossRef]
  46. Shen, Q.; Chen, Y.A.; Tuohy, K.M. A comparative in vitro investigation into the effects of cooked meats on the human faecal microbiota. Anaerobe 2010, 16, 572–577. [Google Scholar] [CrossRef] [PubMed]
  47. Zhao, F.; Zhou, G.; Liu, X.; Song, S.; Xu, X.; Hooiveld, G.; Müller, M.; Liu, L.; Kristiansen, K.; Li, C. Dietary protein sources differentially affect the growth of Akkermansia muciniphila and maintenance of the gut mucus barrier in mice. Mol. Nutr. Food Res. 2019, 63, 1900589. [Google Scholar] [CrossRef]
  48. Hagi, T.; Belzer, C. The interaction of Akkermansia muciniphila with host-derived substances, bacteria and diets. Appl. Microbiol. Biotechnol. 2021, 105, 4833–4841. [Google Scholar] [CrossRef]
  49. Zhu, Y.; Lin, X.; Zhao, F.; Shi, X.; Li, H.; Li, Y.; Zhu, W.; Xu, X.; Li, C.; Zhou, G. Meat, dairy and plant proteins alter bacterial composition of rat gut bacteria. Sci. Rep. 2015, 5, 15220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Heise, J.; Witt, P.; Maneck, C.; Wichmann-Schauer, H.; Maurischat, S. Prevalence and phylogenetic relationship of Clostridioides difficile strains in fresh poultry meat samples processed in different cutting plants. Int. J. Food Microbiol. 2021, 339, 109032. [Google Scholar] [CrossRef] [PubMed]
  51. Tagliabue, A.; Ferraris, C.; Uggeri, F.; Trentani, C.; Bertoli, S.; de Giorgis, V.; Veggiotti, P.; Elli, M. Short-term impact of a classical ketogenic diet on gut microbiota in GLUT1 Deficiency Syndrome: A 3-month prospective observational study. Clin. Nutr. ESPEN 2017, 17, 33–37. [Google Scholar] [CrossRef]
  52. Lopez-Siles, M.; Duncan, S.H.; Garcia-Gil, L.J.; Martinez-Medina, M. Faecalibacterium prausnitzii: From microbiology to diagnostics and prognostics. ISME J. 2017, 11, 841–852. [Google Scholar] [CrossRef]
  53. Falony, G.; Joossens, M.; Vieira-Silva, S.; Wang, J.; Darzi, Y.; Faust, K.; Kurilshikov, A.; Bonder, M.J.; Valles-Colomer, M.; Vandeputte, D. Population-level analysis of gut microbiome variation. Science 2016, 352, 560–564. [Google Scholar] [CrossRef]
  54. Zhernakova, A.; Kurilshikov, A.; Bonder, M.J.; Tigchelaar, E.F.; Schirmer, M.; Vatanen, T.; Mujagic, Z.; Vila, A.V.; Falony, G.; Vieira-Silva, S. Population-based metagenomics analysis reveals markers for gut microbiome composition and diversity. Science 2016, 352, 565–569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Castaner, O.; Goday, A.; Park, Y.-M.; Lee, S.-H.; Magkos, F.; Shiow, S.-A.T.E.; Schröder, H. The gut microbiome profile in obesity: A systematic review. Int. J. Endocrinol. 2018, 2018, 4095789. [Google Scholar] [CrossRef] [PubMed]
  56. Vallianou, N.; Stratigou, T.; Christodoulatos, G.S.; Dalamaga, M. Understanding the role of the gut microbiome and microbial metabolites in obesity and obesity-associated metabolic disorders: Current evidence and perspectives. Curr. Obes. Rep. 2019, 8, 317–332. [Google Scholar] [CrossRef]
  57. Rajkumar, H.; Mahmood, N.; Kumar, M.; Varikuti, S.R.; Challa, H.R.; Myakala, S.P. Effect of probiotic (VSL# 3) and omega-3 on lipid profile, insulin sensitivity, inflammatory markers, and gut colonization in overweight adults: A randomized, controlled trial. Mediat. Inflamm. 2014, 2014, 348959. [Google Scholar]
  58. van den Munckhof, I.C.; Kurilshikov, A.; ter Horst, R.; Riksen, N.P.; Joosten, L.; Zhernakova, A.; Fu, J.; Keating, S.T.; Netea, M.G.; De Graaf, J. Role of gut microbiota in chronic low-grade inflammation as potential driver for atherosclerotic cardiovascular disease: A systematic review of human studies. Obes. Rev. 2018, 19, 1719–1734. [Google Scholar] [CrossRef] [Green Version]
  59. Miller, L.E.; Lehtoranta, L.; Lehtinen, M.J. The Effect of Bifidobacterium animalis ssp. lactis HN019 on Cellular Immune Function in Healthy Elderly Subjects: Systematic Review and Meta-Analysis. Nutrients 2017, 9, 191. [Google Scholar] [CrossRef] [PubMed]
  60. Cani, P.D.; Neyrinck, A.M.; Fava, F.; Knauf, C.; Burcelin, R.G.; Tuohy, K.M.; Gibson, G.; Delzenne, N.M. Selective increases of bifidobacteria in gut microflora improve high-fat-diet-induced diabetes in mice through a mechanism associated with endotoxaemia. Diabetologia 2007, 50, 2374–2383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Fu, J.; Bonder, M.J.; Cenit, M.C.; Tigchelaar, E.F.; Maatman, A.; Dekens, J.A.; Brandsma, E.; Marczynska, J.; Imhann, F.; Weersma, R.K. The gut microbiome contributes to a substantial proportion of the variation in blood lipids. Circ. Res. 2015, 117, 817–824. [Google Scholar] [CrossRef] [PubMed]
  62. Vazquez-Moreno, M.; Perez-Herrera, A.; Locia-Morales, D.; Dizzel, S.; Meyre, D.; Stearns, J.C.; Cruz, M. Association of gut microbiome with fasting triglycerides, fasting insulin and obesity status in Mexican children. Pediatr. Obes. 2021, 16, e12748. [Google Scholar] [CrossRef] [PubMed]
  63. Poretsky, R.; Rodriguez-R, L.M.; Luo, C.; Tsementzi, D.; Konstantinidis, K.T. Strengths and limitations of 16S rRNA gene amplicon sequencing in revealing temporal microbial community dynamics. PLoS ONE 2014, 9, e93827. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Total meat intake across obese and non-obese groups (beta diversity, Bray–Curtis).
Figure 1. Total meat intake across obese and non-obese groups (beta diversity, Bray–Curtis).
Foods 12 00245 g001
Table 1. Characteristics of obese and non-obese groups, n = 92 1.
Table 1. Characteristics of obese and non-obese groups, n = 92 1.
CharacteristicsTotal
(n = 92)
Non-Obese
(BMI 18.50–24.99)
(n = 48)
Obese
(BMI ≥ 30)
(n = 44)
p-Value
Age (years) (mean ± SD)21.10 ± 1.5020.60 ± 1.1021.60 ± 1.70<0.001
BMI (kg/m2)28.50 ± 8.0021.70 ± 1.9036.00 ± 4.70<0.001
Waist to hip ratio (cm)0.70 ± 0.100.70 ± 0.100.80 ± 0.100.01
Fat (kg)42.50 ± 9.4034.80 ± 5.5051.10 ± 3.30<0.001
Skeletal muscle mass (kg)21.00 ± 3.5018.60 ± 2.2023.60 ± 2.60<0.001
Total white meat (g/1000 kcal)32 (19–41)34 (26–41)25 (16–40)0.11
Total red meat (g/1000 kcal)15 (9–19)12 (4–16)16 (13–21)<0.001
LWM (%)3533360.01
HWM (%)656764
LRM (%)625668
HRM (%)384432
Biochemical Measurements
Total Cholesterol (mmol/L)4.10 ± 1.503.60 ± 1.704.50 ± 1.000.01
HDL-Cholesterol (mmol/L)1.00 ± 0.300.90 ± 0.401.00 ± 0.300.24
LDL-Cholesterol (mmol/L)2.90 ± 1.302.60 ± 1.503.30 ± 1.000.01
Total cholesterol/HDL ratio4.30 ± 1.703.90 ± 1.704.70 ± 1.700.05
Triglyceride (mmol/L)0.70 (0.50–1.00)0.50 (0.40–0.70)1.00 (0.80–1.10)<0.001
High-sensitivity C-reactive Protein (ng/mL)1473.90 (787.10–7635.60)1031.20 (553–1514)6385.40 (1235–12,988)<0.001
1 Data are presented as mean ± SD for normal variables, median (1st quartile–3rd quartile) for non-normal variables, and N (%) for categorical variables. BMI: body mass index; HDL: high-density lipoprotein; HRM: high red meat intake; HWM: high white meat intake; LDL: low-density lipoprotein; LRM: low red meat intake; LWM: low white meat intake.
Table 2. Correlations between gut flora by type of meat intake among total participants 1.
Table 2. Correlations between gut flora by type of meat intake among total participants 1.
Gut FloraTotal White Meat (g/1000 kcal)Total Red Meat (g/1000 kcal)
Rp(i/m) QRp(i/m) Q
Flavonifractor plautii−0.040.700.180.31<0.00010.01
Bacteroides (unidentified species)0.210.050.010.230.040.03
Actinobacteria−0.110.310.060.220.050.04
Faecalibacterium Prausnitzii0.220.050.030.180.070.06
Bifidobacterium adolescentis−0.090.420.120.150.150.07
Fusobacteria0.060.610.16−0.110.320.09
Bifidobacterium longum−0.110.310.070.090.410.10
Proteobacteria−0.020.870.210.080.460.12
Clostridium Bolteae0.080.450.150.080.470.13
Firmicutes−0.080.450.13−0.070.490.15
Verrucomicrobia0.010.900.24−0.040.740.16
Akkermansia muciniphila0.010.890.22−0.030.750.18
F:B ratio−0.010.960.25−0.020.830.19
Bacteroides uniformis−0.150.150.04−0.020.850.21
Clostridium difficile0.040.710.190.020.860.22
Bacteria (unidentified phylum)−0.10.360.090.010.940.24
Bacteroidetes0.110.370.100.010.950.25
1 Data are presented as Pearson correlation coefficient adjusted for age. F:B: Firmicutes to Bacteroidetes; (i/m) Q: Benjamini-Hochberg critical value for a false discovery rate of 0.25.
Table 3. Correlations between gut flora and white meat intake in obese and non-obese groups 1.
Table 3. Correlations between gut flora and white meat intake in obese and non-obese groups 1.
Gut FloraNon-ObeseObese
LWMHWMLWMHWM
Rp(i/m) QRp(i/m) QRp(i/m) QRp(i/m) Q
Actinobacteria0.420.050.03−0.210.340.130.060.760.18−0.220.450.12
Bifidobacterium longum0.480.020.020.010.960.25−0.040.840.22−0.280.330.09
Bacteroidetes−0.230.290.110.450.030.05−0.130.530.100.030.910.25
Bacteria (unidentified phylum)−0.090.690.23−0.200.370.140.170.400.030.430.120.06
F:B ratio0.160.460.20−0.250.260.090.120.560.120.060.840.19
Clostridium Bolteae0.250.260.090.60<0.0010.030.040.860.24−0.060.850.21
Firmicutes0.090.670.22−0.430.040.060.120.570.150.050.870.24
Flavonifractor plautii0.270.210.080.60<0.0010.020.020.900.25−0.270.350.10
Clostridium difficile0.230.290.130.030.900.23−0.050.800.210.670.030.03
Bacteroides uniformis0.180.410.190.150.490.20−0.050.790.19−0.430.120.07
Bacteroides (unidentified species)−0.370.090.050.120.580.220.210.290.010.85<0.0010.01
Faecalibacterium Prausnitzii0.190.380.14−0.150.490.17−0.100.630.160.550.040.04
Bifidobacterium adolescentis0.290.170.06−0.230.280.110.120.560.13−0.090.770.16
Proteobacteria0.070.760.25−0.330.130.080.160.430.07−0.150.610.13
Fusobacteria 0.130.510.09−0.050.860.22
Verrucomicrobia0.180.410.16−0.150.490.19−0.170.400.060.080.780.18
Akkermansia muciniphila0.180.410.17−0.180.410.16−0.170.400.040.090.770.15
1 Data are presented as Pearson correlation coefficients adjusted for age. F:B: Firmicutes to Bacteroidetes; HWM: high white meat intake; LWM: low white meat intake; (i/m) Q: Benjamini-Hochberg critical value for a false discovery rate of 0.25. High white meat intake is >25 g/1000 kcal for the obese group, low white meat intake is ≤25 g/1000 kcal for the obese group, high white meat intake is >34 g/1000 kcal for the non-obese group, and low white meat intake is ≤34 g/1000 kcal for the non-obese group. The hyphens in cells indicate that gut flora is not present in these categories or that the number of observations for this group is not high enough to report a correlation coefficient.
Table 4. Correlations between gut flora and red meat intake in obese and non-obese groups 1.
Table 4. Correlations between gut flora and red meat intake in obese and non-obese groups 1.
Gut FloraNon-ObeseObese
LRMHRMLRMHRM
Rp(i/m) QRp(i/m) QRp(i/m) QRp(i/m) Q
Actinobacteria−0.160.480.050.010.950.230.270.140.040.550.100.02
Bifidobacterium longum−0.130.570.06−0.160.440.130.300.100.010.460.180.04
Bacteroidetes0.040.850.190.050.830.190.060.740.19−0.350.320.05
Bacteria (unidentified phylum)0.000.990.230.210.320.080.100.590.180.340.330.07
F:B ratio−0.020.940.22−0.010.980.25−0.150.410.150.310.380.09
Clostridium Bolteae−0.040.870.200.210.320.06−0.140.450.16−0.290.410.11
Firmicutes−0.070.770.17−0.050.820.17−0.160.380.120.200.590.13
Flavonifractor plautii−0.080.730.160.740.000.02−0.060.750.21−0.170.640.14
Clostridium difficile−0.080.720.14−0.190.370.090.010.950.25−0.170.640.16
Bacteroides uniformis0.010.990.25−0.190.380.11−0.040.820.220.140.690.18
Bacteroides (unidentified species)−0.200.380.030.070.740.16−0.040.840.24−0.130.720.20
Faecalibacterium Prausnitzii−0.090.690.13−0.150.490.140.230.200.070.120.750.21
Bifidobacterium adolescentis−0.130.570.080.030.900.220.260.160.060.070.850.23
Proteobacteria0.280.200.02−0.040.840.200.160.400.130.020.960.25
Fusobacteria-- -- −0.300.100.03---
Verrucomicrobia0.090.680.11−0.370.080.05−0.180.340.09---
Akkermansia muciniphila0.090.680.09−0.400.050.03−0.170.370.10---
1 Data are presented as Pearson correlation coefficients adjusted for age. F:B: Firmicutes to Bacteroidetes; HRM: high red meat intake; (i/m) Q: Benjamini-Hochberg critical value for a false discovery rate of 0.25; LRM: low red meat intake. High red meat intake is > 11 g/1000 kcal for the obese group, low red meat intake is ≤ 11 g/1000 kcal for the obese group, high red meat intake is > 16 g/1000 kcal for the non-obese group, and low red meat intake is ≤ 16 g/1000 kcal for the non-obese group. The hyphens in cells indicate that gut flora is not present in these categories or that the number of observations for this group is not high enough to report a correlation coefficient.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Almajed, J.; Al-Musharaf, S.; Abudawood, M.; Sabico, S.; Aljazairy, E.A.; Aljuraiban, G.S. Red and White Meat Intake in Relation to Gut Flora in Obese and Non-Obese Arab Females. Foods 2023, 12, 245. https://doi.org/10.3390/foods12020245

AMA Style

Almajed J, Al-Musharaf S, Abudawood M, Sabico S, Aljazairy EA, Aljuraiban GS. Red and White Meat Intake in Relation to Gut Flora in Obese and Non-Obese Arab Females. Foods. 2023; 12(2):245. https://doi.org/10.3390/foods12020245

Chicago/Turabian Style

Almajed, Jinan, Sara Al-Musharaf, Manal Abudawood, Shaun Sabico, Esra’a A. Aljazairy, and Ghadeer S. Aljuraiban. 2023. "Red and White Meat Intake in Relation to Gut Flora in Obese and Non-Obese Arab Females" Foods 12, no. 2: 245. https://doi.org/10.3390/foods12020245

APA Style

Almajed, J., Al-Musharaf, S., Abudawood, M., Sabico, S., Aljazairy, E. A., & Aljuraiban, G. S. (2023). Red and White Meat Intake in Relation to Gut Flora in Obese and Non-Obese Arab Females. Foods, 12(2), 245. https://doi.org/10.3390/foods12020245

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop