Open-access Comparative analysis of miRNA-mRNA interaction prediction tools based on experimental head and neck cancer data

ABSTRACT

Objective:  Head and neck squamous cell carcinoma (HNSCC) has a poor prognosis largely due to late diagnosis and a lack of reliable biomarkers. MicroRNAs (miRNAs), small non-coding RNAs that regulate gene expression, are promising biomarkers for HNSCC. This study evaluated miRNA-mRNA interactions in HNSCC using conventional computational tools and validated the results using molecular data.

Methods:  We compared three miRNA-mRNA interaction prediction tools, TargetScan, miRDB, and miRWalk, using differentially expressed miRNAs and mRNAs from HNSCC and cancer-free tissues. NanoString nCounter was used to measure miRNA and mRNA expression and the miRTarBase database was used to validate the predicted miRNA-mRNA interactions.

Results:  TargetScan and miRWalk provide a comprehensive overview of potential interactions, whereas miRDB provides functional insights. Our results identified 77 and 154 differentially expressed miRNAs and mRNAs in HNSCC, respectively. miRWalk predicted the highest number of miRNA-mRNA interactions, followed by miRDB and TargetScan. Only 3.3% of interactions were common among the tools. The MiRTarBase analysis confirmed a small subset of the predictions. Biological pathway analysis highlighted the dysregulation of PI3K-Akt and Wnt signaling; miRWalk was the best for elucidating how miRNAs modulate target mRNAs in these key pathways during HNSCC progression.

Conclusion:  miRWalk emerged as the most robust tool for predicting miRNA-mRNA interactions. Our findings highlight the importance of integrating bioinformatics predictions with experimental data to better understand the regulatory networks in HNSCC and identify potential biomarkers for diagnosis and therapy.

Keywords:
Squamous cell carcinoma of head and neck; MicroRNAs; Biomarkers; Gene expression; Target prediction tool; Computational biology; Algorithms

Highlights

  • miRWalk had the highest predicted interactions and validated miRNA networks in HNSCC.

  • Around 3.3% of interactions overlapped across tools, emphasizing the need for multitool approaches.

  • Dysregulated genes and miRNAs were tied to cancerdriving PI3K-Akt and Wnt pathways.

  • The validated approach highlights the importance of integrating computational and molecular data.


In Brief

We evaluated the performance of TargetScan, miRDB, and miRWalk for predicting miRNA-mRNA interactions in HNSCC. Based on clinical tumor and cancer-free tissue data, miRWalk emerged as the most comprehensive tool. Validation using NanoString technology and MiRTarBase confirmed key predictions, highlighting the important roles of the PI3K-Akt and Wnt pathways. This study underscores the importance of integrating bioinformatics and experimental data to better understand HNSCC.

INTRODUCTION

Head and neck cancer is the seventh most common type of cancer worldwide.(1) In Brazil, >10, 000 deaths were caused by head and neck cancer in 2020, and 22,890 new cases have been predicted for 2023-2025. The disease has a higher prevalence among men, among whom it was anticipated to rank among the top 10 cancers in 2023.(2) Head and neck squamous cell carcinoma (HNSCC) accounts for ~90% of this tumor type.(3) Head and neck squamous cell carcinoma comprises a heterogeneous group of tumors originating in the oral cavity, larynx, nasopharynx, oropharynx, and hypopharynx. It has an estimated 5-year survival rate of 40-50%,(4) primarily due to delays in diagnosis and treatment initiation, which is directly impacted by the lack of specific markers.(5)

MicroRNAs (miRNAs) are small non-coding RNA molecules that exert post-transcriptional regulatory functions in conjunction with the RNA-induced silencing complex (RISC). They target the 3’ untranslated region (UTR), 5’ UTR, or promoter region of mRNAs through base pair complementarity.(6) In addition to their biological role, miRNAs are candidate biomarkers owing to their dysregulated expression in tumor tissues; their small size makes them less susceptible to degradation in formalin-fixed paraffin-embedded (FFPE) tissues-the primary sources of clinical specimens.(7) Several studies have suggested that differentially expressed miRNAs act as oncomiRs, tumor suppressors, and prognostic markers in HNSCC.(8-12)

Owing to the heterogeneity in HNSCC tumors, there is currently no consensus on effective biomarkers. Further studies are needed to identify diagnostic and prognostic markers and elucidate tumor development, progression, and dynamics.

Identifying the targets of miRNA-mRNA interactions can aid in the discovery of miRNA regulatory networks. In this context, the most effective approach involves computational prediction followed by the validation of miRNA-mRNA interactions. Bioinformatic analyses have shown that a single miRNA can regulate the expression of several genes, and a gene can be controlled by multiple miRNAs. Therefore, the experimental validation of each potential miRNA target is impractical, costly, and time-consuming. Computational approaches for predicting miRNA targets simplify the process by expediting the selection and reducing the number of potential targets for validation.

Several miRNA-mRNA interaction prediction tools are available.(13) As described by Nazarov and colleagues,(14) some use data-driven methods and other target-based methods. Data-driven methods based on similarity use large transcriptomic datasets and advanced methods like correlation analysis and biclustering to accurately identify and prioritize miRNA-mRNA interactions, accounting for complexities and context-specific relationships and improving the accuracy of interaction predictions. Data-driven methods based on matrix factorization integrate mRNA and miRNA expression data by decomposing the expression matrix into lower-rank matrices to capture data variability and identify shared latent variables. Target-based methods can be split into two categories, the first of which is target prediction, such as TargetScan,(15) DIANA-microT-CDS,(16) miRDB,(17) STarMir,(18) miRGator,(19) miRGate,(20) and DeepMirTar.(21) By combining different features and correlations, these tools employ distinct algorithms and methodologies to predict miRNA targets, considering various features such as sequence complementarity, binding site conservation, thermodynamic stability, and experimental validation. Second, various online repositories have been developed for collecting data on experimentally validated interactions between mRNA and miRNAs, such as Diana-Tarbase,(22) miRTarBase,(23) miRecords,(24) and miR2Disease.(14,25,26)

In this study, we selected the most commonly used prediction tools to analyze potential miRNA targets in HNSCC. We validated the in silico results using molecular barcode hybridization technology (NanoString nCounter) for mRNAs and miRNAs in paired cancer and non-cancer FFPE tissues. This technology is particularly well-suited for situations involving small fragments and a high rate of RNA degradation, such as in FFPE samples. In addition, this platform does not require high computational performance or expertise in data analysis, making it ideal for clinical laboratory applications.

OBJECTIVE

We aimed to compare the performance of miRNA-mRNA interaction prediction tools, considering both experimental and literature evidence, in accurately representing the characteristics of head and neck squamous cell carcinoma tissues.

METHODS

Literature search for miRNA-mRNA interaction prediction tools and comparison using experimental results

We used PubMed Central to search for articles published in English from 2018 to 2024 containing the terms "MicroRNAs," "mRNA-Gene," "Expression Regulation," "carcinoma," "interaction," and "prediction." Tools were chosen based on their online availability and use of an original algorithm or curation method for identifying mRNA-miRNA interactions.

For comparison among tools, we used data on the differentially expressed miRNAs detected between HNSCC and cancer-free tissues, as input. The higher-scoring gene targets for each tool, limited to 10% of the target gene set list (higher scores), were selected and compared with the differentially expressed mRNAs between HNSCC samples and cancer-free samples. Additionally, we also used the MiRTarBase database to search for experimentally validated miRNA-mRNA interactions.

Study population and sample collection

The cohort comprised 18 patients with a confirmed anatomopathological diagnosis of oropharyngeal squamous cell carcinoma and 9 patients with laryngeal squamous cell carcinoma (Table 1). The average patient age was 58 years (range, 39-85 years) and the male/female ratio was approximately 4:1. Most of the patients were smokers or former smokers (Table 1). All samples were obtained from patients treated at the Instituto de Cancer, Dr. Arnaldo Vieira de Carvalho (ICAVC, São Paulo, Brazil). Clinical TNM classification was performed according to the Union for International Cancer Control staging classification system for HNSCC.(27) This study was approved by the Institutional Review Board of Hospital Israelita Albert Einstein (CAAE: 94534318.3.0000.0071; # 3.780.449), and written informed consent was obtained from all participants. Analysis of hematoxylin and eosin (H&E)-stained sections was performed by expert pathologists, who confirmed the presence of ≥50% tumor cells in all HNSCC samples.

Table 1
Clinical and pathological characteristics of the study group

RNA isolation

We used 21 tissue samples from patients with HNSCC (12 oropharynx and 9 larynx) for miRNA profiling, and 12 HNSCC samples were used for mRNA expression analysis. Formalin-fixed paraffin-embedded tissue blocks were cut into 10μm-thick sections using a microtome and adhered to glass slides. Based on the H&E-stained slides, tumor and non-tumor areas were separately scraped using a scalpel and loaded into a 1.5mL tube. Total RNA, including miRNA, was isolated from FFPE tissues using an AllPrep DNA/RNA FFPE kit (Qiagen, Hilden, Germany) according to the manufacturer's protocol.

miRNA expression

An nCounter Human v3 miRNA Expression Assay (NanoString Technologies, Seattle, WA, USA) was used for the miRNA profiling of 21 tissue samples from patients with HNSCC (n=12 oropharynx and 9 larynx) using the total RNA isolated from HNSCC and cancer-free FFPE tissues (60-100ng), according to the manufacturer's recommendations. Each cartridge was scanned with a 555 FOV using an nCounter Digital Analyzer. Data processing was performed using nSolver Analysis Software 4.0 (NanoString Technologies). Background subtraction was performed using the geometric mean of the negative controls plus one standard deviation, normalization to the geometric mean of the positive controls, and normalization to the geometric mean of the top 100 most highly expressed miRNAs. After normalization, we filtered miRNAs with fewer than three counts from >30% of the samples.

mRNA expression

For mRNA, tumor and non-tumor FFPE tissue RNA from 12 patients with oropharyngeal squamous cell carcinoma were analyzed using the nCounter PanCancer Pathways Panel (NanoString Technologies). This tool offers fast and high-throughput multiplex capabilities for analyzing human gene expression. Each panel comes with 770 genes from 13 cancer-associated canonical pathways representing the basic cancer biology [Notch, MAPK, STAT, PI3K, RAS, cell cycle, apoptosis, Hedgehog, APC (Wnt), DNA damage and repair control, transcriptional misregulation, chromatin modification, and TGF-β]. We used ~300ng RNA in the assay, according to the manufacturer's recommendations. Samples with <100 counts for 100% of the target mRNAs were excluded from further analysis. Data were normalized using the geNorm algorithm, which selects the most stable housekeeping genes that exhibit the least variation among the samples.

Statistical analysis

A two-tailed t-test was used to assess the miRNA counts obtained using the nCounter Human v3 assay and mRNA counts derived using the nCounter PanCancer Pathways Panel. nSolver Analysis Software 4.0 was used for statistical analysis, with significance set at p≤0.05 for miRNAs and p<0.05 for mRNAs. Adjusted p-values corrected for multiple testing using the Benjamini-Yekutieli method (assigned as FDR) were estimated for genes involved in key biological processes.

RESULTS

Most commonly used miRNA-mRNA interaction tools

Based on the literature search and accounting for the limitations of the keywords, a total of 247 articles were identified, of which 150 matched our inclusion criteria (Figure 1).

Figure 1
Most used miRNA-mRNA interaction tools from PubMed search for 2018-2024. The x-axis lists the tools that met the study inclusion criteria and y-axis shows the number of publications utilizing each tool

Of the 150 studies, 67% used TargetScan 8.0, 46% miRDB, and 36% miRWalk 2.0. Thus, we selected these tools for our prediction analysis and used MiRTarBase to validate the predicted interactions.

TargetScan 8.0(17) is the most widely used miRNA-mRNA interaction prediction tool, with the current version available since 2021. This tool predicts miRNA targets based on seed region complementarity and evolutionary conservation by applying a context score that reflects the potential of a miRNA-binding site for creating an effective miRNA-mRNA interaction, considering 14 parameters accounting for interactions in canonical and non-canonical sites.

The miRWalk 2.0 database(28) contains results from the TarPmiR algorithm,(29) a random forest-based approach that predicts miRNA-mRNA interactions through the use of "decision trees" with set rules for decision making. It is similar to a fluxogram with "knots" where, if the condition is verified, flow is continued in the same direction, otherwise a different direction is taken towards the next knot. This approach integrates key features such as accessibility, seed match, flanking conservation, folding energy, the length of the largest consecutive pairings, the length of the target site, AU content, stem conservation, difference between stem and flanking conservation, m/e motif, and difference between the numbers of paired positions in the seed region and miRNA 3’ end to predict miRNA target sites.(29)

miRDB(17) is another online database that uses the machine learning algorithm MirTarget(30) for target prediction. Together with the more commonly used miRNA-targeting features used for predicting interactions, this tool combines data from miRNA overexpression and crosslinking and immunoprecipitation binding experiments data to identify transcript targets associated with functional RISC complexes, providing a functional weight for the predictions.

The miRTarBase database(22) relies on experimentally validated interactions, curated data from the literature, and high-throughput experiments. miRTarBase includes interactions confirmed through methods such as reporter assays (western blotting, microarray experiments, and high-throughput sequencing). The database is compiled of data from manually reviewed scientific articles reporting validated miRNA-target interactions.

miRNAs and mRNAs expressed in HNSCC

We identified 77 miRNAs that were differentially expressed between cancer and cancer-free areas; 38 and miRNAs were upregulated and downregulated in tumor tissues, respectively (Figure 2). For mRNAs, 154 were differentially expressed between cancer and cancer-free tissues; 87 were upregulated and 67 were downregulated in tumor tissues. Among the differentially expressed genes, 15 were previously reported to be dysregulated in HNSCC (Table 2). These genes participate in the cell cycle, apoptosis, and PI3K, Wnt, and Hedgehog pathways (Table 2).

Figure 2
MicroRNAs differentially expressed between tumor and non-tumor tissues (p≤0.05). Positive and negative fold-change describe upregulated and downregulated miRNAs in tumor tissue, respectively
Table 2
Genes differentially expressed (p<0.05) in tumor tissue and previously described for HNSCC

miRNA-mRNA interaction predictions and biological implications

We evaluated the potential mRNA targets of the differentially expressed miRNAs using the different prediction tools. Using the 77 differentially expressed miRNAs as input, we identified interactions with 154 differentially expressed mRNAs in tumor tissues: 2,124 interactions using miRWalk, 187 using miRDB, and 82 using TargetScan (Figure 3), of which 42, 9, and 17 were validated using MiRTarBase; however, none were previously associated with HNSCC. Notably, one miRNA can regulate multiple target genes, whereas a single gene can be regulated by several miRNAs.(31)

Figure 3
Differentially expressed genes targeted by differentially expressed miRNAs between tumor and tumor-free tissues. MiRWalk, miRDB, and TargetScan databases were used to predict miRNA targets. The numbers of miRNA-mRNA interactions are presented in the Venn diagram

Only 3.3% (76/2,314) of the interactions were shared among at least two of the prediction tools. Only 2.4% (59/2,463) of the predicted interactions were validated using MiRTarBase. Although the three tools predicted three common interactions (Figure 3), these were not consistent with the experimentally validated database.

Based on our data and the literature, the altered expression of mRNAs in this cohort promotes tumor progression by enhancing proliferative signals and PI3K-Akt and WNT pathway activities. Figure 4 illustrates these alterations based on differentially expressed miRNAs and target genes predicted using at least one of the four tools. In tumor tissues, we identified the upregulation of genes involved in activating the cell cycle and DNA replication, and downregulation of miRNAs predicted to target these genes. miRWalk predicted the most interactions between miRNAs and mRNAs, followed by miRDB and TargetScan, with 111, seven, and five interactions, respectively (Figure 4).

Figure 4
Proposed model of how differently expressed mRNA and miRNA may contribute to tumor progression. The histogram represents the number of interactions predicted by each of the tools for differentially expressed mRNA and miRNA. miRNAs in red and green were downregulated and upregulated in tumor tissue, respectively. Green and red arrows represent upregulated and downregulated genes in tumor tissue, respectively. Genes in blue were differentially expressed (p<0.05). Genes in yellow had an adjusted p<0.05 (FDR)

DISCUSSION

Several miRNA-mRNA interaction prediction tools have been developed, each applying distinct algorithms and methodologies to predict miRNA targets. These algorithms consider various features, such as sequence complementarity, binding site conservation, thermodynamic stability, and experimental validation, to promote the accuracy of these tools.

Although TargetScan, miRDB, and miRWalk share similar features(32) (such as sequence complementarity, total number of paired positions, target site accessibility, and conservation), only ~3% of the miRNA-mRNA interactions was predicted by more than one of them. This is consistent with the 13% similarity previously reported between miRDB and TargetScan.(14) In another study, in which 10 of the top miRNAs with MiRTarBase-validated targets were analyzed using TargetScan, miRanda-mirSVR, Pita, and RNA22, 18% of the validated interactions were not predicted by any tool.(33) This suggests that biological processes may not be fully captured by prediction tools.

The best prediction tool in our study was miRWalk, since it had the highest numbers of predicted and validated interactions. Similarly, a previous study reported better performance for miRWalk among several prediction tools.(29) miRWalk also shared the most interactions with the other tools.

Despite the low similarity between the different tools, they are the most commonly used among researchers in this field. This could be because of their continued updates and community involvement, whereas many other tools have been discontinued.(34) In terms of publication count, TargetScan was the most common, possibly because of the amount of information it provides and its pioneering role in predicting miRNA targets in vertebrate research.(35) In our study, TargetScan identified the fewest validated interactions, likely because of its high false negative rate and stringent criteria for interaction sites.(36) These criteria require high confidence, canonical interactions primarily within the 3’-UTR, and evolutionary conservation. In contrast, miRWalk is more flexible and inclusive, broadening the range of potential interactions to those with mismatches or noncanonical binding patterns across the entire mRNA. Notably, our analysis was limited to the top-scoring interactions from all four tools, which may have contributed to the lower perceived sensitivity of TargetScan.

The nCounter PanCancer Pathways Panel allowed us to identify differences in the expression of genes related to commonly dysregulated pathways in cancer. Genes involved in the PI3K pathway, which are crucial for cellular functions such as transcription, translation, proliferation, and survival, are often overexpressed in HNSCC.(37) We identified several candidate miRNA-mRNA interactions that should be investigated as targets in HNSCC. Notably, integrin subunit alpha 3 (ITGA3) and ITGA6 showed significantly increased expression in patients with HNSCC, with higher levels of ITGA3 linked to poorer overall and relapse-free survival rates.(38) We observed the downregulation of miR-150-5p and upregulation of ITGA6 in tumor tissues, consistent with the findings of Koshizuka et al.(38) This interaction was predicted only by miRWalk. Although it uses a combination of manual curation and machine learning methods to validate interactions with high accuracy (>82%), MiRTarBase did not capture this interaction, indicating possible limitations in curation.

ITGB4 is upregulated in primary and metastatic HNSCC, associated with poor prognosis, and its knockdown reduces EGFR-mediated migration and invasion.(39,40) Overexpression of ITGA6 and laminin subunit gamma 2 (LAMC2) was previously reported in HNSCC; LAMC2 has been associated with poor prognosis and identified as a potential biomarker for cancer invasion.(41-44) Laminin subunit beta 3 (LAMB3) enhances cisplatin sensitivity and reduces cell migration and invasion.(45,46) In our study, these genes were upregulated in tumor tissues and miRNAs predicted to interact with these genes were identified as less expressed in the tumor tissue. Collectively, these findings support the role of these genes in HNSCC tumorigenesis and progression.

Overactivation of WNT-β-catenin signaling promotes cancer cell proliferation and survival.(47) Consistent with our findings, Wnt family member 7 B (WNT7B) is upregulated in oral squamous cell carcinoma (OSCC), and in vitro knockdown of WNT7B in OSCC cell lines inhibits cell invasion.(48) Overexpression of WNT4 in laryngeal cancer cell lines promoted the nuclear accumulation of β-catenin while knockdown of WNT4 had the opposite results.(49) In our study, these genes were upregulated in tumor tissues corroborating the literature, while miRNAs predicted to have these genes as targets were found with less expression in tumor tissue, such as miR-30a and miR-150.

Genes that regulate the cell cycle and apoptosis are crucial for tumor progression and cell proliferation. In our tumor samples, we observed the overexpression of S-phase kinase-associated protein 2 (SKP2) and downregulation of cyclin-dependent kinase inhibitor 1C (CDKN1C). Notably, SKP2 expression has been correlated with shorter overall survival in patients with HNSCC(50) and is often inversely associated with CDKN1C expression in various cancers.(51,52) In OSCC, positive CDKN1C staining was correlated with higher 5-year survival rates. Additionally, CDKN1C mRNA expression decreased from leukoplakia with moderate or severe dysplasia to OSCC.(53) Consistent with other studies, we observed the upregulation of protein kinase membrane-associated tyrosine/threonine 1 (PKMYT1), cyclin A2 (CCNA2), minichromosome maintenance complex component 7 (MCM7), cell division cycle 6 (CDC6),(54-56) cyclin-dependent kinase 4 (CDK4), and CDK6, which are potential immunotherapy targets in HNSCC.(57,58) Additionally, we observed reduced expression of miRNAs (such as miR-145, miR-30a-3p, and let-7c-5p) predicted by miRWalk to target CDK6 and we found these interactions validated in miRTarBase database.(59-61) Further research into these regulatory networks may unveil new therapeutic strategies and biomarkers for the better management and prognosis of this challenging disease.

Changes in the cellular phenotype of tumor cells reflect changes in the expression of several genes caused by different carcinogenic and/or biological processes. Here, we hypothesized that differences in mRNA expression between tumor and non-tumor tissues may be partially caused by differences in miRNA expression, as supported by the miRNA-mRNA interaction tools. Our study provides useful information for elucidating the mechanisms of tumor progression in HNSCC, especially given the current lack of biomarkers.

In summary, depending on the research aim, our findings demonstrated the effective use of TargetScan and miRWalk for gaining a broad overview of potential interactions, miRDB for gaining functional insights, and miRTarBase for validation. Although the focus of miRTarBase on experimentally validated interactions may enhance reliability, it may limit the scope of the results. Combining these four tools may address some of their separate limitations while complementing their strengths. Although these methods are suitable for preliminary investigations, the predictions should be supplemented with experimentally validated data to ensure biological relevance. Through experimental validation, we identified miRWalk as the most appropriate tool for assessing how miRNAs modulate differentially expressed target mRNAs in HNSCC, followed by miRDB and TargetScan. Considering the limitations of each database, miRWalk enables comprehensive exploration by integrating data from the other three databases, offering a well-rounded approach to studying miRNA interactions.

CONCLUSION

Databases and online tools provide excellent resources for the initial review and identification of miRNA and mRNA interactions. However, their output must be interpreted cautiously, considering the limitations of each method.

ACKNOWLEDGEMENTS

The authors acknowledge the contributions of the Instituto de Cancer, Dr. Arnaldo Vieira de Carvalho clinical research team and clinical laboratory during participant recruitment. The authors sincerely thank Dr. Nelson Hamerschlak for his invaluable support and efforts in securing funding for this study. This work was supported by Programa Nacional de Apoio a Atençao Oncologica (PRONON). Larissa Figueiredo Alves Diniz and Bárbara dos Santos Dias are supported by a CAPES fellowship.

REFERENCES

  • 1 Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209-49.
  • 2 Instituto Nacional de Câncer (INCA). Estimativa 2023: Incidência de câncer no Brasil. Rio de Janeiro: INCA; 2022 [citado 2024 Nov 29]. Disponível em: https://www.inca.gov.br/sites/ufu.sti.inca.local/files/media/document/estimativa-2023.pdf
    » https://www.inca.gov.br/sites/ufu.sti.inca.local/files/media/document/estimativa-2023.pdf
  • 3 Curado MP, Hashibe M. Recent changes in the epidemiology of head and neck cancer. Curr Opin Oncol. 2009;21(3):194-200.
  • 4 Lu L, Wu Y, Feng M, Xue X, Fan Y. A novel seven-miRNA prognostic model to predict overall survival in head and neck squamous cell carcinoma patients. Mol Med Rep. 2019;20(5):4340-8.
  • 5 Liao DZ, Schlecht NF, Rosenblatt G, Kinkhabwala CM, Leonard JA, Ference RS, et al. Association of Delayed Time to Treatment Initiation With Overall Survival and Recurrence Among Patients With Head and Neck Squamous Cell Carcinoma in an Underserved Urban Population. JAMA Otolaryngol Head Neck Surg. 2019;145(11):1001-9.
  • 6 Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281-97.
  • 7 von Ahlfen S, Missel A, Bendrat K, Schlumpberger M. Determinants of RNA quality from FFPE samples. PLoS One. 2007;2(12):e1261.
  • 8 Li H, Liu P, Li D, Wang Z, Ding Z, Zhou M, et al. STAT3/miR-130b-3p/MBNL1 feedback loop regulated by mTORC1 signaling promotes angiogenesis and tumor growth. J Exp Clin Cancer Res. 2022;41(1):297.
  • 9 Wang H, Zhou Z, Lin W, Qian Y, He S, Wang J. MicroRNA-21 promotes head and neck squamous cell carcinoma (HNSCC) induced transition of bone marrow mesenchymal stem cells to cancer-associated fibroblasts. BMC Cancer. 2023;23(1):1135.
  • 10 Wang D, Zhao X, Li S, Guo H, Li S, Yu D. The impact of LncRNA-SOX2-OT/let-7c-3p/SKP2 Axis on head and neck squamous cell carcinoma progression: insights from bioinformatics analysis and experimental validation. Cell Signal. 2024;115(218):111018.
  • 11 Gopalakrishnan S, Pandi A, Arumugam P, Jayaseelan VP. MicroRNAs targeting CDKN2A gene as a potential prognostic marker in head and neck squamous cell carcinoma. Mol Biol Res Commun. 2024;13(1):21-7.
  • 12 Huang F, Qi H. MiR-29c-3p/C1QTNF6 Restrains the Angiogenesis and Cell Proliferation, Migration and Invasion in Head and Neck Squamous Cell Carcinoma. Mol Biotechnol. 2023;65(6):913-21.
  • 13 Peterson SM, Thompson JA, Ufkin ML, Sathyanarayana P, Liaw L, Congdon CB. Common features of microRNA target prediction tools. Front Genet. 2014;5(FEB):23.
  • 14 Nazarov PV, Kreis S. Integrative approaches for analysis of mRNA and microRNA high-throughput data. Comput Struct Biotechnol J. 2021;19:1154-62.
  • 15 Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.
  • 16 Reczko M, Maragkakis M, Alexiou P, Grosse I, Hatzigeorgiou AG. Functional microRNA targets in protein coding sequences. Bioinformatics. 2012;28(6):771-6.
  • 17 Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48(D1):D127-31.
  • 18 Rennie W, Liu C, Carmack CS, Wolenc A, Kanoria S, Lu J, et al. STarMir: a web server for prediction of microRNA binding sites. Nucleic Acids Res. 2014;42(Web Server issue):W114-8.
  • 19 Cho S, Jang I, Jun Y, Yoon S, Ko M, Kwon Y, et al. MiRGator v3.0: a microRNA portal for deep sequencing, expression profiling and mRNA targeting. Nucleic Acids Res. 2013;41(Database issue D1):D252-7.
  • 20 Andrés-León E, González Peña D, Gómez-López G, Pisano DG. miRGate: a curated database of human, mouse and rat miRNA-mRNA targets. Database (Oxford). 2015;2015:bav035.
  • 21 Wen M, Cong P, Zhang Z, Lu H, Li T. DeepMirTar: a deep-learning approach for predicting human miRNA targets. Bioinformatics. 2018;34(22):3781-7.
  • 22 Huang HY, Lin YC, Cui S, Huang Y, Tang Y, Xu J, et al. miRTarBase update 2022: an informative resource for experimentally validated miRNA-target interactions. Nucleic Acids Res. 2022;50(D1):D222-30.
  • 23 Sethupathy P, Corda B, Hatzigeorgiou AG. TarBase: A comprehensive database of experimentally supported animal microRNA targets. RNA. 2006;12(2):192-7.
  • 24 Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T. miRecords: an integrated resource for microRNA-target interactions. Nucleic Acids Res. 2009;37(Database issue):D105-10.
  • 25 Jiang Q, Wang Y, Hao Y, Juan L, Teng M, Zhang X, et al. miR2Disease: a manually curated database for microRNA deregulation in human disease. Nucleic Acids Res. 2009;37(Database issue):D98-104.
  • 26 Karagkouni D, Paraskevopoulou MD, Chatzopoulos S, Vlachos IS, Tastsoglou S, Kanellos I, et al. DIANA-TarBase v8: a decade-long collection of experimentally supported miRNA-gene interactions. Nucleic Acids Res. 2018;46 D1:D239-45.
  • 27 Brierley JD, Gospodarowicz MK, O'Sullivan B. TNM Classification of Malignant Tumors. 8th ed. Eighth Edi. Oxford, UK: Wiley Blackwell; 2017. 17p.
  • 28 Dweep H, Sticht C, Pandey P, Gretz N. miRWalk-database: prediction of possible miRNA binding sites by "walking" the genes of three genomes. J Biomed Inform. 2011;44(5):839-47.
  • 29 Ding J, Li X, Hu H. TarPmiR: a new approach for microRNA target site prediction. Bioinformatics. 2016;32(18):2768-75.
  • 30 Liu W, Wang X. Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data. Genome Biol. 2019;20(1):18.
  • 31 Yang H, Dinney CP, Ye Y, Zhu Y, Grossman HB, Wu X. Evaluation of genetic variants in microRNA-related genes and risk of bladder cancer. Cancer Res. 2008;68(7):2530-7.
  • 32 Yu T, Xu N, Haque N, Gao C, Huang W, Huang Z. Popular Computational Tools Used for miRNA Prediction and Their Future Development Prospects. Interdiscip Sci. 2020;12(4):395-413.
  • 33 Oliveira AC, Bovolenta LA, Nachtigall PG, Herkenhoff ME, Lemke N, Pinhal D. Combining results from distinct microRNA target prediction tools enhances the performance of analyses. Front Genet. 2017;8:59.
  • 34 Chen L, Heikkinen L, Wang C, Yang Y, Sun H, Wong G. Trends in the development of miRNA bioinformatics tools. Brief Bioinform. 2019;20(5):1836-52.
  • 35 Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115(7):787-98.
  • 36 Riffo-Campos ÁL, Riquelme I, Brebi-Mieville P. Tools for sequence-based miRNA target prediction: what to choose? Int J Mol Sci. 2016;17(12):1987.
  • 37 Feng H, Zhong L, Yang X, Wan Q, Pei X, Wang J. Development and validation of prognostic index based on autophagy-related genes in patient with head and neck squamous cell carcinoma. Cell Death Discov. 2020;6(1):59.
  • 38 Koshizuka K, Nohata N, Hanazawa T, Kikkawa N, Arai T, Okato A, et al. Deep sequencing-based microRNA expression signatures in head and neck squamous cell carcinoma: dual strands of pre-miR-150 as antitumor miRNAs. Oncotarget. 2017;8(18):30288-304.
  • 39 Schinke H, Shi E, Lin Z, Quadt T, Kranz G, Zhou J, et al. A transcriptomic map of EGFR-induced epithelial-to-mesenchymal transition identifies prognostic and therapeutic targets for head and neck cancer. Mol Cancer. 2022;21(1):178.
  • 40 Chen Y, Liu M, Jin H, Peng B, Dai L, Wang S, Xing H, Wang B, Wu Z. Synthetic Evaluation of MicroRNA-1-3p Expression in Head and Neck Squamous Cell Carcinoma Based on Microarray Chips and MicroRNA Sequencing. Biomed Res Int. 2021;2021:6529255.
  • 41 Zhao L, Chi W, Cao H, Cui W, Meng W, Guo W, et al. Screening and clinical significance of tumor markers in head and neck squamous cell carcinoma through bioinformatics analysis. Mol Med Rep. 2019;19(1):143-54.
  • 42 Chen YJ, Chang JT, You GR, Huang CY, Fan KH, Cheng AJ. Panel biomarkers associated with cancer invasion and prognostic prediction for head-neck cancer. Biomark Med. 2021;15(11):861-77.
  • 43 Jiang P, He S, Li Y, Xu Z. Identification of Therapeutic and Prognostic Biomarkers of Lamin C (LAMC) Family Members in Head and Neck Squamous Cell Carcinoma. Med Sci Monit. 2020;26:e925735.
  • 44 Huang C, Liang Y, Dong Y, Huang L, Li A, Du R, et al. Novel prognostic matrisome-related gene signature of head and neck squamous cell carcinoma. Front Cell Dev Biol. 2022;10:884590.
  • 45 Liu L, Jung SN, Oh C, Lee K, Won HR, Chang JW, et al. LAMB3 is associated with disease progression and cisplatin cytotoxic sensitivity in head and neck squamous cell carcinoma. Eur J Surg Oncol. 2019;45(3):359-65.
  • 46 Kinoshita T, Hanazawa T, Nohata N, Kikkawa N, Enokida H, Yoshino H, et al. Tumor suppressive microRNA-218 inhibits cancer cell migration and invasion through targeting laminin-332 in head and neck squamous cell carcinoma. Oncotarget. 2012;3(11):1386-400.
  • 47 Xie J, Huang L, Lu YG, Zheng DL. Roles of the Wnt Signaling Pathway in Head and Neck Squamous Cell Carcinoma. Front Mol Biosci. 2021;7:590912. Review.
  • 48 Chen C, Luo L, Xu C, Yang X, Liu T, Luo J, et al. Tumor specificity of WNT ligands and receptors reveals universal squamous cell carcinoma oncogenes. BMC Cancer. 2022;22(1):1-14.
  • 49 Wang N, Yan H, Wu D, Zhao Z, Chen X, Long Q, et al. PRMT5/Wnt4 axis promotes lymph-node metastasis and proliferation of laryngeal carcinoma. Cell Death Dis. 2020;11(10):864.
  • 50 Carracedo DG, Astudillo A, Rodrigo JP, Suarez C, Gonzalez MV. Skp2, p27kip1 and EGFR assessment in head and neck squamous cell carcinoma: prognostic implications. Oncol Rep. 2008;20(3):589-95.
  • 51 Yang C, Nan H, Ma J, Jiang L, Guo Q, Han L, et al. High Skp2/low p57Kip2 expression is associated with poor prognosis in human breast Carcinoma: Supplementary issue: Breast cancer detection and screening. Breast Cancer Basic Clin Res. 2015;9:13-21.
  • 52 Pomella S, Cassandri M, D'Archivio L, Porrazzo A, Cossetti C, Phelps D, et al. MYOD-SKP2 axis boosts tumorigenesis in fusion negative rhabdomyosarcoma by preventing differentiation through p57Kip2 targeting. Nat Commun. 2023;14(1):8373.
  • 53 Fan GK, Chen J, Ping F, Geng Y. Immunohistochemical analysis of P57(kip2), p53 and hsp60 expressions in premalignant and malignant oral tissues. Oral Oncol. 2006;42(2):147-53.
  • 54 Cai Y, Yang W. PKMYT1 regulates the proliferation and epithelial-mesenchymal transition of oral squamous cell carcinoma cells by targeting CCNA2. Oncol Lett. 2022;23(2):63.
  • 55 Feng CJ, Li HJ, Li JN, Lu YJ, Liao GQ. Expression of Mcm7 and Cdc6 in oral squamous cell carcinoma and precancerous lesions. Anticancer Res. 2008;28 6A:3763-9.
  • 56 Saloura V, Vougiouklakis T, Zewde M, Kiyotani K, Park JH, Gao G, et al. WHSC1L1 drives cell cycle progression through transcriptional regulation of CDC6 and CDK2 in squamous cell carcinoma of the head and neck. Oncotarget. 2016;7(27):42527-38.
  • 57 Remer E, Badarni M, Hikri E, Dayan A, Levi L, Popovtzer A, et al. CDK 4/6 Inhibition Overcomes Acquired and Inherent Resistance to PI3KKip2 Inhibition in Pre-Clinical Models of Head and Neck Squamous Cell Carcinoma. J Clin Med. 2020;9(10):3214.
  • 58 Duah E, Seligson ND, Persaud AK, Dam Q, Pabla N, Rocco JW, et al. CDK4/6 and autophagy inhibitors synergize to suppress the growth of human head and neck squamous cell carcinomas. Mol Carcinog. 2023;62(8):1201-12.
  • 59 Nakamoto M, Jin P, O'Donnell WT, Warren ST. Physiological identification of human transcripts translationally regulated by a specific microRNA. Hum Mol Genet. 2005;14(24):3813-21.
  • 60 Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D, et al. The let-7 microRNA represses cell proliferation pathways in human cells. Cancer Res. 2007;67(16):7713-22.
  • 61 Chang S, Gao L, Yang Y, Tong D, Guo B, Liu L, et al. miR-145 mediates the antiproliferative and gene regulatory effects of vitamin D3 by directly targeting E2F3 in gastric cancer cells. Oncotarget. 2015;6(10):7675-85.

Edited by

Publication Dates

  • Publication in this collection
    18 Apr 2025
  • Date of issue
    2025

History

  • Received
    02 Sept 2024
  • Accepted
    20 Oct 2024
location_on
Instituto Israelita de Ensino e Pesquisa Albert Einstein Avenida Albert Einstein, 627/701 , 05651-901 São Paulo - SP, Tel.: (55 11) 2151 0904 - São Paulo - SP - Brazil
E-mail: revista@einstein.br
rss_feed Acompanhe os números deste periódico no seu leitor de RSS
Reportar erro