Integrated microRNA, gene expression and transcription factors signature in papillary thyroid cancer with lymph node metastasis

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

Integrated Characterization of MicroRNA and mRNA Transcriptome in Papillary Thyroid Carcinoma

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

Epigenome-wide DNA methylation profiling in colorectal cancer and normal adjacent colon using Infinium Human Methylation 450K

The aims were to profile the DNA methylation in colorectal cancer (CRC) and to explore cancer-specific methylation biomarkers. Fifty-four pairs of CRCs and the adjacent normal tissues were subjected to Infinium Human Methylation 450K assay and analysed using ChAMP R package. A total of 26,093 differentially methylated probes were identified, which represent 6156 genes; 650 probes were hypermethylated, and 25,443 were hypomethylated. Hypermethylated sites were common in CpG islands, while hypomethylated sites were in open sea. Most of the hypermethylated genes were associated with pathways in cancer, while the hypomethylated genes were involved in the PI3K-AKT signalling pathway. Among the identified differentially methylated probes, we found evidence of four potential probes in CRCs versus adjacent normal; HOXA2 cg06786372, OPLAH cg17301223, cg15638338, and TRIM31 cg02583465 that could serve as a new biomarker in CRC since these probes were aberrantly methylated in CRC as well as involved in the progression of CRC. Furthermore, we revealed the potential of promoter methylation ADHFE1 cg18065361 in differentiating the CRC from normal colonic tissue from the integrated analysis. In conclusion, aberrant DNA methylation is significantly involved in CRC pathogenesis and is associated with gene silencing. This study reports several potential important methylated genes in CRC and, therefore, merit further validation as novel candidate biomarker genes in CRC

table_1_Integrated Characterization of MicroRNA and mRNA Transcriptome in Papillary Thyroid Carcinoma.xlsx

<p>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.</p