OBJECTIVE: To explore the key target molecules and potential mechanisms of oridonin against non-small cell lung cancer (NSCLC). METHODS: The target molecules of oridonin were retrieved from SEA, STITCH, SuperPred and TargetPred databases; target genes associated with the treatment of NSCLC were retrieved from GeneCards, DisGeNET and TTD databases. Then, the overlapping target molecules between the drug and the disease were identified. The protein-protein interaction (PPI) was constructed using the STRING database according to overlapping targets, and Cytoscape was used to screen for key targets. Molecular docking verification were performed using AutoDockTools and PyMOL software. Using the DAVID database, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were conducted. The impact of oridonin on the proliferation and apoptosis of NSCLC cells was assessed using cell counting kit-8, cell proliferation EdU image kit, and Annexin V-FITC/PI apoptosis kit respectively. Moreover, real-time quantitative PCR and Western blot were used to verify the potential mechanisms. RESULTS: Fifty-six target molecules and 12 key target molecules of oridonin involved in NSCLC treatment were identified, including tumor protein 53 (TP53), Caspase-3, signal transducer and activator of transcription 3 (STAT3), mitogen-activated protein kinase kinase 8 (MAPK8), and mammalian target of rapamycin (mTOR). Molecular docking showed that oridonin and its key target molecules bind spontaneously. GO and KEGG enrichment analyses revealed cancer, apoptosis, phosphoinositide-3 kinase/protein kinase B (PI3K/Akt), and other signaling pathways. In vitro experiments showed that oridonin inhibited the proliferation, induced apoptosis, downregulated the expression of Bcl-2 and Akt, and upregulated the expression of Caspase-3. CONCLUSION Oridonin can act on multiple targets and pathways to exert its inhibitory effects on NSCLC, and its mechanism may be related to upregulating the expression of Caspase-3 and downregulating the expressions of Akt and Bcl-2.
Diterpenes, Kaurane/chemistry* ; Carcinoma, Non-Small-Cell Lung/pathology* ; Humans ; Network Pharmacology ; Lung Neoplasms/pathology* ; Cell Proliferation/drug effects* ; Apoptosis/drug effects* ; Molecular Docking Simulation ; Protein Interaction Maps/drug effects* ; Cell Line, Tumor ; Signal Transduction/drug effects* ; Gene Expression Regulation, Neoplastic/drug effects* ; Reproducibility of Results ; Gene Ontology

OBJECTIVES: Hepatocellular carcinoma is one of the most common primary malignant tumors with the third highest mortality rate worldwide. This study aims to identify key genes associated with hepatocellular carcinoma prognosis using the Gene Expression Omnibus (GEO) database and provide a theoretical basis for discovering novel prognostic biomarkers for hepatocellular carcinoma. METHODS: Hepatocellular carcinoma-related datasets were retrieved from the GEO database. Differentially expressed genes (DEGs) were identified using the GEO2R tool. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). A protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), and key genes were identified using Cytoscape software. The University of Alabama at Birmingham Cancer Data Analysis Resource (UALCAN) was used to analyze the expression levels of key genes in normal and hepatocellular carcinoma tissues, as well as their associations with pathological grade, clinical stage, and patient survival. The Human Protein Atlas (THPA) was used to further validate the impact of key genes on overall survival. Expression levels of key genes in the blood of hepatocellular carcinoma patients were evaluated using the expression atlas of blood-based biomarkers in the early diagnosis of cancers (BBCancer). RESULTS: A total of 78 DEGs were identified from the GEO database. GO and KEGG analyses indicated that these genes may contribute to hepatocellular carcinoma progression by promoting cell division and regulating protein kinase activity. Sixteen key genes were screened via Cytoscape and validated using UALCAN and THPA. These genes were overexpressed in hepatocellular carcinoma tissues and were associated with disease progression and poor prognosis. Finally, BBCancer analysis showed that ASPM and NCAPG were also elevated in the blood of hepatocellular carcinoma patients. CONCLUSIONS This study identified 16 key genes as potential prognostic biomarkers for hepatocellular carcinoma, among which ASPM and NCAPG may serve as promising blood-based markers for hepatocellular carcinoma.
Humans ; Carcinoma, Hepatocellular/mortality* ; Liver Neoplasms/pathology* ; Prognosis ; Computational Biology/methods* ; Protein Interaction Maps/genetics* ; Biomarkers, Tumor/genetics* ; Gene Expression Regulation, Neoplastic ; Gene Expression Profiling ; Gene Ontology ; Databases, Genetic

OBJECTIVES: Osteoarthritis (OA) and sarcopenia are significant health concerns in the elderly, substantially impacting their daily activities and quality of life. However, the relationship between them remains poorly understood. This study aims to uncover common biomarkers and pathways associated with both OA and sarcopenia. METHODS: Gene expression profiles related to OA and sarcopenia were retrieved from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between disease and control groups were identified using R software. Common DEGs were extracted via Venn diagram analysis. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to identify biological processes and pathways associated with shared DEGs. Protein-protein interaction (PPI) networks were constructed, and candidate hub genes were ranked using the maximal clique centrality (MCC) algorithm. Further validation of hub gene expression was performed using 2 independent datasets. Receiver operating characteristic (ROC) curve analysis was used to evaluate the predictive value of key genes for OA and sarcopenia. Mouse models of OA and sarcopenia were established. Hematoxylin-eosin and Safranin O/Fast Green staining were used to validate the OA model. The sarcopenia model was validated via rotarod testing and quadriceps muscle mass measurement. Real-time reverse transcription PCR (real-time RT-PCR) was employed to assess the mRNA expression levels of candidate key genes in both models. Gene set enrichment analysis (GSEA) was conducted to identify pathways associated with the selected shared key genes in both diseases. RESULTS: A total of 89 common DEGs were identified in the gene expression profiles of OA and sarcopenia, including 76 upregulated and 13 downregulated genes. These 89 DEGs were significantly enriched in protein digestion and absorption, the PI3K-Akt signaling pathway, and extracellular matrix-receptor interaction. PPI network analysis and MCC algorithm analysis of the 89 common DEGs identified the top 17 candidate hub genes. Based on the differential expression analysis of these 17 candidate hub genes in the validation datasets, AEBP1 and COL8A2 were ultimately selected as the common key genes for both diseases, both of which showed a significant upregulation trend in the disease groups (all P<0.05). The value of area under the curve (AUC) for AEBP1 and COL8A2 in the OA and sarcopenia datasets were all greater than 0.7, indicating that both genes have potential value in predicting OA and sarcopenia. Real-time RT-PCR results showed that the mRNA expression levels of AEBP1 and COL8A2 were significantly upregulated in the disease groups (all P<0.05), consistent with the results observed in the bioinformatics analysis. GSEA revealed that AEBP1 and COL8A2 were closely related to extracellular matrix-receptor interaction, ribosome, and oxidative phosphorylation in OA and sarcopenia. CONCLUSIONS AEBP1 and COL8A2 have the potential to serve as common biomarkers for OA and sarcopenia. The extracellular matrix-receptor interaction pathway may represent a potential target for the prevention and treatment of both OA and sarcopenia.
Sarcopenia/genetics* ; Osteoarthritis/genetics* ; Computational Biology/methods* ; Humans ; Protein Interaction Maps/genetics* ; Animals ; Mice ; Gene Expression Profiling ; Gene Ontology ; Transcriptome ; Male ; Signal Transduction/genetics* ; Gene Regulatory Networks

BACKGROUND: PANoptosis has the features of pyroptosis, apoptosis, and necroptosis. Numerous studies have confirmed the diverse roles of various types of cell death in acute liver failure (ALF), but limited attention has been given to the crosstalk among them. In this study, we aimed to explore the role of PANoptosis in ALF and uncover new targets for its prevention or treatment. METHODS: Three ALF-related datasets (GSE14668, GSE62029, and GSE74000) were downloaded from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs). Hub genes were identified through intersecting DEGs, genes obtained from weighted gene co-expression network analysis (WGCNA), and genes related to PANoptosis. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), protein‍‒‍protein interaction (PPI) analyses and gene set enrichment analysis (GSEA) were performed to determine functional roles. Verification was performed using an ALF mouse model. RESULTS: Our results showed that expression of seven hub genes (B-cell lymphoma-2-modifying factor (BMF), B-cell lymphoma-2-interacting protein 3-like (BNIP3L), Caspase-1 (CASP1), receptor-interacting protein kinase 3 (RIPK3), uveal autoantigen with coiled-coil domains and ankyrin repeats protein (UACA), uncoordinated-5 homolog B receptor (UNC5B), and Z-DNA-binding protein 1 (ZBP1)) was up-regulated in liver samples of patients. However, in the ALF mouse model, the expression of BNIP3L, RIPK3, phosphorylated RIPK3 (P-RIPK3), UACA, and cleaved caspase-1 was up-regulated, while the expression of CASP1 and UNC5B was down-regulated. The expression of ZBP1 and BMF increased only during the development of ALF, and there was no significant change in the end stage. Immunofluorescence of mouse liver tissue showed that macrophages expressed all seven markers. Western blot results showed that pyroptosis, apoptosis, and necroptosis were always involved in lipopolysaccharide (LPS)/ d-galactosamine (d-gal)‍-induced ALF mice. The ALF cell model showed that bone marrow-derived macrophages (BMDMs) form PANoptosomes after LPS stimulation. CONCLUSIONS: Our results suggest that PANoptosis of macrophages promotes the development of ALF. The seven new ALF biomarkers identified and validated in this study may contribute to further investigation of diagnostic markers or novel therapeutic targets of ALF.
Animals ; Liver Failure, Acute/genetics* ; Computational Biology ; Mice ; Pyroptosis/genetics* ; Humans ; Protein Interaction Maps ; Apoptosis/genetics* ; Necroptosis/genetics* ; Gene Regulatory Networks ; Gene Ontology ; Gene Expression Profiling ; Disease Models, Animal

OBJECTIVES: To identify cuproptosis- and ferroptosis-related genes involved in nonalcoholic fatty liver disease and to determine the diagnostic value of hub genes. METHODS: The gene expression dataset GSE89632 was retrieved from the Gene Expression Omnibus database to identify differentially expressed genes (DEGs) between the non-alcoholic steatohepatitis (NASH) group and the healthy group using the 'limma' package in R software and weighted gene co-expression network analysis. Gene ontology, kyoto encyclopedia of genes and genomes pathway, and single-sample gene set enrichment analyses were performed to identify functional enrichment of DEGs. Ferroptosis- and cuproptosis-related genes were obtained from the FerrDb V2 database and available literatures, respectively. A combined signature for cuproptosis- and ferroptosis-related genes, called CRF, was constructed using the STRING database. Hub genes were identified by overlapping DEGs, WGCNA-derived key genes, and combined signature CRF genes, and validated using the GSE109836 and GSE227714 datasets and real-time quantitative polymerase chain reaction. A nomogram of NASH diagnostic model was established utilizing the 'rms' package in R software based on the hub genes, and the diagnostic value of hub genes was assessed using receiver operating characteristic curve analysis. In addition, immune cell infiltration in NASH versus healthy controls was examined using the CIBERSORT algorithm. The relationships among various infiltrated immune cells were explored with Spearman's correlation analysis. RESULTS: Analysis of GSE89632 identified 236 DEGs between the NASH group and the healthy group. WGCNA highlighted 8 significant modules and 11,095 pivotal genes, of which 330 genes constituted CRF. Intersection analysis identified IL6, IL1B, JUN, NR4A1, and PTGS2 as hub genes. The hub genes were all downregulated in the NASH group, and this result was further verified by the NASH validation dataset and real-time quantitative polymerase chain reaction. Receiver operating characteristic curve analysis confirmed the diagnostic efficacy of these hub genes with areas under the curve of 0.985, 0.941, 1.000, 0.967, and 0.985, respectively. Immune infiltration assessment revealed that gamma delta T cells, M1 macrophages, M2 macrophages, and resting mast cells were predominantly implicated. CONCLUSIONS Our investigation underscores the significant association of cuproptosis- and ferroptosis-related genes, specifically IL6, IL1B, JUN, NR4A1, and PTGS2, with NASH. These findings offer novel insights into the pathogenesis of NASH, potentially guiding future diagnostic and therapeutic strategies.
Non-alcoholic Fatty Liver Disease/pathology* ; Humans ; Ferroptosis/genetics* ; Copper/metabolism* ; Gene Ontology ; Gene Expression Profiling

OBJECTIVES: To explore ferroptosis-related genes in osteoporosis through bioinformatic analysis and in vitro study. METHODS: Osteoporosis-related genes were identified from dataset GSE35958 in the Gene Expression Omnibus database; and the ferroptosis-related genes were identified from the FerrDb database. These were intersected with the differentially expressed genes in GSE35958 to obtain ferroptosis-related genes in osteoporosis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed for the differentially expressed genes. And Spearman correlation and protein-protein interaction network analysis were performed. Then, the hub genes of ferroptosis in osteoporosis were screened by Degree, MNC, EPC, MCC and DMNC in Cytoscape software CytoHubba plugin; and analyzed with receiver operating characteristic (ROC) curves. The bone marrow mesenchymal stem cells from osteoporosis patients (osteoporosis group) and non-osteoporosis patients (control group) were subjected to quantitative reverse transcription polymerase chain reaction to detect the messenger RNA expression of ferroptosis hub genes in both groups. RESULTS: A total of 32 differentially expressed genes related to ferroptosis in osteoporosis were identified, including 26 up-regulated genes and 6 down-regulated genes. GO enrichment analysis showed that the identified genes were mainly involved in intercellular adhesion, lipid metabolism and cytokine response. KEGG enrichment analysis showed that the genes were mainly involved in signaling pathways of adhesive plaques, MAPK, PI3K-Akt, and Wnt. Spearman correlation analysis showed correlation among differentially expressed genes. Six hub genes for ferroptosis in osteoporosis were obtained, namely MAPK3, CDKN1A, MAP1LC3A, TNF, RELA, and TGF-β1. ROC curve analysis showed that these hub genes had good diagnostic performance in osteoporosis and may become potential biomarkers of osteoporosis. In vitro experiments confirmed significant differences in these hub genes between the control group and the osteoporosis group (all P<0.05). CONCLUSIONS This study has identified six ferroptosis-related hub genes in osteoporosis, which may be used as novel biomarkers for the early diagnosis and treatment of osteoporosis.
Osteoporosis/genetics* ; Humans ; Computational Biology ; Ferroptosis/genetics* ; Protein Interaction Maps/genetics* ; Gene Ontology ; Mesenchymal Stem Cells/metabolism* ; Gene Expression Profiling ; Databases, Genetic

Objective To identify the possibility of IgG Fc binding protein (FCGBP) acting as a prognostic marker of low-grade glioma (LGG) and its correlation with immune infiltration. Methods The expression of FCGBP was analyzed in pan-cancer using The Cancer Genome Atlas (TCGA), Genotypic tissue expression (GTEX), and China Glioma Genome Atlas (CGGA) database. Then, GSE15824 and GSE68848 datasets were selected for further verification. And gene expression Profile Interaction analysis (GEPIA) database and R language were used to analyze the relationship between FCGBP and survival prognosis. Metascape and GSEA were used for functional annotation and enrichment analysis. Finally, the expression of FCGBP gene in LGG immune microenvironment and its correlation with immune cells were analyzed by TIMER database. Results FCGBP was highly expressed in LGG tissues, indicating poor prognosis of LGG patients. Receiver operating characteristic (ROC) curve analysis and COX analysis showed that FCGBP was an independent risk factor for the prognosis of LGG. Moreover, Gene Ontology (GO) demonstrated that FCGBP was involved in cell metabolism, localization, positive, and negative regulation of biological processes, as well as biological adhesion, response to viral and microbial stimulation, and inflammation. GSEA pathway enrichment analysis showed that FCGBP was significantly correlated with Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway, Toll-like receptor (TLR) pathway, chemokine pathway, and P53 pathway. In addition, FCGBP expression was positively correlated with the expression of most immune cells in the immune microenvironment of LGG. Conclusion The high expression of FCGBP in LGG is a risk factor for survival and prognosis, and it is positively correlated with the expression of immune cells.
Humans ; Prognosis ; Glioma/genetics* ; China ; Gene Ontology ; Immunoglobulin G ; Tumor Microenvironment ; Cell Adhesion Molecules

This study explored the anticoagulant material basis and mechanism of Trichosanthis Semen and its shell and kernel based on spectrum-effect relationship-integrated molecular docking. High performance liquid chromatography(HPLC) fingerprints of Trichosanthis Semen and its shell and kernel were established. Prothrombin time(PT) and activated partial thromboplastin time(APTT) in mice in the low-and high-dose(5, 30 g·kg~(-1), respectively) Trichosanthis Semen, the shell, and kernel groups were determined as the coagulation markers. The spectrum-effect relationship and anticoagulant material basis of Trichosanthis Semen and its shell and kernel were analyzed with mean value calculation method of Deng's correlation degree(MATLAB) and the common effective component cluster was obtained. Then the common targets of the component cluster and coagulation were retrieved from TCMSP, Swiss-TargetPrediction, GenCLiP3, GeneCards, and DAVID, followed by Gene Ontology(GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment of the targets. The main anticoagulant molecular mechanism of the component cluster was verified by SYBYL-X 2.1.1. The spectrum-effect relationship of Trichosanthis Semen and its shell and kernel was in positive correlation with the dosage. The contribution of each component to anticoagulation was not the same, suggesting that the material basis for anticoagulation was different, but they have common effective components(i.e. common material basis), such as adenine(peak 3), uracil(peak 4), hypoxanthine(peak 6), xanthine(peak 9), and adenosine(peak 11). Network pharmacology showed that these components can act on multiple target proteins such as NOS3, KDR, and PTGS2, and exert anticoagulant effect through multiple pathways such as VEGF signaling pathway. They involved the biological functions such as proteolysis, cell component such as cytosol, and molecular functions. The results of molecular docking showed that the binding free energy of these components with NOS3(PDB ID: 1 D0 C), KDR(PDB ID: 5 AMN), and PTGS2(PDB ID: 4 COX) was ≤-5 kJ·mol~(-1), and the docking conformations were stable. Spectrum-effect relationship-integrated molecular docking can be used for the optimization, virtual screening, and verification of complex chemical and biological information of Chinese medicine. Trichosanthis Semen and its shell and kernel have the common material basis for anticoagulation and they exert the anticoagulant through multiple targets and pathways.
Animals ; Anticoagulants/pharmacology* ; Drugs, Chinese Herbal/pharmacology* ; Gene Ontology ; Mice ; Molecular Docking Simulation ; Semen

Proteomic tools were used to identify the key proteins that might be associated with bronchial asthma(BA). Firstly, the serum samples from healthy adults and asthmatic patients were collected. Tandem Mass Tag~(TM)(TMT), which removes high-abundance structures and nonspecific proteins, was employed to identify the differentially expressed proteins between asthmatic patients and healthy adults. Gene Ontology(GO) annotation and Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment analysis were carried out for the differentially expressed proteins. The core proteins in the asthma group were screened out by protein-protein interaction(PPI) analysis. Then the core proteins were verified by Western blot for 3 patients with bronchial asthma and 3 healthy adults. A total of 778 differentially expressed proteins were screened out, among which 32 proteins contained quantitative information, including 18 up-regulated proteins and 14 down-regulated proteins. The differentially expressed proteins were enriched in 28 KEGG signaling pathways. The PPI analysis showed that 10 proteins(GDN, 1433 Z, VWF, HEMO, CERU, A1 AT, TSP1, G3 P, IBP7, and KPYM) might be involved in the pathogenesis of bronchial asthma. Compared with those in healthy adults, the expression levels of SLC25 A4, SVEP1, and KRT25 in the sera of asthmatic patients were up-regulated(P<0.05). Therefore, it is hypothesized that a variety of immune signaling pathways and differentially expressed proteins play a role in the pathogenesis of BA, which provides potential target information for the treatment of BA.
Adult ; Humans ; Proteomics ; Gene Ontology ; Proteins ; Disease Susceptibility ; Asthma/genetics*

OBJECTIVE: To screen differentially expressed gene (DEG) related to myelodysplastic syndrome (MDS) based on Gene Expression Omnibus (GEO) database, and explore the core genes and pathogenesis of MDS by analyzing the biological functions and related signaling pathways of DEG. METHODS: The expression profiles of GSE4619, GSE19429, GSE58831 including MDS patients and normal controls were downloaded from GEO database. The gene expression analysis tool (GEO2R) of GEO database was used to screen DEG according to | log FC (fold change) |≥1 and P<0.01. David online database was used to annotate gene ontology function (GO). Metascape online database was used to enrich and analyze differential genes in Kyoto Encyclopedia of Genes and Genomes (KEGG). The protein-protein interaction network (PPI) was constructed by using STRING database. CytoHubba and Mcode plug-ins of Cytoscape were used to analyze the key gene clusters and hub genes. R language was used to diagnose hub genes and draw the ROC curve. GSEA enrichment analysis was performed on GSE19429 according to the expression of LEF1. RESULTS: A total of 74 co-DEG were identified, including 14 up-regulated genes and 60 down regulated genes. GO enrichment analysis indicated that BP of down regulated genes was mainly enriched in the transcription and regulation of RNA polymerase II promoter, negative regulation of cell proliferation, and immune response. CC of down regulated genes was mainly enriched in the nucleus, transcription factor complexes, and adhesion spots. MF was mainly enriched in protein binding, DNA binding, and β-catenin binding. KEGG pathway was enriched in primary immunodeficiency, Hippo signaling pathway, cAMP signaling pathway, transcriptional mis-regulation in cancer and hematopoietic cell lineage. BP of up-regulated genes was mainly enriched in type I interferon signaling pathway and viral response. CC was mainly enriched in cytoplasm. MF was mainly enriched in RNA binding. Ten hub genes and three important gene clusters were screened by STRING database and Cytoscape software. The functions of the three key gene clusters were closely related to immune regulation. ROC analysis showed that the hub genes had a good diagnostic significance for MDS. GSEA analysis indicated that LEF1 may affect the normal function of hematopoietic stem cells by regulating inflammatory reaction, which further revealed the pathogenesis of MDS. CONCLUSION Bioinformatics can effectively screen the core genes and key signaling pathways of MDS, which provides a new strategy for the diagnosis and treatment of MDS.
Computational Biology ; Gene Expression Profiling ; Gene Expression Regulation, Neoplastic ; Gene Ontology ; Humans ; Myelodysplastic Syndromes/genetics*