Abstract
Liver cancer is the fourth leading cause of cancer-related deaths with a steadily increasing rate worldwide, as a well-known hallmark of liver cancer, metabolic alterations are related to liposomal changes, a common characteristic of primary liver cancers based on recent lipidomics studies. Peroxisome proliferator-activated receptor α (PPARα) is a ligand-activated transcription factor with important lipid homeostasis function, therefore we aimed to understand the molecular mechanisms and pathways that activate PPARα after using PPAR-α agonist WY-14643 and identify candidate biomarkers related to PPARα activity and evaluate their effects in liver cancer. The data from differently expressed genes (DEGs) between liver cancer tissue from obese subjects alone and liver tissue after treatment were evaluated by DESeq2 and module genes were analyzed using weighted gene co-expression network analysis (WGCNA). Final candidate genes were identified by intersecting genes among highly ranked DEGs and the brown module, which demonstrated a significant negative correlation with drug treatments. We conducted a protein-protein interaction network, and KEGG enrichment analysis, and core hub genes (CD40, CXCL9, CXCL10, TNFSF14, GBP2, GBP3, APOL3, CLDN1) were identified using the cyto-hubba plugin, among them we focused on GBP2 that plays key roles in oncogenesis and evaluate its expressional with clinical outcomes. In conclusion, the WGCNA-based co-expression network identified GBP2 as one of the hub genes with a negative relation with PPARα agonist treatments. higher expression of GBP2 was closely associated with HCC progression. Therefore, GBP2 might be a potential candidate for the study of PPARα activity in HCC.