Skip to main content

Bioinformatics analysis of the expression and role of microRNA-221-3p in head and neck squamous cell carcinoma



Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer worldwide, associated with a high rate of morbidity and mortality. However, the target genes of miR-221-3p and the underlying mechanism involved in HNSCC are still not clear. Therefore, in the current study, we studied the role of miR-221-3p in the HNSCC.


Tissues collected from 48 control and 21 HNSCC patients were processed to check the differential expression of miR-221-3p by RT-qPCR. Overexpression of microRNA-221-3p (miR-221-3p) is significantly correlated to the onset and progression of HNSCC. We also conducted the meta-analysis of the cancer literature from the cancer genome atlas (TCGA) and the Gene Expression Omnibus (GEO) database to estimate the expression of miR-221-3p in HNSCC. The miR-221-3p target genes in the HNSCC were predicted with the miRWalk and TCGA databases, and functionally annotated via the Gene Ontology. Finally, Spearman’s analysis was used to determine the role of the related target genes in important pathways involved in the development of HNSCC.


We observed a significantly higher expression of miR-221-3p in HNSCC compared to the normal with a summary receiver operating characteristic (sROC) of 0.86(95% Cl: 0.83,0.89). The KEGG and GO comprehensive analysis predicted that miR-221-3p might be involved in the development of HNSCC through the following metabolic pathways, viz. Drug metabolism - cytochrome P450 UGT1A7 and MAOB may be important genes for the role of miR-221-3p.


Based on bioinformatics analysis, our results indicate that miR-221-3p may be used as a non-invasive and hypersensitive biomarker in the diagnosis. Thus, it can be concluded that miR-221-3p may be an extremely important gene locus involved in the process of the deterioration and eventual tumorigenesis of HNSCC. Hopefully, additional work will validate its usefulness as a target for future clinical research.

Peer Review reports


Head and neck squamous cell carcinoma (HNSCC), encompassing oral squamous cell carcinoma, oropharyngeal squamous cell carcinoma, hypopharyngeal squamous cell carcinoma and laryngeal squamous cell carcinoma, is the sixth most common cancer worldwide, accounting for approximately 1–2% of all cancer deaths, with approximately 600,000 new cases detected each year globally [1,2,3]. The treatment of patients at the early stage is relatively successful, but the overall survival rate of recurrent or metastatic HNSCC remains low and has barely seen any improvements for decades [4, 5]. Although the past three decades have seen several improvements in diagnostic tools and treatment regimens, the overall survival rate for advanced (stage III-IV) HNSCC is only 65% [6, 7]. The current standard of care for such patients is surgery along with radiation and chemotherapy, but these have not significantly improved the 5-year survival of HNSCC patients. The main reason for the decline in patient survival is the lack of effective therapeutic targets for the development of HNSCC [8, 9].

Recent years have seen an increasing concentration of research on microRNAs (miRNAs) [10, 11], which are a type of highly conserved single-stranded noncoding RNA containing 17 ~ 22 nucleotides and are involved in the process of tumorigenesis, cell survival, and chemosensitivity [12,13,14]. They bind to the 3′ untranslated region (3’UTR) of different target mRNA (messenger RNA) genes to degrade or inhibit the mRNA translation of target genes associated with a tumor suppressor function [13, 15, 16]. In normal healthy individuals, miR-221-3p is observed to play a role in the process of vascular proliferation [17], while the tumor promoter microRNA-221 is involved in the process of regulating apoptosis of tumor cells [18,19,20] and is associated with a variety of cancer types, including hepatocellular cancer [21], cutaneous melanoma [19], prostate cancer [20], and non-small cell lung cancer [22].

Studies have shown that specific miRNA profiles can be identified between tumor tissue and adjacent healthy tissue in HNSCC patients [2, 23]. For instance, studies have reported a link between miR-221 and vascular invasion in HNSCC [14, 15, 24,25,26]. However, the definite target gene of miR-221-3p and its biological mechanism of action are still not clear. Thus, in the current study, we investigated the expression of miR-221-3p in HNSCC and attempted to explore the correlation between the two.



Tissue specimens were collected from 48 control and 21 HNSCC patients from the Department of Pathology, First Affiliated Hospital of Guangxi Medical University (Nanning, Guangxi Medical University). Total RNA was extracted from formalin-fixed, paraffin-embedded (FFPE) tissues using an FFPE RNA kit (Omega Bio-Tek, Norcross, GA, USA). RNA was reverse-transcribed into cDNA using the Mir-X miRNA qRT-PCR SYBR Kit (Takara, California, USA). The cDNA samples were processed for qPCR with SYBR-Green Master Mix (Takara, Tokyo, Japan) on an ABI 7500 cycler (Applied Biosystems) under the following conditions: initial denaturation at 95 °C for 30 s, followed by 40 cycles of denaturation at 95 °C for 5 s and annealing at 60 °C for 34 s. The expression of miR-221-3p in HNSCC tissues relative to negative control (NC) tissues was calculated using the 2- ΔΔCT method, with U6 as the internal control. The primers for miR-221-3p and U6 were synthesized by TaKaRa (Dalian, Liaoning, China), and the sequences were as follows: miR-221-3p: forward 5′-AGCUACAUUGUCUGCUGGGUUUC − 3′ and reverse 5′ mRQ 3′; U6: forward 5′-GGAACGATACAGAGAAGATTAGC-3′ and reverse 5′-TGGAACGCTTCACGAATTTGCG-3′. All the experiments were repeated three times.

Literature search and selection strategy

A literature search was conducted in PubMed, Chinese Biomedical Literature Database (CBM), Science Direct, China National Knowledge Infrastructure (CNKI) database, Web of Science, Wiley Online Library, EMBASE, China Science and Technology Journal Database (VIP) and Wanfang Database for this study. We searched the databases from the earliest available data to October 1, 2019. The following keywords were used: (HNSCC OR SCC OR “squamous cell cancer” OR “squamous cell carcinoma”) AND (oropharynx OR oropharyngeal OR “head and neck” OR nose OR nasopharynx OR “nasal sinus” OR “nasal cavity” OR “oral cavity” OR hypopharynx OR oral OR laryngopharynx OR larynx OR laryngopharyngeal OR laryngeal OR pharyngeal OR tongue OR tonsil OR tonsillar OR cheek OR palatal OR “paranasal sinuses” OR buccal OR lip) AND (microRNA-221-3p OR miRNA-221-3p OR “miR 221-3p” OR “miRNA 221-3p” OR miRNA221-3p OR miR221-3p).

Selection criteria and data extraction

The databases were searched independently by two researchers who selected the studies based on the following inclusion criteria: (1) comparison of HNSCC and noncancerous tissues; (2) validation of miR-221-3p expression levels via reverse transcription quantitative PCR (RT-qPCR); (3) evaluation of the association between miR-221-3p expression and clinical outcomes; and (4) availability of sufficient data to calculate the mean, standard deviation (SD) and 95% confidence intervals (95% CIs). The articles were eliminated if they met any of the following exclusion criteria: (1) irrelevant to the research focus; (2) inclusion of unqualified data; (3) publication language other than English or Chinese; (4) overlapping or duplicate publications; and (5) letters, reviews, comments, editorials, conference articles, laboratory studies or case reports. The reviewers independently appraised the quality of data in each eligible study and extracted the first author name, year of publication, country of origin, sample type, sample size and analysis method. For articles with incomplete information, the authors were contacted to obtain relevant information.

Microarray data collection from GEO

The microarray data of HNSCC samples uploaded until October 1, 2019, were obtained from the GEO database ( using the following keywords: (HNSCC OR SCC OR “squamous cell cancer” OR “squamous cell carcinoma”) AND (oropharynx OR oropharyngeal OR “head and neck” OR nose OR nasopharynx OR “nasal sinus” OR “nasal cavity” OR “oral cavity” OR hypopharynx OR oral OR laryngopharynx OR larynx OR laryngopharyngeal OR laryngeal OR pharyngeal OR tonsil OR tonsillar OR tongue OR cheek OR palatal OR “paranasal sinuses” OR buccal OR lip) AND (microRNA OR miRNA OR “miR” OR “miRNA”). The inclusion criteria for the microarray datasets were as follows: (1) comprising data from HNSCC and noncancerous tissues; (2) evaluation of the association between miR-221-3p expression and clinical outcomes; and (3) availability of sufficient data to calculate the mean, SD and 95% CI. The exclusion criteria were as follows: (1) unrelated to this study; (2) unqualified data; and (3) publication language other than English.

RNA sequencing data selection from TCGA

The miR-221-3p expression data of HNSCC and normal tissues and the clinicopathological parameters of patients were downloaded and extracted from the OncoLnc website (

Predicting the target genes of miR-221-3p in HNSCC

We use the MiRWalk 2.0 [27] ( and GEPIA ( databases to retrieve differentially expressed genes of miR-221-3p in HNSCC. Gene Ontology (GO) [28] analysis, which was used to define the biological processes (BPs), molecular functions (MFs), and cellular components (CCs) of the target genes, and Kyoto Encyclopedia of Genes and Genomes (KEGG) [29] pathway enrichment analyses were conducted using R (version 3.6.1), and the ClusterProfiler package of R was used to visualize the results. The Search Tool for the Retrieval of Interacting Genes (STRING) ( was then used to establish a protein-protein interaction (PPI) network of the miR-221-3p target genes associated with the significantly enriched pathways. The expression levels of the miR-221-3p target genes in HNSCC and nontumor tissues were determined with UALCOULD [30] ( Finally, the LinkedOmics [31] ( Spearman’s analysis tool was applied to determine the correlation between the expression levels of miR-221-3p and the potential target genes involved in notable signaling pathways.

Statistical analysis

Independent-Samples T Test for the RT-qPCR analysis of two groups. SPSS 24.0 software (IBM, Somers, NY) was used for the statistical analysis. The mean ± SD (X ± s) were calculated, and t test was used for comparisons between groups for all the measurement data. Data that were used for our meta-analysis as well as TCGA data (Nanning, Guangxi Zhuang Autonomous Region, China) were analyzed using STATA 15.0 software. The standard mean difference (SMD) and 95% confidence interval were used to estimate the expression value of miR-221-3p. Scatter plots of the data of both the normal and HNSCC tissues were prepared with GraphPad Prism 6 (Nanning, Guangxi Zhuang Autonomous Region, China) and used to assess the expression level of miR-221-3p. Evidence of bias was assessed by visual funnel plots and Egger’s regression asymmetry test. In addition, the I2 index was used to evaluate the potential heterogeneity of the selected data. When I2 > 50% or P < 0.05, we used the fixed effect model; otherwise, the random effect model was used.. The count rate was expressed as a percentage (%), and we used the X2 test for comparisons. The rates of true positives (TPs), false positives (FPs), false negatives (FNs) and true negatives (TNs) were used to determine the diagnostic value of miR-221-3p. We used Meta-disc software to calculate the specificity, sensitivity, negative and positive likelihood ratios, and summary receiver operating characteristic (sROC) curve. We also performed sensitivity analysis to assess the differences between the sample sizes. Moreover, all the miR-221-3p expression data, which included the TCGA sequencing data, were normalized to log2 to improve the normality of the measurements. P < 0.05 was taken as a statistically significant difference for all analyses.


miR-221-3p is highly expressed in the tumor tissues

The RT-qPCR data clearly showed that miR-221-3p is highly expressed in the tissues of nasopharyngeal carcinoma patients compared to the controls (p < 0.05) (Fig.1).

Fig. 1

MiR-221-3p levels in HNSCC and NC tissue specimens. RT-qPCR was used to detect miR-221-3p was significantly up-regulated in the tumor tissues of HNSCC (n = 21) relative to control tissues (n = 48)(p<0.01)

Selection of relevant literature and microarray data extraction

A total of 17 articles were retrieved from the initial search, and the following paper was selected after the full-text review based on the inclusion criteria: Zhou Cheng. et al. (PMID:30928631) [32] (Fig.2). Sixteen out of 424 GEO microarray datasets, including GSE11163, GSE103931, GSE107591, GSE113956, GSE28100, GSE31277, GSE34496, GSE45238, GSE51829, GSE58911, GSE69002, GSE70289, GSE73171, GSE75630, GSE82064, and GSE98463, met our inclusion criteria. Information for each included dataset is summarized in Table 1.

Fig. 2

The flow chart of literature search and selection of relevant studies. Sixteen out of 424 GEO microarray datasets met our inclusion criteria. At the same time, a total of 17 articles were retrieved from the initial search, and one papers (Zhou Cheng. et al) were selected after the full-text review based on the inclusion criteria

Table 1 The main features of included studies for this meta-analysis

Risk of bias assessment and meta-analysis of miR-221-3p in HNSCThe results of the visual funnel plots show a symmetrical shape, which indicates the absence of any significant publication bias as a whole. Sensitivity analysis performed to assess the heterogeneity of the samples showed that there were no significant differences between the studies (Fig.3). There were 911 HNSCC samples and 289 normal samples from the PubMed, GEO, and TCGA datasets. The meta-analysis performed on all data to explore the expression level of miR-221-3p in HNSCC indicated a high degree of heterogeneity between these studies. Therefore, the random-effects model was selected, and the combined standard mean difference (SMD) was observed to be 0.72 (95% CI: 0.28, 1.15) (Fig.4). The meta-analysis data indicate that miR-221-3p expression is significantly upregulated in HNSCC, which is consistent with the results obtained by RT-qPCR. In addition, miR-221-3p expression was significantly increased in HNSCC in the GSE11163, GSE103931, GSE113956, GSE31277, GSE45238, GSE82064, TCGA and PMID30928631 datasets, and the expression levels of miR-221-3p in HNSCC and normal tissues from each included dataset are shown in Fig.5.

Fig. 3

Bias analysis. a Funnel plots of published data bias for meta-analysis of miR-21. b Sensitivity analysis for meta-analysis of miR-221-3p

Fig. 4

Forest plots. The combined SMD of miR-221-3p expression levels in HNSCC samples was 0.717(95% Cl:0.283, 1.152)

Fig. 5

The expression data of miR-221-3p in HNSCC from TCGA, each included GEO and the literature

Diagnostic value of miR-221-3p in HNSCC

The random-effects model used to analyze the diagnostic value of miR-221-3p for HNSCC showed significant heterogeneity in the likelihood ratio (negative and positive) sensitivity and specificity analyses. The diagnostic meta-analysis results indicated that the pooled specificity, sensitivity, likelihood ratio (negative and positive) and diagnostic odds ratio were 0.78 (95% CI: 0.72–0.82), 0.57 (95% CI: 0.54–0.60), 0.42 (95% CI: 0.31–0.56) and 3.04 (95% CI: 1.72–5.40), 9.59(95% CI: 4.77–19.29) respectively (Fig.6). The receiver operating characteristic (ROC) curve of GSE103931, GSE11163, GSE113956, GSE31277, GSE34496, GSE45238, GSE82064, TCGA and PMID30928631 is statistically significant (p < 0.05, Fig.7). Based on the ROC curve of each study, the overall ROC curve indicated that the area under the sROC curve (AUC) was 0.86 (95% Cl: 0.83, 0.89) (Fig.8).

Fig. 6

The diagnostic value of miR-221-3p in HNSCC. a Sensitivity; (b) Specificity; (c) Negative LR; (d) Positive LR; (e) Diagnostic odds ratio

Fig. 7

Receiver operating characteristic (ROC) curves showing the diagnostic value of miR-221-3p in HNSCC

Fig. 8

Summary ROC (SROC) curve showing the diagnostic performance of miR-221-3p in HNSCC. The sROC curve AUC was 0.86 (95% CI: 0.83–0.89), indicating high predictive power

Clinical features of HNSCC based on the TCGA

Based on the HNSCC sample data extracted from the TCGA database, we analyzed 527 samples, and the results indicated the absence of any statistically significant difference between miR-221-3p expression and age, lymphatic invasion, neoplasm histologic grade and any other clinicopathological features. However, the results indicated statistically significant differences in the expression of miR-221-3p in different tissues, sexes, and even alcohol consumption statuses (Table 2).

Table 2 The correlation between miR-221-3p expression levels and clinic characteristics based on TCGA database

The prospective target genes of miR-221-3p in HNSCC

By using the MiRWalk 2.0 and GEPIA databases, we retrieved 5311 and 467 differentially expressed genes, respectively. Bioinformatics analysis showed a total of 117 overlapping genes that were involved (Fig.9, Table 3). These overlapped 117 genes could be considered as the target genes that miR-221-3p might play a role in HNSCC. The enriched GO and KEGG pathway categories of the overlapped 117 genes with p < 0.05 are shown in Fig.10, Table 4 and Table 5. For the cellular components, the identified target genes were mostly enriched in the actin cytoskeleton, sarcomeres, and contractile fiber part; for the molecular functions, the target genes were mainly enriched in aromatase activity, oxidoreductase activity, and cofactor binding. KEGG enrichment analysis showed that miR-221-3p plays a significant role in HNSCC through a variety of pathways, including the drug metabolism pathways of the cytochrome P450 signaling pathway. The drug metabolism-cytochrome P450-related genes include GSTA1, UGT1A7, CYP3A5, FMO2, MAOB, ADH1C and ADH1B, and the protein-protein interaction (PPI) network of these genes is shown in Fig.11.

Fig. 9

Venn diagram showing overlapping miR-221-3p target genes

Table 3 The promising target genes of miR-221-3p in HNSCC with TCGA and MiRWalk
Fig. 10

Enrichment analysis of GO and KEGG pathways of miR-221-3p target genes. a Gene ontology (GO) of candidate genes; (b) Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis

Table 4 Predictive target genes of miR-221-3p with GO analysis
Table 5 KEGG pathway of validated target genes of miR-221-3p
Fig. 11

PPI network. The nodes represent target genes and the lines are associations between genes, which showing the connections among the seven miR-221-3p targets associated with the Drug metabolism-cytochrome P450 pathway

Validation of the miR-221-3p target genes in the drug metabolism-cytochrome P450 signaling pathway

The expression values of each target gene are shown in Fig. 12, from which we learned that the 5 target genes (UGT1A7, CYP3A5, FMO2, MAOB and ADH1B) related to miR-221-3p were downregulated in the drug metabolism-cytochrome P450 pathway (p < 0.05). Spearman correlation analysis showed that GSTA1, UGT1A7 and MAOB, the target genes of the drug metabolism-cytochrome P450, were correlated with miR-221-3p in HNSCC (p < 0.05, Fig.13).

Fig. 12

The expression levels of miR-221-3p target genes in HNSCC. a-g represents the expression of the following genes: GSTA1, UGT1A7, CYP3A5, FMO2, MAOB, ADH1C, ADH1B

Fig. 13

Spearman analysis. a-g show the association between miR-221-3p expression levels and the seven identified target genes involved in the Drug metabolism - cytochrome P450 pathway


The drug metabolism-cytochrome P450 signaling pathway is an essential signal transduction pathway in cells. It is an important biological function of cell survival, proliferation and apoptosis. For example, Liping Wang and his colleagues suggest [33] that the inhibition of CYP3A5 in cytochrome P450 drug metabolism can drive the migration, proliferation, and invasion of HNSCC cells. Dongfang Wang et al. showed that cytochrome P450 is correlated with the overall survival and vascular invasion of hepatocellular carcinoma (HCC) patients. Additionally, inhibiting the P450 pathway in drug metabolism-cytochrome can increase the efficacy of HNSCC [34]. Among the target genes of miR-221-3p, MAOB and UGT1A7 are associated with the miR-221-3p and drug-cytochrome P450 signaling pathway and are downregulated in HNSCC (p < 0.05). MAOB is regarded as a novel biomarker for accurate prostate cancer diagnosis and treatment [35]. In addition, a large number of valuable epidemiological studies have shown that UGT1A7 affects individuals’ susceptibility to various cancers, such as pancreatic cancer [36] and gastrointestinal carcinomas [37]. Our study suggests that miR-221-3p may play an important role as a gene that promotes the development of HNSCC by reducing the expression of the MAOB and UGT1A7 pathways in HNSCC.

Understanding the pathogenesis of HNSCC and identifying new gene therapy programs have been the focus of recent studies. Many recent studies have reported that miRNAs with different expression patterns in different tumors [38,39,40,41,42] control the progression of tumors. For example, miR-221-3p has been regarded as a tumor biomarker that can be used to assess the clinical prognosis of breast cancer [15]. In healthy humans, it has been shown that miR-221-3p plays a role in the process of vascular proliferation, while the tumor promoter miR-221-3p regulates the apoptosis of tumor cells. The current study focused on the expression pattern of miR-221-3p in HNSCC, identifying the exact target genes of this miRNA and its biological mechanism of action through bioinformatics analysis. The biological pathways relevant to the actions of miR-221-3p will guide further understanding of the mechanisms of HNSCC.

Of course, the limitations of our study should also be mentioned. On the one hand, this study is restricted in that the analysis of differential miRNAs was only based on HNSCC and noncancerous tissues, and other samples (e.g., blood) were not assessed. On the other hand, the target genes of miR-221-3p have not been experimentally confirmed, and our conclusions will need to be confirmed by clinical or molecular biological methods in the future. In addition, only English studies were used for the basis of our meta-analysis, which did not include other potentially relevant studies that were in other languages.


In our study, we found that miR-221-3p levels were apparently upregulated in HNSCC compared to normal tissues. MAOB and UGT1A7 are potentially important targets of miR-221-3p. MiR-221-3p may be used as a noninvasive and hypersensitive biomarker for the diagnosis of HNSCC and is an extremely important gene locus involved in the process of the deterioration and eventual tumorigenesis of HNSCC. Finally, the biological pathways relevant to the actions of miR-221-3p will provide insights into its potential molecular mechanisms in HNSCC. Larger-scale studies will be needed to validate its diagnostic promise. Hopefully, additional work will validate its usefulness as a target for future clinical research.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.



Head and neck squamous cell carcinoma




protein-protein interaction


Gene Ontology


Kyoto Encyclopedia of Genes and Genomes


Gene Expression Omnibus


The cancer genome atlas


Database for Annotation, Visualization and Integrated Discovery


Biological Networks Gene Oncology tool


biological process


cellular component


molecular function


  1. 1.

    Avissar M, Christensen BC, Kelsey KT, Marsit CJ. MicroRNA expression ratio is predictive of head and neck squamous cell carcinoma. Clin Cancer Res. 2009;15(8):2850–5.

    Google Scholar 

  2. 2.

    Allen B, Schneider A, Victoria B, Nunez Lopez YO, Muller M, Szewczyk M, Pazdrowski J, Majchrzak E, Barczak W, Golusinski W, Golusinski P, Masternak MM. Blood serum from head and neck squamous cell carcinoma patients induces altered MicroRNA and target gene expression profile in treated cells. Front Oncol. 2018;8:217.

    Google Scholar 

  3. 3.

    Peitzsch C, Nathansen J, Schniewind SI, Schwarz F, Dubrovska A. Cancer stem cells in head and neck squamous cell carcinoma: identification, Characterization and Clinical Implications. Cancers. 2019;11(5):616.

    Google Scholar 

  4. 4.

    Ausoni S, Boscolo-Rizzo P, Singh B, Da Mosto MC, Spinato G, Tirelli G, Spinato R, Azzarello G. Targeting cellular and molecular drivers of head and neck squamous cell carcinoma: current options and emerging perspectives. Cancer Metastasis Rev. 2016;35(3):413–26.

    Google Scholar 

  5. 5.

    Britt EL, Raman S, Leek K, Sheehy CH, Kim SW, Harada H. Combination of fenretinide and ABT-263 induces apoptosis through NOXA for head and neck squamous cell carcinoma treatment. PLoS One. 2019;14(7):e0219398.

    Google Scholar 

  6. 6.

    Koole K, Clausen MJ, van Es RJ, van Kempen PM, Melchers LJ, Koole R, Langendijk JA, van Diest PJ, Roodenburg JL, Schuuring E, et al. FGFR family members protein expression as prognostic markers in Oral cavity and Oropharyngeal squamous cell carcinoma. Molecular diagnosis therapy. 2016;20(4):363–74.

    Google Scholar 

  7. 7.

    Tang XX, He J, Li B, Zheng Y, Li KJ, Zou S, Chen L. Efficacy and safety of Gefitinib in patients with advanced head and neck squamous cell carcinoma: a meta-analysis of randomized controlled trials. J Oncol. 2019;2019:6273438.

    Google Scholar 

  8. 8.

    Cao W, Liu J, Xia R, Lin L, Wang X, Xiao M, Zhang C, Li J, Ji T, Chen W. X-linked FHL1 as a novel therapeutic target for head and neck squamous cell carcinoma. Oncotarget. 2016;7(12):14537–50.

    Google Scholar 

  9. 9.

    Ahmad P, Sana J, Slavik M, Slampa P, Smilek P, Slaby O. MicroRNAs involvement in Radioresistance of head and neck Cancer. Dis Markers. 2017;2017:8245345.

    Google Scholar 

  10. 10.

    Orosz E, Gombos K, Riedling T, Afiakurue P, Kiss I, Pytel J, Gerlinger I, Szanyi I. Comparative miRNA expression profile analysis of squamous cell carcinoma and Peritumoral mucosa from the Meso- and Hypopharynx. Cancer Genomics Proteomics. 2017;14(4):285–92.

    Google Scholar 

  11. 11.

    Qureshi R, Sacan A. A novel method for the normalization of microRNA RT-PCR data. BMC Med Genomics. 2013;6(Suppl 1):S14.

    Google Scholar 

  12. 12.

    Irani S. miRNAs signature in head and neck squamous cell carcinoma metastasis: a literature review. J Dent. 2016;17(2):71–83.

    Google Scholar 

  13. 13.

    Wang H, Tang Y, Yang D, Zheng L. MicroRNA-591 functions as a tumor suppressor in hepatocellular carcinoma by lowering drug resistance through inhibition of far-upstream element-binding protein 2-mediated Phosphoinositide 3-kinase/Akt/mammalian target of Rapamycin Axis. Pharmacology. 2019:1–14.

  14. 14.

    Hussein S, Mosaad H, Rashed HE, El-Anwar MW. Up-regulated miR-221 expression as a molecular diagnostic marker in laryngeal squamous cell carcinoma and its correlation with Apaf-1 expression. Cancer Biomark. 2017;19(3):279–87.

    Google Scholar 

  15. 15.

    Zhou L, Jiang F, Chen X, Liu Z, Ouyang Y, Zhao W, Yu D. Downregulation of miR-221/222 by a microRNA sponge promotes apoptosis in oral squamous cell carcinoma cells through upregulation of PTEN. Oncol Lett. 2016;12(6):4419–26.

    Google Scholar 

  16. 16.

    Deng X, Ma L, Wu M, Zhang G, Jin C, Guo Y, Liu R. miR-124 radiosensitizes human glioma cells by targeting CDK4. J Neurooncol. 2013;114(3):263–74.

    Google Scholar 

  17. 17.

    Rippe C, Blimline M, Magerko KA, Lawson BR, LaRocca TJ, Donato AJ, Seals DR. MicroRNA changes in human arterial endothelial cells with senescence: relation to apoptosis, eNOS and inflammation. Exp Gerontol. 2012;47(1):45–51.

    Google Scholar 

  18. 18.

    Gillies JK, Lorimer IA. Regulation of p27Kip1 by miRNA 221/222 in glioblastoma. Cell Cycle. 2007;6(16):2005–9.

    Google Scholar 

  19. 19.

    Felicetti F, Errico MC, Bottero L, Segnalini P, Stoppacciaro A, Biffoni M, Felli N, Mattia G, Petrini M, Colombo MP, Peschle C, Care A. The promyelocytic leukemia zinc finger-microRNA-221/−222 pathway controls melanoma progression through multiple oncogenic mechanisms. Cancer Res. 2008;68(8):2745–54.

    Google Scholar 

  20. 20.

    Ambs S, Prueitt RL, Yi M, Hudson RS, Howe TM, Petrocca F, Wallace TA, Liu CG, Volinia S, Calin GA, Yfantis HG, Stephens RM, Croce CM. Genomic profiling of microRNA and messenger RNA reveals deregulated microRNA expression in prostate cancer. Cancer Res. 2008;68(15):6162–70.

    Google Scholar 

  21. 21.

    Li J, Wang Y, Yu W, Chen J, Luo J. Expression of serum miR-221 in human hepatocellular carcinoma and its prognostic significance. Biochem Biophys Res Commun. 2011;406(1):70–3.

    Google Scholar 

  22. 22.

    Ergun S, Arman K, Temiz E, Bozgeyik I, Yumrutas O, Safdar M, Dagli H, Arslan A, Oztuzcu S. Expression patterns of miR-221 and its target Caspase-3 in different cancer cell lines. Mol Biol Rep. 2014;41(9):5877–81.

    Google Scholar 

  23. 23.

    Wald AI, Hoskins EE, Wells SI, Ferris RL, Khan SA. Alteration of microRNA profiles in squamous cell carcinoma of the head and neck cell lines by human papillomavirus. Head Neck. 2011;33(4):504–12.

    Google Scholar 

  24. 24.

    Yang CJ, Shen WG, Liu CJ, Chen YW, Lu HH, Tsai MM. Lin SC: miR-221 and miR-222 expression increased the growth and tumorigenesis of oral carcinoma cells. J Oral Pathol Med. 2011;40(7):560–6.

    Google Scholar 

  25. 25.

    Liu CJ, Kao SY, Tu HF, Tsai MM, Chang KW, Lin SC. Increase of microRNA miR-31 level in plasma could be a potential marker of oral cancer. Oral Dis. 2010;16(4):360–4.

    Google Scholar 

  26. 26.

    Wong TS, Liu XB, Wong BY, Ng RW, Yuen AP, Wei WI. Mature miR-184 as potential oncogenic microRNA of squamous cell carcinoma of tongue. Clin Cancer Res. 2008;14(9):2588–92.

    Google Scholar 

  27. 27.

    Dweep H. Gretz N: miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12(8):697.

    Google Scholar 

  28. 28.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25(1):25–9.

    Google Scholar 

  29. 29.

    Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.

    Google Scholar 

  30. 30.

    Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, Varambally S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia. 2017;19(8):649–58.

    Google Scholar 

  31. 31.

    Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–d963.

    Google Scholar 

  32. 32.

    Cheng Z, Xing D. Ginsenoside Rg3 inhibits growth and epithelial-mesenchymal transition of human oral squamous carcinoma cells by down-regulating miR-221. Eur J Pharmacol. 2019;853:353–63.

    Google Scholar 

  33. 33.

    Wang LP, Chen WH, Zha J, Yan YY, Wei YX, Chen XL, Zhu XX, Ge LH. miR543 acts as a novel oncogene in oral squamous cell carcinoma by targeting CYP3A5. Oncol Rep. 2019;42(3):973–90.

    Google Scholar 

  34. 34.

    Ruwali M, Dhawan A, Pant MC, Rahman Q, Khurana SM, Parmar D. Clinical Management of Head and Neck Cancer Cases: role of Pharmacogenetics of CYP2 and GSTs. Oncol Res Treat. 2016;39(4):221–6.

    Google Scholar 

  35. 35.

    Xu N, Wu Y-P, Ke Z-B, Liang Y-C, Cai H, Su W-T, Tao X, Chen S-H, Zheng Q-S, Wei Y, et al. Identification of key DNA methylation-driven genes in prostate adenocarcinoma: an integrative analysis of TCGA methylation data. J Transl Med. 2019:17(1).

  36. 36.

    Yilmaz L, Borazan E, Aytekin T, Baskonus I, Aytekin A, Oztuzcu S, Bozdag Z, Balik A. Increased UGT1A3 and UGT1A7 expression is associated with pancreatic Cancer. Asian Pac J Cancer Prev. 2015;16(4):1651–5.

    Google Scholar 

  37. 37.

    Cengiz B, Yumrutas O, Bozgeyik E, Borazan E, Igci YZ, Bozgeyik I, Oztuzcu S. Differential expression of the UGT1A family of genes in stomach cancer tissues. Tumor Biol. 2015;36(8):5831–7.

    Google Scholar 

  38. 38.

    Chen W, Fu W, Deng Q, Li Y, Wang K, Bai Y, Wu X, Li G, Wang G, Huang J, He M, Zhang X, Wu T, Wei S, Guo H. Multiple metals exposure and chromosome damage: exploring the mediation effects of microRNAs and their potentials in lung carcinogenesis. Environ Int. 2019;122:291–300.

    Google Scholar 

  39. 39.

    Dong XT, Liu YD. Expression and significance of miR-24 and miR-101 in patients with advanced gastric cancer. Oncol Lett. 2018;16(5):5769–74.

    Google Scholar 

  40. 40.

    Jin FF, Yang R, Wei Y, Wang D, Zhu YN, Wang XH, Lu YS, Wang YB, Zen K, Li LM. HIF-1alpha-induced miR-23a approximately 27a approximately 24 cluster promotes colorectal cancer progression via reprogramming metabolism. Cancer Lett. 2019;440-441:211–22.

    Google Scholar 

  41. 41.

    Jing PY, Zhao N, Xie NL, Ye MX, Zhang Y, Zhang ZP, Li MY, Lai XF, Zhang J, Gu ZP. miR-24-3p/FGFR3 Signaling as a novel Axis is involved in epithelial-Mesenchymal transition and regulates lung adenocarcinoma progression. J Immunol Res. 2018;2018:2834109.

    Google Scholar 

  42. 42.

    Lang B, Shang C, Meng LR. Targeted silencing of S100A8 gene by miR-24 to increase chemotherapy sensitivity of endometrial carcinoma cells to paclitaxel. Med Sci Monit. 2016;22:1953–8.

    Google Scholar 

Download references


Not applicable.


This work was supported by grants from the National Natural Science Foundation of China (No. 81460460,81760542), The Research Foundation of the Science and Technology Department of Guangxi Province, China (grant No. 2016GXNSFAA380252, 2018AB61001 and 2014GXNSFBA118114), the Research Foundation of the Health Department of Guangxi Province, China (No. S2018087). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing this manuscript.

Author information




MK designed the study and accessed the relevant information, YZZ designed the study, analyzed data, and wrote the manuscript. And YZZ and LWW collected and analyzed and interpreted the data. XJL, XZX and CL were involved in statistical analysis. QQL, JMS, XRQ and SS critically revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Min Kang.

Ethics declarations

Ethics approval and consent to participate

All procedures were approved by the Ethics Committee of Guangxi Medical University (Nanning, China). Written informed consent was obtained from all patients or their families prior to enrolment in the present study.

Consent for publication

Not applicable.

Competing interests

The authors declared that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhou, Z., Wu, W., Li, J. et al. Bioinformatics analysis of the expression and role of microRNA-221-3p in head and neck squamous cell carcinoma. BMC Cancer 21, 395 (2021).

Download citation


  • miR-221-3p
  • RT-qPCR
  • TCGA
  • Bioinformatics