Comparative analysis of gene expression between mice and humans in acetaminophen-induced liver injury by integrating bioinformatics analysis

Abstract

Abstract Objective Mice are routinely utilized as animal models of drug-induced liver injury (DILI), however, there are significant differences in the pathogenesis between mice and humans. This study aimed to compare gene expression between humans and mice in acetaminophen (APAP)-induced liver injury (AILI), and investigate the similarities and differences in biological processes between the two species. Methods A pair of public datasets (GSE218879 and GSE120652) obtained from GEO were analyzed using “Limma” package in R language, and differentially expressed genes (DEGs) were identified, including co-expressed DEGs (co-DEGs) and specific-expressed DEGS (specific-DEGs). Analysis of Gene Set Enrichment Analysis (GSEA), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed analyses for specific-DEGs and co-DEGs. The co-DEGs were also used to construct transcription factor (TF)-gene network, gene-miRNA interactions network and protein-protein interaction (PPI) network for analyzing hub genes. Results Mouse samples contained 1052 up-regulated genes and 1064 down-regulated genes, while human samples contained 1156 up-regulated genes and 1557 down-regulated genes. After taking the intersection between the DEGs, only 154 co-down-regulated and 89 co-up-regulated DEGs were identified, with a proportion of less than 10%. It was suggested that significant differences in gene expression between mice and humans in drug-induced liver injury. Mouse-specific-DEGs predominantly engaged in processes related to apoptosis and endoplasmic reticulum stress, while human-specific-DEGs were concentrated around catabolic process. Analysis of co-regulated genes reveals showed that they were mainly enriched in biosynthetic and metabolism-related processes. Then a PPI network which contains 189 nodes and 380 edges was constructed from the co-DEGs and two modules were obtained by Mcode. We screened out 10 hub genes by three algorithms of Degree, MCC and MNC, including CYP7A1, LSS, SREBF1, FASN, CD44, SPP1, ITGAV, ANXA5, LGALS3 and PDGFRA. Besides, TFs such as FOXC1, HINFP, NFKB1, miRNAs like mir-744-5p, mir-335-5p, mir-149-3p, mir-218-5p, mir-10a-5p may be the key regulatory factors of hub genes. Conclusions The DEGs of AILI mice models and those of patients were compared, and common biological processes were identified. The signaling pathways and hub genes in co-expression were identified between mice and humans through a series of bioinformatics analyses, which may be more valuable to reveal molecular mechanisms of AILI.