Revealing immune infiltrate characteristics and potential immune-related genes in hepatic fibrosis: based on bioinformatics, transcriptomics and q-PCR experiments

Background The occurrence and progression of hepatic fibrosis (HF) is accompanied by inflammatory damage. Immune genes play a pivotal role in fibrogenesis and inflammatory damage in HF by regulating immune cell infiltration. However, the immune mechanisms of HF are inadequately studied. Therefore, this research aims to identify the immune genes and biological pathway which involved in fibrosis formation and inflammatory damage in HF and explore immune target-based therapeutics for HF. Methods The expression dataset GSE84044 of HF was downloaded from the GEO database. The crucial module genes for HF were screened according to weighted gene co-expression network analysis (WGCNA). The crucial module genes were mapped to immune-related genes obtained from the ImmPort database to obtain the hepatic fibrosis immune genes (HFIGs). In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses were performed on HFIGs. Then, the protein-protein interaction (PPI) network was conducted on HFIGs and hub genes were identified from the PPI network. Moreover, immune infiltration analysis was performed to identified correlation between hub gene and immune cell infiltration. To verify the reliability of the GSE84044 expression profile data analysis, a rat model of CCl4-induced HF was established, followed by transcriptome sequencing and immunofluorescence analysis and quantitative reverse transcription (q-PCR) experiments were performed in HF rats and normal rat liver tissues. Finally, CMAP platform was used to explore immune target-based therapeutics for HF. Results In the bioinformatics analysis of GSE84044 data, 98 HFIGs were screened. These genes were mainly involved in inflammation-related biological pathways such as NOD-like receptor signaling pathway, NF-kappa B signaling pathway, Toll-like receptor signaling pathway and PI3K-Akt signaling pathway. From the PPI network, 10 hub genes were identified, including CXCL8, IL18, CXCL10, CD8A, IL7, PTPRC, CCL5, IL7R, CXCL9 and CCL2. Immune infiltration analysis showed that immune cells like neutrophils, natural killer (NK) cells, macrophages M1 and macrophages M2 were significantly correlated with the hepatic fibrosis process and hub gene expression was significantly correlated with these immune cells. Notably, most of the biological pathways HFIGs riched and all the hub gene expression except CXCL8 were validated in subsequent transcriptome and qRCR experiments. Finally, 15 small molecule compounds with the potential to reverse the high expression of hub genes were screen out as potential therapeutic agents for HF. Conclusion The immune genes CXCL8, IL18, CXCL10, CD8A, IL7, PTPRC, CCL5, IL7R, CXCL9 and CCL2 may play an essential role in the fibrosis formation and inflammatory damage in HF. The outcomes of this research provide a basis for the study of the immune mechanisms of HF and contribute to the diagnosis and prevention and treatment of HF in clinical practice.