Next Article in Journal
A Comparative Analysis between Subpectoral versus Prepectoral Single Stage Direct-to-Implant Breast Reconstruction
Previous Article in Journal
Atypical Enostoses—Series of Ten Cases and Literature Review
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Medicina
Volume 56
Issue 10
10.3390/medicina56100535
medicina-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Prognosis Associated microRNAs in HNSCC Subtypes Based on TCGA Dataset

by
Cintia M. Chamorro Petronacci
1,
Abel García García
1,
Elena Padín Iruegas
2,
Berta Rivas Mundiña
3,
Alejandro I. Lorenzo Pouso
1 and
Mario Pérez Sayáns
1,*
1
Faculty of Medicine and Dentistry, Health Research, Institute of Santiago de Compostela (Instituto de Investigación Sanitaria de Santiago, IDIS), Oral Medicine, Oral Surgery and Implantology University, Santiago de Compostela University, 15782 Santiago de Compostela, Spain
2
Department of Functional Biology and Health Sciences, Faculty of Physiotherapy, Human Anatomy and Embryology Area, Vigo University, 36001 Pontevedra, Spain
3
Pathology and Therapeutic Unity, Faculty of Medicine and Dentistry, University of Santiago de Compostela, 15782 Santiago de Compostela, Spain
*
Author to whom correspondence should be addressed.
Medicina 2020, 56(10), 535; https://doi.org/10.3390/medicina56100535
Submission received: 22 September 2020 / Revised: 7 October 2020 / Accepted: 9 October 2020 / Published: 13 October 2020
(This article belongs to the Section Genetics and Molecular Medicine)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background and Objectives: Head and Neck Squamous Cell Carcinoma (HNSCC) includes cancers from the oral cavity, larynx, and oropharynx and is the sixth-most common cancer worldwide. MicroRNAs are small non-coding RNAs for which altered expression has been demonstrated in pathological processes, such as cancer. The objective of our study was to evaluate the different expression profile in HNSCC subtypes and the prognostic value that one or several miRNAs may have. Materials and Methods: Data from The Cancer Genome Atlas Program-Head and Neck Squamous Cell Carcinoma (TCGA-HNSCC) patients were collected. Differential expression analysis was conducted by edge R-powered TCGAbiolinks R package specific function. Enrichment analysis was developed with Diana Tool miRPath 3.0. Kaplan-Meier survival estimators were used, followed by log-rank tests to compute significance. Results: A total of 127 miRNAs were identified with differential expression level in HNSCC; 48 of them were site-specific and, surprisingly, only miR-383 showed a similar deregulation in all locations studied (tonsil, mouth, floor of mouth, cheek mucosa, lip, tongue, and base of tongue). The most probable affected pathways based on miRNAs interaction levels were protein processing in endoplasmic reticulum, proteoglycans in cancer (p < 0.01), Hippo signaling pathway (p < 0.01), and Transforming growth factor-beta (TGF-beta) signaling pathway (p < 0.01). The survival analysis highlighted 38 differentially expressed miRNAs as prognostic biomarkers. The miRNAs with a greater association between poor prognosis and altered expression (p < 0.001) were miR-137, miR-125b-2, miR-26c, and miR-1304. Conclusions: In this study we have determined miR-137, miR-125b-2, miR-26c, and miR-1304 as novel powerful prognosis biomarkers. Furthermore, we have depicted the miRNAs expression patterns in tumor patients compared with normal subjects using the TCGA-HNSCC cohort.

1. Introduction

Head and Neck Squamous Cell Carcinoma (HNSCC) includes cancers from the oral cavity, larynx, and oropharynx, and is the sixth-most common cancer worldwide [1]. Despite technological and biological advances, prognosis for patients diagnosed with HNSCC is still low, and its incidence has shown an increasing trend, especially in developed countries [2]. HNSCC therapy implies a combination of surgery, radiotherapy, and/or chemotherapy depending on the cancer stage. However, five-year survival rate has not improved in recent years, remaining at less than 50% for patients diagnosed with HNSCC [3]. Many studies have demonstrated gene expression profiles associated with HNSCC formation and progression [4,5].
MicroRNAs (miRNAs) are small non-coding RNAs (18–25 nucleotides) that may target also non-protein-coding genes [6]. MiRNA expression levels varies in different physiological processes, such as cell development, proliferation, differentiation, apoptosis, immune response, and angiogenesis [7]. In pathological processes, such as cancer, changes in the expression of miRNAs have been observed. This has aroused interest in how the deregulated expression of one or more miRNAs can determine the prognosis and/or diagnosis of cancer-affected patients [8]. This has arisen great interest to understand how altered miRNA expression may impact the prognosis of HNSCC and/or aid the diagnosis of these cancer-affected patients. So far, few studies have profiled miRNAs expression on large clinically-annotated oral cavity cohorts (GSE3524: https://pubmed.ncbi.nlm.nih.gov/15381369/; GSE2280: https://pubmed.ncbi.nlm.nih.gov/15558013/; GSE31056: https://pubmed.ncbi.nlm.nih.gov/21989116/; GSE30784: https://pubmed.ncbi.nlm.nih.gov/18669583/) [9]. Moreover, miRNA expression patterns across different studies may vary considerably due to different patients’ clinical-pathological characteristics (sample type, site of origin), and technology used. Novel non-invasive biomarkers such as miRNAs endowed with prognostic or diagnostic potential are necessary for better HNSCC management [10,11].
The identification of the deregulated miRNAs is the first step to design individualized therapies, with reduced side effects, with miRNAs. Currently, two miRNAs are in the clinical phase for the treatment of liver cancer and hepatitis C, miR-34 and miR-122 (Miravirsen), respectively [12,13]. However, there are not yet miRNAs in any clinical phase for the treatment of oral cancer.
The Cancer Genome Atlas (TCGA) is a database managed by the National Cancer Institute and the National Human Genome Research Institute that provides information about genomic, epigenomic, transcriptomic, and proteomic expression profiles of patients who have or have had different types of cancer (including HNSCC) [14]. New analyses are needed to assess the underlying molecular mechanisms altered in HNSCC. It is necessary to find reliable markers to predict patients’ risk and to improve their prognosis, selecting an appropriate treatment.
The objective of our study is to evaluate the different expression profile in HNSCC subtypes and the prognostic value that one or several miRNAs may have using the TCGA database. A secondary objective is to assess the underlying molecular mechanisms in which these miRNAs may be involved.

2. Materials and Methods

Data from TCGA-HNSCC patients were collected in a database that was specifically designed for this purpose, with repeated verification. This work was developed following the Strengthening the Reporting of Observational studies in Epidemiology (STROBE) guide recommendations [15]. The clinical variables used were age, sex, tobacco consumption, alcohol consumption, tumor stage, localization, vital status (alive or deceased), date of death. This article does not contain any studies with human participants or animals performed by any of the authors.
TCGA-HNSC cohort patient and sample data were collected form the Genomic Data Commons (GDC) server. This work was developed following the STROBE guide recommendations [15]. The clinical variables available were age, sex, tobacco consumption, alcohol consumption, tumor stage, localization, vital status (alive or deceased), and date of death. All are primary tumors without any treatment.

2.1. MiRNAs Differential Expression Analysis (DEA)

TCGA-HNSC RNAseq raw count matrix was downloaded from GDC server using TCGAbiolinks R package (version 4.0) [16]. Then, the ComBat-seq function (sva R package) was applied to remove batch effects due to Tissue Source Site (TSS) and plate. The corrected count matrix was filtered and only the 25% top expressed miRNAs were passed to further steps. Next, counts were within and between lane normalized following EDASeq R package instructions. Finally, a TCGAbiolinks edgeR-powered function was applied to determine differentially expressed miRNAs when comparing (1) whole TCGA-HNSC dataset tumors vs normal samples and (2) site-specific resected tumors vs normal samples, as well. This step had four sub-processes powered by EdgeR package: (1) it converted the count filtered matrix into an edgeR DGElist object, (2) it estimated the common dispersion and each miRNA was assigned the same dispersion estimate, (3) it performed the exact test pair-wise for differential expression between the two groups (tumor and normal) and, finally, the process returned log2 (FC), log counts per million (logCPM), p-value, and False discovery rate (FDR) adjusted p-value for every differentially expressed gene. In this particular step, we established a |Log2 (FC)| > 1 and FDR < 0.01 as cutoffs.
Enrichment analysis for molecular pathways was performed in transcriptionally deregulated miRNAs in TCGA-HNSC cohort. Kyoto Encyclopedia of Gene and Genomes (KEGG) analysis was carried out through Diana Tools mirPath v3.0 (free access). They were calculated by Fisher’s exact test and the p-value adjusted by FDR is reported [17].

2.2. Survival Analysis

Based on normalized RNAseq counts of each differentially expressed miRNA, patients were stratified in two groups of high and low expression. The stratification followed a percentile scale-up method, ranging from percentile 2 to 98. In each iteration, the associated p-value was calculated using Kaplan-Meier survival estimators and the log-rank tests. For every miRNA, the percentile associated with the lowest p-value is reported, since it is the point in which survival differences are maximized between both groups. Due to the fact that TCGA-HNSC tumors are very heterogeneous (i.e., different resection sites) and to report potential transversal prognostic biomarkers to the whole cohort, we decided to prioritize those miRNAs with a lowest p-value associated percentile nearby the median. This would reduce the possibility to report as a potential prognostic biomarker an miRNA particularly associated with one resection site or any other variable.

3. Results

A total of 528 samples were included from TCGA-HNSCC data base. Clinical data have been previously described elsewhere [14]. The majority of the samples were of race catalogued as “white” (N = 452, 85.6%), followed by “black or African American race” (N = 48, 9.1%). In terms of ethnicity, most were classified as “non-Hispanic or Latino” (N = 465, 88.1%), followed by “Hispanic or Latino” (N = 26, 4.9%). A total of 127 miRNAs were identified with differential expression level in HNSCC. On the other hand, site-specific differential expression analysis (DEA) reached out 63 Differential expressed miRNAs (DEMs) and only 1 (hsa-miR-383) was found to be deregulated in all locations studied (downregulated in tonsil, mouth, floor of mouth, cheek mucosa, tongue, base of tongue; upregulated in LOP lesions (lip cancer)).
In the DianaTools tool, 63 down-regulated and 92 up-regulated miRNAs were entered. The p-value threshold was 0.05, but the top candidates have been reported. miRNA target genes (by Diana Tools) allowed to gain their possible functional roles through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis in downregulated miRNAs revealed that protein processing in endoplasmic reticulum (p < 0.001), proteoglycans in cancer (p < 0.01), and the Hippo signaling pathway (p < 0.01) were the most probable pathways based on miRNAs interaction levels. In upregulated miRNAs, most probable pathways were the TGF-beta signaling pathway (p < 0.01) and Hippo signaling pathway (p < 0.01) (Figure 1 and Figure S1 Supplementary Materials).
We found 127 DEMs when comparing all TCGA-HNSC tumors vs normal samples, from which only 48 were also found during site-specific DEA. This last analysis reached out 63 transcriptionally altered miRNAs among the 7 analyzed locations (Figure 2A). We see clearly that the tonsil and LOP lesions have a large number of high-level (|log2 (FC)| > 4) expressed miRNAs that are not shared with any location. Furthermore, the percentage of deregulated miRNAs (under- or over-expressed) is also different in the different areas affected. Mouth, floor of mouth, tongue, and base of tongue samples showed a larger number of repressed miRNAs than upregulated, whereas LOP lesions and tonsil samples showed the opposite. Cheek mucosa showed the lowest proportion of deregulated miRNAs among all the studied resection sites.
Previously identified differentially-expressed miRNAs were submitted to survival analysis in order to assess their potential as HNSC prognostic biomarkers. To be selected, a miRNA should match the following criteria: (1) lower survival groups should show higher miRNA levels when it is upregulated, (2) lower survival group should show lower miRNA levels when it is downregulated, (3) the largest significant miRNA-associated K-M curve should be nearby the 50th quantile (median) and, (4) the K-M curve is associated with p-value < 0.01.
Our results show that 38 out of the 127 identified DEMs significantly associated with TCGA-HNSC patients’ survival. From them, 20 upregulated and 6 downregulated follow 1 and 2 criteria, as well. Finally, only hsa-miR-137, hsa-miR-125b-2, hsa-miR-26b, and hsa-miR-1304 show the most significant association in a median-nearby quantile (Figure 3A,B).

4. Discussion

A total of 528 patients were included in this study. Although HNSCC subtypes show molecular heterogeneity and it seems difficult to find a single prognostic signature, some studies have reported individual biomarkers for the different HNSCC subtypes [18].
In 2016, Nathan et al. identified four miRNA signatures to predict the overall survival in oral cancer. Among these miRNAs, miR-26b and miR-142 were positively associated with survival. In our analysis we have associated lower expression of both miRNAs with poorer prognosis, according with Nathan et al. [19]. However, we demonstrated that miR-26b, but not miR-142, is differentially expressed in HNSCC and could be applied as a biomarker regardless of tumor resection site.
Enrichment analysis in downregulated and upregulated miRNAs revealed that protein processing in endoplasmic reticulum, proteoglycans in cancer, the Hippo signaling pathway, and the TGF-beta signaling pathway are the most probable pathways based on miRNAs interaction levels. This goes in accordance with our previous analysis of miRNA expression in oral cancer [20], and other enrichment analysis [21,22]. The alteration of the Hippo signaling mechanism has been reported in different types of cancer, and the altered expression of some of its components has been associated with cell migration and invasion. Its role in HNSCC has been demonstrated in vitro by numerous authors [23]. The mechanism of TGF-beta signaling is responsible for controlling cell proliferation, angiogenesis, and immune function of epithelial cells. The proteins of this pathway often show alterations in the expression in different malignant tumors, including HNSCC [24]. Proteoglycans, such as perlecan, heparin sulphate, or sulphatase 2, in the tumor microenvironment have been previously demonstrated to participate in proliferation, angiogenesis, and metastasis [25].
The association between miRNAs and their targets, as well as the molecular mechanisms in which they are involved, is complicated, since the same miRNA can influence the expression of several proteins, and in turn, the same protein can be regulated by different miRNAs. This may explain the fact that the results between enrichment analysis studies are varied [26].
The use of databases such as TCGA has allowed progress in cancer research, helping to explore molecular mechanisms affected in cancer and to propose new treatment strategies [27]. Although theoretically different locations (for example, tonsil and tongue) in the oral cavity are exposed to the same risk factors, such as tobacco, alcohol, or Human papillomavirus (HPV), clinically these tumors behave differently (in terms of recurrence and survival) [28] and the results of our analysis also demonstrate a different biological alteration in terms of miRNA expression profile. This highlights the fact that perhaps the approach based on therapeutic targets should also be different, or exclusive to certain locations, although initially by proximity they may seem susceptible to the same treatment plan.
The results obtained in relation to miRNA-signature in the HNSCC are similar to other bioinformatics and meta-analysis studies, both in the type of miRNA (such as miR-210 or miR-375) and in its expression (over- or under-expressed) [18,29]. We must emphasize that miR-21 has been proposed as a prognostic marker on numerous occasions [30,31,32,33,34], especially in squamous cell carcinoma of the tongue, but it has not been highlighted in our analysis or other previous bioinformatics [22]. Other very well studied miRNAs, such as miR-375, whose down expression has always been associated with a poor prognosis and has been proposed as a marker in the HNSCC, have been reflected in our analysis [33,35,36,37] due to its significant downregulation, exclusively.
Regarding miRNAs’ over- or under-expressed proportions, in previous microarray studies [20], we have obtained a greater number of under-expressed than over-expressed miRNAs, and in the results of this study we see that in some locations this trend is maintained (such as in the mouth, floor of the mouth, or tongue), but in other locations the opposite is true (such as in the tonsil, cheek mucosa, lip, or base of the tongue). In other bioinformatics studies we also confirmed this trend [21]. We cannot forget, however, that in cancer we have other mechanisms of genetic regulation, such as mutations, and epigenetics, such as methylations [38] or the expression of other non-coding RNAs such as lncRNAs (long non coding RNAs) [39] or snoRNAs (small nucleolar RNAs) [40], that influence the formation and progression of the HNSCC.
Among the limitations of this study include the fact that it was an in silico analysis and that our results must be validated with the most significant deregulated miRNAs to determine the prognostic value. The results obtained in this study should be confirmed when future new cohorts of these tumor types with a size, composition, and analyzed with a technology that must be equal to or superior to that which the TCGA has employed. Also, further analysis onto different samples and conditions, so different datasets, with data from different sources (RNA sequencing and microarray), would be truly beneficial to underpin our results.
On the other hand, the validation of these miRNAs must be carried out with functional studies, in cell lines and animal models. In the case of under-expressed miRNAs, the restitution of the expression of these miRNAs by plasmids or similar vehicles as should demonstrate a decrease in proliferation, migration, or tumor size. In over-expressed microRNAs, blocking their expression by siRNAs or similar vehicles [41]. Another limitation of this study is that the presence of HPV has not been taken into account in the analysis, since the positivity carried out among the samples was low (n = 36).
This study provides a global view of the differential expression profile of HNSCC with a considerable sample size associated with patient survival. The identification of this profile allows the design of treatment strategies and clinical studies to determine the predictive value of miRNAs.

5. Conclusions

In this study, we found eight miRNAs that were commonly deregulated in the same way, in all the localizations of HNSCC. We also determined miR-3689f and miR-142 as two powerful prognosis biomarkers. The determination of miRNA expression profile and its association with prognosis provides useful information for treatment strategies and the development of new clinical studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/1648-9144/56/10/535/s1, Figure S1. Chondroitin sulfate/Dermatan Sulfate (CSPG/DSPG) proteoglycans molecular pathway and miRNAs deregulated in this study. Proteins target by, at least, one of the 8 deregulated miRNAs are colored in yellow.

Author Contributions

C.M.C.P. and A.G.G. contributed to conception of the article; E.P.I. and A.I.L.P. contributed to its design; and B.R.M. and M.P.S. contributed to the interpretation and analysis of data. All authors have read and agreed to the published version of the manuscript.

Funding

No funding sources.

Acknowledgments

We would like to thank Rubén Fernández Caloto for his involvement in this work and his contribution to the analysis of the results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alsahafi, E.; Begg, K.; Amelio, I.; Raulf, N.; Lucarelli, P.; Sauter, T.; Tavassoli, M. Clinical update on head and neck cancer: Molecular biology and ongoing challenges. Cell Death Dis. 2019, 10, 540. [Google Scholar] [PubMed] [Green Version]
  2. Bray, F.; Me, J.F.; Soerjomataram, I.; Siegel, R.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA A Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Taghavi, N.; Yazdi, I. Prognostic factors of survival rate in oral squamous cell carcinoma: Clinical, histologic, genetic and molecular concepts. Arch. Iran. Med. 2015, 18, 314–319. [Google Scholar] [PubMed]
  4. Tuna, M.; Amos, C.I.; Mills, G.B. Genome-Wide Analysis of Head and Neck Squamous Cell Carcinomas Reveals HPV, TP53, Smoking and Alcohol-Related Allele-Based Acquired Uniparental Disomy Genomic Alterations. Neoplasia 2019, 21, 197–205. [Google Scholar] [CrossRef] [PubMed]
  5. Leemans, C.R.; Snijders, P.J.F.; Brakenhoff, R.H. The molecular landscape of head and neck cancer. Nat. Rev. Cancer 2018, 18, 269–282. [Google Scholar]
  6. Mohr, A.M.; Mott, J.L. Overview of MicroRNA Biology. Semin. Liver Dis. 2015, 35, 3–11. [Google Scholar] [CrossRef] [Green Version]
  7. Sayáns, M.P.; Pilar, G.-D.; Barros-Angueira, F.; Suárez-Peñaranda, J.M.; Fernández, A.C.; Gándara-Rey, J.M.; García, A.G. Current trends in miRNAs and their relationship with oral squamous cell carcinoma. J. Oral Pathol. Med. 2011, 41, 433–443. [Google Scholar] [CrossRef]
  8. Guo, W.; Chen, X.; Zhu, L.; Wang, Q. A six-mRNA signature model for the prognosis of head and neck squamous cell carcinoma. Oncotarget 2017, 8, 94528–94538. [Google Scholar]
  9. Zeljic, K.; Jovanovic, I.; Jovanovic, J.; Magic, Z.; Stankovic, A.; Supic, G. MicroRNA meta-signature of oral cancer: Evidence from a meta-analysis. Upsala J. Med Sci. 2018, 123, 43–49. [Google Scholar] [CrossRef]
  10. Janiszewska, J.; Szaumkessel, M.; Szyfter, K. microRNAs are important players in head and neck carcinoma: A review. Crit. Rev. Oncol. 2013, 88, 716–728. [Google Scholar] [CrossRef]
  11. Min, A.; Zhu, C.; Peng, S.; Rajthala, S.; Costea, D.E.; Sapkota, D. MicroRNAs as Important Players and Biomarkers in Oral Carcinogenesis. BioMed Res. Int. 2015, 2015, 1–10. [Google Scholar] [CrossRef] [Green Version]
  12. Van Der Ree, M.H.; Van Der Meer, A.J.; Van Nuenen, A.C.; De Bruijne, J.; Ottosen, S.; Janssen, H.L.; Kootstra, N.A.; Reesink, H.W. Miravirsen dosing in chronic hepatitis C patients results in decreased microRNA-122 levels without affecting other microRNAs in plasma. Aliment. Pharmacol. Ther. 2015, 43, 102–113. [Google Scholar] [CrossRef] [Green Version]
  13. Beg, M.S.; Brenner, A.J.; Sachdev, J.; Borad, M.; Kang, Y.-K.; Stoudemire, J.; Smith, S.; Bader, A.G.; Kim, S.; Hong, D.S. Phase I study of MRX34, a liposomal miR-34a mimic, administered twice weekly in patients with advanced solid tumors. Investig. New Drugs 2016, 35, 180–188. [Google Scholar]
  14. Pérez-Sayáns, M.; Chamorro-Petronacci, C.M.; Lorenzo-Pouso, A.I.; Iruegas, E.P.; Carrión, A.B.; Peñaranda, J.M.S.; García-García, A. Comprehensive Genomic Review of TCGA Head and Neck Squamous Cell Carcinomas (HNSCC). J. Clin. Med. 2019, 8, 1896. [Google Scholar] [CrossRef] [Green Version]
  15. Von Elm, E.; Altman, D.G.; Initiative, F.T.S.; Pocock, S.J.; Gøtzsche, P.C.; Vandenbroucke, J.P. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: Guidelines for reporting observational studies. J. Clin. Epidemiol. 2008, 61, 344–349. [Google Scholar] [CrossRef] [Green Version]
  16. Colaprico, A.; Silva, T.C.; Olsen, C.; Garofano, L.; Cava, C.; Garolini, D.; Sabedot, T.S.; Malta, T.M.; Pagnotta, S.M.; Castiglioni, I.; et al. TCGAbiolinks: An R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2015, 44, e71. [Google Scholar] [CrossRef] [PubMed]
  17. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-miRPath v3.0: Deciphering microRNA function with experimental support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef] [PubMed]
  18. Jamali, Z.; Aminabadi, N.A.; Attaran, R.; Pournagiazar, F.; Oskouei, S.G.; Ahmadpour, F. MicroRNAs as prognostic molecular signatures in human head and neck squamous cell carcinoma: A systematic review and meta-analysis. Oral Oncol. 2015, 51, 321–331. [Google Scholar] [CrossRef] [PubMed]
  19. Nathan, C.A.; Stern, W. 2016 Oral, Head and Neck Cancer Awareness Week: April 10–16. Ear Nose Throat J. 2016, 95, 52–55. [Google Scholar]
  20. Petronacci, C.M.C.; Pérez-Sayáns, M.; Iruegas, M.E.P.; Peñaranda, J.M.S.; Pouso, A.I.L.; Carrión, A.B.; García, A.G. miRNAs expression of oral squamous cell carcinoma patients. Medicine 2019, 98, e14922. [Google Scholar] [CrossRef]
  21. Amiri-Dashatan, N.; Koushki, M.; Jalilian, A.; Ahmadi, N.A.; Tavirani, M.R. Integrated Bioinformatics Analysis of mRNAs and miRNAs Identified Potential Biomarkers of Oral Squamous Cell Carcinoma. Asian Pac. J. Cancer Prev. 2020, 21, 1841–1848. [Google Scholar] [CrossRef] [PubMed]
  22. Wong, N.; Khwaja, S.S.; Baker, C.M.; Gay, H.A.; Thorstad, W.L.; Daly, M.D.; Lewis, J.S.; Wang, X. Prognostic microRNA signatures derived from The Cancer Genome Atlas for head and neck squamous cell carcinomas. Cancer Med. 2016, 5, 1619–1628. [Google Scholar] [CrossRef] [PubMed]
  23. Santos-De-Frutos, K.; Segrelles, C.; Lorz, C. Hippo Pathway and YAP Signaling Alterations in Squamous Cancer of the Head and Neck. J. Clin. Med. 2019, 8, 2131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. White, R.A.; Malkoski, S.P.; Wang, X.-J. TGFβ signaling in head and neck squamous cell carcinoma. Oncogene 2010, 29, 5437–5446. [Google Scholar] [PubMed] [Green Version]
  25. Flowers, S.A.; Zhou, X.; Wu, J.; Wang, Y.; Makambi, K.; Kallakury, B.V.; Singer, M.S.; Rosen, S.D.; Davidson, B.; Goldman, R. Expression of the extracellular sulfatase SULF2 is associated with squamous cell carcinoma of the head and neck. Oncotarget 2016, 7, 43177–43187. [Google Scholar]
  26. Solomon, M.C.; Radhakrishnan, R.A. MicroRNA’s—The vibrant performers in the oral cancer scenario. Jpn. Dent. Sci. Rev. 2020, 56, 85–89. [Google Scholar] [CrossRef]
  27. Sarode, S.; Sarode, S.; Maniyar, N.; Anand, R.; Patil, S. Oral cancer databases: A comprehensive review. J. Oral Pathol. Med. 2017, 47, 547–556. [Google Scholar]
  28. Caponigro, F.; Ionna, F.; Scarpati, G.D.V.; Longo, F.; Addeo, R.; Manzo, R.; Muto, P.; Pisconti, S.; Leopaldi, L.; Perri, F. Translational Research: A Future Strategy for Managing Squamous Cell Carcinoma of the Head and Neck? Anti-Cancer Agents Med. Chem. 2019, 18, 1220–1227. [Google Scholar] [CrossRef]
  29. Lin, N.; Lin, Y.; Fu, X.; Wu, C.; Xu, J.; Cui, Z.; Lin, D. MicroRNAs as a Novel Class of Diagnostic Biomarkers in Detection of Oral Carcinoma: A Meta-Analysis Study. Clin Lab. 2016, 62, 451–461. [Google Scholar]
  30. Danielsson, K.; Wahlin, Y.B.; Gu, X.; Boldrup, L.; Nylander, K. Altered expression of miR-21, miR-125b, and miR-203 indicates a role for these microRNAs in oral lichen planus. J. Oral Pathol. Med. 2011, 41, 90–95. [Google Scholar] [CrossRef]
  31. Manikandan, M.; Rao, A.K.D.M.; Rajkumar, K.S.; Rajaraman, R.; Munirajan, A.K. Altered levels of miR-21, miR-125b-2*, miR-138, miR-155, miR-184, and miR-205 in oral squamous cell carcinoma and association with clinicopathological characteristics. J. Oral Pathol. Med. 2014, 44, 792–800. [Google Scholar] [CrossRef] [PubMed]
  32. Odar, K.; Boštjančič, E.; Gale, N.; Glavač, D.; Zidar, N. Differential expression of microRNAs miR-21, miR-31, miR-203, miR-125a-5p and miR-125b and proteins PTEN and p63 in verrucous carcinoma of the head and neck. Histopathology 2012, 61, 257–265. [Google Scholar] [CrossRef] [PubMed]
  33. He, Q.; Chen, Z.; Cabay, R.J.; Zhang, L.; Luan, X.; Chen, D.; Yu, T.; Wang, A.; Zhou, X. microRNA-21 and microRNA-375 from oral cytology as biomarkers for oral tongue cancer detection. Oral Oncol. 2016, 57, 15–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Zheng, J.; Xue, H.; Wang, T.; Jiang, Y.; Liu, B.; Li, J.; Liu, Y.; Wang, W.; Zhang, B.; Sun, M. miR-21 downregulates the tumor suppressor P12CDK2AP1 and Stimulates Cell Proliferation and Invasion. J. Cell. Biochem. 2011, 112, 872–880. [Google Scholar] [CrossRef] [PubMed]
  35. Cen, W.-N.; Pang, J.-S.; Huang, J.-C.; Hou, J.; Bao, W.-G.; He, R.; Ma, J.; Peng, Z.-G.; Hu, X.-H.; Ma, F.-C. The expression and biological information analysis of miR-375-3p in head and neck squamous cell carcinoma based on 1825 samples from GEO, TCGA, and peer-reviewed publications. Pathol. Res. Pr. 2018, 214, 1835–1847. [Google Scholar] [CrossRef]
  36. Wang, J.; Huang, H.; Wang, C.; Liu, X.; Hu, F.; Liu, M. MicroRNA-375 sensitizes tumour necrosis factor-alpha (TNF-α)-induced apoptosis in head and neck squamous cell carcinoma in vitro. Int. J. Oral Maxillofac. Surg. 2013, 42, 949–955. [Google Scholar] [CrossRef]
  37. Hu, A.; Huang, J.J.; Xu, W.H.; Jin, X.J.; Li, J.P.; Tang, Y.J.; Huang, X.-F.; Cui, H.-C.; Sun, G.-B.; Li, R.-L.; et al. MiR-21/miR-375 ratio is an independent prognostic factor in patients with laryngeal squamous cell carcinoma. Am. J. Cancer Res. 2015, 5, 1775–1785. [Google Scholar]
  38. Sannigrahi, M.; Sharma, R.; Singh, V.; Panda, N.K.; Rattan, V.; Khullar, M. DNA methylation regulated microRNAs in HPV-16-induced head and neck squamous cell carcinoma (HNSCC). Mol. Cell. Biochem. 2018, 448, 321–333. [Google Scholar] [CrossRef]
  39. Affinito, O.; Pane, K.; Smaldone, G.; Orlandella, F.M.; Mirabelli, P.; Beneduce, G.; Parasole, R.; Ripaldi, M.; Salvatore, M.; Franzese, M. lncRNAs–mRNAs Co–Expression Network Underlying Childhood B–Cell Acute Lymphoblastic Leukaemia: A Pilot Study. Cancers 2020, 12, 2489. [Google Scholar] [CrossRef]
  40. Petronacci, C.M.C.; Sayáns, M.P.; Padín-Iruegas, M.E.; Marichalar-Mendia, X.; Torreira, M.G.; García, A.G. Differential expression of snoRNAs in oral squamous cell carcinomas: New potential diagnostic markers. J. Enzym. Inhib. Med. Chem. 2018, 33, 424–427. [Google Scholar] [CrossRef] [Green Version]
  41. Wagner, E. Tumor-targeted Delivery of Anti-microRNA for Cancer Therapy: pHLIP is Key. Angew. Chem. Int. Ed. 2015, 54, 5824–5826. [Google Scholar] [CrossRef]
Medicina 56 00535 g001 550
Figure 1. Most relevant molecular pathways revealed in Diana Tools results for over-expressed miRNAs (N = 92) and under-expressed (N = 63) in Head and Neck Squamous Cell Carcinoma (HNSCC) analyzed samples. Analysis was performed by Fisher’s exact test and the p-value adjusted by False discovery rate (FDR). For the value of p (x-axis), its logarithm has been previously calculated to represent its significance.
Figure 1. Most relevant molecular pathways revealed in Diana Tools results for over-expressed miRNAs (N = 92) and under-expressed (N = 63) in Head and Neck Squamous Cell Carcinoma (HNSCC) analyzed samples. Analysis was performed by Fisher’s exact test and the p-value adjusted by False discovery rate (FDR). For the value of p (x-axis), its logarithm has been previously calculated to represent its significance.
Medicina 56 00535 g001
Medicina 56 00535 g002 550
Figure 2. Differential expression of miRNAs in HNSCC subtypes. Panel (A) hierarchical cluster. The number of over-expressed, under-expressed, and no variation in expression can be observed by comparing the cases with healthy controls. The left column indicates the origin of the samples. Panel (B) represents the differences in miRNA expression by location in percentage. Under-expressed miRNAs are represented in blue and over-expressed in red. Those miRNAs without differential expression are represented by the color grey.
Figure 2. Differential expression of miRNAs in HNSCC subtypes. Panel (A) hierarchical cluster. The number of over-expressed, under-expressed, and no variation in expression can be observed by comparing the cases with healthy controls. The left column indicates the origin of the samples. Panel (B) represents the differences in miRNA expression by location in percentage. Under-expressed miRNAs are represented in blue and over-expressed in red. Those miRNAs without differential expression are represented by the color grey.
Medicina 56 00535 g002
Medicina 56 00535 g003 550
Figure 3. Evaluation of differently expressed miRNAs and prognosis. Panel (A): colored squares indicate miRNA expression level (red: upregulated; blue: downregulated). The intensity of green squares represents how significant the survival curves are. The circle shows the differential survival with which they are associated: blue—less expression, worse prognosis; red—more expression, worse prognosis. * Represents those microRNAs up- or down-regulated in HNSCC and whose higher or lower expression, respectively, is associated with worse survival. Numbers are the quantile at which each miRNA shows the largest survival association. Panel (B): Kaplan–Meier survival curves (log-rank test P) represent the survival of HNSCC patients in those four miRNAs with the highest level of prognostic significance according to the described criteria. MiR-137: N low expression = 242, N high expression = 275; miR-26b: N low expression = 274, N high expression = 243; miR-1340: N low expression = 314, N high expression = 203; miR-125b-2: N low expression = 316, N high expression = 201.
Figure 3. Evaluation of differently expressed miRNAs and prognosis. Panel (A): colored squares indicate miRNA expression level (red: upregulated; blue: downregulated). The intensity of green squares represents how significant the survival curves are. The circle shows the differential survival with which they are associated: blue—less expression, worse prognosis; red—more expression, worse prognosis. * Represents those microRNAs up- or down-regulated in HNSCC and whose higher or lower expression, respectively, is associated with worse survival. Numbers are the quantile at which each miRNA shows the largest survival association. Panel (B): Kaplan–Meier survival curves (log-rank test P) represent the survival of HNSCC patients in those four miRNAs with the highest level of prognostic significance according to the described criteria. MiR-137: N low expression = 242, N high expression = 275; miR-26b: N low expression = 274, N high expression = 243; miR-1340: N low expression = 314, N high expression = 203; miR-125b-2: N low expression = 316, N high expression = 201.
Medicina 56 00535 g003

Share and Cite

MDPI and ACS Style

Chamorro Petronacci, C.M.; García García, A.; Padín Iruegas, E.; Rivas Mundiña, B.; Lorenzo Pouso, A.I.; Pérez Sayáns, M. Identification of Prognosis Associated microRNAs in HNSCC Subtypes Based on TCGA Dataset. Medicina 2020, 56, 535. https://doi.org/10.3390/medicina56100535

AMA Style

Chamorro Petronacci CM, García García A, Padín Iruegas E, Rivas Mundiña B, Lorenzo Pouso AI, Pérez Sayáns M. Identification of Prognosis Associated microRNAs in HNSCC Subtypes Based on TCGA Dataset. Medicina. 2020; 56(10):535. https://doi.org/10.3390/medicina56100535

Chicago/Turabian Style

Chamorro Petronacci, Cintia M., Abel García García, Elena Padín Iruegas, Berta Rivas Mundiña, Alejandro I. Lorenzo Pouso, and Mario Pérez Sayáns. 2020. "Identification of Prognosis Associated microRNAs in HNSCC Subtypes Based on TCGA Dataset" Medicina 56, no. 10: 535. https://doi.org/10.3390/medicina56100535

Article Metrics

Back to TopTop