Abstract
BACKGROUND: Functions of miRNAs involved in tumorigenesis are well reported, yet, their roles in normal cell lineage commitment remain ambiguous. Here, we investigated a specific "transcription factor (TF)-miRNA-Target" regulatory network during the lineage maturation of biliary tree stem cells (BTSCs) into adult hepatocytes (hAHeps). METHOD: Bioinformatic analysis was conducted based on our RNA-seq and microRNA-seq datasets with four human hepatic-lineage cell lines, including hBTSCs, hepatic stem cells (hHpSCs), hepatoblasts (hHBs), and hAHeps. Short time-series expression miner (STEM) analysis was performed to reveal the time-dependent dynamically changed miRNAs and mRNAs. GO and KEGG analyses were applied to reveal the potential function of key miRNAs and mRNAs. Then, the miRDB, miRTarBase, TargetScan, miRWalk, and DIANA-microT-CDS databases were adopted to predict the potential targets of miRNAs while the TransmiR v2.0 database was used to obtain the experimentally supported TFs that regulate miRNAs. The TCGA, Kaplan-Meier Plotter, and human protein atlas (HPA) databases and more than 10 sequencing data, including bulk RNA-seq, microRNA-seq, and scRNA-seq data related to hepatic development or lineage reprogramming, were obtained from both our or other published studies for validation. RESULTS: STEM analysis showed that during the maturation from hBTSCs to hAHeps, 52 miRNAs were downwardly expressed and 928 mRNA were upwardly expressed. Enrichment analyses revealed that those 52 miRNAs acted as pluripotency regulators for stem cells and participated in various novel signaling pathways, including PI3K/AKT, MAPK, and etc., while 928 mRNAs played important roles in liver-functional metabolism. With an extensive sorting of those key miRNAs and mRNAs based on the target prediction results, 23 genes were obtained which not only functioned as the targets of 17 miRNAs but were considered critical for the hepatic lineage commitment. A "TF-miRNA-Target" regulatory network for hepatic lineage commitment was therefore established and had been well validated by various datasets. The network revealed that the PI3K/AKT pathway was gradually suppressed during the hepatic commitment. CONCLUSION: A total of 17 miRNAs act as suppressors during hepatic maturation mainly by regulating 23 targets and modulating the PI3K/AKT signaling pathway. The regulatory network uncovers possible signatures and guidelines enabling us to identify or obtain the functional hepatocytes for future study.