Functional characteristics of YAP-positive hepatocytes expression in an early stage of NASH with transcriptome sequence analysis

Abstract

Objective: To analyze and compare the differentially expressed genes (DEGs) of Yes-associated protein (YAP)-positive and negative hepatocytes and further understand the preliminary functional characteristics of YAP-positive hepatocytes in an early mouse model of nonalcoholic steatohepatitis (NASH) with transcriptome sequence (RNA-Seq). Methods: C57BL/6 mice were fed with methionine-choline deficiency (MCD) diet for 2 weeks to establish an early NASH model, and the control group was fed with normal diet. Liver tissue was stained with hematoxylin-eosin (HE) and Sirius red, and the pathological score was recorded. The expression of YAP and P-YAP were determined by immunohistochemistry (IHC) in liver tissues. Primary hepatocytes with viability greater than 90% were isolated and purified by collagenase perfusion combined with Percoll density gradient centrifugation. YAP-positive and negative hepatocytes were assessed by YAP antibody, flow cytometry and RNA-Seq analyses. Sequencing results were screened by GO, KEGG and interaction network analysis methods. RT-PCR was used to verify the expression levels of YAP and some DEGs in liver tissue model group. Two samples mean was compared by independent samples t-test. Results: Compared with the control group, the HE-stained liver tissue of MCD-induced mice at 2 weeks showed steatosis (pathological score 1.07±0.21), accompanied by lobular inflammation (pathological score 1.13±0.32) and ballooned hepatocyte (pathological score 0.80) ±0.20). Sirius red staining showed non-significant liver fibrosis (pathological score 0.40±0.40). IHC showed partial YAP-positive hepatocytes expression in an early stage of NASH. RNA-Seq analysis showed that clean reads of YAP-positive and negative hepatocytes were 49 310 604 and 5 4820 036, respectively. Compared with YAP-negative hepatocytes, YAP-positive hepatocytes had differential expression of 5 565 genes, including 1 662 up-regulated genes and 3 903 down-regulated genes. GO analysis of up-regulated genes showed that the metabolic processes related to mitochondrial functions, such as purine nucleoside triphosphate and nucleoside triphosphate were significantly enriched in biological processes (BP), while down-regulated gene analysis showed that olfactory-related receptor were significantly enriched in BP. KEGG analysis showed that DEGs were enriched in 292 pathways, and oxidative phosphorylation (OXPHOS) pathway was significantly enriched in signaling pathway. RT-PCR validated that inflammatory factors (interleukin-1β, interleukin-6), YAP and its target genes (Cyr61, Ankrd1), and Cox5b and Sdhc genes were significantly up-regulated in the OXPHOS pathway, which was consistent with the sequencing results. In addition, eight key genes with interaction network analysis were predicted. Conclusion: Changes in hepatocyte metabolic levels may be associated with increased YAP activity in an early stage of NASH.