miRNA biomarkers for predicting overall survival outcomes for head and neck squamous cell carcinoma

Single-cell transcriptomics reveals heterogeneity in esophageal squamous epithelial cells and constructs models for predicting patient prognosis and immunotherapy

Background

Esophageal squamous cell carcinoma (ESCC), characterized by its high invasiveness and malignant potential, has long been a formidable challenge in terms of treatment.

Methods

A variety of advanced analytical techniques are employed, including single-cell RNA sequencing (scRNA-seq), cell trajectory inference, transcription factor regulatory network analysis, GSVA enrichment analysis, mutation profile construction, and the inference of potential immunotherapeutic drugs. The purpose is to conduct a more comprehensive exploration of the heterogeneity among malignant squamous epithelial cell subgroups within the ESCC microenvironment and establish a model for predicting the prognosis and immunotherapy outcomes of ESCC patients.

Results

An analysis was conducted through scRNA-seq, and three Cluster of malignant epithelial cells were identified using the infer CNV method. Cluster 0 was found to exhibit high invasiveness, whereas Cluster 1 displayed prominent characteristics associated with epithelial-mesenchymal transition. Confirmation of these findings was provided through cell trajectory analysis, which positioned Cluster 0 at the initiation stage of development and Cluster 1 at the final developmental stage. The abundance of Cluster 0-2 groups in TCGA-LUAD samples was assessed using ssGSEA and subsequently categorized into high and low-expression groups. Notably, it was observed that Cluster 0-1 had a significant impact on survival (p<0.05). Furthermore, GSVA enrichment analysis demonstrated heightened activity in hallmark pathways for Cluster 0, whereas Cluster 1 exhibited notable enrichment in pathways related to cell proliferation. It is noteworthy that a prognostic model was established utilizing feature genes from Cluster 0-1, employing the Lasso and stepwise regression methods. The results revealed that in TCGA and GSE53624 cohorts, the low-risk group demonstrated significantly higher overall survival and increased levels of immune infiltration. An examination of four external immunotherapy cohorts unveiled that the low-risk group exhibited improved immunotherapeutic efficacy. Additionally, more meaningful treatment options were identified for the low-risk group.

Conclusion

The findings revealed distinct interactions between malignant epithelial cells of ESCC and subgroups within the tumor microenvironment. Two cell clusters, strongly linked to survival, were pinpointed, and a signature was formulated. This signature is expected to play a crucial role in identifying and advancing precision medicine approaches for the treatment of ESCC.

Recent advances on high-efficiency of microRNAs in different types of lung cancer: a comprehensive review

Carcinoma of the lung is among the most common types of cancer globally. Concerning its histology, it is categorized as a non-small cell carcinoma (NSCLC) and a small cell cancer (SCLC) subtype. MicroRNAs (miRNAs) are a member of non-coding RNA whose nucleotides range from 19 to 25. They are known to be critical regulators of cancer via epigenetic control of oncogenes expression and by regulating tumor suppressor genes. miRNAs have an essential function in a tumorous microenvironment via modulating cancer cell growth, metastasis, angiogenesis, metabolism, and apoptosis. Moreover, a wide range of information produced via several investigations indicates their tumor-suppressing, oncogenic, diagnostic assessment, and predictive marker functions in different types of lung malignancy. miRNA mimics or anti-miRNAs can be transferred into a lung cancer cell, with possible curative implications. As a result, miRNAs hold promise as targets for lung cancer treatment and detection. In this study, we investigate the different functions of various miRNAs in different types of lung malignancy, which have been achieved in recent years that show the lung cancer-associated regulation of miRNAs expression, concerning their function in lung cancer beginning, development, and resistance to chemotherapy, also the probability to utilize miRNAs as predictive biomarkers for therapy reaction.

Screening microRNAs as potential prognostic biomarkers for lung adenocarcinoma

OBJECTIVE

To screen and identify microRNAs (miRNAs) associated with the prognosis of lung adenocarcinoma (LUAD) using clinical samples and construct a prediction model for the prognosis of LUAD.

METHODS

160 patient samples were used to screen and identify miRNAs associated with the prognosis of LUAD. Differentially expressed miRNAs were analyzed using gene chip technology. The selected miRNAs were validated using samples from the validation sample group. Cox proportional hazards regression was used to construct the model and Kaplan-Meier was used to plot survival curves. Model power was assessed by testing the prognosis of the constructed model using real-time polymerase chain reaction (RT-PCR) data.

RESULTS

The data showed that miR-1260b, miR-21-3p and miR-92a-3p were highly expressed in the early recurrence and metastasis group, while miR-2467-3p, miR-4659a-3p, miR-4514, miR-1471 and miR-3621 were lowly expressed. It was further confirmed that miR-21-3p was significantly highly expressed in the early recurrence and metastasis group (p = 0.02). Receiver operating characteristic (ROC) curve results showed cut-off point value of 0.0172, sensitivity of 88.2% and specificity of 100%. The predictive results of the constructed model were in good agreement with the actual prognosis of patients by using the validation sample test (Kappa = 0.426, p < 0.001), with a model sensitivity of 74.4%, a specificity of 68.3%, and an accuracy of 71.3%.

CONCLUSION

miRNAs associated with the prognosis of patients with stage I LUAD were screened and validated, and a risk model for predicting the prognosis of patients was constructed. This model has good consistency with the actual prognosis of patients.

Mitochondrial-related genes markers that predict survival in patients with head and neck squamous cell carcinoma affect immunomodulation through hypoxia, glycolysis, and angiogenesis pathways

Mitochondria play a crucial role in the occurrence and development of tumors. We used mitochondria-related genes for consistent clustering to identify three stable molecular subtypes of head and neck squamous cell carcinoma (HNSCC) with different prognoses, mutations, and immune characteristics. Significant differences were observed in clinical characteristics, immune microenvironment, immune cell infiltration, and immune cell scores. TP53 was the most significantly mutated; cell cycle-related pathways and tumorigenesis-related pathways were activated in different subtypes. Risk modeling was conducted using a multifactor stepwise regression method, and nine genes were identified as mitochondria-related genes affecting prognosis (DKK1, EFNB2, ITGA5, AREG, EPHX3, CHGB, P4HA1, CCND1, and JCHAIN). Risk score calculations revealed significant differences in prognosis, immune cell scores, immune cell infiltration, and responses to conventional chemotherapy drugs. Glycolysis, angiogenesis, hypoxia, and tumor-related pathways were positively correlated with the RiskScore. Clinical samples were subjected to qPCR to validate the results. In this work, we constructed a prognostic model based on the mitochondrial correlation score, which well reflects the risk and positive factors for the prognosis of patients with HNSCC. This model can be used to guide individualized adjuvant and immunotherapy in patients with HNSCC.

Differentially Expressed Genes, miRNAs and Network Models: A Strategy to Shed Light on Molecular Interactions Driving HNSCC Tumorigenesis

Head and neck squamous cell carcinoma (HNSCC) is among the most common cancer worldwide, accounting for hundreds thousands deaths annually. Unfortunately, most patients are diagnosed in an advanced stage and only a percentage respond favorably to therapies. To help fill this gap, we hereby propose a retrospective in silico study to shed light on gene-miRNA interactions driving the development of HNSCC. Moreover, to identify topological biomarkers as a source for designing new drugs. To achieve this, gene and miRNA profiles from patients and controls are holistically reevaluated using protein-protein interaction (PPI) and bipartite miRNA-target networks. Cytoskeletal remodeling, extracellular matrix (ECM), immune system, proteolysis, and energy metabolism have emerged as major functional modules involved in the pathogenesis of HNSCC. Of note, the landscape of our findings depicts a concerted molecular action in activating genes promoting cell cycle and proliferation, and inactivating those suppressive. In this scenario, genes, including VEGFA, EMP1, PPL, KRAS, MET, TP53, MMPs and HOXs, and miRNAs, including mir-6728 and mir-99a, emerge as key players in the molecular interactions driving HNSCC tumorigenesis. Despite the heterogeneity characterizing these HNSCC subtypes, and the limitations of a study pointing to relationships that could be context dependent, the overlap with previously published studies is encouraging. Hence, it supports further investigation for key molecules, both those already and not correlated to HNSCC.

RNA-Based Liquid Biopsy in Head and Neck Cancer

Head and neck cancer (HNC) is a prevalent and diverse group of malignancies with substantial morbidity and mortality rates. Early detection and monitoring of HNC are crucial for improving patient outcomes. Liquid biopsy, a non-invasive diagnostic approach, has emerged as a promising tool for cancer detection and monitoring. In this article, we review the application of RNA-based liquid biopsy in HNC. Various types of RNA, including messenger RNA (mRNA), microRNA (miRNA), long non-coding RNA (lncRNA), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), circular RNA (circRNA) and PIWI-interacting RNA (piRNA), are explored as potential biomarkers in HNC liquid-based diagnostics. The roles of RNAs in HNC diagnosis, metastasis, tumor resistance to radio and chemotherapy, and overall prognosis are discussed. RNA-based liquid biopsy holds great promise for the early detection, prognosis, and personalized treatment of HNC. Further research and validation are necessary to translate these findings into clinical practice and improve patient outcomes.

Six-MicroRNA Prognostic Signature in Patients With Locally Advanced Head and Neck Squamous Cell Carcinoma

PURPOSE

MicroRNAs (miRNAs) have been evaluated as biomarkers in cancers. Therefore, we aimed to identify a prognostic miRNA signature from The Cancer Genome Atlas (TCGA) database and validate it in the Ramathibodi (RA) locally advanced head and neck squamous cell carcinoma (LA-HNSCC) cohort.

METHODS

The correlation between candidate miRNAs and the survival of patients with LA-HNSCC in TCGA database was analyzed. A prognostic miRNA signature model was generated that classified patients into high-risk and low-risk groups. This candidate miRNA signature was further validated in the independent RA cohort using droplet-digital polymerase chain reaction.

RESULTS

In TCGA database, we compared the expression of 277 miRNAs between 519 head and neck squamous cell carcinoma tissues and 44 normal tissues. The expression of hsa-miR-10b, hsa-miR-148b, hsa-miR-99a, hsa-miR-127, hsa-miR-370, and hsa-miR-500a was independently associated with overall survival (OS). Thus, we established the miRNA signature risk score from these six miRNAs and categorized patients into low-risk and high-risk groups. The median OS of TCGA patients was significantly shorter in the low-risk group than in the high-risk group (P < .001). The six-miRNA signature was further validated in the RA cohort (N = 87). The median OS of the low-risk group was significantly shorter compared with the high-risk group (P = .03). In multivariate analysis, the six-miRNA signature was an independent prognostic factor for OS in the RA cohort (HR, 1.958; 95% CI, 1.006 to 3.812; P = .048).

CONCLUSION

We identified a prognostic six-miRNA signature for patients with LA-HNSCC from TCGA cohort and validated it in our independent cohort. However, larger studies are warranted to confirm these results.

Application of miRNA Biomarkers in Predicting Overall Survival Outcomes for Lung Adenocarcinoma

Background

With the development of research, the importance of microRNAs (miRNAs) in the occurrence, metastasis, and prognosis of lung adenocarcinoma (LUAD) has attracted extensive attention. This study is aimed at predicting overall survival (OS) results through bioinformatics to identify novel miRNA biomarkers and hub genes.

Materials and Methods

The data of LUAD-related miRNA and mRNA samples was downloaded from The Cancer Genome Atlas (TCGA) database. Upon screening and pretreatment of initial data, TCGA data were analyzed using R platform and a series of analytical tools to identify biomarkers with high specificity and sensitivity.

Results

7 miRNAs and 13 hub genes that had strong relation to the overall surviving status were identified in patients with LUAD. The expression of seven miRNAs (hsa-miR-19a-3p, hsa-miR-126-5p, hsa-miR-556-3p, hsa-miR-671-5p, hsa-miR-937-3p, hsa-miR-4664-3p, and hsa-miR-4746-5p) could apparently improve the OS rate of patient with LUAD. The 13 hub genes, namely, CCT6A, CDK5R1, CEP55, DNAJB4, EGLN3, HDGF, HOXC8, LIMD1, MKI67, PCP4L1, PPIL1, SCAI, and STK32A, showed a correlation with the OS status.

Conclusion

7 miRNAs were identified as novel biomarkers for the prognosis of patients with LUAD. This study offered a deeper comprehension of LUAD treatment and prognosis from the molecular level and helped enhance the understanding of the pathogenesis and potential molecular events of LUAD.

Transcriptome analysis reveals the prognostic and immune infiltration characteristics of glycolysis and hypoxia in head and neck squamous cell carcinoma

BACKGROUND

This study aims to construct a new prognostic gene signature in survival prediction and risk stratification for patients with Head and neck squamous cell carcinoma (HNSCC).

METHOD

The transcriptome profiling data and hallmark gene sets in the Molecular Signatures Database was used to explore the cancer hallmarks most relevant to the prognosis of HNSCC patients. Differential gene expression analysis, weighted gene co-expression network analysis, univariate COX regression analysis, random forest algorithm and multiple combinatorial screening were used to construct the prognostic gene signature. The predictive ability of gene signature was verified in the TCGA HNSCC cohort as the training set and the GEO HNSCC cohorts (GSE41613 and GSE42743) as the validation sets, respectively. Moreover, the correlations between risk scores and immune infiltration patterns, as well as risk scores and genomic changes were explored.

RESULTS

A total of 3391 differentially expressed genes in HNSCC were screened. Glycolysis and hypoxia were screened as the main risk factors for OS in HNSCC. Using univariate Cox analysis, 97 prognostic candidates were identified (P < 0.05). Top 10 important genes were then screened out by random forest. Using multiple combinatorial screening, a combination with less genes and more significant P value was used to construct the prognostic gene signature (RNF144A, STC1, P4HA1, FMNL3, ANO1, BASP1, MME, PLEKHG2 and DKK1). Kaplan-Meier analysis showed that patients with higher risk scores had worse overall survival (p < 0.001). The ROC curve showed that the risk score had a good predictive efficiency (AUC > 0.66). Subsequently, the predictive ability of the risk score was verified in the validation sets. Moreover, the two-factor survival analysis combining the cancer hallmarks and risk scores suggested that HNSCC patients with the high hypoxia or glycolysis & high risk-score showed the worst prognosis. Besides, a nomogram based on the nine-gene signature was established for clinical practice. Furthermore, the risk score was significantly related to tumor immune infiltration profiles and genome changes.

CONCLUSION

This nine-gene signature associated with glycolysis and hypoxia can not only be used for prognosis prediction and risk stratification, but also may be a potential therapeutic target for patients with HNSCC.

A Pan-Cancer Analysis of the Oncogenic Role of Cell Division Cycle-Associated Protein 4 (CDCA4) in Human Tumors

Purpose

To unravel the oncogenic role of CDCA4 in different cancers from the perspective of tumor immunity.

Methods

Raw data on CDCA4 expression in tumor samples and paracancerous samples were obtained from TCGA and GTEX databases. In addition, we investigated pathological stages and the survival analysis of CDCA4 in pan-cancer across Gene Expression Profiling Interactive Analysis (GEPIA) database. Cox Proportional Hazards Model shows that high CDCA4 levels are associated with several vital indicators in oncology. On the one hand, we explored the correlation between CADA4 expression and tumor immune infiltration by the TIMER tool; On the other hand, we utilized the methods of CIBERSORT and ESTIMATE computational to evaluate the proportion of tumor infiltrating immune cells (TIIC) and the amounts of stromal and immune components based on TCGA database. The use of antineoplastic drugs and the expression of CDCA4 also showed a high correlation via linear regression. Protein–Protein Interaction analysis was performed in the GeneMANIA database, and enrichment analysis was performed and predicted signaling pathways were identified by using Gene Ontology and Kyoto Encyclopedia of Genes. The correlation between CDCA4 expression with Copy number variations (CNV) and methylation is detailed, respectively. Molecular biology experiments including Western blotting, flow cytometry, EDU staining, Transwell and Wound Healing assay to validate the cancer promoting role of CDCA4 in hepatocellular carcinoma (HCC).

Results

Most tumors highly expressed CDCA4. Elevated CDCA4 expression was associated with poor OS and DFS. There was a significant correlation between CDCA4 expression and TITCs. Moreover, markers of TIICs exhibited distinct patterns of CDCA4 associated immune infiltration. In addition, we pay attention to the association between the expression of CDCA4 and the use of the anti-tumor drugs. CDCA4 is related to biological progress (BP), cellular component (CC) and molecular function (MF). Dopaminergic Synapse, AMPK, Sphingolipid, Chagas Disease, mRNA Surveillance were significantly enriched pathways in positive and negative correlation genes with CDCA4. CNV is thought to be a positive correlation with CDCA4 expression. Conversely, methylation is negative correlation with CDCA4 expression. Molecular biology experiments confirm a cancer promoting role for CDCA4 in HCC

Conclusion

CDCA4 may serve as a biomarker for cancer immunologic infiltration and poor prognosis, providing a new way of thinking for cancer treatment.

Development of a novel necroptosis-associated miRNA risk signature to evaluate the prognosis of colon cancer patients

BACKGROUND

Necroptosis is a recently discovered caspase-independent form of cell death which plays an important role in the occurrence and development of cancer. As an important regulatory factor in necroptosis, microRNAs (miRNAs) are important for the development of colon cancer. This study established a novel necroptosis-related miRNA risk signature to evaluate the prognosis of patients with colon adenocarcinoma (COAD).

METHODS

The necroptosis-related miRNAs were selected by assessing the differential expression of miRNAs in 459 COAD patient samples and 8 control samples from The Cancer Genome Atlas (TCGA). Selection operator Cox analyses and survival analyses were used to establish the risk signature of 7 miRNAs related to necroptosis. Functional enrichment analysis and nomograms were used to explore the potential effects of necroptosis-related miRNAs on prognosis and metastasis. The target genes of the necroptosis-related miRNAs were predicted using online databases and the genes related to overall survival (OS) were screened.

RESULTS

The risk signature was based on 7 necroptosis-related miRNAs. Nomograms showed that the risk signature was effective at predicting the prognosis and TNM stage of COAD patients. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses demonstrated that these miRNAs play an important role in cancer development, metastasis, and prognosis. A total of 38 target genes for these miRNAs were found to be associated with the OS in COAD patients.

CONCLUSIONS

This study provided novel evidence that necroptosis-related miRNAs are associated with the prognosis of COAD patients. A risk signature established based on these miRNAs could effectively predict the prognosis and metastasis of COAD in patients.

Multiple microRNA signature panel as promising potential for diagnosis and prognosis of head and neck cancer

MicroRNAs are small non-coding RNA that regulate gene expressions of human body. To date, numerous studies have reported that microRNAs possess great diagnostic and prognostic power in head and neck cancer and had governed a lot of attention. The factor for the successfulness of miRNAs in these aspects is due to cancer being fundamentally tied to genetic changes, which are regulated by these miRNAs. Head and neck cancer, leading the world record for cancer as number sixth, is caused by multiple risk factors such as tobacco consumption, alcohol consumption, dietary factors, ethnicity, family history, and human papilloma virus. It derives at locations such as oral cavity, pharynx, larynx, paranasal sinus and salivary gland and have high rate of mortality with high recurrence rate. Besides, head and neck cancer is also usually having poor prognosis due to its asymptomatic nature. However, this diagnostic and prognostic power can be further improved by using multiple panels of miRNA as a signature or even combined with TNM staging system to obtain even more remarkable results. This is due to multiple factors such as tumour heterogeneity and components of the tumour which may affect the composition of miRNAs. This review covers the examples of such miRNA signatures, compare their diagnostic and prognostic powers, discuss some controversial roles of unreported miRNAs, and the molecular mechanisms of the miRNAs in gene targeting and pathways.

Aberrant ROS Mediate Cell Cycle and Motility in Colorectal Cancer Cells Through an Oncogenic CXCL14 Signaling Pathway

Background: Reactive oxygen species (ROS) act as signal mediators to induce tumorigenesis.

Objective: This study aims to explore whether chemokine CXCL14 is involved in the proliferation and migration of ROS-induced colorectal cancer (CRC) cells.

Methods: The proliferative and migratory capacities of CRC cells treated with or without H₂O₂ were measured by various methods, including the CKK-8 assay, colony formation assay, flow cytometry, wounding healing assay, and migration assay.

Results: The results revealed that H₂O₂ promoted the proliferation and migration of CRC cells by regulating the cell cycle progression and the epithelial to mesenchymal transition (EMT) process. Furthermore, we noted that the expression level of CXCL14 was elevated in both HCT116 cells and SW620 cells treated with H₂O₂. An antioxidant N-Acetyl-l-cysteine (NAC) pretreatment could partially suppress the CXCL14 expression in CRC cells treated with H₂O₂. Next, we constructed CRC cell lines stably expressing CXCL14 (HCT116/CXCL14 and SW620/CXCL14) and CRC cell lines with empty plasmid vectors (HCT116/Control and SW620/Control) separately. We noted that both H₂O₂ treatment and CXCL14 over-expression could up-regulate the expression levels of cell cycle-related and EMT-related proteins. Moreover, the level of phosphorylated ERK (p-ERK) was markedly higher in HCT116/CXCL14 cells when compared with that in HCT116/Control cells. CXCL14-deficiency significantly inhibited the phosphorylation of ERK compared with control (i.e., scrambled shNCs). H₂O₂ treatment could partially restore the expression levels of CXCL14 and p-ERK in HCT116/shCXCL14 cells.

Conclusion: Our studies thus suggest that aberrant ROS may promote colorectal cancer cell proliferation and migration through an oncogenic CXCL14 signaling pathway.

MicroRNA-375: potential cancer suppressor and therapeutic drug

MiR-375 is a conserved noncoding RNA that is known to be involved in tumor cell proliferation, migration, and drug resistance. Previous studies have shown that miR-375 affects the epithelial-mesenchymal transition (EMT) of human tumor cells via some key transcription factors, such as Yes-associated protein 1 (YAP1), Specificity protein 1 (SP1) and signaling pathways (Wnt signaling pathway, nuclear factor κB (NF-κB) pathway and transforming growth factor β (TGF-β) signaling pathway) and is vital for the development of cancer. Additionally, recent studies have identified microRNA (miRNA) delivery system carriers for improved in vivo transportation of miR-375 to specific sites. Here, we discussed the role of miR-375 in different types of cancers, as well as molecular mechanisms, and analyzed the potential of miR-375 as a molecular biomarker and therapeutic target to improve the efficiency of clinical diagnosis of cancer.

Constructe a novel 5 hypoxia genes signature for cervical cancer

Background

Hypoxia, which affects the development, metastasis and prognosis of cancer, represents a key feature of cancer. This study describe a hypoxia risk factor model, with predicting the prognosis of cervical cancer.

Methods

Based on hypoxia pathway related genes, we divided cervical cancer samples into high and low expression groups. A cox analysis was then performed. Genes from these cervical cancer samples showing a significant impact on OS were selected for cluster analysis to obtain two subtypes. The TPM dataset of TCGA was divided into training and validation sets. For the training set, a lasso analysis was conducted as based on cox analysis of meaningful genes and a risk factor model was constructed. The constructed model was verified in internal and external data sets. Finally, RT-PCR, immunohistochemistry were used to detect the expression of relative genes or proteins and functional assays were used to evaluate the biological function of signature genes.

Results

Two molecular subtypes were obtained, Cluster2 vs Cluster1.These subtypes were obtained by clustering with a total of 149 DEGs (Differential expressed genes) being in line with this standard, of which 27 were up-regulated and 122 were down-regulated. The five genes with lambda = 0.0571 were selected to construct the model, the RiskScore = AK4*0.042 + HK2*0.021 + P4HA1*0.22 + TGFBI*0.1 + VEGFA*0.077. Further, in order to verify the signature, we used TCGA-test and GSE44001 chip datasets to test, and finally got a good risk prediction effect in those datasets. Moreover, the result of RT-PCR and immunohistochemistry demonstrated that AK4, HK2, P4HA1, TGFBI and VEGFA were all highly expressed in these cervical cancer tissue samples. The functional study shown that expression of AK4, HK2, P4HA1, TGFBI and VEGFA can regulate the proliferation, migration, and invasion ability of cervical cancer cells in vitro.

Conclusions

In summary, we developed a 5-gene signature prognostic hierarchical system based on the hypoxic pathway of cervical cancer, which is independent of clinical characteristics. And also conducted experimental verifications on these signature gene. Therefore, we propose that use of this classifier as a molecular diagnostic test can provide an effective means for evaluating the prognostic risk of cervical cancer patients, and provide potential targets for the treatment of cervical cancer patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02050-3.

A Super-Enhancer Driven by FOSL1 Controls miR-21-5p Expression in Head and Neck Squamous Cell Carcinoma

MiR-21-5p is one of the most common oncogenic miRNAs that is upregulated in many solid cancers by inhibiting its target genes at the posttranscriptional level. However, the upstream regulatory mechanisms of miR-21-5p are still not well documented in cancers. Here, we identify a super-enhancer associated with the MIR21 gene (MIR21-SE) by analyzing the MIR21 genomic regulatory landscape in head and neck squamous cell carcinoma (HNSCC). We show that the MIR21-SE regulates miR-21-5p expression in different HNSCC cell lines and disruption of MIR21-SE inhibits miR-21-5p expression. We also identified that a key transcription factor, FOSL1 directly controls miR-21-5p expression by interacting with the MIR21-SE in HNSCC. Moreover, functional studies indicate that restoration of miR-21-5p partially abrogates FOSL1 depletion-mediated inhibition of cell proliferation and invasion. Clinical studies confirmed that miR-21-5p expression is positively correlated with FOSL1 expression. These findings suggest that FOSL1-SE drives miR-21-5p expression to promote malignant progression of HNSCC