Abstract
Long non-coding RNA (lncRNA) can act as ceRNA to regulate the expression of target genes by sponging miRNAs, and therefore plays an essential role in tumor initiation and progression. However, functional roles and regulatory mechanisms of lncRNAs as ceRNAs in head and neck squamous cell carcinoma (HNSCC) remain to be determined. We downloaded RNA sequence profiles from The Cancer Genome Atlas (TCGA) database, and identified the differential RNAs by bioinformatics. Then we analyzed the biological processes of differential expressed RNAs (DER), and established their interaction networks and pathway analysis to find out potential biological effects of these DERs. Besides, we also explored the relationship between the DERs and prognosis of HNSCC patients. We obtained 525 tumor samples and 44 paracancerous controls, and there were 1081 DElncRNAs, 1889 DEmRNAs, and 145 DEmiRNAs. GO and KEGG pathways analysis of these DEmRNAs were mainly involved in “Protein digestion and absorption,” “Calcium signaling pathway,” and “ECM-receptor interaction.” The analysis of the ceRNA network identified 61 DElncRNAs as functional ceRNAs whose dysregulated expression may affect the expression of oncogenes/tumor suppressor genes. Furthermore, univariate and multivariate Cox regression analysis revealed that 4 DElncRNAs, 3 EDmiRNAs, and 6 DEmRNAs can predict survival with high accuracy. Survival analysis found that 4 lncRNAs was related to prognostic, including overexpressed RP11-366H4.1, HOTTIP, RP11-865I6.2, and RP11-275N1.1 patients had a worse survival. In conclusion, through constructing the ceRNA network in HNSCC patients, we identified key lncRNA-miRNA-mRNA network in HNSCC. All the DERs might participate in varieties of pathways in the initiation, progression, and invasion of HNSCC. Furthermore, some miRNAs (hsa-mir-99a, hsa-mir-337, and hsa-mir-137) and mRNAs (NOSTRIN, TIMP4, GRB14, HOXB9, CELSR3, and ADGRD2) may be the prognostic genes of HNSCC. This study provided a new target and theoretical basis for further research on molecular mechanisms and biomarkers.
Introduction
Head and neck cancer refers to malignant tumors derived from the nasal cavity, paranasal sinuses, nasopharynx, oral cavity, pharynx, and larynx; the majority of these tumors are squamous cell carcinomas and their carcinogenesis has been associated with cigarette smoking (Hashibe et al., 2009), alcoholism (Hashibe et al., 2009) and HPV (Herrero et al., 2003). Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer worldwide with an incidence of over 650 000 new cases each year (Torre et al., 2015). Despite recent advances in molecular pathology and targeted therapies have remarkably improved the prognosis of patients with HNSCC, there are no specific biomarkers or comparable effective targeted molecular therapies available to screen and treat HNSCC patients, and the 5-year survival for patients with HNSCC still remains low (Argiris et al., 2008). There is an urgent requirement for the investigation of prognostic biomarkers and treatment options based on the current genomic approaches for HNSCC.
The development of high-throughput sequencing endowed a powerful tool to expand the understanding of gene expression mechanisms. It has been demonstrated that non-coding RNAs play important roles in regulating gene expression (Wang and Chang, 2011). Non-coding RNAs can be divided into microRNA (miRNA, 20–23 nucleotides) and long non-coding RNA (lncRNA, >200 nucleotides) according to the transcript size (Derrien et al., 2012). Plenty of studies have discovered that miRNAs can post-transcriptionally inhibit translation or repress mRNAs expression through binding to the 3′-untranslated region (3′-UTR) of the target genes (Chan and Wang, 2015), leading to cleavage of target mRNAs and/or inhibition of their translation. Currently, it is widely accepted that miRNAs may act as oncogenes or suppressor genes during tumor initiation and progression.
In addition to the well-annotated protein-coding genes and miRNA genes, lncRNAs have also emerged as important regulatory molecules of a variety of physiological and pathological process (Kopp and Mendell, 2018). Studies identified that lncRNA post-transcriptionally regulates the expression of mRNA by competing with miRNA (Feng et al., 2018). Interference of dysregulated lncRNA and mRNAs influence the expression of target mRNAs through miRNAs. This miRNA-regulated lncRNA and mRNA network is a part of the “competing endogenous RNA (ceRNA) hypothesis” (Thomson and Dinger, 2016). The ceRNA crosstalk described that ceRNAs, as miRNA sponges, harbored the same miRNA response elements (MREs) and communicated with each other by competing for shared miRNAs (Thomson and Dinger, 2016). It has proved that lncRNAs could disturb the balance of ceRNA network, thus resulting in the initiation and progress of cancers (Ding et al., 2018; Zhang et al., 2018).
In the current study, to identify the aberrant expression profile of ncRNAs in HNSCC patients and further study the roles of lncRNAs in HNSCC progress, genome-wide profile of ncRNAs was screened by using high-throughput sequencing in HNSCC patients. Further, we analyzed and predicted the functions of the aberrant lncRNAs following the lncRNA-miRNA-mRNA network.
Materials and Methods
Study Population and RNA-Sequencing Data Processing
RNA-sequencing data of 546 samples with HNSCC were retrieved from the TCGA data portal1. The RNA expression data (level 3) and clinical data of 546 HNSCC cases were downloaded from the TCGA data portal. The RNA- and miRNA-sequencing data from the 546 samples which were free to download were derived from the IlluminaHiSeq_RNASeq and IlluminaHiSeq_miRNASeq sequencing platforms. The sequencing data of the 546 samples contained the corresponding RNA-seq and miRNA-seq data and were divided into 502 tumor samples and 44 normal samples. In the current study, we mainly took the program code written in Perl and R language to analyze RNA data.
Identification of DEIncRNAs, DEmiRNAs, and DEmRNAs
We identified mRNAs and lncRNAs by using the Ensembl database2. Before conducting differential expression analysis, we ruled out all unexpressed RNAs by removing all rows with a mean read of less than or equal to one. We analyzed the DEIncRNAs, DEmiRNAs, and DEmRNAs by using the using “edgeR” package, a bioconductor package via R language procedure. All P-values used false discovery rate (FDR) to correct the statistical significance of the multiple test. |FoldChange| ≧2 and FDR < 0.01 were considered significant. For the obtained DEIncRNAs, DEmiRNAs, and DEmRNAs, we generated heat maps using the heatmap packages in the R software.
Construction of ceRNA Network
RNAs sharing the same miRNA regulators affects each other by competing for the limited number of miRNA sponges based on the ceRNA hypothesis. The construction of the ceRNA network included: (a) We obtained the DElncRNAs, DEmiRNA, and DEmRNAs based on above method. (b) The experimentally validated miRNA-lncRNA interactions were downloaded from miRcode database3 and then the DElncRNAs and DEmiRNAs were matched with them to obtain their connections, an lncRNA-miRNA regulatory network. Expression correlation (pearson) between DElncRNAs and miRNAs was assessed. (c) Human miRNAs and their targets data were collected as miRNA-gene pairs concurrenced based on the miRDB4, miRTarBase5, and TargetScan6. (d) These two networks, lncRNA-miRNA, and a miRNA-mRNA network were integrated into a comprehensive ceRNA regulatory network to demonstrate the interactions of DEmiRNAs with DElncRNAs or DEmRNAs. The ceRNA network was visualized using Cytoscape v3.4.07.
Functional Enrichment Analysis
Functional enrichment analysis among GO and KEGG pathway levels was performed based on DAVID Bioinformatics Resources (8, version 6.8) and KOBAS 3.0 software9.
Construction of the Prognostic Risk Score
DE genes were assessed by univariate Cox proportional hazards regression analysis to obtain prognostic values (P < 0.001), and the most prominent 10 genes were then used to perform a multivariate Cox proportional hazards regression analysis. A risk score formula was established by including each of the prognosis related genes, weighted by their estimated regression coefficients in the multivariable Cox regression analysis. The risk score for each patient was assessed, and patients were classified into high-risk score or low-risk score groups by using the corresponding median risk score as the cut-off. A receiver operating characteristic (ROC) curve was obtained by using R with survival ROC package. Survival difference between the low-risk and high risk group were calculated by the Kaplan Meier and log-rank test.
Survival Analysis
We further study prognostic DERNAs signature by combining the clinical data of those HNSCC patients. Kaplan Meier curve analysis was performed to evaluate the univariate survival analysis. P < 0.05 was regarded as significant unless specifically indicated.
Results
Patient Characteristics
The detailed clinical and pathological characteristics of patients were summarized in Table 1. The exclusion criteria were samples without complete data for further analysis. Finally, for RNA-Sequence data (level 3) we matched a total of 501 tumor tissues and 43 paracancerous controls while for miRNA-Sequence data from the same cohort we matched 525 tumor tissues and 44 paracancerous controls.
Table 1
| Variable | Number (%) | |
|---|---|---|
| Age | ||
| ≤60 | 244 (51.2%) | |
| >60 | 256 (48.8%) | |
| Unknown | 1 (0.2%) | |
| Sex | ||
| Male | 368 (73.5%) | |
| Female | 133 (26.5%) | |
| Stage | ||
| I | 25 (5.0%) | |
| II | 71 (14.2%) | |
| III | 78 (15.6%) | |
| IV | 259 (51.7%) | |
| Unknown | 68 (13.6%) | |
| Vital Status | ||
| Alive | 283 (56.5%) | |
| Dead | 218 (53.5%) |
Characteristics of the patients (n = 501).
Differentially Expressed (DE) lncRNAs, mRNAs, and miRNAs in HNSCC
RNAs expression profiles of HNSCC patients and corresponding clinical data were downloaded from TCGA database. With the cut-off criteria unified |FoldChange| >2 and FDR < 0.01, obtained from TCGA and finally sorted out 1081 DElncRNAs (759 up-regulated and 322 down-regulated) and 1889 DEmRNAs (726 up-regulated and 1163 down-regulated). To further establish an lncRNA-miRNA-mRNA ceRNA network, we also matched miRNA expression profiles in 569 cases of HNSCC patients among the same cohort. As a result, a total of 145 DEmiRNAs (90 up-regulated and 55 down-regulated) were sorted using the cut-off unified |FoldChange| >1.5 and FDR < 0.05. We outlined DElncRNAs, DEmRNAs, and DEmiRNAs using heatmap in Figure 1.
FIGURE 1
Reconstruction of the ceRNA Network in HNSCC
To better understand the role of DElncRNAs in HNSCC and to further elucidate the interaction among the DElncRNAs, DEmRNA, and DEmiRNAs, we constructed lncRNA-miRNA-mRNA related ceRNA network. Based on miRcode database, which masters the function of transcriptome-wide mircoRNA targeting prediction including lncRNAs, we outlined 71 DElncRNAs and then we matched above DElncRNAs with previous 145 DEmiRNA, finally, we matched 60 DElncRNAs and 13 DEmiRNAs. Using above 13 DEmiRNAs we predict 500 miRNA-target genes based on these different miRNA-mRNA databases, of which 19 mRNAs were included in the 1888 DEmRNAs. Finally, we sorted out the association between DElncRNAs and DEmiRNAs and the link between DEmiRNAs and DEmRNAs. In the ceRNA network, there were 60 DElncRNAs, 13 DEmiRNAs, and 19 DEmRNAs (Figure 2).
FIGURE 2
GO and KEGG Pathway Analysis of DEmRNAs
To further analyze the functional characteristics of DEmRNAsin HNSCC. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed using DAVID Bioinformatics Tool (version 6.8)10 and KOBAS 3.011 websites, respectively. The results were divided into almost three parts, including cellular component, biological process, and molecular function. With the criteria of P-Value < 0.01, 38 GO terms were included in Figure 3A. We also analyzed KEGG pathway of the DEmRNAs using KOBAS 3.0. Almost 12 KEGG pathways were significantly enriched in our analysis. The majority of the enriched pathways were associated with cancers including “Protein digestion and absorption, Calcium signaling pathway, ECM-receptor interaction” (Figure 3B,C).
FIGURE 3
Survival Evaluation of the DElncRNAs, mRNAs, and miRNAs
We obtained 497 cases which had both lncRNA and mRNA expression data and 521 cases which had miRNA expression data for the prognostic evaluation. 496 overlapped cases had lncRNA, mRNA, and miRNA expression data. Using the univariate Cox proportional hazards regression method, 210 lncRNAs, 350 mRNAs, and 41 miRNAs that closely related to the survival of HNSCC with p-value < 0.05 were obtained. We select the most closely related 10 lncRNAs, mRNAs, and miRNAs to perform multivariate Cox proportional hazard regression analysis which reduced the amount to be 4 lncRNAs, 6 mRNAs, and 3 miRNAs that can all predict survival as integral with highest accuracy. The expression heatmap of the lncRNA, mRNA, and miRNA signature is shown in Figure 4–6A–C, which is grouped into low risk and high risk groups based on the median value of risk score of each sample. This analysis can predict 1, 3, 5, and 10 years survival rate accurately for HNSCC patients. The time-dependent receiver operating characteristic (ROC) curves for lncRNAs (Figure 4D–G), mRNAs (Figure 5D–G), and miRNAs (Figure 6D–G) signature have area under curve (AUC) values higher than 0.6.
FIGURE 4
FIGURE 5
FIGURE 6
Next, we performed single gene survival analysis from the above 4-lncRNA, 6-mRNA, and 3-miRNA to see if these genes can act independently as prognostic biomarker, based on low and high expression group which is divided by the median value of gene expression level. Among these genes, 4 lncRNAs, 6 mRNAs, and 3 miRNAs showed significant differential survival between high expression and low expression samples with p-value < 0.05 which is shown in Figure 7–9. The expression profile of the 4 lncRNAs, 6 mRNAs, and 3 miRNAs in tumor sample all showed significant differential expression as compared to normal tissue samples and are shown in Figure 10–12, respectively.
FIGURE 7
FIGURE 8
FIGURE 9
FIGURE 10
FIGURE 11
FIGURE 12
Discussion
The prediction of HNSCC prognosis largely depends on the TNM staging system and histologic grade (Alsaffar et al., 2016). It is very important to investigate novel prognostic biomarkers for HNSCC.
Benefit from the development of sequencing technology, the genome-wide profiling has extensively discovered that almost 98% of the transcriptional outputs are ncRNAs. The recent implication of lncRNAs in many biological functions has established a new scenario to better understand complex processes like cancer (Joung et al., 2017; Bester et al., 2018). Meanwhile, miRNAs have been well investigated in the regulation of transcription, epigenetics modulation, and RNA-protein interaction (Gorbea et al., 2017; Thomou et al., 2017). Indeed, some miRNAs have been proven to be potential diagnostic and therapeutic markers for HNSCC (Koshizuka et al., 2017; Chou et al., 2018). However, the combination of lncRNAs and miRNAs is currently not reported for prognosis. In this study, we developed a cox multivariate regression model with survival-related genes. Afterward, the performance of the lncRNA-miRNA-mRNA based risk score was assessed. Patients in the high risk group significantly had a worse survival than that in the low risk group.
GO analysis and KEGG pathway analysis also suggested that the differential expressed genes were associated with several signal pathways, including extracellular matrix organization, collagen fibril organization, which were associated with the pathogenesis of HNSCC (Kuang et al., 2016; Yang et al., 2017). lncRNAs were recognized to regulate downstream signaling pathways by functioning as ceRNA in tumorigenesis (Chen et al., 2018; Li et al., 2018). In our work, we successfully predicted lncRNA/miRNA/mRNA interactions, which may shed a light to reveal the information for ceRNA research in HNSCC.
Accumulative evidence shows that ncRNAs play a very important regulatory role in tumor progression of HNSCC. For instance, lncRNA HOX transcript antisense RNA (HOTAIR) was found to modulate the progression of HNSCC (Sun et al., 2018). HNSCC samples showed significantly robust expression/activation of HOTAIR compared with normal squamous epithelium. Overexpression of HOTAIR promoted the growth of xenograft tumors in vivo (Sun et al., 2018). In our ceRNA network, the high expression of the HOTAIR and HOTTIP are closely associated with 6 miRNAs (mir-301b, mir-193b, and mir-31 are upregulated, whilst mir-216a, mir-206, and mir-375 are downregulated). The highly homologous miRNA to hsa-mir-301b, hsa-miR-301a-3p, played a role in the emergence and development of laryngeal squamous cell carcinoma by directly regulating the Smad4 (Lu et al., 2015). MiR-193b was over-expressed in HNSCC cell lines. Knockdown of miR-193b in FaDu cancer cells substantially reduced cell proliferation, migration and invasion, along with suppressed tumor formation (Lenarduzzi et al., 2013). HNSCC patients whose tumors expressed high levels of miR-193b experienced a lower disease-free survival compared to patients with low miR-193b expression (Lenarduzzi et al., 2013). The miR-31 was significantly increased in patients with oral carcinoma at all clinical stages, including very small tumors (Liu et al., 2012), miR-31 also facilitate migration and invasion by targeting Numb in HNSCC cells (Chou et al., 2018). miR-206 and miR-375 were down-regulated in HNSCC clinical specimens and cell lines (Nohata et al., 2011; Liu et al., 2017), which is in agreement with our results. Finally, expression level of miR-216-5p in cervical cancer tissues was observably lower than that in corresponding normal tissues (Zhu et al., 2018).
Considering all the identified lncRNAs, HOTTIP’s abnormal expression had the significant impact on the survival of patients (HR = 1.16; P = 0.00132), our data was in agreement with an recent meta-analysis that high HOTTIP expression was significantly correlated with poor OS in cancer patients (Fan et al., 2018). Meanwhile, miRNAs (hsa-mir-99a, hsa-mir-337, and hsa-mir-137) and mRNA (NOSTRIN, TIMP4, GRB14, HOXB9, CELSR3, ADGRD2) were showed to be significant independent factors for poor survival in HNSCC. hsa-mir-99a (Hou et al., 2015), hsa-mir-337 (Wang et al., 2013), and hsa-mir-137 (Langevin et al., 2011) have been reported to be associated with the diagnosis, proliferation, metastasis, and prognosis of various cancers, including HNSCC. Ma et al. (2009) revealed that TIMP4 were differentially expressed and could distinguish cancerous and non-cancerous samples. GRB14 (Morzyglod et al., 2017), HOXB9 (Sun et al., 2017), and CELSR3 (Karpathakis et al., 2016), have been proven to be associated with proliferation and prognosis in multiple cancers, including HNSCC. However, NOSTRIN and ADGRD2 were identified for the first time to be associated with prognosis of HNSCC. More experiments not only in cell lines, but also in animal model, are still needed to investigate the underlying regulatory mechanisms of these ncRNAs.
Taken together, our study identified the lncRNAs expression profile in HNSCC. With the construction of a ceRNA crosstalk network, we constructed a perspective to screen lncRNAs that could be involved in HNSCC tumorigenesis. Moreover, the lncRNA-associated ceRNA network will enable us to better understand the pathogenesis of HNSCC and provide novel lncRNAs and miRNAs as candidate prognosis biomarkers or potential therapeutic targets.
Statements
Author contributions
YP conceived the study and drafted the manuscript. GL acquired, analyzed, and interpreted the data and reviewed the manuscript. YL reviewed the manuscript and supervised the study. DW helped with discussion and comments on the project.
Funding
This work was supported by the National Natural Science Foundation of China (81872200), the Natural Science Foundation of Hubei Province (2018CFB510), the Zhongnan Hospital of Wuhan University Science, Technology and Innovation Seed Fund (CXPY2017029), and the Fundamental Research Funds for the Central Universities (2042018kf0091).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
1.^http://cancergenome.nih.gov/
5.^http://mirtarbase.mbc.nctu.edu.tw/
6.^http://www.targetscan.org/vert_71/
8.^http://david.abcc.ncifcrf.gov/
9.^http://kobas.cbi.pku.edu.cn/index.php
References
1
AlsaffarH. A.GoldsteinD. P.KingE. V.de AlmeidaJ. R.BrownD. H.GilbertR. W.et al (2016). Correlation between clinical and MRI assessment of depth of invasion in oral tongue squamous cell carcinoma.J. Otolaryngol. Head Neck Surg.45:61. 10.1186/s40463-016-0172-0
2
ArgirisA.KaramouzisM. V.RabenD.FerrisR. L. (2008). Head and neck cancer.Lancet3711695–1709. 10.1016/S0140-6736(08)60728-X
3
BesterA. C.LeeJ. D.ChavezA.LeeY. R.NachmaniD.VoraS.et al (2018). An integrated genome-wide CRISPRa approach to functionalize lncRNAs in drug resistance.Cell173649–664e620. 10.1016/j.cell.2018.03.052
4
ChanS. H.WangL. H. (2015). Regulation of cancer metastasis by microRNAs.J. Biomed. Sci.22:9. 10.1186/s12929-015-0113-7
5
ChenX.ChenZ.YuS.NieF.YanS.MaP.et al (2018). Long noncoding RNA LINC01234 functions as a competing endogenous RNA to regulate CBFB expression by sponging miR-204-5p in gastric cancer.Clin. Cancer Res.242002–2014. 10.1158/1078-0432.CCR-17-2376
6
ChouC. H.TuH. F.KaoS. Y.ChiangC. F.LiuC. J.ChangK. W.et al (2018). Targeting of miR-31/96/182 to the Numb gene during head and neck oncogenesis.Head Neck40808–817. 10.1002/hed.25063
7
DerrienT.JohnsonR.BussottiG.TanzerA.DjebaliS.TilgnerH.et al (2012). The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression.Genome Res.221775–1789. 10.1101/gr.132159.111
8
DingJ.YehC. R.SunY.LinC.ChouJ.OuZ.et al (2018). Estrogen receptor beta promotes renal cell carcinoma progression via regulating LncRNA HOTAIR-miR-138/200c/204/217 associated CeRNA network.Oncogene375037–5053. 10.1038/s41388-018-0175-6
9
FanY.YanT.ChaiY.JiangY.ZhuX. (2018). Long noncoding RNA HOTTIP as an independent prognostic marker in cancer.Clin. Chim. Acta482224–230. 10.1016/j.cca.2017.07.031
10
FengL.ShiL.LuY. F.WangB.TangT.FuW. M.et al (2018). Linc-ROR promotes osteogenic differentiation of mesenchymal stem cells by functioning as a competing endogenous RNA for miR-138 and miR-145.Mol. Ther. Nucleic Acids11345–353. 10.1016/j.omtn.2018.03.004
11
GorbeaC.MosbrugerT.CazallaD. (2017). A viral Sm-class RNA base-pairs with mRNAs and recruits microRNAs to inhibit apoptosis.Nature550275–279. 10.1038/nature24034
12
HashibeM.BrennanP.ChuangS. C.BocciaS.CastellsagueX.ChenC.et al (2009). Interaction between tobacco and alcohol use and the risk of head and neck cancer: pooled analysis in the international head and neck cancer epidemiology consortium.Cancer Epidemiol. Biomarkers Prev.18541–550. 10.1158/1055-9965.EPI-08-0347
13
HerreroR.CastellsagueX.PawlitaM.LissowskaJ.KeeF.BalaramP.et al (2003). Human papillomavirus and oral cancer: the international agency for research on cancer multicenter study.J. Nation. Cancer Inst.951772–1783. 10.1093/jnci/djg107
14
HouB.IshinagaH.MidorikawaK.ShahS. A.NakamuraS.HirakuY.et al (2015). Circulating microRNAs as novel prognosis biomarkers for head and neck squamous cell carcinoma.Cancer Biol. Ther.161042–1046. 10.1080/15384047.2015.1045692
15
JoungJ.EngreitzJ. M.KonermannS.AbudayyehO. O.VerdineV. K.AguetF.et al (2017). Genome-scale activation screen identifies a lncRNA locus regulating a gene neighbourhood.Nature548343–346. 10.1038/nature23451
16
KarpathakisA.DibraH.PipinikasC.FeberA.MorrisT.FrancisJ.et al (2016). Prognostic impact of novel molecular subtypes of small intestinal neuroendocrine tumor.Clin. Cancer Res.22250–258. 10.1158/1078-0432.CCR-15-0373
17
KoppF.MendellJ. T. (2018). Functional classification and experimental dissection of long noncoding RNAs.Cell172393–407. 10.1016/j.cell.2018.01.011
18
KoshizukaK.HanazawaT.AraiT.OkatoA.KikkawaN.SekiN. (2017). Involvement of aberrantly expressed microRNAs in the pathogenesis of head and neck squamous cell carcinoma.Cancer Metastasis Rev.36525–545. 10.1007/s10555-017-9692-y
19
KuangJ.ZhaoM.LiH.DangW.LiW. (2016). Identification of potential therapeutic target genes and mechanisms in head and neck squamous cell carcinoma by bioinformatics analysis.Oncol. Lett.113009–3014. 10.3892/ol.2016.4358
20
LangevinS. M.StoneR. A.BunkerC. H.Lyons-WeilerM. A.LaFramboiseW. A.KellyL.et al (2011). MicroRNA-137 promoter methylation is associated with poorer overall survival in patients with squamous cell carcinoma of the head and neck.Cancer1171454–1462. 10.1002/cncr.25689
21
LenarduzziM.HuiA. B.AlajezN. M.ShiW.WilliamsJ.YueS.et al (2013). MicroRNA-193b enhances tumor progression via down regulation of neurofibromin 1.PLoS One8:e53765. 10.1371/journal.pone.0053765
22
LiZ.JiangP.LiJ.PengM.ZhaoX.ZhangX.et al (2018). Tumor-derived exosomal lnc-Sox2ot promotes EMT and stemness by acting as a ceRNA in pancreatic ductal adenocarcinoma.Oncogene373822–3838. 10.1038/s41388-018-0237-9
23
LiuC. J.LinS. C.YangC. C.ChengH. W.ChangK. W. (2012). Exploiting salivary miR-31 as a clinical biomarker of oral squamous cell carcinoma.Head Neck34219–224. 10.1002/hed.21713
24
LiuF.ZhaoX.QianY.ZhangJ.ZhangY.YinR. (2017). MiR-206 inhibits Head and neck squamous cell carcinoma cell progression by targeting HDAC6 via PTEN/AKT/mTOR pathway.Biomed. Pharmacother.96229–237. 10.1016/j.biopha.2017.08.145
25
LuY.GaoW.ZhangC.WenS.HuangfuH.KangJ.et al (2015). Hsa-miR-301a-3p acts as an oncogene in laryngeal squamous cell carcinoma via target regulation of Smad4.J. Cancer61260–1275. 10.7150/jca.12659
26
MaL. J.LiW.ZhangX.HuangD. H.ZhangH.XiaoJ. Y.et al (2009). Differential gene expression profiling of laryngeal squamous cell carcinoma by laser capture microdissection and complementary DNA microarrays.Arch. Med. Res.40114–123. 10.1016/j.arcmed.2008.12.005
27
MorzyglodL.CauzacM.PopineauL.DenechaudP. D.FajasL.RagazzonB.et al (2017). Growth factor receptor binding protein 14 inhibition triggers insulin-induced mouse hepatocyte proliferation and is associated with hepatocellular carcinoma.Hepatology651352–1368. 10.1002/hep.28972
28
NohataN.HanazawaT.KikkawaN.MutallipM.SakuraiD.FujimuraL.et al (2011). Tumor suppressive microRNA-375 regulates oncogene AEG-1/MTDH in head and neck squamous cell carcinoma (HNSCC).J. Hum. Genet.56595–601. 10.1038/jhg.2011.66
29
SunC.HanC.WangP.JinY.SunY.QuL. (2017). HOXB9 expression correlates with histological grade and prognosis in LSCC.Biomed. Res. Int.2017:3680305. 10.1155/2017/3680305
30
SunS.WuY.GuoW.YuF.KongL.RenY.et al (2018). STAT3/HOTAIR signaling axis regulates HNSCC growth in an EZH2-dependent manner.Clin. Cancer Res.242665–2677. 10.1158/1078-0432.CCR-16-2248
31
ThomouT.MoriM. A.DreyfussJ. M.KonishiM.SakaguchiM.WolfrumC.et al (2017). Adipose-derived circulating miRNAs regulate gene expression in other tissues.Nature542450–455. 10.1038/nature21365
32
ThomsonD. W.DingerM. E. (2016). Endogenous microRNA sponges: evidence and controversy.Nat. Rev. Genet.17272–283. 10.1038/nrg.2016.20
33
TorreL. A.BrayF.SiegelR. L.FerlayJ.Lortet-TieulentJ.JemalA. (2015). Global cancer statistics, 2012.CA Cancer J. Clin.6587–108. 10.3322/caac.21262
34
WangK. C.ChangH. Y. (2011). Molecular mechanisms of long noncoding RNAs.Mol. Cell43904–914. 10.1016/j.molcel.2011.08.018
35
WangZ.WangJ.YangY.HaoB.WangR.LiY.et al (2013). Loss of has-miR-337-3p expression is associated with lymph node metastasis of human gastric cancer.J. Exp. Clin. Cancer Res.3276. 10.1186/1756-9966-32-76
36
YangB.ChenZ.HuangY.HanG.LiW. (2017). Identification of potential biomarkers and analysis of prognostic values in head and neck squamous cell carcinoma by bioinformatics analysis.Onco Targets Ther.102315–2321. 10.2147/OTT.S135514
37
ZhangG.LiS.LuJ.GeY.WangQ.MaG.et al (2018). LncRNA MT1JP functions as a ceRNA in regulating FBXW7 through competitively binding to miR-92a-3p in gastric cancer.Mol. Cancer17:87. 10.1186/s12943-018-0829-6
38
ZhuH.ZengY.ZhouC. C.YeW. (2018). SNHG16/miR-216-5p/ZEB1 signal pathway contributes to the tumorigenesis of cervical cancer cells.Arch. Biochem. Biophys.6371–8. 10.1016/j.abb.2017.11.003
Summary
Keywords
head and neck squamous cell carcinoma, ceRNA network, long non-coding RNAs, microRNA, TCGA
Citation
Pan Y, Liu G, Wang D and Li Y (2019) Analysis of lncRNA-Mediated ceRNA Crosstalk and Identification of Prognostic Signature in Head and Neck Squamous Cell Carcinoma. Front. Pharmacol. 10:150. doi: 10.3389/fphar.2019.00150
Received
10 July 2018
Accepted
08 February 2019
Published
04 March 2019
Volume
10 - 2019
Edited by
Francois X. Claret, The University of Texas MD Anderson Cancer Center, United States
Reviewed by
Valerio Costa, Italian National Research Council (CNR), Italy; Dingbo Shi, Sun Yat-sen University, China
Updates
Copyright
© 2019 Pan, Liu, Wang and Li.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guohong Liu, liuguoh@outlook.com Yirong Li, liyirong838@163.com
This article was submitted to Cancer Molecular Targets and Therapeutics, a section of the journal Frontiers in Pharmacology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.