What this is
- This research examines how () affects intestinal health in the context of alcohol consumption.
- disrupts the peripheral circadian clock, potentially worsening the intestinal barrier's response to alcohol-induced damage.
- The study utilizes mouse models to explore the relationship between eating patterns, alcohol intake, and gut microbiota.
Essence
- decreases the resiliency of the intestinal barrier to alcohol-induced injury. This effect is linked to and reduced short-chain fatty acid levels.
Key takeaways
- exacerbates alcohol-induced disruption of intestinal barrier integrity. This is evidenced by decreased expression of tight junction proteins, which are crucial for maintaining gut permeability.
- Alcohol consumption alters the gut microbiota, reducing beneficial short-chain fatty acid-producing taxa. This dysbiosis may contribute to the impaired intestinal barrier function observed in conditions.
- Microbial metabolites from and alcohol-fed mice negatively impact the resiliency of the intestinal epithelial barrier. This suggests that the gut microbiota plays a significant role in mediating the effects of on gut health.
Caveats
- The study's findings are based on animal models, which may not fully translate to human physiology. Further research is needed to confirm these effects in human populations.
- Serum and stool samples were collected at a single time point, which may limit the understanding of temporal changes in markers of intestinal health.
- The caloric intake was not controlled during alcohol treatment, which could introduce variability in the results related to dietary influences on gut health.
Definitions
- wrong-time eating (WTE): Eating during physiological rest periods, misaligning the body's circadian rhythms.
- microbiota dysbiosis: An imbalance in the microbial community in the gut, often linked to health issues.
- short-chain fatty acids (SCFAs): Fatty acids with fewer than six carbon atoms, produced by gut bacteria and important for gut health.
AI simplified
Introduction
Circadian rhythms are endogenous 24-h oscillations controlled by biological clocks that are involved in functions such as the sleep-wake cycle and multiple physiologic activities.1 While the central circadian clock is entrained by light, the peripheral circadian clock has autonomous circadian control and is entrained by other specific factors.2 In the gastrointestinal tract (GIT), the main entrainer, or zeitgeber, is the time of food consumption.3 Therefore, while abnormal exposure to light as in shift work disrupts central circadian rhythms, altered eating patterns during physiologic rest time, defined as wrong-time eating (WTE), disrupts peripheral circadian rhythms and misaligns the peripheral circadian clock of the GIT and the body’s central clock.
Modern lifestyle has led to an increased prevalence of a 24-h society where several aspects of the modern “Western” societies commonly disrupt circadian rhythms. This includes shift work, increased nighttime light exposure (e.g., use of light emitting devices while in bed, bright street lights), social jet lag (different sleep patterns on work days, and free days), frequent long-haul travels over two time zones, irregular eating patterns (e.g., eating at different times each day, eating late at night, eating large meals close to biological rest time), and high-fat diet composition.4–11 While all these lifestyle features contribute to circadian disruption, shift work has the most impactful and dramatic effect on the central circadian clock.12 However, since irregular eating patterns are extremely common then it is highly plausible that “wrong-time” eating could be even more impactful than shift work.
The GI tract is the largest interface between the host and environment and is one of the primary gatekeepers to protect the host from the deleterious effects of environmental factors and maintain a healthy state; conversely, if the “gatekeeping” function of the GI tract is compromised, it could be the site for injurious environmental factors to trigger and/or worsen disease states. The intestinal epithelial barrier is one of the critical elements of this gatekeeping function of the GI tract.13 A network of cell-cell junctions including tight junction proteins, adherens junctions, and desmosomes maintains this barrier. Disruption of this barrier can result in increased permeability and “leakiness” of the intestinal barrier which allows luminal toxins, bacteria, and their metabolites such as endotoxins (i.e., lipopolysaccharide [LPS]) to leave the gut lumen and enter the bloodstream to promote inflammation-mediated diseases and mediate organ damage.14
One environmental factor that can disrupt intestinal barrier integrity and promote endotoxemia is unhealthy alcohol consumption.15 Unhealthy alcohol consumption has been associated with a wide range of harmful health effects. Indeed, alcohol use disorder is associated with chronic liver disease, pancreatitis, muscle/bone disorders, and brain disorders with cognitive decline.16,17 Inflammation is the common underlying mechanism of alcohol-induced organ damage. We and others have shown that the primary source of this inflammatory state in the gut is alcohol-induced disruption of intestinal barrier integrity.18 Intriguingly, not all patients with alcohol use disorder develop organ damage suggesting unhealthy alcohol consumption is required but not sufficient to cause organ damage and additional factor(s) is required for unhealthy alcohol consumption to cause severe enough disruption of the intestinal barrier to initiate and sustained enough inflammatory state to cause organ damage.19 One such environmental factor is disrupted circadian homeostasis because we and others have shown that disrupted circadian can cause microbiota dysbiosis and intestinal leak.20,21
Studies evaluating the potential impact of disrupted circadian on the GI tract and intestinal barrier have primarily focused on central circadian disruption. In rodents, environmental central circadian disruption models (e.g., chronic phase shifts of the light/dark cycle).21 Similarly, genetic disruption through clock gene mutations (e.g., ClockΔ19/Δ19) in mice has worsened intestinal permeability that exacerbates alcohol-induced gut leakiness as well as dysbiosis (resulting in lower microbial diversity) in the presence of a stressor such as alcohol.22,23
In addition, knock-out models of the core clock machinery with Per2, Reverb-a, and Bmal1 have all shown the circadian clock mediates colonic inflammation, dysbiosis, and decreased production of bacterial-derived metabolites like short-chain fatty acids.24–26 In humans, alcohol consumption induces a systemic pro-inflammatory response and harms gut barrier integrity in night workers with circadian misalignment as well as causes dysbiosis of the gut microbiota.12,27
Yet, disruption of the intestinal circadian clock might be more relevant in the regulation of intestinal barrier integrity. However, there are very limited studies to determine whether disruption of intestinal circadian homeostasis by wrong-time eating can negatively impact intestinal barrier function and whether it can worsen alcohol-induced gut leak and endotoxemia. Our study aimed to fill this knowledge gap by utilizing both in vivo animal models and ex vivo organoid models.
Results
Previously, we have shown that central circadian disruption through chronic LD shifting can increase the variability of the circadian period and phase in colonic tissue and organoids before altering food timing to impact peripheral colonic circadian rhythms.21 Food timing may elicit peripheral circadian phase-shifting effects and act as an environmental entrainer.28
In this study, to establish a phenotype of peripheral circadian disruption, we sought to administer the PER2:LUC mice under a food restriction paradigm defined with either day- (rest-phase= “wrong-time”, 7 AM) or nighttime (active-phase= “right-time”, 7 PM) access to food for 12 h for 5 consecutive weekdays followed by a 2-day ad-libitum access to food to mimic human conditions on the weekend. This food timing paradigm represents abnormal eating patterns that have been seen to shift the peripheral clocks in organs such as the liver and colon.29 In addition, the mice were maintained in a 12:12 cycle that followed the standard 12 h light (during the active phase, 7PM) and 12 h dark (during the rest phase). Together, these components formed the basis of our right-time, wrong-time experimental paradigm.
Impact of wrong-time eating (WTE) and alcohol on peripheral (colon) circadian rhythm
Circadian rhythms of colon tissue. (a) Bioluminescence of PER2:LUC in RTE and WTE colon tissue with control (HO). (b) Bioluminescence of PER2:LUC in RTE and WTE colon tissue with alcohol (EtOH). (c) The period of WTE mice has increasing variability than the periods of RTE in colon tissue. There is a significant difference in the period between RTE and WTE colon tissue. Experiments were performed in triplicate (= 3). Purple bars represent H₂O groups and blue bars represent EtOH groups. Two-way ANOVA (results in box) was conducted, and effects are indicated on each graph when significant: *< 0.05, **< 0.01, and ***< 0.001. 2 n p p p
Impact of wrong-time eating (WTE) and alcohol on food/liquid consumption and weight in mice
Average weekly food and liquid intake were significantly affected by treatment and sex effects. Males consumed more food, as expected, but consumed less liquid compared to female mice. Body weight was assessed and was significantly affected by alcohol treatment and sex effects. Male mice were more affected by alcohol treatment, resulting in decreased average body weight. Average weekly alcohol intake was measured and there was a significant impact of the food timing paradigm on alcohol consumption (Supplemental Figure S1).
Impact of wrong-time eating (WTE) on alcohol consumption in mice
The blood alcohol levels were assessed with significantly elevated alcohol effects noted in the right-time eating alcohol-fed group (RTE EtOH), consistent with blood collected following 12 h of food intake and 20% EtOH (v/v) availability. It is significantly increased compared to the H₂O groups. However, there was not a significantly elevated blood alcohol in the wrong-time alcohol-fed (WTE EtOH) group due to a study constraint. The time of sacrifice occurred in the morning although we anticipate this trend to have continued with significantly increased levels if this group was sacrificed at night when the mice are coming off 12 h of food intake and 20% EtOH (v/v) availability (Supplemental Figure S2).
Impact of wrong-time eating (WTE) and alcohol on intestinal permeability
Effect of altered food timing and alcohol consumption on urinary sugar excretion to assess intestinal barrier integrity. (a) Urinary sucralose exhibited a significant effect of food timing and sex but no interaction. (b) urinary lactulose exhibited a significant effect on food timing but no interaction. (c) Sucralose:lactulose ratio was significantly impacted by alcohol treatment, food timing, and sex. (d) urinary sucrose was significantly impacted by food timing and sex, but there was no interaction. (e) Urinary mannitol exhibited a significant effect on food timing and sex but no interaction. (f) Lactulose:mannitol (LM) ratio was significantly impacted by sex, but there was no difference of food timing or alcohol treatment nor was there an interaction. Between= 7–13 mice/treatment group. Three-way ANOVA (results in box) was conducted, and effects are indicated on each graph when significant: *< 0.05, **< 0.01, ***< 0.001, and ****< 0.0001. n p p p p
Impact of wrong-time eating (WTE) and alcohol on AJC protein expression is decreased in colon tissue
Previous studies of central circadian misalignment resulted in colonic hyperpermeability.20,21 The apical junctional complex (AJC) of intestinal epithelial cell proteins contains the tight junction protein zonula occludens protein 1 (ZO-1), occludin, and the adherens junction protein E-cadherin that are key regulators of gut leakiness. We sought to determine whether our altered food timing paradigm and alcohol treatment resulted in changes in the expression of these key AJC proteins that regulate permeability in colon tissue.
Effect of altered food timing and alcohol consumption on AJC proteins. (a-c) section of colonic tissue labeled with DAPI (blue) and ZO-1 (green). ZO-1 in colon tissue is decreased with altered food timing and alcohol treatment. Between= 19–21 mice/treatment group. (d-f) section of colonic tissue labeled with DAPI (blue) and occludin (green). Occludin in colon tissue decreases with altered food timing and alcohol treatment. (d) Immunofluorescent staining of occludin in colon tissue (left) and cell fluorescent measurement (right), in which alcohol treatment significantly affected occludin expression. (e) Western blot analysis of occludin in its cytoplasmic fraction. Food timing and alcohol treatment effects are significant. (f) Western blot analysis of occludin in its membrane fraction. Food timing, treatment, and their interaction effects are significant. Between= 18–21 mice/treatment group. Three-way ANOVA (results in box) was conducted, and effects are indicated on each graph when significant: *< 0.05, **< 0.01, ***< 0.001, and ****< 0.0001. (g) Section of colonic tissue labeled with DAPI (blue) and E-cadherin (green). E-cadherin in colon tissue is decreased with alcohol treatment. Alcohol treatment effects were significant. Between= 18–21 mice/treatment group. Two-way ANOVA (results in box) was conducted, and effects are indicated on each graph when significant: *< 0.05 and **< 0.01. n n p p p p n p p
Impact of wrong-time eating (WTE) and alcohol on colonic inflammation.
Colon tissue from mice was assessed for any histological evidence of inflammation or tissue damage (Supplemental Table S1). Male mice from the RTE EtOH group experienced increased scores of histological evidence of inflammation and tissue damage. Overall, there is a sex and alcohol treatment effect (p = 0.0195).
Impact of wrong-time eating (WTE) and alcohol on stool microbiota composition
However, we also evaluated the impact of alcohol treatment effect on stool microbial communities (i.e., RTE H₂O + WTE H₂O vs. RTE EtOH + WTE EtOH). While there were no significant differences in alpha diversity observed (i.e., variation within each sample), the analysis of beta diversity (i.e., variation between samples/ groups) revealed a significant impact of alcohol on the microbiota.
PERMANOVA analysis showed no differences between male and female groups. Subsequent microbiota analysis was run with combined male and female data. Alcohol-fed groups (RTE EtOH and WTE EtOH) had significantly different effects. Alcohol treatment drove the global changes to the microbiota profile seen in our group comparisons (Table 2).
PERMANOVA and PERMDISP results are based on Aitchison matrix distances for the multi-amplicon (V4) sequence data at the feature level. Significant values are based on 9,999 permutations and corrected for multiple testing using the Benjamini–Hochberg method (bolded, q < 0.05). Total mice sample sizes: RTE H₂O (n = 18); RTE EtOH (n = 17); WTE H₂O (n = 18); WTE EtOH (n = 14).
We analyzed the distance measure of the samples using Non-metric multidimensional scaling (NMDS) (Supplemental Figure S3). The distance between each group is indicative of each group’s dissimilarity and similarity to each other. The RTE H₂O and WTE H₂O groups are more like each other compared to both the RTE EtOH and WTE EtOH groups.
Phylum and genus taxonomic level differential abundance and compositional analysis of the microbial profiles between mice groups were analyzed with both DESeq2 and centered log-ratio Kruskal–Wallis (CLR-KW). These were compared to Boruta to identify important taxa. CLR-KW identified the same and more taxa compared to DESeq2. Of interest, significant taxa that were identified between H₂O and EtOH groups that were associated with alcohol consumption included: Clostridium sensu stricto 1, Muribaculaeceae, Bifidobacterium, Turicibacter, and Atopobiaceae (all q < 0.05) (Supplemental Excel sheet S1).
Taken together, the richness (alpha diversity) of the intestinal microbiota was not impacted by either food timing or alcohol treatment. However, the differences between the microbial communities (beta diversity) were driven mainly by alcohol treatment. The identified taxa of significance are associated with alcohol consumption, supporting how alcohol alters the microbial communities in our mouse groups.
Impact of wrong-time eating (WTE) and alcohol on the microbiota function
We also assessed the biological pathways present between the RTE H₂O and WTE EtOH cecal content supernatant data using Reactome to predict equivalent human pathways.33 A total of 55 pathways were identified (all p < 0.05) and a few pathways of interest to note are adherens junctions interactions, G-protein-coupled receptors (GPCR) signaling and binding signaling by GPCR, and γ-aminobutyric acid (GABA) signaling (supplemental Excel sheet S2). These pathways are decreased in the WTE EtOH samples which are implicated in alcohol consumption.34
Short-chain fatty acid (SCFA) content in cecal content supernatant of mice.= 4 pooled samples from each respective group, RTE H₂O (blue bars) and WTE EtOH (red bars). Two-way ANOVA was conducted, and the results of multiple comparison tests are indicated on the graph when significant: *< 0.05. n p
Relative abundance of short-chain fatty acid (SCFA)-producing taxa in the microbiota of RTE H₂O and WTE EtOH mice. Left image: overall percent abundance of SCFA-producing taxa is decreased in the WTE EtOH group (red bar) compared to the RTE H₂O (blue bar). Right image: overall percent abundance of SCFA-producing taxa is decreased in the WTE EtOH group (left) and much of the changes are occurring in, an important butyrate producer (right). Lachnospiraceae
Relationship between intestinal microbiota and intestinal permeability and barrier integrity outcomes in mice
Significant correlations were noted between sucralose:lactulose ratio and Turicibacter (p-value = 0.034269, R-value = 0.394355), Occludin relative density and Clostridium sensu stricto 1 (p-value = 0.016971, R-value = −0.50322), ZO-1 relative density and Clostridium sensu stricto 1 (p-value = 0.035408, R-value = −0.42979). These results indicate that Turicibacter and Clostridium sensu stricto 1 may be important mediators of microbiota and epithelial barrier integrity in the context of alcohol consumption and serve as potential targets in future studies (supplemental excel sheet S3).
Effect of food timing and alcohol treatment on serum and stool inflammatory markers
Food timing and alcohol treatment can induce intestinal inflammation that can subsequently lead to peripheral inflammation in organs throughout the body. We assessed a few inflammatory markers in the serum and stool. First, the pro-inflammatory cytokine IL-6 was examined in the serum for systemic inflammation. While there was no significant difference between groups, there was a slight increase in IL-6 levels seen in the alcohol-treated groups. Second, lipopolysaccharide (LPS) is a component in the outer membrane of gram-negative bacteria, and LPS-binding protein (LBP) is a type 1 acute-phase protein that binds to LPS to facilitate an immune response and is a well-accepted marker of intestinal barrier integrity and endotoxemia. There was a relative increase of serum LBP between food timing in both WTE groups. Third, we assessed calprotectin which is present within neutrophils and throughout the body. Increased stool calprotectin levels are an indication of neutrophil migration into gastrointestinal tissue due to inflammatory processes. Although there was not a significant difference with food timing or alcohol treatment, the levels of calprotectin increased with food timing effects (Supplemental Figure S5).
Impact of wrong-time eating (WTE) and alcohol on the organoid barrier integrity.
Organoid permeability increases with alcohol and cecal supernatant treatment. (a) FITC-dextran (green) was analyzed at 20x and analyzed for net fluorescence. The presence of green dye inside the lumen of the organoids is indicative of increased paracellular permeability. (b) Net fluorescence levels are significant when comparing the RTE H₂O control (blue bar) to the treatment groups: RTE H₂O cecal supernatant (green bar,= 0.0102), 0.2% EtOH (red bar,= 0.0006), and WTE EtOH cecal supernatant (purple bar,<0.0001). P-values are indicated on each graph when significant: *< 0.05, **< 0.01, ***< 0.001, and ****< 0.0001. p p p p p p p
AJC protein expression is decreased in colonic organoids.
Previous studies with central circadian misalignment have resulted in colonic hyperpermeability in tissue as well as colonic organoids.20,21 The apical junctional complex (AJC) of intestinal epithelial cell proteins contains the tight junction protein zonula occludens protein 1 (ZO-1), occludin, and the adherens junction protein E-cadherin that are key regulators of gut leakiness. The phenotype demonstrated in our in vivo mouse model was recapitulated in our ex vivo organoid model to compare the colon tissue to the colonic organoids generated from them. We sought to determine whether our altered food timing paradigm and alcohol treatment resulted in changes in the expression of these key AJC proteins that regulate permeability in colonic organoids.
Effect of food timing and alcohol treatment on AJC protein expression in organoids. (a) Immunofluorescent staining of tight junction protein ZO-1 and adherens junction protein E-cadherin (green). (b) No differences in ZO-1 or E-cadherin total fluorescence, however, there is a decrease in the WTE EtOH group.= 5 mice/treatment group. (c) Western blot analysis of ZO-1 relative density. No difference due to low n but a similar trend of decreased ZO-1 expression in the WTE EtOH group. Representative data with= 2 mice/treatment group. (d) Western blot analysis of occludin relative density. No difference due to low n but a similar trend of decreased occludin expression in the WTE EtOH group. The blue bar represents the RTE H₂O group, while the red bar represents the WTE EtOH group. Representative data with= 2 mice/treatment group. N n n
Assessing organoid cellular proliferation
To assess cellular proliferation in the organoids, we stained for the cell proliferation marker Ki-67 (Supplemental Figure 6). Previous organoid studies have shown that alcohol-treated organoids have increased levels of Ki-67.37 Alcohol is known to promote cellular proliferation and the migration of intestinal epithelial cells. This data demonstrates an increase in the number of Ki-67 organoids present in the WTE EtOH group.
Discussion
The key findings of the present study are: 1) wrong-time eating (WTE) induces peripheral circadian disruption in the intestine; 2) peripheral circadian misalignment by WTE decreases alcohol intestinal barrier resiliency; 3) WTE-induced decreased resiliency of intestinal barrier to alcohol injurious effects is mediated through disruption of tight junction proteins – ZO-1 and Occludin; 4) alcohol, not food timing, is the primary trigger for microbiota dysbiosis in our model; 5) WTE and alcohol caused decreased relative abundance of SCFA producers and lower butyrate levels in the colonic lumen; and 6) cecal contents from WTE mice increased sensitivity of organoid barrier to cytokines-induced damage (i.e. recapitulate WTE-induced decreased resiliency of intestinal barrier disruption by alcohol in mice) suggesting that luminal factors like bacterial products mediate WTE effects on intestinal barrier function.
Circadian rhythms in the intestines play a crucial role in regulating various physiological processes including digestion, nutrient absorption, gut motility, and cell proliferation. While these peripheral circadian rhythms synchronize with the central circadian clock, peripheral clocks can behave independently based on certain cues. Feeding and fasting are the main drivers of peripheral circadian rhythms, oscillating with rest and active cycles. However, genetic studies of clock genes—Cry, Per1/2, Rev-erbα, and Bmal—demonstrate their effects on intestinal health. For instance, Bmal mutants which have a loss of rhythms, negatively influence cellular regeneration in the colon, leading to worsened colitis in murine DSS-colitis models.38 In addition to the intestine, perturbations to the circadian clocks have resulted in alterations of the gut microbiota, leading to decreased production of bacterial-derived metabolites like short-chain fatty acids.24–26
Previously, we demonstrated how a central circadian misalignment model can promote an increase in colonic permeability and impair barrier integrity.20–22 In this study, we wanted to see if peripheral circadian misalignment also promotes an increase in colonic permeability and impaired barrier integrity. To disrupt peripheral (colon) circadian homeostasis, we used an altered food paradigm (WTE) with and without alcohol to assess if we successfully disrupted circadian rhythms in the colon by the WTE paradigm.9 Then, we utilized a PER2:LUC gene reporter mouse that expresses a Per2-luciferase fusion protein that allows us to measure both colon tissue and colonic organoid circadian clock rhythmicity.39,40 We found that colon tissue had increased variability in the WTE group. Alcohol did not significantly affect circadian rhythmicity suggesting that altered food timing alone is sufficient to affect peripheral circadian rhythms.
To assess the extent of the impact that peripheral circadian misalignment has on the resiliency of the colon, we assessed measures of intestinal permeability, barrier integrity, gut microbiota profiles, as well as the potential effect of the cecal supernatant of these mice on potentially ameliorating and/or exacerbating intestinal outcomes. Altered food timing and alcohol consumption resulted in increased intestinal permeability and exhibited decreased expression of AJC proteins like ZO-1, occludin, and E-cadherin in colon tissue. These AJC proteins are markers of intestinal permeability and are under circadian regulation and this data supports the findings in previous studies.9,41–43
However, despite the disruption of intestinal AJC proteins and increased intestinal barrier by WTE and alcohol, we could not find evidence of systemic inflammation (normal serum cytokine) or endotoxemia (normal serum LBP). It is possible that the duration of our study was not adequate to cause long enough disrupted barrier function to cause enough intestinal leak to endotoxin to be identified by changes in serum LBP or serum cytokines. Further studies with a longer duration of WTE and alcohol are required to answer this question.
One possible mechanism for disruption of intestinal barrier integrity by WTE/alcohol is altered microbiota composition (dysbiosis). Indeed, we identified significant taxa that are implicated in alcohol consumption such as Clostridium sensu stricto 1, Muribaculaceae, Bifidobacterium, Turicibacter, and Atopobiaceae. We also observed a decrease in the relative abundance of putative anti-inflammatory SCFA-producing taxa, like Lachnospiraceae, an important butyrate producer.23,44 We then used the organoid model to determine whether WTE/alcohol disrupts intestinal barrier integrity more directly and whether luminal dysbiotic bacteria products are responsible for the WTE-induced decreased resiliency of the intestinal barrier to the injurious effect of alcohol. Previous work demonstrated that colon organoids are an excellent model to elucidate the mechanism of the alcohol-induced gut leak or disrupted circadian-induced intestinal epithelial injury because colonic epithelial organoids retain the phenotype of the host.20,21,42 Intestinal organoids are self-renewing three-dimensional tissue cultures that are derived from intestinal stem cells and can be used to test specific hypotheses that are not possible in vivo. They are advantageous compared to traditional cells used for cell culture (e.g., Caco-2 cells) because the organoids contain various primary cell types present in vivo. Indeed, PER2:LUC mice have been a useful model for the characterization of the circadian clock because the insertion of the luciferase reporter in the Period gene (a critical component of the molecular circadian clock) allows for easy assessment of Period gene expression across 24 h.39,40 Previous studies have studied organoid circadian phenotypes in the small intestine, but few studies have examined colonic organoids. In the current study, we used colonic organoids and found AJC proteins were significantly disrupted in organoids from WTE/alcohol-fed mice. This study complements our previous study which found that central circadian disruption disturbs the circadian regulation of apical junction proteins.21 Furthermore, using colonic organoids from control mice (RTE/H2O mice), we found that filtered cecal content from WTE alcohol-fed mice induces marked intestinal leak suggesting that colonic lumen in WTE/alcohol-fed mice contain compound(s) that, at least in part, may be responsible for alcohol-induced disruption of intestinal barrier integrity and perhaps WTE-induced decrease resiliency of intestinal barrier to the injurious effect of alcohol.
We found that WTE/alcohol-fed mice have a dysbiotic microbiota community in their stool characterized by a relative abundance of SCFA-producing bacteria. In our prior studies, we also found stool microbiota dysbiosis characterized by a low relative abundance of SCFA producers in central circadian disrupted humans and mice.12,20–23 It is noteworthy that multiple studies in rodents and humans have shown the important role of SCFA in maintaining normal intestinal barrier function with low SCFA being associated with disrupted intestinal barrier integrity.20,23,45 We therefore measured SCFA levels in the filtered cecal content and found that SCFA, butyrate was significantly decreased in the cecal content of WTE/alcohol-fed mice. Thus, it is highly plausible that low SCFA levels in the cecal content are responsible for intestinal organoid leak induced by the cecal content from WTE/alcohol-fed mice. However, other metabolites/compounds can still be contributing factors to this leak. Further studies including untargeted and targeted metabolomics of the cecal content and direct treatment of organoids with these compounds are required to identify the culprit metabolite(s).
The following are several important limitations of the present study. Firstly, we assessed serum and stool samples at one singular time point in the morning, instead of conducting two (one morning and one nightly collection), at the end of the study. This was a byproduct of the study design which might have impacted the difference between groups in various serum and stool markers (e.g., serum blood alcohol, IL-6, serum LBP, and stool calprotectin). Secondly, while we can recapitulate barrier dysfunction in RTE organoids with cecal contents from WTE mice (demonstrating a loss of function), the assessment of cecal contents from RTE mice on WTE organoids (“gain of function” approach) was not performed. This was in part due to limited WTE organoid availability. Thirdly, the caloric intake was not controlled in our alcohol-in-water model with an isocaloric substitution such as maltodextrin. This decision was based on our prior finding that maltodextrin can also affect the intestinal barrier inflammation readouts and the intestinal microbial community.46 Lastly, we performed a global analysis of the cecal content and did not look at specific bacteria or metabolites that should be examined in additional studies. This research paradigm offers valuable insights into how abnormal, late-night eating and alcohol consumption interact with the gut. Although our study focused on mice, these findings may reflect similar trends in humans, such as night shift workers and in groups where individuals tend to stay up late on a regular basis. Understanding the biological effects of these behaviors could help address potential long-term health risks, such as digestive disorders or metabolic dysfunction. Future research could explore targeted interventions, such as rescue treatment with bacterial-derived metabolites like butyrate from non-circadian disrupted individuals or lifestyle adjustments to mitigate the negative consequences of disrupted eating patterns and alcohol intake.
Conclusion
In summary, our study indicates that: (1) Eating during biological rest time (wrong-time eating [WTE], a model of peripheral circadian misalignment) is sufficient to influence circadian rhythmicity, (2) Circadian misalignment by WTE decreases the resiliency of the colon to the injurious alcohol effect and this effect might be mediated through intestinal microbiota dysbiosis that results in low SCFA in the colonic lumen. Unhealthy alcohol consumption in circadian misaligned hosts is a major contributor to intestinal microbiota dysbiosis and gut leak and identifying the microbial metabolites that mediate gut leak and associated intestinal/systemic inflammation may lead to new microbial and circadian targets to prevent organ damage in patients with alcohol use disorders.
Materials and methods
Animal model
All experiments were approved by the Rush University Institutional Animal Care and Use Committee (IACUC). Mice with a Per2 (mPer2Luciferase [mPer2Luc]) knock-in mutation, obtained from Jackson Labs (Bar Harbor, ME), were bred in-house.39 A total of 89 mice (51 male and 38 female) were used for experiments, with 12–13 males and 9–10 females. All mice were at least 8-weeks old when the experiments began. The animals were maintained on a 12-h light:12-h dark (12∶12 LD) cycle with lights on at 7 AM and lights off at 7 PM Central Standard Time.22 Animals were individually housed in ventilated, light-tight cabinets with 40W GE “Cool White” fluorescent bulbs individually controlled by an electronic timer.43
Food timing paradigm
Mice were placed on a food-restricted paradigm. As mice are nocturnal and consume their food in the dark, food timing in this experiment is defined with either day- (rest-phase= “wrong-time”, 7 AM) or nighttime (active-phase= “right-time”, 7 PM). Mice received access to food for 12 h for 5 consecutive weekdays followed by a 2-day ad-libitum access to food to mimic human conditions on the weekend.9 Mice were fed Envigo 2018 standard rodent chow (Teklad, Madison, WI) and given either normal drinking water or alcohol in drinking water (described below). This paradigm lasted a total duration of 10 weeks (Supplemental Figure S7). The general condition and health of the mice were monitored by daily observation, weekly body weight, food intake, and liquid intake measurements.47
Alcohol treatment
The control groups received drinking water (H₂O). The mice in alcohol-fed groups were given 20% alcohol v/v treatment in drinking water.48 The 20% alcohol was only provided during the mice’s food intake periods, so it was available ad-libitum during the weekends when the mice had unrestricted access to both food and liquid.49 Rodents reportedly metabolize alcohol much more rapidly than humans, at around 5.5 times the human rate.50 To achieve blood alcohol levels (BAL) in mice that are pharmacologically relevant and comparable to BAL in alcohol-consuming human (~1.0 g/L), a 20% alcohol dose in drinking water was required.51 To acclimate the mice to alcohol consumption, a 14-day ramp-up was used, starting at 3% on days 1–3, increasing to 5% on days 4–7, 10% on days 8–10, 15% on days 11–14, and then reaching the target treatment of 20% on day 15. The volume was replaced as needed, with a maximum fill of 30 mL, and daily intake of both water and alcohol was measured. After 10 weeks, spontaneously voided stool samples were collected with sterile forceps to assess the intestinal microbial communities, followed by a urine test for intestinal barrier integrity. Finally, at the end of the study, the mice were euthanized via guillotine. Then, trunk blood and tissue samples were collected (Supplemental Figure S8).
Intestinal barrier (permeability) testing
The mice were fasted overnight prior to the urine permeability test, which is performed the next morning.22,30 A 200 µl solution containing lactulose (3.2 mg), sucrose (0.45 mg), sucralose (0.45 mg), and mannitol (0.9 mg) is administered via oral gavage, followed by a subcutaneous injection of 2 mL of 0.9% saline to promote urine production. The mice were then placed into a metabolic chamber with water available ad-libitum for 5 h, after which the urine produced was collected and the total volume was recorded. Urine samples were stored at −80°C until further experimental analysis. This 5-h timeframe allows for intestinal transit, enabling assessment of barrier integrity across the entire intestine, including the colon. Intestinal barrier integrity is determined by measuring urinary concentrations of the administered sugars using gas chromatography, which enables calculation of the percentage of the oral dose that was exerted in the urine over 5 h. Then, intestinal permeability is assessed at the end of the 10-week study period.30
Intestinal microbiota analysis
Stool collection
To collect stool samples for microbiota analysis, individual mice were housed in a bedding-free cage for 12 h. After this 12-h period, the spontaneously produced stool pellets were gathered using sanitized forceps and stored at −80°C until further experimental analysis.
DNA extraction
The total genomic DNA was extracted from end-of-study stool pellets using the FastDNA SPIN Kit. The extracted DNA was then quantified using fluorometric quantitation on a Qubit 3.0 instrument (Life Technologies, Grand Island, NY, USA). To decrease run effects, all samples were extracted together using the same DNA extraction kit, and library preparation was conducted concurrently in 96-well plates.
The variable region 4 (V4) of microbial 16S ribosomal RNA (rRNA) genes was targeted for PCR amplification using the modified Earth Microbiome Project primers 515F/806 R (515F: GTGTGYCAGCMGCCGCGGTAA; 806 R: CCGGACTACNVGGGTWTCTAAT). The PCR method, as previously described done in a two-stage fashion, was used to prepare the amplicons for high-throughput sequencing.52,53 The amplicon libraries were sequenced on an Illumina MiniSeq, utilizing a V2 kit and generating paired-end (150) base reads, at Rush University Medical Center utilizing the Rush Genomics and Microbiome Core Facility.
16S rRNA V4 sequencing analysis
PEAR (paired-end read merger) (v0.9.11) software package (Dalhousie University, Halifax, Nova Scotia, Canada) was used to analyze merged raw sequences.54 A DADA2 algorithm from the QIIME2 (v 2020.8.0) workflow was utilized to remove sequences shorter than 240 bases, with denoising, were then removed before the merged were processed.55,56 The resulting amplicon sequence variants (ASVs) were used for all subsequent analyses. Each ASV was allotted a taxonomy using an unbiased Bayes classifier and the SILVA 138 99% database reference for OTUs.57,58 In total, 4,229,677 sequencing clusters were generated, with an average depth of 44,059 sequences per sample (median = 41,906; min = 1; max = 113,563). Three reagent contaminant ASVs (Rhodococcus, Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium, and Pseudomonas) were identified and removed using the decontam algorithm based on their prevalence in the reagent negative blank controls (n = 4).59 Unassigned and host-associated taxa, such as mitochondrial ASVs, chloroplast, and eukaryote were also removed from the datasets prior to statistical analyses.60 After analysis, the raw sequence data were deposited in the NCBI Sequence Read Archive under BioProject PRJNA881333.
Blood collection and tissue harvest
At the end of the study, the mice were euthanized via guillotine. Blood was collected by trunk collection, then centrifuged at 2,000 RPM for 35 min at 22°C. The serum was collected and stored at −80°C until experimental analysis. Murine tissue samples were collected: taking sections of the liver and the large intestine. The large intestine was flushed clean with cold phosphate-buffered saline (PBS) at least three times to remove excess debris. The tissue segments were snap-frozen and stored at −80°C, or formalin-fixed and paraffin-embedded (FFPE), as previously described.42,61
Per2-luciferase bioluminescence
Rhythmic expression of the clock gene PER2 (PER2 luciferase, PER2:LUC) persists in these colonic tissues and organoids, allowing constant monitoring of the robustness of these oscillations.40 Freshly collected mouse colonic tissue was immediately transferred onto Matrigel-coated 35 mm dishes to minimize tissue movement in the media-filled dishes. Floating colonic organoids were also assessed in these dishes separately. The Kronos Dio apparatus (AB-2500, Atto Co., Tokyo, Japan) quantified bioluminescence levels of the PER-2 luciferase protein (E1602, Promega) in luciferin-containing media. Measurements were taken for 1 min every 15 min, at a controlled temperature of 37°C and 5% CO2, for a minimum of 5 days.
Mouse organoid preparation and culture
The mouse intestinal epithelial organoid cultures from the colon were prepared using a previously described method.42,62–64 Briefly, the colon was washed with ice-cold PBS lacking Ca++ and Mg++ (Thermo Fisher Scientific, #100-10-023) to remove any loose debris. The samples were cut longitudinally, and excess fat was gently scraped off. After two PBS washes, the samples were cut into small pieces and transferred into a 50 mL conical tube. The tissue was then washed five times with PBS, transferred to a separate 50 mL conical with 2.5 mm EDTA chelating buffer, and incubated on a shaker at 4°C for 1 h. The supernatant was removed, and the tissue pieces were washed in PBS, discarding the first fraction. Three mL of PBS was added, and the crypts were gently resuspended and filtered through a 70 μm cell strainer into a fresh conical tube. This washing and filtering was repeated three times. The crypt suspension was centrifuged at 300 g for 5 min after adding 10% FBS. The supernatant was discarded, and the cell pellet was resuspended in 15 mL of organoid media minus growth factors. Base organoid media consisted of: Advanced DMEM/F-12 (Thermo Fisher Scientific), Glutamax (Thermo Fisher Scientific), Penicillin/Streptomycin (Sigma Aldrich), HEPES (Sigma Aldrich), B27-supplement (Thermo Fisher Scientific), N2-supplement (Thermo Fisher Scientific), and N-Acetylcysteine (Sigma Aldrich).
The cell suspension was then centrifuged at 150 g for 3 min. Then, this washing/centrifugation step was repeated twice. After the final centrifugation, the crypts were counted, resuspended in cold Matrigel (60 μL per well), and plated on a warm 24-well culture plate. Then, the plate was incubated at 37°C for 30 min to allow the Matrigel to polymerize before adding 500 μL of organoid medium with the growth factors to each well. Growth factors for colonic organoids included: recombinant mouse EGF (Sigma Aldrich), recombinant mouse Noggin (Sigma Aldrich), recombinant human R spondin 1 (PeproTech), recombinant mouse WNT-3A (PeproTech), A-83-01 (Tocris Bioscience), SB202190 (Sigma Aldrich), Y-27632 (Sigma Aldrich), and nicotinamide (Sigma Aldrich). The organoids were incubated at 37°C/5% CO2. Organoids were observed daily to ensure proper organoid health and robustness. Then, the organoids were passaged at least three times (roughly 3 weeks from the initial preparation) before any experimentation was performed.
Apical out intestinal organoids
We generated apical out organoids from these initial organoids using a previously described method.65 Organoid growth was assessed under the microscope daily. To ensure proper health and robustness, the organoids were passaged at least three times (approximately 3 weeks from the initial preparation) before any experimentation was performed.
Cecal content supernatant
The cecal content supernatant sample preparations were adapted from previously published methods.66 The cecal contents were pooled together (n = 4) for each group, weighed on ice, and 500 μL of cold PBS was added 100 mg cecal content. The sample tubes were then vortexed for 5 min and centrifuged at 13,000 g for 15 min at 4°C. The resulting supernatant was filtered through a 0.22 μM membrane with a 1 ml syringe before the sample was stored at −80°C until experimental use.
Sample preparation
The cecal content supernatant samples were prepared with the automated MicroLab STAR® system from Hamilton Company from previously described methods from Metabolon.67 Several recovery standards were added prior to the first extraction step in for quality control purposes. To extract chemically diverse metabolites while removing proteins, the cecal content supernatant samples were precipitated with methanol under vigorous shaking for 2 min, using a Glen Mills GenoGrinder 2000, followed by centrifugation. This helped dissociate small molecules bound to or trapped in the precipitated protein matrix. The resulting sample extract was divided into five fractions: two for analysis by positive ion mode electrospray ionization (ESI) reverse phase (RP)/ultra-performance liquid chromatography (UPLC)-tandem mass spectrometry (MS/MS), one for negative ion mode ESI RP/UPLC-MS/MS, one for negative ion mode ESI hydrophilic interaction liquid chromatography (HILIC)/UPLC-MS/MS, and one sample was reserved as a backup. Finally, organic solvent was removed from the samples then stored overnight under nitrogen before preparation for analysis.
Ultrahigh performance liquid chromatography-tandem mass spectroscopy (UPLC-MS/MS)
The analysis utilized a Waters ACQUITY UPLC system coupled with a Thermo Scientific Q-Exactive high resolution/accurate mass spectrometer using previously described methods.68 The mass spectrometer was equipped with a heated electrospray ionization (HESI-II) source and Orbitrap mass analyzer, operated at 35,000 mass resolution. The cecum content supernatant sample extract was dried and reconstituted in solvents compatible with analytical standards at fixed concentrations to ensure consistent injection and chromatographic. Four different ionization methods were used to analyze the samples. One aliquot was analyzed using acidic positive ion conditions using a C18 column (Waters UPLC BEH C18–2.1×100 mm, 1.7 µm) and a gradient of water and methanol, containing 0.1% formic acid (FA) and 0.05% perfluoropentanoic acid (PFPA). This method was optimized for more hydrophilic compounds. The second aliquot was also analyzed using acidic positive ion conditions using the same C18 column, but with a gradient of acetonitrile, methanol, water, 0.05% PFPA, and 0.01% FA. This method was optimized for more hydrophobic compounds and operated at a higher overall organic content. A third aliquot was analyzed under basic negative ion optimized conditions using a separate dedicated C18 column. Utilizing methanol and water with 6.5 mm ammonium bicarbonate at pH 8 the basic extracts were gradient eluted. The final, fourth aliquot was determined by negative ionization by elution with a HILIC column (Waters UPLC BEH Amide 2.1 × 150 mm, 1.7 µm). About 10 mm ammonium formate at pH 10.8 with a gradient of water and acetonitrile was utilized. A scan range across 70–1000 m/z alternated between MS and data-reliant MS scans using dynamic prohibition was performed.
Fluorescein isothiocyanate (FITC)-dextran
To assess epithelial barrier integrity, apical-out organoids were treated for 3 h with either control media, 0.2% EtOH, or a 1:100 dilution of cecal content supernatant (RTE H₂O or WTE EtOH) with control media. The organoids were then pelleted, and washed in PBS (4°C, 5 min, 5,000 RPM) before being incubated for 30 min in a solution (2 mg/ml) of 4 kDa fluorescein isothiocyanate (FITC)-dextran (Sigma, FD4-1 G) in a 37° incubator.65 Organoids are pelleted and washed (4°C, 5 min, 5,000 RPM) in a live cell imaging solution (Thermo Fisher, #A14291DJ). Then, the organoids were placed in 8-well chamber slides (Nunc, Rochester, NY) and subsequently imaged using a Zeiss LSM 700 confocal microscope. All images were taken at × 20 magnification with 20–25 enteroids per group. Then, all experiments were performed in triplicate. Data are reported as mean ± standard error of the mean (SEM).
Immunofluorescent staining of colon tissue and organoids
Colon tissue
The colon tissue immunofluorescence (IF) staining was performed on formalin-fixed and paraffin-embedded mouse tissues cut into 5 μm sections. After deparaffinization and rehydration antigen retrieval was carried out using 1X Dako antigen retrieval solution (Agilent # S1699). The slides were then stained with adherens and tight junction protein antibodies for Occludin (mouse, Invitrogen #33–1500), ZO-1 (rabbit, Invitrogen #40–2200), and E-cadherin (mouse, Cell Signaling #14472S). Tissue staining data was collected from n = 9–10 mice per group, with images taken from at least 10 stained tissues. The relative expression of each marker was quantified using ImageJ software (ImageJ, U. S. National Institutes of Health, Bethesda, MD) and representative images were selected. All staining was evaluated by two blinded independent observers using a Zeiss LSM 700 confocal microscope. All images were taken at × 40 magnification. The immunofluorescence of each antibody was quantified by calculating the mean difference between background and overall fluorescence of each region. The fluorescence values for each image and section were then averaged to determine the immunofluorescence value for each mouse, using GraphPad Prism (v9.1) software (GraphPad Software, La Jolla, CA).
Organoids
The organoid immunofluorescent staining was completed on 8-well chamber slides. The staining occurred after three weeks or longer after colon isolation and organoid generation, allowing optimal colonic organoid growth. The organoids in Matrigel were washed three times with cold PBS after the media was removed. Then, 500 μL of Corning Cell Recovery Solution was added to each well, and the Matrigel suspension was moved to a 15 mL conical tube. An additional 500 μL of Cell Recovery Solution was then added to wash the well. To fully dissolve the Matrigel, the samples were kept on ice for 1 h, with regular inversion. The cells were then spun at 300 g for 5 min. The cells were washed twice more with ice-cold PBS after the supernatant was removed.
The organoids were fixed with 4% paraformaldehyde (PFA) on chamber slides for 1 h at room temperature, and then permeabilized with 1% Triton-X100. The samples were blocked with a solution of 1% bovine serum albumin (BSA), 3% goat serum, and 0.2% Triton X-100 for 1 h. The primary staining for the organoids was performed for 24 h at room temperature in a humidified chamber. On the second day, the slides were washed and incubated for 2 h with their relevant secondary antibodies conjugated to Alexa Fluor 488. The samples were further stained for DAPI and mounted using Fluoromount aqueous mounting medium (Sigma Aldrich), after washing.
A Zeiss LSM 700 confocal microscope was used to image the tissue and organoid slides. The specific antibodies for staining are detailed in Supplemental Table S2. All organoid staining data are from n = 1 mice with images from at least 10 stained organoids. Section areas were used to determine the relative expression of each marker using ImageJ software and to pick the representative images for analysis. All staining was evaluated by two blinded independent observers using a Zeiss LSM 700 confocal microscope. All images were taken at × 40 magnification. The immunofluorescence of each antibody was quantified by calculating the mean difference between the background and overall fluorescence of each region. The fluorescence values for each image and section were then averaged to determine the immunofluorescence value for each mouse, using GraphPad Prism (v9.1) software.
Western blot analysis
Tissue samples were immediately frozen in liquid nitrogen and stored at −80°C until used for protein extraction and concentration prior to western blot analysis. To ensure organoid growth and robustness for experimentation, the organoid protein was collected at least 21 days (approximately 3 weeks) after colon isolation and organoid generation. The commercially available kit, NE-PER Nuclear and Cytoplasmic Extraction Reagents (Thermo Scientific 78,833), was used to obtain cellular material from both colon tissue and organoid samples.
2-mercaptoethanol (Bio-Rad) in a laemmli sample buffer with was used to prep the samples. For 2 h each gel was electrophoresed at 100 V with 20 μg of protein per lane. A 4%/7.5% stacking acrylamide Tris gel was used. As previously described, the membrane transference, blocking, and antibody incubation of both primary and secondary were performed.20,21 Autoradiography film (HyBlot CL, Denville Scientific, Metuchen, NJ) with chemiluminescent substrate (ECL, GE Healthcare) was then placed on the membrane, and protein imaging was attained. Image J software was used to examine optical density by densitometric analysis. The specific antibodies are detailed in Supplemental Table S3 (below).
Hematoxylin and eosin staining
To assess histologic changes within the colonic tissue, 4- to 5-µm thick sections of fixed colon tissue samples were stained with hematoxylin and eosin, as previously described.61,69 A gastrointestinal pathologist blinded to the sample identities, scored the tissues across seven categories to measure colon inflammation, tissue injury, and repair: inflammatory cells (0–4), goblet cells (0–4), mucosa thickening (0–4), submucosa cell infiltration (0–4), destruction of architecture (0–4), ulcers (0–1), and crypt abscess (0–1).61
Lipopolysaccharide-binding protein assay
Lipopolysaccharide (LPS) is a component found in the outer membrane of Gram-negative bacteria. LPS-binding protein (LBP) is a type 1 acute-phase protein that binds to LPS, facilitating an immune response. LBP is a well-accepted marker of intestinal barrier integrity and endotoxemia. Serum collected at the time of euthanasia was stored at −80°C until experimental analysis to measure LBP levels using an LBP ELISA kit (Hycult Biotech, HK205), following the manufacturer’s instructions.
Serum interleukin-6 analysis
Serum levels of the cytokine interleukin-6 (IL-6) were measured using a Meso Scale V-PLEX Pro-inflammatory Mouse Kit (Cat. #K152QXG–1, Meso Scale Diagnostics, Rockville, MD). All samples were analyzed in duplicate on a QuickPlex SQ120 instrument (Meso Scale Diagnostics), following the manufacturer’s protocol.
Stool calprotectin
At the end of the study, mouse fecal samples were collected and stored at −80°C until experimental analysis. The fecal samples were subjected to only a single freeze-thaw cycle to maintain stool integrity. The calprotectin levels were measured using an ELISA kit (Cat. #KR6936, Eagle Biosciences, Amherst, NH).
Blood alcohol levels
Blood samples (25 µl each) were removed from storage at −80°C and thawed at 4°C for 15 min. To each sample, a standard (0.1%, or 0.804 mg/ml of n-propanol [nPrOH]; 5 μL) was then added to each sample. Each sample was then transferred and mixed to a 2 mL GC vial by pipette, which was capped right away. The sealed samples were placed on the GC autosampler and allowed to equilibrate at room temperature (22–23°C) for a minimum of 15 min prior to injection. The 25 μl headspace vapor from each sample was injected into a Trace 1310 GC coupled to a Thermo ISQ-LT MS, with a split ratio of 10:1 using a gastight 100 μl syringe. The GC inlet was held at 250°C. A 30 m DB-WAX UI column (0.25 mm ID, and 0.25 μm film thickness) was used to achieve optimal separation. The oven temperature was kept at 40°C for 4 min, and then ramped at 40°C/min to 120°C and held for 1 min. Helium carrier gas flow was maintained at 1.2 mL/min. 250°C was used for both transfer line and ion source. The SIM mode was used to scan ions, at m/z 31 and 45 (EtOH), m/z 31 and 42 (nPrOH), and m/z 59 and 43 (iPrOH [IPA]).
Statistical analysis
Tissue and organoid analyses
Data are reported as mean ± standard error of the mean (SEM). A two-way analysis of variance (ANOVA) was used to evaluate the main effects of food timing (RTE vs. WTE), treatment (H₂O vs. EtOH), and their interaction. Planned a priori comparisons between groups were conducted using a Tukey test to control the type I error.70 A separate two-way ANOVA assessed sex examining the main effects of treatment (H₂O vs. EtOH), food timing (RTE vs. WTE), and their interaction. Statistical significance was set at p < 0.05. All analyses were performed using GraphPad Prism (v9.1) software.
Spearman’s correlation was used to evaluate the relationships between the relative abundances of specific taxa (species) and phenotype, identifying potential targets for future investigation. A significant threshold was set at p < 0.05 and R > 0.30.71 This approach has been applied in prior studies to uncover associations between bacteria and alcohol consumption.
Microbiota analysis
Analyses of alpha- and beta-diversity were used to compare the stool microbial community structure between groups (e.g., EtOH vs. H₂O; WTE vs. RTE) with feature (ASV) counts. Alpha-diversity metrics (Shannon Index, Simpson’s Index, Observed Features (number of taxa), and Pielou’s Evenness (relative abundance of those taxa) were calculated on rarefied datasets of 25,000 sequences/sample, with significance level set at p < 0.05. These analyses were conducted using GraphPad Prism (v9.1) software.
Then, Permutation Multivariate Analysis of Variance (PERMANOVA) and Permutational Analysis of Multivariate Dispersions (PERMDISP) were used to compare microbial community structure, with significance determined using 9,999 permutations and corrected using the Benjamini–Hochberg method (q < 0.05).72,73 Principal coordinates analysis (PCoA) based on a Bray-Curtis dissimilarity distance matrix was used to visualize baseline microbial community structure of group comparisons (RTE vs. WTE and H₂O vs. EtOH) within the software package QIIME2.56 Non-metric multidimensional scaling (NMDS) plots were also generated to visualize bacterial species community data for these group comparisons (RTE H₂O, RTE EtOH, WTE H₂O, and WTE EtOH). Each sample was connected to a centroid representing the mean value of the group.
Boruta feature selection was performed in R.74 Both Centered Log-Ratio Kruskal Wallis (CLR-KW) and DESeq2 testing were performed to identify significantly differentially abundant features between mice groups, with Benjamini–Hochberg method correction (q < 0.05).75–77 Features with an average relative abundance below 0.1% were removed from the analysis. Finally, individual taxa percent mean relative abundances (˃1%) ± standard deviations (SD) were calculated and depicted as stacked histograms.
Serum gas chromatography–Mass spectrometry analysis
The GC/MS data was analyzed using Chromeleon software. This allowed for visual inspection of the retention time and peak area integration for each target analyte (EtOH, nPrOH, and iPrOH). The peak areas for EtOH and nPrOH were then extracted from the data for each sample. Absolute quantitation (in mg/mL and percentage) was calculated using the linear regression equation generated from the calibration curve for each compound. Then, peaks were quantified using the area-under-the-curve method.
Cecal content pathway analysis
The sample data was analyzed with Reactome software.78 Mouse data identifiers were converted to their human equivalents before analysis was performed. Reactome curates human pathways and infers their existence in other species using orthology information from the ENSEMBL Compara database.33
Circadian rhythmicity
Analysis of variance (ANOVA) was used to detect the differences in circadian period between right-time and wrong-time mouse groups. All analyses set significance levels at p < 0.05. The statistics were performed using GraphPad Prism (v9.1), SPSS (v26) (IBM Incorporation), and circular data was measured by Oriana (v4). To analyze the luciferase time-series data, Biodare 2 was used. Periodogram was performed using the Fast Fourier Transformation Nonlinear Least Square algorithm (FFT-NLLS) to measure period, acrophase, and amplitude.79 FFT-NLLS is a variation of the cosinor method and is less sensitive to missing or noisy data. The cosinor fit equation is: Yt=M+Acos2πt/τ+ϕ where M is the MESOR (Midline Statistic of Rhythm, a rhythm-adjusted mean), A is the amplitude (a measure of half the extent of the variation within the cycle), ϕ is the acrophase (a measure of the time of overall highest value), and τ is the period.80,81 The circular variance is calculated as [V = 1-r], where r = length of the mean vector, and V = Circular variance. Circular variance has been found to be an important marker in genetic models of circadian disruption.82 All circular statistics were calculated using Oriana (v4).