A number of recent studies have highlighted the findings that certain lncRNAs are associated with alternative splicing (AS) in tumorigenesis and progression. Although existing work showed the importance of linking certain misregulations of RNA splicing with lncRNAs, a primary concern is the lack of genome-wide comprehensive analysis for their associations.We analyzed an extensive collection of RNA-seq data, quantified 198,619 isoform expressions, and found systematic isoform usage changes between hepatocellular carcinoma (HCC) and normal liver tissue. We identified a total of 1375 splicing switched isoforms and further analyzed their biological functions.To predict which lncRNAs are associated with these AS genes, we integrated the coexpression networks and epigenetic interaction networks collected from text mining and database searching, linking lncRNA modulators such as splicing factors, transcript factors, and miRNAs with their targeted AS genes in HCC. To model the heterogeneous networks in a single framework, we developed a multi-graphic random walk (RWMG) network method to prioritize the lncRNAs associated with AS in HCC. RWMG showed a good performace evaluated by ROC curve based on cross-validation and bootstrapping strategy.As a summary, we identified 31 AS-related lncRNAs including MALAT1 and HOXA11-AS, which have been reported before, as well as some novel lncRNAs such as DNM1P35, HAND2-AS1, and DLX6-AS1. Survival analysis further confirmed the clinical significance of identified lncRNAs.Keywords: lncRNA, alternative splicing, random walk, network analysis in tumorigeneses. Therefore, a key open question is: how many novel lncRNAs are associated with AS modulations in tumorigenesis, genome-wisely.Here we utilized a novel network propagation technology, random walk-based multigraphic model (RWMG), to simultaneously integrate complicated biological connections among lncRNA -effectors (TF, ASF, and miRNA) -AS interaction networks and co-expression networks in a single analysis framework. This method is an extended application inspired by Random walk with restart algorithm to prioritize important lncRNAs that are involved in AS based on the hypothesis that more important genes are likely to receive more links from other networks. In comparison with traditional random walk algorithms, which treat all genes equally, our flexible, scalable method can be formulated to rank a subset of vertices (e.g., PCGs, lncRNAs), based on pre-knowledge as the starting walking vertices. This method is more accurate than other traditional "shortest path" network-based integrative methods, as it can overcome the "noisy" and "incomplete" highly dimensional heterogeneous data.In addition, previous tumor and normal comparison studies are limited to normal adjacent to tumor (NAT) tissues. However, these tissues are not truly 'normal' as they usually surrounded by tumor contaminations. Therefore, many potential cancer biomarkers involved in AS may be missed. By combining the 'pure' healthy liver organs from GTEx with TCGA expression...