Gut microbial metabolism of dietary flavonoids leads to a diverse array of bioactive products that are closely associated with human health. Combining enzyme promiscuity prediction, metabolomics, and in vitro model systems, we identified a chalcone-synthase-like bacterial polyketide synthase that can initiate the metabolism of naringenin by catalyzing the C-ring cleavage. This was validated using a mutant strain of the model organism Bacillus subtilis (ATCC 23857). Our prediction–validation methodology could be used to systematically characterize the products of gut bacterial flavonoid metabolism and identify the responsible enzymes and species. In vitro experiments with Caco-2 cells revealed that naringenin and its bacterial metabolites differentially engage the aryl hydrocarbon receptor (AhR) and orphan nuclear receptor 4A (NR4A). These results suggest that metabolism by gut bacterial species could directly impact the profile of bioactive flavonoids and influence inflammatory responses in the intestine. These results are significant for understanding gut-microbiota-dependent physiological effects of dietary flavonoids.

BACKGROUND:Concurrent consumption of dietary fiber and n-3 polyunsaturated fatty acids reduces colon tumor formation. However, their combined effects on colorectal cancer risk remain unexplored in human trials. OBJECTIVES:This study investigated the synergistic effects of fish oil (FO) and fermentable fiber on the gut transcriptional profiles and microbiome composition in older adults. METHODS:In a randomized controlled crossover pilot study, 30 adults (ages 50-75 y), received fermentable fiber (33 g/d soluble corn fiber; SCF) plus eicosapentaenoic acid and docosahexaenoic acid (EPA + DHA as FO, 7.7 g/d) or a comparator (similar doses of maltodextrin plus corn oil; MD + CO) for 30 d, followed by a 60-d washout period before crossing over to the alternate intervention. Serum phospholipid fatty acids, stool exfoliome [ribonucleic acid sequencing (RNAseq)], microbiome (16S ribosomal ribonucleic acid gene sequencing), butyrate kinase (but) gene abundance (digital droplet polymerase chain reaction), and fecal short-chain fatty acids were analyzed. Linear mixed models were used for the majority of outcome analyses. Differential expression and pathway enrichment analyses were applied to RNAseq data, whereas microbiome diversity was assessed using α and β diversity. RESULTS:Serum EPA and DHA concentrations were higher after SCF + FO than MD + CO supplementation [EPA: β = 0.51; 95% confidence interval: 0.31, 0.72; DHA: β = 0.18; 95% confidence interval: 0.10, 0.27; P < 0.0001]. Analysis of host gut transcriptional networks revealed that SCF + FO supplementation inhibited the glucose-insulin receptor-phosphatidylinositol-3 kinase-signaling axis. Microbiome analysis revealed significant intervention differences in β-diversity (F = 4.4, R2 = 0.08, P = 0.001), and 27 of 73 genera analyzed, several known short-chain fatty acid producers, differed between the 2 interventions (false discovery rate <0.05). Abundance of the but gene from Roseburia sp (P < 0.001) and the genera Roseburia (P = 0.006) were lower in the SCF + FO compared to MD + CO intervention, although fecal butyrate concentrations did not differ. CONCLUSIONS:Thirty-day supplementation of SCF + FO compared with MD + CO showed significant shifts in intestinal cell pathways relevant to colorectal cancer with concomitant differences in gut microbial community structure and butyrate-producing taxa.

An effective therapeutic strategy to treat oncogenic Wnt signaling in the context of colorectal cancer (CRC) remains elusive. A new study from Cho and colleagues describes a novel mechanistic link between the loss of canonical adenomatous polyposis coli (APC) function, membrane cholesterol, and an innovative drug target to specifically suppress the cholesterol-Dvl-β-catenin signaling axis.

Single-cell RNA sequencing (scRNA-seq) has revolutionized our understanding of cellular variability by capturing gene expression profiles of individual cells. The importance of cell-to-cell variability in determining and shaping cell function has been widely appreciated. Nevertheless, differential expression (DE) analysis remains a cornerstone method in analytical practice. Current computational analyses overlook the rich information encoded by variability within the single-cell gene expression data by focusing exclusively on mean expression. To offer a deeper understanding of cellular systems, there is a need for approaches to assess data variability rather than just the mean. Here we present spline-DV, a statistical framework for differential variability (DV) analysis using scRNA-seq data. The spline-DV method identifies genes exhibiting significantly increased or decreased expression variability among cells derived from two experimental conditions. Case studies show that DV genes identified using spline-DV are representative and functionally relevant to tested cellular conditions, including obesity, fibrosis, and cancer.

Polyfluoroalkyl substances (PFAS) are widely used industrial compounds that have been identified as contaminants in almost every component of the global ecosystem, and in human studies, higher levels of PFAS have been correlated with increased incidence of multiple diseases. Based on the results of human and laboratory animal studies, we hypothesize that the orphan nuclear receptor 4A1 (NR4A1) may be a critical target for some PFAS such as the legacy linear polyfluorooctanesulfonate (PFOS) and other sulfonates. We show that PFOS and related compounds bound the ligand binding domain (LBD) of NR4A1 and induced the growth of several cancer cell lines and enhanced tumor growth in an athymic nude mouse model. Using NR4A1-responsive rhabdomyosarcoma Rh30 cells as a model, PFOS induced NR4A1-dependent cell proliferation and Rh30 cell migration and invasion. Moreover, in Rh30 cells, PFOS also induces several NR4A1-regulated genes including the PAX3-FOXO1 oncogene and downstream gene products, and in a chromatin immunoprecipitation assay, PFOS does not decrease NR4A1 binding to the promoter. These results demonstrate that PFOS is an NR4A1 ligand and enhances tumorigenesis through the activation of this receptor.

BACKGROUND:Several studies have indicated the benefits of human milk feeding to infants however, mechanisms behind positive health outcomes have not been investigated. OBJECTIVES:The study aimed to characterize circulating immune cell subpopulation gene expression in human milk-fed (HMF) compared with cow milk formula-fed (FF) infants using single-cell transcriptomics. METHODS:Peripheral blood mononuclear cells (PBMCs) were isolated from healthy HMF (n = 6), and FF (n = 3) infants who were 3-3.5 mo old and enrolled in a nonrandomized clinical trial. Single-cell RNA sequencing was used to generate a PBMC atlas and evaluate gene expression in immune cell subsets. Differential expression analysis was performed on each cell type independently after clustering the cells by similar marker gene expression using the scGEAToolbox. Differentially expressed genes were subjected to pathway analyses using an online functional enrichment analysis program. RESULTS:The relative abundance (%) of T and B lymphocytes, natural killer (NK) cells, and plasmacytoid dendritic cells were similar, whereas monocytes were higher in FF infants than in HMF infants (22.6 ± 10.7 compared with 8.3 ± 5.6; P = 0.0314). In addition, innate and adaptive immune cells from FF infants exhibited a higher activation state compared with HMF infants. We identified 16 distinct cell subsets from the major immune cell types: 3 monocyte subsets, 4 NK subsets, 2 B cell subsets, and 7 T cell subsets. Transcriptional profiles of each peripheral innate and adaptive immune cell subtype varied between HMF and FF infants. Pathway enrichment analysis of cell-specific transcriptional changes within subsets of major cell types revealed that the interleukin (IL)-4/IL-13 signaling pathways were upregulated in FF infants relative to HMF infants. CONCLUSIONS:These findings suggest that human milk downregulates peripheral immune cell cytokine transcriptional signatures linked to allergic inflammation and infection relative to formula feeding.

Single-cell RNA sequencing (scRNA-seq) data simulation is limited by classical methods that rely on linear correlations, failing to capture the intrinsic, nonlinear dependencies. No existing simulator jointly models gene-gene and cell-cell interactions. We introduce qSimCells, a novel quantum computing-based simulator that employs entanglement to model intra- and inter-cellular interactions, generating realistic single-cell transcriptomes with cellular heterogeneity. The core innovation is a quantum kernel that uses a parameterized quantum circuit with CNOT gates to encode complex, nonlinear gene regulatory network (GRN) as well as cell-cell communication topologies with explicit causal directionality. The resulting synthetic data exhibits non-classical dependencies: standard correlation-based analyses (Pearson and Spearman) fail to recover the programmed causal pathways and instead report spurious associations driven by high baseline gene-expression probabilities. Furthermore, applying cell-cell communication detection to the simulated data validates the true mechanistic links, revealing a robust, up to 75-fold relative increase in inferred communication probability only when quantum entanglement is active. These results demonstrate that the quantum kernel is essential for producing high-fidelity ground-truth datasets and highlight the need for advanced inference techniques to capture the complex, non-classical dependencies inherent in gene regulation.

OBJECTIVE:Fatty Acid Desaturase 1 (FADS1) is a rate-limiting enzyme controlling the bioproduction of long-chain polyunsaturated fatty acids (PUFAs). Increasing studies suggest that FADS1 is a potential cancer target. Our previous research has demonstrated the significant role of FADS1 in cancer biology and patient survival, especially in kidney cancers. We aim to explore the underlying mechanism in this study. METHOD AND RESULTS:We found that pharmacological inhibition or knockdown of the expression of FADS1 significantly reduced the intracellular conversion of long-chain PUFAs, effectively inhibits renal cancer cell proliferation, and induces cell cycle arrest. The stable knockdown of FADS1 also significantly inhibits tumor formation in vivo. Mechanistically, we showed that while FADS1 inhibition induces endoplasmic reticulum (ER) stress, FADS1 expression is augmented by ER-stress inducer, suggesting a necessary role of PUFA production in response to ER stress. FADS1-inhibition sensitized cellular response to ER stress inducers, leading to cell apoptosis. Also, FADS1 inhibition-induced ER stress leads to activation of the PERK/eIF2α/ATF4/ATF3 pathway. Inhibiting PERK or knockdown of ATF3 rescued FADS1 inhibition-induced ER stress and cell growth suppression, while ATF3-overexpression aggravates the FADS1 inhibition-induced cell growth suppression and leads to cell death. Metabolomic analysis revealed that FADS1 inhibition results in decreased level of UPD-N-Acetylglucosamine, a critical mediator of the unfolded protein response, as well as impaired biosynthesis of nucleotides, possibly accounting for the cell cycle arrest. CONCLUSION:Our findings suggest that PUFA desaturation is crucial for rescuing cancer cells from persistent ER stress, supporting FADS1 as a new therapeutic target.

Feature selection is a machine learning technique for identifying relevant variables in classification and regression models. In single-cell RNA sequencing (scRNA-seq) data analysis, feature selection is used to identify relevant genes that are crucial for understanding cellular processes. Traditional feature selection methods often struggle with the complexity of scRNA-seq data and suffer from interpretation difficulties. Quantum annealing presents a promising alternative approach. In this study, we implement quantum annealing-empowered quadratic unconstrained binary optimization (QUBO) for feature selection in scRNA-seq data. Using data from a human cell differentiation system and an anticancer drug resistance study, we demonstrate that QUBO feature selection effectively identifies genes whose expression patterns reflect critical cell state transitions associated with differentiation and drug resistance development. Our findings indicate that quantum annealing-powered QUBO reveals complex gene expression patterns potentially missed by traditional methods, thereby enhancing scRNA-seq data analysis and interpretation.

Bis-indole-derived compounds including 1,1-bis(3'-indolyl)-1-(3,5-disubstitutedphenyl)methane (DIM-3,5) analogs bind both orphan nuclear receptors 4A1 (NR4A1) and NR4A2, and DIM-3,5 compounds act as dual receptor inverse agonists and inhibit both NR4A1- and NR4A2-regulated responses. Chromatin immunoprecipitation assays show that β1-integrin and the methyltransferase gene G9a are regulated by both NR4A1 and NR4A2 acting as cofactors for Sp1- and Sp4-dependent gene expression. DIM-3,5 treatment results in the loss of one or more of these nuclear factors from the β1-integrin and G9a promoters. Single-cell and RNAseq analyses show that both receptors regulate common (<10%) and unique genes in SW480 colon cancer cells; however, functional enrichment analysis of the differentially expressed genes converges to several common pathways and gene ontology terms.

Parkinson’s disease (PD), characterized by death of dopaminergic neurons in the substantia nigra, is the second most prevalent progressive neurodegenerative disease. However, the etiology of PD is largely elusive. This study employed the 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) rodent model to examine the effectiveness of 1,4-dihydroxy-2-naphthoic acid (1,4-DHNA), an aryl hydrocarbon receptor (AhR) active gut bacteria-derived metabolite, in mitigating MPTP’s motoric deficits, and the role of AhR in mediating these effects. Male C57BL/6 mice were fed daily with vehicle, 20 mg/kg 1,4-DHNA, or AhR-inactive isomer 3,7-DHNA, for 3 weeks before administration of 80 mg/kg MPTP or vehicle. Four weeks later, mice were assessed for motoric functions. Both 1,4-DHNA and 3,7-DHNA prevented MPTP-induced deficits in the motor pole test and in the adhesive strip removal test. Additionally, 1,4-DHNA improved balance beam performance and completely prevented MPTP-induced reduction in stride length. In contrast, 3,7-DHNA, an AhR-inactive compound, did not improve balance beam performance and had only a partial effect on stride length. This study suggests that natural metabolites of gut microbiota, such as 1,4-DHNA, could be beneficial to counteract the development of motor deficits observed in PD. Thus, this study further supports the hypothesis that pathological and mitigating processes in the gut could play an essential role in PD development. Moreover, this indicates that 1,4-DHNA’s ability to combat various motor deficits is likely mediated via multiple underlying molecular mechanisms. Specifically, AhR is involved, at least partially, in control of gait and bradykinesia, but it likely does not mediate the effects on fine motor skills.

Monounsaturated fatty acids (MUFAs) play a pivotal role in maintaining endoplasmic reticulum (ER) homeostasis, an emerging hallmark of cancer. However, the role of polyunsaturated fatty acid (PUFAs) desaturation in persistent ER stress driven by oncogenic abnormalities remains elusive. Fatty Acid Desaturase 1 (FADS1) is a rate-limiting enzyme controlling the bioproduction of long-chain PUFAs. Our previous research has demonstrated the significant role of FADS1 in cancer survival, especially in kidney cancers. We explored the underlying mechanism in this study. We found that pharmacological inhibition or knockdown of the expression of FADS1 effectively inhibits renal cancer cell proliferation and induces cell cycle arrest. The stable knockdown of FADS1 also significantly inhibits tumor formation in vivo . Mechanistically, we show that while FADS1 inhibition induces ER stress, its expression is also augmented by ER-stress inducers. Notably, FADS1-inhibition sensitized cellular response to ER stress inducers, providing evidence of FADS1’s role in modulating the ER stress response in cancer cells. We show that, while FADS1 inhibition-induced ER stress leads to activation of ATF3, ATF3-knockdown rescues the FADS1 inhibition-induced ER stress and cell growth suppression. In addition, FADS1 inhibition results in the impaired biosynthesis of nucleotides and decreases the level of UPD-N-Acetylglucosamine, a critical mediator of the unfolded protein response. Our findings suggest that PUFA desaturation is crucial for rescuing cancer cells from persistent ER stress, supporting FADS1 as a new therapeutic target.

The etiology of major depressive disorder (MDD) remains poorly understood. Our previous studies suggest a role for the aryl hydrocarbon receptor (AhR) in depression. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is a toxic environmental contaminant, with a high AhR binding affinity, and an established benchmark for assessing AhR activity. Therefore, this study examined the effect of TCDD on depression-like behaviors. Female mice were fed standard chow or a high-fat diet (HFD) for 11 weeks, and their weight was recorded. Subsequently, they were tested for baseline sucrose preference and splash test grooming. Then, TCDD (0.1µg/kg/day) or vehicle was administered orally for 28 days, and mice were examined for their sucrose preference and performances in the splash test, forced swim test (FST), and Morris water maze (MWM) task. TCDD significantly decreased sucrose preference, increased FST immobility time, and decreased groom time in chow-fed mice. HFD itself significantly reduced sucrose preference. However, TCDD significantly increased FST immobility time and decreased groom time in HFD-fed mice. A small decrease in bodyweight was observed only at the fourth week of daily TCDD administration in chow-fed mice, and no significant effects of TCDD on bodyweights were observed in HFD-fed mice. TCDD did not have a significant effect on spatial learning in the MWM. Thus, this study demonstrated that TCDD induces a depression-like state, and the effects were not due to gross lethal toxicity. This study further suggests that more studies should examine a possible role for AhR and AhR-active environmental pollutants in precipitating or worsening MDD.

After birth, the development of secondary lymphoid tissues (SLTs) in the colon is dependent on the expression of the aryl hydrocarbon receptor (AhR) in immune cells as a response to the availability of AhR ligands. However, little is known about how AhR activity from intestinal epithelial cells (IECs) may influence the development of tertiary lymphoid tissues (TLTs). As organized structures that develop at sites of inflammation or infection during adulthood, TLTs serve as localized centers of adaptive immune responses, and their presence has been associated with the resolution of inflammation and tumorigenesis in the colon. Here, we investigated the effect of the conditional loss of AhR activity in IECs in the formation and immune cell composition of TLTs in a model of acute inflammation. In females, loss of AhR activity in IECs reduced the formation of TLTs without significantly changing disease outcomes or immune cell composition within TLTs. In males lacking AhR expression in IECs, increased disease activity index, lower expression of functional-IEC genes, increased number of TLTs, increased T-cell density, and lower B- to T-cell ratio were observed. These findings may represent an unfavorable prognosis when exposed to dextran sodium sulfate (DSS)-induced epithelial damage compared with females. Sex and loss of IEC AhR also resulted in changes in microbial populations in the gut. Collectively, these data suggest that the formation of TLTs in the colon is influenced by sex and AhR expression in IECs.NEW & NOTEWORTHY This is the first research of its kind to demonstrate a clear connection between biological sex and the development of tertiary lymphoid tissues (TLT) in the colon. In addition, the research finds that in a preclinical model of inflammatory bowel disease, the expression of the aryl hydrocarbon receptor (AhR) influences the development of these structures in a sex-specific manner.

The equine gastrointestinal (GI) microbiome has been described in the context of various diseases. The observed changes, however, have not been linked to host function and therefore it remains unclear how specific changes in the microbiome alter cellular and molecular pathways within the GI tract. Further, non-invasive techniques to examine the host gene expression profile of the GI mucosa have been described in horses but not evaluated in response to interventions. Therefore, the objectives of our study were to (1) profile gene expression and metabolomic changes in an equine model of non-steroidal anti-inflammatory drug (NSAID)-induced intestinal inflammation and (2) apply computational data integration methods to examine host-microbiota interactions. Twenty horses were randomly assigned to 1 of 2 groups (n = 10): control (placebo paste) or NSAID (phenylbutazone 4.4 mg/kg orally once daily for 9 days). Fecal samples were collected on days 0 and 10 and analyzed with respect to microbiota (16S rDNA gene sequencing), metabolomic (untargeted metabolites), and host exfoliated cell transcriptomic (exfoliome) changes. Data were analyzed and integrated using a variety of computational techniques, and underlying regulatory mechanisms were inferred from features that were commonly identified by all computational approaches. Phenylbutazone induced alterations in the microbiota, metabolome, and host transcriptome. Data integration identified correlation of specific bacterial genera with expression of several genes and metabolites that were linked to oxidative stress. Concomitant microbiota and metabolite changes resulted in the initiation of endoplasmic reticulum stress and unfolded protein response within the intestinal mucosa. Results of integrative analysis identified an important role for oxidative stress, and subsequent cell signaling responses, in a large animal model of GI inflammation. The computational approaches for combining non-invasive platforms for unbiased assessment of host GI responses (e.g., exfoliomics) with metabolomic and microbiota changes have broad application for the field of gastroenterology.

Tertiary lymphoid tissues (TLTs) are adaptive immune structures that develop during chronic inflammation and may worsen or lessen disease outcomes in a context-specific manner. Immune cell activity governing TLT formation in the intestines is dependent on immune cell aryl hydrocarbon receptor (AhR) activation. Homeostatic immune cell activity in the intestines is further dependent on ligand activation of AhR in intestinal epithelial cells (IECs), yet whether AhR activation and signaling in IECs influences the formation of TLTs in the presence of dietary AhR ligands is not known. To this end, we used IEC-specific AhR deletion coupled with a mouse model of dextran sodium sulfate (DSS)-induced colitis to understand how dietary AhR ligand 3, 3′-diindolylmethane (DIM) influenced TLT formation. DIM consumption increased the size of TLTs and decreased T-cell aggregation to TLT sites in an IEC-specific manner. In DSS-exposed female mice, DIM consumption increased the expression of genes implicated in TLT formation (Interleukin-22, Il-22; CXC motif chemokine ligand 13, CXCL13) in an IEC AhR-specific manner. Conversely, in female mice without DSS exposure, DIM significantly reduced the expression of Il-22 or CXCL13 in iAhRKO mice, but this effect was not observed in WT animals. Our findings suggest that DIM affects the immunological landscape of TLT formation during DSS-induced colitis in a manner contingent on AhR expression in IECs and biological sex. Further investigations into specific immune cell activity, IEC-specific AhR signaling pathways, and dietary AhR ligand-mediated effects on TLT formation are warranted.

Objectives: Preclinical studies show that the concurrent consumption of dietary fiber and n-3 polyunsaturated fatty acids from fish oil reduce colon tumor formation. However, no controlled dietary intervention has explored the combined effects of these two dietary components on CRC risk factors in humans. This pilot study investigated the synergistic effects of fish oil and fermentable fiber on the gut microbiome composition in adult humans. Methods: In a randomized controlled crossover pilot study in 30 men and women ages 50-75, we tested the combined effects of supplemental fermentable fiber (35 g/d soluble corn fiber) and eicosapentaenoic acid, and docosahexaenoic acid (EPA+DHA as fish oil, 6 g/d) compared to similar doses of maltodextrin and corn oil. Blood samples were collected at baseline and end of each 30-day intervention period. Serum was used to measure phospholipid fatty acids. Participants collected stool samples at the beginning, middle, and end of each period, which were used to assess the microbiome via (16S rRNAseq) and short chain fatty acids (SCFA) via GC-MS. Results: Serum DHA and EPA concentrations were statistically significantly higher after fish oil and fiber vs the comparison supplements (P < 0.0001), suggesting good adherence to the intervention. Microbiome analysis revealed significant differences in beta-diversity (P = 0.001). Twenty-six of 116 genera analyzed, were significantly different between the two intervention (FDR < 0.05) and are known SCFA producers. Although fecal butyrate concentrations were not statistically significant different between intervention periods (P = 0.9), the bacteria identified play a key role in fermenting dietary fiber and producing butyrate. Conclusions: Intervention differences in the microbial ecosystem underscore the role of dietary components in modulating the gut microbiota. The insights from this study lay a foundation for future research focused on dietary interventions for CRC prevention, emphasizing the importance of diet-gut microbiome interactions. Funding Sources: Study supported by NCI. Supplements provided by Tate & Lyle Solutions USA LLC and Omega3 Innovations.

We consider the problem of clustering grouped data with possibly non-exchangeable groups whose dependencies can be characterized by a known directed acyclic graph. To allow the sharing of clusters among the non-exchangeable groups, we propose a Bayesian nonparametric approach, termed graphical Dirichlet process, that jointly models the dependent group-specific random measures by assuming each random measure to be distributed as a Dirichlet process whose concentration parameter and base probability measure depend on those of its parent groups. The resulting joint stochastic process respects the Markov property of the directed acyclic graph that links the groups. We characterize the graphical Dirichlet process using a novel hypergraph representation as well as the stick-breaking representation, the restaurant-type representation, and the representation as a limit of a finite mixture model. We develop an efficient posterior inference algorithm and illustrate our model with simulations and a real grouped single-cell dataset.