Bioinformatics analyses of significant prognostic risk markers for thyroid papillary carcinoma

SOX4 is a pivotal regulator of tumorigenesis in differentiated thyroid cancer

The SOX family consists of about 20 transcription factors involved in embryonic development, reprogramming, and cell fate determination. In this study, we demonstrated that SOX4 was significantly upregulated in differentiated thyroid cancer. Immunohistochemical analysis revealed that high SOX4 expression was associated with papillary histology, extrathyroidal extension, lymph node metastasis, and advanced disease stage. Patients whose tumors exhibited high SOX4 expression had a shorter recurrence-free survival, though significance was lost in multivariate Cox regression analysis. SOX4 silencing in thyroid cancer cells slowed cell growth, attenuated clonogenicity, and suppressed anoikis resistance. Additionally, SOX4 knockdown impeded xenograft tumor growth in nude mice. Knockdown of SOX4 expression was accompanied by reduced phosphorylation of AKT and ERK. Furthermore, CRABP2 expression correlated with SOX4 expression, and SOX4 silencing decreased CRABP2 expression and its downstream effectors such as integrin β1 and β4. These results indicate that SOX4 has both prognostic and therapeutic implications in differentiated thyroid cancer, and targeting SOX4 may modulate tumorigenic processes in the thyroid.

The role and regulation of Maf proteins in cancer

The Maf proteins (Mafs) belong to basic leucine zipper transcription factors and are members of the activator protein-1 (AP-1) superfamily. There are two subgroups of Mafs: large Mafs and small Mafs, which are involved in a wide range of biological processes, such as the cell cycle, proliferation, oxidative stress, and inflammation. Therefore, dysregulation of Mafs can affect cell fate and is closely associated with diverse diseases. Accumulating evidence has established both large and small Mafs as mediators of tumor development. In this review, we first briefly describe the structure and physiological functions of Mafs. Then we summarize the upstream regulatory mechanisms that control the expression and activity of Mafs. Furthermore, we discuss recent studies on the critical role of Mafs in cancer progression, including cancer proliferation, apoptosis, metastasis, tumor/stroma interaction and angiogenesis. We also review the clinical implications of Mafs, namely their potential possibilities and limitations as biomarkers and therapeutic targets in cancer.

Comprehensive Pan-Cancer Analysis and the Regulatory Mechanism of ASF1B, a Gene Associated With Thyroid Cancer Prognosis in the Tumor Micro-Environment

Background

The incidence of thyroid cancer, whose local recurrence and metastasis lead to death, has always been high and the pathogenesis of papillary thyroid carcinoma (PTC) has not been clearly elucidated. Therefore, the research for more accurate prognosis-related predictive biomarkers is imminent, and a key gene can often be a prognostic marker for multiple tumors.

Methods

Gene expression profiles of various cancers in the TCGA and GTEx databases were downloaded, and genes significantly associated with the prognosis of THCA were identified by combining differential analysis with survival analysis. Then, a series of bioinformatics tools and methods were used to analyze the expression of the gene in each cancer and the correlation of each expression with prognosis, tumor immune microenvironment, immune neoantigens, immune checkpoints, DNA repair genes, and methyltransferases respectively. The possible biological mechanisms were also investigated by GSEA enrichment analysis.

Results

656 differentially expressed genes were identified from two datasets and 960 DEGs that were associated with disease-free survival in THCA patients were screened via survival analysis. The former and the latter were crossed to obtain 7 key genes, and the gene with the highest risk factor, ASF1B, was selected for this study. Differential analysis of multiple databases showed that ASF1B was commonly and highly expressed in pan-cancer. Survival analysis showed that high ASF1B expression was significantly associated with poor patient prognosis in multiple cancers. In addition, ASF1B expression levels were found to be associated with tumor immune infiltration in THCA, KIRC, LGG, and LIHC, and with tumor microenvironment in BRCA, LUSC, STAD, UCEC, and KIRC. Further analysis of the relationship between ASF1B expression and immune checker gene expression suggested that ASF1B may regulate tumor immune patterns in most tumors by regulating the expression levels of specific immune checker genes. Finally, GSEA enrichment analysis showed that ASF1B high expression was mainly enriched in cell cycle, MTORC1 signaling system, E2F targets, and G2M checkpoints pathways.

Conclusions

ASF1B may be an independent prognostic marker for predicting the prognosis of THCA patients. The pan-cancer analysis suggested that ASF1B may play an important role in the tumor micro-environment and tumor immunity and it has the potential of serving as a predictive biomarker for multiple cancers.

hsa_circ_0001018 promotes papillary thyroid cancer by facilitating cell survival, invasion, G1/S cell cycle progression, and repressing cell apoptosis via crosstalk with miR-338-3p and SOX4

We identified a novel interactome, circ_0001018/miR-338-3p/SOX4, in papillary thyroid cancer (PTC), and we intended to confirm the regulatory relationship between the three and to study the effects of the three in PTC. The bioinformatics method was used to screen out the circular RNA and mRNA of interest. A cellular fractionation assay and fluorescence in situ hybridization (FISH) assay were conducted to prove that circ_0001018 and CCT4 (the host gene of circ_0001018) mRNA primarily localized in the cytoplasm of PTC cell lines. By qRT-PCR analysis, the expression of circ_0001018 and SOX4 mRNA was found upregulated while the expression of miR-338-3p was found downregulated in PTC tissues and cells. circ_0001018 silence significantly inhibited the tumor growth in xenografted nude mice. A series of cytological experiments such as a Cell Counting Kit-8 (CCK-8) assay, a 5-ethynyl-2′-deoxyuridine (EdU) assay, cell cycle profiling, wound healing, a transwell assay, and cell apoptosis were conducted and showed that circ_0001018 and SOX4 promoted cell proliferation, migration, and invasion, inhibited cell apoptosis, and reduced the cell cycle arrest at the G₁ phase in PTC cells. Compared with circ_0001018 and SOX4, miR-338-3p held the opposite function. The regulatory relationship between circ_0001018 and miR-338-3p, and between miR-338-3p and SOX4 mRNA, was validated using a luciferase reporter gene assay and/or RNA immunoprecipitation (RIP assay). Our findings showed that circ_0001018 acted as the tumor promoter via sponging miR-338-3p to elevate SOX4 expression level in PTC. Importantly, this novel circ_0001018/miR-338-3p/SOX4 axis has the potential to be considered as a therapy target for PTC.

NF-κB-Activated miR-574 Promotes Multiple Malignant and Metastatic Phenotypes by Targeting BNIP3 in Thyroid Carcinoma

Thyroid cancer is the most common endocrine malignancy, and miR-574 is significantly upregulated in thyroid cancer. However, the role and underlying mechanism of miR-574 in thyroid cancer development are poorly understood. In this study, we showed that NF-κB/p65 signaling pathway was activated and miR-574 was upregulated in thyroid cancer cells. p65 directly bound to the promoter of miR-574 and activated miR-574 transcription. Functionally, miR-574 inhibited apoptosis, promoted proliferation and migration of thyroid cancer cells, and stimulated thyroid cancer-induced tube formation of endothelial cells. On the molecular level, miR-574 inhibited the expression of BCL2/adenovirus E1B 19 kDa protein-interacting protein 3 (BNIP3) by binding to 3'-UTR of BNIP3. miR-574 also downregulated the expression of apoptosis-inducing factor (AIF), while elevated the levels of MMP2, MMP9, and VEGFA. In vivo, miR-574 promoted xenograft growth, which was associated with reduced apoptosis and enhanced angiogenesis. NF-κB/miR-574 signaling presents multiple oncogenic activities on thyroid cancer development by directly regulating the BNIP3/AIF pathway. Therefore, targeting NF-κB/miR-574 signaling may reduce the aggressiveness of thyroid cancer. IMPLICATIONS: miR-574, directly regulated by NF-κB/p65, promotes tumorigenesis of thyroid cancer via inhibiting BNIP3/AIF pathway.

Bioinformatics analysis to screen key genes in papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common type of thyroid carcinoma, and its incidence has been on the increase in recent years. However, the molecular mechanism of PTC is unclear and misdiagnosis remains a major issue. Therefore, the present study aimed to investigate this mechanism, and to identify key prognostic biomarkers. Integrated analysis was used to explore differentially expressed genes (DEGs) between PTC and healthy thyroid tissue. To investigate the functions and pathways associated with DEGs, Gene Ontology, pathway and protein-protein interaction (PPI) network analyses were performed. The predictive accuracy of DEGs was evaluated using the receiver operating characteristic (ROC) curve. Based on the four microarray datasets obtained from the Gene Expression Omnibus database, namely GSE33630, GSE27155, GSE3467 and GSE3678, a total of 153 DEGs were identified, including 66 upregulated and 87 downregulated DEGs in PTC compared with controls. These DEGs were significantly enriched in cancer-related pathways and the phosphoinositide 3-kinase-AKT signaling pathway. PPI network analysis screened out key genes, including acetyl-CoA carboxylase beta, cyclin D1, BCL2, and serpin peptidase inhibitor clade A member 1, which may serve important roles in PTC pathogenesis. ROC analysis revealed that these DEGs had excellent predictive performance, thus verifying their potential for clinical diagnosis. Taken together, the findings of the present study suggest that these genes and related pathways are involved in key events of PTC progression and facilitate the identification of prognostic biomarkers.

Prediction of novel target genes and pathways involved in tall cell variant papillary thyroid carcinoma

BACKGROUND

Tall cell variant papillary thyroid carcinoma (TCPTC) is reportedly associated with aggressive clinicopathological parameters and poor outcomes; however, the molecular mechanisms underlying TCPTC remain poorly understood.

METHODS

The gene mutation types and mRNA expression profiles of patients with TCPTC were obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed genes (DEGs) were identified. Pathways in the interaction network and the diagnostic approaches of candidate markers for TCPTC were investigated.

RESULTS

BRAF mutation was particularly prevalent in TCPTC with a mutation frequency of 78%. TCPTC was associated with a patient age >45 years, tumor multifocality, extrathyroidal extension, a higher T stage, advanced AJCC TNM stages, BRAF V600E mutation, and poor disease-free survival. We identified 4138 TCPTC-related DEGs and 301 TCPTC-specific DEGs. Intriguingly, the gene expression pattern revealed that the dysregulated levels of both putative oncogenes and tumor suppressors in TCPTC were higher than those in classical/conventional variant PTC (cPTC). Functional enrichment analyses revealed that these DEGs were involved in several cancer-related pathways. A protein-protein interaction (PPI) network was constructed from the 301 TCPTC-specific DEGs, and 3 subnetworks, and 8 hub genes were verified. Receiver operating characteristic (ROC) analyses revealed that 6 hub genes, including COL5A1, COL1A1, COL10A1, COL11A1, CCL20, and CXCL5, could be used not only for the differential diagnosis of PTC from normal samples, but also for the differential diagnosis of TCPTC from cPTC samples.

CONCLUSIONS

Our study might provide further insights into the investigations of the tumorigenesis mechanism of TCPTC and assists in the discovery of novel candidate diagnostic markers for TCPTC.

Genetic expression profile-based screening of genes and pathways associated with papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common subtype of thyroid cancer; however, the specific genes and signaling pathways involved in this cancer remain largely unclear. The present study analyzed three profile datasets, GSE6004, GSE29265 and GSE60542, which were comprised of 47 PTC and 41 normal thyroid tissue samples, to identify key genes and pathways associated with PTC. Initially, differentially-expressed genes (DEGs) between PTC and normal thyroid tissue were screened using R 3.4.0 (2017-04-21, R Foundation, Vienna, Austria, http://www.R-project.org/). These DEGs were then clustered by gene ontology functional terms and representative signaling pathways. Additionally, specific key gene nodes were filtered out from a constructed protein-protein interaction (PPI) network. The results identified a total of 423 shared DEGs associated with PTC, including 211 upregulated and 212 downregulated genes. These 423 genes were primarily enriched in glycosaminoglycan binding, sulfur compound binding, heparin binding, enzyme activator activity, peptidase activator activity and hsa04512: Extracellular matrix (ECM)-receptor interaction. A total of 21 central node genes were identified as key genes in the PTC disease process including complement factor D (CFD), Collagen Type I α 1 Chain (COL1A1), Extracellular Matrix Protein 1 (ECM1) and Fibronectin 1 (FN1). These genes are involved in protease binding, G-protein coupled receptor binding, extracellular matrix structural constituent and peptidase regulator activity. To conclude, using bioinformatics analysis, the present study identified candidate DEGs and critical pathways in PTC that may improve the current understanding regarding the underlying mechanisms of PTC. These genes and pathways may be used as potential therapeutic targets of PTC in the future.

Identification of key genes of papillary thyroid cancer using integrated bioinformatics analysis

OBJECTIVE

To identify novel clinically relevant genes in papillary thyroid carcinoma from public databases.

METHODS

Four original microarray datasets, GSE3678, GSE3467, GSE33630 and GSE58545, were downloaded. Differentially expressed genes (DEGs) were filtered from integrated data. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed, followed by protein-protein interaction (PPI) network construction. The CentiScape pug-in was performed to scale degree. The genes at the top of the degree distribution (≥ 95% percentile) in the significantly perturbed networks were defined as central genes. UALCAN and The Cancer Genome Atlas Clinical Explorer were used to verify clinically relevant genes and perform survival analysis.

RESULT

225 commonly changed DEGs (111 up-regulated and 114 down-regulated) were identified. The DEGs were classified into three groups by GO terms. KEGG pathway enrichment analysis showed DEGs mainly enriched in the PI3K-Akt signaling pathway, pathways in cancer, focal adhesion and proteoglycans in cancer. DEGs' protein-protein interaction (PPI) network complex was developed; six central genes (BCL2, CCND1, FN1, IRS1, COL1A1, CXCL12) were identified. Among them, BCL2, CCND1 and COL1A1 were identified as clinically relevant genes.

CONCLUSION

BCL2, CCND1 and COL1A1 may be key genes for papillary thyroid carcinoma. Further molecular biological experiments are required to confirm the function of the identified genes.

Inhibitory roles of miR‑9 on papillary thyroid cancer through targeting BRAF

MicroRNA‑9 (miR‑9) is reported to be underexpressed in papillary thyroid carcinoma (PTC) tissues; however, the molecular mechanisms underlying the implication of miR‑9 in PTC have yet to be elucidated. The present study aimed to explore the potential roles of miR‑9 in PTC. PTC tissue samples and paired non‑cancerous adjacent tissues were collected from 60 patients with PTC. The human TPC‑1 thyroid gland papillary carcinoma cell line was used to investigate the molecular mechanisms underlying the roles of miR‑9 in PTC. The levels of miR‑9 and its downstream target gene BRAF were detected through reverse transcription‑quantitative polymerase chain reaction. MTT assay and flow cytometry were performed to evaluate cell viability and apoptosis, respectively. A mouse xenograft tumor model was established to observe the effects of miR‑9 on thyroid gland tumorigenesis in vivo. The present study revealed that the expression of miR‑9 was significantly reduced in PTC tissues compared with paired normal tissues. In addition, miR‑9 upregulation suppressed the expression of BRAF in TPC‑1 cells in vitro. Luciferase reporter assay demonstrated that BRAF may be a direct target gene of miR‑9 in TPC‑1 cells. In addition, following transfection with miR‑9 mimics, the viability of TPC‑1 cells was suppressed and their apoptosis was enhanced; conversely, transfection with miR‑9 inhibitor exerted the opposite effects in vitro. miR‑9 overexpression or downregulation also affected in vivo PTC tumorigenesis in athymic mice. The present findings suggested that miR‑9 may suppress the viability of PTC cells and inhibit tumor growth through directly targeting the expression of BRAF in PTC.

Reactivity of thyroid papillary carcinoma cells to thyroid stimulating hormone-dominated endocrine therapy

This study investigated the effect of thyroid stimulating hormone (TSH) on the proliferation of papillary thyroid carcinoma (PTC) cells and the therapeutic effect of levothyroxine sodium (TH). PTC cells (TPC-1) were cultured using 0.1, 1.0 and 10 U/l TSH and 10⁻², 10⁻⁴ and 10⁻⁶ mol/l TH. After the appropriate concentration was screened, TPC-1 cells were further divided into control group, TSH group, TH group and TSH+TH group. The cell proliferation was detected via methyl thiazolyl tetrazolium (MTT) method, TPC-1 cell cycle was detected via flow cytometer, and the mRNA and protein expression of cyclin D1 were detected via real-time polymerase chain reaction (PCR) and enzyme-linked immunosorbent assay (ELISA). Compared with control group, TSH significantly promoted the proliferation of TPC-1 cells (P<0.05 or P<0.01), obviously promoted the transition of TPC-1 cells from G1 phase to S phase (P<0.01) and remarkably increased the mRNA and protein expression of cyclin D1 (P<0.01); but TH had a significant inhibitory effect on these results of TSH (P<0.05 or P<0.01). TSH can promote the proliferation of PTC cells, and the appropriate complement of TH can inhibit its proliferation.

ZNRF3 is downregulated in papillary thyroid carcinoma and suppresses the proliferation and invasion of papillary thyroid cancer cells

Zinc and ring finger 3 (ZNRF3) is a transmembrane E3 ubiquitin ligase that has emerged as an important regulator of cancer development; however, its cancer-related function remains controversial. Here, we investigated the possible role of ZNRF3 in thyroid carcinoma (TC). We found that ZNRF3 is downregulated in papillary thyroid carcinoma (PTC) compared to normal thyroid tissues and inversely correlated with the degree of cell differentiation. Overexpression of ZNRF3 significantly suppressed cell malignant behaviors, including cell proliferation, migration, and invasion in vitro, as well as tumor growth in vivo. Consistent with recent studies showing that ZNRF3 is involved in the Wnt/β-catenin pathway, ZNRF3 overexpression negatively regulated β-catenin activation, modulating PTC cell behaviors. Clinical specimens revealed a significant inverse correlation between ZNRF3 and β-catenin mRNA levels. Taken together, these results provide insight into a potential tumor suppressor role of ZNRF3 in PTC progression, and may have potential clinical relevance for the prognosis and treatment of PTC.