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

OBJECTIVES: Polycystic ovary syndrome (PCOS) is a common endocrine disorder that affects women's health. This study aims to investigate gene and transcription factor (TF) expression differences between PCOS patients and healthy individuals using bioinformatics approaches, and to verify the function of key transcription factors, with the goal of providing new insights into the pathogenesis of PCOS. METHODS: Differentially expressed genes (DEGs) and differentially expressed transcription factors (DETFs) between PCOS patients and controls were identified from the RNA sequencing dataset GSE168404 using bioinformatics methods. Functional enrichment analysis was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The expression and function of core transcription factors were further validated in ovarian tissues of PCOS model mice and control mice using Western blotting and reverse transcription quantitative polymerase chain reaction (RT-qPCR). RESULTS: A total of 332 DEGs were identified between PCOS patients and controls, including 259 upregulated and 73 downregulated genes in the PCOS group. 19 DETFs were further screened, of which 16 were upregulated and 3 were downregulated in PCOS. The upregulated DETFs (including TFCP2L1, DACH1, ESR2, AFF3, SMAD9, ZNF331, HOPX,ATOH8, HIF3α, DPF3, HOXC4, HES1, ID1, JDP2, SOX4, and ID3) were primarily associated with lipid metabolism, development, and cell adhesion. Protein and mRNA expression analysis in PCOS model mice revealed significantly decreased levels of hypoxia-inducible factor (HIF) 1α and HIF2α, and significantly increased expression of HIF3α compared to control mice (all P<0.001). CONCLUSIONS Significant differences in gene and TF expression exist between PCOS patients and healthy individuals. HIF-3α may play a crucial role in PCOS and could serve as a novel biomarker for diagnosis and a potential therapeutic target.
Polycystic Ovary Syndrome/metabolism* ; Female ; Humans ; Animals ; Mice ; Transcription Factors/metabolism* ; Computational Biology ; Gene Expression Profiling ; Adult

OBJECTIVES: Multiple myeloma (MM) is a hematologic malignancy characterized by clonal proliferation of plasma cells and remains incurable. Patients with primary refractory multiple myeloma (PRMM) show poor response to initial induction therapy. This study aims to develop a machine learning-based model to predict treatment response in newly diagnosed multiple myeloma (NDMM) patients, in order to optimize therapeutic strategies. METHODS: NDMM and post-treatment MM patients hospitalized in the Department of Hematology, Third Xiangya Hospital, Central South University, between August 2022 and July 2023 were enrolled. Post-treatment MM patients were categorized into PRMM patients and treatment-responsive MM (TRMM) patients based on therapeutic efficacy. Serum metabolites were detected and analyzed via metabolomics. Based on the metabolomics analysis results and combined with transcriptomic sequencing data of NDMM patients from databases, differentially expressed amino acid metabolism-related genes (AAMGs) among post-treatment NDMM patients with varying therapeutic outcomes were screened. Using bioinformatics analyses and machine learning algorithms, a predictive model for treatment response in NDMM was constructed and used to identify patients at risk for PRMM. RESULTS: A total of 61 patients were included: 22 NDMM, 23 TRMM, and 16 PRMM patients. Significant differences in metabolite levels were observed among the 3 groups, with differential metabolites mainly enriched in amino acid metabolism pathways. Follow-up data were available for 16 of the 22 NDMM patients, including 12 treatment responders (ND_TR group) and 4 with PRMM (ND_PR group). A total of 23 differential metabolites were identified between these 2 groups: 6 metabolites (e.g., tryptophan) were upregulated and 17 (e.g., citric acid) were downregulated in the ND_TR group. Transcriptomic data from 108 TRMM and 77 PRMM patients were analyzed to identify differentially expressed AAMGs, which were then used to construct a prediction model. The area under the receiver operating characteristic curve (AUC) for the model exceeded 0.8, and AUC values in 3 external validation cohorts were all above 0.7. CONCLUSIONS This study delineated the metabolic alterations in MM patients with different treatment response, suggesting that dysregulated amino acid metabolism may be associated with poor treatment response in PRMM. By integrating metabolomics and transcriptomics, a machine learning-based predictive model was successfully established to forecast treatment response in NDMM patients.
Humans ; Multiple Myeloma/drug therapy* ; Machine Learning ; Male ; Female ; Metabolomics/methods* ; Middle Aged ; Aged ; Treatment Outcome ; Transcriptome ; Computational Biology ; Adult ; Multiomics

Dent disease is a rare X-linked recessive inherited renal tubular disorder characterized by low molecular weight proteinuria (LMWP), hypercalciuria, nephrocalcinosis, and other clinical features, and can lead to progressive renal failure. It is primarily caused by mutations in the CLCN5 gene. This article reports the case of a 10-year-old male patient of Chinese descent who was incidentally found to have asymptomatic proteinuria during a routine health examination. Comprehensive biochemical testing and clinical evaluation revealed significant LMWP and hypercalciuria, while renal biopsy showed mesangial cell and matrix proliferation. Whole exome sequencing identified a novel deletion mutation in the CLCN5 gene (NM_001127899.4, c.1158delC, p.F387Lfs*42) causing a frameshift and premature termination, which is likely to disrupt its role in chloride/hydrogen ion exchange and endosomal acidification. Bioinformatic analysis indicated the variant is pathogenic. Genetic testing plays an important role in diagnosing rare kidney diseases. Early identification of pathogenic mutations is essential for facilitating timely intervention and appropriate management, potentially enhancing patient outcomes. This report expands the CLCN5 mutation spectrum and contributes to understanding the genetic and molecular mechanisms of Dent disease.
Humans ; Chloride Channels/genetics* ; Dent Disease/genetics* ; Male ; Child ; Computational Biology ; Mutation ; Proteinuria/genetics* ; Hypercalciuria/genetics*

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 potential molecular targets of quercetin in the treatment of clear cell renal carcinoma (ccRCC). METHODS: The therapeutic targets of quercetin were screened from multiple databases by network pharmacology analysis, and the targets significantly correlated with ccRCC were screened from 4907 plasma proteins using a Mendelian randomization method. The drug-disease network model was constructed to screen the potential key targets. The functions of these targets were evaluated via bioinformatics analysis, and the screened targets were verified in cultured ccRCC cells. RESULTS: Network pharmacology analysis combined with Mendelian randomization identified TP53 (OR=3.325, 95% CI: 1.805-6.124, P=0.0001), ARF4 (OR=0.173, 95% CI: 0.065-0.456, P=0.0003), and DPP4 (OR=0.463, 95% CI: 0.302-0.711, P=0.0004) as the core targets in quercetin treatment of ccRCC. Bioinformatics analysis showed that TP53 was highly expressed in ccRCC, and patients with high TP53 expressions had worse survival outcomes. Molecular docking studies showed that the binding energy between quercetin and TP53 was -5.83 kcal/mol. In cultured 786-O cells, CCK-8 assay and wound healing assay showed that treatment with quercetin significantly inhibited cell proliferation and migration. Quercetin treatment also strongly suppressed the expression of TP53 at both the mRNA and protein levels in 786-O cells as shown by RT-qPCR and Western blotting. CONCLUSIONS TP53 may be the key target of quercetin in the treatment of ccRCC, which sheds light on potential molecular mechanism that mediate the therapeutic effect of quercetin.
Humans ; Quercetin/pharmacology* ; Carcinoma, Renal Cell/genetics* ; Cell Proliferation/drug effects* ; Kidney Neoplasms/genetics* ; Cell Movement/drug effects* ; Tumor Suppressor Protein p53/metabolism* ; Cell Line, Tumor ; Computational Biology

OBJECTIVES: To explore the role of oxidative stress and immune infiltration in rheumatoid arthritis (RA). METHODS: RA datasets GSE55235 (10 RA vs 10 normal samples) and GSE55457 (13 RA vs 10 normal samples) from the GEO database were merged as the test set to identify the differentially expressed genes (DEGs) in RA using R. The DEGs were intersected with oxidative stress-related genes to obtain oxidative stress-associated DEGs. KEGG and GO enrichment analyses of the DEGs were performed, and the RA-related pathways and biological processes were analyzed using GSEA. A protein-protein interaction (PPI) network was constructed using STRING and Cytoscape, and the top 10 key genes were obtained using the Degree algorithm. The validation dataset GSE1919 from GEO database was used for ROC analysis of the key genes to obtain the core genes, and their correlations with infiltrating immune cells were analyzed using CIBERSORT. The results were verified by RT-qPCR for detecting expression levels of the core genes in RA and normal joint samples. RESULTS: We identified 89 oxidative stress-associated DEGs. Enrichment analysis suggested that these DEGs were involved in the biological processes including oxidative stress, chemical stress response, reactive oxygen species response, and lipopolysaccharide response. ROC analysis showed that the 5 core genes (STAT1, MMP9, MYC, CCL5, and JUN) all had AUC values >0.7, indicating their high diagnostic sensitivity and specificity for RA. These genes were closely correlated with immune cells, particularly T cells. RT-qPCR confirmed significant differential expressions of the core genes between RA and normal samples. CONCLUSIONS Oxidative stress and diverse immune responses are features of RA, and the immune responses contribute to activation of oxidative stress. The identified core genes can potential serve as new diagnostic markers for RA.
Arthritis, Rheumatoid/genetics* ; Oxidative Stress/genetics* ; Humans ; Computational Biology ; Protein Interaction Maps ; Gene Expression Profiling ; Gene Regulatory Networks

OBJECTIVES: To explore the potential of pyrroline-5-carboxylate reductase 1 (PYCR1) as a pan-cancer biomarker and investigate its expression, function, and clinical significance in bladder cancer (BLCA). METHODS: Bioinformatics analysis was conducted to evaluate the associations of PYCR1 with prognosis, immune microenvironment remodeling, tumor mutation burden (TMB), and microsatellite instability (MSI) in cancer patients. Using the TCGA-BLCA dataset, univariate and multivariate regression analyses were performed to assess the potential of PYCR1 as an independent prognostic risk factor for BLCA, and a clinical decision model was constructed. The IMvigor210 cohort was utilized to evaluate the potential of PYCR1 for independently predicting the efficacy of immunotherapy. The pRRophetic was employed to screen candidate chemotherapeutic agents for treating BLCA with high PYCR1 expression. The CMap-XSum algorithm and molecular docking techniques were used to explore and validate small molecule inhibitors of PYCR1. RESULTS: A high expression of PYCR1 was significantly associated with poor prognosis, immune cell infiltration, TMB and MSI in various tumors (r>0.3). PYCR1 was overexpressed in BLCA, and high PYCR1 expression was closely related to poor prognosis in BLCA patients (HR: 1.14, 95% CI: 1.02-1.68, P=0.006). The IC₅₀ of the anti-cancer drugs cetuximab, 5-fluorouracil, and doxorubicin increased significantly in BLCA cell lines with high PYCR1 expressions (P<0.0001). CONCLUSIONS High PYCR1 expression is an independent risk factor for poor prognosis in BLCA patients and can serve as a significant indicator for clinical decision-making as well as a marker for predicting sensitivity to chemotherapeutic agents and the efficacy of immunotherapy.
Humans ; Urinary Bladder Neoplasms/genetics* ; Immunotherapy ; Prognosis ; Pyrroline Carboxylate Reductases/metabolism* ; Biomarkers, Tumor/genetics* ; delta-1-Pyrroline-5-Carboxylate Reductase ; Microsatellite Instability ; Tumor Microenvironment ; Mutation ; Computational Biology ; Molecular Docking Simulation

OBJECTIVES: To establish a prognostic model for primary liver cancer (PLC) using bioinformatics methods. METHODS: Based on the data from 404 patients in the Cancer Genome Atlas (TCGA) database, we constructed a prognostic model integrating the differentially expressed genes, anoikis, and immune-related genes (DAIs) using univariate Cox regression and the LASSO-Cox approach. The predictive ability of the model was evaluated using Kaplan-Meier method and receiver-operating characteristic curves, and a nomogram was developed to facilitate its clinical applications. Gene set enrichment analysis (GSEA) was performed to explore the associated pathways and relationship between the DAIs and the tumor immune microenvironment, and the half-maximal inhibitory concentration (IC₅₀) of liver cancer drugs was calculated using the "pRRophetic" R package. We also detected the expression of SEMA7A in paired tumor and adjacent tissues from liver cancer patients. RESULTS: We constructed and validated a prognostic model based on 7 DAIs (NR4A3, SEMA7A, IL11, AR, BIRC5, EGF, and SPP1), and obtained consistent results in both the TCGA training cohort and GEO validation cohort (GSE14520), where the patients in the low-risk group were characterized by more favorable clinical outcomes and immune status. By integrating this prognostic signature with clinical information, a composite nomogram was generated. Somatic mutation analysis showed that TTN, TP53, and CTNNB1 mutations accounted for the largest proportion of total mutations, and the patients in the low-risk-low-TMB group had higher survival rate. Drug sensitivity analysis revealed differences in sensitivity to chemotherapeutic agents between high- and low-risk groups and between TP53 mutations and non-mutations. In clinical tissue specimens, SEMA7A expression was significantly higher in liver cancer tissues than in the adjacent tissues. CONCLUSIONS We established a new prognostic model based on DAIs for predicting clinical outcomes and therapeutic response of patients with primary liver cancer.
Humans ; Liver Neoplasms/diagnosis* ; Prognosis ; Anoikis ; Nomograms ; Computational Biology ; Tumor Microenvironment ; Semaphorins/metabolism*