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*

This study employed bioinformatics to screen the feature genes related to efferocytosis in diabetic kidney disease(DKD) and explores traditional Chinese medicine(TCM) regulating these feature genes. The GSE96804 and GSE30528 datasets were integrated as the training set, and the intersection of differentially expressed genes and efferocytosis-related genes(ERGs) was identified as DKD-ERGs. Subsequently, correlation analysis, protein-protein interaction(PPI) network construction, enrichment analysis, and immune infiltration analysis were performed. Consensus clustering was conducted on DKD patients based on the expression levels of DKD-ERGs, and the expression levels, immune infiltration characteristics, and gene set variations between different subtypes were explored. Eight machine learning models were constructed and their prediction performance was evaluated. The best-performing model was evaluated by nomograms, calibration curves, and external datasets, followed by the identification of efferocytosis-related feature genes associated with DKD. Finally, potential TCMs that can regulate these feature genes were predicted. The results showed that the training set contained 640 differentially expressed genes, and after intersecting with ERGs, 12 DKD-ERGs were obtained, which demonstrated mutual regulation and immune modulation effects. Consensus clustering divided DKD into two subtypes, C1 and C2. The support vector machine(SVM) model had the best performance, predicting that growth arrest-specific protein 6(GAS6), S100 calcium-binding protein A9(S100A9), C-X3-C motif chemokine ligand 1(CX3CL1), 5'-nucleotidase(NT5E), and interleukin 33(IL33) were the feature genes of DKD. Potential TCMs with therapeutic effects included Astragali Radix, Trionycis Carapax, Sargassum, Rhei Radix et Rhizoma, Curcumae Radix, and Alismatis Rhizoma, which mainly function to clear heat, replenish deficiency, activate blood, resolve stasis, and promote urination and drain dampness. Molecular docking revealed that the key components of these TCMs, including β-sitosterol, quercetin, and sitosterol, exhibited good binding activity with the five target genes. These results indicated that efferocytosis played a crucial role in the development and progression of DKD. The feature genes closely related to both DKD and efferocytosis, such as GAS6, S100A9, CX3CL1, NT5E, and IL33, were identified. TCMs such as Astragali Radix, Trionycis Carapa, Sargassum, Rhei Radix et Rhizoma, Curcumae Radix, and Alismatis Rhizoma may provide a new therapeutic strategy for DKD by regulating efferocytosis.
Humans ; Computational Biology ; Diabetic Nephropathies/physiopathology* ; Protein Interaction Maps ; Medicine, Chinese Traditional ; Drugs, Chinese Herbal ; Phagocytosis/genetics* ; Efferocytosis

BACKGROUND: Diabetic foot is a complex condition with high incidence, recurrence, mortality, and disability rates. Current treatments for diabetic foot ulcers are often insufficient. This study was conducted to identify potential therapeutic targets for diabetic foot. METHODS: Datasets related to diabetic foot and diabetic skin were retrieved from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified using R software. Enrichment analysis was conducted to screen for critical gene functions and pathways. A protein interaction network was constructed to identify node genes corresponding to key proteins. The DEGs and node genes were overlapped to pinpoint target genes. Plasma and chronic ulcer samples from diabetic and non-diabetic individuals were collected. Western blotting, immunohistochemistry, and enzyme-linked immunosorbent assays were performed to verify the S100 calcium binding protein A9 (S100A9), inflammatory cytokine, and related pathway protein levels. Hematoxylin and eosin staining was used to measure epidermal layer thickness. RESULTS: In total, 283 common DEGs and 42 node genes in diabetic foot ulcers were identified. Forty-three genes were differentially expressed in the skin of diabetic and non-diabetic individuals. The overlapping of the most significant DEGs and node genes led to the identification of S100A9 as a target gene. The S100A9 level was significantly higher in diabetic than in non-diabetic plasma (178.40 ± 44.65 ng/mL vs. 40.84 ± 18.86 ng/mL) and in chronic ulcers, and the wound healing time correlated positively with the plasma S100A9 level. The levels of inflammatory cytokines (tumor necrosis factor-α, interleukin [IL]-1, and IL-6) and related pathway proteins (phospho-extracellular signal regulated kinase [ERK], phospho-p38, phospho-p65, and p-protein kinase B [Akt]) were also elevated. The epidermal layer was notably thinner in chronic diabetic ulcers than in non-diabetic skin (24.17 ± 25.60 μm vs. 412.00 ± 181.60 μm). CONCLUSIONS S100A9 was significantly upregulated in diabetic foot and was associated with prolonged wound healing. S100A9 may impair diabetic wound healing by disrupting local inflammatory responses and skin re-epithelialization.
Calgranulin B/therapeutic use* ; Diabetic Foot/metabolism* ; Humans ; Datasets as Topic ; Computational Biology ; Mice, Inbred C57BL ; Animals ; Mice ; Protein Interaction Maps ; Immunohistochemistry

BACKGROUND: Gluconeogenesis is a critical metabolic pathway for maintaining glucose homeostasis, and its dysregulation can lead to glycometabolic disorders. This study aimed to identify hub biomarkers of these disorders to provide a theoretical foundation for enhancing diagnosis and treatment. METHODS: Gene expression profiles from liver tissues of three well-characterized gluconeogenesis mouse models were analyzed to identify commonly differentially expressed genes (DEGs). Weighted gene co-expression network analysis (WGCNA), machine learning techniques, and diagnostic tests on transcriptome data from publicly available datasets of type 2 diabetes mellitus (T2DM) patients were employed to assess the clinical relevance of these DEGs. Subsequently, we identified hub biomarkers associated with gluconeogenesis-related glycometabolic disorders, investigated potential correlations with immune cell types, and validated expression using quantitative polymerase chain reaction in the mouse models. RESULTS: Only a few common DEGs were observed in gluconeogenesis-related glycometabolic disorders across different contributing factors. However, these DEGs were consistently associated with cytokine regulation and oxidative stress (OS). Enrichment analysis highlighted significant alterations in terms related to cytokines and OS. Importantly, osteomodulin ( OMD ), apolipoprotein A4 ( APOA4 ), and insulin like growth factor binding protein 6 ( IGFBP6 ) were identified with potential clinical significance in T2DM patients. These genes demonstrated robust diagnostic performance in T2DM cohorts and were positively correlated with resting dendritic cells. CONCLUSIONS Gluconeogenesis-related glycometabolic disorders exhibit considerable heterogeneity, yet changes in cytokine regulation and OS are universally present. OMD , APOA4 , and IGFBP6  may serve as hub biomarkers for gluconeogenesis-related glycometabolic disorders.
Machine Learning ; Humans ; Computational Biology/methods* ; Biomarkers/metabolism* ; Diabetes Mellitus, Type 2/genetics* ; Animals ; Mice ; Gluconeogenesis/physiology* ; Gene Expression Profiling ; Transcriptome/genetics* ; Gene Regulatory Networks/genetics* ; Clinical Relevance

With the rapid development of artificial intelligence, computational pathology has been seamlessly integrated into the entire clinical workflow, which encompasses diagnosis, treatment, prognosis, and biomarker discovery. This integration has significantly enhanced clinical accuracy and efficiency while reducing the workload for clinicians. Traditionally, research in this field has depended on the collection and labeling of large datasets for specific tasks, followed by the development of task-specific computational pathology models. However, this approach is labor intensive and does not scale efficiently for open-set identification or rare diseases. Given the diversity of clinical tasks, training individual models from scratch to address the whole spectrum of clinical tasks in the pathology workflow is impractical, which highlights the urgent need to transition from task-specific models to foundation models (FMs). In recent years, pathological FMs have proliferated. These FMs can be classified into three categories, namely, pathology image FMs, pathology image-text FMs, and pathology image-gene FMs, each of which results in distinct functionalities and application scenarios. This review provides an overview of the latest research advancements in pathological FMs, with a particular emphasis on their applications in oncology. The key challenges and opportunities presented by pathological FMs in precision oncology are also explored.
Humans ; Precision Medicine/methods* ; Medical Oncology/methods* ; Artificial Intelligence ; Neoplasms/pathology* ; Computational Biology/methods*

This study aims to investigate the therapeutic effect and mechanism of Daotan Xixin Decoction on APP/PS1 mice. Twelve APP/PS1 male mice were randomized into four groups: APP/PS1 and low-, medium-, and high-dose Daotan Xixin Decoction. Three C57BL/6 wild-type mice were used as the control group. The learning and memory abilities of mice in each group were examined by the Morris water maze test. The pathological changes of hippocampal nerve cells were observed by hematoxylin-eosin staining and Nissl staining. Immunohistochemistry was employed to detect the expression of β-amyloid(Aβ)_(1-42) in the hippocampal tissue. The high-dose Daotan Xixin Decoction group with significant therapeutic effects and the model group were selected for high-throughput sequencing. The differentially expressed gene(DEG) analysis, Gene Ontology(GO) analysis, Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment analysis, and Gene Set Variation Analysis(GSVA) were performed on the sequencing results. RT-qPCR and Western blot were conducted to determine the mRNA and protein levels, respectively, of some DEGs. Compared with the APP/PS1 group, Daotan Xixin Decoction at different doses significantly improved the learning and memory abilities of APP/PS1 mice, ameliorated the neuropathological damage in the CA1 region of the hippocampus, increased the number of neurons, and decreased the deposition of Aβ_(1-42) in the brain. A total of 1 240 DEGs were screened out, including 634 genes with up-regulated expression and 606 genes with down-regulated expression. The GO analysis predicted the biological processes including RNA splicing and protein folding, the cellular components including spliceosome complexes and nuclear spots, and the molecular functions including unfolded protein binding and heat shock protein binding. The KEGG pathway enrichment analysis revealed the involvement of neurodegenerative disease pathways, amyotrophic lateral sclerosis, and splicing complexes. Further GSVA pathway enrichment analysis showed that the down-regulated pathways involved nuclear factor-κB(NF-κB)-mediated tumor necrosis factor-α(TNF-α) signaling pathway, UV response, and unfolded protein response, while the up-regulated pathways involved the Wnt/β-catenin signaling pathway. The results of RT-qPCR and Western blot showed that compared with the APP/PS1 group, Daotan Xixin Decoction at different doses down-regulated the mRNA and protein levels of signal transducer and activator of transcription 3(STAT3), NF-κB, and interleukin-6(IL-6) in the hippocampus. In conclusion, Daotan Xixin Decoction can improve the learning and memory abilities of APP/PS1 mice by regulating the STAT3/NF-κB/IL-6 signaling pathway.
Animals ; Drugs, Chinese Herbal/administration & dosage* ; Mice ; Male ; Alzheimer Disease/metabolism* ; Computational Biology ; Mice, Inbred C57BL ; High-Throughput Nucleotide Sequencing ; Amyloid beta-Protein Precursor/metabolism* ; Hippocampus/metabolism* ; Mice, Transgenic ; Presenilin-1/metabolism* ; Humans ; Memory/drug effects* ; Maze Learning/drug effects* ; Amyloid beta-Peptides/genetics* ; Disease Models, Animal

With the advancement in big data and artificial intelligence technologies, the extensive application of omics technologies in traditional Chinese medicine(TCM) research has generated large experimental datasets, enabling the exploration of cross-scale correlations among massive data and thereby resulting in the shift toward a data-intensive research paradigm. The emerging approach of multi-omics data fusion analysis, emphasizing technical and computational tools, presents a potential breakthrough in this field. The holistic perspective of TCM aligns with the concept of multi-omics data fusion, yet the data types encountered exhibit high dimensionality with small sample sizes, necessitating data processing techniques such as dimensionality reduction. The current challenge lies in selecting suitable analytical methods for these data to enhance the systematic understanding of physiological functions and disease diagnosis/treatment processes. This paper explores the theories and frameworks of multi-omics data fusion, analyzes methods for fusing high-dimensional, small-sample multi-omics data in TCM, and aims to provide insights for advancing TCM research.
Medicine, Chinese Traditional/methods* ; Humans ; Computational Biology/methods* ; Genomics/methods* ; Sample Size ; Artificial Intelligence ; Multiomics