Identification of prognostic and metastasis-related alternative splicing signatures in hepatocellular carcinoma

As the most common neoplasm in digestive system, hepatocellular carcinoma (HCC) is one of the most important leading cause of cancer deaths worldwide. Its high-frequency metastasis and relapse rate lead to the poor survival of HCC patients. However, the mechanism of HCC metastasis is still unclear. Alternative splicing events (ASEs) have a great effect in cancer development, progression and metastasis. We downloaded RNA sequencing and seven types of ASEs data of HCC samples, in order to explore the mechanism of ASEs underlying tumorigenesis and metastasis of HCC. The data were taken from the The Cancer Genome Atlas (TCGA) and TCGASpliceSeq databases. Univariate Cox regression analysis was used to determine a total of 3197 overall survival-related ASEs (OS-SEs). And based on five OS-SEs screened by Lasso regression, we constructed a prediction model with the Area Under Curve of 0.765. With a good reliability of the model, the risk score was also proved to be an independent predictor. Among identified 390 candidate SFs, Y-box protein 3 (YBX3) was significantly correlated with OS and metastasis. Among 177 ASEs, ATP-binding cassette subfamily A member 6 (ABCA6)-43162-AT and PLIN5-46808-AT were identified both associated with OS, bone metastasis and co-expressed with SFs. Then we identified primary bile acid biosynthesis as survival-related (KEGG) pathway by Gene Set Variation Analysis (GSVA) and univariate regression analysis, which was correlated with ABCA6-43162-AT and PLIN5-46808-AT. Finally, we proposed that ABCA6-43162-AT and PLIN5-46808-AT may contribute to HCC poor prognosis and metastasis under the regulation of aberrant YBX3 through the pathway of primary bile acid biosynthesis.