Spicy Food Intake and Dietary Factors Shape the Gut Microbiome and Metabolism of Mucin and Short-Chain Fatty Acids in Healthy Adults
Uigi Min 1
Jihyun Kim 1
Jinwook Kim 1
Hyunmi Jin 1
Hyunji Oh 1
Soyeon Ahn 1
Hyeonseok Shin 1
Wonhee Lee 1✉ Email
1
A
Research Center, Mitomics Institute Samyang Foods Inc Seoul Republic of Korea
Uigi Min1,+, Jihyun Kim1,+, Jinwook Kim1, Hyunmi Jin1, Hyunji Oh1, Soyeon Ahn1,
Hyeonseok Shin1, and Wonhee Lee1,*
1Research Center, Mitomics Institute, Samyang Foods Inc., Seoul, Republic of Korea
*Correspondence: whlee@roundsquare.ai
Uigi Min and Jihyun Kim contributed equally to this work.
Abstract
The intestinal mucus layer provides a primary barrier between the gut microbiota and the epithelium. Capsaicin, a major component of spicy foods, has been reported to enhance muc2 expression and mucus secretion; however, evidence from preclinical and clinical studies remains inconsistent. How spicy food interacts with external factors such as alcohol to influence mucus-layer homeostasis also remains insufficiently understood. Using shotgun metagenomics, we investigated the associations between spicy food intake, alcohol consumption, intestinal fatty acid–binding protein (I-FABP), liver fatty acid–binding protein (L-FABP), gut microbial composition, and functional pathways in 229 healthy Korean adults. Alcohol intake was positively correlated with urinary I-FABP (ρ = 0.26), indicating mild epithelial stress, whereas spicy food intake was not associated with either FABP biomarker. Individuals who consumed highly spicy food exhibited an increased abundance of short-chain fatty acids (SCFA)-producing and mucin-metabolizing taxa, including Blautia, Coprococcus, and Ruminococcus, along with higher activity of mucin-degradation and SCFA-production pathways. Individuals who consumed a high level of alcohol showed stronger enrichment of mucin-degrading taxa, with reduced SCFA flux and increased abundance of Proteobacteria and Fusobacteriota. Cross-classified dietary groups revealed distinct mucin and SCFA gene activity patterns. Notably, the DHSH group displayed concurrent elevation of mucin turnover and SCFA production with indications of dysbiosis. These findings suggest that spicy food may modulate mucus layer metabolism in a context-dependent manner, whereas alcohol more consistently perturbs mucin–SCFA networks and epithelial integrity.
Keywords:
Gut microbiome
Mucin metabolism
Spicy food intake
Alcohol consumption
Shotgun metagenomics
Short-chain fatty acids
Intestinal epithelial integrity
Intestinal-fatty acid–binding protein
Liver-fatty acid–binding protein
A
A
A
Introduction
The mucus layer, constituting the primary interface between the gut microbiota and host epithelium, performs multiple physiological functions beyond acting as a simple physical barrier, including providing a habitat for microorganisms, regulating the exchange of metabolic products, and coordinating immune responses1,2. This layer is a key determinant of the stability of the gut ecosystem, and structural or chemical perturbations of the mucus layer have been reported as early indicators of increased intestinal permeability, chronic inflammation, and metabolic disorders2. Therefore, understanding the factors that regulate mucus layer homeostasis is essential for characterizing the health of the gut environment.
A
Dietary factors are recognized as major external variables that modify the composition and function of the mucus layer. Capsaicin, the principal bioactive compound in spicy foods, has emerged as a modulator capable of directly influencing mucus layer ecology. Experimental studies evaluating the physiological effects of spicy foods have shown that capsaicin increases muc2 expression in intestinal epithelial cells and promotes the expansion of mucin-utilizing microbial taxa, thereby supporting mucus-layer regeneration and stability3,4. Capsaicin has also been reported to enhance short-chain fatty acid (SCFA) biosynthesis pathways and regulate host energy metabolism5, suggesting a role as a metabolic mediator that strengthens interactions between the host and gut microbiota. However, the findings across various studies are inconsistent68. In rodents given high doses of capsaicin (40 mg/kg), mucosal irritation and barrier disruption were evident9, whereas human intervention studies have shown discomfort at low doses (0.75 mg/day) and metabolic improvements at higher doses (135 mg/day)10,11, indicating that the physiological response to capsaicin is not a linear dose-dependent relationship.
Similarly, epidemiological studies have reported associations between spicy food intake and reduced all-cause mortality and metabolic disease risk1214, although these associations remain inconsistent across populations. This inconsistency is likely attributable to the co-consumption of spicy foods with other dietary factors, such as high-fat meals or alcohol, which may modify or obscure the intrinsic effects of capsaicin15. Thus, the influence of spicy food consumption on the mucus layer and microbial function requires evaluation within the broader context of real-world dietary patterns, emphasizing the need for analytical approaches that can quantitatively capture such contextual variation.
Most prior studies of the gut microbiota have been based on 16S rRNA sequencing, which provides insights into microbial composition and diversity but offers limited resolution for interpreting functional regulation at the metabolic pathway level, particularly for mucus-related gene clusters or SCFA biosynthetic pathways that are directly linked to gut barrier homeostasis16. Consequently, these approaches are insufficient for evaluating the ecological significance of dietary exposure. In contrast, shotgun metagenomics enables the direct profiling of microbial genes and functional characterization of metabolic networks modulated by dietary factors, including spicy food intake17,18.
In this study, we performed shotgun metagenomic profiling of a cohort of more than 200 healthy Korean adults to investigate the mucus-associated metabolic pathways and functional features of the gut microbiota in relation to the frequency of consumption of spicy food. The goal was to characterize how food intake influences the microbial networks involved in mucus utilization, SCFA production, and intestinal barrier maintenance. By interpreting the relationship between spicy food consumption and the gut ecosystem at both structural and functional levels, this study provides a basis for a more refined understanding of the physiological mechanisms through which dietary habits contribute to gut barrier homeostasis.
RESULTS
Sample Collection and Baseline Characteristics
A
A total of 240 healthy adults were initially recruited for this study. However, 11 participants withdrew consent or were lost to follow-up, resulting in a final cohort of 229 individuals.
A
All participants completed a structured questionnaire on dietary and lifestyle habits and provided fecal, urine, and oral epithelial samples. Urine samples were used to quantify inflammation- and epithelial injury-related biomarkers, including TNF-α, IL-6, I-FABP, and L-FABP.
A
Oral samples were used to measure TRPV1 expression, and fecal samples were subjected to shotgun metagenomic sequencing (Fig. 1A). The demographic and lifestyle characteristics of the participants are summarized in Supplementary Table 1.
Fig. 1
Study design and metadata correlation. A) Overview of the study design. B) Correlation plot showing Spearman’s correlation coefficients between dietary factors, meta info, and intestinal epithelial injury biomarkers (I-FABP, L-FABP).
Click here to Correct
The study group consisted of 87 male (38%) and 142 female participants (62%). The mean age was 38.3 ± 11.1 years, and the mean body mass index (BMI) was 22.9 ± 3.4 kg/m², indicating that most individuals were within the normal weight range. The average frequency of spicy food consumption was 1.9 ± 0.7 points, corresponding to approximately “2–4 times per month.” Consumption of high-fat foods (1.9 ± 0.8), instant foods (1.7 ± 0.8), and fast foods (1.4 ± 0.7) fell within the “rarely to < 1 time per week” range, indicating generally healthy dietary habits. Vegetable preference relative to meat was moderate (4.2 ± 2.0 on a 1–10 scale), suggesting a tendency toward vegetable-based diets. Vitamin (1.3 ± 0.5) and probiotic intake scores (1.6 ± 0.5) indicated infrequent use. Alcohol drinking frequency (1.7 ± 1.0) and drinking volume per occasion (1.6 ± 1.0) corresponded to “2–4 times per month, 1–3 drinks per occasion.” Smoking scores were low (0.2 ± 0.5), reflecting a predominantly non-smoking population.
Correlations Between Biomarkers and Metadata
Spicy food is known to stimulate gastrointestinal physiological responses and mucus secretion, producing both beneficial and adverse health effects68, whereas alcohol consumption has been reported as a representative disruptor of the intestinal barrier, inducing epithelial injury and suppressing mucus secretion19,20. Considering these prior findings, we designed the present study with the aim of evaluating how these two dietary factors are associated with markers of epithelial injury and the composition of the gut microbiota.
We first assessed, through analysis of Spearman’s correlation, the relationships between survey-based demographic, dietary, and lifestyle variables and a number of biomarkers: I-FABP and L-FABP, reflecting intestinal barrier damage; inflammatory markers (TNF-α and IL-6); and the sensory receptor for pungency (TRPV1). In the biomarker assessment, TNF-α, IL-6, and TRPV1 were detected at low concentrations and were therefore excluded from subsequent analyses (Supplementary Fig. 1).
Analysis of the Spearman’s correlation between demographic variables and epithelial injury markers showed that I-FABP exhibited weak positive correlations with age (ρ = 0.22), BMI (ρ = 0.19), and stool frequency (ρ = 0.24), and a negative correlation with sex (ρ = − 0.27) (Fig. 1B). L-FABP showed a positive correlation with I-FABP (ρ = 0.31), suggesting that the two markers reflect similar epithelial injury responses. Examination of correlations between dietary factors and biomarkers revealed that alcohol intake (g/day) was positively correlated with I-FABP (ρ = 0.26), whereas no correlation was observed with L-FABP (ρ = 0.05). Spicy food intake showed little association with I-FABP (ρ = 0.01) or L-FABP (ρ = 0.05), and the same was true for other dietary variables.
These patterns suggest that alcohol consumption represents a dietary factor exerting a more direct and sensitive influence on intestinal epithelial injury. Accordingly, the interpretation of gut microbial composition and functional pathway changes requires analytical strategies that incorporate alcohol intake as a principal independent variable while simultaneously adjusting for sex, age, BMI, and stool frequency as confounders. This approach enables a clearer separation of the effects of spicy food versus alcohol consumption and provides a more accurate framework for explaining the physiological variations reflected by epithelial injury markers.
Besides, alcohol exhibited some correlation (although weak) with high-fat and instant food intake, indicating that it may be part of a broader cluster of co-occurring dietary behaviors. Alcohol intake also showed a weak positive correlation with consumption of instant foods (ρ = 0.22) and preference for high-fat foods (ρ = 0.20). Lifestyle variables such as smoking, physical activity, and sleep quality were not correlated with either I-FABP or L-FABP.
Microbiome Composition and Diversity
Gut microbiome profiles were generated by shotgun metagenomic sequencing of fecal samples collected from 229 participants (Supplementary Fig. 2A). A total of 1,527 species were detected across all samples, and the two phyla commonly reported to be dominant in the human gut microbiome, Bacillota and Bacteroidota, were the most abundant21. When evaluating the extent to which environmental and dietary factors contributed to the taxonomic variation, sex, stool frequency, and age showed the highest explanatory power, followed by BMI, I-FABP levels, smoking, and alcohol consumption (Fig. 2A). In contrast, other factors, including meat–vegetable preference, vitamin intake, and sleep quality, explained less than 1% of the variation in microbial community composition.
Fig. 2
Gut microbiota composition in relation to spicy food and alcohol intake. A) PERMANOVA assessment of variation in microbial taxonomy explained by different environmental or dietary variables. The x-axis indicates the percentage of variance explained (R2), while point size reflects the significance level as -log10(p). Variables with P < 0.05 are highlighted in deeppink, and non-significant variables are shown in gray. B) Boxplots showing the expression levels of epithelial injury markers (I-FABP, L-FABP) according to spicy food and alcohol groups. P-values were calculated using the Wilcoxon rank sum test. C–D) Differentially abundant microbial taxa associated with spicy food (C) and alcohol intake (D), identified by MaAsLin2. Bubble plots show taxa with significant associations (FDR < 0.25) between the low- and high-intake groups in each dietary category. The x-axis represents the MaAsLin2 coefficient (effect size), where positive values indicate higher relative abundance in the SH or DH groups, and negative values indicate higher abundance in the SL or DL groups. Bubble size corresponds to the absolute value of the coefficient, and colors indicate the direction of association (red: enriched in the high-intake group; blue: enriched in the low-intake group). Taxa are ordered by effect size.
Click here to Correct
Group Classification by Spicy Food Intake and Comparison of Gut Microbiome Profiles
A
To determine the gut microbial composition as it relates to spicy food consumption, participants were classified as high-intake (SH; n = 40) or low-intake (SL; n = 189) based on consumption frequency. The SH group showed significantly higher spicy food intake scores than did the SL group (p < 0.001), whereas fast-food consumption was significantly lower (p = 0.003). There were no differences between the groups in terms of demographic or lifestyle variables, including age, BMI, stool frequency, or alcohol consumption. No significant differences were observed for other health-related factors, such as use of vitamins, probiotics, or antibiotics (Supplementary Table 1). Shotgun metagenome–based diversity analyses showed no distinct separation in alpha diversity, and Bray–Curtis distance–based principal coordinates analysis (PCoA) also indicated no clear separation of taxonomic composition between the groups (PERMANOVA, p = 0.09) (Supplementary Fig. 2C).
Species-level differential abundance was assessed using MaAsLin2 with a linear mixed-effects model, which identified 21 species with significant group differences (Fig. 2C). The SH group showed significantly higher abundances of the following microbial taxa: Dialister invisus, Blautia_A caecimuris, Limisoma sp000437795, Ruminococcus_E sp003526955, Dorea_A longicatena_B, Faecalimonas phoceensis, Blautia hansenii, Anaerotignum sp001304995, Blautia sp001304935, Amedibacillus dolichus, Eubacterium_R sp000436835, Blautia_A sp900549015, Phocaeicola sp900542985, Butyribacter sp001916135, Butyricicoccus pullicaecorum, CAG-83 sp000435975, and Fusobacterium_A mortiferum. Among these taxa, species of the genera Blautia, Limisoma, Ruminococcus, Butyribacter, and Butyricicoccus are associated with SCFA production and mucus layer metabolism4,22–25, suggesting that spicy food intake may contribute positively to gut metabolic homeostasis. However, species of Fusobacterium have been linked to intestinal inflammation26 and genera such as Dialister and Dorea have been shown to have context-dependent associations with health in metagenomic studies27–31. In contrast, the SL group showed significantly higher abundances of Dialister hominis, Lawsonibacter asaccharolyticus, UMGS856 sp900546265, and Gordonibacter pamelaeae. Lawsonibacter is a known butyrate producer32, and Gordonibacter is associated with ellagic acid metabolism33.
Group Classification by Alcohol Consumption and Comparison of Gut Microbiome Profiles
Next, participants were classified according to their daily alcohol intake (g/day) into a high-consumption group (DH; n = 109) and a low-consumption group (DL; n = 120). The DH group showed a markedly higher alcohol intake than did the DL group (p < 0.001), and the proportion of smokers was also significantly higher (p < 0.001). Among the dietary variables, consumption of spicy food (p = 0.034), high-fat food (p = 0.014), and instant food (p = 0.003) was higher in the DH group. No significant differences were observed between the groups in terms of age, BMI, stool frequency, Bristol stool score, abdominal pain, sleep quality, or other demographic and lifestyle variables (Supplementary Table 1). Shotgun metagenome–based diversity analyses showed no distinct separation in either alpha or beta diversity between the groups (PERMANOVA, p = 0.241) (Supplementary Fig. 2C).
MaAsLin2 analysis of species-level differences identified 18 microbial taxa that were significantly associated with alcohol consumption (Fig. 2D). In the DH group, Faecalimonas umbilicata, Amedibacillus dolichus, Lawsonibacter sp900545895, Evetepia excrementipullorum, Weissella koreensis, Butyricicoccus pullicaecorum, Fimisoma sp000435715, and Frisingicoccus caecimuris were significantly enriched. Species belonging to Lawsonibacter and Butyricicoccus have been reported to be butyrate producers32,34, and species of Weissella are lactic acid bacteria originating from kimchi35. The presence of certain beneficial taxa in the DH group may reflect ecological selection under alcohol-induced stress conditions, wherein species capable of surviving or adapting to such environments persist or show an increase in their relative abundance36,37. In the DL group, the following strains were relatively abundant: CAG-83 sp900548615, CAG-95 sp900066375, CAG-103 sp900543625, CAG-170 sp900553545, Blautia_A sp900066355, Alistipes_A ihumii, Ventricola sp900542395, Limivivens sp900543575, Lawsonibacter sp014287875, and Gemmiger sp900539695. Notably, species of Alistipes and Blautia possess SCFA-producing and mucin-degrading capacities that may support the maintenance of mucosal homeostasis38. The distinction in I-FABP levels between the spicy food- and alcohol-defined groups further supports the notion that dietary factors may differentially influence intestinal epithelial injury (Fig. 2B). I-FABP and L-FABP levels did not differ significantly between the SH and SL groups; in contrast, I-FABP levels were significantly higher in the DH group than in the DL group (p = 0.004).
Comparison of Epithelial Injury Markers and Gut Microbiome Profiles Across Combined Spicy Food–Alcohol Intake Groups
As spicy food consumption and alcohol intake are often co-occurring behaviors, as noted in previous studies, it is necessary to evaluate their combined effects by examining interaction groups15. Based on the dietary questionnaire responses regarding the frequency of spicy food consumption and daily alcohol intake (g/day), the participants were classified into four groups: DLSL (Drink-Low-Spicy-Low; n = 101), DHSL (Drink-High-Spicy-Low; n = 88), DLSH (Drink-Low-Spicy-High; n = 19), and DHSH (Drink-High-Spicy-High; n = 21). Among the demographic and lifestyle characteristics, the sex distribution differed significantly across the groups (p = 0.001), whereas age, BMI, stool frequency, and Bristol Stool Scale scores did not show group-level differences (Supplementary Table 3). Other lifestyle variables, including smoking, physical activity, sleep quality, and hygiene awareness, did not differ significantly between the groups. Comparison of epithelial injury biomarkers showed that I-FABP concentrations were significantly higher in the DHSL group than in the DLSL group (p = 0.0271) (Fig. 3C). Analysis of shotgun metagenomic data revealed no distinct differences in gut microbiome diversity among the four groups. Neither the Shannon index nor the PERMANOVA-based beta-diversity analysis indicated significant differences in taxonomic composition (PERMANOVA, p = 0.218) (Supplementary Fig. 2B-C).
Fig. 3
Group-specific differences in gut microbial species composition. A) Heatmap showing the abundances of 53 microbial species that are significantly different among groups. The top annotated color bars represent metadata variables, with corresponding color scales shown on the right. Species abundances were log10-transformed and z-scored across all samples. B) Phylogenetic distribution of 53 microbial species showing significant group-wise differences identified by MaAsLin2 analysis (FDR < 0.25). Colors represent the regression coefficients from the MaAsLin2 model, with red indicating positive and blue indicating negative associations relative to the other groups. C) Boxplots showing the expression levels of epithelial injury markers (I-FABP, L-FABP) according to different groups. P-values were calculated by the Kruskal–Wallis rank sum test.
Click here to Correct
Next, we examined the differences in species-level abundance associated with the four groups using MaAsLin2. A total of 53 microbial species were significantly associated with at least one group (Fig. 3A). The DHSH group displayed an overall enrichment in Bacillota taxa, most notably multiple Blautia species, including Blautia_A sp900549015, Blautia hansenii, Blautia sp900541955, and Blautia sp001304935. Additional enriched taxa included Roseburia sp003470905, Butyricicoccus pullicaecorum, Agathobaculum sp900544475, Dialister invisus, and Faecalimonas phoceensis. In the DLSH group, a diverse set of Bacillota species was also enriched, particularly members of the Ruminococcaceae and Lachnospiraceae families, such as RuminiclostridiumE siraeum, Anaerobutyricum soehngenii, Butyribacter intestini, Coprococcus eutactus_A, Coprococcus_A catus, Blautia stercoris, Dorea_A longicatena_B, and Lachnoclostridium_B sp900066555. Notably, Coprococcus species are characteristic of this group. Limisoma sp000437795 and Turicibacter sp001543345 were also distinctly enriched, representing only the Muribaculaceae and Turicibacteraceae taxa observed among the groups. The DHSL group was characterized by the enrichment of Frisingicoccus caecimuris, Clostridium_Q sp000435655, and Angelakisella massiliensis, and, unlike other groups, did not show an increased abundance of Blautia species. The DLSL group, representing a low intake of both alcohol and spicy food, showed higher levels of typical commensals, such as Faecalibacterium prausnitzii, Blautia_A luti, Dysosmobacter sp900550685, and Lawsonibacter asaccharolyticus. The 53 associated species spanned multiple families within Bacillota, and to a lesser extent, Bacteroidota, demonstrating phylogenetically diverse responses to spicy food and alcohol intake. The genera Blautia, Coprococcus, and Limisoma were consistently more abundant in the DHSH and DLSH groups (Fig. 3B).
Global Metabolic Pathways Based on Functional Gene Profiles
A
To explore the global metabolic potential of the gut microbiome, unstratified KEGG and MetaCyc pathway abundances were computed using HUMAnN 3.0, followed by MaAsLin2 analysis (Fig. 4). In total, 38 major pathways exhibited differential distributions among the four groups. Based on the MetaCyc ontology, these pathways were categorized into amino acids, carbohydrates, lipids, cofactors, nucleotides, aromatic compounds, secondary metabolites, and energy metabolism.
Fig. 4
Significantly enriched functional metabolic pathways of the gut microbiota. The bar plot shows 38 MetaCyc pathways that differ significantly among the groups. The x-axis represents -log10(P-value) derived from the MaAsLin2 model, and the color of each bar indicates the group with the corresponding enrichment.
Click here to Correct
The DHSH group was characterized by carbohydrate metabolic pathways, such as hexitol fermentation to lactate, formate, ethanol, and acetate (P461-PWY) and acetyl-CoA fermentation to butanoate II (PWY-5676), along with amino acid metabolic pathways, including L-isoleucine biosynthesis III (PWY-5103) and the super-pathway of L-serine and glycine biosynthesis I (SER-GLYSYN-PWY). The formaldehyde assimilation III (P185-PWY) pathway was also enriched in this group.
The DLSH group showed a prominent enrichment of amino acid metabolism pathways, such as L-methionine biosynthesis I (HOMOSER-METSYN-PWY) and L-alanine fermentation to propanoate and acetate (PROPFERM-PWY), but also exhibited a broader diversity of pathways related to lipid, energy, and cofactor metabolism than was evident in the other groups. Notably, aromatic compound degradation pathways, including catechol degradation I (meta-cleavage pathway) (PWY-5415) and 4-hydroxyphenylacetate degradation (3-HYDROXYPHENYLACETATE-DEGRADATION-PWY) were detected.
The DHSL group displayed relatively few distinct features, with L-proline biosynthesis II (from arginine) (PWY-4981) emerging as the dominant amino acid pathway. The DLSL group showed higher abundance of cofactor-related pathways, such as coenzyme A biosynthesis I (COA-PWY) and super-pathway of thiamine diphosphate biosynthesis I (THISYN-PWY), as well as nucleotide metabolism pathways, including inosine-5'-phosphate biosynthesis I (PWY-6123).
Comparison of SCFA-Targeted Gene Family Distributions among Groups
A
To assess the activation status of short-chain fatty acid (SCFA)-producing pathways, species-level stratified KEGG ortholog (KO) profiles were compared between the groups (Fig. 5A). SCFA-related genes were analyzed based on a previously established list of SCFA-producing-associated KOs39. Group-wise comparison revealed that the DHSH group exhibited the highest overall abundance of gene families associated with acetate, propionate, and butyrate production (K00925, K01738, K16363, K01646, K19697, K01034, K01443, and K02535). The major contributing species were Anaerostipes hadrus, Blautia hansenii, Escherichia coli, and Fusobacterium mortiferum. Similarly, the DLSH group showed elevated abundance of acetate-, propionate-, and butyrate-producing gene families (K01026, K00929, and K00138, respectively), with Blautia obeum, Blautia producta, and Klebsiella pneumoniae identified as key contributors. In contrast, the DLSL group displayed a high abundance of the acetate-producing KO K16363 associated with Bacteroides nordii, whereas the propionate- and butyrate-producing pathways were comparatively reduced. The DHSL group exhibited only a limited presence of SCFA production–related genes, suggesting relatively weak activation of these metabolic pathways.
Fig. 5
Comparison of SCFA and mucin-related gene families of the gut microbiota. A) Heatmap illustrating the associations among short-chain fatty acid (SCFA)-related KEGG Orthologs (KOs), microbial species, and compounds (especially, acetate, butyrate, and propionate). Significantly associated KO-species obtained from HUMAnN3 were identified using MaAsLin2 analysis (FDR < 0.25). Lower panel: group-specific KO abundances; upper panel: SCFA compounds associated with each KO and the corresponding microbial species expressing them. B) Heatmap of mucin degradation-related KO and glycoside hydrolases (GHs), and their associated microbial species.
Click here to Correct
Comparison of Mucin-Degrading Gene Family Distributions among Groups
To investigate the mucin degradation–related gene families, both KEGG ortholog (KO) and Glycoside Hydrolase (GH) family annotations were analyzed based on the glycan utilization gene set proposed by Labarthe et al. (2023)40. Distinct combinations of mucin-derived glycan degradation and downstream metabolic pathways were observed in relation to the level of consumption of spicy food and alcohol (Fig. 4B). The DHSH group was enriched in gene families related to starch degradation (GH13), cellulose/β-glucan degradation (GH5), sulfate reduction (K00394, K00395, K00956), and galactose metabolism (K00965, K00849). Moreover, characteristic enrichment of acetyl-CoA-derived SCFA-producing flux genes (K00625, K00925, K00656, K01034, and K00626) was observed, along with glycolytic enzymes (K01689, K01810, and K00927) and C1 metabolism–associated enzymes (K01491 and K00297). The major associated species were Anaerostipes hadrus and Fusobacterium mortiferum. In the DLSH group, numerous glycoside hydrolase families were enriched, including those involved in xylan degradation (GH10, GH120, GH30, and GH39), cellulose degradation (GH16 and GH26), and broad glycan hydrolysis (GH3), indicating a strong capacity for complex polysaccharide degradation. Additional enrichment in the galactose-containing glycan degradation (GH35) and β-galactosidase/β-mannosidase families (GH2) was observed. Metabolic pathways associated with fructose, mannose, and fucose (K01628, K01813, and K01818), sulfate assimilation (K00957), and butyrate production (K00929) were also identified. The dominant contributing species was Blautia producta, with additional contributions from Coprococcus eutactus and Fusobacterium mortiferum. The DHSL group exhibited characteristic enrichment of glycolytic pathways (K00134 and K00927), acetate formation (K00925), and sulfate assimilation (K00957). Within glycoside hydrolases, α-fucosidase families (GH29, GH95), broad glycosidase activity (GH3), and glucuronoxylan xylanohydrolase (GH30) were predominant. Furthermore, genes associated with the uronate pathway (xylulokinase, K00854) and C1 metabolism (methylenetetrahydrofolate dehydrogenase, K01491) were detected and were mainly attributed to Parabacteroides distasonis. In contrast, the DLSL group showed enrichment of glycolytic enzymes such as enolase (K01689) and pyruvate kinase (K01810) related to Bacteroides nordii, but lacked mucin- or diet-derived glycan–degrading enzymes and downstream SCFA-producing genes.
Network Analysis
Integrated network analysis revealed that host-related factors (alcohol, spicy food, I-FABP, L-FABP, and TRPV1) exhibited multiple correlations with key microbial species and metabolic gene clusters (Fig. 6). Consumption of spicy food was positively correlated with abundance of Blautia producta and Fusobacterium mortiferum, which in turn was associated with genes related to SCFA synthesis (K00925, K00929, K01026, and K01034). The abundance of these species also showed positive correlations with mucin-degrading genes (K00965 and K00626), clustering with Blautia obeum, Coprococcus eutactus, Anaerostipes hadrus, Klebsiella pneumoniae, and Escherichia coli.
Fig. 6
Integrated correlation network linking host factors, microbial species, and functional gene families. Correlation network plot showing the relationships among metadata, microbial species, and gene family (KO or GH). Each node represents the gene families (KO/GH), microbial species, and metadata variable. The shapes and colors of the nodes indicate these three categories, while KOs associated with three compounds (acetate, butyrate, and propionate) are highlighted in distinct colors. Edge colors represent correlation coefficients; edge widths correspond to -log10(P-values) derived from Spearman’s or Pearson correlation analysis.
Click here to Correct
Conversely, alcohol consumption was positively correlated with abundance of Parabacteroides distasonis and Bacteroides nordii and was associated with mucin degradation–related KOs (K00854, K00134) and GH families (GH2, GH3, GH5, GH30, and GH95). These taxa formed a cluster with Blautia hansenii but showed relatively weak associations with SCFA-related gene clusters. At the metadata level, I-FABP was positively correlated with L-FABP and alcohol intake but not with spicy food intake. Spicy food intake was positively correlated with alcohol consumption, and the two factors were jointly associated with abundance of Fusobacterium mortiferum.
DISCUSSION
In this study, we evaluated how spicy food consumption and alcohol intake differentially influence intestinal epithelial injury, dietary and lifestyle factors, and gut microbial composition in healthy adults. Among the epithelial injury markers, I-FABP exhibited significant correlations with age, BMI, stool frequency, and sex, indicating that subtle epithelial damage signals are strongly influenced by demographic and physiological factors (Fig. 1B). In contrast, spicy food consumption frequency showed no association with either I-FABP or L-FABP, and no linkage was observed between FABP-based injury markers and spicy sensitivity or pain responses. These findings suggest that typical levels of spicy food consumption in healthy adults do not induce noticeable structural intestinal damage.
However, alcohol intake showed a clear positive correlation with I-FABP, consistent with previous reports linking alcohol to epithelial cytotoxicity, reduced mucus secretion, and weakened barrier function41,42. In contrast, the L-FABP levels did not correlate with alcohol consumption. This pattern suggests that alcohol consumption may directly trigger the release of I-FABP, which is exclusively expressed in small intestinal enterocytes, thereby providing a more specific indication of alcohol-induced enterocyte disruption. As L-FABP is expressed across multiple organs, its association with alcohol may be less readily captured in a cross-sectional setting, potentially limiting its sensitivity to alcohol-specific epithelial disturbances19,43.
Urinary inflammatory cytokine and oral epithelial TRPV1 expression were assessed to characterize systemic inflammation and receptor expression44,45, but both were detected at low concentrations across samples (Supplementary Fig. 1A). This is consistent with previous reports that healthy adults have low basal cytokine levels and limited detectability45,46, reinforcing I-FABP as a relatively stable epithelial injury marker under the conditions of this study.
A
Based on shotgun metagenomic analyses, distinct species-level differences were observed when participants were grouped separately by spicy food or alcohol intake (Fig. 2A–B). In the high-spicy food group (SH), increases were observed in Blautia, Limosima, Ruminococcus, Butyribacter, and Butyricicoccus, taxa associated with SCFA production and mucin metabolism4,22–25. These findings suggest that spicy food intake may selectively influence the mucus layer-associated metabolic pathways, supporting hypotheses from studies showing capsaicin-induced stimulation of mucin production4. Conversely, increases in Fusobacterium_A mortiferum, Dialister, and species of Dorea, taxa with context-dependent or bidirectional health associations26–31, demonstrate that spicy food may act as a complex stimulus capable of eliciting both beneficial and adverse signals, in line with the mixed findings of previous research68.
Changes observed in the high-alcohol-intake group (DH) differed markedly. Although some taxa typically considered beneficial, including Butyricicoccus, Lawsonibacter, and Weissella, were enriched32,34,35, this was in contrast to the DL group, which showed a greater abundance of SCFA-producing and mucin-degrading taxa such as Alistipes and Blautia38. Notably, the epithelial injury marker I-FABP was not different between spicy food groups, but was significantly elevated only in the high-alcohol-intake group (p = 0.004), supporting the interpretation that alcohol exerts a more direct effect on epithelial injury than spicy food. Thus, the presence of beneficial taxa in the DH group may reflect ecological persistence under alcohol-derived stress, rather than improved gut function36,37. Taken together, our results indicate that the intake of spicy food induces variable and mixed shifts within metabolic and mucus-related microbial networks, whereas alcohol intake produces distinct alterations, particularly in pathways linked to mucin metabolism. These findings indicated that spicy foods and alcohol may influence the gut ecosystem through mechanistically divergent routes.
A
When participants were categorized into four combined spicy food–alcohol groups (DLSL, DLSH, DHSL, and DHSH), the overall microbial diversity did not differ significantly across groups, and Bray–Curtis PCoA analyses showed no distinct clustering (Supplementary Fig. 1B–C).
A
Nonetheless, functional pathways related to mucin metabolism and SCFA production were clearly differentiated according to spicy food intake, suggesting that although the global taxonomy remains stable, metabolic network configurations diverge more strongly when spicy food and alcohol exposure intersect (Fig. 3, Fig. 4, Fig. 5). The high-spicy groups (DHSH and DLSH) were enriched in the genera Blautia, Roseburia, Coprococcus, and Limisoma, taxa known to support combined fiber–mucin utilization and SCFA production4,22–24. Their presence aligns with the previously defined fiber–mucin coupling metabotypes 1–2, which represent homeostatic states sustaining both mucin turnover and epithelial energy supply40. This interpretation was reinforced by functional metagenomic analyses, which showed differential activation of mucin-degrading pathways and SCFA-related gene families in high-spicy groups (Figs. 4 and 5A–B). These results support the hypothesis that spicy foods exert compensatory or modulatory effects on gut metabolic networks that are otherwise perturbed by alcohol exposure. Given that mucin turnover critically relies on cross-feeding among SCFA-producing bacteria47, the observed patterns suggest that spicy food intake, regardless of alcohol consumption, modulates the abundance of genes related to butyrate production, including buk (K00929) and bup (K01034) (Fig. 5B). In particular, the Blautia producta-associated butyrate kinase gene buk (K00929) was enriched in the DLSH group, whereas the Fusobacterium mortiferum-associated butyrate phosphotransferase gene bup (K01034) was specific to the DHSH group. In the DHSH group, enhanced carbohydrate metabolism (PWY-5676, P461-PWY), sulfate reduction (K00394, K00395, K00956), and glycosidase pathways (GH13–GH5 families), along with β-galactosidase (K01491) and galactokinase (K00965), suggested an active mucin degradation–SCFA production flux. Blautia hansenii and Anaerostipes hadrus appeared to sustain the acetate–propionate flux (Fig. 5A–B), indicating the active turnover of sulfo-mucin-derived sulfate, mucin glycan degradation, and SCFA production40. In contrast, the DHSL group was characterized by enrichment of Clostridium_Q species, known H₂S producers, whereas the DHSH group contained H₂S consumers and SCFA producers such as Anaerostipes hadrus (Fig. 3A, Fig. 5B). This pattern implies a disrupted cross-feeding flux in DHSL, marked by H₂S accumulation and a reduction in H₂S-consuming taxa (Faecalibacterium, Roseburia, Anaerostipes), a condition linked to proinflammatory metabolic processes48. Network analysis consistently showed that alcohol-associated taxa, such as Parabacteroides and Bacteroides, formed mucin-degrading-dominant networks with weak connectivity to the mucin–SCFA module (Fig. 6). These contrasting network structures (Fig. 6) indicate that spicy food may support mucin homeostasis via mucin–SCFA coupling, whereas alcohol promotes mucin degradation under reduced SCFA flux through cross-feeding loss and H₂S activation20,36,49. Furthermore, the negative correlation between TRPV1 and α-fucosidase gene families (GH29, GH95)—which were disconnected from the SCFA gene cluster—aligns with previous reports suggesting that TRPV1 activation indirectly promotes mucin synthesis4,50. Collectively, these findings imply that the maintenance of mucin homeostasis by spicy food may serve as a determinant of resilience or disruption of gut metabolic health in the presence of alcohol exposure.
Nevertheless, the enrichment of Proteobacteria and Fusobacteriota taxa in the SCFA- and mucin-targeted gene families of the DHSH and DLSH groups (Fig. 5A, B) warrants cautious interpretation. Fusobacterium mortiferum, Escherichia coli, and Klebsiella pneumoniae are well-known dysbiosis-associated taxa that can induce epithelial injury via endotoxin production, and are linked to metabolic and inflammatory disorders26,51,52. Genes such as galK (K00849) and galT (K00965) have been suggested to contribute indirectly to LPS biosynthesis53,54, and incomplete reduction of H₂S-generating pathways (K00394, K00395, K00956, and K00957) by consumers in the human digestive tract may represent an additional risk factor55,56. Notably, Proteobacteria are key constituents of mucin profile 3, reflecting disrupted mucin homeostasis40.
In the DHSH group, the mucin- and SCFA-producing pathways were relatively upregulated (Fig. 5B). Although an appropriate level of mucin turnover can strengthen the protective barrier, excessive degradation has been reported to induce inflammatory injury, indicating the need for interpretation based on prior reviews47. The amino acid metabolic pathways enriched in the DLSH group, including methionine, arginine, ornithine, and alanine metabolism (Fig. 4), may represent terminal metabolic routes associated with mucin degradation40. However, hyperactivation of these pathways has also been linked to disturbances in host amino acid homeostasis, excessive production of H₂S and polyamines, and dysregulation of the serotonin–glucose axis, all of which may contribute to gut and metabolic diseases such as inflammatory bowel disease and insulin resistance57,58. Furthermore, acetate-generating and stratified Klebsiella pneumoniae, a species associated with fatty liver risk, may reflect a metabolically inflammatory transitional state59,60. These findings are consistent with previous reports that high concentrations of capsaicin can negatively impact the gut environment9; however, because this study did not include systematic quantification of capsaicin content in spicy foods, additional evidence is required. Nevertheless, given that the study population comprised healthy adults, epithelial injury markers showed alcohol-specific patterns, and inflammatory markers were detected at low levels, the spicy food intake observed here is likely within a range capable of maintaining mucin homeostasis. However, spicy food consumption may exert bidirectional modulation on mucus metabolism and microbial ecology, acting in either protective or detrimental directions, depending on consumption intensity.
This study had several limitations. The participants were healthy adults without clinical diseases; thus, I-FABP fluctuations should be interpreted as physiological variations rather than acute inflammatory responses. The cross-sectional design precludes causal inference, and classification based on consumption frequency rather than on quantitative intake limits the estimation of actual exposure. Additionally, because this was a remote study, blood and tissue analyses were not performed, and mechanistic interpretations related to mucin homeostasis relied on information in the literature rather than direct measurements. An imbalance in group sample sizes may also introduce statistical constraints; however, permutation-based tests and FDR corrections were applied to minimize potential bias. Future studies should incorporate quantitative dietary assessments and longitudinal designs to more clearly define the dose-response relationships between spicy food or alcohol intake and mucus layer function. With respect to spicy foods in particular, integrating quantitative intake data with refined assessments of individual sensitivity and sensory responses may reveal subtle physiological changes not captured by frequency-based classification alone. Moreover, integrated metabolomic and transcriptomic analyses using blood or tissue samples are required to directly verify whether the microbial functional changes observed in this study translate into epithelial injury signals or mucin regulatory pathways. Such approaches would provide a more precise understanding of the mechanisms by which dietary factors modulate gut barrier homeostasis.
CONCLUSIONS
With this study, we offer an analysis of the effects of alcohol and spicy food consumption on intestinal epithelial injury, gut microbial composition, and metabolic pathways in healthy adults. Alcohol intake was associated with elevated I-FABP levels, reflecting microepithelial injury, whereas spicy food intake promoted mucin metabolism and SCFA-linked compensatory pathways, supporting mucus layer stability. However, enriched amino acid metabolism in the DLSH group represents terminal mucin degradation processes that may predispose individuals to metabolic inflammation. In addition, the concurrent enrichment of Proteobacteria and Fusobacteriota in the DHSH group suggests a potential risk of dysbiosis through mucin degradation and endotoxin activation. Together, these findings suggest that spicy food exerts bidirectional regulation on gut mucus metabolism, maintaining mucin homeostasis at moderate levels but potentially aggravating mucus degradation and metabolic imbalance when combined with excessive alcohol exposure. Therefore, the intensity of spicy food intake and the level of alcohol exposure should be considered key determinants of mucin homeostasis and gut metabolic health.
METHODS
Study Population and Design
A
A
This cross-sectional observational study was conducted between 2024 and 2025 and included 240 healthy Korean adults aged 19 years or older, recruited from Seoul and the surrounding metropolitan area.
A
The study protocol was approved by the Institutional Review Board of Korea National Institute for Bioethics Policy (KoNIBP) (IRB No. P01-202411-02-004) and registered with the Clinical Research Information Service of Korea ( trial number KCT0010676).
A
All participants were provided with a detailed explanation of the study objectives and procedures, and gave written informed consent prior to participation.
A
A
All methods were performed in accordance with the relevant guidelines and regulations.
Healthy adults were specifically recruited to exclude potential disturbances in the gut microbiome caused by diseases or medication and to clearly evaluate functional variations associated with dietary factors in a physiologically stable population. The exclusion criteria included a history of clinically diagnosed disease, current medication use, and no prior consumption of Buldak hot chicken flavor ramen (Samyang Foods Inc., Seoul, Korea), as this particular product was used as the benchmark for perception of spiciness. Eligibility screening was performed using self-report questionnaires and verification steps. To minimize sample degradation, recruitment was limited to participants residing within one–two hours of the study site.
Questionnaire Assessment
To control for individual factors that could influence the gut microbiome composition, participants completed a structured questionnaire comprising 24 items grouped into three categories: demographic information (six items), health information (five items), and dietary lifestyle factors (13 items). The questionnaire was expressly developed for this study and was partly based on published survey instruments used in previous microbiome and dietary studies involving Korean populations61,62.
Demographic information included sex (male or female), age group (20s–60s and older), height (cm), weight (kg), and contact information. Health information covered disease diagnosis and medication status within the past three years (including gastrointestinal, hepatic, endocrine, cardiovascular, musculoskeletal, and neurological disorders) using a four-tier scale (no diagnosis, past diagnosis, under treatment, and currently medicated). Participants were also asked about recent antibiotic use (within the past month) and the intake of dietary supplements such as vitamins or probiotics. Stool-related questions included bowel movement frequency (≤ 2 times per week to > 1 time per day), stool form (Bristol Stool Scale, types 1–7), and the presence or absence of recurrent abdominal pain.
Dietary and lifestyle items encompassed spicy food preference (10-point extended Likert scale), frequency of spicy food consumption (from never to ≥ 5 times/week), intensity of post-consumption abdominal discomfort (visual analog scale, 0–10), onset of abdominal pain after consumption (from none to > 1 d), and pain duration (from none to ≥ 2–3 days). Additional items assessed dietary tendencies such as meat/vegetable preference (1–10, lower scores indicating meat preference), consumption frequencies of high-fat, instant, or fast foods (from never to ≥ 5 times/week), alcohol consumption (from never to ≥ 5 times/week; drinking amount: 1 to ≥ 7 glasses per occasion), smoking (none to ≥ 1 pack/day), physical activity (0 to ≥ 5 times/week), and self-perceived health, mood, sleep quality, lifestyle regularity, and hygiene (5-point Likert scales).
A
To standardize the perception of “spicy food” across participants, Buldak hot chicken flavor ramen (Samyang Foods Inc., Seoul, Korea) was explicitly referenced in all related questions as the benchmark of spiciness. In accordance with IRB recommendations, individuals with no prior consumption experience of this product were excluded from participation.
Alcohol Intake Assessment
Daily alcohol intake (g/day) was calculated based on the self-reported drinking frequency and the average amount of alcohol consumed per occasion obtained from the questionnaire. The calculation was standardized using soju, the most commonly consumed alcoholic beverage in Korea (16% alcohol by volume), as a reference. The amount of pure ethanol contained in one glass (50 mL) was calculated as (50 mL × 0.16 × 0.789 g/mL [ethanol density]), equivalent to 6.3 g of alcohol. Accordingly, the daily alcohol intake of each participant was calculated as: 6.3 g × (the average number of glasses consumed per day).
Drinking frequency reported on a weekly or monthly basis was converted to daily frequency by dividing by 7 or 30, respectively, and multiplying by the average number of drinks per occasion to yield g/day. For nondrinkers, pseudo-counts were applied to prevent distortion by zero values during statistical analysis. This calculation method was adapted from the standardized procedures used in the Korea National Health and Nutrition Examination Survey (KNHANES)63.
Group Classification
For spicy food consumption, participants whose intake frequency exceeded the median value (once per week) were classified into the high engagement group (SH), whereas those at or below the median were assigned to the low engagement group (SL). Alcohol intake groups were defined using the same approach: participants whose estimated alcohol consumption exceeded the median value (2.91 g/day) were categorized as the high-engagement group (DH) and those at or below the median as the low-engagement group (DL).
The classification of spicy food and alcohol intake groups was based on the quartile distribution of each survey variable, with the median used as the cutoff. Participants whose values exceeded the median were placed in the high group, and those with corresponding values at or below the median were placed in the Low group. Accordingly, participants were assigned to one of the following four groups: DLSL (Drink-Low–Spicy-Low), DHSL (Drink-High–Spicy-Low), DLSH (Drink-Low–Spicy-High), and DHSH (Drink-High–Spicy-High).
Sample Collection and Storage
The participants self-collected fecal (approximately 3 g), urine (approximately 25 mL), and oral epithelial cell (approximately 0.1 mL) samples using sterile collection kits provided by the study team.
A
Sampling was performed at home or in the workplace, and specimens were immediately transported via rapid local delivery to the research facility. All samples were stored at − 80°C until analysis to minimize degradation.
DNA Extraction and Shotgun Metagenomic Sequencing
Genomic DNA was extracted using the DNeasy PowerSoil Pro Kit (Qiagen, Hilden, Germany). DNA quality and concentration were verified before shotgun metagenomic sequencing on the Illumina NovaSeq platform with 150-bp paired-end reads. The quality of the raw sequencing reads was initially assessed using FastQC to ensure the overall read quality before downstream processing.
All steps for metagenome sequencing and preprocessing followed the bioBakery workflow (https://github.com/biobakery/biobakery/wiki/biobakery_workflows). Subsequent preprocessing was performed using KneadData (https://github.com/biobakery/kneaddata), which consisted of three main steps: (1) Quality trimming was performed using Trimmomatic, applying a sliding window approach to remove reads when four consecutive bases had Phred scores below 20. Reads shorter than 60 bp were excluded from analysis. (2) Adapter sequences were removed using Trimmomatic. (3) Host-derived and common sequencing components were decontaminated using Bowtie2 by aligning reads against the human reference genome (hg38).
Taxonomic and Functional Profiling
Microbial taxonomic and functional profiles were generated using MetaPhlAn 4 (v4.1.1) and HUMAnN 3.6, respectively. Functional annotations included MetaCyc pathways and KEGG Orthologs (KOs). In the taxonomic analysis, the species-level genome bin (SGBs) output by MetaPhlAn 4 was converted into a GTDB-based taxonomy using the sgb_to_gtdb_profile.py script provided by the MetaPhlAn 4 suite. Downstream analysis retained only taxa or functional features (MetaCyc pathway or KO) that were present at a minimum relative abundance of 0.005% in all samples. Alpha diversity was calculated using the Shannon index and the phyloseq R package. For beta diversity, Bray–Curtis dissimilarity was calculated using the phyloseq R package based on relative abundance data. Principal coordinate analysis (PCoA) was performed using the Bray–Curtis distance matrix. To test for class-level differences in community composition, a permutational multivariate analysis of variance (PERMANOVA) was conducted using 999 permutations, with sample class (DLSL, DLSH, DHSL, and DHSH) as the main factor and BMI, sex, and age as covariates.
Linear Mixed Model Analysis
To identify microbial taxa and functional features associated with metadata, linear mixed models were fitted using the MaAsLin2 R package. The fixed effects included the sample class, BMI, sex, and age. Only taxa and functional features (pathways or gene families) with a prevalence of at least 5% across all samples were included. P-values were corrected for multiple testing using the Benjamini-Hochberg method, as implemented in MaAsLin2.
Microbial Network and Correlation Analysis
To explore the associations between microbial species, functional features, and metadata, Spearman’s rank correlations (for metadata) and Pearson’s correlations (for continuous microbial profiles) were computed. Microbial interaction networks were visualized using Cytoscape (v3.10.3).
Measurement of Urinary Biomarkers
Urinary samples were analyzed to quantify intestinal epithelial injury markers, including intestinal fatty acid–binding protein (I-FABP) and liver fatty acid–binding protein (L-FABP), as well as inflammatory markers interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α). Intestinal epithelial injury markers were measured using the Human I-FABP and L-FABP ELISA Kits (R&D Systems, MN, USA), and inflammatory markers were quantified using the Human IL-6 and TNF-α ELISA Kits (R&D Systems, MN, USA).
Urine samples were thawed gradually at 4°C overnight prior to analysis, and a single freeze–thaw cycle was performed to ensure consistency across all measurements. Before quantification, samples were centrifuged at 500 g for 15 min at 4°C, and the supernatant was collected to remove debris. All measurements and quantification were performed according to the manufacturer’s ELISA protocol.
Analysis of TRPV1 Expression in Oral Epithelial Cells
TRPV1 expression in oral epithelial cells was measured by collecting tongue epithelial samples and quantifying TRPV1 and GAPDH levels using the Human TRPV1 ELISA Kit (Antibodies.com, Cambridge, UK) and Human GAPDH ELISA Kit (Invitrogen, CA, USA). The relative TRPV1 expression was calculated as the TRPV1-to-GAPDH ratio for intergroup comparisons.
The collected epithelial cells were dissolved in 1 mL cold radioimmunoprecipitation assay (RIPA) lysis buffer (Thermo Fisher Scientific, MA, USA) containing protease inhibitors. Cell disruption and protein extraction were performed using a bead-beating grinder (MP Biomedicals, CA, USA) at 4 m/s for 20 s and repeated twice using Lysing Matrix E tubes. Samples were centrifuged at 500 g for 15 min at 4°C to remove debris, and all procedures were conducted on ice. TRPV1 and GAPDH levels were quantified using ELISA, according to the manufacturer’s protocols, and TRPV1 values were normalized to GAPDH levels to account for sample variability.
Statistical Analysis
All statistical analyses were conducted in R (The R Project for Statistical Computing; v4.4.2; https://www.r-project.org) using RStudio (v2024.12.1.563). For multiple comparisons involving microbiome and functional gene analyses, false discovery rate (FDR) correction was performed using the Benjamini–Hochberg method. Following the analytical threshold applied by Huang et al. (2024)64, an adjusted FDR < 0.25 was considered statistically significant for exploratory multi-omics association analyses.
ABBREVIATIONS
SCFA Short
chain fatty acid
I
FABP Intestinal-fatty acid-binding protein
L
FABP Liver-fatty acid-binding protein
TNF
αTumor necrosis factor-α
IL
6 Interleukin-6
TRPV1 Transient receptor potential vanilloid 1
KO KEGG ortholog
GH Glycoside hydrolase
BMI Body mass index
PERMANOVA Permutational multivariate analysis of variance
Acknowledgements
This work was supported by the Research Center of the Mitomics Institute, Samyang Foods Inc.
A
Author Contribution
U. M., Jinwook Kim, and W. L. conceptualized the study. Jihyun Kim developed the analytical pipeline, performed the statistical analyses, and prepared all visualizations. U.M., H.J., and Jihyun Kim curated and validated the metagenomic data. Jinwook Kim, H.O., and S.A. conducted the biomarker analyses. U.M. and Jihyun Kim interpreted the results and wrote the manuscript with input from all the authors. H.S. and W.L. supervised the study. All authors reviewed and approved the final version of the manuscript.
A
Data Availability
The metagenomic sequencing data generated in this study have been deposited in the NCBI Sequence Read Archive (SRA) under the accession number PRJNA1366038. During the peer-review process, the data are accessible to editors and reviewers via the following confidential reviewer access link provided by NCBI:https://dataview.ncbi.nlm.nih.gov/object/PRJNA1366038?reviewer=qnor0m05t9ccvactnnah81gcui.Due to the ethical restrictions imposed by the Institutional Review Board of the Korea National Institute for Bioethics Policy (KoNIBP) (IRB No. P01-202411-02-004), the individual-level demographic and questionnaire data cannot be made publicly available.
A
Anonymized data may be obtained from the corresponding author upon reasonable request.
A
Funding
This work was supported by internal research funding from Samyang Foods Inc. The funder had no role in the study design, data collection, analysis, interpretation, manuscript preparation, or the decision to submit the article for publication.
A
Competing Interests
This work was supported by internal research funding from Samyang Foods Inc. The funder had no role in the study design, data collection, analysis, interpretation, manuscript preparation, or the decision to submit the article for publication. All authors are employees of Samyang Foods Inc.
Electronic Supplementary Material
Below is the link to the electronic supplementary material
References
1.
Paone, P. & Cani, P. D. Mucus barrier, mucins and gut microbiota: the expected slimy partners? Gut 69, 2232–2243 (2020).
2.
Luis, A. S. & Hansson, G. C. Intestinal mucus and their glycans: A habitat for thriving microbiota. Cell. Host Microbe. 31, 1087–1100 (2023).
3.
Shen, W. et al. Anti-obesity effect of capsaicin in mice fed with high-fat diet is associated with an increase in population of the gut bacterium Akkermansia muciniphila. Front. Microbiol. 8, 272 (2017).
4.
Kumar, V. et al. Mucin secretory action of capsaicin prevents high fat diet-induced gut barrier dysfunction in C57BL/6 mice colon. Biomed. Pharmacother. 145, 112452 (2022).
5.
Wang, Y., Tang, C., Tang, Y., Yin, H. & Liu, X. Capsaicin has an anti-obesity effect through alterations in gut microbiota populations and short-chain fatty acid concentrations. Food Nutr. Res. 64, 10–29219 (2020).
6.
Xiang, Y. et al. Beneficial effects of dietary capsaicin in gastrointestinal health and disease. Exp. Cell. Res. 417, 113227 (2022).
7.
Zhang, W. et al. The effects of capsaicin intake on weight loss among overweight and obese subjects: a systematic review and meta-analysis of randomised controlled trials. Br. J. Nutr. 130, 1645–1656 (2023).
8.
Xiang, Q. et al. Capsaicin—the spicy ingredient of chili peppers: A review of the gastrointestinal effects and mechanisms. Trends Food Sci. Technol. 116, 755–765 (2021).
9.
Xiang, Q. et al. Capsaicin, the spicy ingredient of chili peppers: effects on gastrointestinal tract and composition of gut microbiota at various dosages. Foods 11, 686 (2022).
A
10.
Hammer, J. & Vogelsang, H. Characterization of sensations induced by capsaicin in the upper gastrointestinal tract. Neurogastroenterol Motil. 19, 279–287 (2007).
11.
Lejeune, M. P., Kovacs, E. M. & Westerterp-Plantenga, M. S. Effect of capsaicin on substrate oxidation and weight maintenance after modest body-weight loss in human subjects. Br. J. Nutr. 90, 651–659 (2003).
12.
Lv, J. et al. Consumption of spicy foods and total and cause specific mortality: population based cohort study. BMJ 351 (2015).
13.
Bonaccio, M. et al. Chili pepper consumption and mortality in Italian adults. J. Am. Coll. Cardiol. 74, 3139–3149 (2019).
14.
Alatorre-Cruz, J. M. et al. Traditional Mexican food: Phenolic content and public health relationship. Foods 12, 1233 (2023).
15.
Wang, M., Huang, W. & Xu, Y. Effects of spicy food consumption on overweight/obesity, hypertension and blood lipids in China: a meta-analysis of cross-sectional studies. Nutr. J. 22, 29 (2023).
16.
Knight, R. et al. Best practices for analysing microbiomes. Nat. Rev. Microbiol. 16, 410–422 (2018).
17.
Jovel, J. et al. Characterization of the gut microbiome using 16S or shotgun metagenomics. Front. Microbiol. 7, 459 (2016).
18.
Franzosa, E. A. et al. Species-level functional profiling of metagenomes and metatranscriptomes. Nat. Methods. 15, 962–968 (2018).
19.
Sosnowski, K. & Przybyłkowski, A. Ethanol-induced changes to the gut microbiome compromise the intestinal homeostasis: A review. Gut Microbes. 16, 2393272 (2024).
20.
Bishehsari, F. et al. Alcohol and gut-derived inflammation. Alcohol Res. Curr. Rev. 38, 163 (2017).
1.
21. Structure, function and diversity of the healthy human microbiome. Nature486, 207–214 (2012).
2.
22. Parada Venegas, D. et al. Short chain fatty acids (SCFAs)-mediated gut epithelial and immune regulation and its relevance for inflammatory bowel diseases. Front. Immunol.10, 277 (2019).
3.
23. Singh, V. et al. Butyrate producers, “the sentinel of gut”: Their intestinal significance with and beyond butyrate, and prospective use as microbial therapeutics. Front. Microbiol.13, 1103836 (2023).
4.
24. Zhu, Y. et al. Exploration of the Muribaculaceae family in the gut microbiota: diversity, metabolism, and function. Nutrients16, 2660 (2024).
5.
25. Zou, Y. et al. Butyribacter intestini gen. nov., sp. nov., a butyric acid-producing bacterium of the family Lachnospiraceae isolated from human faeces, and reclassification of Acetivibrio ethanolgignens as Acetanaerobacter ethanolgignens gen. nov., comb. nov. Syst. Appl. Microbiol.44, 126201 (2021).
6.
26. Deng, J. et al. Fusobacterium mortiferum and its metabolite 5-aminovaleric acid promote the development of colorectal cancer in obese individuals through Wnt/β-catenin pathway by DKK2. Gut Microbes17, 2502138 (2025).
7.
27. Rashidi, A. et al. No evidence for colonization of oral bacteria in the distal gut in healthy adults. PNAS118, e2114152118 (2021).
8.
28. Shi, H. et al. The gut microbiome as mediator between diet and its impact on immune function. Sci. Rep.12, 5149 (2022).
9.
29. Anthamatten, L. et al. Stratification of human gut microbiomes by succinotype is associated with inflammatory bowel disease status. Microbiome12, 186 (2024).
10.
30. Petrov, V. et al. Analysis of gut microbiota in patients with Parkinson’s disease. Bull. Exp. Biol. Med.162, 734–737 (2017).
11.
31. Liang, X. et al. Gut microbiome, cognitive function and brain structure: a multi-omics integration analysis. Transl. Neurodegener.11, 49 (2022).
12.
32. Sakamoto, M. et al. Draft genome sequence of Lawsonibacter asaccharolyticus… Genome Announc. 6, 10–1128 (2018).
13.
33. Selma, M. V. et al. Description of urolithin production capacity from ellagic acid of two human intestinal Gordonibacter species. Food Funct. 5, 1779–1784 (2014).
14.
34. Geirnaert, A. et al. Butyricicoccus pullicaecorum, a butyrate producer with probiotic potential, is intrinsically tolerant to stomach and small intestine conditions. Anaerobe30, 70–74 (2014).
15.
35. Moon, Y. et al. Intracellular lipid accumulation inhibitory effect of Weissella koreensis OK1-6… J. Appl. Microbiol.113, 652–658 (2012).
16.
36. Chen, G. et al. Gut microbiota dysbiosis: The potential mechanisms by which alcohol disrupts gut and brain functions. Front. Microbiol.13, 916765 (2022).
17.
37. Gloor, G. B. et al. Microbiome datasets are compositional: and this is not optional. Front. Microbiol.8, 2224 (2017).
18.
38. Glover, J. S., Ticer, T. D. & Engevik, M. A. Characterizing the mucin-degrading capacity of the human gut microbiota. Sci. Rep.12, 8456 (2022).
19.
39. Zhang, J. et al. Altered synthesis of genes associated with short-chain fatty acids in the gut of patients with atrial fibrillation. BMC Genomics22, 634 (2021).
20.
40. Labarthe, S. et al. Four functional profiles for fibre and mucin metabolism in the human gut microbiome. Microbiome11, 231 (2023).
21.
41. de Jong, W.–J. et al. The effect of acute alcohol intoxication on gut wall integrity in healthy male volunteers. Alcohol49, 65–70 (2015).
22.
42. Kokesova, A. et al. The intestinal fatty acid-binding protein as a marker for intestinal damage in gastroschisis. PLoS One14, e0210797 (2019).
23.
43. Voth, M. et al. I-FABP and L-FABP are early markers for abdominal injury… Clin. Chem. Lab. Med.53, 771–780 (2015).
24.
44. Chen, H. et al. Urinary interleukin-6 as a diagnostic and prognostic marker of acute kidney injury in dogs. Vet. Rec.191, no–no (2022).
25.
45. Takahashi, N. et al.
A
Epithelial TRPV1 channels: expression, function, and pathogenicity in the oral cavity. J.
A
Oral Biosci.
62, 235–241 (2020).
26.
46. Biancotto, A. et al. Baseline levels and temporal stability of 27 multiplexed serum cytokine concentrations… PLoS One 8, e76091 (2013).
27.
47. Yamaguchi, M. & Yamamoto, K. Mucin glycans and their degradation by gut microbiota. Glycoconj. J.40, 493–512 (2023).
28.
48. Marcelino, V. R. et al. Disease-specific loss of microbial cross-feeding interactions in the human gut. Nat. Commun.14, 6546 (2023).
29.
49. Herrnreiter, C. J. et al. Bacterial dysbiosis and decrease in SCFA correlate with intestinal inflammation following alcohol intoxication and burn injury. eGastroenterology 3 (2025).
30.
50. Kumar, V. et al. Role of TRPV1 in colonic mucin production and gut microbiota profile. Eur. J. Pharmacol.888, 173567 (2020).
31.
51. Li, Z. et al. Fusobacterium in the microbiome: from health to disease across the oral–gut axis and beyond. npj Biofilms Microbiomes11, 200 (2025).
32.
52. Shin, N.–R., Whon, T. W. & Bae, J.–W. Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 33, 496–503 (2015).
33.
53. Tang, M. et al. Correlation assessment among clinical phenotypes… Hum. Mutat.33, 1107–1115 (2012).
34.
54. Debouck, C. et al. Structure of the galactokinase gene of Escherichia coli… Nucleic Acids Res. 13, 1841–1853 (1985).
35.
55. Braccia, D. J. et al. Hydrogen sulfide (H2S) via cysteine degradation is ubiquitous in the human gut microbiome. Front. Microbiol.12, 705583 (2021).
36.
56. Luo, W. et al. Microbial assimilatory sulfate reduction-mediated H2S: An overlooked role in Crohn’s disease. Microbiome12, 152 (2024).
37.
57. Chen, Y. & Fang, J.–Y. Colonic microbiota amino acid metabolism in gut health regulation. Cell Insight 100227 (2025).
38.
58. Li, T.–T. et al. Microbiota metabolism of intestinal amino acids impacts host nutrient homeostasis and physiology. Cell Host Microbe32, 661–675 (2024).
39.
59. Samuelson, D. R. et al. Alcohol-associated intestinal dysbiosis impairs pulmonary defense against Klebsiella pneumoniae. PLoS Pathog. 13, e1006426 (2017).
40.
60. Shen, M. et al. Alcohol-induced gut microbiome dysbiosis enhances the colonization of Klebsiella pneumoniae. mSystems9, e00052–24 (2024).
41.
61. Falony, G. et al. Population-level analysis of gut microbiome variation. Science352, 560–564 (2016).
42.
62. Shin, J.–H. et al. Differential effects of typical Korean versus American-style diets on gut microbial composition… Nutrients 11, 2450 (2019).
43.
63. Hong, J. W., Noh, J. H. & Kim, D.–J. High-risk alcohol consumption in Korean adults. PLoS One12, e0175299 (2017).
44.
64. Huang, K. D. et al. Dietary selective effects manifest in the human gut microbiota from species composition to strain genetic makeup. Cell Rep. 43 (2024).
Supplementary Figure & Table legends
Supplementary Fig. 1. Biomarkers distribution with different four groups. Boxplots showing the expression levels of inflammatory markers (TNF-α, IL-6) and TRPV1 according to spicy-alcohol intake groups. P-values were calculated using the Kruskal-Wallis test.
Supplementary Fig. 2. Overall composition and diversity of the gut microbiota. A) The phylogenetic tree based on GTDB of 1,527 microbial species detected across all 299 samples. Each color represents a distinct phylum-level taxonomy. B) Alpha diversity represented by the Shannon index according to different dietary groups. P-values were calculated using the Kruskal-Wallis test. C) Principal coordinate analysis (PCoA) based on Bray-Curtis distances illustrating β-diversity among microbial communities across all samples, colored by the different dietary intake groups. P-value was calculated using the PERMANOVA test.
Supplementary Table 1. Baseline demographic, dietary, and lifestyle characteristics of study participants according to spicy-food and alcohol-intake groups. The variables included sex distribution, age, BMI, stool frequency, Bristol Stool Scale score, spicy-food-related responses, dietary preferences, lifestyle habits, and urinary epithelial injury biomarkers (I-FABP and L-FABP). Group differences were evaluated using appropriate statistical tests, and P-values are presented.
Supplementary Table 2. Spearman rank-correlation coefficients between metadata variables and intestinal epithelial-injury biomarkers (I-FABP and L-FABP)
. The table summarizes pairwise associations (ρ) among demographic factors, stool characteristics, spicy-food responses, dietary habits, and biomarker concentrations, along with corresponding P-values.
Total words in MS: 8787
Total words in Title: 21
Total words in Abstract: 227
Total Keyword count: 9
Total Images in MS: 6
Total Tables in MS: 0
Total Reference count: 20