Research | Open | Published:
The effect of cigarette smoking on the oral and nasal microbiota
Microbiomevolume 5, Article number: 3 (2017)
The goal of the study was to investigate whether cigarette smoking alters oral and nasal microbial diversity, composition, and structure. Twenty-three current smokers and 20 never smokers were recruited. From each subject, nine samples including supra and subgingiva plaque scrapes, saliva, swabs from five soft oral tissue sites, and one nasal swab from both the anterior nares were collected. 16S rRNA V3-V4 region was sequenced for microbial profiles.
We found that alpha diversity was lower in smokers than in nonsmokers in the buccal mucosa, but in other sample sites, microbial diversity and composition were not significantly different by smoking status. Microbial profiles differed significantly among eight oral sites.
This study investigates the effect of cigarette smoking on different sites of the oral cavity and shows a potential effect of cigarette smoking on the buccal mucosa microbiota. The marked heterogeneity of the oral microbial ecosystem that we found may contribute to the stability of the oral microbiota in most sites when facing environmental perturbations such as that caused by cigarette smoking.
Cigarette smoke has adverse effects on human health. Smokers have increased risk of developing diseases such as lung and other cancers, chronic obstructive pulmonary disease, cardiovascular disease, and periodontitis . The microbial communities in the mouth and the nose have direct contact with cigarette smoke and may thus be affected by it. Cigarette smoke contains numerous toxicants to which smokers are regularly exposed on a periodic basis. These toxicants can potentially perturb the microbial ecology of the mouth via antibiotic effects, oxygen deprivation, or other potential mechanisms . The current study examined this question by characterizing the microbiota in eight oral sites and a nasal swab (Fig. 1a) and comparing them between never smokers and current smokers.
The characteristics of the study subjects were shown in Additional file 1: Table S1. No difference in age, gender, race, and alcohol drink were found by smoking status. Smokers had slightly higher periodontal screening and recording (PSR) scores than nonsmokers.
By the Wilcoxon rank-sum test, smokers did not differ significantly from nonsmokers at any site on any measure of alpha diversity, beta diversity, or taxa-relative abundance with one exception—PD_whole tree diversity was lower in the smokers’ buccal mucosa (p = 0.05) (Fig. 1b–d, Additional file 2: Table S2). Based on the t test for the buccal mucosa, smokers had marginally lower observed species (p = 0.046), PD_whole tree (p = 0.032), and nearly significant lower Shannon index (p = 0.074). Thus, it appears that smokers had lower alpha diversity than that in the nonsmokers’ buccal mucosa based on both Wilcoxon rank-sum and t tests. We plotted the mean difference between smokers and nonsmokers for observed species and PD_whole tree (Additional file 1: Figure S1). For most sites, we can exclude large differences in means, but for sub and supragingival plaques, the confidence intervals are wide, reflecting small sample size. Thus, we cannot exclude the possibility of a substantial smoking effect at these two sites. Moreover, Additional file 2: Figure S1 indicates that alpha diversity in the buccal mucosa is lower in smokers as described above.
Regardless of the smoking status, microbiota among the oral sites differed greatly (Fig. 2, Additional file 2: Table S3). We found no difference in any microbial measurement by age, gender, race, alcohol consumption, and PSR score (data not shown).
Our study reported an effect of cigarette smoke on the oral microbiota in the buccal mucosa. Although our sample size is modest, it is remarkable that overall, we found no associations between microbial features and smoking status in other oral sites. Our results for the buccal mucosa should be examined in other studies since we tested multiple sites. Previous studies have shown inconsistent results regarding the effect of smoking on the oral microbiota. For example, an altered microbiota composition was noted in studies of 62 swabs from the tonsillar pillars , 200 subgingival samples , 30 marginal and subgingival plaque and gingival crevicular fluid samples , and oral wash samples comparing 13–41 current smokers with 77–194 never smokers in 4 groups  and 64 saliva samples  by sequencing 16S rRNA gene and in a study of 292 stimulated oral samples by Human Oral Microbe Identification Microarray (HOMIM) . These studies, however, reported that different taxa were changed by smoking status. Some studies [9–11] reported no association between oral microbiota and smoking. Heterogeneity in smoking definitions, oral sites sampled, study population, and inclusion/exclusion of participants with related conditions (i.e., periodontitis), and in methods for collecting the specimens, likely account for some of the inconsistency with respect to smoking. In addition, our results suggested that smoking did not affect the oral microbiota with a large effective size. Large studies capable of exploring smaller magnitude effects of smoking on specific oral niches are needed to fully understand subtle smoking-induced alteration in microbiota.
Data for smoking effects on the nasal microbiota are sparser and even less clear. In contrast to our null results, swabs from the left and right nasopharynx in 62 subjects found significant differences in relative abundance of diverse taxa by smoking status, but 55% of these smoking-associated taxa were found only in the left or right nasopharynx, suggesting inconsistency .
Differences of microbiota across oral sites, regardless of smoking, point to distinct microbial niches. Consistent with our findings and with the Human Microbiome Project (HMP) , a recent study of 66 Chinese subjects showed significant differences in microbiota composition among the buccal mucosa, saliva, and dental plaque sites . Different oral sites harbor a distinct microbiota, suggesting that the various oral surfaces might provide different receptors for bacterial adhesion, species interactions, and environmental conditions (e.g., oxygen level) for microbial survival and growth.
Ecological studies have suggested that more diverse and complex communities are generally more resistant and resilient to perturbations . In the HMP, the oral microbiota, despite its heterogeneity across subsites, was more stable over time than the other body site evaluated . As shown here, the oral cavity is a highly heterogeneous ecological system containing distinct microbial niches. In addition, it is well accepted that microbial organisms living on oral surfaces can switch from a free-living state to a sessile mode in biofilms with advantageous properties that include enhanced tolerance to many adverse conditions including antimicrobial agents . Consistent with this view, a recent study of 66 healthy subjects with saliva samples at baseline, immediately and 1, 2, 4, and 12 months after antibiotics usage found that saliva microbial composition remained stable . Likewise, our findings suggest that the microbial communities in the oral cavity might be resistant or resilient to disturbances such as cigarette smoking.
Strengths of the current study include extensive sampling of the oral cavity with a well-developed protocol, careful matching of case and control groups, a contrast of never with heavy smokers, state-of-the-art assays including quality control samples for batch effects, and rigorous statistical analysis. These methods should have identified significant associations with smoking status in the most sampled sites. Weaknesses include the modest sample size, especially for subgingival plaque samples due to sequencing failure. The study only allows us to detect the effect of smoking with a large effect size on the oral microbiota. Subtle effects or temporal effects of smoking could not be examined in this study.
Our study showed that cigarette smoking had a significant effect on the microbiota of the buccal mucosa, but not in other oral sites and nasal cavity. The oral cavity is heterogeneous with distinct communities across sites, which may contribute to its stability in the face of potential perturbing factors such as smoking. Larger studies are needed to further examine oral and nasal microbiota effects produced by behaviors such as cigarette smoking.
Following the approval by the institutional review boards of the National Cancer Institute and University of Rochester, 23 current smokers (median duration, 15 years; median intensity, 15 cigarettes per day) and 20 nonsmokers (<100 cigarettes in a lifetime) were recruited at Eastman Institute of Oral Health, University of Rochester. All the subjects signed informed consent and filled out questionnaires. Individuals with antibiotic usage or professional dental cleaning within the last 3 months or diagnosed with periodontal disease or cancer or losing >1 tooth were excluded. Participants were screened by PSR index estimated at the time of recruitment. Groups were frequency-matched for gender and race (Additional file 1: Table S1).
From each subject, we sought nine samples, including supra and subgingiva plaque, saliva, swabs from five soft tissue sites, and one nasal swab from both the anterior nares. Figure 1a shows the detailed sample collection locations in the oral cavity. The samples were collected by following the procedure of Human Microbiome Project (http://hmpdacc.org/doc/HMP_MOP_Version12_0_072910.pdf).
16S rRNA gene sequence analysis
The DNA was extracted from samples as described previously . The V3-V4 regions of the 16S rRNA gene were amplified and sequenced on an Illumina MiSeq instrument using the 300 paired-end protocol at the Institute of Genome Sciences, Genomic Resource Center, University of Maryland School of Medicine . The sequence data were submitted to NCBI BioProject (http://www.ncbi.nlm.nih.gov/bioproject) under accession number PRJNA316469.
Sequence reads were processed to remove low quality and short reads (see details in ). The remaining reads (18,497 ± 14,130 reads/sample) were clustered into Operational Taxonomy Units (OTUs) at 97% identity using the command pick_open_reference_otus.py and Greengenes database as the reference (version 13_8)  in Quantitative Insights into Microbial Ecology (QIIME 1.8.0) . The default parameters were used except the method of usearch61 and percent_subsample of 0.1. OTUs with only one read were excluded.
Alpha diversity was estimated as the number of OTUs (Observed_species), Shannon’s Index (using information of OTU frequency) , and phylogenetic diversity (using information of phylogenetic relationship of different OTUs) (PD_whole_tree) ) by averaging over 20 rarefied tables (1000 reads/sample). Taxonomic beta diversity was measured as unweighted and weighted UniFrac distance based on the OTU table . Relative abundance of taxa was calculated from unrarefied OTU table. To rule out batch effects, 19 random samples were selected to examine the difference within and between batches, and no difference was found in alpha and beta diversity measures (Additional file 1: Figure S2).
The Wilcoxon rank-sum test was used to examine gastric microbiota alpha diversity and taxa-relative abundance differences between groups. For the buccal mucosa, 95% confidence intervals for differences in mean alpha diversity were based on t-statistics. The Spearman correlation was used to examine the correlation between continuous variables. Bonferroni correction was used to adjust for tests of multiple taxa. Permutational multivariate analysis of variance (PERMANOVA, adonis) was used to examine the association between unweighted/weighted UniFrac distance and smoking status/other demographic variables. P values less than 0.05 were considered significant after adjustment for multiple tests.
Human microbiome project
Operational taxonomical unit
Permutational multivariate analysis of variance
Periodontal screening and recording
Quantitative insights into microbial ecology
Staempfli MR, Anderson GP. How cigarette smoke skews immune responses to promote infection, lung disease and cancer. Nat Rev Immunol. 2009;9:377–84.
Macgregor ID. Effects of smoking on oral ecology. A review of the literature. Clin Prev Dent. 1989;11:3–7.
Charlson ES, Chen J, Custers-Allen R, Bittinger K, Li HZ, Sinha R, et al. Disordered microbial communities in the upper respiratory tract of cigarette smokers. PLoS One. 2010;5.
Mason MR, Preshaw PM, Nagaraja HN, Dabdoub SM, Rahman A, Kumar PS. The subgingival microbiome of clinically healthy current and never smokers. Isme Journal. 2015;9:268–72.
Kumar PS, Matthews CR, Joshi V, de Jager M, Aspiras M. Tobacco smoking affects bacterial acquisition and colonization in oral biofilms. Infect Immun. 2011;79:4730–8.
Wu J, Peters BA, Dominianni C, Zhang Y, Pei Z, Yang L, et al. Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 2016.
Morris A, Beck JM, Schloss PD, Campbell TB, Crothers K, Curtis JL, et al. Comparison of the respiratory microbiome in healthy nonsmokers and smokers. Am J Respir Crit Care Med. 2013;187:1067–75.
Belstrom D, Holmstrup P, Nielsen CH, Kirkby N, Twetman S, Heitmann BL, et al. Bacterial profiles of saliva in relation to diet, lifestyle factors, and socioeconomic status. J Oral Microbiol. 2014;6.
Bostrom L, Bergstrom J, Dahlen G, Linder LE. Smoking and subgingival microflora in periodontal disease. J Clin Periodontol. 2001;28:212–9.
Gomes SC, Nonnenmacher C, Susin C, Oppermann RV, Mutters R, Marcantonio RAC. 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.
Haffajee AD, Socransky SS. Relationship of cigarette smoking to the subgingival microbiota. J Clin Periodontol. 2001;28:377–88.
Zhou Y, Gao H, Mihindukulasuriya KA, Rosa PS, Wylie KM, Vishnivetskaya T, et al. Biogeography of the ecosystems of the healthy human body. Genome Biol. 2013;14:R1.
Xu X, He J, Xue J, Wang Y, Li K, Zhang K, et al. Oral cavity contains distinct niches with dynamic microbial communities. Environ Microbiol. 2015;17:699–710.
McCann KS. The diversity-stability debate. Nature. 2000;405:228–33.
Hoiby N, Bjarnsholt T, Givskov M, Molin S, Ciofu O. Antibiotic resistance of bacterial biofilms. Int J Antimicrob Agents. 2010;35:322–32.
Zaura E, Brandt BW, Teixeira de Mattos MJ, Buijs MJ, Caspers MP, Rashid MU, et al. Same exposure but two radically different responses to antibiotics: resilience of the salivary microbiome versus long-term microbial shifts in feces. MBio. 2015;6:e01693–15.
Fadrosh DW, Ma B, Gajer P, Sengamalay N, Ott S, Brotman RM, et al. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome. 2014;2:6.
Yu G, Fadrosh D, Goedert JJ, Ravel J, Goldstein AM. Nested PCR biases in interpreting microbial community structure in 16S rRNA gene sequence datasets. PLoS One. 2015;10, e0132253.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Shannon CE. The mathematical theory of communication. 1963. MD Comput. 1997;14:306–17.
Faith DP, Baker AM. Phylogenetic diversity (PD) and biodiversity conservation: some bioinformatics challenges. Evol Bioinforma. 2006;2:121–8.
Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011;5:169–72.
We thank Bing Ma in Ravel lab for the initial quality control filtering of the sequencing data and for delivering it to us.
This work was supported by the Intramural Research Program of the National Cancer Institute at the National Institutes of Health.
Availability of data and materials
All sequence data and samples’ information are available at NCBI under BioProject ID PRJNA316469.
GY and NEC designed the study. SP and YR collected the samples. MSH and JR performed the sequencing. GY analyzed the data and drafted the manuscript. MHG supervised the statistical analysis. GY, MHG, JJG, JR, and NEC did the data interpretation and manuscript revision. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
The consent forms signed by each participant included their consent to allow us to publish our findings.
Ethics approval and consent to participate
This study was approved by the institutional review boards of the National Cancer Institute and University of Rochester. All participants signed informed consent to participate in the study.
Summary of study subjects by smoking status. Figure S1. Mean difference between smokers and nonsmokers in observed species (a) and PD_whole_tree (b) with 95% confidence interval based on t-distribution. The black dots are mean difference between smokers and nonsmokers (smokers, nonsmokers). Figure S2. Box plot showing no difference of within- and between-plate/batch variation in alpha diversity (Shannon) and beta diversity (unweighted and weighted UniFrac). Boxes are interquartile range (IQR), median values are bands within the boxes, lines outside the boxes are 1.5-times IQR, and dots are outliers. (DOC 391 kb)
Comparison by smoking status in alpha diversity and relative abundance of top ten abundant genera. Table S3. P values based on the Wilcoxon signed-rank test for the Shannon index differences among eight oral sites. (XLS 50 kb)