Integrated Analysis of Thyroid Cancer Public Datasets Reveals Role of Post-Transcriptional Regulation on Tumor Progression by Targeting of Immune System Mediators

DLK1-DIO3 region as a source of tumor suppressor miRNAs in papillary thyroid carcinoma

BACKGROUND

In previous studies, we demonstrated the downregulation of several miRNAs from the DLK1-DIO3 genomic region in papillary thyroid carcinoma (PTC). Due to the large number of miRNAs within this region, the individual contribution of these molecules to PTC development and progression remains unclear.

OBJECTIVE

In this study, we aimed to clarify the contribution of DLK1-DIO3-derived miRNAs to PTC.

METHODS

We used different computational approaches and in vitro resources to assess the biological processes and signaling pathways potentially modulated by these miRNAs.

RESULTS

Our analysis suggests that, out of more than 100 mature miRNAs originated from the DLK1-DIO3 region, a set of 12 miRNAs accounts for most of the impact on PTC development and progression, cooperating to modulate distinct cancer-relevant biological processes, such as cell migration, extracellular matrix remodeling, and signal transduction. The restoration of the expression of one of these miRNAs (miR-485-5p) in a BRAFT199A-positive PTC cell line impaired proliferation and migration, suppressing the expression of GAB2 and RAC1, validated miR-485-5p targets.

CONCLUSIONS

Overall, our results shed light on the role of the DLK1-DIO3 region, which harbors promising tumor suppressor miRNAs in thyroid cancer, and open prospects for the functional exploration of these miRNAs as therapeutic targets for PTC.

The EPRS-ATF4-COLI pathway axis is a potential target for anaplastic thyroid carcinoma therapy

BACKGROUND

Anaplastic thyroid carcinoma (ATC) is recognized as the most aggressive and malignant form of thyroid cancer, underscoring the critical need for effective therapeutic strategies to curb its progression and improve patient prognosis. Halofuginone (HF), a derivative of febrifugine, has displayed antitumor properties across various cancer types. However, there is a paucity of published research focused on the potential of HF to enhance the clinical efficacy of treating ATC.

OBJECTIVE

In this study, we thoroughly investigated the antitumor effects and mechanisms of HF in ATC, aiming to discover lead compounds for treating ATC and reveal novel therapeutic targets for ATC tumors.

METHODS

A series of assays, including CCK8, colony formation, tumor xenograft models, and ATC tumor organoid experiments, were conducted to evaluate the anticancer properties of HF both in vitro and in vivo. Techniques such as drug affinity responsive target stability (DARTS), western blot, immunofluorescence, and immunohistochemistry were employed to pinpoint HF target proteins within ATC. Furthermore, we harnessed the GEPIA and GEO databases and performed immunohistochemistry to validate the therapeutic potential of the glutamyl-prolyl-tRNA-synthetase (EPRS)- activating transcription factor 4 (ATF4)- type I collagen (COLI) pathway axis in the context of ATC. The study also incorporated RNA sequencing analysis, confocal imaging, and flow cytometry to delve into the molecular mechanisms of HF in ATC.

RESULTS

HF exhibited a substantial inhibitory impact on cell proliferation in vitro and on tumor growth in vivo. The DARTS results highlighted HF's influence on EPRS within ATC cells, triggering an amino acid starvation response (AASR) by suppressing EPRS expression, consequently leading to a reduction in COLI expression in ATC cells. The introduction of proline mitigated the effect of HF on ATF4 and COLI expression, indicating that the EPRS-ATF4-COLI pathway axis was a focal target of HF in ATC. Analysis of the expression levels of the EPRS, ATF4, and COLI proteins in thyroid tumors, along with an examination of the relationship between COLI expression and thyroid tumor stage, revealed that HF significantly inhibited the growth of ATC tumor organoids, demonstrating the therapeutic potential of targeting the EPRS-ATF4-COLI pathway axis in ATC. RNA sequencing analysis revealed significant differences in the pathways associated with metastasis and apoptosis between control and HF-treated cells. Transwell assays and flow cytometry experiments provided evidence of the capacity of HF to impede cell migration and induce apoptosis in ATC cells. Furthermore, HF hindered cell metastasis by suppressing the epithelial-mesenchymal transition (EMT) pathway, acting through the inhibition of FAK-AKT-NF-κB/Wnt-β-catenin signaling and restraining angiogenesis via the VEGF pathway. HF also promoted apoptosis through the mitochondrial apoptotic pathway.

CONCLUSION

This study provided inaugural evidence suggesting that HF could emerge as a promising therapeutic agent for the treatment of ATC. The EPRS-ATF4-COLI pathway axis stood out as a prospective biomarker and therapeutic target for ATC.

MicroRNAs in thyroid cancer with focus on medullary thyroid carcinoma: potential therapeutic targets and diagnostic/prognostic markers and web based tools

This review aimed to describe the inculpation of microRNAs (miRNAs) in thyroid cancer (TC) and its subtypes, mainly medullary thyroid carcinoma (MTC), and to outline web-based tools and databases for bioinformatics analysis of miRNAs in TC. Additionally, the capacity of miRNAs to serve as therapeutic targets and biomarkers in TC management will be discussed. This review is based on a literature search of relevant articles on the role of miRNAs in TC and its subtypes, mainly MTC. Additionally, web-based tools and databases for bioinformatics analysis of miRNAs in TC were identified and described. MiRNAs can perform as oncomiRs or antioncoges, relying on the target mRNAs they regulate. MiRNA replacement therapy using miRNA mimics or antimiRs that aim to suppress the function of certain miRNAs can be applied to correct miRNAs aberrantly expressed in diseases, particularly in cancer. MiRNAs are involved in the modulation of fundamental pathways related to cancer, resembling cell cycle checkpoints and DNA repair pathways. MiRNAs are also rather stable and can reliably be detected in different types of biological materials, rendering them favorable diagnosis and prognosis biomarkers as well. MiRNAs have emerged as promising tools for evaluating medical outcomes in TC and as possible therapeutic targets. The contribution of miRNAs in thyroid cancer, particularly MTC, is an active area of research, and the utility of web applications and databases for the biological data analysis of miRNAs in TC is becoming increasingly important.

MiR-495-3p regulates cell migration and invasion in papillary thyroid carcinoma

Background

Papillary thyroid carcinoma (PTC) is the most prevalent histotype of thyroid cancer and the presence of BRAFV600E mutation in these tumors is related to the malignancy and prognosis of the disease. In recent years attention has been focused on the role of microRNAs in the biology of PTC cells, especially in their role in the modulation of pathways related to tumorigenesis. DLK1-DIO3-derived miRNAs have been shown to play important roles in tumor context and are globally downregulated in PTC.

Methods

Based on a previous in silico target prediction and gene enrichment analysis, we identified miR-495-3p as the candidate with the highest tumor suppressor potential role in PTC among DLK1-DIO3-derived miRNAs. We used bioinformatics and an in vitro model of miR-495-3p overexpression to further understand the influence of this molecule on the tumorigenic processes of PTC.

Results

Overexpression of miR-495-3p impaired cell migration and invasion of PTC cells harboring the BRAFV600E mutation and affected the expression of targets predicted in the bioinformatic analysis, such as TGFB2, EREG and CCND1.

Conclusion

Overall, our results indicate that the loss of miR-495-3p expression during PTC development might play an important role in its progression.

CREB3L1 promotes tumor growth and metastasis of anaplastic thyroid carcinoma by remodeling the tumor microenvironment

Anaplastic thyroid carcinoma (ATC) is an extremely malignant type of endocrine cancer frequently accompanied by extrathyroidal extension or metastasis through mechanisms that remain elusive. We screened for the CREB3 transcription-factor family in a large cohort, consisting of four microarray datasets. This revealed that CREB3L1 was specifically up regulated in ATC tissues and negatively associated with overall survival of patients with thyroid cancer. Consistently, high expression of CREB3L1 was negatively correlated with progression-free survival in an independent cohort. CREB3L1 knockdown dramatically attenuated invasion of ATC cells, whereas overexpression of CREB3L1 facilitated the invasion of papillary thyroid carcinoma (PTC) cells. Loss of CREB3L1 inhibited metastasis and tumor growth of ATC xenografts in zebrafish and nude mouse model. Single-cell RNA-sequencing analysis revealed that CREB3L1 expression gradually increased during the neoplastic progression of a thyroid follicular epithelial cell to an ATC cell, accompanied by the activation of the extracellular matrix (ECM) signaling. CREB3L1 knockdown significantly decreased the expression of collagen subtypes in ATC cells and the fibrillar collagen in xenografts. Due to the loss of CREB3L1, ATC cells were unable to activate alpha-smooth muscle actin (α-SMA)-positive cancer-associated fibroblasts (CAFs). After CREB3L1 knockdown, the presence of CAFs inhibited the growth of ATC spheroids and the metastasis of ATC cells. Further cytokine array screening showed that ATC cells activated α-SMA-positive CAFs through CREB3L1-mediated IL-1α production. Moreover, KPNA2 mediated the nuclear translocation of CREB3L1, thus allowing it to activate downstream ECM signaling. These results demonstrate that CREB3L1 maintains the CAF-like property of ATC cells by activating the ECM signaling, which remodels the tumor stromal microenvironment and drives the malignancy of ATC.

Epigenetic signature associated with thyroid cancer progression and metastasis

Thyroid cancer is not among the top cancers in terms of diagnosis or mortality but it still ranks fifth among the cancers diagnosed in women. Infact, women are more likely to be diagnosed with thyroid cancer than the males. The burden of thyroid cancer has dramatically increased in last two decades in China and, in the United States, it is the most diagnosed cancer in young adults under the age of twenty-nine. All these factors make it worthwhile to fully understand the pathogenesis of thyroid cancer. Towards this end, microRNAs (miRNAs) have constantly emerged as the non-coding RNAs of interest in various thyroid cancer subtypes on which there have been numerous investigations over the last decade and half. This comprehensive review takes a look at the current knowledge on the topic with cataloging of miRNAs known so far, particularly related to their utility as epigenetic signatures of thyroid cancer progression and metastasis. Such information could be of immense use for the eventual development of miRNAs as therapeutic targets or even therapeutic agents for thyroid cancer therapy.

Serous BMP8A has Clinical Significance in the Ultrasonic Diagnosis of Thyroid Cancer and Promotes Thyroid Cancer Cell Progression

OBJECTIVE

This study aims to discover a potential cytokine biomarker for early diagnosis of thyroid cancer.

METHODS

We employed data mining of The Cancer Genome Atlas (TCGA) and experimentally elucidated its mechanistic contributions. The differential expression genes (DEGs) between thyroid cancer and health population were analyzed with TCGA online bioinformatic tools. The relative expression of Bone Morphogenetic Protein 8A (BMP8A) was determined by real-time PCR in ultrasonic diagnosed thyroid cancer both in vivo and in vitro. The serous BMP8A content was quantified with an ELISA kit. Protein levels of BMP8A, OCLN, ZEB1, EZH2 and β-Actin were analyzed by Western blot. Cell viability was measured by the MTT assay, and anchorage-independent growth was measured by the soft agar colony formation assay. Cell migrative and invasive capacities were interrogated with transwell chamber assays.

RESULTS

We identified aberrantly high expression of BMP8A in thyroid cancer, which was associated with unfavorable prognosis and tumor progression. The serous BMP8A was also significantly up-regulated in thyroid cancer patients. Ectopic over-expression of BMP8A remarkably stimulated cell viability and anchorage-independent growth. Meanwhile, the migrative and invasive capacities were greatly increased in response to BMP8A over-expression. Mechanistically, we characterized the positive correlation between BMP8A and TCF7L1, and forced expression of TCF7L1 induced BMP8A expression in TPC-1 cells.

CONCLUSION

In summary, we have identified a novel biomarker for early diagnosis in addition to Ultrasound for thyroid cancer, which is subjected to TCF7L1 regulation.

Integrated bioinformatics analysis for the identification of key genes and signaling pathways in thyroid carcinoma

Thyroid carcinoma (TC) is one of the most common types of endocrine neoplasm with poor prognosis due to its aggressive behavior. Biomarkers for early diagnosis and prevention of TC are in urgent demand. By using a bioinformatics analysis, the present study aimed to identify essential genes and pathways associated with TC. First, the GSE27155 and GSE50901 expression profiles were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were obtained using the two microarray datasets and further subjected to integrated analysis. A gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed 45 common DEGs in the two datasets. GO and KEGG pathway analysis indicated that the biological functions of the DEGs included protein binding, cardiac muscle cell potential involved in contraction, aldehyde dehydrogenase activity, the TGF-β receptor signaling pathway and the canonical Wnt signaling pathway. A protein-protein interaction network was also constructed and visualized to display the nodes of the top 9 up- and 36 downregulated common DEGs. The integrated bioinformatics analysis indicated that potassium inwardly rectifying channel subfamily J member 2 (KCNJ2) was the most significantly upregulated DEG. The transcriptional levels of KCNJ2 were confirmed to be elevated in TC tissues compared with those in normal tissues using reverse transcription-quantitative PCR analysis. Furthermore, the expression level of KCNJ2 was significantly associated with the 5-year survival rate of patients with TC, which was determined using the Kaplan-Meier method. In TC cell lines, KCNJ2 was also upregulated as compared with that in a normal control cell line. Finally, small interfering RNA was used to knock down the expression of KCNJ2, which was demonstrated to inhibit cell proliferation, migration and invasion, while increasing apoptosis in TC cells. In conclusion, in the present study, KCNJ2 was screened as an oncogene with a crucial role in TC development and progression and may represent a promising candidate biomarker and therapeutic target for TC.

Immunological analyses reveal an immune subtype of uveal melanoma with a poor prognosis

Uveal melanoma is an aggressive intraocular malignancy that often exhibits low immunogenicity. Metastatic uveal melanoma samples frequently exhibit monosomy 3 or BAP1 deficiency. In this study, we used bioinformatic methods to investigate the immune infiltration of uveal melanoma samples in public datasets. We first performed Gene Set Enrichment/Variation Analyses to detect immunological pathways that are altered in tumors with monosomy 3 or BAP1 deficiency. We then conducted an unsupervised clustering analysis to identify distinct immunologic molecular subtypes of uveal melanoma. We used CIBERSORT and ESTIMATE with RNA-seq data from The Cancer Genome Atlas and the GSE22138 microarray dataset to determine the sample-level immune subpopulations and immune scores of uveal melanoma samples. The Kaplan-Meier method and log-rank test were used to assess the prognostic value of particular immune cells and genes in uveal melanoma samples. Through these approaches, we discovered uveal melanoma-specific immunologic features, which may provide new insights into the tumor microenvironment and enhance the development of immunotherapies in the future.

Identification of key genes and miRNAs markers of papillary thyroid cancer

OBJECTIVE

In this study, crucial genes and microRNAs (miRNAs) associated with the progression, staging, and prognosis of papillary thyroid cancer (PTC) were identified.

METHODS

Four PTC datasets, including our own mRNA-sequencing (mRNA-seq) dataset and three public datasets downloaded from Gene Expression Omnibus and The Cancer Genome Atlas, were used to analyze differentially expressed genes (DEGs) and miRNAs (DEMs) between PTC tumor tissues and paired normal tissues (control). Gene ontology (GO) terms and pathways associated with these DEGs were identified, and protein-protein interactions (PPIs) were analyzed. Additionally, an miRNA-mRNA regulatory network was constructed and the functions of DEMs were explored. Finally, miRNAs/mRNAs associated with tumor staging and prognosis were identified. The expression levels of several key genes and miRNAs were validated by qRT-PCR.

RESULTS

Numerous DEGs and DEMs were identified between tumor and control groups in four datasets. The DEGs were significantly enriched in cell adhesion and cancer-related GO terms and pathways. In the constructed PPI network, ITGA2, FN1, ICAM1, TIMP1 and CDH2 were hub proteins. In the miRNA-mRNA negative regulatory networks, miR-204-5p regulated the largest number of target genes, such as TNFRSF12A. miR-146b, miR-204, miR-7-2, and FN1 were associated with tumor stage in PTC, and TNFRSF12A and CLDN1 were related to prognosis.

CONCLUSIONS

Our results suggested the important roles of ITGA2, FN1, ICAM1, TIMP1 and CDH2 in the progression of PTC. miR-204-5p, miR-7-2, and miR-146b are potential biomarkers for PTC staging and FN1, CLDN1, and TNFRSF12A may serve as markers of prognosis in PTC.

Down-regulation of 14q32-encoded miRNAs and tumor suppressor role for miR-654-3p in papillary thyroid cancer

Papillary thyroid carcinoma (PTC) is the most prevalent malignant neoplasia of the thyroid gland. A fraction of PTC cases show loss of differentiation and aggressive behavior, with radioiodine therapy resistance and metastasis. Although microRNAs (miRNAs) emerged as promising molecular markers for PTC, their role in the loss of differentiation observed during PTC progression remains to be fully understood. We performed the large-scale analysis of miRNA expression during PTC progression in BRAFT1799A-transgenic animals (Tg-Braf) and thyroid cancer cell lines and identified the marked downregulation of several miRNAs from the region 14q32. Data from The Cancer Genome Atlas (TCGA) confirmed the global downregulation of miRNAs from the 14q32 region in human PTC. The regulatory network potentially suppressed by these miRNAs suggests that key cancer-related biological processes such as cell proliferation, adhesion, migration and angiogenesis. Among the downregulated miRNAs, we observed that miR-654-3p levels decrease with long-term PTC progression in Tg-Braf mice and inversely correlate with EMT. The in vitro restoration of miR-654-3p decreased cell proliferation and migration and induced reprogramming of metastasis-related genes, suggesting a tumor suppressor role for this miRNA. In conclusion, we show global downregulation of 14q32-encoded miRNAs in an in vivo model of PTC progression. The potential circuitry in which these miRNAs are involved suggests that these miRNAs could play a key role in the pathophysiology of PTC and therefore be relevant for the development of new therapeutic strategies.

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.

MicroRNAs in thyroid development, function and tumorigenesis

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that modulate the vast majority of cellular processes. During development, the correct timing and expression of miRNAs in the tissue differentiation is essential for organogenesis and functionality. In thyroid gland, DICER and miRNAs are necessary for accurately establishing thyroid follicles and hormone synthesis. Moreover, DICER1 mutations and miRNA deregulation observed in human goiter influence thyroid tumorigenesis. The thyroid malignant transformation by MAPK oncogenes is accompanied by global miRNA changes, with a marked reduction of "tumor-suppressor" miRNAs and activation of oncogenic miRNAs. Loss of thyroid cell differentiation/function, and consequently iodine trapping impairment, is an important clinical characteristic of radioiodine-refractory thyroid cancer. However, few studies have addressed the direct role of miRNAs in thyroid gland physiology. Here, we focus on what we have learned in the thyroid follicular cell differentiation and function as revealed by cell and animal models and miRNA modulation in thyroid tumorigenesis.

MicroRNAs as Biomarkers in Thyroid Carcinoma

Optimal management of patients with thyroid cancer requires the use of sensitive and specific biomarkers. For early diagnosis and effective follow-up, the currently available cytological and serum biomarkers, thyroglobulin and calcitonin, present severe limitations. Research on microRNA expression in thyroid tumors is providing new insights for the development of novel biomarkers that can be used to diagnose thyroid cancer and optimize its management. In this review, we will examine some of the methods commonly used to detect and quantify microRNA in biospecimens from patients with thyroid tumor, as well as the potential applications of these techniques for developing microRNA-based biomarkers for the diagnosis and prognostic evaluation of thyroid cancers.

Integrated analysis of the potential roles of miRNA‑mRNA networks in triple negative breast cancer

Triple negative breast cancer (TNBC) is a type of breast cancer where the tumor cells are negative for the estrogen, progesterone and human epidermal growth factor 2 receptors. To date, expression profiling of microRNA (miRNA/miR) and mRNA sequences have been widely applied for the diagnosis of TNBC. In the present study, an integrated analysis of miRNA‑mRNA profiling arrays was performed. A total of five dysregulated miRNAs in patients with TNBC were identified, including upregulated miR‑558 expression and downregulated miR‑320d‑1, miR‑548v, miR‑99a and miR‑21 expression. In addition, 49 potential target mRNA sequences were identified. Bioinformatics analyses were performed on the identified miRNAs and mRNAs, including gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes pathway and miRNA‑mRNA network analyses. A total of 31 GO terms and three signaling pathways were identified. The results indicated that the differentially expressed miRNAs and their potential target mRNAs may affect the pathogenesis of TNBC, and may therefore be considered as promising biomarkers for the early diagnosis and targeted therapy of patients with TNBC.

Identification of microRNAs associated with medullary thyroid carcinoma by bioinformatics analyses

The present study aimed to investigate the microRNA (miRNA) profile in human medullary thyroid carcinoma (MTC) tissue. The GSE40807 data profile was downloaded from the Gene Expression Omnibus database. Following preprocessing, differentially expressed microRNAs (DEMs) between MTC and healthy tissues were identified. Based on the obtained DEMs, transcription factor (TF)‑miRNA and miRNA‑target gene regulatory association pairs were predicted. Finally, functional enrichment analysis was performed on target genes of DEMs. Fifteen upregulated and 17 downregulated DEMs were identified. In the constructed TF‑miRNA regulatory network, hsa‑miR‑9‑5p was regulated by 9 TFs and hsa‑miR‑1 was regulated by 8 TFs. TFs of nuclear factor of κ light polypeptide gene enhancer in B‑cells 1 (NF‑κB1) and v‑myc avian myelocytomatosis viral oncogene homolog (MYC) regulated 4 and 3 DEMs, respectively. In the miRNA‑target gene regulatory network, hsa‑miR‑1, hsa‑miR‑9‑5p, hsa‑miR‑96‑5p and hsa‑miR‑590‑5p were most upregulated. The target genes of these 4 miRNAs were primarily enriched in the mitogen activated protein kinase (MAPK) signaling pathway. Therefore, MAPK signaling pathway may serve important roles in MTC progression. In conclusion, the DEMs hsa‑miR‑1 and hsa‑miR‑9‑5p, and TFs of NF‑κB1 and MYC may be used as biomarkers for the diagnosis and treatment of MTC.

Expression of MicroRNAs in Thyroid Carcinoma

MicroRNA (miRNA) are negative regulators of gene expression and subsequent protein production. This method of action translates into regulatory control over cellular processes, including development, signaling, metabolism, and apoptosis. A broad range of miRNA are shown to have abnormal expressions in thyroid cancers which could explain the pathology of tumor oncogenesis and disease progression. A review is conducted of the current research on miRNA dysregulation in thyroid cancers, including papillary thyroid carcinoma (PTC), follicular thyroid carcinoma (FTC), anaplastic thyroid cancer (ATC), and medullary thyroid carcinoma (MTC). Dysregulated miRNA and their associated regulatory pathways are identified and their oncogenic and pathological significance are discussed.

Key genes and pathways predicted in papillary thyroid carcinoma based on bioinformatics analysis

PURPOSE

The aim of this study was to predict the key genes and pathways associated with papillary thyroid carcinoma (PTC).

METHODS

Based on the microarray data of GSE3467 from Gene Expression Omnibus database, we identified the differentially expressed genes (DEGs) between 9 PTC samples and 9 normal controls. With the identified DEGs, functional enrichment analyses were performed. Additionally, a protein-protein interaction network was constructed to screened out some key gene nodes. These key nodes were then performed clustering analysis and pathway enrichment analysis. Furthermore, human PTC-associated network was constructed based on these key genes to investigate the potential relationships between genes and PTC.

RESULTS

A total of 651 up-regulated and 692 down-regulated DEGs were identified in PTC samples compared with controls. The up-regulated DEGs, such as complement component 3 (C3), were mainly enriched in hsa04610:Complement and coagulation cascades. The down-regulated DEGs, including paired box 8 (PAX8), peroxisome proliferator-activated receptor gamma (PPARG), and cadherin 1, type 1 were found enriched in hsa05216:Thyroid cancer. Total 33 DEGs were considered as key genes, such as PAX8, PPARG and Jun proto-oncogene (JUN). Disease-associated network analysis found that 15 key genes such as JUN, PPARG and matrix metallopeptidase 9 (MMP9) were involved in this network.

CONCLUSIONS

DEGs of C3, PPARG, PAX8, JUN and MMP9 were differentially expressed in PTC samples and may be used as potential biomarkers in the diagnosis and treatment of PTC. Additionally, pathways of complement and coagulation cascades and thyroid cancer may also play important roles in the development of PTC.