# Impact of smoking cessation, coffee and bread consumption on the intestinal microbial composition among Saudis: A cross-sectional study

Steve Harakeh, Emmanouil Angelakis, Timokratis Karamitros, Dipankar Bachar, Suhad Bahijri, Ghada Ajabnoor, Sulaiman M. Alfadul, Suha A. Farraj, Turki Al Amri, Ahmed Al-Hejin, et al.

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

## ► To cite this version:

Steve Harakeh, Emmanouil Angelakis, Timokratis Karamitros, Dipankar Bachar, Suhad Bahijri, et al.. Impact of smoking cessation, coffee and bread consumption on the intestinal microbial composition among Saudis: A cross-sectional study. PLoS ONE, 2020, 15 (4), pp.e0230895. 10.1371/journal.pone.0230895. hal-03150169

HAL Id: hal-03150169

https://amu.hal.science/hal-03150169

Submitted on 3 Sep 2021

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

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

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

# Impact of smoking cessation, coffee and bread consumption on the intestinal microbial composition among Saudis: A cross-sectional study

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

Steve Harakeh<sup>1*</sup>, Emmanouil Angelakis<sup>2,3</sup>, Timokratis Karamitros<sup>5</sup>, Dipankar Bachar<sup>2</sup>, Suhad Bahijri<sup>5</sup>, Ghada Ajabnoor<sup>5</sup>, Sulaiman M. Alfadul<sup>6</sup>, Suha A. Farraj<sup>1</sup>, Turki Al Amri<sup>7</sup>, Ahmed Al-Hejin<sup>8</sup>, Abdalla Ahmed<sup>9</sup>, Ahmed A. Mirza<sup>1,10</sup>, Raoult Didier<sup>11</sup>, Esam I. Azhar<sup>1,10</sup>

1 Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia, 2 Aix Marseille Université, IRD, APHM, VITROME, IHU-Méditerranée Infection, Marseille, France, 3 Laboratory of Medical Microbiology, Hellenic Pasteur Institute, Athens, Greece, 4 Unit of Bioinformatics and Applied Genomics, Hellenic Pasteur Institute, Athens, Greece, 5 Clinical Biochemistry Department, College of Medicine, Nutrition unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia, 6 King Abdulaziz city for Science and Technology, Riyadh, Saudi Arabia, 7 Family and Community Medicine department, Faculty of Medicine-Rabigh Branch, King Abdulaziz University, Jeddah, Saudi Arabia, 8 Biological Sciences Department, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia, 9 Department of Medical Microbiology, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia, 10 Department of Laboratory Medical Technology, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah, Saudi Arabia, 11 Aix Marseille Université, IRD, APHM, MEPHI, IHU-Méditerranée Infection, Marseille, France

* sharakeh@gmail.com

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

OPEN ACCESS

**Citation:** Harakeh S, Angelakis E, Karamitros T, Bachar D, Bahijri S, Ajabnoor G, et al. (2020) Impact of smoking cessation, coffee and bread consumption on the intestinal microbial composition among Saudis: A cross-sectional study. PLoS ONE 15(4): e0230895. https://doi.org/10.1371/journal.pone.0230895

**Editor:** Erwin G Zoetendal, Wageningen Universiteit, NETHERLANDS

**Received:** October 12, 2018

**Accepted:** March 12, 2020

**Published:** April 29, 2020

**Copyright:** © 2020 Harakeh et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The data generated and analyzed during this study are available in DRYAD Digital Repository: https://doi.org/10.5061/dryad.qv9s4mwbc.

**Funding:** The authors express their sincere gratitude and gratefulness to King Abdulaziz City for Science and Technology, Saudi Arabia, for funding this work under Grant #AR-34-191 to SH.

**Competing interests:** The authors declared that they have no competing interests.

## Abstract

The gut microbiota is often affected by the dietary and lifestyle habits of the host, resulting in a better efficacy that favors energy harvesting from the consumed food. Our objective was to characterize the composition of gut microbiota in adult Saudis and investigate possible association with lifestyle and dietary practices. Feces from 104 Saudi volunteers (48% males) were tested for microbiota by sequencing the V3-V4 region of bacterial 16S ribosomal RNA (rRNA). For all participants, data were collected related to their lifestyle habits and dietary practices. The relative abundance (RA) of Fusobacteria was significantly higher in normal weight Saudis (P = 0.005, false discovery rate-FDR = 0.014). Individuals who consumed more coffee presented marginally significant more RA of Fusobacteria (P = 0.02, FDR = 0.20) in their gut microbiota compared to those reporting low or no coffee intake, but the RA of Fusobacteria was significantly higher in smokers compared to non-smokers (P = 0.009, FDR = 0.027). The RA of Fusobacteria was also significantly higher in those reporting daily consumption of bread (P = 0.005, FDR = 0.015). At the species level, the gut microbiota of people who consumed coffee was dominated by *Bacteroides thetaiotaomicron* followed by *Phascolarctobacterium faecium* and *Eubacterium rectale*. Similarly, the gut microbiota of smokers was also enriched by *B. thetaiotaomicron* and *Lactobacillus amylovorus*. Smoking cessation, bread and coffee consumption induce changes in the intestinal microbial composition of Saudis. This indicates the significance of diet and lifestyle
practices in the determination of the composition of the gut microbiota, which could possibly
lead later to changes in metabolic profile and weight.

**Introduction**

The gastrointestinal microbiome establishes a stable symbiotic, mutually beneficial relation-
ship with the host. However, it is affected by age, drugs and diet among other factors [1, 2].
Dietary and lifestyle habits positively alter the composition of gut microbiota to harvest energy
from consumed food [3, 4]. Recent studies revealed that coffee has antibiotic effects and its
consumption can modify the gut microbiota's composition [5]. Similarly, smoking modifies
the gut microbiota causing an increase in Firmicutes and Actinobacteria population and leads
to a decrease in microbial diversity [6]. Moreover, the type of food intake by the human host
also influences the gut microbiota composition and diversity [1, 2]. The gut microbiota is
important to the process of harvesting, storing and expending energy attained from the diet. It
has been reported that the microbiota of the duodenal aspirates from obese people was able to
modulate fatty acid and sucrose breakdown pathways, likely due to dietary imbalance and help
the host in harvesting energy and lead to increasing adiposity [7].

Saudi Arabia has rapidly developed economically and socially in the past decades with an impact on the lifestyle of its citizens. Such a rapid lifestyle change has influenced the dietary consumption habits of the whole population and resulted in a trend towards intake of high energy and processed fats with high levels of fat, salt and sugar [8]. As a result, obesity has appeared as an endemic disorder that is quickly emerging as a major problem in Saudi Arabia. According to a local study, around 70% of adult Saudis are either overweight or obese [9].

The body fat percentage among Saudi children and teenagers is increasing and leading to an emergence of obesity [10]. To date, there are few studies that have characterized the gut microbiota composition among Saudis [11]. The objective of this study was to use high-throughput 16S ribosomal RNA (rRNA) gene amplicon sequencing to characterize the gut microbiota of a Saudi population and determine if smoking, coffee and bread consumption had an impact on their microbiome.

**Methods**

**Subject selection criteria**

**Subjects and study design.** This cross-sectional study was conducted between January 2015 and December 2015 on healthy adults of both genders, aged 18–55 years and of different body mass index (BMI), recruited from the student population and others attending King Abdulaziz University Medical campus, as well as members of their families and friends. Exclusion criteria included: history of colon cancer, inflammatory bowel disease, acute or chronic diarrhea in the previous 8 weeks and treatment with antibiotics in the 2 months prior to fecal sampling, and intake of medication or supplements. The study was approved by the Ethics Committee of human Research at King Abdulaziz University under agreement number 014-CEGMR-2-ETH-P. All participants were asked to sign a written informed consent after being informed about the purpose of the study and ensured about confidentiality of the data. They were then requested to fill out a questionnaire covering their socio-demographic information, medical history and lifestyle practices. In addition, a structured food frequency questionnaire (FFQ) was administered to evaluate their dietary practices, and weight and height measurements were taken using standardized techniques. The used questionnaire was
previously described partially or fully and used in other manuscripts [12–15]. Weight and height were used to calculate body mass index (BMI = kg m⁻²) and the WHO criteria [16] were used to classify participants as underweight, normal, overweight and obese. Weight categories were defined according to BMI as follows: normal 20–25 kg m⁻², underweight 18–20 kg m⁻², overweight 25–30 kg m⁻², and obese >30 kg m⁻². Stool samples were collected in aseptic conditions with clean, dry screw-top containers and immediately stored at -20 °C.

**Extraction of DNA from stool samples and 16S rRNA sequencing using MiSeq technology.** All participants' stool samples were extracted using a deglycosylation protocol as follows: 250 µL of each sample was placed in a 2 mL tube containing a mixture of acid-washed glass beads (Sigma, Aldrich) and with two or three 0.5 mm glass beads. Mechanical lysis was performed by bead-beating the mixture using a Fast Prep BIO 101 apparatus (Qbiogene, Strasbourg, France) at maximum speed (6.5) for 3×30 seconds. The supernatant was centrifuged at 12,000 rpm for 10 min and the pellet retained. A mixture containing 2 µL of 10×glycoprotein denaturing buffer EndoHf (New England Biolabs) and 17 µL of H₂O was added and heated at 100 °C for 10 minutes. Deglycosylation was performed adding a mixture of 2 µL of 10×G5 reaction buffer (ref B1702 New England Biolabs), 2 µL of EndoHf (New England Biolabs), 2 µL of cellulase (Sigma) and 16 µL of H₂O. The preparation was then incubated overnight at 37 °C. Finally, DNA was extracted using the NucleoSpin® Tissue Mini Kit (Macherey Nagel, Hoerdt, France) according to a previously described protocol [17]. The quantity, purity, integrity and size of DNA and its amenability to PCR amplification were assessed. The concentration of each DNA extraction was measured by a Qubit assay with the high sensitivity kit (Life technologies, Carlsbad, CA, USA) according to the Nextera XT DNA sample prep kit (Illumina) and diluted to 1 ng aliquots of each metagenome for paired end sequencing analysis. DNA extracts were dispensed into 10- to 20-µL single-use aliquots and frozen at -20 °C to avoid repeat freeze-thaw cycles prior to downstream analyses. Samples were then sequenced targeting the V3–V4 regions of the 16S rRNA gene using MiSeq technology as previously described [18, 19].

**Data processing: Filtering the reads, dereplication and clustering.** Paired end fastq files were assembled using FLASH [20]. A total of 7518258 joined reads were filtered and then analyzed in QIIME by choosing chimera slayer for removing chimera and Uclust [16, 20] for Operational Taxonomic Units (OTU) extraction as described previously [18, 19]. All reads were clustered with a threshold of 97% identity to obtain OTU. Extracted OTUs were blasted against SILVA123 SSU database [21] of release and taxonomy were assigned to a species if they matched one with at least 97% identity, as previously described [22, 23]. Briefly, for each OTU, representative sequences were extracted and were searched against the reference database. For each unique representative sequence, we extracted the best matches from the reference database and sorted them by decreasing percentage of similarity rounded to the nearest integer. We used the reference sequences with >97% similarity (or the highest available) for taxonomic assignments into species. When multiple matches with the same percentage of similarity were present, the taxonomy of each rank was obtained by consensus [16, 24]. OTU not assigned to any species were considered "unidentified". As several OTUs matched identical species, the total number of identified species and the number of unidentified OTU was expected to be smaller than the total number of OTUs.

**Statistical analysis.** A non-parametric Kruskal-Wallis test was used with an adjustment for multiple comparisons using the *post hoc* Benjamini-Hochberg correction from the OMICS package in XLSTAT V.2016.02 (Addinsoft, Paris, France). After normalization, an exploratory Pearson principal component analysis was first performed using the relative abundance at the bacterial phyla and general level as an active variable and age, sex, and coffee consumption, smoking and bread consumption as supplementary elements using XLSTAT v2014.3.07 (Addinsoft). Statistical analyses were performed using GraphPad Prism, V 5.0 (La Jolla,
<table><caption>Table 1. Description of the volunteers.</caption><thead><tr><th></th><th></th><th>Number of volunteers (%)</th><th>Median age ± Interquartile range</th></tr></thead><tbody><tr><td rowspan="4">BMI</td><td>Underweight</td><td>21 (20.2%)</td><td>23.0 ± 10</td></tr><tr><td>Normal weight</td><td>31 (29.8%)</td><td>23.0 ± 6</td></tr><tr><td>Overweight</td><td>28 (26.9%)</td><td>25.0 ± 8</td></tr><tr><td>Obese</td><td>24 (23.1%)</td><td>27.0 ± 7</td></tr><tr><td>Smokers</td><td></td><td>19 (19.2%)</td><td>28.0 ± 8</td></tr><tr><td rowspan="3">Coffee</td><td>Daily consumption</td><td>5 (5.0%)</td><td>20.0 ± 8</td></tr><tr><td>Often consumption</td><td>53 (52.5%)</td><td>25.0 ± 7</td></tr><tr><td>No consumption</td><td>43 (42.5%)</td><td>24.0 ± 9</td></tr><tr><td rowspan="3">Bread</td><td>Daily consumption</td><td>67 (67.0%)</td><td>24.0 ± 8</td></tr><tr><td>Often consumption</td><td>21 (21.0%)</td><td>26.0 ± 8</td></tr><tr><td>No consumption</td><td>12 (12.0%)</td><td>25.0 ± 6</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0230895.t001

California, USA). The Fisher's exact test was used for dichotomous variables, the Mann-Whitney test for continuous variables and the Spearman test for correlations. Principal coordinate analysis (PCoA) was obtained using the weighted unifrac distance after data rarefaction at the depth of 50,000 reads per sample, and the Adonis test was performed in QIIME [20, 25]. Linear discriminant analysis effect size (LEfSe) was performed on relative abundance at the genus and species levels using parameters previously recommended, including per-sample normalization of the sum of the values to 1 M (http://huttenhower.sph.harvard.edu/galaxy/) [23].

## Results

We tested 104 volunteers (48% males) with median age ± interquartile range (IR) was 24 ± 7.7 (Table 1).

### Composition of gut microbiota and body mass index

The analysis of the high-quality trimmed reads revealed that the gut microbiota of the subjects contained sequences mostly belonging to Firmicutes, followed by Bacteroidetes. The relative abundance of Fusobacteria was significantly higher (87%) in normal weight subjects (P = 0.005, FDR = 0.014). Lentisphaerae was marginally significantly higher in the gut microbiota of obese subjects than in the gut microbiota of normal (P = 0.01, FDR = 0.11) and underweight subjects (P = 0.008, FDR = 0.022). Lentisphaerae were not detectable in the gut microbiota of overweight subjects. Euryarchaeota and Synergistetes were not detectable in the gut microbiota of obese individuals. PCoA of the overall composition of the genera communities among the groups did not reveal differences among the microbiomes of underweight, normal weight, overweight and obese individuals (S1 Fig). The microbial richness estimated by the Shannon index did not reveal differences among the microbiomes of underweight, normal weight, overweight and obese individuals.

### Composition of gut microbiota and coffee consumption

We found that coffee consumption influenced the gut microbiome of the subjects. Indeed, the gut microbiota of subjects who consumed coffee daily presented a significantly higher relative abundance of Synergistetes when compared to those who consumed coffee often (P = 0.01, FDR = 0.10) or not at all (P = 0.01, FDR = 0.08). Individuals who often consumed coffee presented marginally significantly more relative abundance of Fusobacteria (86%) than those who
consumed coffee daily or not at all (P = 0.02, FDR = 0.20). Individuals who often consumed coffee daily or often had significantly more Tenericutes and Euryarchaeota than individuals who did not consume coffee. Moreover, PCoA of the overall composition of the genera communities between the groups did not reveal differences among the microbiomes of individuals concerning coffee consumption (S1 Fig). The microbial richness did not reveal differences among the microbiomes of subjects who consumed coffee daily, often or not at all (S2 Fig). *Bifido bacterium* spp relative abundance did not differ among individuals who consumed coffee or not (P = 0.6, FDR >0.2). However, we found that the gut microbiota of subjects who consumed coffee daily was enriched by *Eubacterium rectale* whereas the microbiota of those who did not consume coffee was enriched by *Marinomona sarctica*. The gut microbiota of individuals who consumed coffee often was enriched by *Bacteroides thetaiotaomicron* (S3 Fig). Overall, the gut microbiota of subjects who consumed coffee was dominated by *B. thetaiotaomicron* followed by *Phascolarcto bacterium faecium* and *E. rectale*, whereas the gut microbiota of subjects who did not consume coffee was dominated by *Roseburia faecis* (Fig 1).

<figure><img src="image_5.png" /><figcaption>**Fig 1.** Linear Discriminant Analysis (LDA) scores of differentially abundant species among individuals who consume coffee (green) or not (red). The LDA scores represent the effect size of each abundant species. Species enriched in each group with an LDA score >2 are considered.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0230895.g001
## Composition of gut microbiota and bread consumption

The relative abundance of Fusobacteria (88%) was significantly higher in the gut microbiota of subjects who daily consumed bread (P = 0.005, FDR = 0.015). Fusobacteria were not detectable in the gut microbiota of subjects who did not consume bread. Similarly, the relative abundance of Synergistetes and Lentisphaerae were significantly higher in the gut microbiota of subjects who often consumed bread (P = 0.009, FDR = 0.028 and P = 0.004, FDR = 0.011 respectively). Synergistetes, Cyanobacteria, Euryarchaeota and Lentisphaerae were not detectable in the gut microbiota of subjects who did not consume bread. PCoA of the overall composition of the genera communities among the groups did not reveal differences among subjects concerning bread consumption (S1 Fig). Moreover, bread consumption did not affect the gut microbial richness (S2 Fig).

## Composition of gut microbiota in relation to smoking

The analysis of the high-quality trimmed reads revealed that the gut microbiota of subjects who were smokers contained sequences mostly belonging to Bacteroidetes. Moreover, the gut microbiota of the non-smokers had a significantly higher relative abundance of Fusobacteria (83%) (P = 0.009, FDR = 0.027) and Tenericutes (P = 0.008, FDR = 0.018) than the gut microbiota of smokers. Synergistetes, Lentisphaerae and Euryarchaeota were not detectable in the gut microbiota of the non-smokers. PCoA of the overall composition of the genera communities between the groups did not reveal differences between the microbiomes of smokers and non-smokers (S1 Fig). Similarly, we did not find a difference in the microbial richness, between the gut microbiomes of smokers and non-smokers (S2 Fig). Analysis at the species level revealed that the gut microbiota of smokers was also enriched by *B. thetaiotaomicron* and *Lactobacillus amylovorus*, whereas the gut microbiota of the individuals who did not smoke was enriched by *Dialister invisus* and *Ruminococcus bromii* (Fig 2). Finally, we found that the gut microbiota of subjects who were smokers and who consumed coffee were dominated by *B. thetaiotaomicron* followed by *Bacteroides massiliensis*, whereas the gut microbiota of subjects who did not smoke and consumed coffee by *Klebsiella pneumoniae* (Fig 3).

## Discussion

In this study, we compared the gut microbiomes of Saudis and we provide important information about the impact of different dietary habits on their gut microbiome. To the best of our knowledge, this is the second attempt to evaluate the statistical variability of the gut microbiota among Saudis. Moreover, this is the first study that involved female volunteers as Saudi Arabia has a conserved society and people are reluctant to volunteer for studies involving stool samples, particularly females. A limitation of our study was that we did not measure on our Saudi volunteers the waist circumference, visceral fat and the fat content (%), which are more relevant factors than BMI. Metagenomics studies are the preferential technique for the exploration of the gut microbiota diversity but have considerably low reproducibility due to the differences in sampling, sample conservation, DNA extraction protocol, sequencing method, and data analysis strategy [7, 11, 18, 26]. Although next generation sequencing (NGS) technology has also been used in several studies for the exploration of the gut microbiota to the species level [7, 11, 18], an unequivocal identification on the species level is not always possible due to the high similarities of the 16S rRNA gene sequencing in some bacterial species and due to the small sequence that is obtained that is roughly one third of the full 16S rRNA gene.

Previously, Bacteroidetes have been linked to weight gain [1]. The reasons for the changes of the complex gastrointestinal microbiome ecosystem due to obesity are still controversial. In our previous pilot study, we evaluated changes of the gut microbiota of Saudis and reported
<figure><img src="image_6.png" /><figcaption>**Fig 2.** LDA scores of differentially abundant species according to smokers (green) and individuals who do not smoke (red). The LDA scores represent the effect size of each abundant species. Species enriched in each group with an LDA score >2 are considered.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0230895.g002

that Saudis had a significantly lower diversity in comparison to a group of French people. The
difference may have been due to the food intake with less diversity among Saudis in compari-
son to the French [11]. Moreover, we found that obese Saudis possessed significantly more Fir-
micutes than normal weight Saudis [11]. However, we did not find differences in the gut
microbiota population of Firmicutes or Bacteroidetes among individuals with different weight

<figure><img src="image_7.png" /><figcaption>**Fig 3.** LDA scores of differentially abundant species according to individuals who consumed coffee and were smokers (red) and individuals who do not smoke and consume coffee (green). The LDA scores represent the effect size of each abundant species. Species enriched in each group with an LDA score >2 are considered.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0230895.g003
phenotype in this study. Contradictory results on the microbiome can be noted from different
analyses conducted from the same laboratory [26, 27]. In our study we believe that these dis-
crepancies were due to subject selection. In a recent systematic review based on animal experi-
ments and clinical studies, it was found that a weight loss after a bariatric surgery was
associated with an increase of Fusobacteria in the gut microbiota [28]. This may explain our
findings that Fusobacteria were significantly higher in normal weight Saudis compared to the
overweight or obese ones.

Diet plays an important role on the gut microbiomes and can influence the bacterial diversity from carnivore to omnivore to herbivore [29]. Indeed, there are various reports of the influence of dietary components, as well as lifestyle on gut microbiota in humans from different geographical regions [28-34].

A high-fiber diet has been associated with an enrichment of the gut microbiota [31], and gut microbiota differences were associated with different types of diet [32–34]. Previous reports have indicated that Western populations have lower microbial richness than non-Western populations [32]. In a previous study we found that the gut microbiota of Bedouins, a nomadic population in Saudi Arabia, eating a diet that is mostly based on vegetables, fruits, chicken, dairy products, fermented food and rice presented an increased biodiversity when compared to urban Saudis [18]. Even though bread consumption did not affect the gut microbial richness in our study, the relative abundance of Fusobacteria was significantly higher in the gut microbiota of subjects consuming bread on a daily basis, while they were not detectable in the gut microbiota of those who did not consume bread. Similarly, the relative abundance of Synergistetes and Lentisphaerae were significantly higher in the gut microbiota of subjects who often consumed bread but were not detectable in the gut microbiota of those who did not consume bread. Similar effects on microbial richness were reported in mice fed whole-wheat bread, and the relative abundance of microbial species also differed [35]. However, a recent human study reported absence of a significant species compositional change in the microbiota of persons fed either traditionally made sourdough leavened whole-grain bread or industrially made white bread for one week, even though some clinical parameters and the glycemic index of the subjects were affected [36]. The researchers proposed that the short period of the study was not enough to induce changes in the gut microbiota. However, altering the diet will quickly lead to the modulation of the gut microbiome, thus facilitating the adaptability to diverse human lifestyles diets [35, 36].

We also found that coffee consumption influences the gut microbiota of Saudis. Few studies have evaluated the relationship between coffee consumption and gut microbiota [5]. Gniech-witz *et al*. found that the fecal microbiota did not change after coffee consumption [37]. Similarly, healthy adults who consumed three cups of coffee daily for three weeks did not present important modifications to their gut microbiota [38]. In contrast, some individuals presented an increase in the population of *Bifidobacterium* spp in their gut microbiota [38]. In our studied individuals, we did not find that coffee consumption influenced the population of gut *Bifidobacterium* spp. However, we found that the gut microbiota of people who consumed coffee often was enriched by *E. rectale* and *B. thetaotaomicron*. The same bacterial species also dominated in the gut microbiota of smokers. In a recent study, the fecal microbiota of individuals undergoing smoking was also enriched by *Clostridium coccoides*, *E. rectale*, and *Clostridium leptum* subgroup [6, 39]. In addition, a cross-sectional study using fluorescent *in situ* hybridization targeting selected bacterial groups reported that smokers suffering by Crohn's disease presented higher Bacteroidetes–*Prevotella* comparing to nonsmoking patients [40]. It was previously found that the intestinal bacterial composition and diversity began to change from the fourth week after smoking cessation. Although our study was not designed to investigate potential causative mechanisms of the smoking cessation on the gut microbiota, smoking
seems to affect bowel mucosa and mucin expression. Moreover, it is possible that the absorp-
tion of several toxic chemical compounds in cigarette smoke could change metabolism and
alter the composition of gut microbiota. In addition, immune modification due to smoking
may indirectly affect the composition of gut microbiota.

**Conclusion**

We provide evidence that smoking cessation and coffee consumption induces changes in the intestinal microbial composition of Saudis. Food affects the gut microbiota of individuals and we proved that bread consumption modified that gut microbiota of Saudis. Most of these gut microbiota modifications affect the population of Fusobacteria, a phylum that is not a large member within human gut microbiota. We point out the importance of more intensive research in the future to understand the impact of Fusobacteria as also of diet, coffee and smoking cessation on the gut microbiota.

**Supporting information**

S1 Fig. Principle coordinate analysis (PCoA) of the overall composition of the genera communities among the groups.

(TIF)

S2 Fig. Gut microbiota Chao1 indexes according to a, coffee consumption; b, smoking; and c, bread consumption.

(TIF)

S3 Fig. LDA scores of differentially abundant species among individuals based on daily (red) or often (blue) coffee consumption and no coffee consumption (green). The LDA scores represent the effect size of each abundant species. Species enriched in each group with an LDA score >2 are considered.

(TIF)

**Acknowledgments**

The authors gratefully acknowledge Lubna Al-Shaikh and Fares Al-Shahrani for their help in
sample collection.

**Author Contributions**

**Data curation:** Dipankar Bachar, Ahmed Al-Hejin.

**Formal analysis:** Steve Harakeh, Emmanouil Angelakis, Timokratis Karamitros, Ahmed A. Mirza.

**Investigation:** Suha A. Farraj, Esam I. Azhar.

**Methodology:** Suhad Bahijri, Ghada Ajabnoor, Suha A. Farraj, Turki Al Amri, Abdalla Ahmed.

**Software:** Dipankar Bachar.

**Supervision:** Steve Harakeh, Sulaiman M. Alfadul, Raoult Didier.

**Writing – original draft:** Ghada Ajabnoor.

**Writing – review & editing:** Steve Harakeh, Emmanouil Angelakis, Ahmed A. Mirza, Raoult Didier, Esam I. Azhar.
References

1. Angelakis E, Armougom F, Million M, Raoult D. The relationship between gut microbiota and weight gain in humans. Future microbiology. 2012; 7(1):91–109. Epub 2011/12/24. https://doi.org/10.2217/fmb.11.142 PMID: 22191449.

2. Harakeh SM, Khan I, Kumosani T, Barbour E, Almasaudi SB, Bahijri SM, et al. Gut Microbiota: A Contributing Factor to Obesity. Front Cell Infect Microbiol. 2016; 6:95. https://doi.org/10.3389/fcimb.2016.00095 PMID: 27625997

3. Lynch SV, Pedersen O. The Human Intestinal Microbiome in Health and Disease. N Engl J Med. 2016; 375(24):2369–79. Epub 2016/12/16. https://doi.org/10.1056/NEJMra1600266 PMID: 27974040.

4. Spor A, Koren O, Ley R. Unravelling the effects of the environment and host genotype on the gut microbiome. Nature reviews Microbiology. 2011; 9(4):279–90. Epub 2011/03/17. https://doi.org/10.1038/nrmicro2540 PMID: 21407244.

5. Nakayama T, Oishi K. Influence of coffee (Coffea arabica) and galacto-oligosaccharide consumption on intestinal microbiota and the host responses. FEMS microbiology letters. 2013; 343(2):161–8. Epub 2013/04/05. https://doi.org/10.1111/1574-6968.12142 PMID: 23551139.

6. Biedermann L, Zeitz J, Mwinyi J, Sutter-Minder E, Rehman A, Ott SJ, et al. Smoking Cessation Induces Profound Changes in the Composition of the Intestinal Microbiota in Humans. PloS one. 2013.

7. Angelakis E, Armougom F, Carriere F, Bachar D, Laugier R, Lagier JC, et al. A Metagenomic Investigation of the Duodenal Microbiota Reveals Links with Obesity. PloS one. 2015; 10(9):e0137784. Epub 2015/09/12. https://doi.org/10.1371/journal.pone.0137784 PMID: 26356733

8. Darwish MA A-SG, Albahrani S, Sabra AA. Lifestyle and Dietary Behaviors among Saudi Preschool Children Attending Primary Health Care Centers, Eastern Saudi Arabia. Int J Family Med. 2014; 2014:432732. https://doi.org/10.1155/2014/432732 PMID: 25114804

9. Alqarni SSM. A Review of Prevalence of Obesity in Saudi Arabia. Journal of Obesity & Eating Disorders. 2016; 2(2):25. https://doi.org/10.21767/2471-8203.100025

10. Al-Hazzaa HM, Abahussain NA, Al-Sobayel HI, Qahwaji DM, Musaiger AO. Lifestyle factors associated with overweight and obesity among Saudi adolescents. BMC public health. 2012; 12:354. Epub 2012/05/18. https://doi.org/10.1186/1471-2458-12-354 PMID: 22591544

11. Yasir M, Angelakis E, Bibi F, Azhar EI, Bachar D, Lagier JC, et al. Comparison of the gut microbiota of people in France and Saudi Arabia. Nutrition & diabetes. 2015; 5:e153. Epub 2015/04/29. https://doi.org/10.1038/nutd.2015.3 PMID: 25915742

12. K ET Aljoudi S. B, Alsomali F. A, Haddawi Y, Meer E. M, Binsalman A, Bindajam O, et al. Lifestyle Factors and their Relation to Measures of Obesity Amongst Adults Living in Jeddah- Saudi Arabia: A Cross-Sectional Study. Curr Res Nutr Food Sci. 2015; 3(2). https://doi.org/10.12944/CRNFSJ.3.2.02

13. Bahijri SM, Jambi HA, Al Raddadi RM, Ferns G, Tuomilehto J. The Prevalence of Diabetes and Prediabetes in the Adult Population of Jeddah, Saudi Arabia—A Community-Based Survey. PLoS One. 2016; 11(4):e0152559. Epub 2016/04/02. https://doi.org/10.1371/journal.pone.0152559 PMID: 27035920

14. Al-Raddadi R, Bahijri S, Borai A, AlRaddadi Z. Prevalence of lifestyle practices that might affect bone health in relation to vitamin D status among female Saudi adolescents. Nutrition. 2018; 45:108–13. Epub 2017/11/14. https://doi.org/10.1016/j.nut.2017.07.015 PMID: 29129230.

15. Mydaa K. Aljabri RA-R, Bahijri Suhad M., Al Ahmadi Jawaher, Ajabnoor Ghada, Jambi Hanan A. Factors associated with adherence to Mediterranean diet among Saudi non-diabetic patients attending primary health care centers: A cross-sectional study. Journal of Taibah University Medical Sciences. 2019; 14(2):139–48. https://doi.org/10.1016/j.jtumed.2019.01.006 PMID: 31435404

16. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013; 10(10):996–8. Epub 2013/08/21. https://doi.org/10.1038/nmeth.2604 PMID: 23955772.

17. Zoetendal EG, Booijink CC, Klaassens ES, Heilig HG, Kleerebezem M, Smidt H, et al. Isolation of RNA from bacterial samples of the human gastrointestinal tract. Nature protocols. 2006; 1(2):954–9. Epub 2007/04/05. https://doi.org/10.1038/nprot.2006.143 PMID: 17406329.

18. Angelakis E, Yasir M, Bachar D, Azhar EI, Lagier J-C, Bibi F, et al. Gut microbiome and dietary patterns in different Saudi populations and monkeys. Scientific Reports. 2016; 6:32191. https://doi.org/10.1038/srep32191 PMID: 27578328

19. Angelakis E, Bachar D, Henrissat B, Armougom F, Audoly G, Lagier J-C, et al. Glycans affect DNA extraction and induce substantial differences in gut metagenomic studies. Scientific Reports. 2016; 6:26276. https://doi.org/10.1038/srep26276 PMID: 27188959

20. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nature Methods. 2010; 7(5):335. https://doi.org/10.1038/nmeth.f.303 PMID: 20383131
21. 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 research. 2013; 41 (Database issue):D590–6. Epub 2012/11/30. https://doi.org/10.1093/nar/gks1219 PMID: 23193283

22. Guillou L, Bachar D, Audic S, Bass D, Berney C, Bittner L, et al. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic acids research. 2013; 41 (Database issue):D597–604. Epub 2012/11/30. https://doi.org/10.1093/nar/gks1160 PMID: 23193267

23. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome biology. 2011; 12(6):R60. Epub 2011/06/28. https://doi.org/10.1186/gb-2011-12-6-r60 PMID: 21702898

24. Patin NV, Kunin V, Lidstrom U, Ashby MN. Effects of OTU clustering and PCR artifacts on microbial diversity estimates. Microb Ecol. 2013; 65(3):709–19. Epub 2012/12/13. https://doi.org/10.1007/s00248-012-0145-4 PMID: 23233090.

25. CHAO A. ST. Nonparametric estimation of Shannon's index of diversity when there are unseen species in sample. Environmental and Ecological Statistics. 2003; 10:429–43. https://doi.org/10.1023/A:1026096204727

26. Chatelier EL, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013; 500(7464):541. https://doi.org/10.1038/nature12506 PMID: 23985870

27. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010; 464(7285):59. https://doi.org/10.1038/nature08821 PMID: 20203603

28. Guo Y, Huang ZP, Liu CQ, Qi L, Sheng Y, Zou DJ. Modulation of the gut microbiome: a systematic review of the effect of bariatric surgery. European journal of endocrinology. 2018; 178(1):43–56. Epub 2017/09/17. https://doi.org/10.1530/EJE-17-0403 PMID: 28916564.

29. Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science (New York, NY). 2008; 320(5883):1647–51. Epub 2008/05/24. https://doi.org/10.1126/science.1155725 PMID: 18497261

30. Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010; 464(7285):59–65. Epub 2010/03/06. https://doi.org/10.1038/nature08821 PMID: 20203603

31. De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proceedings of the National Academy of Sciences of the United States of America. 2010; 107(33):14691–6. Epub 2010/08/04. https://doi.org/10.1073/pnas.1005963107 PMID: 20679230

32. Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. Nature. 2012; 486(7402):222–7. Epub 2012/06/16. https://doi.org/10.1038/nature11053 PMID: 22699611

33. Escobar JS, Klotz B, Valdes BE, Agudelo GM. The gut microbiota of Colombians differs from that of Americans, Europeans and Asians. BMC microbiology. 2014; 14:311. Epub 2014/12/17. https://doi.org/10.1186/s12866-014-0311-6 PMID: 25495462

34. Schnorr SL, Candela M, Rampelli S, Centanni M, Consolandi C, Basaglia G, et al. Gut microbiome of the Hadza hunter-gatherers. Nature Communications. 2014; 5:3654. https://doi.org/10.1038/ncomms4654 PMID: 24736369

35. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2013; 505(7484):559. https://doi.org/10.1038/nature12820 PMID: 24336217

36. Cotillard A, Kennedy SP, Kong LC, Prifti E, Pons N, Chatelier EL, et al. Dietary intervention impact on gut microbial gene richness. Nature. 2013; 500(7464):585. https://doi.org/10.1038/nature12480 PMID: 23985875

37. Gniechwitz D, Reichardt N, Blaut M, Steinhart H, Bunzel M. Dietary fiber from coffee beverage: degradation by human fecal microbiota. Journal of agricultural and food chemistry. 2007; 55(17):6989–96. Epub 2007/07/31. https://doi.org/10.1021/jf070646b PMID: 17658822.

38. Jaquet M, Rochat I, Moulin J, Cavin C, Bibiloni R. Impact of coffee consumption on the gut microbiota: a human volunteer study. International journal of food microbiology. 2009; 130(2):117–21. Epub 2009/02/17. https://doi.org/10.1016/j.ijfoodmicro.2009.01.011 PMID: 19217682.

39. Biedermann L, Brulisauer K, Zeitz J, Frei P, Scharl M, Vavricka SR, et al. Smoking cessation alters intestinal microbiota: insights from quantitative investigations on human fecal samples using FISH.
Inflammatory bowel diseases. 2014; 20(9):1496-501. Epub 2014/07/30. https://doi.org/10.1097/MIB.0000000000000129 PMID: 25072500.

40. Benjamin JL, Hedin CR, Koutsoumpas A, Ng SC, McCarthy NE, Prescott NJ, et al. Smokers with active Crohn's disease have a clinically relevant dysbiosis of the gastrointestinal microbiota. Inflammatory bowel diseases. 2012; 18(6):1092-1100. Epub 2011/11/22. https://doi.org/10.1002/ibd.21864 PMID: 22102318.