Background. Papillary thyroid carcinoma (PTC) is the commonest thyroid malignancy originating from the follicle cells in the thyroid. Despite a good overall prognosis, certain high-risk cases as in those with lymph node metastasis (LNM) have progressive disease and poorer prognosis. MicroRNAs are a class of non-protein-coding, 19–24 nucleotides single-stranded RNAs which regulate gene expression and these molecules have been shown to play a role in LNM. The integrated analysis of miRNAs and gene expression profiles together with transcription factors (TFs) has been shown to improve the identification of functional miRNA-target gene-TF relationships, providing a more complete view of molecular events underlying metastasis process.Objectives. We reanalyzed The Cancer Genome Atlas (TCGA) datasets on PTC to identify differentially expressed miRNAs/genes in PTC patients with LNM-positive (LNM-P) versus lymph node negative (LNN) PTC patients and to investigate the miRNA-gene-TF regulatory circuit that regulate LNM in PTC.Results. PTC patients with LNM (PTC LNM-P) have a significantly shorter disease-free survival rate compared to PTC patients without LNM (PTC LNN) (Log-rank Mantel Cox test, p = 0.0049). We identified 181 significantly differentially expressed miRNAs in PTC LNM-P versus PTC LNN; 110 were upregulated and 71 were downregulated. The five topmost deregulated miRNAs were hsa-miR-146b, hsa-miR-375, hsa-miR-31, hsa-miR-7-2 and hsa-miR-204. In addition, 395 miRNAs were differentially expressed between PTC LNM-P and normal thyroid while 400 miRNAs were differentially expressed between PTC LNN and normal thyroid. We found four significant enrichment pathways potentially involved in metastasis to the lymph nodes, namely oxidative phosphorylation (OxPhos), cell adhesion molecules (CAMs), leukocyte transendothelial migration and cytokine–cytokine receptor interaction. OxPhos was the most significantly perturbed pathway (p = 4.70E−06) involving downregulation of 90 OxPhos-related genes. Significant interaction of hsa-miR-301b with HLF, HIF and REL/NFkB transcription factors were identified exclusively in PTC LNM-P versus PTC LNN.Conclusion. We found evidence of five miRNAs differentially expressed in PTC LNM-P. Alteration in OxPhos pathway could be the central event in metastasis to the lymph node in PTC. We postulate that hsa-miR-301b might be involved in regulating LNM in PTC via interactions with HLF, HIF and REL/NFkB. To the best of our knowledge, the roles of these TFs have been studied in PTC but the precise role of this miRNA with these TFs in LNM in PTC has not been investigated.

Lymph node metastasis (LNM) in papillary thyroid cancer (PTC) has been shown to be associated with increased risk of locoregional recurrence, poor prognosis and decreased survival, especially in older patients. Hence, there is a need for a reliable biomarker for the prediction of LNM in this cancer. MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene translation or degradation and play key roles in numerous cellular functions including cell-cycle regulation, differentiation, apoptosis, invasion and migration. Various studies have demonstrated deregulation of miRNA levels in many diseases including cancers. While a large number of miRNAs have been identified from PTCs using various means, association of miRNAs with LNM in such cases is still controversial. Furthermore, studies linking most of the identified miRNAs to the mechanism of LNM have not been well documented. The aim of this review is to update readers on the current knowledge of miRNAs in relation to LNM in PTC.

The incidence rate of papillary thyroid carcinoma (PTC) has rapidly increased in the recent decades, and the microRNA (miRNA) is one of the potential biomarkers in this cancer. Despite good prognosis, certain features such as lymph node metastasis (LNM) and BRAF V600E mutation are associated with a poor outcome. More than 50% of PTC patients present with LNM and BRAF V600E is the most common mutation identified in this cancer. The molecular mechanisms underlying these features are yet to be elucidated. This study aims to elucidate miRNA–genes interaction networks in PTC with or without LNM and to determine the association of BRAF V600E mutation with miRNAs and genes expression profiles. Next generation sequencing was performed to characterize miRNA and gene expression profiles in 20 fresh frozen tumor and the normal adjacent tissues of PTC with LNM positive (PTC LNM-P) and PTC without LNM (PTC LNN). BRAF V600E was genotyped using Sanger sequencing. Bioinformatics integration and pathway analysis were performed to determine the regulatory networks involved. Based on network analysis, we then investigated the association between miRNA and gene biomarkers, and pathway enrichment analysis was performed to study the role of candidate biomarkers. We identified 138 and 43 significantly deregulated miRNAs (adjusted p value < 0.05; log2 fold change ≤ −1.0 or ≥1.0) in PTC LNM-P and PTC LNN compared to adjacent normal tissues, respectively. Ninety-six miRNAs had significant expression ratios of 3p-to-5p in PTC LNM-P as compared to PTC LNN. In addition, ribosomal RNA-reduced RNA sequencing analysis revealed 699 significantly deregulated genes in PTC LNM-P versus normal adjacent tissues, 1,362 genes in PTC LNN versus normal adjacent tissue, and 1,576 genes in PTC LNM-P versus PTC LNN. We provide the evidence of miRNA and gene interactions, which are involved in LNM of papillary thyroid cancer. These findings may lead to better understanding of carcinogenesis and metastasis processes. This study also complements the existing knowledge about deregulated miRNAs in papillary thyroid carcinoma development.

No abstract

MicroRNAs (miRNAs) are small non-coding RNAs, which play a critical regulatory role in papillary thyroid carcinoma (PTC). BRAF V600E is a hotspot mutation occurring in a majority of PTC cases and is proposed to be associatedwith poor clinical outcomes. The relationship between BRAF V600E status and miRNA expression in PTC has not been comprehensively studied. In this study, we aimed to identify the differentially expressed miRNAs in PTCs with and without BRAF V600E in an unbiased manner. Five fresh frozen thyroid cancer tissues paired with their respective adjacent normal tissues from PTC patients were subjected to BRAF V600E genotyping using Sanger sequencing and small RNA deep sequencing (miRNAseq). MiRNAs differentially expressed between BRAF V600E-positive and BRAF V600E-negative PTC tissues were validated in silico using The Cancer Genome Atlas (TCGA) THCA datasets containing 420 samples. MiRNA target prediction and pathway enrichment analysis were performed to identify biological pathways altered in this cancer. We identified 174 differentially expressed miRNAs; 80 were significantly over-expressed, while 94 were underexpressed (adj. p-value < 0.1; log2 fold change ? -1 or ? 1). Fifteen miRNAs were significantly differentially expressed only in BRAF V600E-positive PTC, and eight of these were validated in TCGA THCA dataset (hsa-miR-212, -132, -135b-3p/5p, -200b, -200a-3p/5p, -27a-3p/5p, -29a and -1296). Subsequent analysis revealed significant enrichment of cancer-relatedpathways including proteoglycans in cancer, ECM-receptor interaction and MAPK pathways in BRAF V600E-positive PTC. Using the miRNAseq and in silico validation using TCGA THCA study, we identified eight miRNAs that were differentially expressed in PTC tissues with BRAF V600E. This study also complemented the existing knowledge about deregulated miRNAs in PTC development.