Skip to main content

Revealing oral microbiota composition and functionality associated with heavy cigarette smoking

Abstract

Background

Heavy tobacco smoking, a hallmark feature of lung cancer, is drastically predominant in Middle Eastern populations. The precise links between nicotine dependence and the functional contribution of the oral microbiota remain unknown in these populations.

Methods

We evaluated the composition and functional capabilities of oral microbiota with relation to cigarette smoking in 105 adults through shotgun metagenomics using buccal swabs.

Results

The oral microbiota composition in our study subjects was dominated by the phyla Firmicutes, Proteobacteria, Actinobacteria, and Bacteroidetes, in addition to the genera Prevotella and Veillonella, similar to previously described westernized cohorts. Furthermore, the smoker's oral microbiota represented a significant abundance of Veillonella dispar, Leptotrichia spp. and Prevotella pleuritidis when compared to non-smokers. Within the smoking groups, differential relative abundance testing unveiled relative abundance of Streptobacillus hongkongensis, Fusobacterium massiliense, Prevotella bivia in high nicotine dependent compared to low nicotine dependent profiles based on Fagerström Test for Nicotine Dependence. Functional profiling showed marked differences between smokers and non-smokers. Smokers exhibited an enrichment of Tricarballylate utilization and Lactate racemization when compared to the non-smokers. According to their nicotine dependence, enrichment of Xanthosine utilization, p-Aminobenzoyl-Glutamate utilization, and multidrug efflux pump in Campylobacter jejuni biosynthesis modules were detected in the high nicotine dependent group.

Conclusions

These compositional and functional differences may provide critical insight on how variations in the oral microbiota could predispose to respiratory illnesses and smoke cessation relapse in cigarette smokers. In particular, the observed enrichment of Fusobacterium and Prevotella in the oral microbiota possibly suggests an intriguing linkage to gut and lung cancers.

Introduction

The oral microbiota is the second most complex microbial ecosystem after the gut flora, consisting of a dynamic spectrum of microorganisms residing in the oral cavity and its interaction with host genetics, diet, immune system, and many other factors [1]. The bacterial microbiome is the predominant component, with species consisting mainly of obligate aerobes such as Neisseria and Rothia, facultative aerobes such as Streptococcus and Actinomyces, and obligate anaerobes including Firmicutes, Bacteroidetes, and Spirochaetes [2]. The community composition, although similar amongst the buccal mucosa, gingiva, and hard palate; yet is different from the soft surfaces, saliva, and gingival plaques [3]. Saprophytic protozoa such as Entamoeba gingivalis and Trichomonas tenax and fungi such as Candida albicans and Saccharomyces cerevisiae are also native residents of oral microbiota [1].

Despite the similarities in the core microbial composition existing within oral cavities, the species may vary depending on the host's diet and nutrition, genetic predisposition, hormonal factors, antibiotic exposure, alcohol consumption, and repeated infections by pathogenic bacteria. This variation, if pathogenic, is termed dysbiosis, which can cause several alterations to the host's oral and systemic health through multiple pathophysiological processes [4, 5]. Dysbiosis has been reported to be involved in the etiology of oral diseases such as dental caries, gingivitis, and periodontitis; and systemic diseases spanning from infections to cancers, such as respiratory tract infections, gastric ulcers, irritable bowel disease, rheumatoid arthritis, infective endocarditis, and cancers [1, 4, 6].

Tobacco smoking is a well-known preventable cause of death and affects nearly every organ system of the body [7]. The oral cavity is one of the first regions exposed to cigarette smoke and is at a prime disadvantage for increased carcinogenesis, impaired mucosal immunity, and alteration of the oral microbiome [8,9,10]. In turn, smoking increases colonization of the oral cavity by pathogenic bacteria and reduces colonization by commensal bacteria [11, 12]. Smoking enhances biofilm formation and results in greater epithelial adherence by certain pathogens, including Streptococcus pneumonia, Staphylococcus aureus, Streptococcus mutans; thereby, increasing susceptibility to respiratory infections and dental caries respectively in those smokers [8, 10, 12]. Furthermore, smoking contributes to the alteration in the oxygen tension of the oral and upper gastrointestinal microenvironment that encourages persistence of microaerophilic bacteria replacing the commensal beneficial species [12, 13]. Previous studies have shown an increased prevalence of the genera Atopobium, Campylobacter, and Prevotella among smokers and selective depletion of certain phyla, including Proteobacteria [12, 14,15,16]. Thus, tobacco smoking creates a unique dysbiotic environment in the oral cavity, influencing the microbiota composition with far-reaching consequences in the local and systemic health of the host [8]. In this study, we intend to decipher our understanding of the oral microbiota's composition and its alteration due to tobacco smoking and smoking severity (nicotine dependence level). Further, we evaluated the metabolic capabilities of the oral microbiota using shotgun metagenomic sequencing to determine microbial biodiversity and functional capabilities that associate with tobacco smoking in the oral cavity.

Materials and methods

Study population

In this case–control study, we recruited participants over an eight-month period between June 2019 and February 2020 in the emirates of Dubai, Sharjah, and Ajman in the United Arab Emirates. Participants completed self-administered questionnaires that included comprehensive demographic, social, and medical history, among other lifestyle information. Tobacco smokers were defined as those individuals that reported as exclusively cigarette smokers for 5 years or more. Non-smoker controls were defined as individuals who did not report smoking cigarettes or any other tobacco products and were otherwise healthy. We excluded those who reported antibiotic or prescribed probiotic use in the past three months, and those with a preexisting respiratory illness such as asthma and chronic obstructive pulmonary disease.

We have also assessed nicotine dependence by collecting participants' self-administered Fagerström Test for Nicotine Dependence (FTND) scale, as previously described [17]. Briefly, yes or no items are scored with 0 or 1, and multiple-choice items are scored from 0 to 3. The items are summed to yield a total score of 0–10. Higher FTND scores indicate greater physical dependence on nicotine. For further validation, participants also completed the Short Nicotine Dependence Syndrome Scale [18, 19].

We conducted sample size calculation for unmatched case–control with two-sided confidence level of 95 and power of 80 with ratio of controls to cases of 1. Hypothetical proportion of cases with smoking exposure of 50 resulted in minimum of 12 sample size (cases) and a minimum of 12 sample size (controls), total sample size of 24 [20, 21]. During the data collection phase, we collected 539 buccal swabs from 428 non-smokers and 111 smokers, using Isohelix DNA/RNA Buccal Swabs (Isohelix Ltd. Harrietsham, United Kingdom) following the manufacturer's instruction (Isohelix Ltd.). Case–control matching of tobacco smokers and non-smokers group yielded 105 participants consisting of 50 non-smokers and 55 smokers. The swabs were then collected in a sterile container, stored immediately into liquid nitrogen, and then transferred to a − 80 °C freezer until further analysis. Swabs from these 105 participants were further processed for analysis. All participants in the study read and signed an informed consent, and the Research Ethics Committee at the University of Sharjah approved the study protocol.

DNA extraction and library preparation

DNA was extracted using the Qiagen MagAttract PowerSoil DNA KF kit (Formerly MO Bio PowerSoil DNA Kit) using a KingFisher robot. DNA quality was evaluated visually via gel electrophoresis and quantified using a Qubit 3.0 fluorometer (Thermo-Fischer, Waltham, MA, USA). Libraries were prepared with the Illumina Nextera library preparation kit using an in-house protocol (Illumina, San Diego, CA, USA).

Sequencing, data curation, and sequence processing

Paired-end sequencing (150 bp × 2) was done on a NextSeq 500 in medium-output mode. Next, shotgun metagenomic sequence reads were processed with the Sunbeam pipeline [22]. Initial quality evaluation was done using FastQC v0.11.5 (Bioinformatics Group at the Babraham Institute. Software available at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/. Processing took part in four steps: adapter removal, read trimming, low-complexity reads removal, and host-sequence read removal. First, adapter removal was done using Cutadapt v2.6 [23]. Next, trimming was done with Trimmomatic v0.36 [24] using custom parameters (LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36). Then, low-complexity sequences were detected with Komplexity v0.3.6 [22]. Spurious Operational Taxonomic Units (OTUs) reads were determined and removed because they matched one of the pre-specified host/contaminant genomes or due to low complexity or quality. At the end of quality control, the median and range number of quality-filtered reads per sample was 5,062,550 and 606,459, respectively. The remaining reads were taxonomically classified using Kraken2 with the MiniKraken2_v1 database [25] and with the Genome Taxonomy Database (v. 89).

For functional profiling, high-quality (filtered) reads were aligned against the SEED database via translated homology search and annotated to Subsystems, or functional levels, 1–3 using Super-Focus [26].

Statistical analysis

We assessed the alpha diversity with Shannon and Chao1 indices after filtering out spurious OTUs, and then the significance of diversity changes was tested with the Mann Whitney test. Next, we evaluated the beta diversity, underscoring differences across samples; a non-metric multidimensional scaling analysis was used to visualize microbiome similarities. Permutational analysis of variance (PERMANOVA) was used to test for the significance of overall microbiome differences. To assess possible compositional differences in the bacterial community, binomial models (DESEq2 R package) of the form group dependence of taxonomic and subsystem level 3 features were used. P values were calculated with Likelihood Ratio Tests. False Discovery Rate (FDR) p-value adjustments using the Benjamini and Hochberg method were made to correct for multiple testing. An adjusted p-value of less than 0.05 was considered significant. All analyses were conducted in the R environment.

Results

Bacterial summary taxonomic composition

We analyzed buccal swab samples from 105 participants for taxonomic composition, differential abundance, and functional profiling of their oral microbiota. The subjects' characteristics in this study, such as age, gender, body mass index (BMI), ethnicity, and medical history, have been provided (Table 1).

Table 1 Demographics of Study Cohort

First, we evaluated the taxonomic composition generated from high-quality reads and classified them using the MiniKraken2_v1 database [25] as the reference database for bacteria. We aggregated taxa abundances into genera and plotted the relative abundances of the most abundant ones (Additional file 1: Fig. S1). Furthermore, we plotted the relative abundances of the most abundant taxa within the smokers' group based on their FTND score (nicotine dependence); 1−2 (low dependence), 3−4 (low to moderate dependence), 5 – 7 (moderate dependence), and ≥ 8 (high dependence) (Additional file  1: Fig. S2). Nicotine dependence was further evaluated using the Short Nicotine Dependence Syndrome Scale (NDSS-S) [18, 19]. Pearson correlation suggested a significant positive correlation between FTND and NDDS-S for smokers (r = 0.646) (p-value < 0.01) (Data not shown). Next, we estimated alpha diversity (richness and evenness) from taxonomic profiles using Shannon's diversity index and Chao1 richness estimator. No significant differences across different groups were found (Additional file 1: Fig. S3). Last, to assess the overall microbial community compositional changes, PERMANOVA was used to model the effects of smoking and nicotine dependence on oral microbiota composition. We observed a significant taxonomy difference between smoker and non-smoker groups (p-value < 0.04) and a non-significant difference based on nicotine dependence among the smoker group (p-value < 0.09).

Bacterial differential abundance based on smoking and nicotine dependence levels

In order to further assess possible compositional differences in the bacterial community, as suggested in Additional file 1: Figure S1, we conducted negative binomial models as mentioned in methods. First, comparison of the average relative abundance between smokers and non-smokers groups revealed that profiles obtained from smokers have a statistically significant abundance of Veillonella dispar (Log2FoldChange 2.327, P. adjusted value < 0.0000003), Leptotrichia sp000469385 (Log2FoldChange 1.913, P. adjusted value < 0.0013), and Prevotella pleuritidis (Log2FoldChange 1.896, P. adjusted value < 0.00019). On the other hand, there was a statistically significant under-representation of Haemophilus_A (Log2FoldChange − 2.33, P. adjusted value < 0.00007), Gemella cuniculi (Log2FoldChange − 1.976, P. adjusted value < 0.00019), Neisseria subflava_B (Log2FoldChange − 1.87, P. adjusted value < 0.00006), Gemella haemolysans_B (Log2FoldChange − 1.75, P. adjusted value < 0.00085), Neisseria perflava (Log2FoldChange − 1.73, P. adjusted value < 0.0012), Streptococcus oralis_BA (Log2FoldChange − 1.56, P. adjusted value < 0.0004), and Streptococcus mitis_AT (Log2FoldChange − 1.39, P. adjusted value < 0.0013) in smokers (Fig. 1). We further evaluated the average relative abundance among smokers based on nicotine dependence (Fagerström score), which revealed that profiles obtained from more nicotine dependent smokers have a statistically significant abundance of Streptobacillus hongkongensis (Log2FoldChange 4.78, P. adjusted value < 0.00004), Fusobacterium massiliense (Log2FoldChange 4.63, P. adjusted value < 0.00000004), Prevotella sp000163055 (Log2FoldChange 4.42, P. adjusted value < 0.00008), and Prevotella bivia (Log2FoldChange 2.46, P. adjusted value < 0.00024) (Fig. 2).

Fig. 1
figure 1

Differentially abundant taxa between smokers and non-smokers group. Panel shows relative abundance of normalized counts for the top 10 taxa. Results were calculated by negative binomial models (DESEq2 R package) of the form group for differential abundance testing of taxonomic and subsystem level 3 features. P values were calculated with Likelihood Ratio Tests method. All of the above comparison are significant. Smoker and non-smoker corresponding abundance are colored in blue and red, respectively

Fig. 2
figure 2

Differentially abundant taxa based on FNTD nicotine dependence score. Panel shows relative abundance of normalized counts for the top 10 taxa. Results were calculated by negative binomial models (DESEq2 R package) of the form group for differential abundance testing of taxonomic and subsystem level 3 features. P values were calculated with Likelihood Ratio Tests method. All of the above comparison are significant. Nicotine dependence FTND scores; low, low to moderate, moderate, and high are colored in red, green, blue and pink, respectively

Functional profiling of oral microbiota in smoker vs. non-smokers

We used shotgun metagenomic sequencing to determine the functional contribution of the oral microbiota in smokers vs. non-smokers using the SEED hierarchical categorization. Functional profiling showed significant enrichment of Tricarballylate utilization (Log2FoldChange 2.52, P. adjusted value < 0.0013), Aminoglycoside adenylyltransferases (Log2FoldChange 2.39, P. adjusted value < 0.002), Bacteriocins in Lactobacilli (Log2FoldChange 2.29, P. adjusted value < 0.0012), Lactate racemization (Log2FoldChange 1.003, P. adjusted value < 0.0001), and Methionine salvage (Log2FoldChange 0.7, P. adjusted value < 0.0004) in smokers. It also revealed a significant depletion of Two-component Response Regulator of Virulence ResDE (Log2FoldChange − 1.28, P. adjusted value < 0.0009), Listeria Pathogenicity Island LIPI-1 extended (Log2FoldChange − 0.888, P. adjusted value < 0.00006), and CarD (Log2FoldChange − 0.139, P. adjusted value < 0.0007) in smokers (Fig. 3).

Fig. 3
figure 3

Differentially abundant gene functions of smokers vs. non-smokers group. Panel shows relative abundance of normalized counts for functional genes using SEED hierarchical categorization. All of the above comparison are significant. Smoker and non-smoker corresponding abundance are colored in blue and red, respectively

Functional profiling of oral microbiota based on nicotine dependence severity

Finally, we examined differentially abundant gene functions based on the Fagerström score for nicotine dependence among smokers. Pairwise functional differences determined a significant difference between low and more nicotine dependent groups (p-value < 0.02, p-value FDR < 0.05). For example, we show enrichment of Xanthosine utilization (xap region) (Log2FoldChange 3.38, P. adjusted value < 0.00007), p-Aminobenzoyl-Glutamate utilization (Log2FoldChange 1.33, P. adjusted value < 0.00056), Multidrug efflux pump in Campylobacter jejuni (CmeABC operon) (Log2FoldChange 1.14, P. adjusted value < 0.00007), Glycine biosynthesis (Log2FoldChange 1.02, P. adjusted value < 0.00062), Isoleucine degradation (Log2FoldChange 0.989, P. adjusted value < 0.00021). We also noted depletion of Type VI secretion systems (Log2FoldChange − 1.99, P. adjusted value < 0.00027), Rrf2 family transcriptional regulators (Log2FoldChange − 0.598, P. adjusted value < 0.00067), and ABC transporter oligopeptide (TC 3.A.1.5.1) (Log2FoldChange -0.351, P. adjusted value < 0.00001) in the more nicotine dependence group (Fig. 4).

Fig. 4
figure 4

Differentially abundant gene functions based on FNTD nicotine dependence score. Panel shows relative abundance of normalized counts for functional genes using SEED hierarchical categorization. All of the above comparison are significant. Smoking dependence, low, low to moderate, moderate, and high are colored in red, green, blue and pink, respectively

Discussion

The mouth is a highly heterogeneous ecological system with dynamic interplay between the host and oral microbiome [27]. The collective function of microbial communities is a major determinant of homeostasis or dysbiosis, and host factors such as inflammation and dietary sugars may ultimately favor health or disease such as dental caries and periodontitis [28]. In this report, we attempted to explore oral microbial profiles and functions that influence host homeostasis in the background of cigarette smoking. We explored the oral microbiota of chronic tobacco smokers in the Middle-Eastern population and described, for the first time, the functional contribution of the oral bacterial community based on nicotine dependence assessed by the Fagerström scale [17]. A final study population of 105 subjects, with an average age of 30 years, recruited in northern emirates of UAE was used for shotgun metagenomics analysis. We used buccal swabs, a more specific sampling method for the bio-adherent bacteria as compared to mouth wash sampling previously conducted in a UAE based-study [29]. Consistent with several previous reports, we detected a significant taxonomic difference between smoker and non-smoker groups, but no significant differences in terms of microbial diversity and richness, as shown in Additional file 1: Figure S3 [30,31,32]. Interestingly, a previous study conducted in the UAE determined only a marginal significance of the overall oral microbial differences in smokers compared with non-smokers, underscoring the geographic and ethnic contribution [29]. However, our findings were not consistent with other groups reporting a significant change in richness and diversity [33, 34]. The observed fluctuations in oral microbiota richness and diversity reporting by several groups are not unusual and further assert the high complexity and significant effects of several factors such as diet, geography, ethnicity, and host factors. That said, the oral microbiota in our study exhibit comparable dominance of phyla Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes, and genera Prevotella and Veillonella to that of oral microbiota in other populations across the globe [16, 34, 35].

Differential abundance testing of bacterial communities based on nicotine dependence scores revealed a relative abundance of Streptobacillus hongkongensis among more nicotine dependent smokers (high Fagerström score). Previous studies reported the isolation of S. hongkongensis from patients with quinsy, pneumonia, and septic arthritis [36, 37], which was later reported as part of the human oropharynx natural reservoir [38]. Increased risk of developing serious respiratory illnesses might be partly attributed to more nicotine dependent smokers. That said, we acknowledge that the overall number of reads attributed to this species is generally very low and requires further validation. Furthermore, complications of streptobacilliary infections may include endocarditis, brain abscesses, amnionitis, as well as persistent severe arthritis [39].

Smoking tobacco is the single largest risk factor for the development of lung cancers. Several studies established that Fusobacterium nucleatum plays a major role in colorectal carcinogenesis via Fap2 mediated binding to tumor-overexpressed Gal-GalNAc-binding lectin [40,41,42]. Therefore, F. nucleatum was deemed useful as a microbial biomarker for colorectal cancer detection [43]. Interestingly, we discovered that the phylogenetically similar Fusobacterium massiliense, which exhibited substantial sequence similarity with F. nucleatum, has a significant relative abundance among more nicotine dependent smokers. Furthermore, protein–protein BLAST analysis of the Fap2 surface protein of F. nucleatum ATCC 23,726 produced a significant sequence alignment with pyridoxal phosphate-dependent aminotransferase of F. massiliense [41, 44], the active form of vitamin B6. A previous study examined over 44,000 individuals and evaluated their smoking history and B6 vitamin supplement use over 10 years, this study found that high dosages of vitamin B6 supplements were associated with 3–4 folds increase in lung cancer risk among smokers at baseline, although the exact mechanism of this association is not yet known [45]. Fusobacterium, similar to Bacteroides, Bifidobacterium, Actinobacteria, and Proteobacteria, possess a vitamin B6 biosynthesis pathway. Bacteroidetes and Proteobacteria likely produce vitamin B6, starting from deoxyxylulose 5-phosphate and 4-phosphohydroxy-l-threonine [29]. Several prevailing hypotheses may explain the link. First, several B vitamins, including B6, B9 (folate), and B12 interact with homocysteine and methionine in this complex one-carbon metabolism pathway, and disruption of this process may promote carcinogenesis [46]. Second, a study reported that among B6 metabolism markers, it was the inflammation-related changes in a vitamin B6 catabolism marker, the 4-pyridoxic acid/pyridoxal plus pyridoxal 5′-phosphate ratio, which was linked to increased lung cancer risk [47]. Third, excessive supplementation of folic acid and vitamin B12 was found to be associated with changes in DNA methylation of several genes that could be reactivated or deregulated during carcinogenesis [48]. Altogether, perhaps enrichment of F. massiliense among more nicotine dependent smokers suggest a possible linkage to lung cancer in a pyridoxal phosphate-dependent manner. Tobacco smoking, colorectal cancer, and a high relative abundance of gut Prevotella were linked to each other in an intriguing association [49]. Here, we also noted an increase in the relative abundance of Prevotella sp000163055 and Prevotella bivia in oral microbiota of heavy smokers, thereby suggesting a possible downstream effect on the development of colorectal cancers.

The metabolic capabilities of oral microbiota were evaluated using a shotgun metagenomic sequencing approach to determine microbial biodiversity and functional capabilities associated with tobacco smoking in the oral cavity. Functional profiling showed significant enrichment of Tricarballylate utilization among smokers vs. non-smokers group, a good chelator of magnesium that could lead to magnesium deficiency [50]. Magnesium plays a vital role in tobacco addiction by inhibiting several essential steps of nicotine addiction, such as dopamine secretion, NMDA receptor stimulation by glutamate, and the synthesis of substance P and nitric oxide [51, 52]. A previous study showed a significant decrease in the number of cigarettes smoked and Fagerström scores after 28 days of magnesium therapy [53]. This observation of enriched bacterial genes involved in Tricarballylate utilization among smokers suggests an intriguing role of oral dysbiosis in maintaining nicotine addiction. Moreover, a significant increase in the nickel-dependent lactate racemase enzymes was observed in smokers, consistent with the toxic nickel exposure from tobacco smoking [54, 55].

Finally, we examined the differentially abundant gene functions in correlation with the Fagerström score for nicotine dependence among smokers. Significant enrichment of xanthosine utilization was observed among more nicotine dependent smokers, which is a catabolite of purine nucleotides that leads to caffeine synthesis [56]. This enrichment could be linked to the positive association between smoking and coffee consumption, in which heavy smokers require greater coffee consumption than others to obtain an equivalent satisfactory effect of caffeine, as reported in a study of two European cohorts [57]. Lastly, we noted an enrichment of the Multidrug efflux pump in Campylobacter jejuni (CmeABC operon) biosynthesis module in the heavy smokers' group, an important component of bacterial virulence that can predispose heavy smokers to additional risk of tobacco-related morbidity and mortality [58]. It is important to mention that our findings need further validation on a larger cohort. The data obtained from self-administered questionnaires was subject to self-reporting bias; however, a study staff was available during the questionnaire to answer any questions.

Conclusions

We used the shotgun metagenomics approach to shed new light on the complex functional profiles of the oral microbiota in tobacco smokers from the Middle East. To the best of our knowledge, this is the first report on oral microbiota role in heavy smoking among Middle Eastern populations based on nicotine dependence assessed by the Fagerström test. Our data identified significant compositional and functional variations in oral microbial communities, especially among the more nicotine dependent (heavy smokers) that have been linked to several respiratory illnesses and smoking cessation relapse. We hope this information may help us to understand the oral microbiome compositional changes in smokers and their impact on respiratory health and tobacco control strategies.

Availability of data and materials

The data are all published. The data that support the findings of this study are available on request from the corresponding author.

Abbreviations

BMI:

Body mass index

FTND:

Fagerström Test for Nicotine Dependence

HTN:

Hypertension

MENA:

Middle East and North Africa

NDDS-S:

Short Nicotine Dependence Syndrome Scale

SEED:

Categorization system that organizes gene functional categories into a hierarchy

SUPER-FOCUS:

A Tool for Agile Functional Analysis of Shotgun Metagenomic Data

References

  1. Wade WG. The oral microbiome in health and disease. Pharmacol Res. 2013;69:137–43.

    Article  CAS  PubMed  Google Scholar 

  2. Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner ACR, Yu W-H, Lakshmanan A, Wade WG. The Human Oral Microbiome. J Bacteriol. 2010;192:5002–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Segata N, Haake SK, Mannon P, Lemon KP, Waldron L, Gevers D, Huttenhower C, Izard J. Composition of the adult digestive tract bacterial microbiome based on seven mouth surfaces, tonsils, throat and stool samples. Genome Biol. 2012;13:29.

    Article  CAS  Google Scholar 

  4. Appanna V. Dysbiosis, Probiotics, and Prebiotics: In Diseases and Health. 2018, p. 81–122

  5. Zaura E, Nicu EA, Krom BP, Keijser BJF. Acquiring and maintaining a normal oral microbiome: current perspective. Front Cell Infect Microbiol. 2014;4:85.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Kilian M, Chapple ILC, Hannig M, Marsh PD, Meuric V, Pedersen AML, Tonetti MS, Wade WG, Zaura E. The oral microbiome - an update for oral healthcare professionals. Br Dent J. 2016;221:657–66.

    Article  CAS  PubMed  Google Scholar 

  7. World Health O: WHO report on the global tobacco epidemic 2019: Offer help to quit tobacco use. 2019.

  8. Huang C, Shi G. Smoking and microbiome in oral, airway, gut and some systemic diseases. J Transl Med. 2019;17:85.

    Article  Google Scholar 

  9. West R. Tobacco smoking: health impact, prevalence, correlates and interventions. Psychol Health. 2017;32:1018–36.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Tomar SL, Hecht SS, Jaspers I, Gregory RL, Stepanov I. Oral health effects of combusted and smokeless tobacco products. Adv Dent Res. 2018;30:4–10.

    Article  Google Scholar 

  11. Mason MR, Preshaw PM, Nagaraja HN, Dabdoub SM, Rahman A, Kumar PS. The subgingival microbiome of clinically healthy current and never smokers. ISME J. 2015;9:268–72.

    Article  PubMed  Google Scholar 

  12. Brook I. The impact of smoking on oral and nasopharyngeal bacterial flora. J Dent Res. 2011;90:704–10.

    Article  CAS  PubMed  Google Scholar 

  13. Jaspers I. Cigarette smoke effects on innate immune mechanisms in the nasal mucosa. Potential effects on the microbiome. Ann Am Thorac Soc. 2014;11:38–42.

    Article  CAS  Google Scholar 

  14. Yang Y, Zheng W, Cai Q-Y, Shrubsole MJ, Pei Z, Brucker R, Steinwandel MD, Bordenstein SR, Li Z, Blot WJ, et al. Cigarette smoking and oral microbiota in low-income and African–American populations. J Epidemiol Community Health. 2019;73:1108–15.

    Article  PubMed  Google Scholar 

  15. Vallès Y, Inman CK, Peters BA, Ali R, Wareth LA, Abdulle A, Alsafar H, Anouti FA, Dhaheri AA, Galani D, et al. Types of tobacco consumption and the oral microbiome in the United Arab Emirates Healthy Future (UAEHFS) Pilot Study. Sci Rep. 2018;8:1–11.

    Article  CAS  Google Scholar 

  16. Wu J, Peters BA, Dominianni C, Zhang Y, Pei Z, Yang L, Ma Y, Purdue MP, Jacobs EJ, Gapstur SM, et al. Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 2016;10:2435–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Heatherton TF, Kozlowski LT, Frecker RC, Fagerstrom KO. The Fagerstrom test for nicotine dependence: a revision of the Fagerstrom Tolerance Questionnaire. Br J Addict. 1991;86:1119–27.

    Article  CAS  PubMed  Google Scholar 

  18. Shiffman S, Sayette MA. Validation of the nicotine dependence syndrome scale (NDSS): a criterion-group design contrasting chippers and regular smokers. Drug Alcohol Depend. 2005;79:45–52.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Becona E. Becona E, Fernandez E, Lopez A, MiguezMdel C, Castro J, Nogueiras L, Florez G, Alvarez S, Vazquez D: The Short Nicotine Dependence Syndrome Scale (NDSS-S) in Spanish smokers. Psicothema. 2011;2011(23):126–32.

    Google Scholar 

  20. Fleiss JL, Levin B, Paik MC. Statistical methods for rates and proportions. New York: Wiley; 2013.

    Google Scholar 

  21. Kelsey JL, Whittemore AS, Evans AS, Thompson WD. Methods in observational epidemiology. Monograph Epidemiol Biostat; 1996.

  22. Clarke EL, Taylor LJ, Zhao C, Connell A, Lee JJ, Fett B, Bushman FD, Bittinger K. Sunbeam: an extensible pipeline for analyzing metagenomic sequencing experiments. Microbiome. 2019;7:46.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. 2011;2011(17):3.

    Google Scholar 

  24. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019;20:257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Silva GG, Green KT, Dutilh BE, Edwards RA. SUPER-FOCUS: a tool for agile functional analysis of shotgun metagenomic data. Bioinformatics. 2016;32:354–61.

    Article  CAS  PubMed  Google Scholar 

  27. Xu X, He J, Xue J, Wang Y, Li K, Zhang K, Guo Q, Liu X, Zhou Y, Cheng L. Oral cavity contains distinct niches with dynamic microbial communities. Environ Microbiol. 2015;17:699–710.

    Article  PubMed  Google Scholar 

  28. Lamont RJ, Koo H, Hajishengallis G. The oral microbiota: dynamic communities and host interactions. Nat Rev Microbiol. 2018;16:745–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Vallès Y, Inman CK, Peters BA, Ali R, Wareth LA, Abdulle A, Alsafar H, Al Anouti F, Al Dhaheri A, Galani D. Types of tobacco consumption and the oral microbiome in the United Arab Emirates Healthy Future (UAEHFS) Pilot Study. Sci Rep. 2018;8:1–11.

    Article  CAS  Google Scholar 

  30. Jiang Y, Zhou X, Cheng L, Li M. The impact of smoking on subgingival microflora: from periodontal health to disease. Front Microbiol. 2020;11:66.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Bostrom L, Bergstrom J, Dahlen G, Linder LE. Smoking and subgingival microflora in periodontal disease. J Clin Periodontol. 2001;28:212–9.

    Article  CAS  PubMed  Google Scholar 

  32. Gomes SC, Nonnenmacher C, Susin C, Oppermann RV, Mutters R, Marcantonio RA. The effect of a supragingival plaque-control regimen on the subgingival microbiota in smokers and never-smokers: evaluation by real-time polymerase chain reaction. J Periodontol. 2008;79:2297–304.

    Article  PubMed  Google Scholar 

  33. Joshi V, Matthews C, Aspiras M, de Jager M, Ward M, Kumar P. Smoking decreases structural and functional resilience in the subgingival ecosystem. J Clin Periodontol. 2014;41:1037–47.

    Article  PubMed  Google Scholar 

  34. Takeshita T, Kageyama S, Furuta M, Tsuboi H, Takeuchi K, Shibata Y, Shimazaki Y, Akifusa S, Ninomiya T, Kiyohara Y, Yamashita Y. Bacterial diversity in saliva and oral health-related conditions: the Hisayama Study. Sci Rep. 2016;6:22164.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Contreras M, Costello EK, Hidalgo G, Magris M, Knight R, Dominguez-Bello MG. The bacterial microbiota in the oral mucosa of rural Amerindians. Microbiology. 2010;156:3282–7.

    Article  CAS  PubMed  Google Scholar 

  36. Woo PC, Wu AK, Tsang CC, Leung KW, Ngan AH, Curreem SO, Lam KW, Chen JH, Chan JF, Lau SK. Streptobacillus hongkongensis sp. Nov., isolated from patients with quinsy and septic arthritis, and emended descriptions of the genus Streptobacillus and Streptobacillus moniliformis. Int J Syst Evol Microbiol. 2014;64:3034–9.

    Article  CAS  PubMed  Google Scholar 

  37. Eisenberg T, Ewers C, Rau J, Akimkin V, Nicklas W. Approved and novel strategies in diagnostics of rat bite fever and other Streptobacillus infections in humans and animals. Virulence. 2016;7:630–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Lau SK, Chan JF, Tsang CC, Chan SM, Ho ML, Que TL, Lau YL, Woo PC. Human oropharynx as natural reservoir of Streptobacillus hongkongensis. Sci Rep. 2016;6:24419.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Kerr J. Manual of clinical microbiology. BMJ Publishing Group; 2004.

  40. Mima K, Cao Y, Chan AT, Qian ZR, Nowak JA, Masugi Y, Shi Y, Song M, da Silva A, Gu M, et al. Fusobacterium nucleatum in colorectal carcinoma tissue according to tumor location. Clin Transl Gastroenterol. 2016;7:e200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Abed J, Emgard JE, Zamir G, Faroja M, Almogy G, Grenov A, Sol A, Naor R, Pikarsky E, Atlan KA, et al. Fap2 mediates fusobacterium nucleatum colorectal adenocarcinoma enrichment by binding to tumor-expressed Gal-GalNAc. Cell Host Microbe. 2016;20:215–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Brennan CA, Garrett WS. Fusobacterium nucleatum — symbiont, opportunist and oncobacterium. Nat Rev Microbiol. 2019;17:156–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Peng BJ, Cao CY, Li W, Zhou YJ, Zhang Y, Nie YQ, Cao YW, Li YY. Diagnostic performance of intestinal Fusobacterium nucleatum in colorectal cancer: a meta-analysis. Chin Med J (Engl). 2018;131:1349–56.

    Article  Google Scholar 

  44. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    Article  CAS  PubMed  Google Scholar 

  45. Brasky TM, White E, Chen C-L. Long-term, supplemental, one-carbon metabolism-related vitamin B use in relation to lung cancer risk in the vitamins and lifestyle (VITAL) Cohort. J Clin Oncol. 2017;35:3440–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Kim Y-I. Folate and DNA methylation: a mechanistic link between folate deficiency and colorectal cancer? Cancer Epidemiol Prev Biomark. 2004;13:511–9.

    CAS  Google Scholar 

  47. Zuo H, Ueland PM, Eussen SJ, Tell GS, Vollset SE, Nygård O, Midttun Ø, Meyer K, Ulvik A. Markers of vitamin B6 status and metabolism as predictors of incident cancer: the H ordaland H ealth S tudy. Int J Cancer. 2015;136:2932–9.

    Article  CAS  PubMed  Google Scholar 

  48. Kok DE, Dhonukshe-Rutten RA, Lute C, Heil SG, Uitterlinden AG, van der Velde N, van Meurs JB, van Schoor NM, Hooiveld GJ, de Groot LC. The effects of long-term daily folic acid and vitamin B 12 supplementation on genome-wide DNA methylation in elderly subjects. Clin Epigen. 2015;7:121.

    Article  CAS  Google Scholar 

  49. Stewart CJ, Auchtung TA, Ajami NJ, Velasquez K, Smith DP, De La Garza R, 2nd, Salas R, Petrosino JF, . Effects of tobacco smoke and electronic cigarette vapor exposure on the oral and gut microbiota in humans: a pilot study. PeerJ. 2018;6:e4693.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Lewis JA, Escalante-Semerena JC. The FAD-dependent tricarballylate dehydrogenase (TcuA) enzyme of Salmonella enterica converts tricarballylate into cis-aconitate. J Bacteriol. 2006;188:5479–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Hiranita T, Soto PL, Kohut SJ, Kopajtic T, Cao J, Newman AH, Tanda G, Katz JL. Decreases in cocaine self-administration with dual inhibition of the dopamine transporter and sigma receptors. J Pharmacol Exp Ther. 2011;339:662–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Veeneman MM, Broekhoven MH, Damsteegt R, Vanderschuren LJ. Distinct contributions of dopamine in the dorsolateral striatum and nucleus accumbens shell to the reinforcing properties of cocaine. Neuropsychopharmacology. 2012;37:487–98.

    Article  CAS  PubMed  Google Scholar 

  53. Nechifor M, Chelarescu D, Mandreci I, Cartas N. Magnesium influence on nicotine pharmacodependence and smoking. Magnes Res. 2004;17:176–81.

    CAS  PubMed  Google Scholar 

  54. Desguin B, Goffin P, Viaene E, Kleerebezem M, Martin-Diaconescu V, Maroney MJ, Declercq J-P, Soumillion P, Hols P. Lactate racemase is a nickel-dependent enzyme activated by a widespread maturation system. Nat Commun. 2014;5:3615.

    Article  PubMed  CAS  Google Scholar 

  55. Torjussen W, Zachariasen H, Andersen I. Cigarette smoking and nickel exposure. J Environ Monit. 2003;5:198–201.

    Article  CAS  PubMed  Google Scholar 

  56. Ashihara H, Kato M, Crozier A. Distribution, Biosynthesis and Catabolism of Methylxanthines in Plants. In Methylxanthines. Berlin, Heidelberg: Springer Berlin Heidelberg; 2011: 11–31

  57. Treur JL, Taylor AE, Ware JJ, McMahon G, Hottenga JJ, Baselmans BM, Willemsen G, Boomsma DI, Munafo MR, Vink JM. Associations between smoking and caffeine consumption in two European cohorts. Addiction. 2016;111:1059–68.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Warren GW, Alberg AJ, Kraft AS, Cummings KM. The 2014 Surgeon General’s report: “The Health Consequences of Smoking–50 Years of Progress”: a paradigm shift in cancer care. Cancer. 2014;120:1914–6.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Nouran Nada, M.Besher Adi, Aya Fathy, Rushud Abdulsalam, and Rafla Ali for helping in data and sample collection phases of the study. We also thank Dr Hussam Alshraideh for his valuable suggestions on bioinformatics analysis. Sequencing data provided by microbiome insights, Canada.

Funding

The authors were supported by UAE funding agencies; MTA funded by University of Sharjah (Grant code: 1901090253). QH funded by University of Sharjah (Grant code: 150317).

Author information

Authors and Affiliations

Authors

Contributions

Conception and design: MTA, NRD, MSA, and QH. Acquisition of data: ME, DMHA, IMID. Processing of specimens and generation of data: MTA, ME, DMHA, IMID Analysis and interpretation of data: MTA, NRD, ME, DMHA. Drafting or revising of manuscript: MTA, NRD, MSA, QH, ME, DMHA, IMID. Final approval of manuscript: MTA, NRD, MSA, QH, ME, DMHA, IMID. MTA has access to all study data and takes responsibility for the data integrity and accuracy. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mohammad Tahseen Al Bataineh.

Ethics declarations

Ethics approval and consent to participate

All participants in the study read and signed an informed consent; and the Research Ethics Committee at University of Sharjah approved the study protocol, reference # REC-18–02-18–01.

Consent for publication

There are no individual person identifiers in this manuscript. Consent for publication was not sought.

Competing interests

All authors declare they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1.

Figure S1. Oral microbiota community composition of smokers vs. non-smoker groups. Figure S2. Oral microbiota composition based on Fagerström Test for Nicotine Dependence (FTND) score. Figure S3. Effects of smoking and nicotine dependence on oral microbiota richness and diversity.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Al Bataineh, M.T., Dash, N.R., Elkhazendar, M. et al. Revealing oral microbiota composition and functionality associated with heavy cigarette smoking. J Transl Med 18, 421 (2020). https://doi.org/10.1186/s12967-020-02579-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-020-02579-3

Keywords