INTRODUCTION

Diarrhea is a common gastrointestinal disorder clinically defined by loose stools and more frequent bowel movements1. Despite significant global efforts in improving hygiene conditions and public health awareness, diarrhea remains one of the significant contributors to the global disease burden2. Previous research3 has shown that chronic diarrhea or loose or watery stools affected up to 26.9% of adults in the United States, severely impacting people’s quality of life. Furthermore, according to data from 2016, there were over 4.4 billion cases of diarrhea worldwide, resulting in 1.6 million deaths2, generating substantial medical and healthcare costs, and having a profound impact on socioeconomy4. Additionally, chronic diarrhea is positively associated with the risk of colon cancer, cardiovascular diseases, and all-cause mortality5.

Smoking, a modifiable risk factor, has long been recognized as closely linked to respiratory and cardiovascular diseases6,7. However, its impact on the gastrointestinal tract remains a complex area of research. Previous studies have indicated that smoking is a risk factor for Crohn’s disease (CD) and irritable bowel syndrome (IBS)8,9, while other studies have noted conflicting results10,11. Although the aforementioned conditions may present with diarrhea symptoms12, research directly investigating the association between tobacco exposure and diarrhea remains scarce to date. Meanwhile, the relationship between exposure to SHS and bowel health also remains controversial13-15. Therefore, further research is required to clarify the relationship between smoking and diarrhea, rather than merely within the framework of disease specificity. Given this background, we utilized the available data from NHANES to assess the relationship between different tobacco exposure situations and diarrhea.

METHODS

Study population and data sources

This is a secondary dataset analysis of data from the National Health and Nutrition Examination Survey (NHANES), which represents a continuous cross-sectional surveillance program employing a stratified, multistage probability-clustering sampling strategy to obtain nationally representative data from non-institutionalized United States civilians16. This surveillance system, administered by the National Center for Health Statistics (NCHS), operates under Protocol 2005–2006 with formal authorization from the NCHS Institutional Review Board. Written informed consent was obtained from all participants following ethical review and approval by the National Health Statistics Ethics Review Committee (Protocol No. 2005-06).

Our study participants were derived from NHANES conducted between 2005 and 2010, for which we have pooled the data of three waves of investigation: 2005–2006, 2007–2008, 2009–2010, for a cross-sectional approach. Participants were included in the study if they met the following criteria: aged ≥20 years, had complete tobacco exposure information, completed the Gut Health Questionnaire, and had complete covariate information.

Bowel health questionnaire

The Bowel health questionnaire was completed in the Mobile Examination Center (MEC) survey room using a computerized personal survey system. The BSFS, commonly used in research and clinical practice, was used to define chronic diarrhea. Participants were shown a card with color pictures and descriptions of the seven BSFS (types 1 to 7)17. Participants who indicated their typical or most frequent stool type as BSFS Type 6 or BSFS Type 7 were classified as having diarrhea. Other participants were classified as having normal bowel habits18.

Exposure: Smoking

The tobacco exposure situations were evaluated by trained interviewers using the Computer-Assisted Personal Interviewing (CAPI) system. Based on the classification of smoking status, this study further categorizes participants’ exposure to SHS. Specifically, based on responses to the questions ‘Have you smoked at least 100 cigarettes in your lifetime?’ and ‘Do you currently smoke cigarettes?’, participants who answered ‘yes’ to both questions are considered ‘smokers’, those who answered ‘yes’ to the first question but ‘no’ to the second are considered ‘ex-smokers’, and those who answered ‘no’ to both questions are considered ‘never smokers’. For never smokers, based on their response to the question ‘Does anyone smoke inside your home?’, those who answered ‘yes’ are considered ‘never smokers exposed to SHS’, and those who answered ‘no’ are considered ‘never smokers not exposed to SHS’.

For smokers, the smoking duration is calculated as current age minus age at smoking initiation (age at smoking initiation was based on their answers to ‘Age started smoking cigarettes regularly?’). The average number of cigarettes consumed daily is calculated by their answers to ‘Avg cigarettes/day during past 30 days’. Based on the participants’ responses, they were categorized into three groups: 1–10, 11–20, and >20 cigarettes per day. Additionally, participants were classified as those smoking immediately after waking (within 5 minutes) or non-immediate smokers based on their answers to ‘How soon after waking do you smoke?’.

Measurement of serum cotinine levels

Isotope dilution-high-performance liquid chromatography/atmospheric pressure chemical ionization tandem mass spectrometry was used to measure serum cotinine levels19. As suggested by a previous study20, we created cotinine categories using the newly recommended cutoff point of 3 ng/mL to represent smoking exposure19. Those below the lower detection threshold were considered unexposed.

Covariates

Based on previous studies21,22 and our clinical experience, we preliminarily selected a set of potential factors as covariates to analyze the relationship between tobacco exposure and diarrhea, including gender, age, race (Mexican American, Other Hispanic, Non-Hispanic White, Non-Hispanic Black, and Other), education level (lower than high school, high school, higher than high school), marital status (married, unmarried, other), poverty-to-income ratio (PIR: <1.30, 1.30–3.49, ≥3.50)21, BMI (kg/m2) (normal weight <25, overweight 25–29.9, obese ≥30)21 , hypertension (hypertension, non-hypertension)22 and diabetes(diabetes, non-diabetes, borderline)22. The BMI data were measured in the MEC by trained health technicians. Other covariates were collected by trained interviewers using the CAPI system.

Statistical analysis

All analyses accounted for sample weights. In descriptive analyses, categorical variables were presented through weighted percentages, whereas continuous variables were represented by means and standard deviations (SD). To detect variations in baseline characteristics between diarrhea and non-diarrhea patients for categorical data, the chi-squared test was used. For continuous variables, the t-test was applied to evaluate differences between groups.

Weighted logistic regression models were used to investigate the relationship between different tobacco exposure scenarios and diarrhea. First, with never smokers not exposed to SHS as the reference group, we evaluated the likelihood of having diarrhea among three groups: smokers, ex-smokers, and never smokers exposed to SHS. Additionally, for smokers (not including ex-smokers, never smokers not exposed to SHS, and never smokers exposed to SHS), using three separate regression models, we separately examined the effects of smoking time (immediate smokers and non-immediate smokers), smoking duration, and smoking dose on the likelihood of having diarrhea. For all the logistic regressions, we used three models. Model 1 did not account for any potential covariates. Model 2 was adjusted for age, sex, and race. Model 3 was additionally adjusted for marital status, education level, household income, BMI, hypertension, and diabetes. For smokers (not included ex-smokers, never smokers not exposed to SHS, never smokers exposed to SHS), restricted cubic spline (RCS) analysis with four knots (the 5th, 35th, 65th, and 95th quantiles), controlling for all covariates, was conducted to examine potential nonlinear relationships between smoking duration and the likelihood of diarrhea, and likelihood ratio tests were employed to test for nonlinear relationship. Stratification analyses were used to assess the robustness of the tobacco exposure-diarrhea association and to assess whether the association was modified by specific variables. The variables examined included gender, education level (lower than high school, high school, higher than high school), PIR (<1.30, 1.30–3.49, ≥3.50), BMI (normal weight, overweight, obese), hypertension (hypertension, and non-hypertension), and diabetes (diabetes, non-diabetes, and borderline). Except for the stratification component itself, all covariates were adjusted: age, gender, race, education level, marital status, PIR, BMI, hypertension, and diabetes.

The Baron and Kenny four-step method was employed to examine the mediating role of serum cotinine in the association between tobacco exposure and diarrhea: initially, the association between tobacco exposure (i.e. smokers) and diarrhea was confirmed based on the aforementioned analysis; subsequently, logistic regression analyses were used to investigate the associations between tobacco exposure and serum cotinine concentration; next, the association between serum cotinine concentration and chronic diarrhea was assessed and finally, serum cotinine concentration was included as a covariate to examine whether it could attenuate the association between tobacco exposure and diarrhea. Furthermore, for all multivariable models that may contain multicollinearity in this research, variance inflation factor (VIF) values were calculated to assess multicollinearity, with VIF <10 as the acceptable criterion.

All analyses were performed using R version 3.4.3, Empower software (X&Y Solutions, Inc., Boston, MA), and Stata software version 18 (StataCorp, College Station, TX, USA). A two-sided p<0.05 was considered statistically significant.

RESULTS

Basic characteristics of participants

Participant selection followed a sequential exclusion process: Firstly, individuals without data on tobacco exposure were excluded (n=112). Then, participants lacking covariate information (n=2174) and those lacking data on bowel health (n=1528) were removed. Finally, 13318 participants were included in this study. The specific screening process for participants is illustrated in Figure 1.

Figure 1

The screening process flowchart for the cross-sectional analysis of tobacco exposure and diarrhea in the National Health and Nutrition Examination Survey (NHANES) 2005–2010 (N=13318)

https://www.tobaccoinduceddiseases.org/f/fulltexts/224313/TID-24-144-g001_min.jpg

The average age of the 13318 participants was 46.5 ± 16.5 years, of which 1012 (7.6%) were diagnosed with chronic diarrhea. Table 1 demonstrates that chronic diarrhea exhibited higher prevalence among older adults (p<0.001), females (p<0.001), those with lower level of education (p<0.001), individuals with reduced PIR (p<0.001), obese participants (p<0.001), hypertensive patients (p<0.001), diabetic patients (p<0.001), and tobacco-exposed individuals (p<0.001).

Table 1

Baseline characteristics of participants in a cross-sectional study of tobacco exposure and diarrhea in the NHANES 2005–2010 (N=13318)

CharacteristicsAll n (%)Non-diarrhea n (%)Diarrhea n (%)p
Total, n13318123061012
Age (years), mean ± SD46.5 ± 16.546.2 ± 16.549.9 ± 15.7<0.001
Gender<0.001
Male6490 (48.7)6052 (49.2)428 (42.3)
Female6828 (51.3)6254 (50.8)584 (57.7)
Race0.063
Mexican American1037 (7.8)943 (7.7)98 (9.6)
Other Hispanic544 (4.1)496 (4.0)48 (4.7)
Non-Hispanic White9594 (72.0)8899 (72.3)689 (68.1)
Non-Hispanic Black1433 (10.8)1312 (10.7)123 (12.2)
Other710 (5.3)656 (5.3)54 (5.3)
Education level<0.001
Lower than high school2356 (17.7)2105 (17.1)264 (26.0)
High school3225 (24.2)2969 (24.1)257 (25.4)
Higher than high school7737 (58.1)7232 (58.8)491 (48.5)
Marriage0.002
Married/cohabitating8723 (65.5)8061 (65.5)661 (65.3)
Unmarried2186 (16.4)2049 (16.6)133 (13.1)
Other2409 (18.1)2196 (17.8)218 (21.6)
PIR<0.001
<1.302533 (19.0)2288 (18.6)254 (25.1)
1.30–3.494787 (35.9)4429 (36.0)357 (35.3)
≥3.505999 (45.0)5590 (45.4)401 (39.6)
BMI<0.001
Normal weight4160 (31.2)3893 (31.6)259 (25.6)
Overweight4472 (33.6)4158 (33.8)309 (30.5)
Obese4686 (35.2)4255 (34.6)444 (43.9)
Hypertension<0.001
Hypertension4025 (30.2)3647 (29.6)390 (38.5)
Non-hypertension9293 (69.8)8659 (70.4)622 (61.5)
Diabetes<0.001
Diabetes1048 (7.9)919 (7.5)137 (13.5)
Non-diabetes12041 (90.4)11185 (90.9)846 (83.6)
Borderline229 (1.7)202 (1.6)29 (2.9)
Tobacco exposure<0.001
Never smokers not exposed to SHS6658 (50.0)6218 (50.5)429 (42.4)
Never smokers exposed to SHS383 (2.9)349 (2.8)35 (3.4)
Ex-smokers3296 (24.7)3027 (24.6)272 (26.9)
Smokers2981 (22.4)2712 (22.0)276 (27.3)

[i] PIR: poverty-to-income ratio. SHS: secondhand smoke. BMI: body mass index (kg/m2); normal weight <25.0, overweight 25.0–29.9, obese ≥30.0. As this is a descriptive statistic, all variables in the table are independently included in a separate model for calculation. All analyses were conducted using sample weights to ensure national representativeness. The p-value of continuous variables is obtained through linear regression, while the p-value of categorical variables is obtained through chi-squared test.

The relationship between tobacco exposure and diarrhea

Table 2 shows the logistic regression analyses of different tobacco exposure situations and diarrhea. Compared to never smokers not exposed to SHS, smokers exhibited a high probability of diarrhea. After adjusting for all covariates, the association remained significant (AOR=1.56; 95% CI: 1.25–1.93). However, no significant relationships were observed in either never smokers exposed to SHS (AOR=1.37; 95% CI: 0.87–2.16) or ex-smokers (AOR=1.22; 95% CI: 0.98–1.51) when compared to the never smoker not exposed to SHS.

Table 2

Logistic multivariable regression analysis results of the association between different tobacco exposure situations and diarrhea in the NHANES 2005–2010 (N=13318; smokers=2952)

CharacteristicsModel 1Model 2Model 3
OR (95% CI)pAOR (95% CI)pAOR (95% CI)p
Tobacco exposure situation
Never smokers not exposed to SHS (ref.)111
Never smokers exposed to SHS1.44 (0.93–2.24)0.1051.54 (0.98–2.40)0.0581.37 (0.87–2.16)0.168
Ex-smokers1.30 (1.06–1.60)0.0111.26 (1.01–1.56)0.0381.22 (0.98–1.51)0.079
Smokers1.48 (1.21–1.81)<0.0011.67 (1.15–2.06)<0.0011.56 (1.25–1.93)<0.001
Smokers only
Smoking time
Immediate smokers (ref.)111
Non-immediate smokers1.62 (1.13–2.30)0.0081.63 (1.14–2.34)0.0081.69 (1.17–2.45)0.005
Smoking duration (years)
Continuous (every additional year)1.02 (1.01–1.03)0.0011.02 (1.01–1.03)0.0011.02 (1.00–1.03)0.007
Smoking dose (cigarettes per day)
Continuous (every additional cigarette per day)1.01 (1.00–1.02)0.2141.01 (1.00–1.02)0.2071.01 (0.99–1.02)0.201
≤10 (ref.)111
11–200.97 (0.67–1.38)0.8470.96 (0.65–1.41)0.8420.97 (0.66–1.43)0.890
>201.25 (0.81–1.95)0.3151.29 (0.79–2.10)0.3141.30 (0.79–2.13)0.296

[i] Model 1: no covariates were adjusted. Model 2: adjusted for age, gender, and race. Model 3: adjusted as for Model 2 plus education level, marital status, PIR, BMI, hypertension and diabetes. AOR: adjusted odds ratio. SHS: secondhand smoke. The p-value is determined by the Wald test.

Among smokers (not including ex-smokers, never smokers not exposed to SHS, never smokers exposed to SHS), those who smoked immediately upon waking exhibited a significantly elevated likelihood of diarrhea after adjusting for all covariates (AOR=1.69; 95% CI: 1.17–2.45) (Table 2). Each additional year of smoking duration was associated with a 2% increase in diarrhea likelihood (AOR=1.02; 95% CI: 1.00–1.03). RCS curve (Figure 2) further confirmed a linear relationship between smoking duration and diarrhea (p for overall=0.025; p for nonlinear=0.587). However, whether smoking dose was treated as a continuous variable (AOR=1.01; 95% CI: 0.99–1.02) or a categorical variable (11–20 vs ≤10 cigarettes per day: AOR=0.97; 95% CI: 0.66–1.43; >20 vs ≤10 cigarettes per day: AOR=1.30; 95% CI: 0.79–2.13), no significant association was observed between smoking dose and diarrhea probability.

Figure 2

Restricted cubic spline (RCS) analysis of the relationship between smoking duration and diarrhea (only for smokers) in the NHANES 2005–2010 (N=2952)

https://www.tobaccoinduceddiseases.org/f/fulltexts/224313/TID-24-144-g002_min.jpg

Notably, as the analyses investigating the associations of smoking time, smoking duration, and smoking dose with diarrhea were restricted to the smokers, the sample size for the latter three regression models in Table 2 and Figure 2 corresponds to the number of smokers, specifically n=2952. Besides, given the presence of multicollinearity between age and smoking duration, age was not adjusted for as a covariate in Model 2 and Model 3 when examining the relationship between smoking duration and diarrhea (Table 2). And the VIF values of the above two models are both <2.5.

Stratification analysis (Figure 3) confirmed the stability of this association, with more pronounced tobacco exposure-diarrhea relationships observed in males, individuals with higher level of education, those with elevated PIR values, normal weight individuals, non-hypertensive and non-diabetic populations. However, due to the limited sample size, none of the never smokers exposed to SHS with borderline diabetes had diarrhea, so there is no specific value in the last row of column A. Furthermore, it should be clarified that stratified analysis is classified as exploratory analysis without multiple comparison correction and does not provide confirmatory conclusions. Its findings need to be validated in prospective studies with larger sample sizes and thus cannot be used as a basis for stratified intervention in clinical practice for the time being.

Figure 3

Stratification analysis results on the association between different tobacco exposure situations and diarrhea in the NHANES 2005–2010 (N=13318)

https://www.tobaccoinduceddiseases.org/f/fulltexts/224313/TID-24-144-g003_min.jpg

Mediation analysis

In the above study, we have found that compared to never smokers not exposed to SHS, only smokers showed a significant increase in the likelihood of diarrhea, indicating that smokers passed the first step of the Baron and Kenny four-step method for mediation analysis. Therefore, we conducted the follow-up analysis on smokers.

As demonstrated in Table 3, after adjusting for all covariates, smokers exhibited higher odds of exceeding the serum cotinine threshold compared to never smokers not exposed to SHS. However, smokers with elevated serum cotinine levels did not show a statistically significant increase in diarrhea incidence (AOR=1.19; 95% CI: 0.85–1.68).

Table 3

Mediating effect analysis of serum cotinine concentration on the association between tobacco exposure and diarrhea in NHANES 2005–2010 (N=3273; smokers=985)

CharacteristicsModel 1Model 2Model 3
OR (95% CI)pAOR (95% CI)pAOR (95% CI)p
Tobacco exposure
Never smokers not exposed to SHS (ref.)111
Smokers526.56 (341.68–811.47)<0.001719.21 (444.82–1162.86)<0.001655.02 (396.55–1081.95)<0.001
Cotinine
Under threshold (ref.)111
Above threshold1.12 (0.83–1.52)0.4671.34 (0.97–1.84)0.0741.19 (0.85–1.68)0.306

[i] Model 1 no covariates were adjusted. Model 2 adjusted for age, gender, and race. Model 3 adjusted as for Model 2 plus education level, marital status, PIR, BMI, hypertension and diabetes. AOR: adjusted odds ratio. SHS: secondhand smoke. The table presents results of logistic regression analyses conducted separately for the association between tobacco exposure and serum cotinine levels, as well as the relationship between serum cotinine levels and diarrhea. The two regression analyses are all a part of the mediation analysis. The p-value is determined by the Wald test.

Specifically, Table 3 shows the second and third steps of Baron and Kenny’s four-step method for mediation analysis, respectively. As the research results of the third step were not statistically significant, the mediation pathway was not established, and therefore the fourth step analysis was not conducted. Furthermore, it should be clarified that among the 13318 participants, serum cotinine data were available for only 4518 individuals, comprising 2288 never smokers not exposed to SHS and 985 smokers. Consequently, in the first regression analysis investigating the association between tobacco exposure and serum cotinine concentration presented in this table, the sample size was 3273 (n=3273). For the second regression analysis examining the relationship between serum cotinine concentration and diarrhea restricted to smokers, the sample size corresponded to 985 smokers (n=985).

DISCUSSION

Utilizing pooled data from NHANES 2005–2010, this study represents a systematic investigation of the association between tobacco exposure and diarrhea. We identified that active tobacco exposure was associated with an increased likelihood of diarrhea, while passive exposure and former smoking showed no significant association. Stratification analysis confirmed the robustness of these findings. Both longer smoking duration and smoking immediately after waking were associated with elevated diarrhea likelihood, whereas daily cigarette consumption demonstrated no significant association. Additionally, no mediating effect of serum cotinine was observed in the relationship between tobacco exposure and diarrhea.

As a modifiable lifestyle factor, smoking is one of the earliest identified environmental risk factors for IBD23. Previous study24 has shown that smokers and ex-smokers exhibit a significantly higher prevalence of CD compared to non-smokers. However, the relationship between smoking and UC remains controversial. A study by Mahid et al.11 suggested that former smokers had an increased risk of UC compared to never smokers, while current smokers had a reduced risk. However, a recent Mendelian randomization study9 using genetic data from three independent populations identified smoking as a causal risk factor for UC rather than a protective factor. For IBS, a meta-analysis by Nicholas et al.25 established smoking as a risk factor for diarrhea-predominant IBS. Although diarrhea is not equivalent to these specific diseases, the above studies have to some extent confirmed a significant association between active tobacco exposure and diarrhea, which is also the main conclusion of our study.

Research on the association between passive tobacco exposure and diarrhea remains limited and inconsistent. A Japanese study13 linked passive tobacco exposure to increased UC risk, while van der Heide et al.14 identified passive smoking as a risk factor for CD. Conversely, an Israeli multicenter study15 found no association between passive smoking and IBD. Additionally, research on passive tobacco exposure and IBS is lacking. In our study, never smokers exposed to SHS showed no statistically significant difference in diarrhea likelihood compared to never smokers not exposed to SHS, suggesting that the association between passive tobacco exposure and diarrhea is not significant. However, due to the limited sample size of never smokers exposed to SHS in this study and the uncertain dose of SHS exposure, we cannot currently draw conclusions about the association between passive tobacco exposure and diarrhea based on this study. Larger scale prospective studies are needed in the future to determine this relationship.

Tobacco smoke is a highly complex aerosol containing over 4000 compounds26, and current assessments of tobacco exposure primarily rely on serum cotinine levels, a metabolite of nicotine27. Nicotine has been demonstrated to induce multiple intestinal effects, including hypoperfusion of the rectum and of acutely damaged colonic tissue28 and reduced smooth muscle tone/contractility29. Concurrently, nicotine may modulate the interaction between cytokine profile dysregulation and cell cycle response alterations in UC and CD, thereby influencing disease progression30. However, our study revealed no association between serum cotinine concentrations and diarrhea, suggesting other tobacco constituents instead of nicotine may drive this association. Importantly, given IBS substantially higher prevalence compared to UC and CD in diarrheal populations31,32, this observation does not negate nicotine’s pathogenic roles in UC and CD pathophysiology, but rather implies nicotine-independent mechanisms may underlie IBS and other diarrhea-predominant intestinal disorders.

Interestingly, our study found that individuals who smoked immediately after waking had a significantly higher probability of developing diarrhea compared to other smokers. Although the mechanism remains unclear, this observation may relate to the fasting state. Reduced intestinal content during fasting may allow tobacco components to directly interact with the gut microbiota and mucosal cells, potentially leading to dysbiosis or oxidative damage33. Additionally, fasting-state exposure might enhance the direct stimulation of the vagus nerve and intestinal smooth muscle by tobacco components, triggering contractions and functional diarrhea or IBS34,35. Besides, the association between smoking duration and increased diarrhea likelihood may reflect cumulative inflammatory effects over time; this finding also provides novel insights into subsequent mechanistic exploration.

Strengths and limitations

This study investigated the relationship between different tobacco exposure statuses and diarrhea using a nationally representative adult sample. Additionally, we examined the mediating effect of serum cotinine concentration and discussed the potential mechanisms of action, which provide a theoretical basis for subsequent mechanistic studies and large-scale prospective studies. Furthermore, we adjusted for many potential covariate factors using a broad set of covariates, thereby enhancing the reliability of the study results. However, this study has the following limitations: First, this study is a cross-sectional study; therefore, we cannot infer causality between tobacco exposure and diarrhea, and the results are limited to discussions of association. Second, this study is based on a single NHANES database, so the findings may not be generalizable to other countries, and the extrapolation of the results is limited. Third, although many covariate factors were adjusted for in this study, residual confounding may still exist, which could influence the results. Fourth, a large amount of data in this study was collected based on respondents’ recall, so recall bias and social desirability bias are unavoidable. In addition, this study defined diarrhea using the Bristol Stool Form Scale; although this definition has been recognized in previous studies, according to the World Health Organization, the definition of diarrhea should also consider stool frequency, so the results of this study may have certain biases.

CONCLUSIONS

Active tobacco exposure demonstrates a significant association with elevated diarrhea likelihood, whereas no substantial associations were observed for passive exposure or ex-smokers, a conclusion reinforced by stratification analysis. Both extended smoking duration and immediate post-awakening smoking were associated with increased diarrhea likelihood, though daily cigarette consumption showed no significant association. Furthermore, serum cotinine levels exhibited no mediating effect in this relationship.