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.

Isoflavone is a secondary metabolite of the soybean phenylpropyl biosynthesis pathway with physiological activity and is beneficial to human health. In this study, the isoflavone content of 205 soybean germplasm resources from 3 locations in 2020 showed wide phenotypic variation. A joint genome-wide association study (GWAS) and weighted gene coexpression network analysis (WGCNA) identified 33 single-nucleotide polymorphisms and 11 key genes associated with soybean isoflavone content. Gene ontology enrichment analysis, gene coexpression, and haplotype analysis revealed natural variations in the Glyma.12G109800 (GmOMT7) gene and promoter region, with Hap1 being the elite haplotype. Transient overexpression and knockout of GmOMT7 increased and decreased the isoflavone content, respectively, in hairy roots. The combination of GWAS and WGCNA effectively revealed the genetic basis of soybean isoflavone and identified potential genes affecting isoflavone synthesis and accumulation in soybean, providing a valuable basis for the functional study of soybean isoflavone.

BACKGROUND: Knee osteoarthritis (KOA) is a common degenerative joint disease characterized by cartilage degradation, inflammation, and pain. Traditional Chinese Medicine, including JDJM (a herbal formula derived from the renowned Du Huo Ji Sheng Tang), has been used to alleviate symptoms of KOA, but its underlying mechanisms remain unclear.
OBJECTIVE: This study aims to elucidate the potential therapeutic mechanisms of JDJM in treating KOA through network pharmacology, weighted gene co-expression network analysis (WGCNA), molecular docking, and experimental validation in animal models.
METHODS: The active compounds of JDJM were identified through TCMSP database searches, and their potential targets were predicted using network pharmacology. WGCNA was employed to identify key modules and hub genes associated with KOA. Molecular docking was performed to assess the binding affinities of key compounds to critical inflammatory targets. Molecular dynamics (MD) simulations were used to evaluate the stability of the protein-ligand complexes. An experimental KOA model in rabbits was used to validate the therapeutic effects of JDJM. Histopathological examinations and inflammatory marker analyses were conducted to confirm the findings.
RESULTS: Network pharmacology and WGCNA analyses identified 21 key targets and pathways potentially involved in the therapeutic effects of JDJM. Molecular docking results showed that Glyasperin C had the highest docking scores with EGF and IL-1β, followed by Stigmasterol with IL-6, Myricanone with INS, and Sesamin with VEGFA. MD simulations confirmed the stability of these protein-ligand complexes, indicating strong and stable interactions. In the rabbit KOA model, JDJM treatment significantly improved knee joint morphology and reduced the levels of inflammatory markers, such as IL-6 and TNF-α. Histopathological analysis revealed reduced cartilage degradation and inflammation in the JDJM-treated group compared to controls.
CONCLUSIONS: JDJM exhibits promising anti-inflammatory and cartilage-protective effects, making it a potential therapeutic option for KOA patients. Further experimental and clinical studies are warranted to confirm these findings and translate them into clinical practice.

.

BACKGROUND: Delayed cerebral ischemia (DCI) is a common and serious complication of subarachnoid hemorrhage (SAH). Its pathogenesis is not fully understood. Here, we developed a predictive model based on peripheral blood biomarkers and validated the model using several bioinformatic multi-analysis methods.
METHODS: Six datasets were obtained from the GEO database. Characteristic genes were screened using weighted correlation network analysis (WGCNA) and differentially expressed genes. Three machine learning algorithms, elastic networks-LASSO, support vector machines (SVM-RFE) and random forests (RF), were also used to construct diagnostic prediction models for key genes. To further evaluate the performance and predictive value of the diagnostic models, nomogram model were constructed, and the clinical value of the models was assessed using Decision Curve Analysis (DCA), Area Under the Check Curve (AUC), Clinical Impact Curve (CIC), and validated in the mouse single-cell RNA-seq dataset. Mendelian randomization(MR) analysis explored the causal relationship between SAH and stroke, and the intermediate influencing factors. We validated this by retrospectively analyzing the qPCR levels of the most relevant genes in SAH and SAH-DCI patients. This experiment demonstrated a statistically significant difference between SAH and SAH-DCI and normal group controls. Finally, potential small molecule compounds interacting with the selected features were screened from the Comparative Toxicogenomics Database (CTD).
RESULTS: The fGSEA results showed that activation of Toll-like receptor signaling and leukocyte transendothelial cell migration pathways were positively correlated with the DCI phenotype, whereas cytokine signaling pathways and natural killer cell-mediated cytotoxicity were negatively correlated. Consensus feature selection of DEG genes using WGCNA and three machine learning algorithms resulted in the identification of six genes (SPOCK2, TRRAP, CIB1, BCL11B, PDZD8 and LAT), which were used to predict DCI diagnosis with high accuracy. Three external datasets and the mouse single-cell dataset showed high accuracy of the diagnostic model, in addition to high performance and predictive value of the diagnostic model in DCA and CIC. MR analysis looked at stroke after SAH independent of SAH, but associated with multiple intermediate factors including Hypertensive diseases, Total triglycerides levels in medium HDL and Platelet count. qPCR confirmed that significant differences in DCI signature genes were observed between the SAH and SAH-DCI groups. Finally, valproic acid became a potential therapeutic agent for DCI based on the results of target prediction and molecular docking of the characterized genes.
CONCLUSIONS: This diagnostic model can identify SAH patients at high risk for DCI and may provide potential mechanisms and therapeutic targets for DCI. Valproic acid may be an important future drug for the treatment of DCI.

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.

Eimeria tenella is the main pathogen responsible for coccidiosis in chickens. The life cycle of E. tenella is, arguably, the least complex of all Coccidia, with only one host. However, it presents different developmental stages, either in the environment or in the host and either intracellular or extracellular. Its signaling and metabolic pathways change with its different developmental stages. Until now, little is known about the developmental regulation and transformation mechanisms of its life cycle. In this study, protein profiles from the five developmental stages, including unsporulated oocysts (USO), partially sporulated (7 h) oocysts (SO7h), sporulated oocysts (SO), sporozoites (S) and second-generation merozoites (M2), were harvested using the label-free quantitative proteomics approach. Then the differentially expressed proteins (DEPs) for these stages were identified. A total of 314, 432, 689, and 665 DEPs were identified from the comparison of SO7h vs USO, SO vs SO7h, S vs SO, and M2 vs S, respectively. By conducting weighted gene coexpression network analysis (WGCNA), six modules were dissected. Proteins in blue and brown modules were calculated to be significantly positively correlated with the E. tenella developmental stages of sporozoites (S) and second-generation merozoites (M2), respectively. In addition, hub proteins with high intra-module degree were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway enrichment analyses revealed that hub proteins in blue modules were involved in electron transport chain and oxidative phosphorylation. Hub proteins in the brown module were involved in RNA splicing. These findings provide new clues and ideas to enhance our fundamental understanding of the molecular mechanisms underlying parasite development.

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.