The log2 fold change value of each gene according to RRA result and TCGA-HCC dataset were present in Figure 3E and Figure S1B. The predictive risk model based on the six hub genes expression value and multivariate Cox coefficients was constructed as follow: risk score = (-2.17× the expression value of CD163) + (-1.37× the expression value of EHHADH) + (1.28× the expression value of KIAA0101) + (-0.96× the expression value of SLC16A2) + (1.46× the expression value of SPP1) + (1.57× the expression value of THBS4). Here, THBS4 is linked to hepatocellular carcinoma.