lncRNA TMEM51-AS1 and RUSC1-AS1 function as ceRNAs for induction of laryngeal squamous cell carcinoma and prediction of prognosis

BACKGROUND: Long non-coding RNAs (lncRNAs) can function as competing endogenous RNAs (ceRNAs) to interact with miRNAs to regulate target genes and promote cancer initiation and progression. The expression of lncRNAs and miRNAs can be epigenetically regulated. The goal of this study was to construct an lncRNA-miRNA-mRNA ceRNA network in laryngeal squamous cell carcinoma (LSCC) and reveal their methylation patterns, which was not investigated previously. METHODS: Microarray datasets available from the Gene Expression Omnibus database were used to identify differentially expressed lncRNAs (DELs), miRNAs (DEMs), and genes (DEGs) between LSCC and controls, which were then overlapped with differentially methylated regions (DMRs). The ceRNA network was established by screening the interaction relationships between miRNAs and lncRNAs/mRNAs by corresponding databases. TCGA database was used to identify prognostic biomarkers. RESULTS: Five DELs (downregulated: TMEM51-AS1, SND1-IT1; upregulated: HCP5, RUSC1-AS1, LINC00324) and no DEMs were overlapped with the DMRs, but only a negative relationship occurred in the expression and methylation level of TMEM51-AS1. Five DELs could interact with 11 DEMs to regulate 242 DEGs, which was used to construct the ceRNA network, including TMEM51-AS1-miR-106b-SNX21/ TRAPPC10, LINC00324/RUSC1-AS1-miR-16-SPRY4/MICAL2/ SLC39A14, RUSC1-AS1-miR-10-SCG5 and RUSC1-AS1-miR-7-ZFP1 ceRNAs axes. Univariate Cox regression analysis showed RUSC1-AS1 and SNX21 were associated with overall survival (OS); LINC00324, miR-7 and ZFP1 correlated with recurrence-free survival (RFS); miR-16, miR-10, SCG5, SPRY4, MICAL2 and SLC39A14 were both OS and RFS-related. Furthermore, TRAPPC10 and SLC39A14 were identified as independent OS prognostic factors by multivariate Cox regression analysis. CONCLUSION: DNA methylation-mediated TMEM51-AS1 and non-methylation-mediated RUSC1-AS1 may function as ceRNAs for induction of LSCC. They and their ceRNA axis genes (particularly TMEM51-AS1-miR-106b-TRAPPC10; RUSC1-AS1-miR-16-SLC39A14) may be potentially important prognostic biomarkers for LSCC.