Cancer stem cells (CSCs) are heterogeneous with self-renewal and differentiation ability. The mRNA expression-based stemness index (mRNAsi) described the similarity between tumor cells and CSCs, which is positively associated with the poor prognosis of cancer patients. However, the key prognostic genes related to mRNAsi in hepatocellular carcinoma (HCC) remains unclear. A 9-gene signature related to mRNAsi and HCC prognosis including PSMG3, SNRPD1, DTYMK, PIGU, NME1, TXNL4A, IPO4, PES1, and REXO4 was obtained. High expression of this signature indicates poor prognosis of HCC. PIGU was an independent prognostic factor of HCC, which was significantly associated with progression of HCC. Among them, DTYMK and NME1 enriched in pyrimidine metabolism, SNRPD1 and TXNL4A enriched in spliceosome and PIGU enriched in glycosyl phosphatidylinositol (GPI)-anchor biosynthesis pathways. High levels of IPO4, NME1, PES1, PIGU and SNRPD1 were closely associated with metastasis of HCC, and low levels of IPO4, PIGU and REOX4 were significantly associated with sorafenib resistance of HCC. High expression of the 9-gene signature was negatively correlated with the stromal cell infiltration, and positively correlated with specific immune subtypes-related to angiogenesis, M1/M2 macrophage polarization, and M2 response. The 9-gene signature was negatively correlated with the stroma, and SNRPD1 and TXNL4 were positively correlated with immune infiltrate. NME1 was negatively correlated with tumor purity. Therefore, a 9-gene signature related to mRNAsi and poor prognosis in HCC were identified, which can be used as biomarkers for the diagnosis of HCC and functional mechanism exploration of CSCs in HCC. These genes such as IPO4 and PIGU might drive the transition of tumor cells into CSCs which possibly controls the balance between metastasis and drug resistance in HCC. The challenge on balance between metastasis and drug resistance for tumor therapy was firstly reported by the present study.