Excavating novel diagnostic and prognostic long non-coding RNAs (lncRNAs) for head and neck squamous cell carcinoma: an integrated bioinformatics analysis of competing endogenous RNAs (ceRNAs) and gene co-expression networks

Bioengineered. 2021 Dec;12(2):12821-12838. doi: 10.1080/21655979.2021.2003925.

Abstract

Long non-coding RNAs (lncRNAs) have been demonstrated to fine-tune gene regulations that govern a broad spectrum of oncogenic processes. Nonetheless, our understanding of the roles of lncRNAs and their interactions with miRNAs and mRNAs in HNSCC is still highly rudimentary. Here, we present a comprehensive bioinformatics analysis in which competing endogenous RNA (ceRNA) network construction and weighted gene co-expression network analysis (WGCNA) were combined to explore novel diagnostic and prognostic lncRNAs for HNSCC. Differentially expressed mRNAs (DEGs), miRNAs (DEMs) and lncRNAs (DELs) were identified based on the RNA sequencing data and clinical data retrieved from TCGA database. LncRNA-regulated ceRNA networks were constructed based on the interactive RNA pairs predicted by miRDB, miRcode and TargetScan. WGCNA was conducted to identify lncRNAs that were significantly correlated with patient overall survival (OS) and HNSCC tumor. RT-qPCR was employed to validate the expression of lncRNAs in HNSCC cell lines and patient sera. A ceRNA network consisting of 90 DEGs, 7 DEMs and 67 DELs associated with clinical traits was established. WGCNA and Kaplan-Meier survival analysis revealed that 5 DELs (MIR4435-2 HG, CASC9, LINC01980, STARD4-AS1 and MIR99AHG) were significantly correlated with OS of HNSCC patients, whereas DEL PART1 was most significantly correlated with the HNSCC tumor. The in silico predicted expression patterns of PART1, LINC01980 and MIR4435-2 HG were further validated in HNSCC cell lines and patient sera. Collectively, the present study provided novel insights into the lncRNA-regulated ceRNA networks in HNSCC and identified novel lncRNAs that harbor diagnostic and prognostic potentials for HNSCC.Abbreviations BP, biological process. CC, cellular component. ceRNA, competing endogenous RNA. DEG, differential expressions of mRNA. DEL, differentially expressed lncRNA. DEM, differentially expressed miRNA. ESCC, esophageal squamous cell carcinoma. FPKM, Fragments Per Kilobase Million. GO, Gene Ontology. GS, gene significance. HNSCC, head and neck squamous cell carcinoma. KEGG, Kyoto Encyclopedia of Genes and Genomes. LncRNA, long non-coding RNA. MCC, Maximal Clique Centrality. ME, module eigengenes. MF, molecular functions. MM, module membership. MRE, miRNA-binding site. MYO5A, Myosin-Va. PART1, prostate androgen-regulated transcript 1. RBM3, RNA‑binding motif protein 3. TCGA, The Cancer Genome Atlas. TOM, topological overlap measure. TSCC, tongue squamous cell carcinoma. WGCNA, weighted gene co-expression network analysis.

Keywords: Head and neck squamous cell carcinoma; WGCNA; biomarkers; ceRNA network; lncRNA.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Cell Line, Tumor
  • Computational Biology*
  • Gene Expression Profiling
  • Gene Expression Regulation, Neoplastic*
  • Gene Regulatory Networks*
  • Head and Neck Neoplasms / diagnosis*
  • Head and Neck Neoplasms / genetics*
  • Humans
  • Kaplan-Meier Estimate
  • MicroRNAs / genetics
  • MicroRNAs / metabolism
  • Prognosis
  • RNA, Long Noncoding / genetics
  • RNA, Long Noncoding / metabolism*
  • RNA, Messenger / genetics
  • RNA, Messenger / metabolism
  • Reproducibility of Results
  • Squamous Cell Carcinoma of Head and Neck / diagnosis*
  • Squamous Cell Carcinoma of Head and Neck / genetics*
  • Survival Analysis

Substances

  • MicroRNAs
  • RNA, Long Noncoding
  • RNA, Messenger

Grants and funding

This work was supported by the Guangxi Key Laboratory of Early Prevention and Treatment for Regional High Frequency Tumor [GK2018-08, GKE2019-15, GKE-ZZ202012, GKE-ZZ202128]; the General Program of Guangxi Natural Science Foundation [2018GXNSFAA281032].