Engineering acid-tolerant microbial strains is a cost-effective approach to overcoming acid stress during industrial fermentation. We previously constructed an acid-tolerant strain (Escherichia coli SC3124) with enhanced growth robustness and productivity under mildly acidic conditions by fine-tuning the expression of synthetic acid-tolerance module genes consisting of a proton-consuming acid resistance system (gadE), a periplasmic chaperone (hdeB), and ROS scavengers (sodB, katE). However, the precise acid-tolerance mechanism of E. coli SC3124 remained unclear. In this study, the growth of E. coli SC3124 under mild acid stress (pH 6.0) was determined. The final OD600 of E. coli SC3124 at pH 6.0 was 131% and 124% of that of the parent E. coli MG1655 at pH 6.8 and pH 6.0, respectively. Transcriptome analysis revealed the significant upregulation of the genes involved in oxidative phosphorylation, the tricarboxylic acid (TCA) cycle, and lysine-dependent acid-resistance system in E. coli SC3124 at pH 6.0. Subsequently, a weighted gene coexpression network analysis was performed to systematically determine the metabolic perturbations of E. coli SC3124 with mild acid treatment, and we extracted the gene modules highly associated with different acid traits. The results showed two biologically significant coexpression modules, and 263 hub genes were identified. Specifically, the genes involved in ATP-binding cassette (ABC) transporters, oxidative phosphorylation, the TCA cycle, amino acid metabolism, and purine metabolism were highly positively associated with mild acid stress responses. We propose that the overexpression of synthetic acid-tolerance genes leads to metabolic changes that confer mild acid stress resistance in E. coli. Integrated omics platforms provide valuable information for understanding the regulatory mechanisms of mild acid tolerance in E. coli and highlight the important roles of oxidative phosphorylation and ABC transporters in mild acid stress regulation. These findings offer novel insights to better the design of acid-tolerant chasses to synthesize value-added chemicals in a green and sustainable manner.

Patients with thymoma (THYM)-associated myasthenia gravis (MG) typically have a poor prognosis and recurring illness. This study aimed to discover important biomarkers associated with immune cell infiltration and THYM-associated MG (THYM-MG) development. Gene expression microarray data were downloaded from The Cancer Genome Atlas website (TCGA) and Gene Expression Omnibus (GEO). A total of 102 differentially expressed genes were investigated. According to the immune infiltration data, the distribution of Tfh cells, B cells, and CD4 T cells differed significantly between the THYM-MG and THYM-NMG groups. WGCNA derived 25 coexpression modules; one hub module (the blue module) strongly correlated with Tfh cells. Combining differential genes revealed 21 intersecting genes. LASSO analysis subsequently revealed 16 hub genes as potential THYM-MG biomarkers. ROC curve analysis of the predictive model revealed moderate diagnostic value. The association between the 16 hub genes and infiltrating immune cells was further evaluated in TIMER2.0 and the validation dataset. Draggability analysis identified the therapeutic target genes PTGS2 and ALB, along with significant drugs including Firocoxib, Alclofenac, Pyridostigmine, and Stavudine. This was validated through MD simulation, PCA, and MM-GBSA analyses. The interaction between numerous activated B cells and follicular helper T cells is closely associated with THYM-MG pathogenesis from a bioinformatics perspective. Hub genes (including SP6, SCUBE3, B3GNT7, and MAGEL2) may be downregulated in immune cells in THYM-MG and associated with progression.

BACKGROUND: The development and progression of hepatocellular carcinoma (HCC) have been reported to be associated with immune-related genes and the tumor microenvironment. Nevertheless, there are not enough prognostic biomarkers and models available for clinical use. Based on seven prognostic genes, this study calculated overall survival in patients with HCC using a prognostic survival model and revealed the immune status of the tumor microenvironment (TME).
OBJECTIVE: To develop a novel immune cell-related prognostic model of HCC and depict the basic profile of the immune response in HCC.
METHODS: We obtained clinical information and gene expression data of HCC from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets. TCGA and ICGC datasets were used for screening prognostic genes along with developing and validating a seven-gene prognostic survival model by weighted gene coexpression network analysis and least absolute shrinkage and selection operator regression with Cox regression. The relative analysis of tumor mutation burden (TMB), TME cell infiltration, immune checkpoints, immune therapy, and functional pathways was also performed based on prognostic genes.
RESULTS: Seven prognostic genes were identified for signature construction. Survival receiver operating characteristic curve analysis showed the good performance of survival prediction. TMB could be regarded as an independent factor in HCC survival prediction. There was a significant difference in stromal score, immune score, and estimate score between the high-risk and low-risk groups stratified based on the risk score derived from the seven-gene prognostic model. Several immune checkpoints, including VTCN1 and TNFSF9, were found to be associated with the seven prognostic genes and risk score. Different combinations of checkpoint blockade targeting inhibitory CTLA4 and PD1 receptors and potential chemotherapy drugs hold great promise for specific HCC therapies. Potential pathways, such as cell cycle regulation and metabolism of some amino acids, were also identified and analyzed.
CONCLUSIONS: The novel seven-gene (CYTH3, ENG, HTRA3, PDZD4, SAMD14, PGF, and PLN) prognostic model showed high predictive efficiency. The TMB analysis based on the seven genes could depict the basic profile of the immune response in HCC, which might be worthy of clinical application.

Biomarkers and pathways associated with renal ischemia reperfusion injury (IRI) had not been well unveiled. This study was intended to investigate and summarize the regulatory networks for related hub genes. Besides, the immunological micro-environment features were evaluated and the correlations between immune cells and hub genes were also explored.
GSE98622 containing mouse samples with multiple IRI stages and controls was collected from the GEO database. Differentially expressed genes (DEGs) were recognized by the R package limma, and the GO and KEGG analyses were conducted by DAVID. Gene set variation analysis (GSVA) and weighted gene coexpression network analysis (WGCNA) had been implemented to uncover changed pathways and gene modules related to IRI. Besides the known pathways such as apoptosis pathway, metabolic pathway, and cell cycle pathways, some novel pathways were also discovered to be critical in IRI. A series of novel genes associated with IRI was also dug out. An IRI mouse model was constructed to validate the results.
The well-known IRI marker genes (Kim1 and Lcn2) and novel hub genes (Hbegf, Serpine2, Apbb1ip, Trip13, Atf3, and Ncaph) had been proved by the quantitative real-time polymerase chain reaction (qRT-PCR). Thereafter, miRNAs targeted to the dysregulated genes were predicted and the miRNA-target network was constructed. Furthermore, the immune infiltration for these samples was predicted and the results showed that macrophages infiltrated to the injured kidney to affect the tissue repair or fibrosis. Hub genes were significantly positively or negatively correlated with the macrophage abundance indicating they played a crucial role in macrophage infiltration.
Consequently, the pathways, hub genes, miRNAs, and the immune microenvironment may explain the mechanism of IRI and might be the potential targets for IRI treatments.

UNASSIGNED: The purpose of this study was to identify potential biomarkers in the tubulointerstitium of focal segmental glomerulosclerosis (FSGS) and comprehensively analyze its mRNA-miRNA-lncRNA/circRNA network.
UNASSIGNED: The expression data (GSE108112 and GSE200818) were downloaded from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/). Identification and enrichment analysis of differentially expressed genes (DEGs) were performed. the PPI networks of the DEGs were constructed and classified using the Cytoscape molecular complex detection (MCODE) plugin. Weighted gene coexpression network analysis (WGCNA) was used to identify critical gene modules. Least absolute shrinkage and selection operator regression analysis were used to screen for key biomarkers of the tubulointerstitium in FSGS, and the receiver operating characteristic curve was used to determine their diagnostic accuracy. The screening results were verified by quantitative real-time-PCR (qRT-PCR) and Western blot. The transcription factors (TFs) affecting the hub genes were identified by Cytoscape iRegulon. The mRNA-miRNA-lncRNA/circRNA network for identifying potential biomarkers was based on the starBase database.
UNASSIGNED: A total of 535 DEGs were identified. MCODE obtained eight modules. The green module of WGCNA had the greatest association with the tubulointerstitium in FSGS. PPARG coactivator 1 alpha (PPARGC1A) was screened as a potential tubulointerstitial biomarker for FSGS and verified by qRT-PCR and Western blot. The TFs FOXO4 and FOXO1 had a regulatory effect on PPARGC1A. The ceRNA network yielded 17 miRNAs, 32 lncRNAs, and 50 circRNAs.
UNASSIGNED: PPARGC1A may be a potential biomarker in the tubulointerstitium of FSGS. The ceRNA network contributes to the comprehensive elucidation of the mechanisms of tubulointerstitial lesions in FSGS.

BACKGROUND: Gastric cancer (GC) is a heterogeneous malignancy with variable clinical outcomes. The immune system has been implicated in GC development and progression, highlighting the importance of immune-related gene expression patterns and their prognostic significance.
OBJECTIVE: This study aimed to identify differentially expressed immune-related genes (DEIRGs) and establish a prognostic index for GC patients using comprehensive bioinformatic analyses.
METHODS: We integrated RNA sequencing data from multiple databases and identified DEIRGs by overlapping differentially expressed genes with immune-related genes. Functional enrichment analysis was performed to uncover the biological processes and signaling pathways associated with DEIRGs. We conducted a Weighted Gene Co-expression Network Analysis (WGCNA) to identify key gene modules related to with GC. Cox regression analysis was conducted to determine independent prognostic DEIRGs for overall survival prediction. Based on these findings, we developed an immune-related gene prognostic index (IRGPI) based on these findings. The prognostic value of the IRGPI was validated using survival analysis and an independent validation cohort. Functional enrichment analysis, gene mutation analysis, and immune cell profiling were performed to gain insights into the biological functions and immune characteristics associated with the IRGPI-based subgroups.
RESULTS: We identified 493 DEIRGs significantly enriched in immune-related biological processes and signaling pathways associated with GC. WGCNA analysis revealed a significant module (turquoise module) associated with GC, revealing potential therapeutic targets. Cox regression analysis identified RNASE2, CGB5, CTLA4, and DUSP1 as independent prognostic DEIRGs. The IRGPI, incorporating the expression levels of these genes, demonstrated significant prognostic value in predicting overall survival. The IRGPI-based subgroups exhibited distinct biological functions, genetic alterations, and immune cell compositions.
CONCLUSIONS: Our study identified DEIRGs and established a prognostic index (IRGPI) for GC patients. The IRGPI exhibited promising prognostic potential and provided insights into GC tumor biology and immune characteristics. These findings have implications for guiding therapeutic strategies.

UNASSIGNED: Flavonoids, as secondary metabolites in plants, play important roles in many biological processes and responses to environmental factors.
UNASSIGNED: Apricot fruits are rich in flavonoid compounds, and in this study, we performed a combined metabolomic and transcriptomic analysis of orange flesh (JN) and white flesh (ZS) apricot fruits.
UNASSIGNED: A total of 222 differentially accumulated flavonoids (DAFs) and 15855 differentially expressed genes (DEGs) involved in flavonoid biosynthesis were identified. The biosynthesis of flavonoids in apricot fruit may be regulated by 17 enzyme-encoding genes, namely PAL (2), 4CL (9), C4H (1), HCT (15), C3\'H (4), CHS (2), CHI (3), F3H (1), F3\'H (CYP75B1) (2), F3\'5\'H (4), DFR (4), LAR (1), FLS (3), ANS (9), ANR (2), UGT79B1 (6) and CYP81E (2). A structural gene-transcription factor (TF) correlation analysis yielded 3 TFs (2 bHLH, 1 MYB) highly correlated with 2 structural genes. In addition, we obtained 26 candidate genes involved in the biosynthesis of 8 differentially accumulated flavonoids metabolites in ZS by weighted gene coexpression network analysis. The candidate genes and transcription factors identified in this study will provide a highly valuable molecular basis for the in-depth study of flavonoid biosynthesis in apricot fruits.

BACKGROUND: Regulatory T cells (Tregs) and natural killer (NK) cells play an essential role in the development of bladder urothelial carcinoma (BUC).
OBJECTIVE: To construct a prognosis-related model to judge the prognosis of patients with bladder cancer, meanwhile, predict the sensitivity of patients to chemotherapy and immunotherapy.
METHODS: Bladder cancer information data was obtained from The Cancer Genome Atlas and GSE32894. The CIBERSORT was used to calculate the immune score of each sample. Weighted gene co-expression network analysis was used to find genes that will have the same or similar expression patterns. Subsequently, multivariate cox regression and lasso regression was used to further screen prognosis-related genes. The prrophetic package was used to predict phenotype from gene expression data, drug sensitivity of external cell line and predict clinical data.
RESULTS: The stage and risk scores are independent prognostic factors in patients with BUC. Mutations in FGFR3 lead to an increase in Tregs percolation and affect the prognosis of the tumor, and additionally, EMP1, TCHH and CNTNAP3B in the model are mainly positively correlated with the expression of immune checkpoints, while CMTM8, SORT1 and IQSEC1 are negatively correlated with immune checkpoints and the high-risk group had higher sensitivity to chemotherapy drugs.
CONCLUSIONS: Prognosis-related models of bladder tumor patients, based on Treg and NK cell percolation in tumor tissue. In addition to judging the prognosis of patients with bladder cancer, it can also predict the sensitivity of patients to chemotherapy and immunotherapy. At the same time, patients were divided into high and low risk groups based on this model, and differences in genetic mutations were found between the high and low risk groups.

Background: Kawasaki disease (KD) is an acute vasculitis, that is, the leading cause of acquired heart disease in children, with approximately 10%-20% of patients with KD suffering intravenous immunoglobulin (IVIG) resistance. Although the underlying mechanism of this phenomenon remains unclear, recent studies have revealed that immune cell infiltration may associate with its occurrence. Methods: In this study, we downloaded the expression profiles from the GSE48498 and GSE16797 datasets in the Gene Expression Omnibus database, analyzed differentially expressed genes (DEGs), and intersected the DEGs with the immune-related genes downloaded from the ImmPort database to obtain differentially expressed immune-related genes (DEIGs). Then CIBERSORT algorithm was used to calculate the immune cell compositions, followed by the WGCNA analysis to identify the module genes associated with immune cell infiltration. Next, we took the intersection of the selected module genes and DEIGs, then performed GO and KEGG enrichment analysis. Moreover, ROC curve validation, Spearman analysis with immune cells, TF, and miRNA regulation network, and potential drug prediction were implemented for the finally obtained hub genes. Results: The CIBERSORT algorithm showed that neutrophil expression was significantly higher in IVIG-resistant patients compared to IVIG-responsive patients. Next, we got differentially expressed neutrophil-related genes by intersecting DEIGs with neutrophil-related module genes obtained by WGCNA, for further analysis. Enrichment analysis revealed that these genes were associated with immune pathways, such as cytokine-cytokine receptor interaction and neutrophil extracellular trap formation. Then we combined the PPI network in the STRING database with the MCODE plugin in Cytoscape and identified 6 hub genes (TLR8, AQP9, CXCR1, FPR2, HCK, and IL1R2), which had good diagnostic performance in IVIG resistance according to ROC analysis. Furthermore, Spearman\'s correlation analysis confirmed that these genes were closely related to neutrophils. Finally, TFs, miRNAs, and potential drugs targeting the hub genes were predicted, and TF-, miRNA-, and drug-gene networks were constructed. Conclusion: This study found that the 6 hub genes (TLR8, AQP9, CXCR1, FPR2, HCK, and IL1R2) were significantly associated with neutrophil cell infiltration, which played an important role in IVIG resistance. In a word, this work rendered potential diagnostic biomarkers and prospective therapeutic targets for IVIG-resistant patients.

UNASSIGNED: The aim of our study was to identify key biomarkers of glomeruli in focal glomerulosclerosis (FSGS) and analyze their relationship with the infiltration of immune cells.
UNASSIGNED: The expression profiles (GSE108109 and GSE200828) were obtained from the GEO database. The differentially expressed genes (DEGs) were filtered and analyzed by gene set enrichment analysis (GSEA). MCODE module was constructed. Weighted gene coexpression network analysis (WGCNA) was performed to obtain the core gene modules. Least absolute shrinkage and selection operator (LASSO) regression was applied to identify key genes. ROC curves were employed to explore their diagnostic accuracy. Transcription factor prediction of the key biomarkers was performed using the Cytoscape plugin IRegulon. The analysis of the infiltration of 28 immune cells and their correlation with the key biomarkers were performed.
UNASSIGNED: A total of 1474 DEGs were identified. Their functions were mostly related to immune-related diseases and signaling pathways. MCODE identified five modules. The turquoise module of WGCNA had significant relevance to the glomerulus in FSGS. TGFB1 and NOTCH1 were identified as potential key glomerular biomarkers in FSGS. Eighteen transcription factors were obtained from the two hub genes. Immune infiltration showed significant correlations with T cells. The results of immune cell infiltration and their relationship with key biomarkers implied that NOTCH1 and TGFB1 were enhanced in immune-related pathways.
UNASSIGNED: TGFB1 and NOTCH1 may be strongly correlated with the pathogenesis of the glomerulus in FSGS and are new candidate key biomarkers. T-cell infiltration plays an essential role in the FSGS lesion process.