Computational and Mathematical Methods in Medicine
Volume 2017 (2017), Article ID 9803018, 8 pages
https://doi.org/10.1155/2017/9803018

Research Article

Biomarker MicroRNAs for Diagnosis of Oral Squamous Cell Carcinoma Identified Based on Gene Expression Data and MicroRNA-mRNA Network Analysis

Hui Zhang,¹ Tangxin Li,² Linqing Zheng,¹ and Xiangya Huang¹

¹Guanghua School of Stomatology, Affiliated Stomatological Hospital, Guangdong Province Key Laboratory of Stomatology, Guangdong, China
²Oral and Maxillofacial Center, Kiang Wu Hospital, Macau

Correspondence should be addressed to Xiangya Huang; hxiangya@mail.sysu.edu.cn

Received 22 May 2017; Revised 1 August 2017; Accepted 6 August 2017; Published 17 September 2017

Academic Editor: Xiaofeng Song

Copyright © 2017 Hui Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Oral squamous cell carcinoma is one of the most malignant tumors with high mortality rate worldwide. Biomarker discovery is critical for early diagnosis and precision treatment of this disease. MicroRNAs are small noncoding RNA molecules which often regulate essential biological processes and are good candidates for biomarkers. By integrative analysis of both the cancer-associated gene expression data and microRNA-mRNA network, miR-148b-3p, miR-629-3p, miR-27a-3p, and miR-142-3p were screened as novel diagnostic biomarkers for oral squamous cell carcinoma based on their unique regulatory abilities in the network structure of the conditional microRNA-mRNA network and their important functions. These findings were confirmed by literature verification and functional enrichment analysis. Future experimental validation is expected for the further investigation of their molecular mechanisms.

1. Introduction

Oral squamous cell carcinoma (OSCC) is the sixth most common cancer with more than 300,000 cases worldwide each year [1]. It is the most malignant tumor in the oral and maxillofacial regions and accounts for 90% of oral cancers [2, 3]. The risk factors for OSCC could be tobacco, alcohol consumption, betel quid (BQ) chewing, Bidi smoking, and genetic predisposition [4, 5]. OSCC can metastasize to lymph-node, even to remote organs with high mortality rate. The present diagnosis of OSCC often happened at late stage and the treatment can be unsuccessful due to its local recurrence. The precise early diagnosis is critical and essential to the future prevention and personalized treatment of this disease.

MicroRNA is a family of functional noncoding RNA molecules containing about 22 nucleotides, which play roles in the posttranscriptional gene regulation. Since many key biological processes including the development, differentiation, and cell cycles are regulated by microRNAs, the abnormal expression of microRNAs is often associated with the initialization and progression of many diseases [6, 7]. Thus miRNAs usually could serve as suitable biomarkers for many diseases, such as neurodevelopmental disorders [8], cancer, and cardiovascular disease [9–11].

Previous studies have demonstrated that microRNAs played important roles in OSCC. For example, microRNA-23b/27b cluster is reported as tumor suppressive and regulates the MET oncogene in OSCC [12]. MicroRNA-27a-3p can regulate transition from epithelial to mesenchymal in OSCC by targeting YAP1 [13]. The apoptosis-related protein expression and radiosensitivity in BQ-associated OSCC are regulated by microRNA-17-5p [14]. Metabolic shift in OSCC is mediated by microRNA-340 targeting glucose transporter-1 [15]. Tumor growth and activation of NF-κB signaling were promoted via the regulation of NLK by microRNA-92b in OSCC [16]. In addition, microRNA-17/20a was suggested as a prognostic marker since it can inhibit cell migration in OSCC [2]. Circulating microRNA-21 and PTEN (phosphatase and tensin homolog) are reported as promising biomarkers for detection of OSCC [17].

From the above introduction, we believe that microRNAs are good candidates to act as diagnostic and prognostic biomarkers of OSCC. As we know, OSCC is a complex and heterogeneous disease. For this reason, more precise and personalized biomarkers are needed for the diagnosis, prognosis, and treatment of OSCC. Until now, very few studies have focused on the expression data of OSCC to integrate it with the microRNA-mRNA network structural analysis for biomarker discovery in OSCC, especially the application of bioinformatics and network analysis to the study of the functions of microRNAs in the OSCC initialization and progression.

The experimental methods for biomarker discovery are time-consuming and costly. Bioinformatics screening will be helpful to the efficient biomarker screening. Previously, several models have been developed to infer key and biomarker microRNAs in complex diseases from conditional gene expression data. Differential expression genes (DEGs) are often used to screen biomarker genes, but only few DEGs are validated as biomarkers; therefore integrative analysis of DEGs with other information is very necessary for efficient biomarker discovery. As described in previous work [18–21], these models screen the potential biomarker based on the scrutinizing of the structure of the conditional microRNA-mRNA network. By statistical analysis of the network structure and functions of the biomarker microRNA’s targets, the model can very effectively identify novel putative microRNA biomarkers for the diagnosis of complex diseases. So we here apply the model to the biomarker microRNA discovery for diagnosis of OSCC.

2. Materials and Methods

The schematic pipeline of the present work for the data collection, model construction, biomarker microRNAs prediction, and validation and enrichment analysis of the targets of the predicted microRNAs is presented in Figure 1. The details of the step-by-step procedures for the screening of OSCC diagnostic biomarker microRNAs are described as follows.

Figure 1: Schematic pipeline for the identification of oral squamous cell carcinoma (OSCC) microRNA biomarkers.

2.1. Gene Expression Data Collection

The data for the OSCC gene expression and microRNA expression were extracted from the GEO database [22]. The OSCC associated expression data in the GEO database are listed in Table 1. Eight OSCC associated data sets measured from different microarray platforms are available in the GEO database. After the condition filtering, the final data used for the construction of OSCC-specific microRNA-mRNA network are GSE30784 and GSE28100. The former is the mRNA expression data including 167 OSCC samples and 45 samples as control [23] and the latter is the microRNA expression data with 17 OSCC samples and 3 control samples [24]. The data were normalized and the differentially expressed mRNAs were identified based on linear models in Limma R package [25, 26]. The value and other parameters were calculated with the empirical Bayes (eBayes) method. The Benjamini-Hochberg correction was applied to adjust the values. The adjusted value less than 0.05 was chosen as the cut-off criteria.

Table 1: OSCC gene expression data collected from GEO database.

Accession/ID
PMID
Platform
OSCC
Control
Gene/microRNA

GSE3524
15381369
GPL96

Gene

GSE70604
26700817
GPL2986

Gene

23624915

GSE37991
25204733
GPL6883

Gene

23362108

GSE28100
22761427
GPL10850
n = 17
n = 3
MicroRNA

GSE23558
22072328
GPL6480

Gene

GSE25099
21853135
GPL5175

Gene

GSE30784
18669583
GPL570
n = 167
n = 45
Gene

GSE10121
18472963
GPL6353

Gene

The reported OSCC associated microRNAs were also collected from PubMed citations by the search criteria “(Oral squamous cell carcinoma OR OSCC) AND (miRNA OR microRNA) AND (biomarker OR marker)”. They were checked manually and listed in Table 2.

Table 2: Previously reported microRNA biomarkers for OSCC from PubMed.

MicroRNA
PMID
Biomarker Type
Samples
Expression in OSCC
NOD
(Ref.)
TFP
(Ref.)

miR-155-5p
26307116
24692283
Prognostic biomarker
Tissue
Upregulated
71
0.211

miR-483-5p
26224475
Diagnostic/prognostic biomarker
Serum
Upregulated
4
0.135

miR-216a
25955794
Prognostic biomarker
Tissue
Downregulated
6
0.087

miR-21-3p 
 miR-96-5p 
 miR-141-3p 
 miR-130b-3p
25532855
Prognostic biomarker
Tissue
Upregulated
N/A
8
22
2
N/A
0.159
0.167
0.119

miR-21-5p
24755828
Prognostic biomarker
Tissue
Upregulated
38
0.135

miR-31-5p
22083872
Diagnostic biomarker
Saliva
Upregulated
10
0.142

The NOD and TFP are calculated based on the reference human microRNA-mRNA network, while the measurements in Section 3.2 are calculated based on the OSCC-specific network.

2.2. Prediction of microRNA Biomarkers for Diagnosis of OSCC

As reported in the previous researches in Shen’s group [18–21], two measurements are important for candidate biomarker microRNA. The first one is the novel of degree measurement (NOD). It measures the number of genes solely targeted by certain microRNA [19, 20]. This character is reasonable since the abnormal alteration of this type of interaction cannot be compensated by another microRNA-mRNA interaction pair as most of the microRNA-mRNA interactions are synergic. The other measurement is the transcription factor percentage (TFP), which was defined as the percentage of transcription factor (TF) genes of all targets of the microRNA [18].

According to the above hypothesis, the OSCC-specific microRNA-mRNA network was constructed by mapping the detected differentially expressed microRNAs in OSCC onto the reference human microRNA-mRNA network. The reference network was constructed with the integration of the experimentally validated and computational predicted microRNA-mRNA pairs. The experimentally validated data included information from miRecords [27], TarBase [28], miR2Disease [29], and miRTarBase [30], while the computational predicted microRNA-mRNA pairs are extracted from no fewer than 2 databases among HOCTAR [31], ExprTargetDB [32], and starBase [33]. With the reconstructed conditional network, the above-mentioned measurements, that is, the NOD and TFP, were calculated for each microRNA in the OSCC-specific network. MicroRNAs with significantly large NOD and TFP values (Wilcoxon signed-rank test, value < 0.05) were screened as putative biomarkers.

2.3. Functional Enrichment Investigation of the Targets of Predicted OSCC Diagnostic MicroRNA Biomarkers

Functional enrichment analysis of the genes targeted by the identified candidate biomarker microRNAs from the OSCC-specific microRNA-mRNA networks was performed through three different tools: Gene Ontology Annotation, KEGG Pathway Analysis, and Ingenuity Pathway Analysis (IPA). Here, the first two analyses were conducted on the DAVID (Database for Annotation, Visualization, and Integrated Discovery) online analysis webpage [34]. The significantly enriched pathways and ontologies for OSCC with value less than 0.05 were ranked. The value was calculated based on the hypergeometric test and FDR adjustment was used for multiple test correction.

3. Result and Discussion

3.1. The Characterization of the Previous Reported OSCC Diagnostic Biomarker MicroRNAs

We checked the PubMed citations and the previously reported biomarker microRNAs for OSCC were listed in Table 2. From Table 2, it is clear that all the reported microRNAs have high NOD and TFP values except miR-21-3p, which cannot be extracted from the reconstructed microRNA-mRNA network. This observation confirmed that the model using the NOD and TFP as two measurements for the evaluation of the potential biomarkers is applicable for OSCC biomarker discovery.

3.2. Predicted Diagnostic Biomarker MicroRNAs for OSCC

We first identified 56 dysregulated microRNAs and 3375 differentially expressed genes in OSCC by using the pipeline presented in Figure 1. Five microRNAs were identified through Wilcoxon signed-rank test with value less than 0.05. These microRNAs were predicted to be candidate biomarkers for the diagnosis of OSCC. Their network structural characteristics in the microRNA-mRNA network, including the number of targets and NOD and TFP values calculated based on the conditional OSCC-specific microRNA-mRNA network, are listed in Table 3. Among the five microRNAs, miR-155-5p was reported in previous work as biomarker [35, 36]. The other four microRNAs, that is, miR-148b-3p, miR-629-3p, miR-27a-3p, and miR-142-3p, are the novel putative biomarkers identified for OSCC.

Table 3: Candidate OSCC microRNA biomarkers identified by our model.

MicroRNA
NOD
value
TFP
value
Number of targets

miR-148b-3p
15
5.06E-06
0.168
0.00912
95

miR-155-5p
27
1.25E − 12
0.259
1.53E − 05
58

miR-629-3p
10
0.00845
0.214

14

miR-27a-3p
10
0.00845
0.194

31

miR-142-3p
10
0.00845
0.164
0.0168
55

3.3. Literature-Based Validation of Identified MicroRNA Biomarkers

The targets of the five putative microRNAs are presented in Figure 2. From the figure, we can see that some of the targets of these microRNAs have been reported to be associated with OSCC (genes colored red) or other oral diseases (genes colored yellow) according to the PubMed citations.

Figure 2: Candidate oral squamous cell carcinoma (OSCC) microRNA biomarkers with their target genes. Here, ellipses in black represent 5 microRNA biomarkers. Genes in red are associated with OSCC according to literature reports while genes in yellow are reported to be associated with other oral diseases.

3.4. Functional Enrichment Analysis of Target Genes of Candidate MicroRNA Biomarkers

The functional enrichment analysis was further performed to explore the roles of the uniquely regulated genes of the identified microRNAs in OSCC with DAVID and IPA tools. Figure 3 presented the Gene Ontology (GO) Annotation for targets of the identified microRNA biomarkers in OSCC. The three levels of GO analysis are presented in Figures 3(a)–3(c), respectively, for biological process, cellular component, and molecular function. The top 10 significantly enriched items are listed for each level. Most of the dysregulated biological processes are the positive/negative biological or cellular processes, the regulation of cell cycle, and the response to oxygen-containing compound. The former are well-known popular cancer-associated processes, while the latter is associated specifically with OSCC [37–39]. The most enriched molecular functions are the general cancer-associated items, such as protein binding, protein kinase activity, and receptor signaling protein activity. The functions for carbohydrate binding [40–42] and glycosaminoglycan binding [43, 44] were also discovered in the OSCC studies.

(a)

(b)

(c)

Figure 3: Gene Ontology (GO) Annotation for genes targeted by identified microRNA biomarkers. (a), (b), and (c), respectively, represent three levels of GO: biological process, cellular component, and molecular function. The statistical significance value ( value) has been negative 10-based log-transformed. The top 10 significantly enriched items are listed for each level.

The result of the pathway enrichment analysis of the targeted genes of the putative microRNA biomarkers is displayed in Figure 4. The most common cancer-associated pathways like p53 signaling pathway and cell cycle pathway are enriched in both the DAVID and IPA methods. There are still other pathways such as PI3K-Akt signaling pathways and colorectal cancer metastasis were screened by these two enrichment analyses. The Aryl hydrocarbon receptor [45–47], the HGF [48, 49], ECM receptor interaction [50, 51], Hepatitis B [52, 53], and glucocorticoid receptor signaling [54] are all supported by the PubMed citations.

(a)

(b)

Figure 4: KEGG Pathway Enrichment Analysis and Ingenuity Pathway Analysis for genes targeted by identified microRNA biomarkers. The statistical significance value ( value) has been negative 10-based log-transformed. The top 10 significantly enriched pathways are listed, respectively, in (a) and (b).

4. Conclusions

In this research, we applied the concepts of NOD and TPF to the integrative analysis of OSCC gene expression and the microRNA-mRNA network. We identified five microRNAs that could be putative biomarkers for OSCC. Among them, one has been reported as biomarker and two are reported as associated microRNAs. The other two are the novel finding microRNA biomarkers. As a result, four novel biomarker microRNAs, that is, miR-148b-3p, miR-629-3p, miR-27a-3p, and miR-142-3p, are discovered in our work. The literature checking and the functional enrichment analysis confirmed our finding. Therefore, further experimental verification and clinical testing were suggested for these putative OSCC microRNA biomarkers.

Abbreviations

OSCC: | Oral squamous cell carcinoma
TF: | Transcription factor
NOD: | Novel out degree
TFP: | Transcription factor gene percentage
KEGG: | Kyoto Encyclopedia of Genes and Genomes
DAVID: | Database for Annotation, Visualization, and Integrated Discovery
IPA: | Ingenuity Pathway Analysis.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Xiangya Huang and Linqing Zheng conceived and designed the experiments; Hui Zhang performed the experiments; Hui Zhang and Tangxin Li analyzed the data; Xiangya Huang wrote the paper. Hui Zhang and Tangxin Li contributed equally to this work.

References

Q. Wang, P. Gao, X. Wang, and Y. Duan, “Investigation and identification of potential biomarkers in human saliva for the early diagnosis of oral squamous cell carcinoma,” Clinica Chimica Acta, vol. 427, pp. 79–85, 2014.

C.-C. Chang, Y.-J. Yang, Y.-J. Li et al., “MicroRNA-17/20a functions to inhibit cell migration and can be used a prognostic marker in oral squamous cell carcinoma,” Oral Oncology, vol. 49, no. 9, pp. 923–931, 2013.

Y. Li, B. Li, B. Xu et al., “Expression of p53, p21(CIP1/WAF1) and eIF4E in the adjacent tissues of oral squamous cell carcinoma: establishing the molecular boundary and a cancer progression model,” International Journal of Oral Science, vol. 7, no. 3, pp. 161–168, 2015.

R. Subapriya, A. Thangavelu, B. Mathavan, C. R. Ramachandran, and S. Nagini, “Assessment of risk factors for oral squamous cell carcinoma in Chidambaram, southern India: a case-control study,” European Journal of Cancer Prevention, vol. 16, no. 3, pp. 251–256, 2007.

H.-B. Santos, T.-K. dos Santos, A.-R. Paz et al., “Clinical findings and risk factors to oral squamous cell carcinoma in young patients: A 12-year retrospective analysis,” Medicina Oral, Patologia Oral y Cirugia Bucal, vol. 21, no. 2, pp. e151–e156, 2016.

Y. Tang, W. Yan, J. Chen, C. Luo, A. Kaipia, and B. Shen, “Identification of novel microRNA regulatory pathways associated with heterogeneous prostate cancer,” BMC Systems Biology, vol. 7, supplement 3, article S6, 2013.

Y. W. Kong, D. Ferland-McCollough, T. J. Jackson, and M. Bushell, “microRNAs in cancer management,” The Lancet Oncology, vol. 13, no. 6, pp. e249–e258, 2012.

L. Shen, Y. Lin, Z. Sun, X. Yuan, L. Chen, and B. Shen, “Knowledge-guided bioinformatics model for identifying autism spectrum disorder diagnostic microrna biomarkers,” Scientific Reports, vol. 6, article 39663, 2016.

X.-M. Zhao, K.-Q. Liu, G. Zhu et al., “Identifying cancer-related microRNAs based on gene expression data,” Bioinformatics, vol. 31, no. 8, pp. 1226–1234, 2015.

G.-M. Qin, R.-Y. Li, and X.-M. Zhao, “Identifying Disease Associated miRNAs Based on Protein Domains,” IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 13, no. 6, pp. 1027–1035, 2016.

S. Becker, A. Florian, A. Patrascu et al., “Identification of cardiomyopathy associated circulating miRNA biomarkers in patients with muscular dystrophy using a complementary cardiovascular magnetic resonance and plasma profiling approach,” Journal of Cardiovascular Magnetic Resonance, vol. 18, no. 1, article no. 244, 2016.

I. Fukumoto, K. Koshizuka, T. Hana Zawa et al., “The tumor-suppressive microRNA-23b/27b cluster regulates the MET oncogene in oral squamous cell carcinoma,” International Journal of Oncology, vol. 49, no. 3, pp. 1119–1129, 2016.

G. Zeng, W. Xun, K. Wei, Y. Yang, and H. Shen, “MicroRNA-27a-3p regulates epithelial to mesenchymal transition via targeting YAP1 in oral squamous cell carcinoma cells,” Oncology Reports, vol. 36, no. 3, pp. 1475–1482, 2016.

S.-Y. Wu, A. T. H. Wu, and S.-H. Liu, “Microrna-17-5p regulated apoptosis-related protein expression and radiosensitivity in oral squamous cell carcinoma caused by betel nut chewing,” Oncotarget, vol. 7, no. 32, pp. 51482–51493, 2016.

P. Xu, Y. Li, H. Zhang, M. Li, and H. Zhu, “MicroRNA-340 Mediates Metabolic Shift in Oral Squamous Cell Carcinoma by Targeting Glucose Transporter-1,” Journal of Oral and Maxillofacial Surgery, vol. 74, no. 4, pp. 844–850, 2016.

Z. Liu, C. Diep, T. Mao et al., “MicroRNA-92b promotes tumor growth and activation of NF-κB signaling via regulation of NLK in oral squamous cell carcinoma,” Oncology Reports, vol. 34, no. 6, pp. 2961–2968, 2015.

W. Ren, C. Qiang, L. Gao et al., “Circulating microRNA-21 (MIR-21) and phosphatase and tensin homolog (PTEN) are promising novel biomarkers for detection of oral squamous cell carcinoma,” Biomarkers, vol. 19, no. 7, pp. 590–596, 2014.

W. Yan, L. Xu, Z. Sun et al., “MicroRNA biomarker identification for pediatric acute myeloid leukemia based on a novel bioinformatics model,” Oncotarget, vol. 6, no. 28, pp. 26424–26436, 2015.

W. Zhang, J. Zang, X. Jing et al., “Identification of candidate miRNA biomarkers from miRNA regulatory network with application to prostate cancer,” Journal of Translational Medicine, vol. 12, article 66, 2014.

J. Zhu, S. Wang, W. Zhang et al., “Screening key microRNAs for castration-resistant prostate cancer based on miRNA/mRNA functional synergistic network,” Oncotarget, vol. 6, no. 41, pp. 43819–43830, 2015.

J. Chen, D. Zhang, W. Zhang et al., “Clear cell renal cell carcinoma associated microRNA expression signatures identified by an integrated bioinformatics analysis,” Journal of Translational Medicine, vol. 11, no. 1, article 169, 2013.

T. Barrett, T. O. Suzek, D. B. Troup et al., “NCBI GEO: mining millions of expression profiles—database and tools,” Nucleic Acids Research, vol. 33, pp. D562–D566, 2005.

C. Chen, E. Méndez, J. Houck et al., “Gene expression profiling identifies genes predictive of oral squamous cell carcinoma,” Cancer Epidemiology Biomarkers and Prevention, vol. 17, no. 8, pp. 2152–2162, 2008.

H. M. Jung, B. L. Phillips, R. S. Patel et al., “Keratinization-associated miR-7 and miR-21 regulate tumor suppressor reversion-inducing cysteine-rich protein with kazal motifs (RECK) in oral cancer,” The Journal of Biological Chemistry, vol. 287, no. 35, pp. 29261–29272, 2012.

G. K. Smyth, “Linear models and empirical Bayes methods for assessing differential expression in microarray experiments,” Statistical Applications in Genetics and Molecular Biology, vol. 3, Art. 3, 29 pages, 2004.

M. E. Ritchie, B. Phipson, D. Wu et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, 2015.

F. Xiao, Z. Zuo, G. Cai, S. Kang, X. Gao, and T. Li, “miRecords: an integrated resource for microRNA-target interactions,” Nucleic Acids Research, vol. 37, no. 1, pp. D105–D110, 2009.

P. Sethupathy, B. Corda, and A. G. Hatzigeorgiou, “TarBase: a comprehensive database of experimentally supported animal microRNA targets,” RNA, vol. 12, no. 2, pp. 192–197, 2006.

Q. Jiang, Y. Wang, Y. Hao et al., “miR2Disease: a manually curated database for microRNA deregulation in human disease,” Nucleic Acids Research, vol. 37, no. 1, pp. D98–D104, 2009.

S.-D. Hsu, F.-M. Lin, W.-Y. Wu et al., “MiRTarBase: a database curates experimentally validated microRNA-target interactions,” Nucleic Acids Research, vol. 39, supplement 1, pp. D163–D169, 2011.

V. A. Gennarino, M. Sardiello, R. Avellino et al., “MicroRNA target prediction by expression analysis of host genes,” Genome Research, vol. 19, no. 3, pp. 481–490, 2009.

E. R. Gamazon, H.-K. Im, S. Duan et al., “ExprTarget: an integrative approach to predicting human microRNA targets,” PLoS ONE, vol. 5, no. 10, Article ID e13534, 2010.

J.-H. Yang, J.-H. Li, P. Shao, H. Zhou, Y.-Q. Chen, and L.-H. Qu, “StarBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data,” Nucleic Acids Research, vol. 39, no. 1, pp. D202–D209, 2011.

G. Dennis Jr. et al., “DAVID: database for annotation, visualization, and integrated discovery,” Genome Biology, vol. 4, no. 5, 2003.

M. I. Rather, M. N. Nagashri, S. S. Swamy, K. S. Gopinath, and A. Kumar, “Oncogenic microRNA-155 down-regulates tumor suppressor CDC73 and promotes oral squamous cell carcinoma cell proliferation: implications for cancer therapeutics,” Journal of Biological Chemistry, vol. 288, no. 1, pp. 608–618, 2013.

O. Baba, S. Hasegawa, H. Nagai et al., “MicroRNA-155-5p is associated with oral squamous cell carcinoma metastasis and poor prognosis,” Journal of Oral Pathology and Medicine, vol. 45, no. 4, pp. 248–255, 2016.

E. Sasabe, Z. Yang, S. Ohno, and T. Yamamoto, “Reactive oxygen species produced by the knockdown of manganese-superoxide dismutase up-regulate hypoxia-inducible factor-1α expression in oral squamous cell carcinoma cells,” Free Radical Biology and Medicine, vol. 48, no. 10, pp. 1321–1329, 2010.

K.-F. Hua, P.-C. Liao, Z. Fang et al., “Generation of Reactive Oxygen Species by Polyenylpyrroles Derivatives Causes DNA Damage Leading to G2/M Arrest and Apoptosis in Human Oral Squamous Cell Carcinoma Cells,” PLoS ONE, vol. 8, no. 6, Article ID e67603, 2013.

M. Czesnikiewicz-Guzik, B. Lorkowska, J. Zapala et al., “NADPH oxidase and uncoupled nitric oxide synthase are major sources of reactive oxygen species in oral squamous cell carcinoma. Potential implications for immune regulation in high oxidative stress conditions,” Journal of Physiology and Pharmacology, vol. 59, no. 1, pp. 139–152, 2008.

P. Li, L. Y. Xiao, and H. Tan, “Muc-1 promotes migration and invasion of oral squamous cell carcinoma cells via PI3K-Akt signaling,” International Journal of Clinical and Experimental Pathology, vol. 8, no. 9, pp. 10365–10374, 2015.

Y.-M. Ding, J.-H. Dong, L.-L. Chen, and H.-D. Zhang, “Increased expression of galectin-1 is associated with human oral squamous cell carcinoma development,” Oncology Reports, vol. 21, no. 4, pp. 983–987, 2009.

S. Aggarwal, S. C. Sharma, and S. N. Das, “Galectin-1 and galectin-3: Plausible tumour markers for oral squamous cell carcinoma and suitable targets for screening high-risk population,” Clinica Chimica Acta, vol. 442, pp. 13–21, 2015.

K. Ueda, “Heparin induces apoptosis through suppression of AKt in oral squamous cell carcinoma cells,” Anticancer Research, vol. 29, no. 4, pp. 1079–1088, 2009.

S. Mori, M. Nose, H. Morikawa et al., “A novel evaluation system of metastatic potential of oral squamous cell carcinoma according to the histopathological and histochemical grading,” Oral Oncology, vol. 34, no. 6, pp. 549–557, 1998.

Y. Kimura, A. Kasamatsu, D. Nakashima et al., “ARNT2 regulates tumoral growth in oral squamous cell carcinoma,” Journal of Cancer, vol. 7, no. 6, pp. 702–710, 2016.

E. A. Stanford, A. Ramirez-Cardenas, Z. Wang et al., “Role for the aryl hydrocarbon receptor and diverse ligands in oral squamous cell carcinoma migration and tumorigenesis,” Molecular Cancer Research, vol. 14, no. 8, pp. 696–706, 2016.

X. Xie, Y. Jiang, Y. Yuan et al., “MALDI imaging reveals NCOA7 as a potential biomarker in oral squamous cell carcinoma arising from oral submucous fibrosis,” Oncotarget, vol. 7, no. 37, pp. 59987–60004, 2016.

S. Klosek, K. Nakashiro, S. Hara, C. Li, S. Shintani, and H. Hamakawa, “Constitutive activation of Stat3 correlates with increased expression of the c-Met/HGF receptor in oral squamous cell carcinoma,” Oncology Reports, 2004.

H. Borgiel-Marek, B. Marek, I. Niedzielska, and D. Kajdaniuk, “Serum concentration of fibroblastic growth factor in oral squamous cell carcinoma before and after surgery,” International Journal of Oral and Maxillofacial Surgery, vol. 38, no. 5, p. 537, 2009.

Z. Liu, Y. Niu, C. Li, Y. Yang, and C. Gao, “Integrating multiple microarray datasets on oral squamous cell carcinoma to reveal dysregulated networks,” Head and Neck, vol. 34, no. 12, pp. 1789–1797, 2012.

Y. He, F. Shao, W. Pi et al., “Largescale transcriptomics analysis suggests over-expression of BGH3, MMP9 and PDIA3 in oral squamous cell carcinoma,” PLoS ONE, vol. 11, no. 1, Article ID e0146530, 2016.

M. M. Bakri, H. M. Hussaini, A. Holmes, R. D. Cannon, and A. M. Rich, “Revisiting the association between candidal infection and carcinoma, particularly oral squamous cell carcinoma,” Journal of Oral Microbiology, vol. 2, no. 2010, article no. 5780, 2010.

F.-M. Fang, C.-F. Li, H.-Y. Huang et al., “Overexpression of a chromatin remodeling factor, RSF-1/HBXAP, correlates with aggressive oral squamous cell carcinoma,” American Journal of Pathology, vol. 178, no. 5, pp. 2407–2415, 2011.

B. Li et al., “[Assessment of mental health status in oral squamous cell carcinoma patients and its correlation with catecholamines level],” Shanghai Kou Qiang Yi Xue, vol. 22, no. 6, pp. 671–675, 2013.