Next Article in Journal
Lung Function Assessment as an Early Biomonitor of Mercury-Induced Health Disorders in Artisanal and Small-Scale Gold Mining Areas in Indonesia
Next Article in Special Issue
A Survey of Health Disparities, Social Determinants of Health, and Converging Morbidities in a County Jail: A Cultural-Ecological Assessment of Health Conditions in Jail Populations
Previous Article in Journal
Income Changes Due to Disability Ratings and Participation in Economic Activities Caused by Industrial Accidents: A Population-Based Study of Data from the Fourth Panel Study of Workers’ Compensation Insurance (PSWCI)
Previous Article in Special Issue
Using a Virtual Community (the Health Equity Learning Collaboratory) to Support Early-Stage Investigators Pursuing Grant Funding
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 15
Issue 11
10.3390/ijerph15112479
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cigarette Smoking Modulation of Saliva Microbial Composition and Cytokine Levels

by
Mary Rodríguez-Rabassa
1,2,
Pablo López
1,
Ronald E. Rodríguez-Santiago
1,
Antonio Cases
3,
Marcos Felici
3,
Raphael Sánchez
1,
Yasuhiro Yamamura
1 and
Vanessa Rivera-Amill
1,*
1
AIDS Research Infrastructure Program, Ponce Research Institute, Ponce Health Sciences University, 395 Dr. Luis F. Sala Street, Ponce, PR 00716-2348, USA
2
Clinical Psychology Program, School of Behavioral & Brain Science, Ponce Health Sciences University, Ponce, PR 00716-2348, USA
3
Tobacco Control and Oral Health Division, Department of Health, Commonwealth of Puerto Rico, San Juan, PR 00716-2348, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(11), 2479; https://doi.org/10.3390/ijerph15112479
Submission received: 19 September 2018 / Revised: 26 October 2018 / Accepted: 3 November 2018 / Published: 7 November 2018
(This article belongs to the Special Issue Proceedings of Research Centers at Minority Institutions (RCMI) Translational Science 2017 Conference)
Browse Figures
Versions Notes

Abstract

:
Tobacco use has been implicated as an immunomodulator in the oral cavity and contributes to the development of oral cancer. In the present study, we investigated the effects of cigarette smoking on bacterial diversity and host responses compared to healthy nonsmoking controls. Saliva samples were collected from eighteen smokers and sixteen nonsmoking individuals by passive drool. The 16S rRNA gene was used to characterize the salivary microbiome by using the Illumina MiSeq platform. Cytokine and chemokine expression analyses were performed to evaluate the host response. Significant differences in cytokine and chemokine expression levels of MDC, IL-10, IL-5, IL-2, IL-4, IL-7, adrenocorticotropic hormone (ACTH), insulin, and leptin were observed between smokers and nonsmokers. Taxonomic analyses revealed differences between the two groups, and some bacterial genera associated with the smokers group had correlations with hormones and cytokines identified as statistically different between smokers and nonsmokers. These factors have been associated with inflammation and carcinogenesis in the oral cavity. The data obtained may aid in the identification of the interactions between the salivary microbiome, host inflammatory responses, and metabolism in smokers.

1. Introduction

Although policies to control cigarette smoking in the population have been successfully implemented in Puerto Rico, it is still an important public health problem on the island [1]. In 2016, the prevalence of adults who smoked in Puerto Rico was 10.6% [2]. Tobacco use is a key modifiable behavioral risk factor associated with chronic diseases [3]. Cigarette smoke has been linked to many common medical conditions such as cancer, coronary heart disease, chronic obstructive pulmonary diseases, and periodontitis [4,5,6,7]. Recent studies have established that cigarette smoking has a significant effect on the microbiota and cytokines expression of the buccal mucosa [8,9,10,11,12]. These factors have been associated with inflammation and carcinogenesis in the oral cavity and other body compartments [13,14]. The incidence rate of oral cancer in Puerto Rico is high and may be related in part to tobacco use [15,16]. Interestingly, men are still more likely to smoke (32.1%) than women (25.1%) [17], which coincides with a higher incidence of oral cavity and pharyngeal carcinomas in men in the Puerto Rican population [18]. Oral cancer is the eighth most prevalent cancer worldwide and constitutes a significant global public health problem [19,20]. Oral cancers have multifactorial origins such as smoking cigarettes, heavy alcohol consumption, bacterial translocation, and human papilloma virus infection (HPV), among others [21,22,23].
Several recent studies have suggested that environmental factors (e.g., diet and smoking) and host susceptibility (e.g., stress, anxiety, and depression) affect the equilibrium of the oral microbiome (dysbiosis) and expand the growth of oral pathogens that may have effects on other parts of the body [24,25,26]. For example, cancer in digestive organs (colorectal and pancreatic) has been associated with oral microflora (Porphyromonas gingivalis, Fusobacterium nucleatum), and in many cases with smoking [27,28,29,30,31]. In addition, microbial composition influences the development and function of innate and adaptive branches of the immune system [32].
The principal objective of the study was to analyze the relationship between cigarette smoking, salivary microbiome composition, and host response among adult Puerto Rican smokers by using high-throughput 16S rRNA gene sequencing and Luminex xMAP technology. This information may aid in the identification of the interactions between the salivary microbiome, host immune responses, and metabolic processes in smokers.

2. Methods

The current study was conducted following the Declaration of Helsinki, and the protocol was approved by the Institutional Review Board of the Ponce Medical School Foundation, Inc. (IRB approval no: 150209-YY). All participants signed their informed consent before sample collection and completion of study questionnaires.

2.1. Study Subjects

This cross-sectional study involved the collection of participants’ saliva and depression symptomatology as described below. Thirty-four subjects were recruited for this study, 16 nonsmokers (sex: 6 males, 10 females) and 18 current smokers (sex: 10 males, 8 females). Fifty percent of the current smokers were recruited from the Puerto Rico Department of Health Quitline (n = 9) and 50% from the Puerto Rican general population (n = 9). All nonsmokers were recruited from the Puerto Rican general population.

2.2. Depression Symptomatology

The Patient Heath Questionnaire-9 (PHQ-9) was used to assess depression symptoms. This 9-item self-report scale evaluates how often the described symptom bothered the participant over the 2 weeks prior to the test being filled out. Each of the items can be scored from 0 (not at all) to 3 (nearly every day). The general score can range from 0 to 27 and is interpreted as follows: a score of 0–4 means that the subject has minimal or no symptoms, 5–9 signifies mild depression, 10–14 denotes moderate depression, 15–19 suggests moderately severe depression, and 20–27 signifies that the subject has severe depression. Its validity and reliability as a diagnostic measure, as well as its utility in assessing depression severity and monitoring treatment response, are well-established [33,34,35,36,37].

2.3. Sample Collection, DNA Extraction, and Amplification

Approximately 1 mL of saliva was collected by passive drool into an UltraSal-2 saliva collection device (Oasis diagnostics, Vancouver, WA, USA). The samples were aliquoted and stored at −20 °C until processing. Bacterial DNA was extracted from saliva using the QIAamp DNA Microbiome Kit (Qiagen, Hilden, Germany) according to the manufacturer’s recommendations. The conserved 16S V3-V4 region was amplified using PCR protocol. We performed the 16S rRNA gene amplification by using FastStart master mix from Roche (Roche Diagnostics, Mannhein, Germany). The forward and reverse primers used for the amplification were 16F (5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3′) and 16R (5′-GTCTCGGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVHHHTATCTAATCC-3′), respectively [38]. The PCR conditions used were as follows: A denaturing step at 95 °C for 3 min, followed by 25 cycles at 95 °C for 30 s, 55 °C for 30 s, and 78 °C for 30 s, and a final extension at 72 °C for 5 min. In order to detect possible contamination during the procedures, the amplifications were carried out with an appropriate negative control. Each amplicon was confirmed by 1.5% agarose gel.

2.4. Illumina MiSeq Sequence Determination by Nextera XT

After PCR amplification, we used a Qubit dsDNA HS assay kit and a Qubit fluorometer (Life Technologies, Carlsbad, CA, USA) to quantify the concentration of the amplicons. The Nextera XT DNA Sample Preparation Kit (Illumina, Sand Diego, CA, USA) was used according to the manufacturer’s instructions to generate paired-end DNA libraries. Briefly, the PCR products were fragmented and tagged with an adapter, and the library was purified by using AMPure XP beads (Beckman Coulter, Indianapolis, IN, USA). A normalization procedure was performed to ensure equal library representation during sequencing. Each sample was pooled into a 1.5 mL tube and heated at 96 °C for 2 min. The denatured library was placed on ice for 5 min; PhiX (12.5 pM) library control was added to provide quality control checks. The amplicon library was loaded onto the Illumina MiSeq flow cell. A sequencing run of 2 × 250 bp MiSeq paired-end reads was performed. The data obtained were base-called, and reads with the same barcode were collected and assigned to a sample on the instrument. After the procedure, the software generated analysis output in the FASTQ file format.

2.5. Metagenomics Data Quality Control and Analysis

MiSeq data was extracted, decompressed, and analyzed using the QIIME software (v.1.9.0) (University of Colorado, Boulder, CO, USA) [39]. Forward and reverse reads (FASTQ file) of each sample were joined using the join_paired_ends.py script with QIIME default parameters, and the quality of sequences was assessed using FastQC software (v.10.1, Babraham Institute, Cambridge, UK) (www.bioinformatics.babraham.ac.uk/projects/fastqc/). The resulting file containing all sequences was length and quality trimmed by discarding sequences shorter than 75 nucleotides and a quality threshold of 20 using Trim Galore! software (Babraham Institute, Cambridge, UK) (htps://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). The resulting files were converted to FASTA using the fastq_to_fasta script (FASTX-Toolkit, Hannon Lab, Cold Spring Harbor, NY, USA, http://hannonlab.cshl.edu/fastx_toolkit/license.html) and then used to determine the taxonomy composition of each sample (pick_closed_reference_otu.py script) with a 97% similarity threshold using the Greengenes reference dataset (gg_otus_13_8-release). The beta diversity analyses were produced with the core_diversity_analyses.py from QIIME [40]. The rarefaction of the Operational Taxonomic Unit (OTU) table for alpha diversity was 4394. We also calculated beta diversity with unweighted UniFrac, as well as microbe relative abundance files [41,42]. The OTU table was rarefied to the same number of OTUs using the single_rarefication.py QIIME command. The rarefied OTU table was used as input for Explicet (v.2.10.5, University of Colorado, Boulder, CO, USA) to create stacked bar charts and other statistical data [43].
The sequence data were submitted to NCBI BioProject (http://www.ncbi.nlm.nih.gov/bioproject) under accession number PRJNA407051.

2.6. Cytokine Assay

Proinflammatory and anti-inflammatory cytokine/chemokine concentrations in saliva were calculated by using a Milliplex Human Cytokine/Chemokine Magnetics Bead Panel (TNF-α, IL-12 (p70), MDC, IL-10, IFN-γ, TNF-β, IL-1β, IL-5, IL-2, IL-6, IL-4, IL-1RA, IL-13, IL-17, IL-7, and GM-CSF) and a Human Bone Magnetic Bead Panel (insulin, leptin, and adrenocorticotropic hormone (ACTH)). The procedures were performed according to the manufacturer’s recommendations (Millipore, Billerica, MA, USA). A Luminex MAGPIX instrument with xPonent software version 4.2 (Luminex, TX, USA) and a Milliplex Analyst version 3.5.5.0 (VigeneTech, Carlisle, MA, USA) were used to collect the data and for analyses, respectively. Briefly, a minimum of 50 microsphere events were acquired from each analyte, and cytokine/chemokine concentrations were calculated by measuring the median fluorescent intensity (MFI) in a 5-parameter curve-fitting method [44]. Measurements were performed in duplicate.

2.7. Microbial Biomarkers and Correlation Model

To determine smoker microbial biomarkers, we used LEfSe [45] to compare the bacterial taxonomic composition of the smoking and nonsmoking groups. The taxonomy was evaluated at the genus level, and it was given a p-value of 0.05 to be considered significant and a minimum threshold of 3.0 in the linear discriminant analysis (LDA). To understand correlations between the bacteria, cytokines, and hormones, we generated a correlation network. The correlations and the network were generated using qgraph in R 3.4.0 (University of Amsterdam, Amsterdam, The Netherlands) [46].

2.8. Statistical Analysis

The statistical significance between nonsmoker and smoker variables (depression and cytokine expression) were performed using a nonparametric Mann-Whitney-Wilcoxon test implemented in SPSS version 20 (SPSS Inc., Chicago, IL, USA). The level of statistical significance was set at p < 0.05. The tests of significance of unweighted UniFrac distance were performed using a two-sided Student’s two-sample t-test. The nonparametric p-values (Bonferroni-corrected) were calculated using 999 Monte Carlo permutations implemented by QIIME (University of Colorado, Boulder, CO, USA) [39]. A Spearman correlation test was used for the correlation network.

3. Results

3.1. Demographic Parameters of Study Participants

A total of 34 participants were enrolled in the study, including 16 nonsmokers and 18 current smokers. None of the participants were taking antibiotics. The participants from the smoker group mostly had over three years of smoking history (83.3%). Half of the participants reported light/moderate cigarette consumption (<20 cigarettes per day), whereas the others reported being heavy cigarette users (≥20 cigarettes per day). Demographic data showed that the individuals in the smoker group were predominantly male (55.6%), whereas the nonsmoker group was mostly female (62.5%). The mean age for individuals in the smoker group was 45 years (range 33–49), whereas in the nonsmoker group it was 34 years (range 28–50). The majority of the participants in the smoker group (67.0%), as well as in the nonsmoker group (75.0%), had a college education. The participants from the nonsmoker group showed higher alcohol use (68.8%) than those from the smoker group (38.9%). No statistically significant differences were detected for age, sex, education level, and alcohol use (Table 1). However, the PHQ-9 mean scores showed significantly increased (p = 0.02) depression symptomatology in smokers when compared to nonsmokers (Table 1).

3.2. Microbial Sequencing and Bacterial Diversity

According to the beta diversity of the unweighted UniFrac analysis, the smoker and nonsmoker groups were significantly different (p-value < 0.05). The salivary microbial population between the two groups differed by the most abundant bacteria (Figure 1A). In addition, the intraindividual variability was greater in smokers than in nonsmokers (Figure 1B). The five phyla that were more abundant in nonsmokers’ saliva samples were Firmicutes (66%), Bacteroidetes (16%), Actinobacteria (5%), Fusobacteria (5%), and Proteobacteria (4%), which was 96% of all sequences. In contrast, Proteobacteria (40%), Firmicutes (29%), Bacteroidetes (23%), Fusobacteria (5%), and Actinobacteria (2%) dominated in smoker samples and represented 99% of all sequences (Figure 2).
At the genus level, Streptococcus was most prevalent in both groups, although at a lower percentage in the smoker group (Figure 3). The bacterial compositions (the top ten) in the nonsmoker group samples were as follows: Streptococcus (35%), Veillonella (10%), Prevotella (8%), Porphyromonas (5%), Staphylococcus (5%), Actinomyces (3%), Leptotrichia (2%), Fusobacterium (2%), Granulicatella (2%), and Haemophilus (2%). The bacterial compositions in the smoker group samples differed from the nonsmoker group: Streptococcus (15%), Haemophilus (14%), Prevotella (13%), Neisseria (13%), Porphyromonas (9%), Veillonella (6%), Fusobacterium (5%), Aggregatibacter (2%), Staphylococcus (2%), and Actinobacillus (1%). The highest significance levels (p ≤ 0.002) of Streptococcus, Staphylococcus, Actinomyces, Granulicatella, Fusobacterium, and Leptotrichia were in nonsmoker samples. Actinobacillus was highest in smokers (p = 0.002). A tendency of increased abundance of Porphyromonas, Neisseria, Haemophilus, Prevotella, and Aggregatibacter was observed in smokers. The inner-group variance in smokers was higher than the comparison of the variance between smokers and nonsmokers.

3.3. Cytokines

We used a multiplex assay from Millipore to examine the expression patterns of pro- and anti-inflammatory cytokines in nonsmoker and smoker saliva samples from study participants. According to our data, IL-4, IL-2, and ACTH (p = 1.6 × 10−4, p = 9.07 × 10−10, and p = 5.83 × 10−4, respectively) were significantly higher in smoker samples. However, MDC, IL-10, IL-5, IL-7, insulin, and leptin were downmodulated compared to nonsmoker samples (p = 0.007, p = 4.28 × 10−4, p = 0.001, p = 0.050, p = 3.66 × 10−4, and p = 6.77 × 10−4, respectively). No significant difference was observed in the expression of TNF-α, IL-12 (p70), IFN-γ, TNF-β, IL-1 β, IL-6, IL-1ra, IL-13, IL-17, IL-7, and GM-CSF (Table 2).

3.4. Correlation Network

The enrichment of bacteria was done using LEfSe. The parameters used were p < 0.05, and the linear discriminant analysis was at least 3.0. In contrast to the genera presented in Figure 3, which were the most abundant genera, the variance of the bacteria abundance was compared between the two groups and then associated to one of the groups by linear discrimination. The LEfSe analysis resulted in 26 genera being different. Out of the 26 genera, 12 were associated with the smoker group, and 14 genera were associated with the nonsmoker group (Figure 4A). To observe the relationships between the hormones, cytokines, and bacteria identified as enriched in the smoker group, we generated a correlation network using the qgraph R package. When the smoker group-associated bacteria were combined with the cytokines and hormones in networks, we found five discrete correlation networks (Figure 4B). Interestingly, only one network contained bacteria, cytokines, and hormones, whereas the others were just one or two types of the aforementioned factors. The network that contained the most variables had two bacterial genera (Necropsobacter and Actinobacillus), four cytokines (IL-2, IL-5, IL-10, and MDC), and three hormones (insulin, leptin, and ACTH). It demonstrated that bacteria have direct correlations with other bacteria and cytokines, but indirect correlations with hormones. Thus, cytokines were the factors that were able to interact among themselves and with bacteria and hormones in a direct manner. The only negative correlations were between leptin and insulin, and between IL-10 and IL-5.

4. Discussion

Cigarette smoking likely modulates the microbial composition of the oral cavity. Prior studies have confirmed that tobacco smoke disrupts the homeostasis of commensal microbial composition, gingival diseases, and dental alteration [47,48]. In the present study, the use of tobacco cigarettes had a significant impact on the relative abundance of the three more prevalent bacterial phyla (Firmicutes, Proteobacteria, and Bacteroidetes) found in the saliva of samples collected in Puerto Rico. Recent analyses indicate that Firmicutes are a phylum more prevalent in saliva [49,50]. In our nonsmoker samples, Firmicutes were the most abundant phylum, while in smoker samples, Proteobacteria were the most abundant. Similarly, a significant increase in Bacteroidetes was observed in smoker samples. These findings suggest that smoking has a critical impact on oral cavity homeostasis. The positive effects of the host-microbiome symbiosis, which supports host defense functions and maintains a healthy digestive tract, may be altered by the microenvironment changes caused by the variety of organic and inorganic chemical compounds present in cigarette smoke [51,52]. Members of these phyla include some of the bacteria identified in oral diseases such as caries and periodontitis, among others. While Proteobacteria are more frequently associated with intestinal diseases, some studies indicate an association with cardiovascular disease, chronic obstructive pulmonary disease (COPD), asthma, and metabolic conditions such as diabetes [53,54].
At the genus level, the impact of cigarette smoke was observed as a lower abundance of Streptococcus in smoker samples. In the healthy oral cavity, species of Streptococcus (Firmicutes), Actinomyces (Actinobacteria), and Neisseria (Proteobacteria) are considered to be normal flora that establish a cooperative relationship with the host [55]. In contrast, some species of Actinobacillus (Proteobacteria), Prevotella (Bacteroidetes), and Porphyromonas (Bacteroidetes) are most often involved in the etiology of periodontitis and other oral diseases [56]. Our results showed that the frequency of these three genera in our smoker samples was significantly higher when compared to nonsmokers, suggesting an increased probability of oral infection and inflammation in the smoker participants. Additionally, recent studies have suggested that Actinobacillus, Prevotella, and Porphyromonas are involved in the dysregulation of cytokine networks, damage in the crevicular epithelium, leukocyte killing, the breakdown of periodontal tissues, and an increase in mucosal permeability [57,58]. The normal flora prevents oral diseases by inhibiting the adherence of pathogens onto specific surfaces, in turn preventing significant changes in the relative abundance of microbes and decreasing virulence factors [25]. Understanding that the oral cavity is one of the most important gateways of microbes into the body, we can suggest that the adverse effect of cigarette smoke on the oral microbiome will negatively impact the local and systemic compartments.
Depression symptomatology was significantly increased in the smoker group compared to the nonsmoker group. This association is consistent with previous studies, revealing that smoking and exposure to cigarette smoke can lead to the development of depression [59,60]. In addition, individuals with mental illnesses tend to smoke more than those without mental illnesses [61]. Several common pathways may explain the comorbidity of nicotine dependence and depression, including neurotransmitters such as serotonin and dopamine and the inflammatory pathways (proinflammatory cytokines). Decreased serotonin function has been related to cigarette smoking [62], and reduced serotonergic receptors to depression [63]. The nicotine in cigarettes stimulates the production of dopamine [64], which is a crucial modulator of mood and depression [65]. Additionally, both smoking and depression are associated with inflammation [66,67].
We also wanted to determine the immune response markers within the smoker group. The cytokines identified as different in the smoker group corresponded to a T-cell mediated immune response. The response was not related to one single type of T-cell, as IL-2, IL-5, and IL-10 corresponded to TH1, TH2, and Tregs. The TH1 and TH2 responses have been reported to increase pulmonary inflammation, asthma, and allergic inflammation [68,69,70,71]. It should be noted that even though IL-5, IL-10, MDC, insulin, and leptin had lower means in the smoker group, it is interesting that they maintained the correlations that have been reported in other studies. An interpretation of this could be that there was an inflammatory response, which could then have turned into a severe dysbiosis toward metabolic syndrome, as seen in the kwashiorkor and tumor genesis scenarios [72,73]. A longitudinal study of a similar cohort could detect if the correlations of the variables indicated above may be used as early and long-term markers of metabolic syndrome in the cigarette smoking population.
When comparing the most abundant genera and the LEfSe association of the genera to the smoker and nonsmoker groups, Necropsobacter was identified by LEfSe as associated with the smoker group. A correlation analysis revealed that IL-2, IL-5, IL-10, MDC, insulin, leptin, ACTH, Necropsobacter, and Actinobacillus were associated within the same network. Previous work has identified IL-2, IL-5, IL-10, and MDC to be modulated by smoking [74,75,76,77,78]. Among them, IL-2 is the only cytokine that increased in the smoker group when compared to nonsmokers. A study on insulin-dependent diabetes mellitus patients revealed that there is a decreased production of IL-2 when compared to controls [79], which could indicate that the observed correlation of this cytokine with insulin in the smoker group in the present study was not a sign of insulin-dependent diabetes mellitus. In addition, the positive relationship observed between IL-2 and ACTH has been documented to play a role in depressive disorders [80,81]. Cytokine IL-5 levels were decreased in smokers, and the levels positively correlated with insulin. This observation has been detected in adipose tissue, where insulin resistance was correlated with decreased levels of IL-5 [82]. The model also positively correlated IL-5 with IL-10. The relationship between IL-5 and IL-10 has been documented as a positive one in the scenario of acute stress, in which the inflammatory profile shifts toward an anti-inflammatory state [83]. Another cytokine that had reduced expression when compared to nonsmokers was IL-10. One of the hormones that correlated with IL-10 in our model was leptin. The positive correlation of IL-10 with leptin has been documented in countering the progression of obesity and metabolic syndrome and in ameliorating inflammation [84,85,86].
Previous studies have found that smoking and smoking cessation are associated with insulin resistance [87,88,89,90], which is in accord with lower insulin levels in the smoker group. The quantification of oral insulin is supported by a previous study that measured oral insulin in smokers and demonstrated an increased risk of developing type 2 diabetes [91]. In addition to insulin resistance, tobacco induces leptin resistance [92]. A study by Cnop established the relationship between intra-abdominal fat, insulin resistance, and leptin. It demonstrated that with obesity, there is an increased risk of leptin and insulin resistance [93]. In our correlation model, leptin demonstrated a negative correlation with insulin, which has been shown to play a role in the modulation of insulin resistance in subjects undergoing a weight loss program [94]. Patients with chronic periodontitis also showed lower salivary leptin concentrations when compared to healthy volunteers [95]. Even though leptin appeared to have a lower mean in smokers when compared to nonsmokers, there was a correlation with insulin, which has been demonstrated in smokers and long-term nicotine gum users [96]. We also observed higher levels of adrenocorticotropic hormone in the smoking group. ACTH has been previously used as a marker for stress and has been positively correlated with cortisol and smoking cravings [97,98]. ACTH has also been observed to be elevated due to exposure to cigarette smoke, and could lead to metabolic syndrome [99]. Tobacco and cardiac function have been associated with changes in angiotensin-converting enzyme and ACTH levels [100].
Both Necropsobacter and Actinobacillus belong to the Pasteurellaceae family, meaning that they are closely related in terms of phylogenetic classification [101,102]. Isolation of the Necropsobacter genus in humans has been reported in clinical case studies related to bacteremia and abscesses [103], whereas Actinobacillus has been reported in the pathogenesis of periodontitis [104]. Interestingly, smoking has been linked to abscesses and periodontitis [105,106], suggesting that both Necropsobacter and Actinobacillus may have a role in the development of smoking-related abscesses.
Our study had certain limitations. First, this study had a small sample size. A larger cohort would be essential to provide stronger statistical power and to obtain robust evidence on the associations identified in our study. Second, we did not control for other possible confounders (e.g., oral health, hygiene, and physical properties of saliva). Further studies are needed to determine temporal changes in microbiome communities and immune response. Longitudinal studies will allow for identifying the time dependence of the results and the development of chronic diseases such as cancer and cardiovascular disease, among others.

5. Conclusions

While the study design (cross-sectional) limited the results and did not allow for the temporal association between cigarette smoking, depression, and microbiome composition, our study showed that cigarette smoking altered the expression levels of cytokines, chemokines, and hormones. It also demonstrated that the salivary microbiome differed between nonsmokers and smokers. Such differences may be a result of toxicants found in cigarette smoke, which can have an impact on host defense and metabolism and has implications for chronic diseases. The correlation analysis indicated significant correlations between oral hormones, cytokines, and bacteria, suggesting that the data obtained in this study support the use of noninvasive biomarkers in the identification of high-risk factors for the development of chronic diseases.

Author Contributions

M.R.-R., P.L., A.C., M.F., and Y.Y. designed the study. R.S. acquired the sequencing data. M.R.-R., R.E.R.-S., P.L., and V.R.-A. analyzed and interpreted the data. M.R.-R., R.E.R.-S., P.L., and V.R.-A. wrote the manuscript. V.R.-A. offered advice in the writing and revising of the manuscript. All of the authors approved the final version of the manuscript.

Funding

The study was supported by the AIDS Research Infrastructure Program (NIMHD RCMI Program; NIMHD-G12-MD007579) and 5 NU58DP005055-05-00 from the Centers for Disease Control and Prevention.

Acknowledgments

The authors wish to thank participant’s willingness to take part in this study and the Puerto Rico Quit line staff for their support in conducting this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ortiz, A.P.; Diaz-Toro, E.C.; Calo, W.A.; Correa-Fernandez, V.; Cases, A.; Santos-Ortiz, M.C.; Mazas, C.; Mejia, L.; Wetter, D.W. Characteristics of smokers accessing the Puerto Rico quitline. Puerto Rico Health Sci. J. 2008, 27, 213–219. [Google Scholar]
  2. Centers for Disease Control and Prevention. State Tobacco Activities Tracking and Evaluation (State) System. Available online: https://www.cdc.gov/statesystem/cigaretteuseadult.html (accessed on 9 September 2018).
  3. Puerto Rico Department of Health. Puerto Rico Chronic Disease Action Plan 2014–2020. Available online: http://www.iccp-portal.org/sites/default/files/plans/Puerto Rico Chronic Disease Action Plan English.pdf (accessed on 7 February 2018).
  4. Stampfli, M.R.; Anderson, G.P. How cigarette smoke skews immune responses to promote infection, lung disease and cancer. Nat. Rev. Immunol. 2009, 9, 377–384. [Google Scholar] [CrossRef] [PubMed]
  5. Warren, G.W.; Alberg, A.J.; Kraft, A.S.; Cummings, K.M. 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–1916. [Google Scholar] [CrossRef] [PubMed]
  6. Hozawa, A.; Houston, T.; Steffes, M.W.; Widome, R.; Williams, O.D.; Iribarren, C.; Pletcher, M.J.; Daviglus, M.L.; Carr, J.J.; Jacobs, D.R., Jr. The association of cigarette smoking with self-reported disease before middle age: The coronary artery risk development in young adults (cardia) study. Prev. Med. 2006, 42, 193–199. [Google Scholar] [CrossRef] [PubMed]
  7. Kumar, P.S.; Matthews, C.R.; Joshi, V.; de Jager, M.; Aspiras, M. Tobacco smoking affects bacterial acquisition and colonization in oral biofilms. Infect. Immun. 2011, 79, 4730–4738. [Google Scholar] [CrossRef] [PubMed]
  8. Yu, G.; Phillips, S.; Gail, M.H.; Goedert, J.J.; Humphrys, M.S.; Ravel, J.; Ren, Y.; Caporaso, N.E. The effect of cigarette smoking on the oral and nasal microbiota. Microbiome 2017, 5, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Pei, H.; Yu, Q.; Xue, Q.; Wei, F.; Tao, L. Previous history and cigarette smoking as interfering factors for the effect of vaccine on human papillomavirus infection. J. Infect. Dis. 2014, 209, 1304. [Google Scholar] [CrossRef] [PubMed]
  10. Wu, J.; Peters, B.A.; Dominianni, C.; Zhang, Y.; Pei, Z.; Yang, L.; Ma, Y.; Purdue, M.P.; Jacobs, E.J.; Gapstur, S.M.; et al. Cigarette smoking and the oral microbiome in a large study of American adults. ISME J. 2016, 10, 2435–2446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Belstrom, D.; Holmstrup, P.; Nielsen, C.H.; Kirkby, N.; Twetman, S.; Heitmann, B.L.; Klepac-Ceraj, V.; Paster, B.J.; Fiehn, N.E. Bacterial profiles of saliva in relation to diet, lifestyle factors, and socioeconomic status. J. Oral. Microbiol. 2014, 6. [Google Scholar] [CrossRef] [PubMed]
  12. Lee, J.; Taneja, V.; Vassallo, R. Cigarette smoking and inflammation: Cellular and molecular mechanisms. J. Dent Res. 2012, 91, 142–149. [Google Scholar] [CrossRef] [PubMed]
  13. Prasad, G.; McCullough, M. Chemokines and cytokines as salivary biomarkers for the early diagnosis of oral cancer. Int. J. Dent. 2013, 2013, 813756. [Google Scholar] [CrossRef] [PubMed]
  14. Garcia-Castillo, V.; Sanhueza, E.; McNerney, E.; Onate, S.A.; Garcia, A. Microbiota dysbiosis: A new piece in the understanding of the carcinogenesis puzzle. J. Med. Microbiol. 2016, 65, 1347–1362. [Google Scholar] [CrossRef] [PubMed]
  15. Hayes, R.B.; Bravo-Otero, E.; Kleinman, D.V.; Brown, L.M.; Fraumeni, J.F., Jr.; Harty, L.C.; Winn, D.M. Tobacco and alcohol use and oral cancer in puerto rico. Cancer Causes Control 1999, 10, 27–33. [Google Scholar] [CrossRef] [PubMed]
  16. Psoter, W.J.; Morse, D.E.; Sanchez-Ayendez, M.; Vega, C.M.; Aguilar, M.L.; Buxo-Martinez, C.J.; Psoter, J.A.; Kerr, A.R.; Lane, C.M.; Scaringi, V.J.; et al. Increasing opportunistic oral cancer screening examinations: Findings from focus groups with general dentists in puerto rico. J. Cancer Educ. 2015, 30, 277–283. [Google Scholar] [CrossRef] [PubMed]
  17. Martell, B.N.; Garrett, B.E.; Caraballo, R.S. Disparities in adult cigarette smoking—United States, 2002–2005 and 2010–2013. MMWR Morb. Mortal Wkly. Rep. 2016, 65, 753–758. [Google Scholar] [CrossRef] [PubMed]
  18. Brown, L.M.; Gridley, G.; Diehl, S.R.; Winn, D.M.; Harty, L.C.; Otero, E.B.; Fraumeni, J.F., Jr.; Hayes, R.B. Family cancer history and susceptibility to oral carcinoma in puerto rico. Cancer 2001, 92, 2102–2108. [Google Scholar] [CrossRef]
  19. Perera, M.; Al-Hebshi, N.N.; Speicher, D.J.; Perera, I.; Johnson, N.W. Emerging role of bacteria in oral carcinogenesis: A review with special reference to perio-pathogenic bacteria. J. Oral. Microbiol. 2016, 8, 32762. [Google Scholar] [CrossRef] [PubMed]
  20. Ferlay, J.; Soerjomataram, I.; Dikshit, R.; Eser, S.; Mathers, C.; Rebelo, M.; Parkin, D.M.; Forman, D.; Bray, F. Cancer incidence and mortality worldwide: Sources, methods and major patterns in globocan 2012. Int. J. Cancer 2015, 136, E359–E386. [Google Scholar] [CrossRef] [PubMed]
  21. Petersen, P.E. Oral cancer prevention and control—The approach of the world health organization. Oral. Oncol. 2009, 45, 454–460. [Google Scholar] [CrossRef] [PubMed]
  22. Hooper, S.J.; Crean, S.J.; Fardy, M.J.; Lewis, M.A.; Spratt, D.A.; Wade, W.G.; Wilson, M.J. A molecular analysis of the bacteria present within oral squamous cell carcinoma. J. Med. Microbiol. 2007, 56, 1651–1659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Guerrero-Preston, R.; Godoy-Vitorino, F.; Jedlicka, A.; Rodriguez-Hilario, A.; Gonzalez, H.; Bondy, J.; Lawson, F.; Folawiyo, O.; Michailidi, C.; Dziedzic, A.; et al. 16s rrna amplicon sequencing identifies microbiota associated with oral cancer, human papilloma virus infection and surgical treatment. Oncotarget 2016, 7, 51320–51334. [Google Scholar] [CrossRef] [PubMed]
  24. Kato, I.; Vasquez, A.A.; Moyerbrailean, G.; Land, S.; Sun, J.; Lin, H.S.; Ram, J.L. Oral microbiome and history of smoking and colorectal cancer. J. Epidemiol. Res. 2016, 2, 92–101. [Google Scholar] [CrossRef] [PubMed]
  25. Zarco, M.F.; Vess, T.J.; Ginsburg, G.S. The oral microbiome in health and disease and the potential impact on personalized dental medicine. Oral. Dis. 2012, 18, 109–120. [Google Scholar] [CrossRef] [PubMed]
  26. Schwabe, R.F.; Jobin, C. The microbiome and cancer. Nat. Rev. Cancer 2013, 13, 800–812. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Ahn, J.; Segers, S.; Hayes, R.B. Periodontal disease, porphyromonas gingivalis serum antibody levels and orodigestive cancer mortality. Carcinogenesis 2012, 33, 1055–1058. [Google Scholar] [CrossRef] [PubMed]
  28. Zackular, J.P.; Rogers, M.A.; Ruffin, M.T.T.; Schloss, P.D. The human gut microbiome as a screening tool for colorectal cancer. Cancer Prev. Res. 2014, 7, 1112–1121. [Google Scholar] [CrossRef] [PubMed]
  29. Michaud, D.S.; Izard, J. Microbiota, oral microbiome, and pancreatic cancer. Cancer J. 2014, 20, 203–206. [Google Scholar] [CrossRef] [PubMed]
  30. Michaud, D.S.; Izard, J.; Wilhelm-Benartzi, C.S.; You, D.H.; Grote, V.A.; Tjonneland, A.; Dahm, C.C.; Overvad, K.; Jenab, M.; Fedirko, V.; et al. Plasma antibodies to oral bacteria and risk of pancreatic cancer in a large european prospective cohort study. Gut 2013, 62, 1764–1770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Farrell, J.J.; Zhang, L.; Zhou, H.; Chia, D.; Elashoff, D.; Akin, D.; Paster, B.J.; Joshipura, K.; Wong, D.T. Variations of oral microbiota are associated with pancreatic diseases including pancreatic cancer. Gut 2012, 61, 582–588. [Google Scholar] [CrossRef] [PubMed]
  32. Littman, D.R.; Pamer, E.G. Role of the commensal microbiota in normal and pathogenic host immune responses. Cell Host Microbe 2011, 10, 311–323. [Google Scholar] [CrossRef] [PubMed]
  33. Kroenke, K.; Spitzer, R.L.; Williams, J.B. The phq-9: Validity of a brief depression severity measure. J. Gen. Intern. Med. 2001, 16, 606–613. [Google Scholar] [CrossRef] [PubMed]
  34. Lowe, B.; Unutzer, J.; Callahan, C.M.; Perkins, A.J.; Kroenke, K. Monitoring depression treatment outcomes with the patient health questionnaire-9. Med. Care 2004, 42, 1194–1201. [Google Scholar] [PubMed]
  35. Lowe, B.; Schenkel, I.; Carney-Doebbeling, C.; Gobel, C. Responsiveness of the phq-9 to psychopharmacological depression treatment. Psychosomatics 2006, 47, 62–67. [Google Scholar] [CrossRef] [PubMed]
  36. Spitzer, R.L.; Williams, J.B.; Kroenke, K.; Hornyak, R.; McMurray, J. Validity and utility of the prime-md patient health questionnaire in assessment of 3000 obstetric-gynecologic patients: The prime-md patient health questionnaire obstetrics-gynecology study. Am. J. Obstet. Gynecol. 2000, 183, 759–769. [Google Scholar] [CrossRef] [PubMed]
  37. Lowe, B.; Kroenke, K.; Herzog, W.; Grafe, K. Measuring depression outcome with a brief self-report instrument: Sensitivity to change of the patient health questionnaire (phq-9). J. Affect. Disord. 2004, 81, 61–66. [Google Scholar] [CrossRef]
  38. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glockner, F.O. Evaluation of general 16s ribosomal rna gene pcr primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41, e1. [Google Scholar] [CrossRef] [PubMed]
  39. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Pena, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. Qiime allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [PubMed]
  40. Navas-Molina, J.A.; Peralta-Sanchez, J.M.; Gonzalez, A.; McMurdie, P.J.; Vazquez-Baeza, Y.; Xu, Z.; Ursell, L.K.; Lauber, C.; Zhou, H.; Song, S.J.; et al. Advancing our understanding of the human microbiome using qiime. Methods Enzymol. 2013, 531, 371–444. [Google Scholar] [PubMed]
  41. Vazquez-Baeza, Y.; Pirrung, M.; Gonzalez, A.; Knight, R. Emperor: A tool for visualizing high-throughput microbial community data. Gigascience 2013, 2, 16. [Google Scholar] [CrossRef] [PubMed]
  42. Lozupone, C.A.; Knight, R. The unifrac significance test is sensitive to tree topology. BMC Bioinform. 2015, 16, 211. [Google Scholar] [CrossRef] [PubMed]
  43. Robertson, C.E.; Harris, J.K.; Wagner, B.D.; Granger, D.; Browne, K.; Tatem, B.; Feazel, L.M.; Park, K.; Pace, N.R.; Frank, D.N. Explicet: Graphical user interface software for metadata-driven management, analysis and visualization of microbiome data. Bioinformatics 2013, 29, 3100–3101. [Google Scholar] [CrossRef] [PubMed]
  44. Kuehl, M.N.; Rodriguez, H.; Burkhardt, B.R.; Alman, A.C. Tumor necrosis factor-alpha, matrix-metalloproteinases 8 and 9 levels in the saliva are associated with increased hemoglobin a1c in type 1 diabetes subjects. PLoS ONE 2015, 10, e0125320. [Google Scholar] [CrossRef] [PubMed]
  45. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [PubMed]
  46. Epskamp, S.; Cramer, A.; Waldorp, L.; Schmittmann, V.; Borsboom, D. Qgraph: Network visualizations of relationships in psychometric data. J. Stat. Softw. 2012, 48, 1–18. [Google Scholar] [CrossRef]
  47. Charlson, E.S.; Chen, J.; Custers-Allen, R.; Bittinger, K.; Li, H.; Sinha, R.; Hwang, J.; Bushman, F.D.; Collman, R.G. Disordered microbial communities in the upper respiratory tract of cigarette smokers. PLoS ONE 2010, 5, e15216. [Google Scholar] [CrossRef] [PubMed]
  48. Grover, N.; Sharma, J.; Sengupta, S.; Singh, S.; Singh, N.; Kaur, H. Long-term effect of tobacco on unstimulated salivary ph. J. Oral. Maxillofac. Pathol. 2016, 20, 16–19. [Google Scholar] [CrossRef] [PubMed]
  49. Han, C.S.; Martin, M.A.; Dichosa, A.E.K.; Daughton, A.R.; Frietze, S.; Kaplan, H.; Gurven, M.D.; Alcock, J. Salivary microbiomes of indigenous tsimane mothers and infants are distinct despite frequent premastication. PeerJ 2016, 4, e2660. [Google Scholar] [CrossRef] [PubMed]
  50. Li, J.; Quinque, D.; Horz, H.P.; Li, M.; Rzhetskaya, M.; Raff, J.A.; Hayes, M.G.; Stoneking, M. Comparative analysis of the human saliva microbiome from different climate zones: Alaska, Germany, and Africa. BMC Microbiol. 2014, 14, 316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Kilian, M.; Chapple, I.L.; Hannig, M.; Marsh, P.D.; Meuric, V.; Pedersen, A.M.; Tonetti, M.S.; Wade, W.G.; Zaura, E. The oral microbiome—An update for oral healthcare professionals. Br. Dent J. 2016, 221, 657–666. [Google Scholar] [CrossRef] [PubMed]
  52. Hoffmann, D.; Wynder, E.L. Chemical constituents and bioactivity of tobacco smoke. IARC Sci. Publ. 1986, 145–165. [Google Scholar]
  53. Rizzatti, G.; Lopetuso, L.R.; Gibiino, G.; Binda, C.; Gasbarrini, A. Proteobacteria: A common factor in human diseases. Biomed. Res. Int. 2017, 2017, 9351507. [Google Scholar] [CrossRef] [PubMed]
  54. Larsen, J.M.; Musavian, H.S.; Butt, T.M.; Ingvorsen, C.; Thysen, A.H.; Brix, S. Chronic obstructive pulmonary disease and asthma-associated proteobacteria, but not commensal prevotella spp., promote toll-like receptor 2-independent lung inflammation and pathology. Immunology 2015, 144, 333–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Patil, S.; Rao, R.S.; Amrutha, N.; Sanketh, D. Oral microbial flora in health. World J. Dent. 2013, 4, 262–266. [Google Scholar] [CrossRef]
  56. Costalonga, M.; Herzberg, M.C. The oral microbiome and the immunobiology of periodontal disease and caries. Immunol. Lett. 2014, 162, 22–38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Tribble, G.D.; Lamont, R.J. Bacterial invasion of epithelial cells and spreading in periodontal tissue. Periodontol. 2000 2010, 52, 68–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Hasan, A.; Palmer, R.M. A clinical guide to periodontology: Pathology of periodontal disease. Br. Dent. J. 2014, 216, 457–461. [Google Scholar] [CrossRef] [PubMed]
  59. Stafford, L.; Berk, M.; Jackson, H.J. Tobacco smoking predicts depression and poorer quality of life in heart disease. BMC Cardiovasc. Disord 2013, 13, 35. [Google Scholar] [CrossRef] [PubMed]
  60. Kim, N.H.; Kim, H.C.; Lee, J.Y.; Lee, J.M.; Suh, I. Association between environmental tobacco smoke and depression among Korean women. BMJ Open 2015, 5, e007131. [Google Scholar] [CrossRef] [PubMed]
  61. Smith, P.H.; Mazure, C.M.; McKee, S.A. Smoking and mental illness in the U.S. Population. Tob. Control 2014, 23, e147–e153. [Google Scholar] [CrossRef] [PubMed]
  62. Malone, K.M.; Waternaux, C.; Haas, G.L.; Cooper, T.B.; Li, S.; Mann, J.J. Cigarette smoking, suicidal behavior, and serotonin function in major psychiatric disorders. Am. J. Psychiatry 2003, 160, 773–779. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, L.; Zhou, C.; Zhu, D.; Wang, X.; Fang, L.; Zhong, J.; Mao, Q.; Sun, L.; Gong, X.; Xia, J.; et al. Serotonin-1a receptor alterations in depression: A meta-analysis of molecular imaging studies. BMC Psychiatry 2016, 16, 319. [Google Scholar] [CrossRef] [PubMed]
  64. Mello, N.K. Hormones, nicotine, and cocaine: Clinical studies. Horm. Behav. 2010, 58, 57–71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Diehl, D.J.; Gershon, S. The role of dopamine in mood disorders. Compr. Psychiatry 1992, 33, 115–120. [Google Scholar] [CrossRef]
  66. Taylor, L.; Loerbroks, A.; Herr, R.M.; Lane, R.D.; Fischer, J.E.; Thayer, J.F. Depression and smoking: Mediating role of vagal tone and inflammation. Ann. Behav. Med. 2011, 42, 334–340. [Google Scholar] [CrossRef] [PubMed]
  67. Vargas, H.O.; Nunes, S.O.; de Castro, M.R.; Vargas, M.M.; Barbosa, D.S.; Bortolasci, C.C.; Venugopal, K.; Dodd, S.; Berk, M. Oxidative stress and inflammatory markers are associated with depression and nicotine dependence. Neurosci. Lett. 2013, 544, 136–140. [Google Scholar] [CrossRef] [PubMed]
  68. Nakamura, Y.; Miyata, M.; Ohba, T.; Ando, T.; Hatsushika, K.; Suenaga, F.; Shimokawa, N.; Ohnuma, Y.; Katoh, R.; Ogawa, H.; et al. Cigarette smoke extract induces thymic stromal lymphopoietin expression, leading to t(h)2-type immune responses and airway inflammation. J. Allergy Clin. Immunol. 2008, 122, 1208–1214. [Google Scholar] [CrossRef] [PubMed]
  69. Singh, S.P.; Gundavarapu, S.; Pena-Philippides, J.C.; Rir-Sima-ah, J.; Mishra, N.C.; Wilder, J.A.; Langley, R.J.; Smith, K.R.; Sopori, M.L. Prenatal secondhand cigarette smoke promotes th2 polarization and impairs goblet cell differentiation and airway mucus formation. J. Immunol. 2011, 187, 4542–4552. [Google Scholar] [CrossRef] [PubMed]
  70. Singh, S.P.; Mishra, N.C.; Rir-Sima-Ah, J.; Campen, M.; Kurup, V.; Razani-Boroujerdi, S.; Sopori, M.L. Maternal exposure to secondhand cigarette smoke primes the lung for induction of phosphodiesterase-4d5 isozyme and exacerbated th2 responses: Rolipram attenuates the airway hyperreactivity and muscarinic receptor expression but not lung inflammation and atopy. J. Immunol. 2009, 183, 2115–2121. [Google Scholar] [PubMed]
  71. Mishra, N.C.; Rir-Sima-Ah, J.; Langley, R.J.; Singh, S.P.; Pena-Philippides, J.C.; Koga, T.; Razani-Boroujerdi, S.; Hutt, J.; Campen, M.; Kim, K.C.; et al. Nicotine primarily suppresses lung th2 but not goblet cell and muscle cell responses to allergens. J. Immunol. 2008, 180, 7655–7663. [Google Scholar] [CrossRef] [PubMed]
  72. Kali, A. Human microbiome engineering: The future and beyond. J. Clin. Diagn. Res. 2015, 9, DE01–DE04. [Google Scholar] [CrossRef] [PubMed]
  73. Arthur, J.C.; Perez-Chanona, E.; Muhlbauer, M.; Tomkovich, S.; Uronis, J.M.; Fan, T.J.; Campbell, B.J.; Abujamel, T.; Dogan, B.; Rogers, A.B.; et al. Intestinal inflammation targets cancer-inducing activity of the microbiota. Science 2012, 338, 120–123. [Google Scholar] [CrossRef] [PubMed]
  74. Solovyeva, I.A.; Sobko, E.A.; Ishchenko, O.P.; Kraposhina, A.Y.; Demko, I.V.; Eidemiller, N.S. Impact of tobacco smoking on the clinical and functional indicators and markers of systemic inflammation in patients with severe asthma. Ter. Arkh. 2015, 87, 27–33. [Google Scholar] [CrossRef] [PubMed]
  75. Holownia, A.; Wielgat, P.; Rysiak, E.; Braszko, J.J. Intracellular and extracellular cytokines in a549 cells and thp1 cells exposed to cigarette smoke. Adv. Exp. Med. Biol. 2016, 910, 39–45. [Google Scholar] [PubMed]
  76. Zhao, J.; Li, X.; Xie, F.; Yang, Z.; Pan, X.; Zhu, M.; Shang, P.; Nie, C.; Liu, H.; Xie, J. Immunomodulatory effects of cigarette smoke condensate in mouse macrophage cell line. Int. J. Immunopathol. Pharmacol. 2017, 30, 315–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Dwivedi, S.; Goel, A.; Khattri, S.; Mandhani, A.; Sharma, P.; Pant, K.K. Tobacco exposure by various modes may alter proinflammatory (il-12) and anti-inflammatory (il-10) levels and affects the survival of prostate carcinoma patients: An explorative study in north Indian population. Biomed. Res. Int. 2014, 2014, 158530. [Google Scholar] [CrossRef] [PubMed]
  78. Ritter, M.; Goggel, R.; Chaudhary, N.; Wiedenmann, A.; Jung, B.; Weith, A.; Seither, P. Elevated expression of tarc (ccl17) and mdc (ccl22) in models of cigarette smoke-induced pulmonary inflammation. Biochem. Biophys. Res. Commun. 2005, 334, 254–262. [Google Scholar] [CrossRef] [PubMed]
  79. Zier, K.S.; Leo, M.M.; Spielman, R.S.; Baker, L. Decreased synthesis of interleukin-2 (il-2) in insulin-dependent diabetes mellitus. Diabetes 1984, 33, 552–555. [Google Scholar] [CrossRef] [PubMed]
  80. Musselman, D.; Royster, E.B.; Wang, M.; Long, Q.; Trimble, L.M.; Mann, T.K.; Graciaa, D.S.; McNutt, M.D.; Auyeung, N.S.; Oliver, L.; et al. The impact of escitalopram on il-2-induced neuroendocrine, immune, and behavioral changes in patients with malignant melanoma: Preliminary findings. Neuropsychopharmacology 2013, 38, 1921–1928. [Google Scholar] [CrossRef] [PubMed]
  81. Chen, Y.; Ouyang, J.; Liu, S.; Zhang, S.; Chen, P.; Jiang, T. The role of cytokines in the peripheral blood of major depressive patients. Clin. Lab. 2017, 63, 1207–1212. [Google Scholar] [CrossRef] [PubMed]
  82. Molofsky, A.B.; Nussbaum, J.C.; Liang, H.E.; Van Dyken, S.J.; Cheng, L.E.; Mohapatra, A.; Chawla, A.; Locksley, R.M. Innate lymphoid type 2 cells sustain visceral adipose tissue eosinophils and alternatively activated macrophages. J. Exp. Med. 2013, 210, 535–549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Gutierrez-Meza, J.M.; Jarillo-Luna, R.A.; Rivera-Aguilar, V.; Miliar-Garcia, A.; Campos-Rodriguez, R. Cytokine profile of nalt during acute stress and its possible effect on iga secretion. Immunol. Lett. 2017, 188, 68–78. [Google Scholar] [CrossRef] [PubMed]
  84. Nakata, M.; Yamamoto, S.; Okada, T.; Yada, T. Aav-mediated il-10 gene transfer counteracts inflammation in the hypothalamic arcuate nucleus and obesity induced by high-fat diet. Neuropeptides 2017, 62, 87–92. [Google Scholar] [CrossRef] [PubMed]
  85. Nakata, M.; Yamamoto, S.; Okada, T.; Gantulga, D.; Okano, H.; Ozawa, K.; Yada, T. Il-10 gene transfer upregulates arcuate pomc and ameliorates hyperphagia, obesity and diabetes by substituting for leptin. Int. J. Obes. 2016, 40, 425–433. [Google Scholar] [CrossRef] [PubMed]
  86. Wei, R.; Hu, Y.; Dong, F.; Xu, X.; Hu, A.; Gao, G. Hepatoma cell-derived leptin downregulates the immunosuppressive function of regulatory t-cells to enhance the anti-tumor activity of cd8+ t-cells. Immunol. Cell Biol. 2016, 94, 388–399. [Google Scholar] [CrossRef] [PubMed]
  87. Filozof, C.; Fernandez Pinilla, M.C.; Fernandez-Cruz, A. Smoking cessation and weight gain. Obes. Rev. 2004, 5, 95–103. [Google Scholar] [CrossRef] [PubMed]
  88. Rimm, E.B.; Manson, J.E.; Stampfer, M.J.; Colditz, G.A.; Willett, W.C.; Rosner, B.; Hennekens, C.H.; Speizer, F.E. Cigarette smoking and the risk of diabetes in women. Am. J. Public Health 1993, 83, 211–214. [Google Scholar] [CrossRef] [PubMed]
  89. Perry, I.J.; Wannamethee, S.G.; Walker, M.K.; Thomson, A.G.; Whincup, P.H.; Shaper, A.G. Prospective study of risk factors for development of non-insulin dependent diabetes in middle aged British men. BMJ 1995, 310, 560–564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Kawakami, N.; Takatsuka, N.; Shimizu, H.; Ishibashi, H. Effects of smoking on the incidence of non-insulin-dependent diabetes mellitus. Replication and extension in a Japanese cohort of male employees. Am. J. Epidemiol. 1997, 145, 103–109. [Google Scholar] [CrossRef] [PubMed]
  91. Persson, P.G.; Carlsson, S.; Svanstrom, L.; Ostenson, C.G.; Efendic, S.; Grill, V. Cigarette smoking, oral moist snuff use and glucose intolerance. J. Intern. Med. 2000, 248, 103–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Oda, E.; Kawai, R. Significance of heart rate in the prevalence of metabolic syndrome and its related risk factors in japanese. Circ. J. 2009, 73, 1431–1436. [Google Scholar] [CrossRef] [PubMed]
  93. Cnop, M.; Landchild, M.J.; Vidal, J.; Havel, P.J.; Knowles, N.G.; Carr, D.R.; Wang, F.; Hull, R.L.; Boyko, E.J.; Retzlaff, B.M.; et al. The concurrent accumulation of intra-abdominal and subcutaneous fat explains the association between insulin resistance and plasma leptin concentrations: Distinct metabolic effects of two fat compartments. Diabetes 2002, 51, 1005–1015. [Google Scholar] [CrossRef] [PubMed]
  94. Wang, T.N.; Chang, W.T.; Chiu, Y.W.; Lee, C.Y.; Lin, K.D.; Cheng, Y.Y.; Su, Y.J.; Chung, H.F.; Huang, M.C. Relationships between changes in leptin and insulin resistance levels in obese individuals following weight loss. Kaohsiung J. Med. Sci. 2013, 29, 436–443. [Google Scholar] [CrossRef] [PubMed]
  95. Purwar, P.; Khan, M.A.; Mahdi, A.A.; Pandey, S.; Singh, B.; Dixit, J.; Sareen, S. Salivary and serum leptin concentrations in patients with chronic periodontitis. J. Periodontol. 2015, 86, 588–594. [Google Scholar] [CrossRef] [PubMed]
  96. Eliasson, B.; Smith, U. Leptin levels in smokers and long-term users of nicotine gum. Eur. J. Clin. Investig. 1999, 29, 145–152. [Google Scholar] [CrossRef]
  97. McKee, S.A.; Sinha, R.; Weinberger, A.H.; Sofuoglu, M.; Harrison, E.L.; Lavery, M.; Wanzer, J. Stress decreases the ability to resist smoking and potentiates smoking intensity and reward. J. Psychopharmacol. 2011, 25, 490–502. [Google Scholar] [CrossRef] [PubMed]
  98. Von der Goltz, C.; Koopmann, A.; Dinter, C.; Richter, A.; Rockenbach, C.; Grosshans, M.; Nakovics, H.; Wiedemann, K.; Mann, K.; Winterer, G.; et al. Orexin and leptin are associated with nicotine craving: A link between smoking, appetite and reward. Psychoneuroendocrinology 2010, 35, 570–577. [Google Scholar] [CrossRef] [PubMed]
  99. McDonald, S.D.; Walker, M.; Perkins, S.L.; Beyene, J.; Murphy, K.; Gibb, W.; Ohlsson, A. The effect of tobacco exposure on the fetal hypothalamic-pituitary-adrenal axis. BJOG 2006, 113, 1289–1295. [Google Scholar] [CrossRef] [PubMed]
  100. Garcia, C.M.C.; Garcia, R.J.F.; Mangas, R.A.; Millan, J. Tobacco and arterial pressure (ii.). The acute effects on the angiotensin-converting enzyme. An. Med. Intern. 1990, 7, 392–395. [Google Scholar]
  101. Christensen, H.; Korczak, B.M.; Bojesen, A.M.; Kuhnert, P.; Frederiksen, W.; Bisgaard, M. Classification of organisms previously reported as the sp and stewart-letscher groups, with descriptions of necropsobacter gen. Nov. And of necropsobacter rosorum sp. Nov. For organisms of the sp group. Int. J. Syst. Evol. Microbiol. 2011, 61, 1829–1836. [Google Scholar] [CrossRef] [PubMed]
  102. Myhrvold, V.; Brondz, I.; Olsen, I. Application of multivariate analyses of enzymic data to classification of members of the actinobacillus-haemophilus-pasteurella group. Int. J. Syst. Bacteriol. 1992, 42, 12–18. [Google Scholar] [CrossRef] [PubMed]
  103. Tan, T.Y.; Ng, L.S.; Eng, L.C. Five clinical cases of necropsobacter rosorum bacteremia. J. Clin. Microbiol. 2013, 51, 1337–1340. [Google Scholar] [CrossRef] [PubMed]
  104. Nonnenmacher, C.; Mutters, R.; de Jacoby, L.F. Microbiological characteristics of subgingival microbiota in adult periodontitis, localized juvenile periodontitis and rapidly progressive periodontitis subjects. Clin. Microbiol. Infect. 2001, 7, 213–217. [Google Scholar] [CrossRef] [PubMed]
  105. Lehnerdt, G.; Senska, K.; Fischer, M.; Jahnke, K. smoking promotes the formation of peritonsillar abscesses. Laryngorhinootologie 2005, 84, 676–679. [Google Scholar] [CrossRef] [PubMed]
  106. Schafer, P.; Furrer, C.; Mermillod, B. An association of cigarette smoking with recurrent subareolar breast abscess. Int. J. Epidemiol. 1988, 17, 810–813. [Google Scholar] [CrossRef] [PubMed]
Ijerph 15 02479 g001 550
Figure 1. Clustering of samples from nonsmokers and smokers. (A) for beta diversity, principal coordinate analysis plots (PCoA) based on unweighted UniFrac distance were used to visualize the variation present in all samples by smoking status. Emperor web browser was used to perform the analysis. Significant cluster was observed in the saliva microbial community among smokers (blue dots) and nonsmokers (red dots). (B) unweighted UniFrac plot representing the intraindividual variability within nonsmokers (NS) and smokers (S). The asterisks (*) indicate a significant difference (nonparametric p-value with Bonferroni correction), where p = 0.01.
Figure 1. Clustering of samples from nonsmokers and smokers. (A) for beta diversity, principal coordinate analysis plots (PCoA) based on unweighted UniFrac distance were used to visualize the variation present in all samples by smoking status. Emperor web browser was used to perform the analysis. Significant cluster was observed in the saliva microbial community among smokers (blue dots) and nonsmokers (red dots). (B) unweighted UniFrac plot representing the intraindividual variability within nonsmokers (NS) and smokers (S). The asterisks (*) indicate a significant difference (nonparametric p-value with Bonferroni correction), where p = 0.01.
Ijerph 15 02479 g001
Ijerph 15 02479 g002 550
Figure 2. Taxonomic composition of nonsmokers and smokers at the phylum level. FASTA files were used to determine the taxonomy composition of each sample with a 97% similarity threshold using the Greengenes reference dataset.
Figure 2. Taxonomic composition of nonsmokers and smokers at the phylum level. FASTA files were used to determine the taxonomy composition of each sample with a 97% similarity threshold using the Greengenes reference dataset.
Ijerph 15 02479 g002
Ijerph 15 02479 g003 550
Figure 3. Composition of the major genera observed in nonsmoker (NS; n = 16) and smoker (S; n = 18) saliva samples. (A) comparison of the taxonomic analyses at the genus level by group. (B) composition of the major genera by study participant. The bacterial abundance was established by the analysis of the 16S rRNA gene using a Quantitative Insights into Microbial Ecology (QIIME) pipeline.
Figure 3. Composition of the major genera observed in nonsmoker (NS; n = 16) and smoker (S; n = 18) saliva samples. (A) comparison of the taxonomic analyses at the genus level by group. (B) composition of the major genera by study participant. The bacterial abundance was established by the analysis of the 16S rRNA gene using a Quantitative Insights into Microbial Ecology (QIIME) pipeline.
Ijerph 15 02479 g003
Ijerph 15 02479 g004 550
Figure 4. Bacteria associations and correlation network integrating hormones, cytokines, and bacteria variables enriched in the smoker group. (A) LEfSe-associated bacteria with nonsmoker and smoker groups, with a p-value < 0.05 and linear discriminant analysis (LDA) >3.0. (B) correlation network of bacteria associated by LEfSE in the smoker group, and statistically significant hormones and cytokines. Only significant correlations are shown as determined by qgraph. Gray lines represent negative correlations, whereas black lines represent positive correlations. The line thickness determines the strength of the correlation.
Figure 4. Bacteria associations and correlation network integrating hormones, cytokines, and bacteria variables enriched in the smoker group. (A) LEfSe-associated bacteria with nonsmoker and smoker groups, with a p-value < 0.05 and linear discriminant analysis (LDA) >3.0. (B) correlation network of bacteria associated by LEfSE in the smoker group, and statistically significant hormones and cytokines. Only significant correlations are shown as determined by qgraph. Gray lines represent negative correlations, whereas black lines represent positive correlations. The line thickness determines the strength of the correlation.
Ijerph 15 02479 g004
Table
Table 1. Summary of study subjects by smoking status.
Table 1. Summary of study subjects by smoking status.
VariableNonsmokersSmokersp-Value *
Individuals (n)1618-
Age in years (mean, range)34 (28–50)45 (33–49)0.109
Sex (male/female)6/1010/80.384
RaceHispanicHispanic-
Education level (percentage) 0.695
  ≤High school4 (25%)6 (33%)
  ≥College12 (75%)12 (67%)
Total years smoking-≤3 years = 2 >3 years = 15 NR a = 1-
Cigarettes per day (number)
  Light/moderate (≤20 cigarettes)-9-
  Heavy (≥20 cigarettes)-9-
Alcohol use (yes/no)11/57/110.721
Assessment of depression symptoms mean ± standard deviation6.56 ± 7.3610.89 ± 7.040.020
a NR: not reported; * exact p-values were based on Mann-Whitney U tests, with the significance level set at 0.05.
Table
Table 2. Analytes in saliva samples from smokers and nonsmokers.
Table 2. Analytes in saliva samples from smokers and nonsmokers.
AnalyteNonsmokers Mean ± Standard Deviation (pg/mL)Smokers Mean ± Standard Deviation Mean (pg/mL)p-Value *
TNF-α28.3 ± 46.1911.2 ± 11.330.646
IL-12 (p70)8.4 ± 5.05182.1 ± 749.430.281
MDC ↓311.8 ± 399.44114.8 ± 2460.007
IL-10 ↓7.7 ± 12.201.6 ± 0.234.28 × 10−4
IFN-γ5.2 ± 3.506.9 ± 7.190.746
TNF-β2.5 ± 1.121.8 ± 0.330.126
IL-1 β443.1 ± 1358.49214.4 ± 352.950.384
IL-5 ↓2.7 ± 1.051.7 ± 0.120.001
IL-2 ↑1.1 ± 0.022.4 ± 1.219.07 × 10−10
IL-637.1 ± 92.9426.1 ± 71.880.102
IL-4 ↑1.6 ± 1.4815.6 ± 17.101.6E-4
IL-1ra10,021 ± 2117.6510,296 ± 4958.390.746
IL-131.9 ± 1.131.3 ± 0.900.365
IL-176.4 ± 10.958.5 ± 12.240.164
IL-711.3 ± 8.837.4 ± 9.900.050
GM-CSF30.9 ± 53.9122.5 ± 16.640.721
ACTH ↑1.5 ± 0.809.2 ± 7.445.83 × 10−4
Insulin ↓364.7 ± 87.4027.3 ± 22.233.66 × 10−4
Leptin ↓33.8 ± 16.0319.3 ± 4.276.77 × 10−4
Analytes were measured using the Luminex MAGPIX (with xPONENT 4.1 program) and a cytokine magnetic bead panel (MILLIPLEX HCYTMAG-60K-PX39); (↑): significantly higher levels of analytes in the saliva samples from smokers; (↓): significantly lower levels of analytes in the saliva samples from smokers; * p-values were based on the Mann-Whitney-Wilcoxon test, with a significance level set at 0.05.

Share and Cite

MDPI and ACS Style

Rodríguez-Rabassa, M.; López, P.; Rodríguez-Santiago, R.E.; Cases, A.; Felici, M.; Sánchez, R.; Yamamura, Y.; Rivera-Amill, V. Cigarette Smoking Modulation of Saliva Microbial Composition and Cytokine Levels. Int. J. Environ. Res. Public Health 2018, 15, 2479. https://doi.org/10.3390/ijerph15112479

AMA Style

Rodríguez-Rabassa M, López P, Rodríguez-Santiago RE, Cases A, Felici M, Sánchez R, Yamamura Y, Rivera-Amill V. Cigarette Smoking Modulation of Saliva Microbial Composition and Cytokine Levels. International Journal of Environmental Research and Public Health. 2018; 15(11):2479. https://doi.org/10.3390/ijerph15112479

Chicago/Turabian Style

Rodríguez-Rabassa, Mary, Pablo López, Ronald E. Rodríguez-Santiago, Antonio Cases, Marcos Felici, Raphael Sánchez, Yasuhiro Yamamura, and Vanessa Rivera-Amill. 2018. "Cigarette Smoking Modulation of Saliva Microbial Composition and Cytokine Levels" International Journal of Environmental Research and Public Health 15, no. 11: 2479. https://doi.org/10.3390/ijerph15112479

APA Style

Rodríguez-Rabassa, M., López, P., Rodríguez-Santiago, R. E., Cases, A., Felici, M., Sánchez, R., Yamamura, Y., & Rivera-Amill, V. (2018). Cigarette Smoking Modulation of Saliva Microbial Composition and Cytokine Levels. International Journal of Environmental Research and Public Health, 15(11), 2479. https://doi.org/10.3390/ijerph15112479

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

Article Metrics

Back to TopTop