Identification of key genes, pathways and potential therapeutic agents for liver fibrosis using an integrated bioinformatics analysis

Background Liver fibrosis is often a consequence of chronic liver injury, and has the potential to progress to cirrhosis and liver cancer. Despite being an important human disease, there are currently no approved anti-fibrotic drugs. In this study, we aim to identify the key genes and pathways governing the pathophysiological processes of liver fibrosis, and to screen therapeutic anti-fibrotic agents. Methods Expression profiles were downloaded from the Gene Expression Omnibus (GEO), and differentially expressed genes (DEGs) were identified by R packages (Affy and limma). Gene functional enrichments of each dataset were performed on the DAVID database. Protein–protein interaction (PPI) network was constructed by STRING database and visualized in Cytoscape software. The hub genes were explored by the CytoHubba plugin app and validated in another GEO dataset and in a liver fibrosis cell model by quantitative real-time PCR assay. The Connectivity Map L1000 platform was used to identify potential anti-fibrotic agents. Results We integrated three fibrosis datasets of different disease etiologies, incorporating a total of 70 severe (F3–F4) and 116 mild (F0–F1) fibrotic tissue samples. Gene functional enrichment analyses revealed that cell cycle was a pathway uniquely enriched in a dataset from those patients infected by hepatitis B virus (HBV), while the immune-inflammatory response was enriched in both the HBV and hepatitis C virus (HCV) datasets, but not in the nonalcoholic fatty liver disease (NAFLD) dataset. There was overlap between these three datasets; 185 total shared DEGs that were enriched for pathways associated with extracellular matrix constitution, platelet-derived growth-factor binding, protein digestion and absorption, focal adhesion, and PI3K-Akt signaling. In the PPI network, 25 hub genes were extracted and deemed to be essential genes for fibrogenesis, and the expression trends were consistent with GSE14323 (an additional dataset) and liver fibrosis cell model, confirming the relevance of our findings. Among the 10 best matching anti-fibrotic agents, Zosuquidar and its corresponding gene target ABCB1 might be a novel anti-fibrotic agent or therapeutic target, but further work will be needed to verify its utility. Conclusions Through this bioinformatics analysis, we identified that cell cycle is a pathway uniquely enriched in HBV related dataset and immune-inflammatory response is clearly enriched in the virus-related datasets. Zosuquidar and ABCB1 might be a novel anti-fibrotic agent or target.