Dataset Information

Identifying prognostic lncRNAs based on a ceRNA regulatory network in laryngeal squamous cell carcinoma.

ABSTRACT:

Purpose

Growing evidence demonstrates that long non-coding RNAs (lncRNAs) play a crucial role as competing endogenous RNAs (ceRNAs) in tumor occurrence. The lncRNAs' functions and clinical significance in laryngeal squamous cell carcinoma (LSCC) remain unclear. The study aims to reveal the lncRNA-associated ceRNA regulatory network of LSCC and clarify its clinical relevance.

Methods

Here, we obtained LSCC transcriptome data from The Cancer Genome Atlas (TCGA) database and identified the differential expression profile of lncRNAs, miRNAs, and mRNAs by the EdgeR R package. The function enrichment analysis of mRNAs was performed using clusterProfiler R package and GSEA3.0. Then, we constructed a ceRNA network and prognosis model based on lncRNAs through bioinformatic methods. Moreover, we explored the functions of prognosis-related lncRNA in LSCC by CCK-8 and transwell assay.

Results

1961 lncRNAs, 69 miRNAs, and 2224 mRNAs were identified as differentially expressed genes in LSCC tissues. According to the transcriptome differential expression profile, a ceRNA network containing 61 lncRNAs, 21 miRNAs, and 77 mRNAs was established. Then, four lncRNAs (AC011933.2, FAM30A, LINC02086, LINC02575) were identified from the ceRNA network to build a prognosis model for LSCC patients. And we found that LINC02086 and LINC02575 promoted the proliferation, migration, and invasion of LSCC cells while AC011933.2 and FAM30A inhibited these biological functions in vitro. Furthermore, we validated that LINc02086/miR-770-5p/SLC26A2 axis promoted migration in LSCC.

Conclusion

Four lncRNAs of the ceRNA network were abnormally expressed and related to patient prognosis in LSCC. They played a significant role in the progress of LSCC via affecting the proliferation and metastasis of tumor cells.

SUBMITTER: Shi Y

PROVIDER: S-EPMC8207764 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:ObjectiveThis study aims to construct a systematic mRNA-miRNA-lncRNA network to identify novel lncRNAs and miRNAs biomarkers for laryngeal squamous cell carcinoma (LSCC).MethodsThe mRNA, miRNA and lncRNA expression profiles of LSCC were obtained from Gene Expression Omnibus (GEO) database. The differentially expressed mRNAs, miRNAs and lncRNAs (DEmRNAs, DEmiRNAs and DElncRNAs) were screened between LSCC tissues and controls. Functional analysis of DEmRNAs, DEmRNAs targeted by DEmiRNAs and DEmRNAs targeted by DElncRNAs were respectively performed. The miRWalk, starbase and DIANA-LncBase were respectively used to predict DEmiRNAs-DEmRNAs, DElncRNAs-DEmRNAs and DElncRNAs-DEmiRNAs pairs. ceRNA network was built by DEmiRNAs-DEmRNAs and DElncRNAs-DEmiRNAs pairs. LncRNA subcellular localization was predicted using lncLocator. Using published The Cancer Genome Atlas (TCGA) and external datasets (GSE127165 and GSE133632), we also validated the expression of key DElncRNAs and DEmiRNAs in ceRNA network. The diagnostic and prognostic value of candidate genes was evaluated by ROC curve analysis and survival analysis, respectively.ResultsThere were 5 mRNA datasets, 3 miRNA datasets and 2 lncRNA datasets in this study. Totally, 2957 DEmRNAs, 61 DElncRNAs and 23 DEmiRNAs were identified. Functional analysis of DEmRNAs shows that they were significantly enriched in cancer-related pathways, such as DNA replication and extracellular matrix organization. There were 11 DEmiRNAs, 17 DElncRNAs and 967 DEmRNAs in the ceRNA network. Notably, up-regulated lncRNA DGCR5-down-regulated has-miR-338-3p/has-miR-139-5p pairs in this network were experimentally validated. Moreover, down-regulated AL121839.2, down-regulated LINC02147, up-regulated AC079328.2, up-regulated AC004943.2 and up-regulated HMGA2-AS1 were located in the cytoplasm. AL121839.2 and LINC02147 interacted with has-miR-1246. AC004943.2, AC079328.2 and HMGA2-AS1 targeted has-miR-3185, has-miR-3137 and has-miR-582-5p, respectively. Based on the TCGA and external datasets (GSE127165 and GSE133632), DGCR5 and AC004943.2 were significantly up-regulated while AL121839.2 and LINC02147, has-miR-338-3p, has-miR-139-5p and has-miR-582-5p were significantly down-regulated, which were consistent with our integration analysis. DGCR5, AL121839.2, LINC02147, AC004943.2, has-miR-338-3p, has-miR-139-5p and has-miR-582-5p could predict the occurrence of LSCC. Survival analysis suggested that only, AL121839.2 has potential prognostic value for LSCC.ConclusionThis study provided novel insights into the ceRNA network and uncovered novel lncRNAs and miRNAs with diagnostic value in LSCC.

Project description:A cutaneous squamous cell carcinoma (cSCC) derived from keratinocytes is the second most common cause of non-melanoma skin cancer. The accumulation of the mutational burden of genes and cellular DNA damage caused by the risk factors (e.g., exposure to ultraviolet radiation) contribute to the aberrant proliferation of keratinocytes and the formation of a cSCC. A cSCC encompasses a spectrum of diseases that range from recursor actinic keratosis (AK) and squamous cell carcinoma (SCC) in situ (SCCIS) to invasive cSCCs and further metastatic SCCs. Emerging evidence has revealed that lncRNAs are involved in the biological process of a cSCC. According to the ceRNA regulatory theory, lncRNAs act as natural miRNA sponges and interact with miRNA response elements, thereby regulating the mRNA expression of their down-stream targets. This study was designed to search for the potential lncRNAs that may become potential therapeutic targets or biomarkers of a cSCC. Considering the spirit of the study to be adequately justified, we collected microarray-based datasets of 19 cSCC tissues and 12 normal skin samples from the GEO database (GSE42677 and GSE45164). After screening the differentially expressed genes via a limma package, we identified 24 differentially expressed lncRNAs (DElncRNAs) and 3221 differentially expressed mRNAs (DEmRNAs). The miRcode, miRTarBase, miRDB and TargetScan databases were used to predict miRNAs that could interact with DElncRNAs and DEmRNAs. A total of 137 miRNA-lncRNA and 221 miRNA-mRNA pairs were retained in the ceRNA network, consisting of 31 miRNAs, 11 DElncRNAs and 155 DEmRNAs. For the functional analysis, the top enriched biological process was enhancer sequence-specific DNA binding in Gene Ontology (GO) terms. The FoxO signaling pathway, autophagy and cellular senescence were the top enrichment terms based on a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The combination of a STRING tool and Cytoscape software (plug-in MCODE) identified five core mRNAs and built a core mRNA-associated ceRNA network. The expression for five identified core mRNAs and their related nine lncRNAs was validated using the external dataset GSE7553. Finally, one lncRNA HLA-F-AS1 and three mRNAs named AGO4, E2F1 and CCND1 were validated with the same expression patterns. We speculate that lncRNA HLA-F-AS1 may sponge miR-17-5p or miR-20b-5p to regulate the expression of CCND1 and E2F1 in the cSCC. The present study may provide potential diagnostic and therapeutic targets for cSCC patients.